- Journal List
- Front Plant Sci
- PMC10847352
Learn more: PMC Disclaimer | PMC Copyright Notice
Abstract
Growing evidence underscores the role of nutrients and fertigation systems in soilless production, influencing medicinal cannabis biomass and secondary metabolite content. This study delves into the impact of enhanced nutrient regimes on the ‘ionome’ and its ramifications for biomass and cannabinoid production in medicinal cannabis, comparing two distinct fertigation systems: recirculation and drain-to-waste. Notably, we assess the optimal harvest time for maximizing profitability. In comparing the experimental variant with elevated levels of phosphorus (P), potassium (K), and iron (Fe) in the nutrient solution to the control variant, we observe distinct patterns in element composition across stems, leaves, and flowers, with significant differences between fertigation systems. Total nitrogen content was determined through the Kjeldahl method. Flame atomic absorption spectrometry (FAAS) and inductively coupled plasma optical emission spectrometry (ICP-OES) were employed for elemental analysis. Cannabinoid identification and quantification used high-performance liquid chromatography with a diode-array detector (HPLC/DAD). Followed statistical analyses included ANOVA and Tukey’s HSD test. Although the augmented nutrient regimen does not substantially increase plant biomass, interesting differences emerge between the two fertigation systems. The recirculation fertigation system proves more profitable during the recommended harvest period. Nonetheless, the altered nutrient regime does not yield statistically significant differences in final inflorescence harvest mass or cannabinoid concentrations in medicinal cannabis. The choice of fertigation system influences the quantity and quality of harvested inflorescence. To optimize the balance between the dry biomass yield of flowers and cannabinoid concentration, primarily total THC yield (sum of tetrahydrocannabinolic acid, Δ9-tetrahydrocannabinol, and Δ8-tetrahydrocannabinol), we propose the 11th week of cultivation as the suitable harvest time for the recirculation system. Importantly, the recirculation system consistently outperformed the drain-to-waste system, especially after the ninth week, resulting in significantly higher total THC yields. Enriched nutrition, when compared with control, increased THC yield up to 50.7%, with a remarkable 182% surge in the recirculation system when compared with the drain-to-waste system.
1. Introduction
The pharmaceutical industry faces unique challenges in large-scale cultivation and quality control of plant-based medications, especially with Cannabis sativa L., which is subject to varying international regulations (Chandra et al., 2017; Malík et al., 2021). Indoor cultivation has evolved through techniques like ‘sinsemilla’ (cultivating non-pollinated female plants), cuttings from superior mother plants, and hydroponic systems (Moeller and Lindholst, 2014). In vitro plant tissue culture techniques have emerged as a helpful approach (Hesami et al., 2021; Král et al., 2022; Šenkyřík et al., 2023), allowing for large-scale production of genetically identical plants in controlled laboratory conditions (Lata et al., 2016). These methods, combined with intensive breeding, have increased yields of female inflorescences and their THC (Δ9-tetrahydrocannabinol) levels (Toonen et al., 2006) and improved cannabinoid profile uniformity by altering plant architecture (Danziger and Bernstein, 2021). Indoor cultivation now incorporates automated systems (Malík et al., 2022), with hydroponics gaining popularity for its potential to yield high-quality plant material (Bouchard and Dion, 2009; Vanhove et al., 2011).
Environmental factors, particularly nutrient availability, significantly affect the quantity and quality of secondary metabolites such as cannabinoids and terpenes in cannabis (Janatová et al., 2018). Nutrients like nitrogen (N), phosphorus (P), potassium (K), and iron (Fe) play pivotal roles in plant growth and secondary metabolism (Bernstein et al., 2019b; Bevan et al., 2021; Kumar et al., 2021). Phosphorus, for instance, influences root growth, flower and seed production, and stem strength (Zheng, 2022) and affects terpene profiles in aromatic plants (Rioba et al., 2015). While inadequate P supplementation can lead to deficiency symptoms, excessive supply often accumulates in roots (Shiponi and Bernstein, 2021a; Llewellyn et al., 2023). The impact of increased phosphorus supplementation on medicinal cannabis varies across organs and compounds, with key cannabinoids remaining unaffected in upper flowers (Bernstein et al., 2019b; Shiponi and Bernstein, 2021b).
Potassium regulates water and nutrient transport, cell turgor pressure, disease resistance, stem strength, and inflorescence quality and yield (Oosterhuis et al., 2014; Johnson et al., 2022). It influences various secondary metabolites, such as phenolic compounds (Nguyen et al., 2010; Varo et al., 2022), flavonoids (Chrysargyris et al., 2017; Gaaliche et al., 2019), carotenoids (San Martín-Hernández et al., 2021), and organic acids (Naumann et al., 2020), but excessive potassium can reduce secondary metabolites like cannabinoids and terpenoids (Saloner and Bernstein, 2022). Optimizing K supply is crucial to maintaining medicinal cannabis desired functionality, yield, and secondary metabolite profiles. Iron facilitates oxidation-reduction and electron transfer reactions, activates enzymes, and is essential for photosynthesis and respiration (Briat et al., 2015; Zheng, 2022). Fe deficiency leads to chlorosis in young leaves, impaired root development, yield reduction, and compromised nutritional quality in crops (Kumar et al., 2021; Llewellyn et al., 2023). Fe concentrations in plant nutrition may have varying effects on secondary metabolite production (Shi et al., 2018; Chaouqi et al., 2023), but no specific study has explored the influence of Fe nutrition in indoor medicinal cannabis cultivation. Maintaining appropriate rootzone pH is crucial for nutrient availability, microorganism activity, and root development, affecting water and nutrient uptake (Zheng, 2021). The recommended pH range in hydroponic culture is typically between 5.5–6.0 (Velazquez, 2013), while soilless production suggests 5.5–6.5 (Zheng, 2022). If the pH in the soilless production drops below 5.5, there is a risk of toxicity due to excessive levels of manganese (Mn), while a pH higher than 6.5 can result in limited availability of essential elements such as P, Fe, and Mn for plant uptake (Balliu et al., 2021). pH levels affect nutrient availability and variations between recirculation and drain-to-waste systems, with recirculation systems exhibiting fluctuations (Malík et al., 2023).
This investigation examines the physiological and chemical reactions of medicinal cannabis plants when exposed to varying concentrations of phosphorus (P), potassium (K), and iron (Fe) in the nutrient solution. The study compares the results obtained under two distinct fertigation systems. In light of these objectives, the following hypotheses were put forth: (1) Alterations in the concentrations of P, K, and Fe in the nutrient solution can lead to different inflorescence yields of medicinal cannabis plants; (2) changes in the cannabinoid profile of indoor-cultivated cannabis plants; (3) these induced changes will exhibit a correlation with the content of macro-elements and micro-element (Fe) within different plant organs, including leaves, stems, and flowers; (4) furthermore, the induced changes will display variability between different nutrition systems, namely recirculation and drain-to-waste approaches.
2. Materials and methods
2.1. Parameters of the cultivation room
Cannabis plants were cultivated in a controlled environment within growing tables in a soilless system. A cultivation chamber has an area of 15 m2 (3 × 5 m), with the actual cultivation space delineated by four cultivation tables, covering a total area of 8 m2. The tables, measuring 2 m2 (1 × 2 m), were designated for individual experiments and featured independent 100-litre tanks for nutrient solutions. These tanks were constructed using food-grade inert plastic materials. Each table accommodated a maximum of 55 black conical, square pots made of polypropylene (PP). These pots had a volume of 3.45 liters and dimensions: top – 15 cm × 15 cm, base – 11.5 cm × 11.5 cm, height – 20 cm. Capillaries were used for fertigation, ensuring each plant received individualized watering through a needle applicator. The fertigation system was programmed to provide nine cycles per day, with each cycle lasting 60 seconds. During each cycle, 94 mL of the nutrient solution was delivered to each plant (equivalent to 846 mL per plant per day). In the recirculation system, the plant drainage was piped back to the storage tank. In the drain-to-waste system, the used nutrient solution was directed to a separate waste tank and not mixed with the original solution. Microclimate parameters, including relative humidity, temperature, and CO2 levels, were regulated and monitored by an air ventilation unit. A methane-burning generator facilitated CO2 enrichment in the growth environment (550 ppm). Six double-ended high-pressure sodium lamps provided illumination with a power output of 1,000 W, delivering a suitable light spectrum for plant growth. The lamps generated a photosynthetic photon flux density (PPFD) of 1,029 μmol·m−2·s−1, with a total power of 6,000 W. The light conditions were continuously recorded using a data logger, capturing measurements every minute.
2.2. Plant material and cultivation conditions
The plant material utilized in this study consisted of cuttings (apical parts with at least 3 fully expanded leaves) obtained from C. sativa L. ‘McLove’ mother plants (Atherton and Li, 2023). These plants belong to chemotype I, characterized by a high ratio of Δ9-tetrahydrocannabinolic acid:cannabidiolic acid (THCA : CBDA > 1.0) (Handbook of cannabis, 2016). The mother plants were carefully maintained in a dedicated separate growth chamber under controlled conditions. A total of 220 cuttings were prepared. These cuttings were then cultivated for a period of 21 days in rock-wool cubes measuring 4 × 4 cm. Once the cuttings had developed roots, they were transferred to a separate cultivation room and placed in 3.45-liter pots filled with three liters of expanded clay growing medium. During the vegetative phase, the light regime consisted of 18 hours of light and 6 hours of darkness, with a temperature of 25°C, a relative humidity of 60%, and a CO2 concentration of 550 ppm (1.065 mg·L−1). During the dark phase, the temperature was lowered to 22°C while maintaining the same humidity level. The vegetative phase lasted 7 days, after which the cultivation conditions were adjusted for the generative phase. During the generative phase, the light period was set to 12 hours of light and 12 hours of darkness. The temperature and CO2 concentration remained the same as during the vegetative phase, while the relative humidity was reduced to 40%. Starting from the 10th week, the plants were irrigated only with demineralized (DM) water.
2.3. Treatments
In the cultivation room, the plants were divided into 4 tables, each containing 55 plants (27.5 plants/m2), as shown in Figure 1. Plants were grown in two fertigation systems. The first cycle (1C) involved the recirculation of the nutrient solution, while the second cycle (2C) employed the drain-to-waste system. At the same time, the plants were divided into two groups: control treatment (CN) and enhanced treatment (ET). Compared to CN, the ET underwent an augmented nutritional regime. For the 1C, the nutrient solution was prepared by mixing DM water every 7 days, starting from the experiment’s first day. Since the 10th week, plants were irrigated solely with DM. Solution pH was adjusted to 5.9, and EC values were recorded immediately after preparing the fresh solution and remeasured on the last day before its replacement. A sample of the nutrient solution was collected after each preparation. For the 2C, everything remained the same as for 1C, except the stock solution was prepared after depleting the stock tank, approximately every 2nd to 3rd day. Nutrient content was increased according to the age of plants. Table 1 presents the composition of the nutrient solution used for the CN as determined through measurement, and descriptions of the provided nutrients for ET can be found in Table 2. In contrast to the CN variant, the content of P, K, and Fe increased by an average of 83% P, 39% K, and 860% Fe from the 5th week to the 9th week in the ET variant.
Table 1
Week | 1 | 2 | 3 | 4 | 5 | 6–9 | 10–13 |
---|---|---|---|---|---|---|---|
N | 101 ± 1.6 | 116 ± 1.9 | 130 ± 1.8 | 150 ± 1.9 | 130 ± 1.8 | 150 ± 1.9 | DM |
P | 32. ± 0.8 | 39 ± 0.8 | 44 ± 0.6 | 52 ± 0.8 | 44 ± 0.6 | 52 ± 0.8 | DM |
K | 125 ± 1.9 | 151 ± 1.4 | 173 ± 1.9 | 193 ± 1.6 | 173 ± 1.9 | 193 ± 1.6 | DM |
Ca | 99 ± 1.3 | 119 ± 1.4 | 132 ± 1.4 | 146 ± 1.2 | 132 ± 1.4 | 146 ± 1.3 | DM |
Mg | 25 ± 0.4 | 31 ± 0.4 | 35 ± 0.5 | 39 ± 0.5 | 35 ± 0.5 | 39 ± 0.5 | DM |
S | 22 ± 0.3 | 27 ± 0.3 | 31 ± 0.3 | 35 ± 0.4 | 31 ± 0.3 | 35 ± 0.4 | DM |
Fe | 0.9 ± 0.09 | 1.1 ± 0.09 | 1.2 ± 0.11 | 1.4 ± 0.08 | 1.2 ± 0.11 | 1.4 ± 0.08 | DM |
Mn | 0.7 ± 0.07 | 0.7 ± 0.05 | 0.8 ± 0.08 | 0.9 ± 0.07 | 0.8 ± 0.08 | 0.9 ± 0.07 | DM |
Zn | 0.2 ± 0.03 | 0.3 ± 0.03 | 0.3 ± 0.04 | 0.3 ± 0.03 | 0.3 ± 0.04 | 0.3 ± 0.03 | DM |
Cu | 0.1 ± 0.01 | 0.1 ± 0.01 | 0.1 ± 0.01 | 0.1 ± 0.02 | 0.1 ± 0.01 | 0.1 ± 0.02 | DM |
B | 0.1 ± 0.02 | 0.2 ± 0.01 | 0.2 ± 0.02 | 0.2 ± 0.02 | 0.2 ± 0.02 | 0.2 ± 0.02 | DM |
Mo | 0.01 ± 0.00 | 0.02 ± 0.00 | 0.02 ± 0.00 | 0.02 ± 0.00 | 0.02 ± 0.00 | 0.02 ± 0.00 | DM |
EC | 0.97 ± 0.01 | 1.19 ± 0.01 | 1.46 ± 0.01 | 1.74 ± 0.01 | 1.46 ± 0.01 | 1.74 ± 0.01 | DM |
DM, demineralized water.
Table 2
Week | 1 | 2 | 3 | 4 | 5 | 6–9 | 10–13 |
---|---|---|---|---|---|---|---|
N | 100 ± 0.9 | 115 ± 1.0 | 129 ± 0.1 | 150 ± 1.8 | 129 ± 1.7 | 150 ± 1.8 | DM |
P | 32 ± 0.5 | 39 ± 0.4 | 42 ± 0.5 | 51 ± 0.9 | 92 ± 1.9 | 93 ± 1.9 | DM |
K | 125 ± 1.6 | 143 ± 1.2 | 172 ± 1.9 | 194 ± 1.4 | 258 ± 2.5 | 266 ± 2.7 | DM |
Ca | 98 ± 0.9 | 118 ± 1.2 | 132 ± 1. | 146 ± 1.2 | 133 ± 1.1 | 144 ± 1.4 | DM |
Mg | 25 ± 0.2 | 29 ± 0.4 | 35 ± 0.4 | 39 ± 0.4 | 33 ± 0.4 | 39 ± 0.7 | DM |
S | 22 ± 0.1 | 26 ± 0.9 | 31 ± 0.4 | 34 ± 0.3 | 31 ± 0.4 | 35 ± 0.3 | DM |
Fe | 0.9 ± 0.01 | 1.1 ± 0.06 | 1.3 ± 0.05 | 1.3 ± 0.08 | 12.3 ± 0.4 | 13.8 ± 0.9 | DM |
Mn | 0.7 ± 0.06 | 0.7 ± 0.04 | 0.7 ± 0.09 | 0.9 ± 0.04 | 0.8 ± 0.06 | 0.9 ± 0.07 | DM |
Zn | 0.2 ± 0.04 | 0.3 ± 0.03 | 0.3 ± 0.04 | 0.3 ± 0.02 | 0.3 ± 0.04 | 0.3 ± 0.03 | DM |
Cu | 0.1 ± 0.01 | 0.1 ± 0.02 | 0.1 ± 0.01 | 0.1 ± 0.02 | 0.1 ± 0.01 | 0.1 ± 0.01 | DM |
B | 0.1 ± 0.01 | 0.2 ± 0.01 | 0.2 ± 0.02 | 0.2 ± 0.01 | 0.2 ± 0.01 | 0.2 ± 0.02 | DM |
Mo | 0.01 ± 0.00 | 0.02 ± 0.01 | 0.02 ± 0.00 | 0.02 ± 0.00 | 0.02 ± 0.01 | 0.02 ± 0.01 | DM |
EC | 0.96 ± 0.01 | 1.20 ± 0.01 | 1.45 ± 0.02 | 1.54 ± 0.01 | 2.05 ± 0.02 | 2.34 ± 0.06 | DM |
DM, demineralized water.
2.4. Plant material sampling
For each treatment and cycle, three plants were harvested every 7 days throughout the entire vegetation period. One plant was randomly selected from each highlighted sector (1–3), as shown in Figure 1. Additionally, a random plant from the edge (outside the sectors) was transferred to an empty space within each sector. The sampled plants were weighed fresh as a whole and then divided into leaves, stems, and flowers, which were weighed separately. Subsequently, the materials were dried at a constant moisture level of 8–10% at 25°C and reweighed. To determine the dry matter, a reference amount of each plant part was dried at 105°C until a constant weight was achieved. Just before analysis, the plant parts were homogenized. The flowers (including leaves until the 4th week) were frozen using liquid nitrogen and then ground using a mortar and pestle. The dried leaves (from the 5th week) and stems were ground using a grinder.
2.5. Dry decomposition and elemental analysis
The leaves, stems, and flowers were individually analyzed to determine the content of macroelements (excluding nitrogen), microelements, and trace elements in the plant samples. The plant biomass, which had been weighed and homogenized, was placed in a beaker and covered with a watch glass. The beaker was then placed on a hotplate set at 160°C, and the temperature was gradually increased to 350°C over a period of 4 hours to facilitate decomposition of the samples. Subsequently, the samples were transferred to a muffle furnace and maintained at a temperature of 450–500°C for 12 hours, as described in previous studies (Miholová et al., 1993). After cooling, 1 mL of 65% HNO3 was added to the beakers, which were then placed on a hot plate set at 120°C for 60 minutes. Following this step, the samples were annealed for 90 minutes in an oven at 500°C and suspended in 1.5% HNO3 with stirring using an ultrasonic bath. Elemental analysis of the samples was conducted using flame atomic absorption spectrometry (FAAS) on a Varian 280FS instrument (Varian, Australia), coupled with inductively coupled plasma optical emission spectrometry (ICP-OES) performed on a Varian Vista-PRO (Varian, Australia), (Mester and Sturgeon, 2003).
2.6. Determination of nitrogen content in plant material
The Kjeldahl method was employed to determine the total nitrogen content in the plant material, encompassing both organic and ammonia nitrogen. Approximately 0.50 g of the sample was weighed for analysis, followed by mineralization. The mineralization process was conducted in glass vials, where 2 g of catalyst (a mixture of 100 g K2SO4, 1 g CuSO4, 0.1 g Se) and 10 mL of concentrated sulfuric acid (H2SO4) were added to the sample. Decomposition took place for 90 minutes at a temperature of 420°C. After mineralization, the samples were prepared for distillation. Within the apparatus, 20 mL of distilled water was automatically added to the vial, followed by distillation into H3BO3, allowing for the determination of the total nitrogen content in the sample. The content was ascertained through titration using HCl (0.5 mol·L−1) and subsequently measured using the Gerhardt Vapodest 30s instrument (Königswinter, Germany) (Baker and Thompson, 1992).
2.7. Identification and quantification of cannabinoids
An optimized method of dynamic maceration was employed to identify and quantify cannabinoid content using HPLC/DAD (high-performance liquid chromatography with a diode-array detector), as described by Patel et al. (2017). Subsequently, 0.150 g of the sample was weighed into a 50 mL beaker, and 5 mL of solvent (96% ethanol) was added. The sample was then dynamically macerated for 60 minutes, filtered using a Morton filtration device (P16 porosity), and the filtrate was transferred to 50 mL vials. The same solvent in the same volume was added again to the initial 0.150 g of plant material, and the sample was dynamically macerated for another 60 minutes. This process was repeated three times, resulting in a composite sample derived from the three-phase dynamic maceration process in a 1:100 (V:W) ratio. This sample was diluted 20 times (using 96% ethanol), filtered through a syringe filter (0.22 μm), transferred to vials, and prepared for analysis in this form. Samples of the extracts were introduced into a high-performance liquid chromatography system equipped with diode array detection (HPLC-DAD; Agilent 1,260, Agilent Technologies Inc., United States), utilizing a Luna® 1C8 column (2) with dimensions of 250 × 3 mm and a particle size of 3 μm (Phenomenex, United States). The mobile phase employed was an isocratic mixture of acetonitrile/H2O (31:9, v/v) containing 0.1% HCOOH (v/v) and 0.1 mol/L NH4COOH (without pH adjustment). The flow rate was set at 0.55 mL/min, the temperature at 37°C, the sample injection volume at 8 μl, and UV detection was performed at 275 nm (Križman, 2020). To ensure accuracy, the instrument was externally calibrated using Δ9-tetrahydrocannabinol (Δ9-THC), Δ8-tetrahydrocannabinol (Δ8-THC), cannabigerolic acid (CBGA), cannabigerol (CBG) and cannabinolic acid (CBNA) ranging from 0.3 to 10 mg·L−1, and THCA ranging from 0.3 to 100 mg·L−1 (Sigma-Aldrich, Czech Republic) as reference standards. Data analysis was carried out using OpenLAB CDS software, ChemStation Edition, Rev. C.01.5.
2.8. Statistical analyses
The data were subjected to analysis of variance (one-way ANOVA) followed by Tukey’s Honestly Significant Difference (HSD) test using IBM SPSS Statistics software (version 25, 2017, Chicago, Illinois, United States).
3. Results
Implementing an enhanced nutritional regime and using different fertigation systems (1C and 2C) resulted in alterations in the plant nutrient composition of cannabis plants. Simultaneously, the utilization of two distinct systems resulted in variations in the pH levels of the nutrient solution. In the 2C system, a consistent pH value of 5.9 was upheld throughout the cycle, whereas in the 1C system, pH exhibited fluctuations within the range of 5.9–6.95. The nitrogenous compounds exhibited the lowest concentration in the stems, while the flowers displayed the highest concentration (Figure 2). Comparing the leaves and flowers of CN and ET plants under the augmented nutritional regime in the 1C cycle, statistically significant differences in N concentrations were observed only twice in the leaves (6th and 13th week) and twice in the stems (8th and 13th week). Notably, the most significant variation in N concentrations between CN and ET within the 1C system occurred in the stems during the 13th week, with a difference of 28% (CN: 10.44 mg·g−1, ET: 7.5 mg·g−1; Figure 2A). In contrast, N concentrations in the stems and leaves of CN and ET plants in 2C did not significantly differ in any of the tested weeks. Nevertheless, the concentration of N in flowers was significantly different in weeks 9 and 10 in the 2C regime. N concentrations exhibited over all the highest differences between 1C and 2C of ETs under the augmented nutritional regime, primarily noticeable in the flowers from week 8 to 11. During the last week, significant differences were observed only in the stems and leaves. The most significant variation in N concentrations between fertigation systems in ETs was observed during the 13th week, with a difference of 42% (21% between CNs) in the leaves (1C: 42.18 mg·g−1, 2C: 45.28 mg·g−1; Figure 2C).
The phosphorus (P) content exhibited lower levels in the leaves and stems compared to the flowers, displaying a cumulative pattern over time (Figure 3). In the stems, the P concentration of CN and ET in the 1C cycle significantly differed only in weeks 7 and 13. Significant differences were observed in the leaves in weeks 8, 9, and 13, while in the flowers, differences were found in weeks 8 and 10. Notably, the most significant variation in P concentration between nutritional treatments was observed in the leaves during the 13th week, with a difference of 63% (CN: 63.48 mg·g−1, ET: 38.9 mg·g−1; Figure 3A). In contrast to the 1C cycle, the P concentration of CN and ET in 2C plants showed significant differences only in the 11th week in the stems, with a difference of 108% (CN: 33.3 mg·g−1, ET: 69.15 mg·g−1; Figure 3B). Furthermore, the P concentration in the leaves and flowers of ET plants also varied between the 1C and 2C cycles. Significant differences were observed in the leaves in weeks 7, 9, 10, 11, and 12, while in the flowers, differences were found in weeks 7, 9, and 10. However, when comparing these two regimes, no statistically significant differences were observed in the P concentration of stems. The most notable differences in P concentrations between the 1C and 2C cycles in ETs were observed during the 11th week, with a difference of 31.6% (36% between CNs) in the leaves (1C: 64.87 mg·g−1, 2C: 49.31 mg·g−1; Figure 3C).
The potassium (K) content displayed the lowest levels in the stems and the highest levels in the leaves (Figure 4). In the leaves, the concentration of K for CN and ET with the 1C cycle started to exhibit significant differences from the 6th week until the 12th week. Significant differences were observed in the flowers in weeks 8, 10, and 11, while in the stems, differences were found only in week 9. Notably, the most significant variation in K concentration between CN and ET was observed in the stems during the 9th week, with a difference of 45.5% (CN: 16.12 mg·g−1, ET: 29.57 mg·g−1; Figure 4A). In contrast to the 1C cycle, the concentration of K in the 2C cycle of CN and ET plants showed significant differences only in the stems in weeks 10, 11, and 13. The most substantial differences in K concentrations between nutritional treatments were observed in the stems during the 10th week, with a difference of 77% (CN: 18.02 mg·g−1, ET: 31.88 mg·g−1; Figure 4B). Likewise, the K concentration also varied between the 1C and 2C cycles of ETs, with significant differences occurring in weeks 6, 7, 9, and 10. The most notable differences in K concentrations between 1C and 2C in ETs were observed in the stems during the 10th week, with a difference of 38.8% (20% between CNs) (1C: 23.03 mg·g−1, 2C: 31.88 mg·g−1; Figure 4C).
The iron (Fe) content exhibited the lowest levels in leaves and the highest levels in stems, displaying a cumulative trend over time (Figure 5). In the stems, the concentration of Fe for CN and ET with the 1C cycle showed significant differences in weeks 9, 11, and 13. Notably, the most significant variation in Fe concentration between CN and ET was observed in the stems during the 11th week, with a difference of 188.7% (CN: 337.2 mg·kg−1, ET: 955.8 mg·kg−1; Figure 5A). In contrast to the 1C cycle, the concentration of Fe in the stems of CN and ET plants in the 2C cycle exhibited significant differences only in week 10, in the leaves merely in weeks 6 and 10, and in the flowers from the 7th until the 8th week. The most substantial difference in Fe concentrations between CN and ET of 2C plants was observed in the flowers during the 8th week, with a difference of 55% (CN: 167.3 mg·kg−1, ET: 108.2 mg·kg−1; Figure 5B). Besides, the Fe concentration varied between the 1C and 2C cycles of ETs, with significant differences observed in the stems during weeks 6 and 11, in the leaves during weeks 11 and 12, and in the flowers during week 8. The most notable difference in Fe concentrations between 1C and 2C in ETs was observed in the stems during the 11th week, with a difference of 90.6% (11% between CNs) (1C: 337.2 mg·kg−1, 2C: 642.6 mg·kg−1; Figure 5C).
Implementing an augmented nutritional regime in two distinct nutritional cycles led to noticeable effects on the growth of cannabis plants. The biomass increase was relatively slow during the initial four weeks, but a significant acceleration was observed from week 5 onwards. Notably, the highest weekly dry weight gain was observed in the flowers (Figure 6). In the 1C cycle, the biomass increases in stems, leaves, and flowers for both CN and ET plants was nearly identical until week 4. However, from week 5 to 8, some differences in leaf and flower biomass emerged (Figure 6A), and in the 13th week, a statistically significant disparity in flower dry weight was noted. In contrast, the dry weight of stems, leaves, and flowers in the 2C cycle did not exhibit significant increases for CN and ET plants, except for the dry weight of leaves in the 3rd week and the dry weight of flowers in the 11th week. ET plants reached their maximum dry biomass at week 11, while CN plants achieved it by week 12 (Figure 6B). Biomass variations were also observed between the 1C and 2C cycles in ET plants, commencing from week 6 (Figure 6C).
The implementation of an augmented nutritional regime had statistically significant impact on the concentration of total THC (sum of Δ9-tetrahydrocannabinol, Δ8-tetrahydrocannabinol, and tetrahydrocannabinolic acid) total CBG (sum of cannabigerolic acid and cannabigerol), and CBNA in the flowers of medicinal cannabis plants, with distinct patterns observed for different fertigation systems and nutrition treatments (Figure 7). The total THC and CBNA concentration curves exhibited similar trends within the same growing cycle and treatment. At the same time, total THC yield (g/m2) was evaluated, where there was a clear trend in favor of ET across all comparisons.
Initially, the concentration of total THC in leaves and flowers almost stagnated in both treatments until week 4. However, from week 5 onwards, substantial growth with significant differences in total THC concentration was observed because only flowers were analyzed (Figures 7A-C). In the 1C cycle, significant differences in total THC concentration between CN and ET treatments occurred in weeks 5, 6, 9, and 10. The maximum concentrations were reached at week 11 for both CN (18.5%) and ET (19.4%) (Figure 7A). In contrast, in the 2C cycle, significant differences in total THC concentration between CN and ET treatments were observed only in the 6th week. CN reached its peak concentration at week 8 (15.9%), while ET reached its peak at week 9 (15.7%) (Figure 7B). Significant differences in total THC levels between 1C and 2C of ET treatments were observed from week 10 until the end of the experiment, where the 1C overperformed 2C (Figure 7C). The yield of this dominant cannabinoid (THC) showed significant differences in the 1C regime in favor of the ET treatment. Specifically, there was an increase in yield per square meter in weeks 5, 6 and 10 by 402%, 315% and 46%, respectively (Figure 8A). In drain to waste system (2C), this trend was similar. Still, a statistically significant difference was found only in the eleventh week, when there was an increase in THC yield in the ET variant by 50.7% (Figure 8B). An interesting finding is that when comparing the 1C and 2C regimes with the ET variant, there was an overall lower THC yield with the drain-to-waste system (2C). Specifically, there was a statistically significant difference in weeks 6, 10, 12 and 13 when there was a decrease in the THC yield of the 2C system to 6.6%, 25.3%, 41.5% and 31.1% of the THC yield value of the 1C recirculation system (Figure 8C).
For total CBG, the concentration in leaves and flowers nearly remained stagnant until week 4, with a notably lower concentration compared to THC by more than an order of magnitude. However, a significant increase in total CBG concentration was observed in week 5, focusing only on the flowers (Figures 7D-F). In the 1C cycle, the concentrations of total CBG in CN and ET treatments began to differ significantly from week 7 until week 11 and peaked at week 10 for CN (1.1%) and week 6 for ET (0.9%). The most significant differentiation between CN and ET occurred at week 10, with a difference of 218% in favor of CN. Following the peak concentrations, a foretellable decrease in concentration was observed, followed by a slight increase in the last three weeks. In the 2C cycle, significant differences in total CBG concentration between CN and ET treatments were observed only in weeks 6, 8, and 11. The maximum concentrations were reached at week 9 for both CN (1.3%) and ET (0.8%) (Figure 7E). Significant differences in total CBG levels between 1C and 2C of ET treatments were observed in weeks 1, 2, and from week 6 except the 10th week until the end of the experiment (Figure 7F).
Regarding CBNA concentration, significant variations between CN and ET treatments in the 1C cycle were observed in weeks 5 and 6. CBNA concentration peaked at week 10 for CN (0.06%) and week 13 for ET (0.06%) (Figure 7G). In the 2C cycle, significant differences in CBNA concentration between CN and ET treatments were observed in weeks 6, 7, and 11. CN reached its maximum concentration (0.04%) in 2C at week 11, with a significant difference of 55%, while ET reached its maximum in the 13th week (Figure 7H). CBNA concentrations between 1C and 2C of ET treatments showed significant differences in the first two weeks, followed by nearly identical concentration levels in weeks 3 and 4. From week 5 onwards, a steady growth trend was observed, with 1C outperforming 2C. Significant differences were observed in weeks 1, 2, 5, and 6 and from week 10 until the end of the experiment. As mentioned, the CBNA concentrations in both ET variants reached their maximum at week 13, with a significant difference of 34% (Figure 7I).
4. Discussion
The present study focused on investigating the effects of an enhanced nutritional regime and different fertigation systems on the ionome composition, biomass, and cannabinoid content of indoor-cultivated cannabis plants. The results indicate substantial variations in N, P, K, and Fe concentrations in different plant parts under the augmented nutritional regime. These alterations could be attributed to nutrient uptake, translocation, and allocation modifications resulting from the introduced nutritional variations. Undoubtedly, nutrition plays a pivotal role in shaping the development, functionality, and metabolic processes of various plant organs and tissues (Shiponi and Bernstein, 2021b; De Prato et al., 2022). Existing knowledge has extensively documented the ideal thresholds of specific macronutrients, including N, P and K, along with micronutrients like Fe, necessary for ensuring the normal growth and functioning of both root systems and above-ground biomass (Shiponi and Bernstein, 2021a; Shiponi and Bernstein, 2021b; Malík et al., 2022; Saloner and Bernstein, 2022; Llewellyn et al., 2023), as well as the production of valuable secondary metabolites in medicinal cannabis plants (Caplan et al., 2017; Bernstein et al., 2019b; Shiponi and Bernstein, 2021b; Saloner and Bernstein, 2022). It is crucial to emphasize that nutrient solution fertigation was provided to the plants only up to the ninth week of cultivation. Subsequently, from the ninth week onwards, the plants received fertigation solely with deionized (DM) water. Consequently, the plants relied solely on their accumulated nutrient reserves during this period. Modifications in the nutrient solution’s P, K, and Fe concentrations within two distinct nutrient cycles undeniably exerted discernible impacts on the levels of macro and micro-elements in cannabis plants’ above-ground organs. While the N concentration remained unchanged between the CN and ET systems, some statistically significant differences in N distribution among the above-ground organs of cannabis plants were observed within each fertigation system. The composition and proportional disbalance of nutrients in nutrient solutions have been recognized as factors capable of influencing the profiles of cannabinoids and terpenes in cannabis. These observed variations are likely linked to the complex interplay between nutrient levels in the solution and the growth of the plant. The different fertigation systems employed in this study also exhibited discernible effects on the concentrations of both macro- and microelements within the tissues of cannabis plants. It should be emphasized that the interactions between cations and anions of nutrients during root cell membrane transport have been extensively documented in the scientific literature (Meychik and Yermakov, 2001; Roberts, 2006; Langenfeld et al., 2022).
Nevertheless, variations in the pH of the nutrient solution, such as those observed in the RS (1C), and the increased nutrient quantity provided in the DS (2C), have the potential to influence the availability of specific nutrients and, consequently, the physiological and metabolic reactions within the plants (Malík et al., 2023). Furthermore, a more pronounced variation in N concentrations was evident when comparing the two distinct fertigation systems, which aligns with our study. Nitrogen (N) plays a vital role in organs with high metabolic activity, such as leaves and flowers, where it is a crucial component of proteins and compounds related to energy metabolism. As a result, N tends to accumulate in smaller quantities in structural compartments with lower metabolic activity, such as stems (Marschner and Marschner, 2012). In our investigation comparing the 2C and 1C systems, the 2C system consistently exhibited higher N concentrations in flowers during weeks 8 to 11 (an average 8.7% increase in favor of the 2C system, with values ranging from 6.3% to 11.6%). This contrasts with the findings of Malík et al. (2023), where they observed higher levels of N in plant tissues in the recirculation systems. In our case, we presume that the more frequent nutrient replenishment (2–3 times a week) in the drain-to-waste system resulted in more significant nutrient accumulation in both the above-ground biomass and the roots. Subsequently, when nutrition was stopped, and only DM water was irrigated, there was translocation from the roots to the above-ground biomass. Also, another hypothesis could be that some of the nutrients that were supplied in addition to the ET variant (P, K, or Fe) could be limiting in the recirculation system. Therefore, when it was supplied more in drain-to-waste, other nutrients (in this case, N) could subsequently be more incorporated into tissues (Ågren et al., 2012; Asao, 2012).
As expected, alterations in the concentrations of P, K, and Fe in the nutrient solution in some observations significantly impacted the concentrations of these elements in the above-ground organs of medicinal cannabis plants. Notably, in the 1C regime, differences between CN and ET were more pronounced than in the 2C regimen. This observation can be attributed to the fluctuation and pH increase in the 1C compared with the stable pH in the 2C regime (Malík et al., 2023). Surprisingly, despite doubling the concentration of P in the nutrient solution in ET compared to CN within the 1C fertigation system, P accumulation was notably higher in CN than in ET. This observation might be attributed to the fact that root analysis was not conducted; therefore, it remains plausible that any excess P is encased in the roots. As indicated by the accumulation data, P content within the plant tissues did not exhibit a decline following the cessation of nutrient supplementation. This phenomenon suggests a potential release of P from vacuoles (Shiponi and Bernstein, 2021a; Shiponi and Bernstein, 2021b), serving as a defense mechanism against P toxicity in the shoot (Hawkesford et al., 2012). This mechanism prevents the exposure of shoot cells to damaging P concentrations by compartmentalizing excessive P supply. The compartmentalization of P is recognized as a fundamental mechanism for averting the accumulation of cytoplasmatic P to toxic levels. Under adequate P nutrition, approximately 70–95% of intercellular phosphate (Pi) is sequestered within vacuoles. Consequently, the regulation of transporters under varying P conditions facilitates the maintenance of Pi cellular homeostasis (Liu et al., 2016). So far, there is no information about P transporters in Cannabis sativa L. This information is needed to understand better P remobilization and translocation in the cannabis plant (Shiponi and Bernstein, 2021b). An alternative explanation for this unexpected result could be a potential nutrient lockout effect in the ET due to the elevated P concentration in the solution. Phosphorus has an inclination to form complexes with soil minerals, such as iron Fe, potentially diminishing its overall bioavailability. However, it is essential to note that such an effect was observed by Mardamootoo et al. (2021) exclusively in soil-based cultivation systems. When soluble phosphatic fertilizers are applied to soils, they initially dissolve, causing an immediate increase in soil solution P concentration. Subsequently, P primarily engages in adsorption and precipitation processes (Prasad and Power, 1997). The reactions that transpire among phosphate ions in the soil solution, soil constituents, and non-phosphatic components in the fertilizers primarily sequester P from the solution phase, rendering phosphate less soluble over time (Sample et al., 2015). This phenomenon is commonly referred to as P fixation, adsorption, or retention. Consequently, P becomes markedly immobile in soils and tends to remain in proximity to the point of application (Prasad and Power, 1997; Mardamootoo et al., 2021). In plant cells, P serves as a critical component of nucleic acids, membrane lipids, and phosphorylated intermediates involved in energy metabolism. As a result, maintaining cellular phosphorus homeostasis is essential for ensuring the normal function of various physiological and biochemical processes (Raghothama, 2015). The elevated P concentration observed in the 2C system compared to the 1C system may likely be attributed to an increase and fluctuations in pH within the 1C regime, a phenomenon in line with the findings of Lefever et al. (2017). The higher pH levels, reaching up to 6.95 in the nutrient solution, could lead to the partial precipitation of phosphates by calcium (Ca2+) and magnesium (Mg2+). This precipitation could result in the formation of insoluble and consequently unavailable salts within the hydroponics solution (Lee et al., 2017).
K exhibits a clear accumulation advantage in the ET variant in 1C as a highly bioavailable plant element. In accordance with (Bernstein et al., 2019a), it is noteworthy that cannabis plant stems tend to exhibit a preference for K accumulation, as evidenced by concentrations similar to those found in fan leaves. Typically, this nutrient’s temporary storage occurs in the stem’s xylem or phloem parenchyma. The K concentration in cannabis stems remained notably high at the end of the developmental period, a phase characterized by high nutrient demand, indicating active accumulation. While our study did not measure sugar leaves but entire inflorescences, where K accumulation was slightly higher than in the stem, the results for K suggest that the nutrient concentration does not significantly decrease in tissues after nutrient cessation, indicating its excessive accumulation. This phenomenon, known as the luxury consumption of K and its temporary storage, has been observed in various plant species (Marschner and Marschner, 2012). It is pertinent to add that, based on the experiments conducted by Saloner et al. (2019) and Saloner and Bernstein (2022), there appears to be no competition for nutrient uptake between K and the uptake of N and P at suitable concentrations. Yep and Zheng (2021) studied the yield of Cannabis sativa inflorescence with K fertilizer in aquaponic solutions, similar to hydroponic systems. The added K and micronutrient fertilizers did not affect vegetative growth or leaf physiology. However, a positive linear correlation was observed between K concentration in the nutrient solution and both apical inflorescence yield (g/plant) and total inflorescence yield. Additionally, K fertilizer enhanced the harvest index.
Similarly, Fe demonstrates higher accumulation in the ET in 1C, except for the end of the cultivation, where the trend reverses. In contrast, interestingly, within the 2C fertigation system, significant differences in element accumulation between CN and ET treatments for any of the elements were absent. Generally, the differences in nutrient accumulation were predominantly insignificant when comparing CN and ET. These distinctions can be considered negligible, except for a few isolated instances. Furthermore, the disparities between the 1C and 2C regimens can be attributed to the 2C system’s superior control over nutrient solutions, ensuring stable and precise nutrient concentrations (Malík et al., 2023). Analyzing these variations in nutrient accumulation between the 1C and 2C fertigation systems provides valuable insights into how fertigation systems shape the ‘ionome’ composition of cannabis plants. We found that N concentrations favored the 2C system during later flowering stages, notably in weeks 8 through 11. P concentrations in flowers also showed advantages for the 2C system, particularly in weeks 7, 9, and 10. Conversely, P concentrations in leaves favored the 1C regimen, especially in weeks 9 through 12. K concentrations in leaves followed a similar trend, favoring the 1C system in weeks 6, 7, and 9. Fe concentrations exhibited variations favoring both systems, with Fe in flowers favoring the 1C system at week 8, whereas Fe in leaves and stems favored the 1C system in weeks 11, 12, and 6, respectively. These findings underscore the complexity of nutrient dynamics in fertigation, with the 2C regimen’s precise control over nutrient concentrations offering valuable insights for optimizing nutrient uptake and enhancing plant development. Considering the results regarding all concentrations of the studied elements in all the plant organs examined, our third hypothesis concerning the influence of the changed nutritional regime on the composition of the ‘ionome’ was confirmed.
Regarding the dry weight of above-ground organs, the ET had minimal effects compared to the CN. When comparing CN with ET in the 1C fertigation system, statistically significant changes were primarily observed in the dry weight of leaves. The supplementation of P did not induce any significant alterations in cannabinoid concentration within the flowers, but it led to a decrease in concentration within the inflorescence leaves. It is noteworthy that studies conducted with industrial hemp have demonstrated varying responses to P fertilization depending on growing conditions and cultivars (Bernstein et al., 2019b). In instances where the initial P concentration in the soil is relatively high, the additional P application tends to have limited effects. Phosphorus is indeed a critical nutrient required in relatively substantial amounts by plants. However, conventional practices in the medicinal cannabis industry involve the application of higher P concentrations than those typically used for most other crops. This practice is rooted in the belief that cannabis plants necessitate elevated P concentrations to optimize their functionality and enhance yield (Bevan et al., 2021). It is worth noting that mineral nutrition significantly influences plant morpho-development, and our findings indicate that supplemental P beyond optimal requirements does not contribute to further increases in plant biomass. This outcome aligns with prior research on hemp (Vera et al., 2010). It reinforces our conclusion that medicinal cannabis exhibits a broad optimal P-level range (Shiponi and Bernstein, 2021a). Interestingly, in the ET, plants exhibited greater leaf biomass during the early generative phase of flowering. However, their leaf biomass declined as flower maturation progressed and became statistically more minuscule than the CN. These findings intriguingly suggest that ET group plants in the 1C system displayed a more foliated profile at the onset of flowering, a phenomenon in line with the work of Hawkesford et al. (2012). Shiponi and Bernstein (2021a) recommend a minimum P requirement of 15 mg·L−1 P, with a suggested application rate of 30 mg·L−1, based on the functional physiology and ionome profiling revealing genotypic variability in P sensitivity. Notably, our tested P application rates, ranging from a minimum of 32 mg·L−1 to a maximum of 93 mg·L−1, significantly exceed those employed in the reference study. Nevertheless, fewer leaves were present in the inflorescence during the last two weeks of cultivation in the 1C system. This phenomenon could simplify post-harvest inflorescence arrangement and contribute to a higher quality of harvested inflorescence. In conclusion, when comparing CN and ET at 1C, it becomes evident that ET produces noticeably more biomass in flowers only at week 8, potentially suggesting this is an optimal harvesting time. However, in the subsequent weeks, flower biomass continued to increase. As we will explore in the following discussion, cannabinoid concentrations, mainly THC, were also found to be higher in the later stages of cultivation. Therefore, a later harvest time beyond week 8 seems preferable to optimize economic efficiency while maximizing yield. When comparing CN with ET in the 2C fertigation system, statistically significant differences were observed only during the 11th week of cultivation. Comparing the two fertigation systems, 1C and 2C, statistically significant differences in the dry weight of leaves and flowers favored 1C at weeks 6 and 7 and favored 1C only for flowers at week 10. However, no statistically significant differences were observed at the end of cultivation. In summary, the results of the dry weight yield of biomass suggest that neither the altered nutritional regime of the ET variant nor the type of fertigation system had a statistically conclusive effect on the final inflorescence harvest mass of medicinal cannabis. Therefore, our first hypothesis about the inflorescence yield was only partially confirmed because we recorded a statistically significant difference between CN and ET in the harvest period only at week 11 in the 2C system.
Cannabinoid levels represent a critical aspect in assessing the excellence of medicinal cannabis. Subjects purchasing the resulting cultivated product – the dry female inflorescence of medicinal cannabis – when selecting a supplier are oriented largely according to THC concentrations and, eventually, CBD in dry matter (Zhu et al., 2021). Since cannabinoids are the constituents conferring exceptional value to cannabis, it is imperative to conduct additional research into the influence of mineral nourishment on cannabis productivity and the association between yield and potency (Bevan et al., 2021). Prior research suggests that plants’ mineral nutrition can influence the production of secondary metabolites in cannabis (Caplan et al., 2017; Saloner and Bernstein, 2021; Saloner and Bernstein, 2022). An inverse correlation between cannabis yield and cannabinoid concentrations has been observed in previous studies. These studies have consistently reported a linear decrease in cannabinoid concentrations with increasing yield (the dilution effect) (Caplan et al., 2017; Yep et al., 2020; Shiponi and Bernstein, 2021a). This effect was not observed to a significant extent in our case. Additionally, future investigations should encompass the examination of other compounds acknowledged for their influence on product quality. Treatments involving metals such as Fe and Cu have demonstrated the potential to enhance secondary metabolite production in numerous plant species (Gorelick and Bernstein, 2014). The impact of metals, including Ni, Ag, Fe, and Co, on bioactive compound production is attributed to their ability to modulate various aspects of secondary metabolism Zhao et al., 2001). In the context of hemp cultivation, it has been observed that N supplementation leads to increased plant height and biomass (Papastylianou et al., 2018). Conversely, it is noteworthy that P or K fertigation treatments have shown limited efficacy in eliciting substantial responses (Aubin et al., 2015). It is important to acknowledge that while these findings provide valuable insights into plant growth dynamics, their direct relevance to cannabis cultivation is somewhat limited. In the context of cannabis, the paramount consideration is the concentration of therapeutic cannabinoids, a factor that far outweighs concerns related to overall biomass.
This study correspondingly delved into the intricate relationship between fertigation systems and cannabinoid concentrations in medicinal cannabis, providing insights contributing to the growing body of knowledge in this field. Our findings elucidate critical factors influencing both yield and cannabinoid potency and yield.
These findings, confirming our 4th hypothesis, align with prior research, particularly studies exploring the impact of fertigation systems on cannabis cultivation. For example, a study comparing recirculation and drain-to-waste systems found that recirculation led to higher yields of THCA and CBNA, the prominent cannabinoids in medical cannabis chemotype I, but resulted in lower sesquiterpene concentrations. Drain-to-waste, however, allowed for better control over nutrient delivery but consumed more resources and yielded fewer monoterpenes and THCA (Malík et al., 2023). In our case, there was also a significant increase in THC yield in the recirculation system (up to 182% more in the 10th week).
A previous investigation (Bernstein et al., 2019b) noted that employing mild nutritional treatments, which closely adhered to the optimal range for plant growth, resulted in subtle changes in plant development. However, there emerged an observable influence on the cannabinoid profile. These findings suggest the intriguing possibility that slight adjustments in nutritional status could play a role in modulating secondary metabolism in cannabis. This statement is also confirmed by our results, where the combination of changes in the yield and concentrations of the dominant cannabinoid (THC) resulted in an increase in the yield of total THC (g·m2) by an average of 48% during harvest period for the variant with enriched nutrition.
In accordance with prior research, our findings align with established knowledge, indicating that the lowest recommended P supply for optimal yield is 30 mg·L−1, and optimal yields are maintained at concentrations up to 90 mg·L−1 (Shiponi and Bernstein, 2021b). An earlier study also observed that the impact of increased P supplementation on medicinal cannabis is dependent upon the specific organ and compound. Specifically, the concentrations of key cannabinoids such as THC, CBD, CBN, and CBG remained unaffected by the treatment of enhanced P supplementation. The addition of P did not significantly affect cannabinoid concentrations within the flowers; however, it did lead to a reduction in cannabinoid concentration in the leaves of the inflorescence (Bernstein et al., 2019b). It is worth noting that in our study, the experimental ET variants were exposed to elevated P levels, reaching as high as 93 mg·L−1 P toward the end of cultivation. Importantly, similar to existing research, we did not observe any significant impact of elevated P levels in the nutrient solution on the production of cannabinoids. In contrast, interestingly, an older study from Coffman and Gentner (1977) reported a notable increase in total cannabinoid content per plant with higher P supply. This is consistent with our finding that the nutrient-enhanced variant achieved higher THC yields.
In line with previous research, our findings confirm that increasing K supply beyond 60 mg·L−1 K does not yield beneficial effects. It is noteworthy that variations in K supply had a relatively modest impact on both cannabinoid and terpenoid levels within the plant (Saloner and Bernstein, 2022). Notably, the decrease observed in these compounds was generally mild and less pronounced compared to the significant influence of N on secondary metabolism (Saloner and Bernstein, 2021). It is worth noting that in our study, the ET variant had elevated K levels up to 265 mg·L−1. Similarly, we did not observe any significant impact of increased K levels on the augmentation of cannabinoid concentrations. However, it is necessary to point out again that there was a significant increase in THC yield per square meter in ET variant.
To the best of our knowledge, no prior studies have explored the influence of Fe concentration in the nutrient solution on the content of secondary metabolites in the harvested inflorescence of medicinal cannabis. Our investigation revealed a notable outcome, which examined the effects of elevated macro- and micronutrient levels, including iron, in the ET variant. Contrary to some expectations, the increased iron content did not substantially impact the final cannabinoid content of the harvested medicinal cannabis inflorescence.
Exploring additional variables, such as specific nutrient formulations and variations in environmental conditions, will be essential to comprehensively understand their influence on cannabinoid production in indoor-cultivated cannabis plants.
The complexity of the interactions between various nutrients, environmental factors, and cannabinoid biosynthesis remains a rich area for investigation. Refining our understanding of these dynamics can unlock even more precise methods for optimizing cannabinoid production in indoor-cultivated cannabis. Such advancements hold great promise for both enhancing the therapeutic potential of indoor-cultivated cannabis products and improving their overall cultivation efficiency. This knowledge is valuable for cultivators seeking to maximize both yield and cannabinoid potency in a sustainable and resource-efficient manner.
5. Conclusion
In this study, we explored the intricate relationship between enhanced nutrition, different fertigation systems, and the ionome composition, biomass, and cannabinoid content of indoor-cultivated cannabis plants, comparing the yield of THC in the harvested inflorescence per square meter. Our research reveals multifaceted connections between nutrient supply, plant development, and valuable secondary metabolites production in medicinal cannabis. Significant nutrient variations, including N, P, K, and Fe concentrations in various plant organs under the augmented nutritional regime, result from shifts in nutrient uptake, translocation, and allocation due to changes in nutrient supply. This underscores the pivotal role of nutrition in shaping plant development, functionality, and metabolic processes, aligning with prior studies. Our findings emphasize the need for precise nutrient management strategies in cannabis cultivation. Adjusting nutrient regimens in different fertigation systems yielded significant outcomes. Notably, increased iron content in the ET variant did not substantially impact cannabinoid content, challenging assumptions. The impact of different fertigation systems on nutrient concentrations within plant tissues was evident. Adjusting the nutrient regimen in the 2C system led to a significant 108% increase in P accumulation in the stems in the 11th week, and there was a preference for K accumulation in cannabis stems, with up to a 77% increase in K levels in the ET. The ET increased the Fe concentrations, mainly in stems, up to 189%. The yield of THC in flowers per square meter with enriched nutrition increased by up to 50.7% compared to the control variant. There was even an 182% increase in THC yield in the recirculation system in the case of the ET comparison. As we conclude this study, we invite further investigations into the nuanced relationship between mineral nutrition and cannabis productivity, considering the complex interplay between yield and potency. Our research resonates beyond the scientific realm, offering actionable guidance for cultivators striving to maximize both yield and cannabinoid potency. While our study has yielded valuable insights into the relationship between fertigation systems and cannabinoid concentrations in indoor-cultivated cannabis, it is essential to acknowledge that further investigations should specifically delve deeper into the impact of enhanced nutritional compositions within the nutrient solution.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
JV: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – review & editing. MM: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – review & editing. JŠ: Data curation, Formal analysis, Validation, Visualization, Writing – original draft. PT: Project administration, Resources, Supervision, Validation, Writing – review & editing.
Funding Statement
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by METROFOOD-CZ research infrastructure project (MEYS Grant No: LM2023064) including access to its facilities. This study was also supported by the European Regional Development Fund (project NUTRISK no. CZ.02.1.01/0.0/0.0/16_019/0000845).
Conflict of interest
The 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.2024.1322824/full#supplementary-material
References
- Ågren G. I., Wetterstedt J.Å.M., Billberger M. F. K. (2012). Nutrient limitation on terrestrial plant growth – modeling the interaction between nitrogen and phosphorus. New Phytol. 194, 953–960. doi: 10.1111/j.1469-8137.2012.04116.x [PubMed] [CrossRef] [Google Scholar]
- Asao T. (2012). Hydroponics: a standard methodology for plant biological researches (Rijeka, Croatia: InTech; ). [Google Scholar]
- Atherton H. R., Li P. (2023). Hydroponic cultivation of medicinal plants–plant organs and hydroponic systems: techniques and trends. Horticulturae 9, 349. doi: 10.3390/horticulturae9030349 [CrossRef] [Google Scholar]
- Aubin M., Seguin P., Vanasse A., Tremblay G. F., Mustafa A. F., Charron J. (2015). Industrial hemp response to nitrogen, phosphorus, and potassium fertilization. Crop Forage Turfgrass Mgmt 1, 1–10. doi: 10.2134/cftm2015.0159 [CrossRef] [Google Scholar]
- Baker W. H., Thompson T. L. (1992). Determination of total nitrogen in plant samples by Kjeldahl. Plant Anal. reference procedures South. region United States 1992, 13–16. [Google Scholar]
- Balliu A., Zheng Y., Sallaku G., Fernández J. A., Gruda N. S., Tuzel Y. (2021). Environmental and cultivation factors affect the morphology, architecture and performance of root systems in soilless grown plants. Horticulturae 7, 243. doi: 10.3390/horticulturae7080243 [CrossRef] [Google Scholar]
- Bernstein N., Gorelick J., Koch S. (2019. a). Interplay between chemistry and morphology in medical cannabis (Cannabis sativa L.). Ind. Crops Products 129, 185–194. doi: 10.1016/j.indcrop.2018.11.039 [CrossRef] [Google Scholar]
- Bernstein N., Gorelick J., Zerahia R., Koch S. (2019. b). Impact of N, P, K, and humic acid supplementation on the chemical profile of medical cannabis (Cannabis sativa L). Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00736 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Bevan L., Jones M., Zheng Y. (2021). Optimisation of nitrogen, phosphorus, and potassium for soilless production of cannabis sativa in the flowering stage using response surface analysis. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.764103 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Bouchard M., Dion C. B. (2009). Growers and facilitators: probing the role of entrepreneurs in the development of the cannabis cultivation industry. J. Small Business Entrepreneurship 22, 25–37. doi: 10.1080/08276331.2009.10593440 [CrossRef] [Google Scholar]
- Briat J.-F., Dubos C., Gaymard F. (2015). Iron nutrition, biomass production, and plant product quality. Trends Plant Sci. 20, 33–40. doi: 10.1016/j.tplants.2014.07.005 [PubMed] [CrossRef] [Google Scholar]
- Caplan D., Dixon M., Zheng Y. (2017). Optimal rate of organic fertilizer during the flowering stage for cannabis grown in two coir-based substrates. horts 52, 1796–1803. doi: 10.21273/HORTSCI12401-17 [CrossRef] [Google Scholar]
- Chandra S., Lata H., ElSohly M. A., Walker L. A., Potter D. (2017). Cannabis cultivation: Methodological issues for obtaining medical-grade product. Epilepsy Behav. 70, 302–312. doi: 10.1016/j.yebeh.2016.11.029 [PubMed] [CrossRef] [Google Scholar]
- Chaouqi S., Moratalla-López N., Alonso G. L., Lorenzo C., Zouahri A., Asserar N., et al. (2023). Effect of Soil Composition on Secondary Metabolites of Moroccan Saffron (Crocus sativus L.). Plants 12, 711. doi: 10.3390/plants12040711 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Chrysargyris A., Drouza C., Tzortzakis N. (2017). Optimization of potassium fertilization/nutrition for growth, physiological development, essential oil composition and antioxidant activity of Lavandula angustifolia Mill. J. Soil Sci. Plant Nutr. 291–306. doi: 10.4067/S0718-95162017005000023 [CrossRef] [Google Scholar]
- Coffman C. B., Gentner W. A. (1977). Response of greenhouse grown Cannabis sativa L. to nitrogen, phosphorus and potassium. Agron. J. 69 (5), 832–836. doi: 10.2134/agronj1977.00021962006900050026x [CrossRef] [Google Scholar]
- Danziger N., Bernstein N. (2021). Plant architecture manipulation increases cannabinoid standardization in ‘drug-type’ medical cannabis. Ind. Crops Products 167, 113528. doi: 10.1016/j.indcrop.2021.113528 [CrossRef] [Google Scholar]
- De Prato L., Ansari O., Hardy G. E., St. J., Howieson J., O’Hara G., et al.. (2022). The cannabinoid profile and growth of hemp (Cannabis sativa L.) is influenced by tropical daylengths and temperatures, genotype and nitrogen nutrition. Ind. Crops Products 178, 114605. doi: 10.1016/j.indcrop.2022.114605 [CrossRef] [Google Scholar]
- Gaaliche B., Ladhari A., Zarrelli A., Ben Mimoun M. (2019). Impact of foliar potassium fertilization on biochemical composition and antioxidant activity of fig (Ficus carica L.). Scientia Hortic. 253, 111–119. doi: 10.1016/j.scienta.2019.04.024 [CrossRef] [Google Scholar]
- Gorelick J., Bernstein N. (2014). “Chapter five – elicitation: an underutilized tool in the development of medicinal plants as a source of therapeutic secondary metabolites,” in Advances in agronomy (Amsterdam: Elsevier; ), 201–230. doi: 10.1016/B978-0-12-800138-7.00005-X [CrossRef] [Google Scholar]
- Hawkesford M., Horst W., Kichey T., Lambers H., Schjoerring J., Møller I. S., et al.. (2012). “Functions of macronutrients,” in Marschner’s mineral nutrition of higher plants (Amsterdam: Elsevier; ), 135–189. doi: 10.1016/B978-0-12-384905-2.00006-6 [CrossRef] [Google Scholar]
- Hesami M., Baiton A., Alizadeh M., Pepe M., Torkamaneh D., Jones A. M. P. (2021). Advances and perspectives in tissue culture and genetic engineering of cannabis. IJMS 22, 5671. doi: 10.3390/ijms22115671 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Janatová A., Fraňková A., Tlustoš P., Hamouz K., Božik M., Klouček P. (2018). Yield and cannabinoids contents in different cannabis (Cannabis sativa L.) genotypes for medical use. Ind. Crops Products 112, 363–367. doi: 10.1016/j.indcrop.2017.12.006 [CrossRef] [Google Scholar]
- Johnson R., Vishwakarma K., Hossen Md. S., Kumar V., Shackira A. M., Puthur J. T., et al.. (2022). Potassium in plants: Growth regulation, signaling, and environmental stress tolerance. Plant Physiol. Biochem. 172, 56–69. doi: 10.1016/j.plaphy.2022.01.001 [PubMed] [CrossRef] [Google Scholar]
- Král D., Šenkyřík J. B., Ondřej V. (2022). Expression of genes involved in ABA and auxin metabolism and LEA gene during embryogenesis in hemp. Plants 11, 2995. doi: 10.3390/plants11212995 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Križman M. (2020). A simplified approach for isocratic HPLC analysis of cannabinoids by fine tuning chromatographic selectivity. Eur. Food Res. Technol. 246, 315–322. doi: 10.1007/s00217-019-03344-7 [CrossRef] [Google Scholar]
- Kumar S., Kumar S., Mohapatra T. (2021). Interaction between macro- and micro-nutrients in plants. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.665583 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Langenfeld N. J., Pinto D. F., Faust J. E., Heins R., Bugbee B. (2022). Principles of nutrient and water management for indoor agriculture. Sustainability 14, 10204. doi: 10.3390/su141610204 [CrossRef] [Google Scholar]
- Lata H., Chandra S., Techen N., Khan I. A., ElSohly M. A. (2016). In vitro mass propagation of Cannabis sativa L.: A protocol refinement using novel aromatic cytokinin meta-topolin and the assessment of eco-physiological, biochemical and genetic fidelity of micropropagated plants. J. Appl. Res. Medicinal Aromatic Plants 3, 18–26. doi: 10.1016/j.jarmap.2015.12.001 [CrossRef] [Google Scholar]
- Lee J. Y., Rahman A., Azam H., Kim H. S., Kwon M. J. (2017). Characterizing nutrient uptake kinetics for efficient crop production during Solanum lycopersicum var. cerasiforme Alef. growth in a closed indoor hydroponic system. PloS One 12, e0177041. doi: 10.1371/journal.pone.0177041 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Lefever K., Laubscher C. P., Ndakidemi P. A., Nchu F. (2017). “Effects of pH and Phosphorus Concentrations on the Chlorophyll Responses of Salvia chamelaeagnea (Lamiaceae) Grown in Hydroponics,” in Chlorophyll. Eds. Jacob-Lopes E., Zepka L. Q., Queiroz M. I. (Rijeka: InTech; ). doi: 10.5772/67610 [CrossRef] [Google Scholar]
- Liu T.-Y., Huang T.-K., Yang S.-Y., Hong Y.-T., Huang S.-M., Wang F.-N., et al.. (2016). Identification of plant vacuolar transporters mediating phosphate storage. Nat. Commun. 7, 11095. doi: 10.1038/ncomms11095 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Llewellyn D., Golem S., Jones A. M. P., Zheng Y. (2023). Foliar symptomology, nutrient content, yield, and secondary metabolite variability of cannabis grown hydroponically with different single-element nutrient deficiencies. Plants 12, 422. doi: 10.3390/plants12030422 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Malík M., Praus L., Tlustoš P. (2023). Comparison of recirculation and drain-to-waste hydroponic systems in relation to medical cannabis (Cannabis sativa L.) plants. Ind. Crops Products 202, 117059. doi: 10.1016/j.indcrop.2023.117059 [CrossRef] [Google Scholar]
- Malík M., Velechovský J., Praus L., Janatová A., Kahánková Z., Klouček P., et al.. (2022). Amino acid supplementation as a biostimulant in medical cannabis (Cannabis sativa L.) plant nutrition. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.868350 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Malík M., Velechovský J., Tlustoš P. (2021). The overview of existing knowledge on medical cannabis plants growing. Plant Soil Environ. 67, 425–442. doi: 10.17221/96/2021-PSE [CrossRef] [Google Scholar]
- Mardamootoo T., Du Preez. C.C., Barnard. J. H. (2021). Phosphorus management issues for crop production: A review. Afr. J. Agric. Res. 17, 939–952. doi: 10.5897/AJAR2020.15205 [CrossRef] [Google Scholar]
- Marschner H., Marschner P. (2012). Marschner’s mineral nutrition of higher plants. 3rd ed (Amsterdam: Academic Press; ). [Google Scholar]
- Mester Z., Sturgeon R. (2003). Sample preparation for trace element analysis. 1st ed (Amsterdam: Elsevier; ). [Google Scholar]
- Meychik N. R., Yermakov I. P. (2001). Ion exchange properties of plant root cell walls. Plant Soil 234, 181–193. doi: 10.1023/A:1017936318435 [CrossRef] [Google Scholar]
- Miholová D., Mader P., Száková J., Slámová A., Svatoš Z. (1993). CzechoSlovakian biological certified reference materials and their use in the analytical quality assurance system in a trace element laboratory. Fresenius J. Anal. Chem. 345, 256–260. doi: 10.1007/BF00322606 [CrossRef] [Google Scholar]
- Moeller K., Lindholst C. (2014). Sanctioning large-scale domestic cannabis production – potency, yield, and professionalism. Scandinavian J. Forensic Sci. 20, 20–27. doi: 10.2478/sjfs-2014-0002 [CrossRef] [Google Scholar]
- Naumann M., Koch M., Thiel H., Gransee A., Pawelzik E. (2020). The importance of nutrient management for potato production part II: plant nutrition and tuber quality. Potato Res. 63, 121–137. doi: 10.1007/s11540-019-09430-3 [CrossRef] [Google Scholar]
- Nguyen P. M., Kwee E. M., Niemeyer E. D. (2010). Potassium rate alters the antioxidant capacity and phenolic concentration of basil (Ocimum basilicum L.) leaves. Food Chem. 123, 1235–1241. doi: 10.1016/j.foodchem.2010.05.092 [CrossRef] [Google Scholar]
- Oosterhuis D. M., Loka D. A., Kawakami E. M., Pettigrew W. T. (2014). “The physiology of potassium in crop production,” in Advances in agronomy (Amsterdam: Elsevier; ), 203–233. doi: 10.1016/B978-0-12-800132-5.00003-1 [CrossRef] [Google Scholar]
- Papastylianou P., Kakabouki I., Travlos I. (2018). Effect of nitrogen fertilization on growth and yield of industrial hemp (Cannabis sativa L.). Not Bot. Horti Agrobo 46, 197–201. doi: 10.15835/nbha46110862 [CrossRef] [Google Scholar]
- Patel B., Wene D., Fan Z. T. (2017). Qualitative and quantitative measurement of cannabinoids in cannabis using modified HPLC/DAD method. J. Pharm. Biomed. Anal. 146, 15–23. doi: 10.1016/j.jpba.2017.07.021 [PubMed] [CrossRef] [Google Scholar]
- Prasad R., Power J. F. (1997). Soil fertility management for sustainable agriculture. (Boca Raton: CRC/Lewis Publishers; ). [Google Scholar]
- Raghothama K. G. (2015). Phosphorus and plant nutrition: an overview. Phosphorus: Agriculture and the Environment 46, 353–378. doi: 10.2134/agronmonogr46.c11 [CrossRef] [Google Scholar]
- Rioba N. B., Itulya F. M., Saidi M., Dudai N., Bernstein N. (2015). Effects of nitrogen, phosphorus and irrigation frequency on essential oil content and composition of sage (Salvia officinalis L.). J. Appl. Res. Medicinal Aromatic Plants 2, 21–29. doi: 10.1016/j.jarmap.2015.01.003 [CrossRef] [Google Scholar]
- Roberts S. K. (2006). Plasma membrane anion channels in higher plants and their putative functions in roots. New Phytol. 169, 647–666. doi: 10.1111/j.1469-8137.2006.01639.x [PubMed] [CrossRef] [Google Scholar]
- Saloner A., Bernstein N. (2021). Nitrogen supply affects cannabinoid and terpenoid profile in medical cannabis (Cannabis sativa L.). Ind. Crops Products 167, 113516. doi: 10.1016/j.indcrop.2021.113516 [CrossRef] [Google Scholar]
- Saloner A., Bernstein N. (2022). Effect of potassium (K) supply on cannabinoids, terpenoids and plant function in medical cannabis. Agronomy 12, 1242. doi: 10.3390/agronomy12051242 [CrossRef] [Google Scholar]
- Saloner A., Sacks M. M., Bernstein N. (2019). Response of medical cannabis (Cannabis sativa L.) genotypes to K supply under long photoperiod. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01369 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Sample E. C., Soper R. J., Racz G. J. (2015). “Reactions of phosphate fertilizers in soils,” in ASA, CSSA, and SSSA books. Eds. Khasawneh F. E., Sample E. C., Kamprath E. J. (Muscle Shoals: Madison, WI, USA: American Society of Agronomy, Crop Science Society of America, Soil Science Society of America; ), 263–310. doi: 10.2134/1980.roleofphosphorus.c12 [CrossRef] [Google Scholar]
- San Martín-Hernández C., Gómez-Merino F. C., Saucedo-Veloz C., Quintana-Obregón E. A., Muy-Rangel M. D., Trejo-Téllez L. I. (2021). Nitrogen and potassium supplied by phenological stages affect the carotenoid and nutritive content of the tomato fruit. Not Bot. Horti Agrobo 49, 12320. doi: 10.15835/nbha49212320 [CrossRef] [Google Scholar]
- Šenkyřík J. B., Křivánková T., Kaczorová D., Štefelová N. (2023). Investigation of the Effect of the Auxin Antagonist PEO-IAA on Cannabinoid Gene Expression and Content in Cannabis sativa L. Plants under In Vitro Conditions. Plants 12, 1664. doi: 10.3390/plants12081664 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Shi P., Song C., Chen H., Duan B., Zhang Z., Meng J. (2018). Foliar applications of iron promote flavonoids accumulation in grape berry of Vitis vinifera cv. Merlot grown in the iron deficiency soil. Food Chem. 253, 164–170. doi: 10.1016/j.foodchem.2018.01.109 [PubMed] [CrossRef] [Google Scholar]
- Shiponi S., Bernstein N. (2021. a). Response of medical cannabis (Cannabis sativa L.) genotypes to P supply under long photoperiod: Functional phenotyping and the ionome. Ind. Crops Products 161, 113154. doi: 10.1016/j.indcrop.2020.113154 [CrossRef] [Google Scholar]
- Shiponi S., Bernstein N. (2021. b). The highs and lows of P supply in medical cannabis: effects on cannabinoids, the ionome, and morpho-physiology. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.657323 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Toonen M., Ribot S., Thissen J. (2006). Yield of illicit indoor cannabis cultivation in the Netherlands. J. Forensic Sci. 51, 1050–1054. doi: 10.1111/j.1556-4029.2006.00228.x [PubMed] [CrossRef] [Google Scholar]
- Vanhove W., Van Damme P., Meert N. (2011). Factors determining yield and quality of illicit indoor cannabis (Cannabis spp.) production. Forensic Science International 212, 158–163. doi: 10.1016/j.forsciint.2011.06.006 [PubMed] [CrossRef] [Google Scholar]
- Varo M. A., Martin-Gomez J., Serratosa M. P., Merida J. (2022). Effect of potassium metabisulphite and potassium bicarbonate on color, phenolic compounds, vitamin C and antioxidant activity of blueberry wine. LWT 163, 113585. doi: 10.1016/j.lwt.2022.113585 [CrossRef] [Google Scholar]
- Velazquez L. A. (2013). First advances on the development of a hydroponic system for cherry tomato culture. In 2013 10th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE), 155–159. doi: 10.1109/ICEEE.2013.6676029 [CrossRef] [Google Scholar]
- Vera C. L., Malhi S. S., Phelps S. M., May W. E., Johnson E. N. (2010). N, P, and S fertilization effects on industrial hemp in Saskatchewan. Can. J. Plant Sci. 90, 179–184. doi: 10.4141/CJPS09101 [CrossRef] [Google Scholar]
- Yep B., Gale N. V., Zheng Y. (2020). Comparing hydroponic and aquaponic rootzones on the growth of two drug-type Cannabis sativa L. cultivars during the flowering stage. Ind. Crops Products 157, 112881. doi: 10.1016/j.indcrop.2020.112881 [CrossRef] [Google Scholar]
- Yep B., Zheng Y. (2021). Potassium and micronutrient fertilizer addition in a mock aquaponic system for drug-type Cannabis sativa L. cultivation. Can. J. Plant Sci. 101, 341–352. doi: 10.1139/cjps-2020-0107 [CrossRef] [Google Scholar]
- Zhao J., Zhu W.-H., Hu Q. (2001). Selection of fungal elicitors to increase indole alkaloid accumulation in catharanthus roseus suspension cell culture. Enzyme Microbial Technol. 28, 666–672. doi: 10.1016/S0141-0229(01)00309-X [PubMed] [CrossRef] [Google Scholar]
- Zheng Y. (2021). “Advances in understanding plant root behaviour and rootzone management in soilless culture systems,” in Advances in horticultural soilless culture. Ed. Gruda N. S. (London: Burleigh Dodds Science Publishing; ), 23–44. doi: 10.1201/9781003048206-3 [CrossRef] [Google Scholar]
- Zheng Y. (2022). Handbook of cannabis production in controlled environments. 1st ed (Boca Raton: CRC Press; ). doi: 10.1201/9781003150442 [CrossRef] [Google Scholar]
- Zhu B., Guo H., Cao Y., An R., Shi Y. (2021). Perceived importance of factors in cannabis purchase decisions: A best-worst scaling experiment. Int. J. Drug Policy 91, 102793. doi: 10.1016/j.drugpo.2020.102793 [PMC free article] [PubMed] [CrossRef] [Google Scholar]