Development of Elite Mother Palms from the Best Performing Slow Vertical Growth Oil Palm (Elaeis guineensis Jacq.) Genotypes

1. Introduction

The Oil Palm (Elaeis guineensis Jacq.) is popularly known as golden palm one of the most economically important and largest vegetable oil yielding perennial crop (4-6 t oil/ha/ year), one acre of Oil palm plantations may yield up to 10 times more oil than other major oilseed crops. Indonesia stands as the dominant player, contributing a massive 59% of the global production, translating to 47 million metric tons. Following Indonesia, Malaysia holds a substantial 24% share with 19 million metric tons. India ranks 13 in Palm Oil Production, with 305,000 Metric Tons, 0.4% of global production [27]. Thailand, Colombia, and Nigeria also contribute notable shares, with 4%, 2%, and 2% respectively. Other oil palm growing countries make up the remaining percentage, each contributing between 0.74% to 1% of the global production. This distribution highlights the concentration of palm oil production in Southeast Asia, particularly in Indonesia and Malaysia, which together account for over 80% of the world's palm oil production. The data provides a clear picture of the key players in the palm oil market and their respective contributions. There is a significant disparity between the actual and potential use of oil palm germplasm [20]. The enhancement of oil palm crops through the extensive use of a limited number of closely related parents may lead to a restricted genetic base and inbreeding depression. To address this issue, oil palm breeders have identified genetically diverse, trait-specific germplasm lines through the study and characterization of existing oil palm germplasm. The primary goal for oil palm breeders is to create high-yielding cultivars with a broad genetic foundation [12]. The most crucial characteristic for oil palm production is oil yield, making its optimization the main objective in enhancing trait-specific oil palm varieties [14]. One of the biggest challenges in any hybridization program is finding the optimal pairing of two (or more) parental genotypes to maximize variance within related breeding populations. This maximizes the likelihood of identifying superior transgressive segregants in the segregating populations. The most difficult task for plant breeders in every hybridization attempt is figuring out how to blend multiple parental genotypes effectively to maximize diversity within related breeding populations. While DxD oil palm genotypes are not ideal for direct commercial oil production, they play a crucial role in breeding programs. By assessing and selecting for key traits such as mesocarp content, fruit size, disease resistance, and growth vigor, breeders can develop promising lines that contribute to the production of high-yield tenera palms. The integration of traditional breeding techniques with modern technologies like MAS will further enhance the efficiency and effectiveness of these programs, ultimately leading to more productive and sustainable oil palm cultivars [4,10]. Genetic Coefficient of Variation (GCV) indicates the proportion of variability in traits due to genetic differences among the genotypes. Phenotypic Coefficient of Variation (PCV) includes both genetic and environmental contributions to trait variability [26]. Higher GCV values suggest that there is significant genetic variation for that trait, which could be exploited in breeding programs. Selecting progenies based on one or more descriptors can cause unwanted modifications in others due to negative correlations between them [16]. To address this issue, principal component analysis (PCA) is utilised as a multivariate method to discover and identify the independent descriptors that influence plant traits [1]. To address the harvesting issue, we must develop dwarf or low height increment genotype. Reduced height increments have a long-term impact on harvesting costs [3,22]. The oil palm sector prioritizes height reduction due to the high expense of harvesting tall palms [25]. Identifying dwarf palms with commercial significance is challenging, with few known instances globally. Harvesting is a serious issue in elderly oil palm estates, hence the creation of dwarf palms is extremely beneficial to farmers.

2. Results

The comparison of means using the Tukey test showed that genotypes means performance for BI, BN and BW are 0.24 % , 7.49 and 105.95 Kg/palm/year, respectively. Whereas, for GT, HI, HT, LA, SR and TDM are 271.96, 32.99, 252.37, 7.22, 0.52 and 226.08% (Table 1 and Table 2). It is common to observe a negative correlation between bunch quantity and average bunch weight.

Descriptive statistics were generated from preliminary data analysis to describe the 308 Oil Palm dura genotypes for each attribute under this study. The range of variation and the potential for selection of the genotypes were determined by estimating the mean, maximum, and lowest values of these important traits [17]. In order to confirm the genotype variability of each trait relative to the mean, the Coefficient of Variation (CV) was also computed. The lower the CV, the more homogenous the data set for a particular trait. The x-axis represents Bunch Number and the y-axis represents Bunch Index (Figure 1).

Individual data points are marked with orange circles. The background color represents the density of data points. Blue areas indicate lower density, green areas indicate moderate density, yellow areas indicate higher density, and red areas indicate the highest density of data points. A blue line indicates a linear regression fit, showing the overall trend in the data. The positive slope of the regression line suggests a positive correlation between "Bunch Number" and "Bunch Index." As the bunch number increases, the bunch index also tends to increase. The highest density of data points (red and yellow areas) is concentrated around a certain range of bunch numbers (approximately 1 to 10) and bunch indices (approximately 0.1 to 0.4).There are fewer data points in the blue areas, indicating lower density. A few data points lie outside the main dense cluster, suggesting outliers with unusually high or low values for the given relationship. The scatter plot with the density heat map and regression line provides a clear visualization of the relationship between bunch number and bunch index in the dataset. The positive correlation implies that as the bunch number increases, the bunch index also tends to increase. This information can be useful for understanding growth patterns and selecting traits for breeding programs [18]. Individual data points are marked with orange circles. The background color represents the density of data points. Blue areas indicate lower density, green areas indicate moderate density, yellow areas indicate higher density, and red areas indicate the highest density of data points. A blue line specifies a linear regression fit, showing the overall trend in data. The positive slope of the regression line suggests a positive correlation between "Height Increment" and "Bunch Weight." As the height increment increases, the bunch weight tends to increase as well. The highest density of data points (red and yellow areas) is concentrated around a certain range of height increments (approximately 20 to 50) and bunch weights (approximately 50 to 150). There are fewer data points in the blue areas, indicating lower density. A few data points lie outside the main dense cluster, suggesting outliers with unusually high or low values for the given relationship. The genotypes are scattered based on their values in different traits (Figure 2).

High GCV (62.5 %) and PCV (45.83%) of BI indicate significant genetic variability, this trait appears to have substantial genetic influence. High GCV (66.9%) and PCV (52.47%) of BN suggest significant genetic variability and potential significance in genetic studies. Extremely high GCV (83.5%) and PCV (56.01%) of BW suggest that this trait has a high genetic influence and variability, indicating significance. Lower GCV (13.8%) and PCV (11.16%) of GT indicate less variability, suggesting this trait might be less influenced by genetic factors and more stable. Moderate GCV (40.9%) and PCV (33.59%) of HI indicate genetic influence and variability, which may be significant. Moderate to high GCV (43.4%) and PCV (33.00%) of HT suggest genetic influence and significance. Similar GCV (40 %) and PCV (34.07%) of LA values indicate genetic influence and potential significance. High GCV (88.5%) and PCV (57.69%) of SR values suggest a significant genetic contribution, indicating importance in genetic studies.

Simple correlations between oil palm yield and morphological parameterss, when analysed across 308 DxD genotypes, revealed both positive and negative correlations between above said parameters (Figure 3).

The highlights were the positive and extremely significant (P≤0.05) associations between the BI (0.87) and BW. Similarly, BW was emphatically profoundly critical (P≤0.05) with SR (0.42), HI (0.47), LA (0.15), BN (0.86), TDM (0.67); however, it was negatively significant with GT (-0.04) and HT (-0.03). Height increment (HI) had a significant correlations with SR (0.23), LA (0.25), BN (0.30), BW (0.47), TDM (0.72) BI (0.16) and HT (0.08). There were positive and significant relationships between the BN and other important parameters such as SR (0.45), HI (0.30), LA (0.02), TDM (0.47), and BI (0.84), but not with the GT (-0.02) and HT (-0.06). Furthermore, the BI revealing a optimistic association with SR (0.36) and is highly significant (P≤0.01) with BN (0.84) and BW (0.87).

The heatmap cluster analysis yields two types of dendrograms: genotype dendrograms with a horizontal orientation and character dendrograms with a vertical orientation (Figure 4). Based on the dendrogram of genotype, 308 items are divided into two major clusters, namely the first cluster consisting BN, HT and TDM and the other group consisting of remaining traits. The heatmap represents the clustering of different genotypes and traits (Figure 6). The horizontal dendrogram at the top clusters genotypes, while the vertical dendrogram on the left clusters traits. The color scale signifies the series of values, with blue representing minimum values and red representing maximum values. Each cell in the heatmap corresponds to a specific genotype-trait combination, with the color intensity reflecting the magnitude of the value. The numbers at the bottom (0.63, 0.22, 9.97, etc.) likely represent the mean or scaled values for each trait across the genotypes.

Genotype clustering (Horizontal Dendrogram): Genotypes are grouped based on similarity in their trait profiles. Close branches indicate genotypes with similar characteristics. Trait clustering (Vertical Dendrogram): Traits are grouped based on their correlation across genotypes. Close branches indicate traits that vary similarly across the genotypes. Color Intensity: Blue cells indicate lower trait values, while red cells indicate higher trait values. The intensity of the color provides a visual indication of the relative magnitude of each value. There is a clear distinction in trait values across different genotypes, with some genotypes showing consistently higher values (red) and others showing lower values (blue). Traits such as those on the right side of the heatmap tend to have higher values, indicating they might be positively correlated with each other.

In order to get more trustworthy information on how to recognize genotype groups with desired qualities for oil palm breeding, the PCA was carried out (Figure 5). Out of the nine principal components (PCs), only three had Eigen values greater than one, accounting for 67.19% of the variance across the various genotypes [2,23]. As seen in Figure 2, PC1 had a larger contribution to variability (36.56%), followed by PC 2 (15.87%). With the exception of height and girth, all attributes had positive factor loadings on the PC1.

3. Discussion

The scatter plot with the density heatmap and regression line provides a clear visualization of the relationship between height increment and bunch weight in the dataset. The positive correlation implies that as the height increment of the plants increases, the bunch weight also tends to increase [5]. Families with high bunch weights typically have low numbers, and vice versa; if bunch number is decreased by removing inflorescences, the weight of the remaining bunches increases. Bunch weight increases with palm age, whereas number declines.

Moderate GCV (37.9%) and PCV (45.07%) of TDM indicate genetic influence, making it potentially significant. The GCV indicates the proportion of variability in traits due to genetic differences among the genotypes. PCV includes both genetic and environmental contributions to trait variability. Higher GCV values suggest that there is significant genetic variation for that trait, which could be exploited in breeding programs. The closer the GCV and PCV values are, the lower the environmental influence on the expression of the trait, indicating that selection based on phenotype will be more effective. Traits like SR, BW, and BN show high GCV, indicating substantial genetic variability. Assign appropriate weights to each trait based on their importance in the breeding objectives. Not all traits may be equally important, and some traits may have more significant economic or agronomic value [6]. Traits like GT have lower GCV, suggesting less genetic variability or more environmental influence, similar results was found with [24]. The differences between GCV and PCV indicate the relative influence of environmental factors, with smaller differences suggesting lower environmental impact.

An important yield-determining characteristic for Oil Palm FFB is sex ratio. The progenies differences in sex ratio and other agronomic characteristics were examined. The current results concur with the previous studies by [9] and [15]. The Oil Palm FFB and oil output are the main factors considered in breeding and selection. Consequently, a reduction in male inflorescence raises the SR, which in turn causes a poor fruit set, low fresh fruit bunch yield, and low oil output. An increase in female inflorescences may put at risk the formation of male inflorescences. Based on genotype and male inflorescence output, oil palm sex ratios can be raised or lowered. The oil palm inflorescence sex, according to [8], generated changes throughout time in response to both internal and external factors, consequential in periodic cycles of male and female inflorescences.

Consider the correlation between traits, positive correlations between traits imply that selecting for one trait may inadvertently improve another, while negative correlations may require trade-offs [11]. Correlation analysis was used to identify the variables that had a strong link with one another [21]. One may get the conclusion that two or more variables can both explain the same phenomenon when they have a strong correlation with one another. The negative association indicates that these features have a detrimental influence on the BW [Figure 5]. It is therefore suggested that potential genotypes can be chosen using the traits that have been investigated [22].

To identify and depict pair wise structural differences between items, cluster analysis is utilized. Using fairly similar clustering techniques, the genotypes and significant characteristics are clustered. The data were divided into two directions (genotypes and characteristics) using two-way clustering. One advantage of two-way clustering is that it allows for the simultaneous grouping of specific objects and their variables [7,28]. The heatmap cluster analysis effectively highlights the relationships among genotypes and traits [13], enabling the identification of patterns and correlations that are crucial for breeding programs. The distinct clusters suggest that specific genotypes and traits are closely related, which can inform selective breeding strategies to optimize desirable traits such as yield and protein content.

The present results revealed that PC1contributed maximum variability due to yield related traits. The characters with positive values indicate highest contribution and its importance in divergence, whereas negative values indicate the least contribution to the total divergence [2,29]

4. Materials and Methods

4.1. Experimental Area Description

During the growing season, a field experiment was undertaken at ICAR-Indian Institute of Oil Palm Research, Pedavegi, Eluru district, Andhra Pradesh. The location is located at around 16°48′ N latitude and 81°12′ E longitude. The location receives an average annual rainfall of around 785.6 mm has an average temperature range of 28 °C to 36 °C, and is primarily red sandy loam soil (Figure 6).

Figure 6. Dura improvement Block with 308 genotypes.

Figure 6. Dura improvement Block with 308 genotypes.

4.2. Oil Palm Genotypes

In this investigation, 308 Oil Palm genotypes planted during 2013 were utilized for the development of dwarf dura breeding lines. This collection of Oil Palm genotypes resulted from eleven crosses made between promising selected African germplasm (Zambia and Tanzania) planted in 1998 (Figure 7). Eleven crosses, each with four replicates, make up nine palms in a randomized complete block design (Table 3). As a result, these genotypes were chosen as a benchmark for evaluating new genotypes. The experimental data were draw (average of three years) together from 2021 to 2023.

4.3. Biometric Observations

The following nine traits were assessed in present study: bunch index (BI, %), bunch number (BN), bunch weight (BW, kg/palm/year), Girth (GT, cm), height (HT, cm), height increment (HI, cm), leaf area (LA, sqm), sex ratio (SR) and total dry matter (TDM, kg). Every harvest, the number of fruit bunches produced is recorded for each individual palm to note the FFB (Fresh fruit bunch) yield. During the busiest time of year, harvesting often occurs twice a month; otherwise, it happens just once. It is given annually as the quantity of bunches per palm we considered usually for noting FFB yield. The total weight of all harvested bunches is indicated in kilograms per palm per year (kg/palm/year) and every harvest weight of harvested FFB is recorded in kg throughout the year. For measuring height of the palm, using a long, specifically marked wooden stick, measured the height of the adult palm from ground level to the 41st leaf base using a measuring tape. Every year, the height of the palm was measured in centimeters (cm). The girth of the adult palm measured using a measuring tape at a height of 50 cm above the ground. The palm's circumference was measured in centimeters. The performance of individual oil palm can be known by calculating by using standard formulas and elite palms can be identified for further research programme. Height increment (HI) in cm for adult palm=Height/ (Age of the palm-2), Sex ratio (SR)=Female/(Male + Female + Hermaphrodite), bunch dry weight= 0.5275xFFB (kg/palm/year); Bunch Index (BI)= Bunch Dry Weight /Total Dry Matter in kg; and Leaf Area (LA) of adult palm in Sqm =0.57 x NLL x LLL x LLW /100 /100. The detailed formulae are given in supplementary table. The most promising genomic selection indices were those that correlated with deviations from the shared regression of the relevant characteristics. These indices were especially effective for managing the observed strong negative phenotypic and genomic correlations [19]. The aim was to maximize protein yield, either through elevated grain yield or increased protein content, while maintaining the population average for the other respective trait.

4.4. Data Analysis

To assess the performance of evaluated genotypes, analysis of variance, averages were compared using the HSD test (honestly significant difference) or Tukey’s test. Recurrent reciprocal selection (R.R.S.) has allowed oil palm populations to thrive throughout several mating cycles, resulting in a higher frequency of favorable alleles linked with desirable characteristics by using graphic pad software and OPSTAT.

5. Conclusions

This study investigated and identified potential traits that can be utilized in future oil palm breeding programs. The use of mother palms in breeding programmes improved prediction accuracy for individual traits as well as the genomic selection indices. This improvement was particularly notable when combining phenotypic and genomic information in a genomics-assisted selection approach during preliminary yield trials. The availability of genetic diversity, as well as suitable parent selection for developing seeds that add value to the created cultivars and address crop limiting concerns, are critical variables in breeding project success. This study examined 308 oil palm populations to analyze oil production components, vegetative features, and output. Selected genotypes G8, G300, and G221 had the lowest yearly height increments, measuring 28.98, 29.19, and 30.87 cm, respectively. The outcome of the present study shows that they are slow height increment genotypes with a high FFB yield (> 25 T/Ha). The oil palm industry future will involve in the development of innovative progenies with slow growth, high fresh fruit yields, and bunch components using elite dura-type female parents derived from these appraised populations. Meanwhile, genotype evaluation should focus on enhancing the capacity to choose extraordinary parental source to produce the finest DxP hybrid tenera palms.