Improvement and Prediction of Secondary Metabolites Biosynthesis Under Yeast Extract Elicitation of Neem Cell Suspension Culture Using Response Surface Methodology

Neem is the medicinal plant used as antimalarial, antibacterial, antiviral, insecticide and antimicrobial. The aim of this study was investigate effect of yeast extract and sampling time on cell growth, secondary metabolites synthesis, SQS1 and MOF1 genes expression and prediction of secondary metabolites synthesis by response surface methodology. Highest fresh and dry cell weight were 580.25 g/L and 21.01 g/L obtained 6 days after addition 100 mg/L yeast extract. Highest azadirachtin accumulation and production were 16.08 mg/g DW and 219.78 mg/L obtained at 25 mg/L yeast extract for 2 and 4 days. Maximum mevalonic acid accumulation and production were 1.746 mg/g DW and 23.77 mg/L observed 2 days after addition 50 mg/L yeast extract addition. Highest amount of squalene accumulation (0.215 mg/g DW) and production (4.526 mg/L) were achieved 4 days after application of 50 mg/L yeast extract. Prediction results showed that highest azadirachtin accumulation (13.607 mg/g DW) and production (190.50 mg/L), mevalonic acid accumulation (0.500 mg/g DW) and production (5.57 mg/L) and squalene accumulation (0.301 mg/g DW) obtained by using 245 mg/L yeast extract for 2 days, 71 mg/L yeast extract for 2 days, 200 mg/L yeast extract for 4.96 days, whitout yeast extract for 6.54 days and 4 days, respectively and predicted that highest squalene production is achieved by long-term exposure to high concentrations of yeast extract. qRT-PCR analysis showed the maximum relative gene expression of SQS1 and MOF1 by addition 150 mg/L yeast extract for 4 days and 25 mg/L yeast extract for 2 days, respectively. Ex µl step The qRT-PCR analysis of SQS1 and MOF1 genes performed with two biological and two technical replications. The mean comparison of relative expression of genes carried out using Duncan’s multiple range test at a probability level of 0.01. study of independent including of yeast of ndings. In qRT-PCR analysis the maximum relative gene expression of SQS1 and MOF1 by addition 150 mg/L yeast extract for 4 days and 25 mg/L yeast extract for 2 days, respectively.


Introduction
Plants are a rich source of medicinal compounds used to make medicine. About a quarter of the drugs approved by the FDA and the EMA are produced from plants, which shows their importance (Borkotoky and Banerjee 2020; Thomford et al. 2018). Neem (Azadirachta indica) is a member of meliaceae family and very important in traditional medicine. Studies have shown that neem leaf extract have free radical scavenging and anti-in ammatory activity and inhibits HSV-1 and MHV viruses (Alzohairy 2016;Sarkar et al. 2020;Tiwari et al. 2010). It has been reported that neem leaf extract can increase immunity against HIV/AIDS by increasing CD4 + cell levels (Mbah et al. 2007). Also, it has also been reported that neem extract may be used against the COVID-19 infection (Roy and Bhattacharyya 2020). Neem tree has a variety of secondary metabolites, including azadirachtin, mevalonic acid, squalene, nimbin, nimbiol, polyphenolic avonoids and etc. (Borkotoky and Banerjee 2020;Farjaminezhad and Garoosi 2020).
Secondary metabolites are not necessary to maintain the plant life cycle, but play an important role in environmental adaptation (Park et al. 2020). Secondary metabolites have various uses and are mainly used as drugs, avorings, fragrances, pigments, biopesticides and food additives (Murthy et al. 2014). Studies show that the production of secondary metabolites depends on plant genetics, environmental factors, climate, season, growth period, plant parts, pre-and post-harvest processes, and extraction methods (Açıkgöz 2020). All of these factors create problems in the production of secondary metabolites in the traditional method and increases the cost of production (Ramachandra Rao and Ravishankar 2002). Therefore, plant tissue culture is used to produce secondary metabolites due to its reliability and predictability, independence from geographical, seasonal and environmental factors, modi cation or elimination of unwanted taste and production of high quality and standard product (Abd El-Salam et al. 2015). Cell suspension culture is the best option to increase the production of secondary metabolites and respond to the increasing industrial demand for secondary metabolites (Rani et al. 2020).
There are several strategies to improve the production of secondary metabolites in plant tissue and cell culture (Park et al. 2020). One of these strategies is the use of elicitors. Elicitors are compounds that stimulate the production of secondary metabolites by biochemical changes in the plant (Namdeo 2007). Elicitors do this by activating the signal transductio cascades (Karalija et al. 2020). These signals act as stress factors and regulate enzyme activity and the production of secondary metabolites (Sharma and Agrawal 2018). Elicitors are divided into two types depending on the origin including biotic and abiotic. Elicitors are produced from microbial or plant sources. Yeast extract is one of the elicitors derived from microbial sources (Ramirez-Estrada et al. 2016).
An experimental design is a set of tools for studying a system and includes planning and performing a set of experiments to determine the impact of the studied variables on the system. The obtained valid model is a model that contains information about the effect of experimental conditions on the measurement response rate. The required experiment is performed in such a way that a large amount of information is obtained from a limited number of experiments. Once the appropriate model is obtained, it can be used to predict future observations within the original design range. Therefore, it is appropriate to use the experimental design to develop and optimize a wide range of laboratory studies. The response surface methodology (RSM) was rst used in chemical experiments and includes statistical and mathematical tools for designing and analyzing response surfaces. The experimental design method and the response surface methodology are closely related and the use of the response surface methodology is based on experimental data (Mäkelä 2017). This statistical model is an effective tool for optimizing complex processes and save time during the experimental phase (Menezes Maciel Bindes et al. 2019). In the present study the effect of yeast extract and sampling time investigated on cell growth, azadirachtin, mevalonic acid and squalene accumulation and production and predicted their effects on accumulation and production by response surface methodology. Also, the effect of yeast extract was studied on the squalene synthase 1 (SQS1) and squalene epoxidase 1 (MOF1) genes expression. This is the rst comprehensive study and so far no study has been done on gene expression and prediction on the effect of yeast extract on cell suspension culture of neem.

Plant material and cell suspension culture establishment
The leaves of neem collected from Bandar Abbas city of Iran and surface sterilized with ethanol and sodium hypochlorite and cultured on MS medium containing 1 mg/L picloram and 2 mg/L kinetin and maintained in the growth chamber at 25 ± 2 °C in the dark. The friable calli were transferred to the liquid MS medium with same concentrations of picloram and kinetin and kept on a rotary shaker at 110 rpm with 26 ± 2 °C in the dark and sub-cultured every 12 days (Farjaminezhad and Garoosi 2019).

Elicitation with yeast extract
The cell suspension cultures were transferred to 100 mL Erlenmeyer asks containing 25 mL liquid MS medium supplemented with 1 mg/L picloram and 2 mg/L kinetin with an initial cell density of 2.6 × 10 5 (SCV = 8%). According to the our previous study growth curve (Farjaminezhad and Garoosi 2019), eight days after culture different concentrations of yeast extract including 0, 25, 50, 100, 150 and 200 mg/L were added to cell suspension culture (Merck, Germany). The cultures were kept on a rotary shaker at 110 rpm with 26 ± 2 °C in the dark and sampling were performed after 2, 4, 6, 8, 10 and 12 days.
Fresh and dry cell weight determination Fresh and dry cell weights were measured by Godoy-Hernández and Vázquez-Flota (2006) method with a little modi cation. The cells were collected by Whatman No. 1 lter paper using Büchner funnel under vacuum and retained under vacuum for 30 s and weighed immediately. Then the collected cells transferred to oven with 50 °C temperature and maintained for 72 and weighed immediately for dry cell weight.

Mevalonic acid, squalene and azadirachtin extraction
Azadirachtin, mevalonic acid, and squalene were extracted by Ra q and Dahot (2010) method with a modi cations. One mililier of dichloromethane was added to 100 mg of dried and powdered cells and sonicated for 25 min at room temperature, centrifuged at 7000 rpm for 15 min, and naly supernatant was collected. The procedure was performed in duplicate. The dichloromethane was evaporated at 50 °C in a water bath and samples were dried. The dried samples were redissolved in 1.5 mL HPLC-grade distilled water and maintained at − 20 °C.
RNA extraction, cDNA synthesis and qRT-PCR analysis For RNA extraction, cDNA synthesis and qRT-PCR analysis the samples with highest amount azadirachtin accumulation at each concentrations of yeast extract were used. Total RNA was extracted using RNX-Plus kit (Cinaclon, Iran) based on the producer's guidance. The quantity of extracted RNA was evaluated by Nano-Drop 200C spectrophotometer (Thermo Scienti c, USA). Then, extracted total RNAs were treated with DNase I, RNase-free (Sinaclon, Iran) according to the producer's guidance to eliminate remaining genomic DNA. For single strand cDNA synthesis, the mixture of 5 µg of total RNA, 0.5 µg/µL Oligo (dT) 18 primer (Cinaclon, Iran) and 12.5 μL DEPC-treated water in the tube maintained at 65 °C for 5 min and immediately transferred on ice.

Statistical analysis and experimental design by response surface methodology (RSM)
The treatment of yeast extract was performed in a factorial experiments based on a completely randomized design in triplicate. The rst factor was different concentration of yeast extract and the second factor was different sampling times. Data were analyzed by IBM SPSS Statistics software, Version 24.0 (Armonk, NY, USA).The measured indices means compared by using Duncan's multiple range test at a probability level of 0.01. The qRT-PCR analysis of SQS1 and MOF1 genes performed with two biological and two technical replications. The mean comparison of relative expression of genes carried out using Duncan's multiple range test at a probability level of 0.01.
Response surfece methodology was used to study the effects of independent variables including different concentrations of yeast extract and different sampling time and their interaction on accumulation and production of azadirachtin, mevalonic acid and squalene. The sampling times were selected on the basis of preliminary studies. The central composite design (CCD) with two variables and ve different levels (-2, -1, 0, +1, +2) was used to optimization of yeast extract concentration and sampling time.
Total 13 experiments were conducted to test the ve levels of yeast extract and sampling time with full-factorial CCD. By using coded units, the experimental and predicted values for the azadirachtin, mevalonic acid and squalene accumulation and production in terms of the different variables of yeast extract and sampling time are presented in Tables 2 and 4. The predicted response was calculated using the quadratic polynomial model. The predicted responses were calculated by a second-order polynomial (quadratic) model was shown as follows: Where, Y is the response variable, β 0 is the average response obtained during replicated experiments of the CCD, β i , β ii and β ij are the linear, quadratic and cross-product effects, respectively; X i and X j are the independent coded variables. Response surface regression coe cient and Analysis of Variance (ANOVA) predicted the effects of independent variables on azadirachtin, mevalonic acid and squalene accumulation and production from cell suspension culture of neem. The data were analyzed using Design Expert (12.0.0 version) software.

Fresh and dry cell weight
The results demonstrated that different yeast extract concentrations, sampling time and thier interactions affected both the fresh and dry cell weight ( Table 1). Application of yeast extract decreased the fresh and dry cell weight. The most suitable condition for neem cell suspension growth was control without any concentrations of yeast extract. In this conditions maximized the fresh and dry cell weight,which were 413.41 g/L and 14.47 g/L, respectively. Based on results, by addition 25, 50, 100, 150 and 200 mg/L yeast extract the fresh cell weight reduced 29.61%, 35.61%, 27.73%, 45.52% and 48.66% and dry cell weight reduced 22.53%, 25.57%, 8.64%, 27.99% and 33.86% compared to the control, respectively. The sampling time of 6 and 4 days also gave the highest fresh cell weight and dry cell weight of 410.69 g/L and 16.23 g/L. The fresh cell weight increased from the 2 nd days to the 6 th days of sampling time but decreased from 6 th day to the 12 th days, but the dry cell weight increased from the 2 nd days to the 4 th days of sampling time and decreased from 4 th day to the 12 th days. Therefore, the frech cell weight on the 6 th days of sampling was 47.00%, 1.66%, 88.52%, 111.52% and 104.56% higher than 2 nd , 4 th , 8 th , 10 th and 12 th days of sampling and dry cell weight on 4 th day of sampling was 22.67%, 9.90%, 80.05%, 95.30% and 99.96% higher than 2 nd , 6 th , 8 th , 10 th and 12 th days of sampling. By combination study of the effect of different concentrations of yeast extract and sampling times, it was found the the highest fresh and dry cell weight were 580.25 g/L and 21.01 g/L obtained 6 day after addition 100 mg/L yeast extract and 4 days after using 50 mg/L yeast extract. In this treatments, fresh cell weight increased 9.00%, 50.86%, 23.23%, 109.97% and 146.06% compared to control, 25, 75 and 100 mg/L yeast extract at same sampling time and the dry cell weight increased 24.17%, 34.33%, 21.37%, 40.91% and 81.43% compared to the control, 25, 75 and 100 mg/L yeast extract at same sampling time (Table S1).

Azaidrachtin accumulation and production
HPLC-DAD analysis demonstrated that the applied concentrations of yeast extract, sampling times and their interaction signi cantly stimulated the azadirachtin accumulation and production in treated cells compared to the control (Table 1).
Accumulation and production of azadirachtin showed a dose-dependent response to the yeast extract. The more increase in azadirachtin accumulation and production observed at 25 mg/L yeast extract. Yeast extract with lower cocnetration including 25 mg/L showed the highest azadirachtin accumulation and production which was 9.67 mg/g DW and 118.53 mg/L. The azadirachtin accumulation and production from control treatment to 25 mg/L yeast extract increased and then decreased with increasing yeast extract concentration from 25 to 100 mg/L. At the 25 mg/L yeast extract the azadirachtin accumulation was 161.57%, 19.79%, 42.45%, 36.34% and 59.37% and the azadirachtin production was 119.74%, 33.80%, 31.16%, 51.28% and 106.81% higher than control, 50, 100, 150 and 200 mg/L (Fig. 1a). In terms of sampling time, the highest azadirachtin accumulation and production were 9.20 mg/g DW and 125.65 mg/L observed at 2 days after treatment. On the 2 nd day of sampling, the azadirachtin accumulation increased 17.94%, 42.22%, 75.44%, 58.97% and 34.02% and azadirachtin production increased 0.53%, 31.53%, 160.65%, 207.37% and 142.42% compared to the 2 nd , 4 th , 8 th , 10 th and 12 th days, respectively (Fig. 1b).
The effect of different concentrations of yeast extract along with different sampling times is shown in Table S1. The best conditions for induction of azdiracetin accumulation and production were application of 25 mg/L for 2 days (16.08 mg/g DW) and 4 days (219.78 mg/L), respectively. Under these conditions, the azadirachtin accumulation increased 462.24%, 42.43%, 66.11%, 82.93% and 146.62% compared to the control, 50, 100, 150 and 200 mg/L yeast extract the 2 nd day of sampling and the azadirachtin production increased 302.16%, 54.09%, 96.23%, 51.01% and 191.79% compared to the control, 50, 100, 150 and 200 mg/L at same day, respectively (Table S1).
Model for predicting of azadirachtin accumulation and production The sampling times of 2, 4, 6, 8 and 10 days choosed to RSM analysis and prediction of azadirachtin accumulation and production. The combined effects of different concentrations of yeast extract and sampling times were investigated by the response surface methodology using a central composite design (CCD). The speci c interaction of different concentrations of yeast extract and sampling times with the measured and predicted response values of azadirachtin accumulation and production are shown in Table 2. In this study, the experiment No. 2 in 100 mg/L yeast extract and sampling on the 2 nd day had the highest amount of azadirachtin accumulation (9.89 mg/g DW) and the experiment No. 13 in 150 mg/L yeast extract and sampling on the 8 th day had the lowest azadirachtin accumulation (4.29 mg/g DW). Also, the experiment No. 2 with application of 100 mg/L yeast extract and sampling time of 2 days had the highest azadirachtin production (170.37 mg/L) and the experiment No. 3 with application 50 mg/L yeast extract and sampling time of 8 days had the lowest azadirachtin production (32.44 mg/L).
Analysis of variance (ANOVA) of the results of the response surface model was presented in Table 3. In the analysis of variance of CCD, the coe cient of determination of the model for azadirachtin accumulation and production were 95.76% and 94.08%, respectively; which indicates that 95.76% and 94.08% of the actual levels can correspond to the predicted levels. Also, the p-value of the models were signi cant and the proposed models were appropriate. Therefore, the following formuls were obtained to predict the azadirachtin accumulation and production using yeast extract and sampling times: Azadirachtin accumulation (mg/g DW) = 6.46 + 0.2452 A -1.36 B -1.09 AB -0.4178 A 2 + 0.2269 B 2 Azadirachtin production (mg/L) = 121.32 + 1.97 A -38.14 B + 11.21 AB -17.73 A 2 -3.25 B 2 Optimization of the response surface of azadriachtin accumulation and production The interaction between different yeast extract concentrations and different sampling times is shown in Fig. 1. Increasing the concentration of yeast extract along with increasing the exposure time to 6 days had a signi cant effect on the azadirachtin accumulation. The best yeast extract concentrations and sampling times for maximum azadirachtin accumulation were between 100-200 mg/L and 2-4 days. The highest amount of azadirachtin accumulation in this analysis was obtained with application of 100 mg/L yeast extract and sampling time was 2 days. However, the optimal conditions for maximizing the azadirachtin accumulation (13.607 mg/g DW) predicted two days after application of 245 mg/L yeast extract (Fig 1c and d). Also, the results indicated that azadirachtin production depends on sampling time. According to Fig. 1e and f, at all concentrations of yeast extract, azadirachtin production gradually decreases with increasing sampling time. The highest azadirachtin production is achieved by applying 0-150 mg/L of yeast extract and sampling time of 2-6 days. In this analysis, the highest azadirachtin production was 2 days after application of 100 mg/L yeast extract, but it is predicted that the highest azadirachtin production with 71 mg/L is obtained by culturing of cell for 2 days at medium containing 190.50 mg/L yeast extract.

Mevalonic acid accumulation and production
The mevalonic acid accumulation and production were signi cantly changed under different yeast extract concnetrations, sampling times and their interactions (Table 1). Among the studied concentrations of yeast extract, 50 mg/L increased the accumulation of mevalonic acid compared to the control and 25, 100, 150 and 200 mg/L decreased it compared to the control.
With increasing the concentration of yeast extract in the culture medium from 0 to 25 mg/L, the accumulation of mevalonic acid decreased and then with increasing the concentration of yeast extract to 50 mg increased and with increasing the concentration of yeast extract form 50 mg/L to 200 mg/L its amount decreased. The highest amount of mevalonic acid accumulation was 0.632 mg/g DW obtained at 50 mg/L yeast extract, which was 17.31% higher than the control and 48.71%, 381.59%, 265.15% and 353.86% higher than the 25, 100, 150 and 200 mg/L yeast extract, respectively. Therefore, higher concentrations of yeast extract in the culture medium prevented the mevalonic acid accumulation. With increasing the yeast extract concentration from 0 to 25 mg/L the production of malonic acid decreased and with increasing the concentration of yeast extract from 25 to 50 mg/L increased; however, the difference between control and application of 50 mg/L yeast extract was not statistically signi cant. The highest mevalonic acid production (8.45 mg/L) was obtained at 50 mg/Lyeast extract. By adding 50 mg/L yeast extract to culture medium, the mevalonic acid production increased 8.94%, 47.06%, 381.84%, 302.82% and 424.19% compared to the control and 25, 100, 150 and 200 mg/L yeast extract, respectively (Fig. 2a). Therefore, the application of moderate concentrations of yeast extract had a positive effect on the accumulation and production of mevalonic acid. By investigation of the effect of sampling times, we founded that prolonged exposure of neem cell suspension culture with yeast extract reduces the mevalonic acid accumulation. The highest amount of mevalonic acid accumulation with an average of 0.778 mg/g DW was observed on the second day of sampling. In general, two days after treatment the amount of mevalonic acid accumulation was 55.88%, 79.76%, 213.19%, 915.48% and 5525.94% higher than the 4 th , 6 th , 8 th , 10 th and 12 th days, respectively. Also, between different sampling times the highest mevalonic acid production (9.95 mg/L) was observed on the 2 th day of sampling. By culturing of neem cells for 2 days, the mevalonic acid production was increased 23.63%, 50.01%, 359.17%, 1724.97% and 11214.56% compared to days 4, 6, 8, 10 and 12, respectively (Fig. 2b). The interactions of yeast extract concentrations and sampling times showed that the highest mevalonic accumulation (1.746 mg/g DW) and production (23.77 mg/L) obtained two days after application of 50 mg/L yeast extract. In this conditions, mevalonic acid accumulation was 8.92%, 212.90%, 1954.11%, 334.27% and 2716.13% higher then the control, 25, 100, 150 and 200 mg/L yeast extract and mevalonic acid production was 29.82%, 216.93%, 1801.60%, 370.69% and 525.53% higher than the 25, 50, 75 and 100 mg/L yeast extract on same day (Table S1).

Model for predicting of mevalonic acid accumulation and production
Based on the results, the sampling times of 2, 4, 6, 8 and 10 days selected to RSM analysis and prediction of mevalonic acid accumulation and production. The speci c interaction of different concentrations of yeast extract and sampling times with the measured and predicted response values of mevalonic acid accumulation and production is shown in Table 2. The experiment No.
10 (six days after addition of 200 mg/L yeast extract) had the highest amount of mevalonic acid accumulation (0.508 mg/g DW) and the experiment No. 2 (two days after application of 100 mg/L yeast extract) had the lowest of mevalonic acid accumulation. Also, the experiment No. 12 with application 100 mg/L yeast extract for 6 days had the highest amount of mevalonic acid production (6.97 mg/g DW) and the experiment No. 2 with using 100 mg/L yeast extract for two days produced the lowest of mevalonic acid production. Analysis of variance results showed that the coe cient of determination (R 2 ) of the model for mevalonic acid accumulation and production were 91.67% and 87.53%, respectively; which indicates that 91.67% and 87.53% of the actual levels can correspond to the predicted levels. Also, the p-value of the models were signi cant and the proposed models were appropriate (Table 3). Therefore, the following formuls were obtained to predict the mevalonic acid accumulation and production using yeast extract and sampling times: The interaction between different yeast extract concentrations and different sampling times is shown in Fig. 2. The results showed that yeast extract concentration is the most important factors. Accordingly, gradual increases in the concentration of yeast extract leads to a decrease in the accumulation of mevalonic acid in the cell suspension culture of neem. Therefore, in this analysis the highest amount of accumulation of mevalonic acid was observed at their highest level. Depending on the concentration of yeast extract, sampling time increased or decreased the accumulation of mevalonic acid. According to the Fig. 2c and d, the highest amount of mevalonic acid accumulation was obtained 6 days after applying 200 mg/L of yeast extract, but it is predicted that the highest amount of mevalonic acid accumulation with the amount of 0.500 mg/g DW was obtained by culture of the cells for 4.96 days in medium containing 200 mg/L yeast extract. Also, according to the results of this analysis the production of mevalonic acid depends on the sampling time. At different concentrations of yeast extract, increasing the duration of cell suspension culture exposure to yeast extract increased the amount of mevalonic production from 2 days to 6 days and then decreased. The highest mevalonic acid production obtained by exposing of the cell suspension with 100 mg/L yeast extract for 6 days. The optimal conditions for maximizing the mevalonic acid production predicted at cell suspension culture without yeast extract after 6.54 days, which can produce 5.57 mg/L mevalonic acid.

Squalene accumulation and production
The results showed that different concentrations of yeast extract, sampling times and their intractions had signi cant effect on squalene accumulation and production ( Table 1). The study of squalene accumulation in cell suspension culture in the presence of yeast extract showed that yeast extract inhibited squalene accumulation. Increasing the concentration of yeast extract in the culture medium reduced the accumulation of squalene. With increasing the yeast extract concentration from 25 to 50 mg/L, the squalene accumulation increased and then decreased. The highest squalene accumulation between different concentrations of yeast extract was obtained (0.067 mg/g DW) at control condition. With using 25, 50, 100, 150 and 200 mg/L yeast extract sqaulene accumulation decreased 59.70%, 5.97%, 35.82%, 42.27% and 64.17% than control, respectively. Application of 50 mg/L yeast extract increased squalene production compared to the control and application of 25, 100, 150 and 200 mg/L decreased it.
Among different concentrations of yeast extract, the highest production of squalene (1.132 mg/L) was obtained at 50 mg/L, which was 14.85%, 227.71%, 66.62%, 123.72% and 398.77% higher compared to the control and 25, 100, 150 and 200 mg/L, respectively (Fig. 3a). Between different sampling times, the highest amount of squalene accumulation was 0.098 mg/g DW obtained on the 4 th day of sampling. In general, 4 days after treatment the amount of squalene accumulation compared to 2, 6, 8, 10 and 12 days was 547.87%, 44.95%, 363.76%, 149.67% and 361.09% higher, respectively. Investigation of different sampling times showed that during the rst 4 days, the production of squalene in the cell suspension culture of neem increased and then decreased. The highest amount of squalene production (1.738 mg/L) was observed on the 4 th day of sampling. So, on the 4 th day of sampling the amount of squalene production compared to the 2 nd , 6 th , 8 th , 10 th and 12 th days was 732.07%, 68.63%, 748.53%, 343.38% and 474.63% higher, respectively (Fig. 3b). Further study of squalene accumulation and production was performed using simultaneous examination of different concentrations of yeast extract and sampling times. The results showed that 4 days after application of 50 mg/L yeast extract the highest amount of squalene accumulation (0.215 mg/g DW) and production (4.526 mg/L) obtained, which were 283.93%, 1333.3%, 69.29%, 41.45% and 923.81% and 377.43%, 1817.79%, 105.54%, 99.47% and 1739.84% higher compared to the control, 25, 100, 150 and 200 mg/L, respectively (Table S1).

Model for predicting of squalene accumulation and production
The sampling times of 4, 6, 8, 10 and 12 days selected to RSM analysis and prediction of squalene accumulation and production.
The experimental and predicted values of yeast extract-induced squalene accumulation are shown in Table 4. The experiment No. 6 (six days after application of 50 mg/L yeast extract) had the highest squalene accumulation (0.1351 mg/g DW) and the experiment No. 10 (150 mg/L yeast extract after 10 days) had the lowest squalene production (0.012 mg/g DW). The experiment No. 3 with application of 100 mg /L yeast extract and sampling time of 4 days had maximum squalene production (2.21 mg/L) and the experiment No. 13 with application of 50 mg/L yeast extract and sample time of 10 days had the lowest squalene production. The results of analysis of variance of central composite design showed that the coe cient of determination (R 2 ) of the model is 94.20% for squalene accumulation and 94.48% for squalene production. This indicates that 94.20% and 94.48% of the actual levels of squalene accumulation and production correspond to the predicted levels. Also, the p-value of the models were signi cant for squalene accumulation and production, which indicates the suitability of the models (Table 3). Therefore, mathematical models were obtained to predict the squalene accumulation and production under yeast extract elicitation: Squalene accumulation (mg/g DW) = 0.0199 -0.0169 A -0.0304 B + 0.0274 AB + 0.0058 A 2 + 0.0135 B 2 Squalene production (mg/L) = 0.2150 -0.2509 A -0.5438 B + 0.4834 AB + 0.0725 A 2 + 0.2466 B 2 Optimization of the response surface of squalene accumulation and production Therefore, depending on the concentration of yeast extract, sampling time can have a positive or negative effect on squalene accumulation. Based on this analysis results, the highest amount of squalene accumulation was obtained 6 days after application of 50 mg/L yeast extract, but the prediction results showed that highest squalene accumulation obtained at free-yeast extract medium after 4 days, which can accumulate 0.301 mg/g DW of squalene. Also, increasing the yeast extract concentration during low exposure periods reduced squalene production, but increasing yeast extract concentration at high exposure times enhanced squalene production. According to Fig. 3e and f, it can be said that using 0-100 mg/L yeast extract for 4-8 days can produce an acceptable amount of squalene. It is predicted that the highest amount of squalene production is achieved by longterm exposure to high concentrations of yeast extract.
qRT-PCR analysis of SQS1 and MOF1 genes expression qRT-PCR analysis showed relative gene expression of SQS1 was signi cantly increased 68.96% at 4 days application of 150 mg/L yeast extract compared with the control cells after 12 days. But, after addition of 25 mg/L yeast extract for 2 days, 50 mg/L yeast extract for 2 days, 100 mg/L yeast extract for 2 days and 200 mg/L yeast extract for 12 days, the relative expression of SQS1 gene decreased 73. 26%, 74.49%, 38.29% and 73.53% compared to the control after 12 day. Also, the application of yeast extract signi cantly up-regulated MOF1 gene. The highest relative expression of MOF1 gene obtained by addition 25 mg/L yeast extract for 2 days, which was 52.41%, 14.46%, 26.02%, 22.51% and 39.90% higher than control at 12 days, 50 mg/L yeast extract at 2 days, 100 mg/L yeast extract at 2 days, 150 mg/L yeast extract at 4 days and 200 mg/L yeast extract at 12 days (Fig. 4).

Discussion
Many researchers have been reported that the use of elicitor is effective in cell growth and increasing bioactive compounds. The age of the cell culture, the duration of exposure to the elicitor and the type of elicitor have a very effective role in increasing the effect of the elicitors (Açıkgöz 2020;Nazir et al. 2019). Elicitor concentration is important factor for elicitation process and optimum level depend on the plant species (Vasconsuelo and Boland 2007). The studies showed that high concentrations of elicitors triggers the hypersensitive response and leading to cell death (Park et al. 2020). Yeast exreact is applied to enhance growth and secondary metabolites production in cell and hairy root cultures of different plants such as Artemisia annua (Putalun et al. 2007) and Salvia miltiorrhiza (Yan et al. 2006). In this study, the highest fresh and dry cell weight obtained 6 day after addition 100 mg/L yeast extract. The positive effect of using yeast extract on cell and plant growth and increasing biomass has been reported in studies. Krstić-Milošević et al. (2017) reported that biotic eleicitors such as yeast extract improved root growth and biomass production in Gentiana dinarica. Bayraktar et al. (2016) enhanced biomass production of Stevia rebaudiana by yeast extract.
The highest azadirachtin accumulation and production were obtained at 25 mg/L yeast extract for 2 and 4 deys, maximum mevalonic acid accumulation and production observed 2 days after addition 50 mg/L yeast extract addition and highest amount of squalene accumulation and production were achieved 4 days after application of 50 mg/L yeast extract. The stimulating effect of yeast extract on the iso avonoid production in Pueraria candollei hairy root cultures was reported (Udomsuk et al. 2011). Rani et al. (2020) showed that utility of yeast extract enhance diadzein and genistein production in cell culture of Pueraria candollei. Sharma and Agrawal (2018) produced 544.60 µg/g DW plumbagin in root cullus of plumbago zeylanica by using 100 mg/L yeast extract. Elicitation of Helianthus tuberosus L. with 0.25 mg/L yeast extract increased inulin content 1.18-fold compared to control (Ma et al. 2017). Park et al. (2016) observed the highest amount of rosmarinic acid at 4.98 mg/g by using 750 mg/L yeast extract in Agastache rugosa cell culture. Szopa et al. (2018) increased lignan production in microshoot cultures of Schisandra chinensis by elicitation with 5000 mg/L yeast extract on the rst day of the growth period and with 1000 and 3000 mg/L on the 20th day.
Elicitors activate the plant's defense response, causing consecutive cellular and molecular events and activating biosynthetic genes involved in the production of secondary metabolites (Jiao et al. 2016). Various studies have shown the upregulated accumulation of different secondary metabolites by different elicitors such as yeast extract (van der Heijden et al. 1989;Yoon et al. 2000). For example, in the hairy root cultures of Artemisia annua L. application fungal elicitors enhanced artemisinin production and upregulate expression pro le of mevalonate and methylerythritol phosphate biosynthetic genes (Ahlawat et al. 2014). These results are consistent with our ndings. In this study, qRT-PCR analysis showed the maximum relative gene expression of SQS1 and MOF1 by addition 150 mg/L yeast extract for 4 days and 25 mg/L yeast extract for 2 days, respectively. In cell cutlures of Uncaria tomentosa and Tobernamontana divaricate, the activities of enzymes involved in tritepenoid biosynthesis such as IDI and SS stimulated by elicitors (Flores-Sánchez et al. 2002;Fulton et al. 1994). Ge and Wu (2005) reported that yeast elicitor stimulated HMGR activity in Salvia miltiorrhiza hairy root cultures. In cell culture of A. rugosa the transcript levels of HPPR under yeast extract treatment were higher than the control (Park et al. 2016 Table 3 Analysis of variance of azadirachtin, mevalonic acid and squalene accumulation and production in the central composite design under the in uence of yeast extract and sampling time.