Abstract
Nitrogen (N) and sulfur (S) were believed to be important nutrient management tools for the production of sweet basil (Ocimum basilicum L. ‘German’) with desirable oil content and composition and acceptable herbage yields. A multilocation research study was initiated to evaluate the effect of N (0, 60, 120, and 180 kg·ha−1 N) and S (0, 20, 40, and 80 kg·ha−1 S) rates on biomass production, oil content, and oil composition for sweet basil. The three locations in Mississippi (Stoneville, Poplarville, and Verona) were selected based on soil type, geographic, and climatic variation. Location, N rate, and their interaction were significant on basil dry herbage yields. The herbage yield means were 4967 kg·ha−1, 2907 kg·ha−1, and 2122 kg·ha−1 for Poplarville, Verona, and Stoneville, respectively. Oil content was significantly affected by location with means of 0.69%, 0.80%, and 0.64% for Stoneville, Poplarville, and Verona, respectively. Location, N, and S had significant effects on oil yields with means of 14.7, 38.7, and 18.5 kg·ha−1 for Stoneville, Poplarville, and Verona, respectively. The significant quadratic response of essential oil yields to N fertilization rates showed oil yields were maximized for fertilization between 50 and 60 kg·ha−1 N. In contrast, the response to S fertilization appeared to continue to increase beyond the maximum fertilization rate evaluated of 80 kg·ha−1 S. Location and N application rates had a significant effect on the yields of the major basil oil constituents (–)-linalool, eugenol, (–)-bornyl acetate, and eucalyptol, whereas S had a significant effect on eucalyptol yield only. Eucalyptol concentration was positively correlated to the concentration of (–)-bornyl acetate. This is the first study to quantify (in real concentration) the response of the major sweet basil oil constituents (–)-linalool, eugenol, (–)-bornyl acetate, and eucalyptol as a function of N and S application rates. Also, it is the first study to demonstrate a strong response of basil oil yield to S. The results from this study demonstrated that N and S applications can be used as management tools with respect to sweet basil production, oil content, and oil composition.
Sweet basil (Ocimum basilicum L.) is the most widely grown species of the genus Ocimum L. (Javanmardi et al., 2002; Simon et al., 1990, 1999). Recent research demonstrated that sweet basil could be a feasible essential oil crop for Mississippi and possibly the southeastern United States (Zheljazkov et al., 2008a, 2008b); however, basil plant nutrition has not been investigated in the southeastern United States. Prior research in other parts of the world demonstrated a significant effect of plant nutrition in general (Golcz et al., 2006; Sifola and Barbieri, 2006; Singh et al., 2004a; Topalov, 1962; Zheljazkov, 1998) and N fertilization (Arabaci and Bayram, 2004; Golcz et al., 2006; Sifola and Barbieri, 2006; Singh et al., 2004a, 2004b; Yassen et al., 2003) on sweet basil productivity and oil composition.
Sulfur is another macronutrient with important functions in plants. It is a part of plant proteins as a component of the amino acids cysteine and methionine and coenzymes (R-SH) (Marschner, 1999). One of the plant nutrient deficiencies identified in both Europe and North America was found to be from a lack of S in the plants. Traditionally, growers have not applied S because: 1) a large S deposition as sulfate (SO4) or mild solution of sulfuric acid (H2SO4) was available from the atmosphere as a result of the significant SO2 emissions; 2) S was found as “impurities” in mineral fertilizers; and 3) S-SO4 was a byproduct in the formulation of other nutrient fertilizers such as ordinary phosphate. Hence, soils have been receiving S with the rain or with the application of N, P, or K mineral fertilizers for many years. As a result of the use of high purity fertilizers in recent decades and changes in processing, continuous cropping with high-yielding cultivars and a significant reduction of anthropogenic SO2 emission (MacGrath et al., 1996; Riley et al., 2002), these S input modes have been diminished or eliminated, whereas the S removal from agricultural land continued. As a result, many soils have been depleted of available S (Hedge and Babu, 2007; Morris, 1994; Srinivasarao et al., 2008). In addition, S-SO4 can be leached from the surface soil, especially in regions with high precipitation (above 1200 mm/year) such as in Mississippi. According to Stats Canada and the U.S. Environmental Protection Agency (Statistics Canada, 2003; USEPA, 2006), SO2 emissions in North America (Canada and the United States) have been reduced by ≈50% for the period 1980 to 2000. These factors have resulted in a significant reduction of wet S deposition on agricultural land.
Sulfur deficiency could reduce N uptake, potential crop yield, and the quality of plants (Marschner, 1999). Previous research has found that SO4 additions increased yields in alfalfa (Medicago sativa), wheat (Triticum aestivum), corn (Zea mays), canola or oilseed rape (Brassica napus), and potato (Solanum tuberosum) (Millins and Mitchell, 1989; Pavlista, 2005; Rehm, 2005; Stewart and Porter, 1969; Weil and Mughogho, 2000). Agricultural crops vary widely with respect to their tolerance to low solution S and different crops may have different critical tissue S concentration (Hitsuda et al., 2005). Sweet basil requirement for S are unknown or have not been reported in the literature.
The hypothesis to be evaluated in this study was that N and S rates would have a significant effect on sweet basil productivity and oil content and would significantly impact the concentration and yield of (–)-linalool, eugenol, (–)-bornyl acetate, and eucalyptol. Furthermore, N and S fertilization of sweet basil could be used as a tool to target the production of basil with increased essential oil yields or oil with a desirable composition (an elevated concentration of a compound of interest). The chemical composition of basil is of a great interest to various industries that are using basil essential oil. Hence, in systems in which basil is grown for essential oil production, it is the oil composition and yield that are more important than herbage (shoot) biomass. The objective of this study was to evaluate the effects of N and S application rates on total basil production, oil content, dry herbage yield, essential oil composition, and yield of major oil constituents. The research has focused on the quantification of (–)-linalool, eugenol, (–)-bornyl acetate, and eucalyptol, because these are the major important basil oil constituents and previous research indicated that these constituents were relatively conservative (stable) traits. The minor constituents fluctuate depending on chemotype or environmental factors (Kruger et al., 2002; Topalov, 1962; Zheljazkov et al., 2008a). Herewith, the actual (real) concentrations of the four major constituents in basil oil as a function of N and S fertilization are reported.
Material and Methods
Plant material and experimental design.
The experiment was conducted in 2006 at three locations in Mississippi: the North Mississippi Research and Extension Center in Verona (lat. 34°43′22′′ N, long. −88°43′22′′ W), the Delta Research and Extension Center in Stoneville (lat. 33°25′18″ N, long. −90°54′28′′ W), and at the South Mississippi Branch Experiment Station in Poplarville (lat. 30°50′11″ N, long. −89°32′44″ W) using sweet basil (O. basilicum L. ‘German’). Basil seedlings were produced from certified basil seeds (donated by the Research Institute for Roses and Medicinal Plants, Kazanluk, Bulgaria). All the basil seedling transplants for the three locations were produced in a greenhouse at the Verona location during May through June. Seeds were placed in 48-cell plastic trays filled with Metromix 300 (The Scotts Co., Marysville, OH) growth medium. Transplants were grown in the greenhouse for 40 d with a day/night temperature regime of 22 to 25 °C and 18 °C, respectively. Seedlings were fertilized weekly with 1.8 g of 20N–20P2O5–20K2O dissolved in 300 mL of water and irrigated daily as needed. This particular fertilizer did not contain S. Basil seedlings (≈12 cm in height) were transplanted into the field at the three locations in 19 June 2006. To follow the common horticultural practices in Mississippi and the southeastern United States, basil plants were transplanted in previously prepared raised beds. The irrigation was provided by a drip tape irrigation system (1.29 L·min−1/30 m with emitter spacing of 0.3 m) placed 3 to 4 cm into the soil in the middle of the bed at two locations and as overhead sprinkler irrigation at the Stoneville location.
Before soil preparation, soil samples were collected from the basil experimental sites and analyzed for extractable nutrients following the Mississippi State University Soil Testing Laboratory analytical procedures, which included the Lancaster extraction method (Cox, 2001). The three locations had different soil types and concentration of available plant nutrients (Table 1).
Selected initial soil characteristics (0 to 15 cm), the concentration of extractable nutrients in the soil, and the average daily minimum and maximum temperatures of the three locations in Mississippi.


The experiment at each location was a randomized complete block design with four replicates. The experimental units were single-bed plots (2.1 m × 6 m) with 40 plants per plot. Basil was planted in an offset manner with two rows on each bed 30 cm on centers. Phosphorus and K fertilizers (free of S and at the same rate across the experimental site at every location) were added, if needed, based on the soil tests results before the land preparation and subsequently incorporated. Both N (as ammonium nitrate) and S (as sulfur bentonite, 90% S) were applied 1 to 3 d before transplanting. Single N and S fertilization events were applied to this relatively short-season crop (Zheljazkov, 1998; Zheljazkov et al., 2008a). Basil plants at the three locations were harvested at the same time, on 10 Aug. (55 d after transplanting), at full bloom (the optimal commercial stage for oil production), when the content and the composition of the essential oil would be optimal (Simon et al., 1990; Topalov, 1962; Zheljazkov, 1998). Plants were harvested by hand cutting the herbage at ≈12 cm above the soil surface. Herbage fresh weight was recorded at the time of harvest. Because of the large number of samples and the fact that it was impossible to extract the fresh samples from all plots at the same time, the plants from all plots were dried at 40 °C for 72 h and the dry weight recorded.
Essential oil extraction and analyses.
The essential oils from all basil samples were extracted with steam distillation using a modified Clevenger collector apparatus (Furnis et al., 1989) (Scientific Glass, CA) using representative subsamples of 150 g of dry herbage (stems, leaves, and flowers) from each plot with a distillation time of 120 min. Basil essential oil yield was calculated as the weight (g) of oil per weight (g) of dry basil tissue.
Gas chromatography–mass spectrophotometry analyses.
Chemical standards and basil oil from the field experiment were analyzed by gas chromatography–mass spectrophotometry (GC-MS) on a Varian CP-3800 GC (Palo Alto, CA) coupled to a Varian Saturn 2000 MS/MS. The GC was equipped with a DB-5 fused silica capillary column (30 m × 0.25 mm, with film thickness of 0.25 μm) operated using the following conditions: injector temperature, 240 °C, column temperature, 60 to 240 °C at 3 °C/min then held at 240 °C for 5 min; carrier gas, He; injection volume, 1 μL (splitless). The MS mass ranged from 40 to 650 m/z, filament delay of 3 min, target total ion chromatogram (TIC) of 20,000, a prescan ionization time of 100 μsec, an ion trap temperature of 150 °C, manifold temperature of 60 °C, and a transfer line temperature of 170 °C.
Quantitative analysis.
Commercial standard eugenol was purchased from Aldrich (St. Louis, MO), whereas (–)-linalool, eucalyptol, and (–)-bornyl acetate were purchased from Fluka (Buchs, Switzerland). With five concentration points, an external standard least squares regression for quantification was used. All four analytes were used to formulate separate calibration curves. All calculations were performed by generation of standard curves within Varian's Saturn GC/MS Workstation software package Version 6.40 (Table 2). The chromatograms of each of the oils from the field experiments (Fig. 1 shows a sample chromatogram) were compared with the standard injections. The target peaks (Fig. 1) were confirmed by both retention time and mass spectral data. Confirmed integrated peaks were then used to determine the percentage of each chemical constituent in the essential oil based on the equation for the line (Table 3). There were three commercially unavailable analytes that were constituents in sweet basil samples (Table 3). Identification of these components was obtained using retention times, Kovats Indices, and mass spectra. Kovats Indices were calculated using the equation: KI (x) = 100[(log RT (x) – log Pz)/(log RT (Pz+1) – log RT (Pz)], where: RT(Pz) ≤ RT(x) ≤ RT(Pz+1), and P4 …P25 are n-paraffins.
Standard curve data of the commercial standards run on Varian CP-3800 GC coupled to a Varian Saturn 2000 mass spectrometry/mass spectrometry.


Commercially unavailable analytes identified by retention time and Kovats Index.



Representative basil chromatogram (Verona).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Representative basil chromatogram (Verona).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Representative basil chromatogram (Verona).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Statistical analysis.
Covariate analysis (Proc Mixed; SAS Institute, 2003) was executed on the 11 response variables to determine which response variables were significantly (P = 0.05) affected by location, S fertilization, N fertilization, and by their interactions with class variables location and replication. When response variables were affected by fertilization, polynomial regression parameters for S fertilization and N fertilization rates were estimated using Proc GLM (SAS Institute, 2003).
Results
Dry herbage yields.
Location, N rate, and their interaction were significant on basil dry herbage yields (DHY) (Table 4). Mean yields were 4967, 2907, and 2122 kg·ha−1 for Poplarville, Verona, and Stoneville, respectively, and were all significantly different from each other. Herbage yields at Poplarville were approximately twice the yields at the other two locations, indicating that location (environment) could be a major modifier of basil herbage yields. Nitrogen application rates greatly affected herbage yields (Fig. 2). In general, N application rates of up to 60 kg·ha−1 increased yields at each location. Further increases of N rates did not bring a corresponding significant yield increase. Sulfur application rates had no effect on herbage yields.
Analysis of variance P values for the main and interaction effects of location, N rate, and S rate on basil dry herbage yields (DHY), oil content, oil yields, (−)-linalool, eugenol, (−)-bornyl acetate (Bornyl A), and eucalyptolz.



Herbage production as affected by nitrogen fertilization at three locations (fitted curves have R 2 = 0.72).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Herbage production as affected by nitrogen fertilization at three locations (fitted curves have R 2 = 0.72).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Herbage production as affected by nitrogen fertilization at three locations (fitted curves have R 2 = 0.72).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Oil content and yields.
Oil content (as percent in dry herbage) was affected by location (Table 4) with means of 0.69%, 0.80%, and 0.64% for Stoneville, Poplarville, and Verona, respectively. These three means were significantly different from each other. Analysis of variance found an S fertilization by location interaction effect on basil oil content (Table 4). Regression on S was statistically significant but biologically insignificant with R 2 = 0.09. Thus, the oil content was affected by unknown location effects but was little affected by the fertilization. Overall, basil essential oil content at all locations varied from 0.39% to 0.97%, a typical range for oil content in sweet basil (Topalov, 1962; Zheljazkov et al., 2008a).
Basil essential oil yield was a function of herbage yields and essential oil content and was calculated as the weight (g) of oil per weight (g) of dried herbage for each treatment. Oil yields ranged from 5.5 to 56 kg·ha−1. Location, N fertilization, and S fertilization had significant effects on basil essential oil yields (Table 4). For the reason that the oil content was little affected by N fertilization but herbage yields (DHY) were, the oil yields showed similar (to herbage yield) responses to N fertilization. The average oil yields (means) for Stoneville, Poplarville, and Verona were 14.7, 38.7, and 18.5 kg·ha−1, respectively, and were different among locations. The high oil yields at Poplarville resulted from the overall greater average herbage yields and greater oil content relative to other locations. Sulfur fertilization had a significant effect on oil yields (R 2 = 0.62) with yields increasing after increased S application rates (Fig. 3). Overall, essential oil yields were maximized at an intermediate level of N application rates (at ≈50 to 60 kg·ha−1) (Fig. 4), whereas maximum response to S fertilization appeared to be beyond 80 kg·ha−1 S (Fig. 3).

Oil production as a function of sulfur fertilization at three locations (fitted curves have R 2 = 0.62).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Oil production as a function of sulfur fertilization at three locations (fitted curves have R 2 = 0.62).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Oil production as a function of sulfur fertilization at three locations (fitted curves have R 2 = 0.62).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Oil production as a function of nitrogen fertilization and location (fitted curves have R 2 = 0.54).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Oil production as a function of nitrogen fertilization and location (fitted curves have R 2 = 0.54).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Oil production as a function of nitrogen fertilization and location (fitted curves have R 2 = 0.54).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Concentrations of (–)-linalool, eugenol, (–)-bornyl acetate, and eucalyptol.
The concentration of (−)-linalool in basil oil ranged from 15% to 50% of the total oil. The means of (−)-linalool concentrations at Poplarville and Stoneville (31.0% and 31.1%), respectively, were almost equal, whereas the mean of (−)-linalool concentration at Verona was significantly greater at 44%. A regression of (−)-linalool concentration on S application rate was not significant. Regression on N application rate was significant but caused little change in (−)-linalool concentrations (Fig. 5), indicating N fertilization was of little biological significance to (−)-linalool concentration of basil essential oil.

Concentrations of (−)-linalool as a function of nitrogen fertilization and location (fitted curves have R 2 = 0.52 and lines for Stoneville and Poplarville appear as one line).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Concentrations of (−)-linalool as a function of nitrogen fertilization and location (fitted curves have R 2 = 0.52 and lines for Stoneville and Poplarville appear as one line).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Concentrations of (−)-linalool as a function of nitrogen fertilization and location (fitted curves have R 2 = 0.52 and lines for Stoneville and Poplarville appear as one line).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
In general, eugenol concentration in basil oil ranged from 2.2% to 5.4%, whereas (−)-bornyl acetate concentration ranged from 0.3% to 1.8% of the oil. Only location had a significant effect on eugenol or (−)-bornyl acetate concentrations, suggesting that the concentrations of these constituents might be relatively stable traits and could be influenced by the environment and climate but not by N or S fertilization. All location averages were statistically different (P = 0.05). The means of eugenol concentrations at the three locations were 2.5%, 4.1%, and 3.2%, whereas the means of (−)-bornyl acetate concentrations were 0.7%, 0.5%, and 1.3% for Stoneville, Poplarville, and Verona, respectively. The concentration of eucalyptol in basil oil varied from 2.8% to 9.2% and was significantly affected by location and by N fertilization. Locations were significantly different from each other with mean eucalyptol concentrations of 5.7%, 4.3%, and 7.7% for Stoneville, Poplarville, and Verona, respectively. Regression on N rate was statistically significant but caused a relatively small change in eucalyptol concentration (Fig. 6).

Concentrations of eucalyptol as affected by location and by nitrogen fertilization (fitted curves have an R 2 = 0.57).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Concentrations of eucalyptol as affected by location and by nitrogen fertilization (fitted curves have an R 2 = 0.57).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Concentrations of eucalyptol as affected by location and by nitrogen fertilization (fitted curves have an R 2 = 0.57).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Effect on yields per area basis (kg·ha−1) of the four major basil oil constituents.
Location and N application rates had significant effect on the yields (per area basis) of (−)-linalool, eugenol, (−)-bornyl acetate, and eucalyptol (Table 5). Sulfur application rate had an effect on eucalyptol yield, but not on the yield of the other constituents. The relation of eucalyptol yield to S fertilization was quadratic.
Tests for significance of linear and nonlinear effects of N fertilization rate and S fertilization rate at three locations and their interactions on yield of the four major constituents in basil essential oilz.


Yields of (−)-linalool were significantly affected by location and N fertilization (Table 5) with intermediate N rates (50 to 60 kg·ha−1) maximizing (−)-linalool yields per area basis (Fig. 7). Overall, (–)-linalool yield means were 4.91, 11.57, and 8.23 kg·ha−1 for Stoneville, Poplarville, and Verona, respectively. Eugenol yields were much greater at Poplarville than in the other two locations (Fig. 8) with yield means of 0.37, 1.54, and 0.61 kg·ha−1 for Stoneville, Poplarville, and Verona, respectively. Intermediate N fertilization rates (50 to 60 kg·ha−1) maximized eugenol yields per area basis (Fig. 8). The yields of (−)-bornyl acetate at the three locations were 0.10, 0.18, and 0.23 kg·ha−1 for Stoneville, Poplarville, and Verona, respectively. Like with the other constituents, intermediate N application rates near 50 to 60 kg·ha−1 resulted in maximum yields of (−)-bornyl acetate per area basis (Fig. 9). Eucalyptol yields were 0.83, 1.60, and 1.46 kg·ha−1 for Stoneville, Poplarville, and Verona, respectively. Like with the other two constituents, N application rates at ≈60 kg·ha−1 resulted in maximum yields of eucalyptol per area basis (Fig. 10). Although the effect of S rates on eucalyptol was significant, this resulted in relatively small yield change (Fig. 11). Interestingly, the maximum response to S seemed to be beyond the highest S rate of 80 kg·ha−1 (Fig. 11).

Yields of (−)-linalool as a function of location and nitrogen fertilization (fitted curves have R 2 = 0.65).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yields of (−)-linalool as a function of location and nitrogen fertilization (fitted curves have R 2 = 0.65).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Yields of (−)-linalool as a function of location and nitrogen fertilization (fitted curves have R 2 = 0.65).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of eugenol as a function of location and nitrogen fertilization (fitted curves have an R 2 = 0.43).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of eugenol as a function of location and nitrogen fertilization (fitted curves have an R 2 = 0.43).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Yield of eugenol as a function of location and nitrogen fertilization (fitted curves have an R 2 = 0.43).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of (−)-Bornyl acetate as a function of location and nitrogen fertilization. (Fitted curves have R2 = 0.43).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of (−)-Bornyl acetate as a function of location and nitrogen fertilization. (Fitted curves have R2 = 0.43).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Yield of (−)-Bornyl acetate as a function of location and nitrogen fertilization. (Fitted curves have R2 = 0.43).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of eucalyptol as a function of location and nitrogen fertilization (fitted curves have and R 2 = 0.39).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of eucalyptol as a function of location and nitrogen fertilization (fitted curves have and R 2 = 0.39).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Yield of eucalyptol as a function of location and nitrogen fertilization (fitted curves have and R 2 = 0.39).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of eucalyptol as a function of sulfur fertilization and location (fitted curves have an R 2 = 0.30).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415

Yield of eucalyptol as a function of sulfur fertilization and location (fitted curves have an R 2 = 0.30).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Yield of eucalyptol as a function of sulfur fertilization and location (fitted curves have an R 2 = 0.30).
Citation: HortScience horts 43, 5; 10.21273/HORTSCI.43.5.1415
Correlations between the measured responses.
The yields of the four oil constituents [(−)-linalool, eugenol, (−)-bornyl acetate, and eucalyptol] are a function of the percentage oil, oil yield, and herbage yield to the extent that they were calculated from those values. Therefore, those correlations (marked in bold) are expected to be significant (Table 6). There were some significant correlations among the measured responses (Table 6). The concentration of (−)-bornyl acetate in basil oil was negatively correlated to basil herbage yields, which might be an indication that greater amounts of this compound are synthesized and accumulated when the plant growth conditions are less favorable (resulting in lower herbage yields). The concentration of eucalyptol was correlated to (−)-linalool concentration and to oil content, i.e., the greater the oil content or (−)-linalool content, the greater the concentration of eucalyptol. Also, eucalyptol concentration was positively correlated to oil yield, presumably because it was positively correlated to oil content to which oil yields were a function. As expected, oil yields were positively correlated to oil content. The yields (per area basis) of the four major basil oil constituents [(−)-linalool, eugenol, (–)-bornyl acetate, and eucalyptol] were positively correlated to dry herbage yields, to the oil content and oil yields, which should be self-explanatory. (−)-Linalool yields were positively correlated to eucalyptol concentration, hence confirming the positive correlation (indicated previously) between the concentrations of these two constituents. Also, the yields of eugenol, (−)-bornyl acetate, and eucalyptol were positively correlated to (−)-linalool yields. Furthermore, yields of (−)-bornyl acetate and eucalyptol were positively correlated to the eugenol yields, and eucalyptol yields were correlated to the (−)-bornyl acetate yields.
Correlations and their P value statistical significances for all measurements.


Discussion
Overall, basil herbage yields in this study were relatively high and were similar to sweet basil yields reported in the literature or reported recently for Mississippi (Topalov, 1962; Zheljazkov et al., 2008a, 2008b). Basil essential oil content and yield in this study were within the typical range for basil reported in the literature for other countries (Anwar et al., 2005; Bowes and Zheljazkov, 2004; Juliani and Simon, 2002) and comparable to the results from recent research in Mississippi (Zheljazkov et al., 2008a, 2008b). The concentration of major oil constituents (−)-linalool, eugenol, (−)-bornyl acetate, and eucalyptol were within the range of these constituents in European sweet basil chemotypes (Marotti et al., 1996; Simon et al., 1990; Topalov, 1962; Zheljazkov et al., 2008a).
Although the concentration of the four major oil constituents is important with relation to the oil quality, the actual yield of each constituent may significantly affect the overall economic importance of the basil crops. It may be possible to use N and S fertilization as a tool to target the production of basil essential oil with a specific oil profile, i.e., constituent-specific oil. Some industries require (−)-linalool; others make use of eugenol or eucalyptol or (−)-bornyl acetate. For example, (−)-linalool (produced by fractional distillation of basil oil) is widely used as an aromatic agent in soaps, shampoos, detergents, as a biochemical pesticide, and even for the production of vitamin E. Furthermore, (−)-linalool can be used in the production of linalyl esters with different aromas. Hence, different industries may require special profiles of sweet basil oil, which could render premium prices for primary producers.
This is the first report to find such a strong response of basil essential oil yields to S fertilization, confirming part of our hypothesis for a possible significant effect of S fertilization on basil essential oil yields. In addition, to our knowledge, this is the first study to quantify (in real concentration) the response of the major sweet basil oil constituents (−)-linalool, eugenol, (−)-bornyl acetate, and eucalyptol as a function of N and S application rates in a variety of environments. In most cases, N rate, S rate, or N by S interaction resulted in altered oil composition of sweet basil but not in oil content. These results suggest that sweet basil reacts differently to the environmental conditions in different parts of Mississippi with respect to productivity and oil composition in general. In general, this and previous studies (Zheljazkov et al., 2008a, 2008b) demonstrated sweet basil could be grown as a viable essential oil crop in Mississippi and possibly other areas or regions in the southeastern United States.
This study partially confirmed the hypothesis that N and S rates, and probably their interaction, would have a significant effect on sweet basil productivity, essential oil content, and composition. However, the assumption that S addition may improve N use efficiency, which may have a direct effect on herbage yields, was not supported and this part of the hypothesis should be rejected. The lack of a strong yield response to N may be the result of: 1) relatively late transplanting time and the relatively short vegetation period of sweet basil in Mississippi; 2) relatively high fertility soils; 3) frequent heavy rainfalls and high temperatures; and 4) relatively shallow roots system of this plant. The relatively short period from transplanting to harvest of basil may not be sufficient for the complete N or S utilization from the applied fertilizer, especially for S, which is less soluble compared with most N sources. Previous research has demonstrated that sweet basil may require relatively low N application rates, as low as even 5 kg·ha−1 N, to provide good yields (Arabaci and Bayram, 2004). However, other studies found increased yield response with increasing N application rates of up to 80 kg·ha−1 N (Singh et al., 2004a), 80 to 120 kg·ha−1 N (Yassen et al., 2003), 200 kg·ha−1 N, or even up to 300 kg·ha−1 N (Sifola and Barbieri, 2006). These differences in sweet basil response to N application rates could be the result of differences in soil organic matter, soil organic N, and its mineralization rates among various soils. Hence, this and previous studies suggested that sweet basil yield response to N application rates depends largely on environmental and management conditions. Basil is a frost-sensitive species and is transplanted relatively late, when larger pools of inorganic N would be readily available as a result of mineralization of organic residues. On the other hand, with the extreme rain events that occur with some regularity in Mississippi, water may transport necessary nutrients beyond reach of the basil roots and may contribute to denitrification nullifying a potential yield response to N application rates.
Conclusions
From this work at three locations it has been demonstrated that location, N, and S had significant effects on sweet basil oil yields. Basil essential oil yields were maximized at N fertilization of ≈50 to 60 kg·ha−1 N, whereas maximum response to S fertilization appeared to be beyond 80 kg·ha−1 S. Location and N had significant effects on the yields of the major basil oil constituents (−)-linalool, eugenol, (−)-bornyl acetate, and eucalyptol, whereas S had a significant effect on eucalyptol yield only. We found that eucalyptol concentration was positively correlated to the concentration of (−)-bornyl acetate. This is the first study to quantify the response of the major sweet basil oil constituents (−)-linalool, eugenol, (−)-bornyl acetate, and eucalyptol as a function of N and S application rates. In addition, it is the first study to demonstrate a strong response of basil oil yield to S. Further research with higher S application rates and intermediate N background may be needed to estimate the optimal sweet basil yield and compositional response to N and S. In addition, further research is needed to explain the strong effect of location on basil productivity.
Literature Cited
Anwar, M. , Patra, D.D. , Chand, S. , Kumar, A. , Naqvi, A.A. & Khanuja, S.P.S. 2005 Effect of organic manures and inorganic fertilizer on growth, herb and oil yield, nutrient accumulation, and oil quality of French basil Commun. Soil Sci. Plant Anal. 36 1737 1746
Arabaci, O. & Bayram, E. 2004 The effect of nitrogen fertilization and different plant densities on some agronomic and technologic characteristic of Ocimum basilicum L. (basil) J. Agron. 3 255 262
Bowes, K. & Zheljazkov, V. 2004 Ocimum sanctum L. and Ocimum basilicum L. grown in Nova Scotia, Canada show potential as essential oil crops J. Amer. Soc. Hort. Sci. 129 789 794
Cox, M.S. 2001 The Lancaster soil test method as an alternative to the Mehlich 3 soil test method Soil Sci. 166 484 489
Furnis, B.S. , Hannaford, A.J. , Smith, P.W.G. & Tatchell, A.R. 1989 Vogel's textbook of practical chemistry 5th Ed Longman Scientific & Technical New York, NY 171 175
Golcz, A. , Politycka, B. & Seidler-Lozykowska, B. 2006 The effect of nitrogen fertilization and stage of plant development on the mass and quality of sweet basil leaves (Ocimum basilicum L.) Herba Polonica. 52 22 30
Hedge, D.M. & Babu, S.N.S. 2007 Correcting sulfur deficiencies in soils and crops Ind. J. Fertilizers. 3 65 79
Hitsuda, K. , Yamada, M. & Klepker, D. 2005 Sulfur requirement of eight crops at early stages of growth Agron. J. 97 155 159
Javanmardi, J. , Khalighi, A. , Kashi, A. , Bais, H.P. & Vivanco, J.M. 2002 Chemical characterization of basil (Ocimum basilicum L.) found in ‘Local’ accessions and used in traditional medicines in Iran J. Agr. Food Chem. 50 5878 5883
Juliani, H.R. & Simon, J.E. 2002 Antioxidant activity of basil 575 579 Janick J. & Whipkey A. Trends in new crops and new uses ASHS Press Alexandria, VA
Kovats, E. 1965 Gas chromatographic characterization of organic substances in the retention index system Giddings J.C. & Keller R.A. Advances in chromatography Vol. I Marcel Dekker, Inc New York, NY 229 247
Kruger, H. , Wetzel, S.B. & Zeiger, B. 2002 The chemical variability of Ocimum species J. Herbs Spices Med. Plants 9 335 344
MacGrath, S.P. , Zhao, F.J. & Withers, P.J.A. 1996 Development of sulphur deficiency in crops and its treatment Proc. Fertilizer Soc. No. 379 The Fertilizer Society London, UK
Marotti, M. , Piccaglia, R. & Giovanelli, E. 1996 Differences in essential oil composition of basil (Ocimum basilicum L.) Italian cultivars related to morphological characteristics J. Agr. Food Chem. 44 3926 3929
Marschner, H. 1999 Mineral nutrition of higher plants Academic Press London, UK
Millins, G.L. & Mitchell C.C. Jr 1989 Wheat varietal response to sulfur J. Plant Nutr. 12 909 922
Morris, R.J. 1994 Sulfur: Deficiencies spawn demand Farm Chem. Int. 8 20 21
Pavlista, A.D. 2005 Early-season applications of sulfur fertilizers increase potato yield and reduce tuber defects Agron. J. 97 599 603
Rehm, G.W. 2005 Sulfur management for corn growth with conservation tillage Soil Sci. Soc. Amer. J. 69 709 717
Riley, N.G. , Zhao, F.J. & McGrath, S.P. 2002 Leaching losses of sulphur from different forms of sulfphur fertilizers: A field lysimeter study Soil Use Manage. 18 120 126
SAS Institute 2003 SAS for Windows, v 9.1.3 SAS Inst Cary, NC
Sifola, M.I. & Barbieri, G. 2006 Growth, yield and essential oil content of three cultivars of basil grown under different levels of nitrogen in the field Sci. Hort. 108 408 413
Simon, J.E. , Morales, M.R. , Phippen, W.B. , Vieira, R.F. & Hao, Z. 1999 Basil: A source of aroma compounds and popular culinary and ornamental herb 499 505 Janick J. & Simon J.E. Perspectives on new crops and new uses ASHS Press Alexandria, VA
Simon, J.E. , Quinn, J. & Murray, R.G. 1990 Basil: A source of essential oils 484 489 Janick J. & Simon J.E. Advances in new crops Timber Press Portland, OR
Singh, K. , Singh, P.P. , Beg, S.U. , Kumar, D. & Patra, D.D. 2004a Effect of harvesting time and nitrogen application on essential oil yield and quality of a methyl chavicol rich cultivar of Indian basil (Ocimum basilicum L.) J. Medic. Arom. Plant Sci. 26 293 296
Singh, K. , Singh, P.P. , Beg, S.U. , Kumar, D. & Patra, D.D. 2004b Effect of NPK fertilizers on growth, oil yield and quality of French basil (Ocimum basilicum L.) J. Spices Arom. Crops. 13 52 54
Srinivasarao, C. , Wani, S.P. , Sahrawat, K.L. , Rego, T.J. & Pardhasarandhi, G. 2008 Zinc, boron, and sulfur deficiencies are holding back the potential of rainfed crops in semi-arid India: Experiences from participatory watershed management Int. J. Plant Prod. 2 89 99
Statistics Canada 2003 Agriculture Census 2003: Canadian agriculture at a glance Feb. 2006 <http://www.statcan.ca/english/Pgdb/agricul.htm>.
Stewart, B.A. & Porter, L.K. 1969 Nitrogen – sulfur relationships in wheat (Triticum aestivum L.) and corn (Zea mays) and beans, (Phaseolus vulgaris) Agron. J. 61 267 271
Topalov, V.D. 1962 Essential oil and medicinal plants Hr. G. Danov Press Plovdiv, Bulgaria
USEPA 2006 The acid rain challenge: Achieving progress through committed actions. SO2 emission reduction requirements Sept. 2006 <http://www.epa.gov/aimarkets/usca/brochure/acid_chart1htm>.
Weil, R.R. & Mughogho, S.K. 2000 Sulfur nutrition of maize in four regions of Malawi Agron. J. 92 649 656
Yassen, M. , Ram, P. , Yadav, A. & Singh, K. 2003 Response of Indian basil (Ocimum basilicum) to irrigation and nitrogen schedule in Central Uttar Pradesh Annals Plant Phys. 17 177 181
Zheljazkov, V.D. 1998 Basil 317 320 Phehlivanor M. , Moskov G. , Jankov B. , Terziev J. , Zheljazkov V. & Yantcheva H. Plant production edition 1 Academic edition of Higher Institute of Agriculture Plovdiv, Bulgaria
Zheljazkov, V.D. , Callahan, A. & Cantrell, C.L. 2008a Yield and oil composition of thirty-eight basil (Ocimum basilicum L.) accessions grown in Mississippi J. Agr. Food Chem. 56 241 245
Zheljazkov, V.D. , Cantrell, C.L. , Tekwani, B. & Khan, S. 2008b Content, composition, and bioactivity of the essential oil of three basil genotypes as a function of harvesting J. Agr. Food Chem. 56 380 385