Allometric Expression of Sesame Cultivars in Response to Various Environments and Nutrition

: Sesame cultivars perform differently in various environments comprising speciﬁc locations and years. Micronutrient fertilizers increase crop yields and can enhance resistance to environmental hazards and contribute to potential yield. For assessing the effect of environment and micronutrients, ﬁeld experiments were carried out at three different locations (BARI, URF and NARC) of Pothwar region, Pakistan, in two succeeding years (2014, 2015). Four sesame cultivars and three micronutrients (Zn, B and Mn) with a control treatment were used in RCB design with a split-plot procedure and four replications. The results showed signiﬁcant differences in the main effects of all tested factors (cultivars, micronutrients, locations and years). Maximum mean values of plant height, capsules plant − 1 , biomass yield and seed yield were recorded for SG-30, while TH-6 was the lowest for these parameters. Maximum mean values of capsule length, capsule weight, seeds capsules − 1 and the weight of 1000 seeds were recorded for TH-6 but vice versa for SG-30. Maximum capsule length and seeds capsule − 1 were observed at NARC, while other growth and yield attributes were maximum at BARI and URF, including biomass yield and seed yield. The interactive effect of cultivar × location × year was highly signiﬁcant for all growth and yield traits except capsule length. revealed a differential response for all traits at three locations during both years. Application of micronutrients signiﬁcantly augmented all growth and yield features of sesame. Micronutrient fertilizers (i.e., ZnSO 4 at 10 kg ha − 1 , borax at 10 kg ha − 1 and MnSO 4 at 5 kg ha − 1 ) would increase yield and yield-attributing traits in low- or medium-precipitation areas where suitable cultivars have been designated.


Introduction
Pakistan is greatly deficient in edible oil production. During 2020-2021 (July-March), the country's demand for edible oil was 3.291 million tons, whereas 0.374 million tons (11.4%) were locally produced, while the remaining 2.917 million tons (88.6%) of edible oil was covered by imports [1]. Therefore, an increase in the local production of oilseeds is needed. The Pothwar plateau is a key part of the subtropical dry land zone, which extends across latitude 32 • 10 to 34 • 9 N and longitude 71 • 10 to 73 • 55 E, with 1.82 million hectares in Pakistan. The parent material of Pothwar plateau soils is varied in nature, and the parent material is loess, alluvium, colluviums and mixed by nature. The soils have variable parent material and chemical composition, and are prone to erosion due to undulated landforms.
The rainfall is erratic and differs greatly from 400 mm in the southwest to 1500 mm in the northeast part of the region. Pothwar is categorized as semi-arid to humid [2] and more than 70% of annual precipitation falls in the summer months [3]. The people of Pothwar have conventionally looked to rainfed agriculture for their livelihood. Traditional oilseeds such as sesame, which has been grown for centuries, need to be more focused on and explored to minimize the dependency on imports of edible oil.
Sesame (Sesamum indicum L.) seed comprises oil (50%), protein (up to 25%), crude fiber (3.2%), carbohydrate (16-18%) and ash (5.7%), and is rich source of vitamin E, phosphorus and calcium [4,5]. Sesame is known as the queen of oilseeds because of its high level of polyunsaturated fatty acids such as oleic acid (29.3 to 41.4%) and linoleic acid (40.7 to 49.3%) [6]. There is a huge gap between the potential production and actual yield of sesame in this country. Several reasons, such as meagre input management (little or no fertilization), macro-and micronutrient deficiencies in the soils, and biotic and abiotic stresses, cause low production of sesame.
Several research studies on sesame and other oilseed crops have specified that there is a diverse range of differences between genotypes and environments [7]. The yield of genotypes is uneven and differs over locations and environmental variability; specifically, location effects, seasonal fluxes and their interaction strongly impact the phenology of genotypes concerning optimum production [8]. Optimizing the yield and quality of the crop and more economic profits can be obtained through nutrient fertilization, which can also minimize the environmental risks. The idea of integrated plant nutrition is to maintain the soil fertility for high crop productivity by using combined and cheap sources such as fertilizers with the plant-required nutrients at the optimum level [9,10]. Farmers do not often use micronutrient fertilizers for crop production under rainfed conditions; therefore, micronutrient deficiencies may occur, and lower soil moisture reduces the nutrients' availability to plants [11]. Deficiencies of micronutrients in crops have increased owing to soil erosion, exhaustive farming, micronutrient losses through leaching and using marginal lands [12]. The main obstacle lies in the physiological problems associated with micronutrient deficiencies, such as hormonal imbalances that lead to reduced sesame yields [13]. Furthermore, boron, zinc and manganese help plants in chlorophyll formation and increase the photosynthetic activity [14]. Several studies have documented the increase in the growth and yield attributes of sesame due to micronutrients such as boron, zinc and manganese [15][16][17]. B, Zn, and Mn deficiencies were observed in the soils of the Pothwar plateau [18].
This study aimed to examine the response of sesame cultivars to the different environments of Pothwar Pakistan and to evaluate the effect of Zn, B and Mn on the growth and yield features of sesame.

Field Experimentation and Treatments
Three experiments were carried out: one each at Barani Agricultural Research Institute (BARI) Chakwal (latitude 32.9 • N; longitude 72.8 • E; elevation, 523 m; annual rainfall, 450-550 mm); University Research Farm (URF), Chakwal Road, Rawalpindi (latitude 33.1 • N; longitude 73.0 • E; elevation, 510 m; annual rainfall, 650-850 mm) and the National Agricultural Research Center, Islamabad (NARC) (latitude 33.7 • N; longitude 73.1 • E; elevation, 504 m; annual rainfall, 1000-1200 mm). The experiments were performed by using 4 sesame cultivars (SG-30, TS-3, TS-5 and TH-6) and 3 micronutrients, namely zinc, boron and manganese, with 1 control treatment. Experiments were established in a randomized complete block design (RCBD) with four replicates. All plots were randomized in each replicate, with cultivars allotted in the main plots and micronutrients in the subplots. Crops were planted in July 2014 and 2015 with single-row hand drill in a plot size of 2.7 m × 5 m, with six rows in each plot placed 45 cm apart. A 5 kg ha −1 seed rate was kept, with a seeding depth of 3-4 cm. Recommended doses of fertilizers (N at 90 kg ha −1 and P 2 O 5 at 60 kg ha −1 ) were applied as urea and di-ammonium phosphate 3 of 16 (DAP) during the last plowing at all locations. All cultivation operations remained the same at each site. ZnSO 4 (10 kg ha −1 ), borax (10 kg ha −1 ) and MnSO 4 (5 kg ha −1 ) were applied as the Zn, B and Mn sources, respectively. The side banding method was used for micronutrient application just after sowing across all locations. Manual thinning was carried out after the complete emergence to maintain the distance between plants at 10 cm.

Soil and Weather Data
Composite soil samples at two depths (0-15 cm and 15-30 cm) from all locations were collected before sowing and the samples were analyzed for physiochemical properties (Table 1). Zn and Mn were examined by DTPA extraction [19], while B was measured colorimetrically by means of azomethine-H [20] after hot water extraction [21]. Weather data of the experimental sites in both years were obtained from Pakistan's Meteorological Department, Islamabad (Figure 1).

Crop Data
Crop data were collected on the growth and yield-attributing characteristics of sesame. Data for plant height (PH) (using a measuring tape), number of primary

Crop Data
Crop data were collected on the growth and yield-attributing characteristics of sesame. Data for plant height (PH) (using a measuring tape), number of primary branches (NBP), secondary branches (NSB) and number of capsules (NC) were obtained from 10 randomly selected plants in each plot at the time of maturity, and the per-plant average was calculated. Data for capsule length (CL) and capsule weight (CW) were recorded from 20 capsules, which were randomly collected from plants in each plot. Firstly, the length of each capsule was measured by using a ruler scale, and the average capsule length was worked out. Secondly, the weight of a sample of 20 capsules was measured by using an electric balance, and the mean capsule weight was calculated. One hundred capsules of mature plants in each plot were collected randomly. After manual threshing of the capsules obtained, the total number of seeds was counted with the help of a seed counter, and the mean number of seeds per capsule (SPC) was calculated. Mature plants in the 2 central rows of a 2 m row length were manually harvested from each plot, and the tagged bundles of plants were kept vertical in the field. After few days of sun drying, the biomass yield (BioY) of each plant bundle was recorded using a weight balance, and the results were converted into kg ha −1 . Later, the plant bundles were threshed manually, and clean seeds were obtained by removing plant debris. The weight of each seed sample was recorded using an electric balance, and the seed yield (SY) in kg ha −1 was calculated. Three batches of 1000 seeds were isolated from the bulk seed sample obtained from each plot by using a seed counter, the weight of each batch was recorded utilizing an electric balance, and the average weight of 1000 seeds (SW1000) was calculated.

Statistical Analysis
A general linear model was adopted for analysis of variance (ANOVA) using a computer program (statistix 8.1). An RCB design (factorial) was used to determine the significance level of the main experimental factors (cultivars, micronutrients, locations and years) and their possible interactions. Furthermore, to identify significant variation between the means, the least significance test (LSD) was applied with a probability level of 5%. The linear association between the agronomic traits were determined by Pearson's correlation, and the correlation coefficient (r) values [22] along with the p-values were computed using statistix 8.1. Additionally, principal component analysis (PCA) [23] was executed using R software (Version 4.0.0) to classify numerous agronomic traits of sesame in response to the treatment interactions (cultivar × location × year).

Cultivars
The main effects of sesame cultivars for observed growth and yield characteristics were found to be highly significant at p < 0.01 (Table 2). Maximum values of PH, NC, biomass yield and seed yield (725.0 kg ha −1 ) were observed for the cultivar SG-30, but this was lowest for CL, CW, SPC and 1000 SW. Maximum NPB and NSB were observed for TS-3. On the other hand, TH-6 had the highest values for CL, CW, SPC and 1000 SW. However, it was lowest for the remaining attributes including seed yield (363.4 kg ha −1 ). The cultivar TS-5 was statistically on par with SG-30 for seed yield, while TS-3 was statistically parallel with SG-30 for seeds capsule −1 , and capsule length and weight ( Table 2).
** Significant at p ≤ 0.01. * Significant at p ≤ 0.05. NS = Non-significant. Values in columns having the different letter are significantly different at p ≤ 0.05.

Locations
Location means revealed statistically significant differences (p < 0.01) for all studied traits ( Table 2). BARI achieved the maximum 1000 SW and seed yield (660 kg ha −1 ), while URF produced sesame plants with the maximum PH, NPB, NSB, capsules plant −1 (47.4), CW and biomass yield (3161 kg ha −1 ). Although URF and BARI were statistically on par for the locational means of capsules plant −1 , CW, SPC and 1000 SW, NARC had the maximum results for CL and seeds capsule −1 , but this location produced sesame plants with minimum height, NPB, NSB, CW, 1000 SW, biomass yield (2071 kg ha −1 ) and seed yield (597.6 kg ha −1 ).

Year
Year means showed differential responses for the observed parameters, except CL, SPC and seed yield (Table 2), where both years revealed statistically parallel (p > 0.05) results for these traits. However, PH, NPB, NSB and biomass yield (2657 kg ha −1 ) were significantly higher during the first year (2014). However, capsules plant −1 , CW and 1000 SW were significantly (p < 0.01) superior during second year (2015).

Cultivar × Location × Year Interactions
The three-way interactions of year, location and cultivar for the yield traits of sesame were highly significant (p < 0.01) for the recorded growth and yield parameters (except capsule length) (

Cultivar × Location × Year Interactions
The three-way interactions of year, location and cultivar for the yield traits of sesame were highly significant (p < 0.01) for the recorded growth and yield parameters (except capsule length) ( Table 2). Graphical presentations of the interactions (showing significant differences) for quantitative characteristics of sesame are shown in Figures 2-10.         The data represented in Figure 2 indicated that plant height was maximum for TS-5 (140 cm), which was statistically parallel with TS-3 and SG-30 at URF in 2014, while minimum plant height (78.5 cm) was noted for TH-6 at NARC in 2014.
It is evident from the data that SG-30 produced the maximum (10.02) number of primary branches at BARI during 2014 (Figure 3), while TH-6 produced the minimum number of primary branches (0.14) at NARC during 2015. The maximum number of secondary branches per plant was recorded for TS-3 (8.81) at URF in 2014. However, TH-6 had similarly low results for secondary branches (0.03) at BARI and NARC during 2015 ( Figure 4). The results revealed that SG-30 produced the maximum primary branches across all locations during both years. Being a single-stem cultivar, TH-6 showed the minimum results for primary and secondary branches per plant across all locations during both years. However, the location URF in 2015 was found to be more conducive for producing the maximum secondary branches.
The maximum number of capsules plant −1 was recorded for TS-5 (74.2) at URF during 2015, which was statistically parallel with SG-30 (71.8) at BARI during the second year ( Figure 5), whereas TH-6 produced the fewest capsules plant −1 (14) at BARI during 2014. In addition, TH-6 produced the minimum capsules plant −1 across all locations during both years, while the rest of the cultivars revealed differential responses at three locations and two years.
The data for capsule weight (g) revealed that the maximum capsule weight was observed for TH-6 at BARI (0.427 g), which was statistically on a par with URF during the second year ( Figure 6). However, the minimum capsule weight (0. The data represented in Figure 2 indicated that plant height was maximum for TS-5 (140 cm), which was statistically parallel with TS-3 and SG-30 at URF in 2014, while minimum plant height (78.5 cm) was noted for TH-6 at NARC in 2014.
It is evident from the data that SG-30 produced the maximum (10.02) number of primary branches at BARI during 2014 (Figure 3), while TH-6 produced the minimum number of primary branches (0.14) at NARC during 2015. The maximum number of secondary branches per plant was recorded for TS-3 (8.81) at URF in 2014. However, TH-6 had similarly low results for secondary branches (0.03) at BARI and NARC during 2015 ( Figure 4). The results revealed that SG-30 produced the maximum primary branches across all locations during both years. Being a single-stem cultivar, TH-6 showed the minimum results for primary and secondary branches per plant across all locations during both years. However, the location URF in 2015 was found to be more conducive for producing the maximum secondary branches.
The maximum number of capsules plant −1 was recorded for TS-5 (74.2) at URF during 2015, which was statistically parallel with SG-30 (71.8) at BARI during the second year ( Figure 5), whereas TH-6 produced the fewest capsules plant −1 (14) at BARI during 2014. In addition, TH-6 produced the minimum capsules plant −1 across all locations during both years, while the rest of the cultivars revealed differential responses at three locations and two years.
The data for capsule weight (g) revealed that the maximum capsule weight was observed for TH-6 at BARI (0.427 g), which was statistically on a par with URF during the second year ( Figure 6). However, the minimum capsule weight (0.209 g) was noted for TS-3 at NARC during the second year. It is apparent from Figure 6 that the capsule weight of TH-6 was considerably higher across all locations during both years.
Similarly, the cultivar TH-6 achieved statistically parallel high numbers of seeds capsule −1 at NARC, i.e., 67.4 and 65.7 in 2014 and 2015, respectively. Likewise, seeds capsule −1 was higher for TS-5 at NARC during both years, but the lowest number of seeds capsule −1 (52.5) was recorded for SG-30 at NARC during the first year. Other interaction means were not significantly different (Figure 7).
In the same way, TH-6 produced the maximum 1000 SW at BARI (3.67 g), followed by URF during the second year. The minimum 1000 SW was observed for SG-30 and TS-5 (statistically similar) at NARC during 2015 (1.97 g and 2.00 g, respectively). The data revealed that BARI and URF were more favorable for the seed indices of sesame (Figure 8).
Three-way interaction means for biomass yield revealed that TS-3 produced the highest biomass yield (4264 kg ha −1 ) at URF in 2014, which was statistically on par with 2015. The lowest value (885 kg ha −1 ) of biomass yield was reckoned for TH-6 at BARI during the first year. However, TH-6 produced the minimum biomass yield across all three locations during both years, while the rest of the cultivars revealed differential responses for biomass production (Figure 9).
The results in Figure 10 indicated that the maximum seed yield (904.3 kg ha −1 ) was recorded for SG-30 at BARI in 2015, which was statistically on par with TS-5 at URF in the first and second year. However, the lowest seed yield (294.6 kg ha −1 ) was observed for TH-6 at BARI during 2014. Similar to biomass yield, TH-6 produced the lowest seed yield across all locations during both years.

Correlations among Agronomic Traits
Seed yield revealed a significant (p < 0.01) and positive correlation with PH, NPB, NSB, NC and biomass yield (Table 3). Similarly, biomass yield was significantly and positively correlated with these traits. However, biomass and seed yield were negatively correlated with CL, CW, SPC and seed indices. Usually, yield-attributing traits positively contribute towards final seed yield, but in the present case, the cultivar TH-6 revealed some unique characteristics, such as longer and heavier capsules, and higher SPC and seed indices, although, because of the single stem, the shortest plants and the lowest number of capsules, it had the lowest results for biomass and seed yield. Therefore, other growth and yield-attributing traits were negatively correlated with CL, CW, SPC and seed indices.

Principal Component Analysis
The total eigenvalues represent the quantity of characteristics under evaluation. The first two factors denote the information contained in 5.85% and 1.78% of the original quantity of attributed traits. From the interaction of cultivars, location and years, the first three important components were basic, controlling 85.9% of the total variability ( Table 3). The initial two components had the greatest influence with an aggregate contribution to the total variation of 76.3%, with an eigenvalue greater than 1. Agronomic traits including BioY, NPB and SY had the highest positive estimation in PC1, followed by 1000-SW, PH and NC in the second component's segment, whereas CL and CW had the highest negative stacking in PC1, followed by SPC in the second principal component. Among all PCs, PC1, PC2 and PC3 contributed 58.4%, 17.8% and 9.6%, respectively, to the total diversity (Table 4).
The associations among the various traits and sesame cultivars are additionally represented by the biplot in Figure 11. The biplot illustrates that the mean performance of the treatment interaction explained 76.3% of the total variation. The vast majority of the treatment interactions (cultivar × location × year) were dispersed in all coordinates of the principal components. The sesame parameter BioY was positively associated with SY, NBP and NSB, and negatively connected with CL, CW and SPC; most of these characteristics had higher magnitudes, represented by the length of vector, and made a correspondingly greater contribution to the overall variation. A few variables displayed smaller magnitudes, with small or medium vector lengths, such as SPC, PH and NC. The traits with close vector angles presented positive correlations with each other, as exhibited by BioY, SY, NBP and NSB. 3). The initial two components had the greatest influence with an aggregate contribution to the total variation of 76.3%, with an eigenvalue greater than 1. Agronomic traits including BioY, NPB and SY had the highest positive estimation in PC1, followed by 1000-SW, PH and NC in the second component's segment, whereas CL and CW had the highest negative stacking in PC1, followed by SPC in the second principal component. Among all PCs, PC1, PC2 and PC3 contributed 58.4%, 17.8% and 9.6%, respectively, to the total diversity (Table 4). The associations among the various traits and sesame cultivars are additionally represented by the biplot in Figure 11. The biplot illustrates that the mean performance of the treatment interaction explained 76.3% of the total variation. The vast majority of the treatment interactions (cultivar × location × year) were dispersed in all coordinates of the principal components. The sesame parameter BioY was positively associated with SY, NBP and NSB, and negatively connected with CL, CW and SPC; most of these characteristics had higher magnitudes, represented by the length of vector, and made a correspondingly greater contribution to the overall variation. A few variables displayed smaller magnitudes, with small or medium vector lengths, such as SPC, PH and NC. The traits with close vector angles presented positive correlations with each other, as exhibited by BioY, SY, NBP and NSB.    The distance between the treatment interactions (cultivar × location × year) and the biplot origin is an exclusive proportion. In this manner, treatment interactions such as TH6 × BARI × 2015, TH6 × URF × 2015, TH6 × BARI × 2014, TH6 × NARC × 2014, TH6 × NARC × 2015 and TS3 × URF × 2014 are those that were present far away from the origin of the biplot, and were considered most diverse treatment interactions. Such treatment interactions may or may not be greater, but they may be useful for some important agronomic attributes.

Micronutrients
Zinc, boron and manganese treatments considerably amplified (p < 0.01) the growth and yield-related traits of sesame (Table 2). However, the increase over the control was 12.2%, 17.4% and 8.4% in capsules plant −1 ; 18.7%, 10.3% and 5.4% biomass yield; and 14%, 16.8% and 7.9% in seed yield for Zn, B and Mn applications, respectively. Application of zinc resulted in the maximum increase in vegetative growth traits such as PH, NPB and NSB plus biomass yield. Likewise, the maximum increase in all yield attributes including seed yield was detected for boron treatment followed by zinc treatment (without significant differences, except for capsules plant −1 ).

Discussion
The sesame cultivars revealed differences in the expression of their morphological characteristics due to genetic variability and the prevailing environmental factors. Similar explanations were suggested by the previous findings [24,25]. Among the branched cultivars, SG-30 and TH-5 produced the highest seed yield because of their dominant genetic ability and yield stability in all growing conditions, mainly due to producing a higher number of capsules plant −1 . However, TH-6 is an unbranched and short stature cultivar which was found to produce the lowest number of capsules, biomass yield and seed yield, but other traits concerning the yield of this cultivar were dramatically higher. Although the yield potential of TH-6 was moderately higher than other cultivars, if grown in favorable environmental conditions, it was highly susceptible to excessive soil moisture and high relative humidity, and reasonably tolerant to pests and diseases [26].
The results are the pooled means of sesame cultivars across three locations with different environmental conditions, particularly the amount and pattern of rainfall received during the cropping period. In addition, the sesame genotypes are unstable and diverse in their trait expression with changes in the growing environment and seasons (years) [27]. Therefore, sesame responded differentially in two growing seasons in terms of vegetative growth traits. The results of the current research are in agreement with previous work [24]. In the present study, the first year received comparatively higher rainfall during the growing period across all locations (Figure 1). Therefore, on average, growth attributes including biomass yield were higher during first year, which might be due to optimum soil moisture availability, adequate nutrient uptake and appropriate assimilate translocation into plant parts, which might have increased the photosynthetic capacity and resulted in higher vegetative growth.
The distinct effects of sesame in different environmental conditions could be attributed to irregular rainfall, soil and the phenotypic expression of the cultivars. Optimum soil moisture with proper soil aeration guaranteed an adequate oxygen supply, which regulated root respiration and the absorbance of nutrients. Similarly, optimum soil moisture availability owing to moderate and evenly distributed rainfall during the cropping period at BARI and URF ensured the higher growth and yield of sesame. Besides, the additional supply of macronutrients sustained the nutrient flow and high photosynthetic capacity with maximum assimilate supply to the sink that contributed to development, which might have ensured maximum growth and yield. However, sesame is highly prone to drought in several physiological stages but is principally vulnerable to heavy rainfall and waterlogging [28], which could hinder crop growth and yield. It is usually grown in overly dry areas with precipitation of 300-600 mm [29]. Excessive soil moisture limits the soil aeration and oxygen supply, which restricts root respiration and nutrient uptake [30]. At NARC, during later growth stages, high soil moisture due to heavy rainfall in the months of September 2014 and October 2015 (Figure 1) might have reduced the growth and development of sesame and resulted in drastic yield reductions at this site. Such stress conditions might cause restricted plant metabolism, poor plant-water-nutrients relationships, limited photosynthetic capacity, and a slower crop growth rate and dry matter accumulation. In addition, heavy precipitation can cause nutrient leaching, leading to low productivity of the sesame crop. A considerable decline in the growth and yield of sesame due to excess soil moisture was also reported by [31].
Zinc and boron were found to be deficient in the soils of the experimental sites, while manganese was at a marginal level (Table 1). Therefore, the application of deficient micronutrients along with an adequate supply of macronutrients enhanced sesame growth and yield. Moreover, the increase in growth characteristics might be due to the critical role of zinc in hormone metabolism, protein and enzyme synthesis, and enzymatic activities [32], also increasing plants' tolerance to biotic and abiotic stresses, thus improving plant growth [33] and ultimately higher biomass. Additionally, boron is essentially required for the growth and development of a plant's reproductive parts [15,34]. Furthermore, boron improves the chances of fruit setting during pollen tube growth and germination [35]. Consequently, the highest results for NC, CL CW and 1000 SW were observed for boron application. Similarly, manganese is fundamentally crucial for plant metabolic activities and development, and it aids in chlorophyll formation and increases the accessibility of calcium and phosphorus [36]. Therefore, manganese application might have increased the growth and yield of sesame. The results of current study are parallel to the findings of previous research [15,17,37,38]. On the contrary, restricted growth and yield in the control treatments was due to a deficiency in micronutrients. For example, internode shortening and short height occurred as a result of reduced growth hormone development due to Zn deficiency [15]. Sesame is also very sensitive to boron deficiency. Its absence during the flowering stage decreases pollen viability and causes stamen and pistil abortion, resulting in poor seed set and low seed yield due to malformed branching [34]. Moreover, owing to boron deficiency, activation of enzymatic and non-enzymatic oxidation increased quinine concentrations and polyphenol oxidase, which are perilous for plant growth and development [39].

Conclusions
It can be concluded that TS-5 produced the maximum yield at URF, and SG-30 produced the maximum yield at BARI and NARC. These cultivars are recommended for general cultivation in the Pothwar region. Micronutrient fertilizer application produced significant increases in the productivity of sesame. Application of micronutrient fertilizers (10 kg ZnSO 4 ha −1 , 10 kg borax ha −1 and 5 kg ha −1 MnSO 4 ) along with the recommended dose of macronutrients significantly improved the productivity of sesame. The reduced production of sesame under high rainfall areas showed the low suitability of sesame under such conditions. Hence, low-and medium-rainfall areas with medium texture and well-drained soils are best suited for sesame cultivation.