Morphological, Physiological and Quality Performances of Basil Cultivars under Different Fertilization Types

: The prospect of replacing traditional chemical fertilization with organic and microorganism-based fertilization meets the current demand for more sustainable cropping systems and healthy food. In this respect, research was carried out to evaluate the effects of the factorial combination between four basil cultivars (‘Aromat de Buzau’, ‘Macedon’, ‘Cuisoare’ and ‘Seraﬁm’) and three types of fertilization, namely chemical fertilization (with a solid chemical fertilizer), organic fertilization (with chicken manure formulate) and microorganisms’ fertilization (with microorganisms formulate), on basil yield, biochemical and physiological parameters and essential oil composition. The results showed that the biometric parameters (plant height, number of stems and leaves and leaf area) were signiﬁcantly inﬂuenced by the cultivar; ‘Macedon’ obtained the highest values of plant height (64.7 cm) and number of stems (20.33) and leaves (618.3) and ‘Seraﬁm’ the largest leaf area (4901.7 cm 2 per plant), while the type of fertilization did not affect these parameters. Regarding the biomass, the inﬂuence of the cultivar was not signiﬁcant on fresh biomass but was signiﬁcant on dry biomass, with ‘Macedon’ showing the highest value (56.4 g · plant − 1 dry biomass). The mentioned parameters were signiﬁcantly inﬂuenced by the type of fertilization, with the highest values recorded with chemical fertilization. Both the cultivar and the fertilization type signiﬁcantly inﬂuenced the physiological parameters (the total content of assimilatory pigments and photosynthesis). Five phenolic compounds were quantiﬁed from leaf extracts by HPLC-MS (caffeic acid, hyperoside, isoquercitrin, rutin and quercitrin). Hyperoside was identiﬁed only in ‘Macedon’, while the rest of the compounds were found in all the cultivars and varied depending on the cultivar and fertilization type. Regarding the composition of the essential oil, variation was found depending on the cultivar and fertilization type. In ‘Aromat de Buzau’, the main compounds were methyl chavicol and β -linalool; in ‘Macedon’, geranial and neral; and in ‘Cuisoare’ and ‘Seraﬁm’, β -linalool. Moreover, the PCA showed that the ‘Seraﬁm’ cultivar has exclusive properties compared to the other cultivars. Our results highlight that identifying the most effective interaction between genotype and fertilization type allows to optimize yield and quality targets for sweet basil


Introduction
Sweet basil (Ocimum basilicum L.) is a medicinal, culinary and ornamental species of tropical origin that is characterized by a high ecological plasticity, being cultivated worldwide [1,2].This species can be grown successfully both in open fields and in greenhouses, as well as in pots indoors [3].Green-leafed cultivars, i.e., 'Genovese', are used in the preparation of pesto, a typical green sauce in Italian cuisine [4].Moreover, basil can be consumed fresh as a salad or dry in the preparation of some Mediterranean dishes and drinks [5].Due to its popularity, a number of cultivars adapted to various local conditions with different phenotypes and chemotypes have been introduced in the market [6,7].The main phenotypic characteristics varying among the different cultivars are plant height, leaf shape and color and flower color [8].In a general review on basil, Simon et al. [8] described the plant characteristics of basil cultivars belonging to the Ocimum basilicum species that are found on the North American market and showed that plant height can vary from 29 cm ('Green Ruffles') up to 65-70 cm ('Sweet Dani').Moreover, the color of leaves can be green ('Genovese'), green-purple ('Anise') or purple ('Red Rubin Purple Leaf'), while the color of the flowers can be white ('Genovese'), pink ('Dark Opal') or bright purple ('Purple Ruffles').Other characteristics that can vary are spread (cm), stem and spike color and the number of days to flowering [8].Basil is rich in essential oil, and its composition determines its specific aroma and chemotype [9,10].The most common chemotype found within the European market is considered to have the best aroma and quality due to the high content of linalool and methyl chavicol.Other chemotypes are those from 'Reunion', with a high content of methyl chavicol; from tropics, which have a high content of methyl cinnamate; and from Eastern Europe, Russia and many parts of Asia and North Africa, which have a high content of eugenol [11,12].The content and composition of phenolic compounds and essential oil depend on the cultivar and the cultivation technology.Zheljazkov et al. [13] evaluated the essential oil content depending on the cultivar and the growing location and found that the 'Mesten' cultivar had 0.067% at Beaumont, Mississippi and 0.481% at Verona, Mississippi, while the 'German' cultivar had 0.236% at Beaumont and 0.389% at Verona.Regarding the cultivation technology, Baczek et al. [14] found that the content of linalool was higher in plants grown in an open field compared to those grown in a polytunnel.Basil is also rich in phenolic compounds such as rosmarinic acid, chicoric acid and caffeic acid, which give it bioactivities such as antioxidant, antimicrobial or insecticidal activity [15,16].
Fertilization with organic fertilizers has become more and more attractive for farmers as consumer are willing to pay premium price for organic produce.Indeed, in the EU alone, the land area under certified organic management has increased from 9.5 million hectares in 2012 to 14.7 million hectares in 2020 [17].Microorganism-based products (i.e., arbuscular mycorrhizal fungi and plant growth-promoting rhizobacteria) can be used in organic agriculture to stimulate growth and control pests [18].Moreover, due to the complex mechanisms of action, such as increasing the bioavailability of nutrients by solubilizing macronutrients such as phosphorus or inducing systemic resistance in plants, they can influence the synthesis of compounds with a defense role such as polyphenols [19,20].Currently, the interest in healthy food rich in bioactive compounds has increased, and in this respect, farming management allows for improving the quality of products by increasing the content of these compounds.Organic fertilization may stimulate beneficial microorganisms and subsequently could stimulate the synthesis and accumulation of bioactive compounds [21][22][23], while chemical fertilization can have the opposite effect [24,25].Moreover, in conventional cultivation systems, the use of chemical fertilizers and pesticides can cause the accumulation of some chemical residues both in soil and plant products, in contrast to certified organic crops where the use of chemical fertilizers and pesticides is prohibited [26].Caruso et al. [21], in a comparative study, showed that in sweet pepper (Capsicum annuum L. ssp.annuum), the total content of polyphenols increased with microorganism-enriched conventional fertilization and the total flavonoid content increased with microorganism-enriched organic fertilization compared to conventional fertilization.
Due to the shortage of literature reports about the prospect of replacing traditional chemical fertilization with organic and microorganism-based fertilization types which better meet the current demands for more sustainable crop systems and healthy food, the present investigation aimed to assess the interaction effect between cultivar and fertilization type on the yield, biochemical and physiological parameters and essential oil composition of basil.

Experimental Site
This research was carried out on sweet basil (Ocimum basilicum L.) in 2019 and 2020 at the experimental field of the Didactic and Experimental Station V. Adamachi within the Iasi University of Life Sciences, Romania.The soil was anthropic chambic chernozem with the following characteristics: 61% sand, 33% clay and 6% silt; pH 7.1; EC 495 µS•cm −1 ; 2.79% organic matter; 2.8 g•kg −1 N, 32 mg•kg −1 available P (Olsen method), 218 mg•kg −1 available K (ammonium acetate method) and 4.1 g•kg −1 CaCO 3 ; C/N 5.93.The main meteorological conditions during the research are presented in Table 1.The experiment was established in mid-April by direct sowing in polystyrene multicell trays, with 31.3 cm 3 alveoli.The seedlings were planted in mid-May in the open field, spaced 15 cm to 45 cm in row/between row spacing, resulting in a density of 14.8 plants per square meter.
During the cropping season, the following practices were performed: drip irrigation; manual weeding twice; harvesting was performed at the beginning of flowering (BBCH 61) [27], at the beginning of August.No phytosanitary treatments were necessary [28].

Experimental Design
The experimental protocol was based on the combination of two factors, and a split plot design was arranged for the treatment distribution in the field, with three repetitions: factor A was the cultivar [a 1 -'Aromat de Buzau' (AB); a 2 -'Macedon' (M), a 3 -'Cuisoare' (C); a 4 -'Serafim' (S)], assigned to the plots; factor B was the type of fertilization [b 1 -chemical (Ch); b 2 -organic (O); b 3 -microorganisms (Mo)] assigned to the sub-plots.
The cultivars used in the experiment were developed and commercially propagated in Romania, at the Buzau Vegetable Research Development Station; three of them have already been homologated (AB, M and S), while the fourth is a genotype under test, with the aroma of clove oil (Syzygium aromaticum L.) (C).These cultivars have different morphological and phytochemical characteristics: Aromat de Buzau is a cultivar with green leaves and white flowers, and methyl chavicol and linalool are the main essential oil constituents; Serafim has red leaves and pink flowers, and linalool and eugenol are the main essential oil constituents; Macedon has green leaves and white flowers, and geranial and neral are the main essential oil constituents; Cuisoare has green leaves and purple flowers, and linalool and eugenol are the main essential oil constituents [18].
The dose of N active substance (a.s.) per hectare from the organic fertilization with Orgevit represented approximately 70% of the dose of N a.s.per ha associated with the chemical fertilization with Cristaland, because it was taken into account that the N anion from the oxidation of urea and ammonium is not adsorbed by the surface of the soil colloids, and thus, a leaching loss of a N-NO 3 − fraction is expected.The application of beneficial microorganisms served to evaluate their potential in stimulating plant nutrient absorption in the absence of fertilization [21].

Biometric and Agroproductivity Characteristics Determination
In order to determine the biometric characteristics, the height of the plants was evaluated by measuring them with a ruler and expressing the values obtained in cm, followed by determinations regarding the number of lateral stems per plant [29].
To assess the number of leaves per plant and the leaf surface per plant (cm 2 •plant −1 ), the basil plants were harvested by cutting them 5 cm above the ground.The leaf area index (LAI) was determined using the Li-3100 Area Meter, (LICOR, inc., Lincoln, NE, USA) [30].
The amount of fresh biomass (leaves and stems), expressed in grams per plant, was determined immediately after harvesting by weighing with a Kern analytical balance, with a precision of 0.01 g.The amount of dry biomass was determined after drying the plants in a sheltered, naturally ventilated place for 30 days [31].

Total Chlorophyll Content Determination
The total chlorophyll content was determined with a CCM-200 plus non-destructive portable chlorophyll content meter (Opti-Sciences, ADC BioScientific Ltd., Hoddesdon, Hertfordshire, UK); the recorded values were expressed as CCI units (Chlorophyll Content Index).The measurements were taken one day before harvest.For each experimental treatment, 30 readings with 20 plants were performed.Fully developed leaves at the middle plant height were selected [32].

Photosynthesis Determination
Photosynthesis was determined using an LCi system (ADC Bioscientific UK Ltd., Hoddesdon, Hertfordshire, UK).A broad leaf chamber (6.4 cm 2 ) was used, and the measurements were performed between 9 and 11 a.m.The results were expressed as µmol m −2 s −1 [32].

The Color of Leaves
The color parameters of leaves (L, a and b) were assessed using a MiniScan XE Plus color meter (HunterLab, Reston, VA, USA).The value of L indicates lightness, a indicates the degree of red (+a) or green (−a) and b denotes yellow (+b) or blue (−b) color of leaves.C is the chroma [33].

Extraction and Determination of Phenolic Compounds
The phenolic compounds (caffeic acid, hyperoside, isoquercitrin, rutin and quercitrin) were determined from 10% leaf extract in 70% ethanol by ultrasonication for 30 min at room temperature.An Agilent 1100 HPLC system by Agilent Technologies Inc., Santa Clara, CA, USA, and a Zorbax SB-C18 column were used according to the method described by Mocan et al. [34].

Extraction and Analysis of the Essential Oil Composition
The essential oil was extracted from fresh whole aboveground plant material by steam distillation for three hours and the results were expressed as %.
A GC/FID-GC/MS system (Agilent 5975C MSD coupled to Agilent 7890A GC by Agilent Technologies Inc., Santa Clara, CA, USA) was used to analyze the composition of the essential oil.The complete method was described by Teliban et al. [18], Burducea et al. [35] and Adams [36].

Statistical Analysis
The results were reported as means ± standard errors of the two-year experiment (2019 and 2020), after raw data processing by ANOVA and mean separation through the Duncan multiple range test (p < 0.05) using SPSS v21 software (IBM Corp, Armonk, NY, USA).The principal component analysis (PCA) was performed using OriginPro 2020 Academic by OriginLab Corporation, Northampton, MA, USA.This analysis aimed to reduce the number of variables to determine correlations and interactions between different inputs [37].

Results
The effect of the cultivar and the fertilization type on the biometric parameters is shown in Table 2. 'Macedon' had the highest values of plant height (+65% compared to 'Serafim'), ramifications (+56% compared to 'Aromat de Buzau') and number of leaves (+96% compared to 'Cuisoare'), while 'Serafim' displayed the highest values of leaf area (+44% compared to 'Aromat de Buzau').The fertilization type did not induce significant differences with reference to the mentioned parameters (Table 2).The effect of the cultivar and the fertilization type on biomass parameters is shown in Table 3. Regarding the fresh biomass, 'Serafim' had the highest value of leaves weight (+45% compared to 'Aromat de Buzau'), and 'Macedon' had the highest weight of stems (+68% compared to 'Serafim') and of total plant (+14% compared to 'Serafim').Regarding the dry biomass, 'Cuisoare' had the highest leaves weight (+47% compared to 'Aromat de Buzau'), and 'Macedon' had the highest weight of stems (+162% compared to 'Serafim') and of total plant (+60% compared to 'Serafim').Chemical fertilization showed the highest significant influence on the fresh and dry weights of stems and total plant.From the significant interaction between cultivar and fertilization type, it arose that the highest values of plant height and leaf number were recorded in 'Macedon' under the chemical fertilization, the highest ramification number was in organically fertilized 'Macedon' and the largest leaf area was in 'Serafim' supplied with organic fertilizer (Table 4).From the significant interaction between cultivar and fertilization type on fresh biomass (Table 5), it can be observed that the chemical fertilization led to the highest fresh and dry biomass of leaves in 'Cuisoare' and of stems and total plant in 'Macedon'.
The effect of cultivar and fertilization type on physiological and color parameters is shown in Table 6.The highest content of assimilatory pigments, expressed in CCI (Chlorophyll Content Index), was recorded in 'Serafim' (35.63 CCI), 161% higher than that in 'Aromat de Buzau' which had the lowest value (13.6 CCI).The highest value of photosynthesis was detected in 'Aromat de Buzau', 202% higher than that in 'Serafim'.The chemical fertilization elicited the highest content of assimilatory pigments, 10% higher than the microorganism formulation, with the latter leading to the highest value of photosynthesis, 31% higher than that with chemical fertilization.Among the color parameters (L lightness-darkness, a redness-greenness and b yellowness-blueness), L was significantly higher in 'Aromat de Buzau' and under microorganism fertilization; the highest value of a and the lowest value of b were recorded in 'Serafim', a red-leafed basil cultivar.From the significant interaction between cultivar and fertilization type on the content of assimilatory pigments (Table 7), it arose that the highest value was recorded in 'Serafim' under the chemical fertilization, 176% higher than that in 'Aromat de Buzau' fertilized with microorganisms, and the highest value of photosynthesis was recorded in 'Aromat de Buzau' supplied with the microorganism formulation, 357% higher compared to the chemically fertilized Serafim.
To highlight the correlations and interactions between the experimental factors and the variables examined, a PCA was performed.The two principal components shown in the biplot graph (Figure 2) overall contributed to 72.7% of the total variability (48.48% and 24.22% for PC1 and PC2, respectively).
The extracted eigenvectors' values highlight that both the cultivars and fertilization types are based on PC1, in different ways.The cultivars 'Aromat de Buzau' and 'Macedon' have positive values, while 'Serafim' has a negative value and 'Cuisoare' is near the origin (Figure 2).
The results showed that the cultivar 'Serafim' was closely connected with the leaf area index, and 'Macedon' was connected to the other morphological variables analyzed.The leaves' weight was influenced by chemical fertilization in the cultivars 'Cuisoare' and 'Serafim' and by microorganism fertilization in 'Cuisoare'.The other yield parameters were affected by both chemical and organic fertilizations in the cultivar 'Macedon'.As shown in Figure 1, the basil cultivar 'Macedon' under microorganism treatment had the highest oil content (0.22%), though not significantly different from the chemical and organic fertilization, while the organically fertilized 'Serafim' accumulated the lowest oil amount (0.07%).Generally, the organic fertilization resulted in a lower oil content compared to the other two fertilization types.To highlight the correlations and interactions between the experimental factors and the variables examined, a PCA was performed.The two principal components shown in the biplot graph (Figure 2) overall contributed to 72.7% of the total variability (48.48% and 24.22% for PC1 and PC2, respectively).The leaves of 'Serafim' depend on the Chlorophyll Content Index (CCI), especially in the red band, and 'Aromat de Buzau' is sensitive to the photosynthesis process (Figure 2).
The effects of the fertilization types on the evaluated variables were analyzed separately, and the results showed that the chemical fertilization influenced the CCI, number of leaves, plant height and all yield parameters, whereas fertilization with microorganisms increased the number of stems.
The PCA identified that phenolic compounds were specifically connected with the cultivar; hyperoside was found in 'Macedon', quercitrin in 'Aromat de Buzau' and rutin, isoquercitrin and quercitrin in 'Cuisoare'.
The statistical analysis indicated that the data obtained for all cultivars were clustered on the 2D PC diagram, which proves that the analyzed variables did not depend on the fertilization type, but only on the cultivar.
The cultivar 'Serafim' showed different characteristics compared to the other cultivars, which suggests that it is a cultivar with exclusive properties.These differences could be explained by the fact that 'Serafim' is a cultivar with purple leaves, thus having a higher pigment content due to the presence of anthocyanins compared to the cultivars with green leaves, where they were not detected.The correlation matrix between the analyzed variables of basil supports the mentioned findings, and the complete description of these variables is presented in the supplementary materials (Supplementary Table S1).
In the cultivar 'Aromat de Buzau', 30 compounds were identified in the essential oil (Table 9).Methyl-chavicol was found in the largest amount with values between 42.99% under the chemical fertilization up to 49.29% in the case of microorganism fertilization.The next compound found in large amount was β-linalool, with values between 13.07% with the chemical fertilization and 25.16% with the organic one.Other important compounds detected were β-elemene, germacrene D and epi-α-cadinol.RIcalc-calculated Kovats index; RIlit-Kovats Index by literature data [27]; tr-traces.
In the cultivar 'Macedon', 26 compounds were identified in the essential oil (Table 10).Geranial was found in the largest amount with values ranging from 26.19% in the case of chemical fertilization up to 32.20% under organic fertilization.The next compound found in large amounts was neral, with values between 20.52% with the chemical fertilization and 25.94% with the organic one.Other compounds detected in remarkable amounts were nerol, β-caryophyllene and (E)-γ-bisabolene.In the cultivar 'Cuisoare', 36 compounds were identified in the essential oil (Table 11).β-Linalool was found in the largest amount with values from 30.42% when applying microorganisms up to 40.17% in the case of chemical fertilization.The next important compound was epi-α-cadinol, with values between 9.92% with the chemical fertilization and 13.52% under the microorganism formulation.Other important compounds were eugenol, α-trans-bergamotene, γ-cadinene and germacrene D. RIcalc-calculated Kovats index; RIlit-Kovats index by literature data [27]; tr-traces.
In the cultivar 'Serafim', 30 compounds were identified in the essential oil (Table 12).β-Linalool was found in the largest amount with values from 49.52% when applying the microorganism treatment up to 60.80% in the case of organic fertilization.The next important compound was eugenol, with values between 6.81% under the organic fertilization and 10.37% with the microorganism application.Other main compounds were β-elemene, germacrene D, camphor, α-trans-bergamotene, γ-cadinene and α-guaiene.

Discussion
Currently, the consumer interest towards healthy foods rich in bioactive compounds has increased [38][39][40].A strategy to increase these compounds is the application of farming management able to ensure a balance between the quantity and quality of agricultural products [41,42].The aim of this study was to evaluate the effect of the interaction between cultivar and fertilization type on the morphology, physiology and synthesis of bioactive compounds in basil cultivated in the field.To this end, four basil cultivars ('Aromat de Buzau', 'Macedon', 'Cuisoare' and 'Serafim') were evaluated in combination with three types of fertilization, i.e., chemical fertilization (with a solid chemical fertilizer), organic fertilization (with a chicken manure formulate) and microorganism fertilization (with a microorganisms formulate).As expected, the morphological parameters (plant height, number of stems and leaves and leaf area) were significantly influenced by the cultivar, with 'Macedon' showing the highest plant height (64.67 cm) and number of stems (20.33) and leaves (618.33) and 'Serafim' showing the largest leaf area (4901.67 cm 2 per plant) and the smallest height (39.00 cm).This is due to the great diversity among the existing phenotypically different basil cultivars.For example, Svecova and Neugebauerová [43] investigated 34 cultivars of basil and showed that plant height varied from 143 to 570 mm, while Juske-viciene et al. [44], analyzing ten cultivars of basil, showed that height ranged from 44.0 to 77.6 cm in the greenhouse and from 37.2 to 63.4 cm in the open field.In this study, the fertilization type did not affect the biometric characteristics but significantly influenced the fresh and dry biomass, which attained the highest values with the chemical fertilization.Basil reacts positively to both organic and chemical fertilization.Matlok et al. [45] showed that both plant height and biomass were higher in the Genovese and Violetto cultivars grown on a substrate containing neutral peat (70%), extract of common nettle (10%), horse manure (20%) and organic controlled-release fertilizer Bioilsa N 12.5 compared to plants grown on peat (100%) and under mineral fertilization with ammonium nitrate as a result of the higher nutrient content (N, P, K and Mg) available from horse manure.In a comparison between two basil cultivars, Burducea et al. [46] found that the values of yield and morphological parameters were the highest with chemical fertilization (chemical > AMF > organic > 40 t ha −1 biosolids > 20 t ha −1 biosolids > control).In the present study, the physiological parameters, the total content of pigments and photosynthesis were significantly affected by the cultivar, thus confirming the results of previous research [13].The pigment content was stimulated by chemical fertilization and photosynthesis was stimulated by microorganism fertilization.Similarly, photosynthesis and other associated parameters (stomatal conductance and water use efficiency) increased in Corylus avellana after inoculation with AMF (Trichoderma harzianum and Glomus intraradices) [47].
In addition to its use as an aromatic spice within the food and beverage industries, basil is also known as a medicinal plant due to its antimicrobial, antiseptic, antioxidant and anti-inflammatory effects [48].The chemical compounds which make basil a valuable plant are phenolic compounds, mainly caffeic and rosmarinic acid, rutin and isoquercitrin, and essential oil constituents such as linalool and methyl chavicol [49,50].Many factors affect the content and composition of phenolic compounds and essential oil profile, such as the cultivar, climate, season, sampling period or plant part used for extraction [51][52][53].The fertilization type can also influence the phenolic or essential oil composition [54,55].For instance, in this study, it was observed that depending on the fertilizer used, the content of each phenolic compound varied.In general, fertilization based on microorganisms enhanced all the phenolic compounds analyzed.For example, the highest value of caffeic acid was recorded in the cultivar 'Cuisoare', and hyperoside and isoquercitrin were accumulated to a more remarkable extent by 'Macedon' and 'Cuisoare' and rutin by 'Cuisoare', while quercetin was accumulated more by 'Aromat de Buzau'.The differences between the cultivars with regard to the phenolic compounds analyzed upon microorganism fertilization suggest the important influence of the genotype on the phenolic profile.The same observation was made by Cruz et al. [56] in a study regarding three basil cultivars and the effect of nitrogen input on different parameters, including phenolic compounds.Additionally, the influence of the cultivar on the synthesis of different phenolic compounds, regardless of the fertilizer applied, was observed for hyperoside, which was detected only in trace amounts in the cultivars 'Aromat de Buzau', 'Cuisoare' and 'Serafim'.
As in the case of phenolic compounds, the fertilization type influenced the qualitative and quantitative composition of the essential oil.For instance, some components were produced only when specific fertilizer types were applied: cis/trans-muurola-3,5-diene in 'Aromat de Buzau' with chemical fertilization; cis-thujone and geranyl acetate in 'Macedon' with organic fertilization; trans-muurola-3,5-diene in 'Cuisoare' with microorganism fertilization; β-myrcene, cis-β-ocimene and terpinolene in 'Serafim' with chemical and organic fertilization (Tables 9-12).As for the essential oil composition, there can be variations depending on the fertilization type but also on the cultivar; β-linalool was produced to the highest extent under organic fertilization in the cultivars 'Aromat de Buzau' and 'Serafim' and with microorganism fertilization in 'Macedon' or chemical fertilization in 'Cuisoare'.Moreover, by analyzing the essential oil composition, it was observed that the main components differed depending on the cultivar; β-linalool and methyl chavicol in 'Aromat de Buzau', neral and geranial in 'Macedon' and β-linalool in 'Cuisoare' and 'Serafim'.The highest values of most of the main components, regardless of the cultivar, were obtained under organic and microorganism fertilization, except β-linalool in 'Cuisoare' which showed the highest content with chemical fertilization.Knowledge of the qualitative and quantitative composition, depending on the fertilization and cultivar, is a very important aspect to obtain essential oil rich in specific important components for different medical purposes.For example, methyl chavicol, which has antioxidant and anti-lipase activities [57], was best produced by the cultivar 'Aromat de Buzau' under microorganism fertilization; β-linalool, with antimicrobial (e.g., Candida albicans, Staphylococcus aureus and Escherichia coli), antioxidant, anti-inflammatory and anticancer activities [58], had the highest content in 'Cuisoare' under chemical fertilization; and neral was only synthesized by 'Macedon', with the most remarkable production under organic fertilization.
The PCA and, more specifically, the eigenvector values related to this study revealed that both cultivar and fertilization are based on PC1, which produced the greatest effect on the parameters examined.PCA is a powerful statistical technique that can highlight, for example, the influence of different fertilization types on plants, as was shown in the case of basil fertilized with biosolids [19] or Chinese chives (Allium tuberm) under the action of slow-release fertilizer [59].The fertilization type-for example, chemical or organic (manure-based)-directly influences the microorganism communities in the substrate and the enzymatic activity in the soil with the role of plant growth stimulation or protection [60].On the other hand, 'Serafim' showed completely different results compared to the other cultivars, which suggests that it is a cultivar with exclusive properties.'Serafim' is a purpleleafed cultivar, which makes the pigment content higher than that in green cultivars due to the additional presence of anthocyanin compounds.For example, Šamec et al. [61], through the PCA of the physical, chemical and phytochemical parameters of four cultivars of strawberry, were able to highlight specific cultivar properties by grouping in the left side of the PCA plot the color parameters L* and C* and in the right side the polyphenolic compounds, which indicated that the cultivars with a higher polyphenolic content are darker and more colorful.

Conclusions
In this study, the four basil cultivars examined showed different biometrics and growth parameters in terms of plant height, number of stems and leaves, leaf area and dry biomass, whereas the fertilization type only affected the fresh and dry biomass, with the highest amounts obtained with chemical fertilization.Either the cultivar or the fertilization type significantly influenced the physiological parameters, such as the total content of assimilatory pigments and photosynthesis, the phenolic compounds investigated (caffeic acid, hyperoside, isoquercitrin, rutin and quercitrin) and the essential oil composition.Fertilization with microorganisms led to the production of beneficial phenolic compounds and essential oil components in larger amounts compared to organic and chemical fertilization.The latter enhanced the biomass yield, whereas organic fertilization in Serafim elicited a large leaf surface, which is desirable for food or decoration purposes.
In the present research, the genotype proved to be a factor showing a major influence, regardless of the fertilization type, which is essential to achieve specific targets such as a larger amount of a certain component of the essential oil (microorganism fertilization), a higher yield (chemical fertilization) or a larger leaf surface (organic fertilization).

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/agronomy12123219/s1,Table S1: The correlation matrix between the analyzed variables of basil.Funding: The authors wish to thank "Ion Ionescu de la Brad" Iasi University of Life Sciences for the financial support of the experiments.
Data Availability Statement: Not applicable.

Figure 2 .
Figure 2. The 2D principal subspace for different cultivars of basil and fertilization treatments.With red color are the experimental factors (Cultivar × Fertilization type), and with blue color are the parameters and bioactive compounds analyzed at sweet basil.

Table 1 .
Meteorological conditions during the study.

Table 2 . Influence of cultivar and fertilization type on the biometric characteristics. Treatment Plant Height (cm) Ramifications (No. per Plant) No. of Leaves per Plant Leaf Area Index (LAI) (cm 2 per Plant)
Within each column: n.s.-no statistically significant difference; values associated with the same lowercase letters are not statistically different at p < 0.05 according to Duncan's test.

Table 3 .
Influence of cultivar and fertilization type on basil yield characteristics.Within each column: n.s.-no statistically significant difference; values associated with the same lowercase letters are not statistically different at p < 0.05 according to Duncan's test.

Table 4 .
Interaction between cultivar and fertilization type on the biometric characteristics.

per Plant) No. of Leaves per Plant Leaf Area Index (LAI, cm 2 Per Plant)
Within each column, values associated with the same lowercase letters are not statistically different at p < 0.05 according to Duncan's test.AB-'Aromat de Buzau'; M-'Macedon'; C-'Cuisoare'; S-'Serafim'; Ch-chemical; O-organic; Mo-microorganisms.

Table 5 .
Interaction between cultivar and fertilization type on basil yield characteristics.

Table 6 .
Influence of cultivar and fertilization type on physiological and color parameters.
Within each column: n.s.-no statistically significant difference; values associated with the same lowercase letters are not statistically different at p < 0.05 according to Duncan's test.CCI-Chlorophyll Content Index; L-lightness-darkness; a-redness-greenness; b-yellowness-blueness.

Table 7 .
Interaction between cultivar and fertilization type on physiological and color parameters.Within each column, values associated with the same lowercase letters are not statistically different at p < 0.05 according to Duncan's test.