Effects of Temperature, Relative Humidity, and Carbon Dioxide Concentration on Growth and Glucosinolate Content of Kale Grown in a Plant Factory

The growth of plants and their glucosinolate content largely depend on the cultivation environment; however, there are limited reports on the optimization of ambient environmental factors for kale grown in plant factories. This study was conducted to investigate the effects of temperature, relative humidity, and the carbon dioxide (CO2) concentration on kale growth and glucosinolate content in different growth stages of cultivation in a plant factory. Kale was grown under different temperatures (14, 17, 20, 23, and 26 °C), relative humidities (45, 55, 65, 75, and 85%), and CO2 concentrations (400, 700, 1000, 1300, and 1600 ppm) in a plant factory. Two and four weeks after transplantation, leaf samples were collected to evaluate the physical growth and glucosinolate contents. The statistical significance of the treatment effects was determined by two-way analysis of variance, and Duncan’s multiple range test was used to compare the means. A correlation matrix was constructed to show possible linear trends among the dependent variables. The observed optimal temperature, relative humidity, and CO2 range for growth (20–23 °C, 85%, and 700–1000 ppm) and total glucosinolate content (14–17 °C, 55–75%, and 1300–1600 ppm) were different. Furthermore, the glucosinolate content in kale decreased with the increase of temperature and relative humidity levels, and increased with the increase of CO2 concentration. Most of the physical growth variables showed strong positive correlations with each other but negative correlations with glucosinolate components. The findings of this study could be used by growers to maintain optimum environmental conditions for the better growth and production of glucosinolate-rich kale leaves in protected cultivation facilities.


Introduction
Kale (Brassica oleracea var. alboglabra Bailey) is a salad species that is one of the most versatile and commercially valuable vegetables due to its short growth period, various uses, and desirable metabolic and nutritional profiles [1][2][3]. This crisp and hearty vegetable is often consumed raw in salads and smoothies but can also be consumed in steamed, sautéed, or cooked states. Kale originates from China and has since gained particular attention in other countries due to its constituent cancer-preventive and human-healthpromoting phytochemicals (i.e., glucosinolates, carotenoids, phenols, and vitamins) [4][5][6]. Glucosinolates are amino-acid-derived, active secondary metabolites that mainly contain sulfur-and nitrogen-related compounds (i.e., β-D-thioglucose, tryptophan, phenylalanine, sulfonated oxime moiety). They can be classified into aliphatic, aromatic, and indole groups [7], where each group consists of several chemical constituents. Progoitrin, sinigrin, glucoraphanin, and gluconapin are the major constituents of the aliphatic group. Similarly, 4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin and neoglucobrassicin, and gluconasturtiin are the major indole and aromatic glucosinolate constituents, respectively [1]. Glucosinolates are composed of relatively few amino acids and chain-elongated homologs through an independent metabolic pathway ( Figure 1) and are available in all parts of almost all varieties of plants of the Brassicales order; however, the content is higher in the reproductive tissues (i.e., flowers and seeds) than in vegetative tissues [8]. The breakdown products of glucosinolates have a significant amount of anticarcinogenic activity for decreasing the risk of developing lung, stomach, colon, and rectum cancers; helping to maintain low blood pressure and reducing the risk of developing type 2 diabetes [6,9,10]. groups [7], where each group consists of several chemical constituents. Progoitrin, s nigrin, glucoraphanin, and gluconapin are the major constituents of the aliphatic grou Similarly, 4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin and neoglu cobrassicin, and gluconasturtiin are the major indole and aromatic glucosinolate consti uents, respectively [1]. Glucosinolates are composed of relatively few amino acids an chain-elongated homologs through an independent metabolic pathway ( Figure 1) and a available in all parts of almost all varieties of plants of the Brassicales order; however, th content is higher in the reproductive tissues (i.e., flowers and seeds) than in vegetativ tissues [8]. The breakdown products of glucosinolates have a significant amount of ant carcinogenic activity for decreasing the risk of developing lung, stomach, colon, and re tum cancers; helping to maintain low blood pressure and reducing the risk of developin type 2 diabetes [6,9,10].  [11].
Kale growth and the formation of glucosinolates depend on crop genetic factors, ti sue type, crop health, agronomic factors (i.e., water supply and fertigation), cultivatio facilities (i.e., plant factory, greenhouse, and open field), and environmental factors suc as temperature, relative humidity, carbon dioxide (CO2), light type, intensity, photope iod, and cultivation methods [12][13][14]. The physical development stage is also a major d terminant of the glucosinolates composition in kale [15]. Although kale can be easily cu tivated in open fields using traditional methods, the quality and quantity of the grow and glucosinolate content cannot be ensured, as they are extremely sensitive to climat and field conditions [16]. In recent years, farmers have produced kale in protected cult vation facilities, such as plant factories and greenhouses, due to the possibility of adjustin the growth environment and achieving fast and sustainable growth rates, functional com ponent-rich and high-quality yield, lower rates of disease and pest infestation, and low labor costs in addition to the possibility of year-round production with minimum influ ence from geological and climatic conditions [17][18][19]. Moreover, hydroponic cultivatio systems with ion-specific (ISE-sensor-based) nutrient management could enhance th growth and nutritional profile of kale by 15 to 60% [20][21][22][23][24][25]. However, major environme tal factors (i.e., temperature, relative humidity, and CO2) have to be specifically optimize according to crop to ensure sustainable kale growth and glucosinolate formation.
The physical growth of kale can be easily determined by measuring its physical prop erties such as plant length, width, weight, number of leaves, and stem diameter, where the glucosinolate content needs to be identified by laboratory analysis. The deposition  [11].
Kale growth and the formation of glucosinolates depend on crop genetic factors, tissue type, crop health, agronomic factors (i.e., water supply and fertigation), cultivation facilities (i.e., plant factory, greenhouse, and open field), and environmental factors such as temperature, relative humidity, carbon dioxide (CO 2 ), light type, intensity, photoperiod, and cultivation methods [12][13][14]. The physical development stage is also a major determinant of the glucosinolates composition in kale [15]. Although kale can be easily cultivated in open fields using traditional methods, the quality and quantity of the growth and glucosinolate content cannot be ensured, as they are extremely sensitive to climatic and field conditions [16]. In recent years, farmers have produced kale in protected cultivation facilities, such as plant factories and greenhouses, due to the possibility of adjusting the growth environment and achieving fast and sustainable growth rates, functional component-rich and high-quality yield, lower rates of disease and pest infestation, and lower labor costs in addition to the possibility of year-round production with minimum influence from geological and climatic conditions [17][18][19]. Moreover, hydroponic cultivation systems with ion-specific (ISE-sensor-based) nutrient management could enhance the growth and nutritional profile of kale by 15% to 60% [20][21][22][23][24][25]. However, major environmental factors (i.e., temperature, relative humidity, and CO 2 ) have to be specifically optimized according to crop to ensure sustainable kale growth and glucosinolate formation.
The physical growth of kale can be easily determined by measuring its physical properties such as plant length, width, weight, number of leaves, and stem diameter, whereas the glucosinolate content needs to be identified by laboratory analysis. The deposition of glucosinolates in growing plants and their distribution to plant organs are significantly affected by environmental factors [26], with temperature being one of the key factors. Several studies have been conducted to determine the process and effects of temperature on seed germination, physical development, flower formation, and yield [27][28][29][30]. However, physiological processes and their integration are sped up under higher temperatures with both positive and negative effects. For example, high temperatures promote faster growth and greater fruit production of plants, especially in cereal crops, but they also remove functional components from leaves through high transpiration rates [31]. Generally, elevated temperatures affect the structural components of chloroplasts significantly, causing effects such as variation in thylakoids, granum stacking, and swelling with photosystem II reduction, resulting in disruption to the cellular cytoplasm, cell breakdown and, ultimately, cell death. In addition, rising temperatures interrupt protein mechanisms, RNA synthesis, enzymatic interactions, and cell function. As a result, these imbalances and abnormal cell functioning affect the growth and accumulation of glucosinolate synthesis [32,33].
The relative humidity of the ambient environment also directly affects plant growth by resisting water and nutrient consumption. During transpiration, the relative humidity level becomes saturated. As a result, plants halt transpiration and nutrient uptake from the soil or growing media at high relative humidity levels where there is a lack of air circulation, resulting in gradual rotting in cases of long-term humidity saturation [34][35][36]. The maintenance of optimum relative humidity is essential for better growth and glucosinolate accumulation. Several researchers have reported that the photosynthesis rate is proportional to the relative humidity level as a higher range of relative humidity lowers water stress in the leaves and increases stomatal conductance. Although higher relative humidity increases the nutrient concentration, the nutrient solution supply and plant transpiration rate need to be monitored carefully [34,37].
The CO 2 concentration influences the photosynthetic rate, metabolism, and physiological and chemical defense of plants [13,38]. A lack of CO 2 would not only result in a lower biomass but the plants would also be of inferior quality and strength. As an essential substrate of the photosynthesis process, CO 2 is directly absorbed by plants. CO 2 also influences the transpiration process of plants. A meta-analysis was conducted, and it was reported that elevated CO 2 could reduce transpiration by up to 22% in different plant species [39]. CO 2 also preserves the essential nutrient components along with water by reducing the transpiration rate [40,41]. La et al. [38] investigated the effects of CO 2 elevation at different nitrogen levels on the growth and glucosinolate content of Chinese kale and reported that all physical growth variables significantly increased with the elevation of CO 2 at each nitrogen level; however, total glucosinolate content was only increased under low nitrogen level and elevated CO 2 concentration.
The temperature, relative humidity, and CO 2 concentration are the basic environmental factors that affect kale growth and, especially, glucosinolate formation. As they are interrelated, these factors should not be studied in isolation. The proper combination of these factors needs to be specifically confirmed for each crop to ensure optimal growth, a favorable nutritional profile, and identification of the ideal harvesting time. To date, very few studies have investigated the effects of these environmental factors on kale, especially when grown in plant factories using hydroponic cultivation methods. Therefore, the objective of this study was to investigate the effects of temperature, relative humidity, and CO 2 on the growth and glucosinolate content at different stages of kale growth based on cultivation in a plant factory.

Plant Factory and Seedling Preparation
Plant factories are fully-closed crop cultivation systems that are fitted with artificial lights and used to grow high-value vegetables and medicinal plants throughout the year by utilizing artificially controlled ambient environmental factors [17,19]. In this study, five small chambers were prepared, as shown in Figure 2, to implement five different treatment conditions with varied temperature, relative humidity, and CO 2 concentrations. The targeted environmental factors (i.e., temperature, relative humidity, and CO 2 ) and other environmental factors (i.e., light sources, light intensity, photoperiod, and nutrient solution (EC and pH)) were maintained according to the experimental plan (Table 1). A wireless sensor network (XBee-Pro, Digi, Hopkins, MN, USA) was used to monitor the ambient environmental parameters and control the relevant actuators, as detailed by Chung et al. [42]. Three plant beds were placed vertically in each cultivation chamber, and a nutrient solution tank was kept at the bottom (floor). Each plant bed had 24 planting positions and 6 mist spray nozzles for spraying the nutrient solution onto plant roots as a fine mist for a duration of 2 min at 13-min intervals. Commercial nutrient solutions A and B (Daeyu Co., Ltd., Seoul, Korea) were used, and the target nutrient level was monitored and managed once a day using EC and pH sensors. small chambers were prepared, as shown in Figure 2, to implement five different treatment conditions with varied temperature, relative humidity, and CO2 concentrations. The targeted environmental factors (i.e., temperature, relative humidity, and CO2) and other environmental factors (i.e., light sources, light intensity, photoperiod, and nutrient solution (EC and pH)) were maintained according to the experimental plan (Table 1). A wireless sensor network (XBee-Pro, Digi, Hopkins, MN, USA) was used to monitor the ambient environmental parameters and control the relevant actuators, as detailed by Chung et al. [42]. Three plant beds were placed vertically in each cultivation chamber, and a nutrient solution tank was kept at the bottom (floor). Each plant bed had 24 planting positions and 6 mist spray nozzles for spraying the nutrient solution onto plant roots as a fine mist for a duration of 2 min at 13-min intervals. Commercial nutrient solutions A and B (Daeyu Co., Ltd., Seoul, Korea) were used, and the target nutrient level was monitored and managed once a day using EC and pH sensors.   A commercial kale variety with green smooth leaves and a hard and fibered stem was cultivated in the experiments using a recycle-type aeroponic nutrient management system. Kale seeds were sown in a hydroponic germination sponge, covered with wet paper (until germination), kept in the plant factory under a controlled environment for germination, and grown until transplantation ( Figure 3). The maintained temperature, relative humidity, CO 2 concentration, light type, and photoperiod were 25 ± 3 • C, 65 ± 5%, 1000 ± 100 ppm, fluorescent, and 16/8 (day/night hours), respectively. Nutrient-rich water was provided into the root zone, and the EC and pH of the nutrient solution were 1.2 ± 1.00 (dS m −1 ) and 6.5 ± 0.5, respectively. After three weeks of germination, healthy seedlings with true leaves were transplanted into the plant bed with the sponge. A commercial kale variety with green smooth leaves and a hard and fibered stem was cultivated in the experiments using a recycle-type aeroponic nutrient management system. Kale seeds were sown in a hydroponic germination sponge, covered with wet paper (until germination), kept in the plant factory under a controlled environment for germination, and grown until transplantation ( Figure 3). The maintained temperature, relative humidity, CO2 concentration, light type, and photoperiod were 25 ± 3 °C, 65 ± 5%, 1000 ± 100 ppm, fluorescent, and 16/8 (day/night hours), respectively. Nutrient-rich water was provided into the root zone, and the EC and pH of the nutrient solution were 1.2 ± 1.00 (dS m −1 ) and 6.5 ± 0.5, respectively. After three weeks of germination, healthy seedlings with true leaves were transplanted into the plant bed with the sponge.

Experimental Design
Different separate experiments were conducted to investigate the influences of temperature, relative humidity, and CO2 on kale growth and glucosinolate content. Five treatments with various environmental factors were applied in each experiment. For example, temperatures of 14, 17, 20, 23, and 26 °C were varied while other factors were kept constant. Similarly, five relative humidity levels and CO2 concentrations were implemented in experiments 2 and 3 to evaluate the effects of relative humidity and CO2, respectively. The targeted and monitored levels of temperature, relative humidity, and CO2 along with other growth factors are summarized in Table 1. The light source and ratio, intensity, photoperiod, and pH and EC levels were selected following the findings of Zhang et al. [43], Lefsrud et al. [44], Naznin et al. [45], and Jones [24], respectively.

Experimental Design
Different separate experiments were conducted to investigate the influences of temperature, relative humidity, and CO 2 on kale growth and glucosinolate content. Five treatments with various environmental factors were applied in each experiment. For example, temperatures of 14, 17, 20, 23, and 26 • C were varied while other factors were kept constant. Similarly, five relative humidity levels and CO 2 concentrations were implemented in experiments 2 and 3 to evaluate the effects of relative humidity and CO 2 , respectively. The targeted and monitored levels of temperature, relative humidity, and CO 2 along with other growth factors are summarized in Table 1. The light source and ratio, intensity, photoperiod, and pH and EC levels were selected following the findings of Zhang et al. [43], Lefsrud et al. [44], Naznin et al. [45], and Jones [24], respectively.

Sample Collection and Data Acquisition
Two and four weeks after transplantation (Figure 4), sample collection was performed in two steps. First, mature and healthy plants were visually selected and collected from the plant beds for physical growth evaluation. Three plants from each bed and nine plants from three beds (replicates) of each cultivation chamber were collected randomly among 72 plants (24 plants/bed × 3 plant beds). To analyze the glucosinolate content, three normal-sized, mature, healthy leaves were harvested from each plant bed (one leaf from each collected plant), and a total of nine leaves were collected from three plant beds (as a replication) from each cultivation chamber. The measured values for each growth parameter and the glucosinolate content were averaged to represent one data point. As a result, nine data points for growth parameters and one data point for the glucosinolate content were recorded from each plant bed. In total, 270 data points were collected for growth evaluation (9 data points/bed × 2 sampling times × 3 replications × 5 treatments) and 30 data points were collected to assess the glucosinolate content (1 data point/bed × 2 sampling times × 3 replications × 5 treatments) for each experiment. The physical growth variables, namely, the plant height, width, weight, number of leaves, stem diameter, chlorophyll level, leaf length, width, and weight were measured, and the leaves were transferred to the chemical laboratory immediately (to minimize the degradation) for glucosinolate analysis using a commercial high-performance liquid chromatography (HPLC) machine (model: 1200 series, Agilent Technologies, Santa Clara, CA, USA). The chlorophyll concentration was also measured using a commercial device (model: SPAD 502DL, Spectrum Technology Inc., Aurora, IL, USA).

Sample Collection and Data Acquisition
Two and four weeks after transplantation (Figure 4), sample collection was performed in two steps. First, mature and healthy plants were visually selected and collected from the plant beds for physical growth evaluation. Three plants from each bed and nine plants from three beds (replicates) of each cultivation chamber were collected randomly among 72 plants (24 plants/bed × 3 plant beds). To analyze the glucosinolate content, three normal-sized, mature, healthy leaves were harvested from each plant bed (one leaf from each collected plant), and a total of nine leaves were collected from three plant beds (as a replication) from each cultivation chamber. The measured values for each growth parameter and the glucosinolate content were averaged to represent one data point. As a result, nine data points for growth parameters and one data point for the glucosinolate content were recorded from each plant bed. In total, 270 data points were collected for growth evaluation (9 data points/bed × 2 sampling times × 3 replications × 5 treatments) and 30 data points were collected to assess the glucosinolate content (1 data point/bed × 2 sampling times × 3 replications × 5 treatments) for each experiment. The physical growth variables, namely, the plant height, width, weight, number of leaves, stem diameter, chlorophyll level, leaf length, width, and weight were measured, and the leaves were transferred to the chemical laboratory immediately (to minimize the degradation) for glucosinolate analysis using a commercial high-performance liquid chromatography (HPLC) machine (model: 1200 series, Agilent Technologies, Santa Clara, CA, USA). The chlorophyll concentration was also measured using a commercial device (model: SPAD 502DL, Spectrum Technology Inc., Aurora, IL, USA).

Estimation of Glucosinolate Content
The glucosinolates of the freshly harvested kale leaves were extracted and analyzed as described by Doheny-Adams et al. [46]. The whole process was conducted according to the ISO 9167:2019 [47], and the process is divided into four major steps: (a) tissue disruption, (b) extraction in methanol, (c) purification and desulfation, and (d) separation and identification of glucosinolates by HPLC analysis (1200 series, Agilent Technologies, CA, USA). The collected leaf samples were stored in an airtight box, taken to the chemical laboratory immediately to freeze in liquid nitrogen, and stored at −80 °C for 48 h to reduce the activity of myrosinase. For freeze-drying, samples were lightly wrapped with aluminum foil and transported on dry ice to load into the freeze drier (Lyotrap, LTE scientific Ltd., Oldham, UK) within 30 s. The freeze-dried leaf samples were ground to make a homogenized fine powder using a grinder (EK2311, Salter, Tonbridge, UK). Then, 100 mg of the freeze-dried samples was preheated for 3 min at 75 °C and 4.5 mL of preheated 70% methanol at 75 °C was added. The sample was incubated for 10 min at 75 °C (with manual shaking every 2 min) and then centrifuged by a rotor at 4000 rpm (B 3.11, Jouan, Nantes, France) for 10 min. In the purification step, 25 mg of sulfatase and 1ml of 40% ethanol were mixed and centrifuged for 1 min at 8000 rpm. The supernatant was shifted to a new Eppendorf tube and 1 mL of pure ethanol was injected for precipitating the sulfatase before the second centrifugation. Finally, the sulfatase pellet was air dried after separating

Estimation of Glucosinolate Content
The glucosinolates of the freshly harvested kale leaves were extracted and analyzed as described by Doheny-Adams et al. [46]. The whole process was conducted according to the ISO 9167:2019 [47], and the process is divided into four major steps: (a) tissue disruption, (b) extraction in methanol, (c) purification and desulfation, and (d) separation and identification of glucosinolates by HPLC analysis (1200 series, Agilent Technologies, Santa Clara, CA, USA). The collected leaf samples were stored in an airtight box, taken to the chemical laboratory immediately to freeze in liquid nitrogen, and stored at −80 • C for 48 h to reduce the activity of myrosinase. For freeze-drying, samples were lightly wrapped with aluminum foil and transported on dry ice to load into the freeze drier (Lyotrap, LTE scientific Ltd., Oldham, UK) within 30 s. The freeze-dried leaf samples were ground to make a homogenized fine powder using a grinder (EK2311, Salter, Tonbridge, UK). Then, 100 mg of the freeze-dried samples was preheated for 3 min at 75 • C and 4.5 mL of preheated 70% methanol at 75 • C was added. The sample was incubated for 10 min at 75 • C (with manual shaking every 2 min) and then centrifuged by a rotor at 4000 rpm (B 3.11, Jouan, Nantes, France) for 10 min. In the purification step, 25 mg of sulfatase and 1ml of 40% ethanol were mixed and centrifuged for 1 min at 8000 rpm. The supernatant was shifted to a new Eppendorf tube and 1 mL of pure ethanol was injected for precipitating the sulfatase before the second centrifugation. Finally, the sulfatase pellet was air dried after separating from the supernatant and diluted in 2 mL of water. For desulfation, 0.5 cc of Sephadex slurry was used to prepare the columns and 2 mL of imizadole formate (6 M) was added on each for activation. The columns were cleaned twice with 1 mL of water each time. The columns were washed again using 1 mL of 20 mM sodium acetate, and 75 µL of purified sulfatase (0.05-0.3 U/mL) was injected. After that, columns were incubated for 24 h at 28 • C before desulfoglucosinolates were eluted with two 1 mL volumes of water. After 24 h of incubation, elution of desulfoglucosinolates was performed thrice using 1.5 mL of distilled water and filtered through 0.45-µm polytetrafluoroethylene (PTFE) syringe filters (Millipore, Bedford, MA, USA) into an HPLC vial. A reverse phase C18 column (150 × 3.0 mm, 3 µm, Inertsil ODS-3, GL Sciences, Tokyo, Japan) was used, which was equilibrated for 30 min using ultrapure water (solvent A) and 100% acetonitrile (solvent B) with detection at 227 nm. The flow rate was 0.4 mL min −1 , and separation was performed according to the default program. As an external standard, sinigrin (0.1 mg/mL; Sigma, St. Louis, MO, USA) was utilized. The identification and quantification of individual glucosinolate components was performed by comparing the sinigrin retention time and using their HPLC areas and response factor, respectively. In this study, the obtained retention time for progoitrin, sinigrin, glucobrassicin, 4-methoxyglucobrassicin, and neoglucobrassicin were 5.97, 7.13, 21.93, 24.68, and 30.37 min, respectively.

Statistical Analysis
All the presented physical growth parameters and glucosinolate content values are the means of independent measurements for different treatments of each environmental factor. The significance of differences between mean values was determined by two-way analysis of variance (ANOVA). Data were analyzed considering 95% confidence levels and two-sided confidence intervals. Duncan's multiple range test was used to simultaneously compare means (SAS Institute Inc, Campus drive Cary, NC, USA). A correlation matrix recording correlation coefficients was created to show the inter-relationships between variables.

ANOVA of the Environmental Factors
The effects of ambient environmental factors (temperature, relative humidity, and CO 2 ) on plant physical growth variables and total glucosinolate content were analyzed using two-way ANOVA analysis. Five different treatment conditions for each environmental factor and two sampling times were considered when conducting the ANOVA analysis for each growth variable and the glucosinolate content. The results of the two-way ANOVA analysis for the plant height, width, weight, and total glucosinolate content are shown in Table 2 out of nine physical variables and five identified glucosinolate components. The F-values of the treatments and sampling times were higher than the F crit values, except for some growth and glucosinolate variables under the CO 2 treatments, which confirms the adequacy of the hypothesis. This ANOVA analysis indicates that the treatments and sampling times had significant impacts (p < 0.05) on the growth and glucosinolate content (except for some CO 2 treatments). However, some P-values under the CO 2 treatments were higher than 0.05, which also indicates that those growth or glucosinolate variables were not notably affected by the CO 2 treatments. The overall results show that a single unit change of each environmental factor will affect the plant growth and glucosinolate content.  Table 3 shows the magnitude, direction, and linear pairwise relationship between the identified glucosinolate variables under the considered ambient environmental factors (temperature, relative humidity, and CO 2 ). Among the five identified glucosinolate variables under the temperature experiments, sinigrin and glucobrassicin were strongly and positively correlated among them and identified glucosinolate variables under the relative humidity experiments, and they were strongly and negatively correlated with the variables identified in the CO 2 experiments. The correlations were statistically significant at a 0.1% level (except for some variables). Although, progoitrin had a significant positive correlation with each of the five identified glucosinolate variables under the CO 2 experiments, no significant correlations were observed with other variables. Except for some strong correlations, 4-methoxyglucobrassicin, and neoglucobrassicin were also not significantly correlated with other identified glucosinolate variables. Strong negative correlations with a 0.1% significance level were observed for most of the identified glucosinolate variables under the relative humidity and CO 2 experiments. However, identified glucosinolate variables under the CO 2 experiments were strongly and positively correlated. They were statistically significant at a 0.1% level (except for the C_Sin). The multicollinearity issue can also be predicted from the correlation matrix. A highly correlated value (>0.7) hinders the evaluation of the true effects of the predictor variables. According to Table 3, some of the glucosinolate variables had notable evidence of strong correlation. For example, C_Pro showed positive correlations of 0.90, 0.90, and 0.99 with C_Glu, C_4-met, and C_Neo, respectively, and C_4-met showed negative correlations of −0.91 and −0.93, with T_Sin and T_Glu, respectively. Variance inflation factor (VIF) was also investigated and the values varied from 1 to 3 for most of the variables, indicating that the variables were slightly explained by other independent variables. However, the VIF values of T_Sin (4.72) and T_Glu (6.46) were relatively high [48].

Evaluation of Temperature Effects
A statistical analysis was conducted to evaluate the effects of temperature on kale growth, and the results are shown in Table 4. Regarding kale physical properties, an overall high growth rate was observed at 20-23 • C, and the lowest growth rate occurred at 14 • C. However, some physical parameters showed greater numbers at 17 and 26 • C. They were plant height (17 • C) and width (26 • C) after two weeks of transplantation, and chlorophyll level (17 • C), leaf length and weight (26 • C) after four weeks of transplantation. The data points of no. of leaves, stem diameter, and leaf parameters (length, width, and weight) were very close to the mean (low standard deviation); however, the data points of other growth variables, especially plant height and width, were spread out over a wide range of values. Standard deviation was greater in samples collected after four weeks of cultivation, compared to the two weeks. According to Duncan's range test results, significant differences were observed for the plant width, weight, and leaf parameters (length, width, and weight) at 2-week sampling time, and the plant height, and leaf parameters at 4-week sampling time, depending on the temperature levels. Contrariwise, the rest of the growth variables (specifically the number of leaves, stem diameter, and chlorophyll level) did not show statistical significance regarding the temperature variations.

Evaluation of Temperature Effects
A statistical analysis was conducted to evaluate the effects of temperature on kale growth, and the results are shown in Table 4. Regarding kale physical properties, an overall high growth rate was observed at 20-23 °C, and the lowest growth rate occurred at 14 °C. However, some physical parameters showed greater numbers at 17 and 26 °C. They were plant height (17 °C) and width (26 °C) after two weeks of transplantation, and chlorophyll level (17 °C), leaf length and weight (26 °C) after four weeks of transplantation. The data points of no. of leaves, stem diameter, and leaf parameters (length, width, and weight) were very close to the mean (low standard deviation); however, the data points of other growth variables, especially plant height and width, were spread out over a wide range of values. Standard deviation was greater in samples collected after four weeks of cultivation, compared to the two weeks. According to Duncan's range test results, significant differences were observed for the plant width, weight, and leaf parameters (length, width, and weight) at 2-week sampling time, and the plant height, and leaf parameters at 4-week sampling time, depending on the temperature levels. Contrariwise, the rest of the growth variables (specifically the number of leaves, stem diameter, and chlorophyll level) did not show statistical significance regarding the temperature variations.   Figure 5 shows the effects of temperature on the glucosinolate content, based on various components, of harvested kale leaves after two and four weeks of transplantation. Glucobrassicin was found to be a dominant glucosinolate component in both cases. However, an inverse relationship was observed between the contents of all glucosinolate components and increased ambient temperature levels and cultivation period. The total glucosinolate content became lower at each increased temperature level, and the lowest total glucosinolate of kale leaves in each sampling time was observed at 26 • C. A high standard deviation trend was observed for each glucosinolate component due to the low sampling number. Among the five levels of temperature, the total glucosinolate content was higher at 14-17 • C in both cultivation periods. According to Duncan's range test, the concentrations of sinigrin, glucobrassicin, and 4-methoxyglucobrassicin were significantly different for each temperature level in samples collected after two weeks of cultivation; however, no significant differences were observed among the other glucosinolate components (except sinigrin at 4th week) for different temperature levels and cultivation periods.  Figure 5 shows the effects of temperature on the glucosinolate content, based on various components, of harvested kale leaves after two and four weeks of transplantation. Glucobrassicin was found to be a dominant glucosinolate component in both cases. However, an inverse relationship was observed between the contents of all glucosinolate components and increased ambient temperature levels and cultivation period. The total glucosinolate content became lower at each increased temperature level, and the lowest total glucosinolate of kale leaves in each sampling time was observed at 26 °C. A high standard deviation trend was observed for each glucosinolate component due to the low sampling number. Among the five levels of temperature, the total glucosinolate content was higher at 14-17 °C in both cultivation periods. According to Duncan's range test, the concentrations of sinigrin, glucobrassicin, and 4-methoxyglucobrassicin were significantly different for each temperature level in samples collected after two weeks of cultivation; however, no significant differences were observed among the other glucosinolate components (except sinigrin at 4 th week) for different temperature levels and cultivation periods.  The interactions of each growth variable and glucosinolate component under the different temperature treatments after four weeks of transplantation were analyzed using the correlation matrix and the results are summarized in Table 5, where the level of significance and VIF are also mentioned. In many cases, strong positive and negative correlation coefficients were observed. The physical growth variables were found to be strongly correlated with each other (around 0.50-0.98). However, the chlorophyll level showed negative correlations with all the growth variables, and positive correlations were observed with progoitrin, sinigrin, and glucobrassicin. The progoitrin, sinigrin, and glucobrassicin were strongly positively correlated with each other and negatively correlated with 4-methoxyglucobrassicin and neoglucobrassicin. Moreover, the detected glucosinolate components (except 4-methoxyglucobrassicin) were negatively correlated with most of the physical growth variables. According to the VIF analysis, most of the variables were moderately correlated; however, some variables (i.e., leaf length and glucobrassicin) were highly correlated, which might adversely affect other variables.

Evaluation of Relative Humidity Effects
The effects of relative humidity on kale growth properties are summarized in Table 6. The growth status was evaluated at two different stages (two and four weeks after transplantation). Most of the physical growth variables were prominent at the 85% relative humidity level at both sampling periods, except for the number of leaves, stem diameter, chlorophyll level (55%), and leaf length (75%) in the second week, and the number of leaves, leaf weight (65%), and chlorophyll level (45%) in the fourth week of cultivation. Except for some growth variables, the overall lowest growth performance was observed at the 45% relative humidity level in both sampling times. The data points of some growth variables, such as plant height and width, leaf length and width, were spread out over a wide range compared to other growth variables, and greater standard deviations were observed in samples collected after four weeks of cultivation, compared to the two weeks of cultivation. According to Duncan's range test, all the growth variables (except the chlorophyll level) were significantly different at the 2-week sampling time depending on the relative humidity levels. A similar result was observed (except for the plant weight, chlorophyll level, and leaf width) at the 4-week sampling time. the correlation matrix and the results are summarized in Table 5, where the level of significance and VIF are also mentioned. In many cases, strong positive and negative correlation coefficients were observed. The physical growth variables were found to be strongly correlated with each other (around 0.50-0.98). However, the chlorophyll level showed negative correlations with all the growth variables, and positive correlations were observed with progoitrin, sinigrin, and glucobrassicin. The progoitrin, sinigrin, and glucobrassicin were strongly positively correlated with each other and negatively correlated with 4-methoxyglucobrassicin and neoglucobrassicin. Moreover, the detected glucosinolate components (except 4-methoxyglucobrassicin) were negatively correlated with most of the physical growth variables. According to the VIF analysis, most of the variables were moderately correlated; however, some variables (i.e., leaf length and glucobrassicin) were highly correlated, which might adversely affect other variables.

Evaluation of Relative Humidity Effects
The effects of relative humidity on kale growth properties are summarized in Table  6. The growth status was evaluated at two different stages (two and four weeks after transplantation). Most of the physical growth variables were prominent at the 85% relative humidity level at both sampling periods, except for the number of leaves, stem diameter, chlorophyll level (55%), and leaf length (75%) in the second week, and the number of leaves, leaf weight (65%), and chlorophyll level (45%) in the fourth week of cultivation. Except for some growth variables, the overall lowest growth performance was observed at the 45% relative humidity level in both sampling times. The data points of some growth   The results of the glucosinolate analysis for different relative humidity treatments and cultivation periods are shown in Figure 6. The aliphatic glucosinolates (i.e., progoitrin, sinigrin) and indole glucosinolate (i.e., glucobrassicin) were the most prominent components at both of the sampling times. The overall glucosinolate concentrations decreased slightly in the samples collected after the fourth week of cultivation. A high standard deviation was observed, especially for the glucobrassicin, as the sample number was low and sometimes all glucosinolate components were not detected in some samples. According to Duncan's range test, no significant differences were observed among the glucosinolate components (except for the progoitrin, glucobrassicin, and 4-methoxyglucobrassicin at the 2-week sampling time) for different relative humidity treatments and cultivation periods. Table 7 shows the correlation matrix of physical and glucosinolate properties for different relative humidity treatments after four weeks of cultivation. All physical growth variables, except for the number of leaves and the chlorophyll level, showed strong positive correlations with one another. A fairly good correlation (both positive and negative) was observed between the physical variables and glucosinolate components. However, most of the glucosinolate components were negatively correlated with each other. Moreover, the VIF values were also determined. Except for the leaf length (VIF: 11.93), other predictors were moderately correlated, which resulted in a low influence on other independent variables.
The results of the glucosinolate analysis for different relative humidity treatments and cultivation periods are shown in Figure 6. The aliphatic glucosinolates (i.e., progoitrin, sinigrin) and indole glucosinolate (i.e., glucobrassicin) were the most prominent components at both of the sampling times. The overall glucosinolate concentrations decreased slightly in the samples collected after the fourth week of cultivation. A high standard deviation was observed, especially for the glucobrassicin, as the sample number was low and sometimes all glucosinolate components were not detected in some samples. According to Duncan's range test, no significant differences were observed among the glucosinolate components (except for the progoitrin, glucobrassicin, and 4-methoxyglucobrassicin at the 2-week sampling time) for different relative humidity treatments and cultivation periods.    Table 7 shows the correlation matrix of physical and glucosinolate properties for different relative humidity treatments after four weeks of cultivation. All physical growth variables, except for the number of leaves and the chlorophyll level, showed strong positive correlations with one another. A fairly good correlation (both positive and negative) was observed between the physical variables and glucosinolate components. However, most of the glucosinolate components were negatively correlated with each other. Moreover, the VIF values were also determined. Except for the leaf length (VIF: 11.93), other predictors were moderately correlated, which resulted in a low influence on other independent variables.

Evaluation of CO2 Effects
A summary of the effects of CO2 treatments and cultivation periods on kale growth is given in Table 8. The overall growth performance was higher under 700-1000 ppm CO2. However, a notable growth rate of some parameters (i.e., chlorophyll level and leaf length at the 2-week sampling time, and plant height and chlorophyll level at the 4-week sampling time) was observed under 400 ppm of CO2. Relatively low growth performance was observed under 1300 and 1600 ppm of CO2 in both sampling periods. Besides this, the spread of standard deviations of the growth variables was almost similar for both of the sampling periods. Comparatively high standard deviations were observed for the plant parameters (height, width, and weight) and chlorophyll level compared to other growth variables. Based on Duncan's range test results, there were no significant differences for the growth parameters (except the plant height at the 4-week sampling time) under different CO2 concentrations and cultivation periods.
The effects of different CO2 concentrations on the glucosinolate content are shown in Figure 7. The optimal CO2 level in relation to the total glucosinolate content was 1300 ppm at both the second and fourth weeks of cultivation. The progoitrin, sinigrin, and neoglu-

Evaluation of CO 2 Effects
A summary of the effects of CO 2 treatments and cultivation periods on kale growth is given in Table 8. The overall growth performance was higher under 700-1000 ppm CO 2 . However, a notable growth rate of some parameters (i.e., chlorophyll level and leaf length at the 2-week sampling time, and plant height and chlorophyll level at the 4-week sampling time) was observed under 400 ppm of CO2. Relatively low growth performance was observed under 1300 and 1600 ppm of CO 2 in both sampling periods. Besides this, the spread of standard deviations of the growth variables was almost similar for both of the sampling periods. Comparatively high standard deviations were observed for the plant parameters (height, width, and weight) and chlorophyll level compared to other growth variables. Based on Duncan's range test results, there were no significant differences for the growth parameters (except the plant height at the 4-week sampling time) under different CO 2 concentrations and cultivation periods.
The effects of different CO 2 concentrations on the glucosinolate content are shown in Figure 7. The optimal CO 2 level in relation to the total glucosinolate content was 1300 ppm at both the second and fourth weeks of cultivation. The progoitrin, sinigrin, and neoglucobrassicin contents decreased after the two weeks of cultivation. The low sampling number caused high standard deviations of the detected glucosinolate components. Glucobrassicin was found to be a dominant component in the samples collected after 4-weeks. The results of Duncan's range test showed a significant difference in glucosinolate components (except for the 4-methoxyglucobrassicin and neoglucobrassicin at the 2-week sampling, and sinigrin, glucobrassicin, and neoglucobrassicin at the 4-week sampling time) under different CO 2 concentrations and both cultivation periods.   Table 9 shows the interactions of each physical and functional parameter with one another, along with the significance levels for different CO2 treatments after four weeks of transplantation. Weak correlations (both positive and negative) were detected among most of the physical growth variables. However, the glucosinolate components (except sinigrin) were strongly positively correlated with each other and mostly negatively correlated with physical growth variables. According to the VIF analysis, the VIF values of the physical variables varied from 1 to 2, except plant height (4.23), indicating low correlations. Contrariwise, the VIF values of the glucosinolate components (except progoitrin) were comparatively high, which indicated highly correlated relationships and influences on other predictors.  Table 9 shows the interactions of each physical and functional parameter with one another, along with the significance levels for different CO 2 treatments after four weeks of transplantation. Weak correlations (both positive and negative) were detected among most of the physical growth variables. However, the glucosinolate components (except sinigrin) were strongly positively correlated with each other and mostly negatively correlated with physical growth variables. According to the VIF analysis, the VIF values of the physical variables varied from 1 to 2, except plant height (4.23), indicating low correlations. Contrariwise, the VIF values of the glucosinolate components (except progoitrin) were comparatively high, which indicated highly correlated relationships and influences on other predictors.

Discussion
There is an interaction between plant growth and glucosinolate concentration, which strongly depends on the environmental conditions and water-nutrient consumption rate, along with the plant species, growth method, cultivation period, and cultivation facilities used [12][13][14]16,49,50]. In this study, the growth rate of kale increased with the cultivation period. The overall maximum growth rate was observed at 20-23 • C, around 85% relative humidity, and 700-1000 ppm CO 2 (Tables 4, 6 and 8). The optimal temperature, relative humidity, and CO 2 range for total glucosinolate content were 14-17 • C, 55-75%, and 1300-1600 ppm. However, the glucosinolate content of kale decreased notably as cultivation period, temperature, and relative humidity level increased (Figures 5 and 6). Contrariwise, it increased with increased CO 2 concentration ( Figure 7). All biological processes of plants speed up at higher temperatures [51]. However, the sensitivity of plants to the atmospheric temperature depends on the growth stage. Plants always seek to maintain a balance between the plant-body temperature and air temperature. If the plant is heated up, the transpiration rate increases to cool down plants, which increases water and nutrient uptake, resulting in phonological changes in plants [52,53]. This assimilation process occurs quickly in the early growth stage. We observed a high growth rate at 23-26 • C in the 2nd week, which was reduced along with the temperature range (to 20-23 • C) in the 4th week. However, the rapid transpiration process also ejects many nutrient components, which lowers the concentration of glucosinolate components, as shown in Figure 5 [33]. Conversely, a high level of accumulation of functional components (i.e., glucosinolates) occurs at low temperatures. Steindal et al. [29] explained that low temperatures activate cold acclimatization processes, including many biochemical and physiological changes, to improve the cold tolerance capacity. These procedures reduce the growth and accumulation of osmolytes and the functional component composition. In this study, the lowest rate of physical growth and the highest concentration of glucosinolates were also observed at 14 • C (Table 4, Figure 5). Velasco et al. [49] reported an inverse relationship between low temperatures and the total glucosinolate content. Relative humidity is directly related to CO 2 acclimation through the stomata response, which is connected with plant growth and nutritional levels. Ahmed et al. [54] reviewed several studies and reported that a relative humidity of lower than 40% and higher than 85% causes stomatal malfunctioning, inhibiting the plant growth rate and photosynthesis. They also mentioned that the optimal range of relative humidity for leafy vegetables (i.e., lettuce) is 70-80%. In this study, maximal growth was found at a relative humidity range of 75-85%, and no significant difference was observed at a relative humidity range of 45-85% for the glucosinolate components ( Figure 6), which matches the findings of previous studies. In addition to the effects of temperature and relative humidity, a significant impact of the CO 2 concentration was observed on the accumulation of glucosinolate components rather than the growth rate of kale. In open environments, the concentration of CO 2 remains constant (300-400 ppm), but this concentration can be increased in protected cultivation facilities (i.e., greenhouses and plant factories). Usually, the demand for CO 2 increases with the increment of plant growth parameters and biomass [55]. In this study, CO 2 concentrations of 700 to 1000 ppm were associated with better growth performance (Table 8), and higher glucosinolate formation was observed under 1300 to 1600 ppm range of CO 2 (Figure 7). Higher concentrations of CO 2 help to synthesize larger amounts of carbohydrates and other functional components through photosynthesis [13,38,40]. Moreover, lower reduction of the photosynthetic ingredients under elevated CO 2 concentrations improves the glucosinolate content [39,41]. An overaccumulation of glucosinolates was observed under experiment-1 (temperature) compared to experiments-2 and 3 (relative humidity and CO 2 ). We know glucosinolates are significantly affected by the variety, genetics, plant growth stage, irrigation level, growing media, and environmental variables (i.e., temperature, humidity, CO 2 , and light conditions). For example, Chen et al. [56] investigated the variation of glucosinolates in Chinese Brassica campestris vegetables (Chinese cabbage, purple cai-tai, choysum, pakchoi, and turnip) and reported that total glucosinolates varied from 14-130 mg/100 g fresh weight (FW), where He et al. [57] observed the minimum (28.9 µmol/100 g FW) in broccoli and maximum (278 µmol/100 g FW) in Chinese kale. From seedling to early flowering, the total glucosinolate content increased with plant age in B. oleracea leaves. After that point, the aliphatic glucosinolate content decreased dramatically over time as the glucosinolates transferred in the flower buds [49]. Qian et al. [58] investigated the effect of light quality on glucosinolate composition and content of Chinese kale sprouts under 23 • C temperature, 80% relative humidity, 16/8 h photoperiod, and red: blue: white light condition, and observed 167.32-288.70 and 72.66-87.48 µmol/g DW of total glucosinolates in shoots and roots, respectively. Similarly, temperature, humidity, and CO 2 have an individual effect on glucosinolate components and accumulation. Rosa and Rodrigues [59] reported that the amount of glucosinolates increases 4-35% in the Brassica species in summer compared to winter seasons. They also observed 386 ± 71 µmol/100 g DW of total glucosinolates in the Chinese cabbage leaves under 20 • C, which increased up to 409 ± 104 µmol/100 g DW under 30 • C. The possible reason behind this increment is the proportional relationship between temperature and the photosynthesis rate. However, glucosinolate components and contents are degraded under both very hot and cold temperatures. Although the light types, intensity, photoperiod, and EC-pH were kept constant in this study (for experiments 1, 2, and 3), the variation of glucosinolate levels was observed due to the individual effect of temperature, humidity, and CO 2 . As the experiment 1, 2, and 3 were conducted separately, the overaccumulation of glucosinolates under temperature treatments might have occurred due to the overall growing condition; however, it is very important to maintain consistency between experiments. To minimize the inconsistency between experiments, the following measurements could be considered to handle and minimize the variations. First, similar seedlings could be prepared as much as possible, so that pre-transplanting cannot affect the final harvested product. Moreover, the number of samples could be increased by cultivating kale in bigger and multiple plant beds. Finally, maintenance of the same cultivation condition through more accurate and precise control of the environmental variables is necessary.
In the correlation matrixes, strong positive correlations were observed among all the physical growth variables, except for the chlorophyll level. Negative correlations with other growth variables were shown. The most likely reason for this phenomenon is that chlorophyll is an indicator of the health of the photosynthetic apparatus, and the concentration (amount per mass) is a function of the leaf area. As the midrib and petiole of kale (depending on the cultivar) are large, the midrib might become enlarged, diluting the concentration of chlorophyll in the lamina during the growth period, which results in negative correlations with other growth variables [60]. Moreover, the efficiency of chlorophyll varies over time due to the engagement-disengagement of assorted photoprotective mechanisms under fluctuating light conditions. This results in energy loss (absorbed by chlorophyll as heat) and affects carbohydrate (glucose) accumulation. This might be another reason for the negative correlation with kale growth [61,62]. Besides this, glucosinolate components were strongly positively correlated with each other under elevated CO 2 concentrations, because glucosinolate synthesis is proportionally related to photosynthesis [38]. The mechanisms of biochemical reactions are, in fact, very complex, and in many processes, the biochemical pathways are only hypothetical or assumed, and the intermediate reactions and products are not fully known. At any stage in the biochemical chain, double bonds, which are very reactive, may be affected by temperature, relative humidity, and CO 2 as well as by free radicals in the environment. Particularly, glucosinolate synthesis can be illustrated according to the following steps: (1) radical substitution and the addition of water occurs at elevated temperatures and/or in the presence of radicals; (2) addition reactions to carbonnitrogen double bonds, resulting in carbonic acid esters; (3) electrophilic addition of water to double bonds that creates two new sigma bonds, resulting in the formation of alcohol; and (4) the occurrence of rearrangement, transposition, and isomerization involving double bonds, allyl radicals, and the glucose cycle. The increased reactivity of double bonds makes them very susceptible to environmental factors, specifically temperature, relative humidity, and CO 2 .
Cartea and Velasco [6] reported that the concentrations of glucosinolate components vary depending on genetics and environmental factors, along with the crop cultivation methods, harvest, storage, and even the processes of meal preparation. Velasco et al. [49] specifically showed that the concentrations of aliphatic glucosinolate components gradually increase in vegetative tissues (i.e., leaves) and are transferred to the reproductive tissues (i.e., flowers and seeds) during the flowering period. In addition, the indole glucosinolate components of leaves and flower buds gradually decrease after a certain period of cultivation. However, the concentrations of aromatic glucosinolates do not vary significantly with the cultivation period. In this study, kale was cultivated in the plant factory using an aeroponic method (one type of hydroponics). The fast growth rate due to proper ambient environment and nutrient management might be a reason for high glucosinolate accumulation in the early stage (two weeks after transplantation), and it gradually declined with the cultivation period (four weeks after transplantation). Determination of the proper harvesting time of brassicaceous plants has been investigated in several studies [63,64]. Based on the environmental factors and cultivation methods used in this study, early harvesting (2-3 weeks after transplantation) is suggested as a possible strategy to achieve glucosinolate-rich kale.

Conclusions
This study was conducted to investigate the effects of temperature, relative humidity, and CO 2 on the growth and glucosinolate content of kale plants hydroponically grown in a plant factory, where five different treatments of each environmental variable were applied separately, and samples were collected after two different periods of cultivation. According to the results, the optimal temperature, relative humidity, and CO 2 range for growth and total glucosinolate content were 20-23 • C, 85%, and 700-1000 ppm, and 14-17 • C, 55-75%, and 1300-1600 ppm, respectively. The glucosinolate content of kale was high in the early growth stage, with low temperature and humidity levels, and elevated CO 2 concentrations. Strong positive correlations were observed among the physical growth variables, and weak correlations were found between the growth and glucosinolate parameters, which indicated that high physical growth might not ensure the high concentration of glucosinolates. According to the findings of this study, early harvesting (i.e., after 2 weeks of transplantation) could be preferred. As the optimum level of temperature, humidity, and CO 2 was different in two-and four-week sampling times, dynamic ambient environment management might be adopted. Farmers could maintain the optimum range of each environmental variable separately based on their target (growth or glucosinolate level), or preferred combined management of the temperature, relative humidity, and CO 2 during kale cultivation within protected cultivation facilities.