Have a personal or library account? Click to login
Multi-Trait Composite Performance Index for the Identification of Elite Pistachio (Pistacia vera L.) Genotypes with Superior Kernel Efficiency Cover

Multi-Trait Composite Performance Index for the Identification of Elite Pistachio (Pistacia vera L.) Genotypes with Superior Kernel Efficiency

Open Access
|Dec 2025

Full Article

Introduction

Pistachio (Pistacia vera L.), a member of the Anacardiaceae family, represents one of the most economically significant nut crops globally, with cultivation spanning semi-arid regions across multiple continents. The global production of pistachios has increased dramatically over the past few decades, from around 50 thousand tons in 1970 to more than 1 million tons in 2020 (FAOSTAT, 2023). Based on FAOSTAT data, the most significant increases in area of cultivation occurred in Jordan, Madagascar, Kyrgyzstan, California, Turkey, Uzbekistan, Azerbaijan, and Mexico (FAOSTAT, 2023). However, Iran, Turkey, and the USA continue to maintain the highest pistachio production globally. The edible seed is fixed within a drupe, protected by a thin endocarp, and has attracted human attention for millennia (Eshankulov, 2017). The pistachio is sometimes called the golden tree in Turkey, Iran, and Syria, referencing its economic value (Eshankulov, 2017). Pistacia vera is a diploid (2n = 30), xerophytic, subtropical woody perennial; it is wind-pollinated, avian-dispersed, deciduous, and dioecious (Hormaza and Wünsch, 2007; Khanazarov et al., 2009). The species lineage dates back to between 1.82 and 3.38 million years ago (Xie et al., 2014). It is both heat tolerant and cold resistant (Khanazarov et al., 2009), growing in semi-arid deserts and dry mountain slopes (Vavilov, 1967). Its deeply penetrating roots allow it to weather arid periods (Whitehouse, 1957), and jasmonic acid production provides salt tolerance (Zeng et al., 2019). Studies conclude that the crop originated in Central Asia, where the only extant wild populations occur (Whitehouse, 1957; Zohary, 1952; Khanazarov et al., 2009). The current wild range covers foothill zones in Turkmenistan, Uzbekistan, Tajikistan, Kyrgyzstan, southeastern Kazakhstan, northern Afghanistan, and northeastern Iran (Zohary, 1952; Khanazarov et al., 2009). Wild pistachios grow along a narrow ecoline between 500 and 800 masl in Central Asia (Khanazarov et al., 2009), though they may be found at extremes of 450 or 2000 masl (Zhukovsky, 1971). The northern Tian-Shan Region, covering around 30,000 ha, is mainly distributed in northwestern Kyrgyzstan on foothills (800–1600 masl) of the Fergana and Chingir-Tash mountains (Kayimov et al., 2011). Talebi et al. (2012) linked high genetic diversity among P. vera to different geographic distribution zones, highlighting the importance of these natural populations as genetic resources.

The domestication and breeding history of pistachio differs fundamentally from that of annual grain crops, as it illustrates very different concepts of land tenure, trans-generational forethought, and complex cultivation practices, including grafting and artificial pollination (Miller, 2008). Historically, pistachio cultivation has primarily relied on grafting, suggesting that, as with many Eurasian tree crops, domestication resulted from genetically locking hybrids or favored individuals in place (Zohary and Spiegel-Roy, 1975). The most commercially available pistachio cultivars worldwide are the result of the selection of chance seedlings (Parfitt et al., 1995), with prominent varieties such as Kerman, Ohadi, and Siirt obtained through selection breeding programs in their respective countries. The delay in dispersal of pistachios may be due to its reproductive biology, notably a need for mating pairs, large closely clustered populations to facilitate wind pollination, and knowledge of grafting and artificial propagation for maintaining clonal variants (Vargas et al., 2008). As a dioecious and wind-pollinated plant, early pistachio cultivators would have had to transplant a full population of multiple male and female plants into a suitable new area rather than simply collect a few seeds. Today, only 300,000 ha of Central Asia is covered with wild P. vera (Chernova, 2004), but a few centuries ago, much larger territories were populated with pistachio groves. In the 1950s, Whitehouse (1957) reported that "the pistachio is seen everywhere from mountainous Bokhara through the region of Samarkand and the whole of Fergana up to the western Tian-Shan", indicating the historical abundance of these genetic resources in the region.

Kyrgyzstan possesses substantial pistachio genetic resources within its natural wild populations, yet these remain largely unexploited for systematic improvement. Despite ranking tenth globally in pistachio production with approximately 921 tonnes in 2023 (FAOSTAT, 2023), Kyrgyzstan does not have any registered improved pistachio varieties. Approximately 40,000 hectares of pistachio trees grow from seed without grafting across the mountainous regions of the country, representing a vast untapped genetic resource. Selection breeding, being the oldest and most straightforward breeding method, offers a practical approach for genetic improvement in countries where sophisticated breeding infrastructure may be limited. The natural populations in Kyrgyzstan's Tian-Shan region provide an ideal setting for selection programs, as they harbor significant phenotypic and genetic variation shaped by diverse environmental conditions. Therefore, the aim of this study was to evaluate the morphological and pomological diversity among naturally growing pistachio genotypes in the mountainous regions near Kuduk-Sai village, Nooken district, Jalal-Abad region, Kyrgyzstan, and to identify promising superior genotypes using multivariate statistical approaches including a Composite Performance Index (CPI) based on key traits such as nut size, inshell weight, split ratio, shell thickness, and kernel yield.

Materials and Methods
Study Area and Environmental Conditions

This research was conducted during 2021 and 2024, with a two-year periodicity observed between the study years. The study site is located at coordinates 41.217889°N, 72.238854°E, at an elevation ranging from 634 to 685 meters above sea level, encompassing a total area of 80 hectares (Fig. 1). The research was carried out in a naturally established pistachio (Pistacia vera L.) forest, where trees originated from spontaneously dispersed seeds through natural vectors including birds and grazing animals such as goats. Trees of varying ages (5–150 years) were sampled to maximize genetic diversity across natural regeneration cohorts. Age estimation was based on trunk diameter and stand history records. While tree age and microsite conditions may influence phenotypic expression, this study focused on preliminary screening of genetic variability rather than isolating age or environmental effects. The region is characterized by arid climatic conditions with annual precipitation ranging from 300 to 350 mm. The soil exhibits high calcium carbonate content, typical of calcareous substrates in semi-arid regions. No other tree species naturally coexist in this monospecific pistachio forest ecosystem. The irregular spatial distribution of trees, with inter-tree distances varying from 1.0 to 30 meters, reflects the stochastic nature of seed dispersal and establishment, resulting in an average stand density of approximately 200 trees per hectare.

Fig. 1:

An image from the area under investigation (original)

Stand Characteristics and Management

The forest contains approximately 17,500 pistachio trees, with a sex ratio of approximately 60% female and 40% male trees. Among the female population, 80-90% exhibited fruit-bearing capacity during the study period. The trees receive no anthropogenic management interventions and develop entirely under natural conditions without supplementary irrigation, fertilization, pruning, or pest management practices. The forest serves dual ecological and socioeconomic functions, providing erosion control on steep slopes while supporting local livelihoods through fruit harvest.

Selection Methodology and Pomological Evaluation

A comprehensive field survey was conducted across the entire 80-hectare study area to identify superior genotypes based on fruit-bearing performance and preliminary visual assessment. Male trees were systematically excluded from the selection process. Among fruit-bearing female trees, 169 promising genotypes were identified and marked in the initial year (2021) based on apparent fruit quality characteristics. From these marked genotypes, 50 representative fruit samples were collected for detailed laboratory analysis. Pomological evaluation included comprehensive morphometric and quality assessments of the following parameters: fruit length, thickness, and width (mm); in-shell weight (g); 100-nuts weight (g); kernel color characteristics; splitting percentage (%); shell thickness (mm); kernel weight (g); and kernel ratio (%). These quantitative and qualitative evaluations were designed to identify genotypes with superior characteristics suitable for potential domestication or germplasm conservation, thereby contributing to the genetic improvement of pistachio cultivation in the region while preserving valuable genetic diversity present in this naturally evolved population (Fig. 2).

Fig. 2:

Appearances of 19 promising pistachio genotypes (original images)

Data Standardization and Trait Selection

A total of 169 pistachio (Pistacia vera L.) genotypes were evaluated based on five morphological and pomological traits: nut size (mm), in-shell weight (g), split ratio (%), shell thickness (mm), and kernel yield (%). Prior to multivariate analysis, all quantitative data were standardized using the z-score transformation to eliminate the effects of differing measurement scales and units: Zij=(XijX¯j)/Sj{Z_{ij}} = \left( {{X_{ij}} - {{\bar X}_j}} \right)/{S_j} where Zij is the standardized value of trait j for genotype i; Xij is the observed value; X¯j{{\bar X}_j} and Sj are the mean and standard deviation of trait j, respectively.

This standardization ensured that each variable contributed equally to the subsequent multivariate analyses.

Composite Performance Index

A Composite Performance Index (CPI) was developed to integrate multiple trait dimensions into a single quantitative score for ranking genotypes. The CPI was calculated as a weighted linear combination of standardized trait values, expressed as follows: CPIi=( 0.10×ZSize +(0.15×ZIn-shell weight )+(0.20 ×ZSplit Ratio  )(0.25×ZShell Thickness )+(0.30×ZKernel Yield )CP{I_i} = \left( {0.10 \times {Z_{{\rm{Size }}}} + \left( {0.15 \times {Z_{{\rm{In - shell weight }}}}} \right) + (0.20} \right.\;\left. { \times {Z_{{\rm{Split Ratio }}}}} \right) - \left( {0.25 \times {Z_{{\rm{Shell Thickness }}}}} \right) + \left( {0.30 \times {Z_{{\rm{Kernel Yield }}}}} \right) where Z represents the standardized value (Z-score) of each trait, calculated as Z = (X - μ)/σ, with X being the observed value, μ the population mean, and σ the standard deviation. Weights were assigned based on economic importance and breeding objectives. Kernel yield received the highest weight (0.30) as it directly determines commercial productivity and revenue. Shell thickness was assigned a negative weight (-0.25) due to its inverse relationship with kernel recovery and processing efficiency. Split ratio (0.20) was weighted for its impact on marketability of whole kernels, while in-shell weight (0.15) and nut size (0.10) received moderate weights as secondary selection criteria. These weights reflect current market preferences and industry standards for, where kernel quality and yield predominate over physical dimensions. Each genotype's CPI score was subsequently ranked from highest to lowest, facilitating the identification of high-yielding and morphologically desirable genotypes.

Statistical Analysis

Data collected from pomological evaluations were subjected to comprehensive statistical analysis to assess genetic variability and identify superior genotypes. Descriptive statistics including mean, standard deviation, and standard error were calculated for all measured traits. Analysis of variance (ANOVA) was performed to determine significant differences among genotypes for size, in-shell weight, split ratio, shell thickness, kernel yield, and kernel weight, with statistical significance determined at p<0.05 level. Duncan's Multiple Range Test was employed as post-hoc analysis to distinguish statistically different genotypes at α = 0.05 significance level, with each genotype evaluated using three replications to ensure data reliability. To facilitate multi-trait selection, z-score standardization was applied to normalize all traits using the formula Z = (X - μ) / σ where X represents individual observation, μ is the population mean, and σ is the standard deviation. Genotypes were ranked according to their CPI values, with higher positive values indicating superior overall performance. Pearson's correlation network analysis was performed using R software, and all graphical outputs were prepared using Origin software version 2021.

Results
Composite Performance Index and Genotype Ranking

The Composite Performance Index (CPI) effectively integrated multiple trait dimensions into a single quantitative metric, enabling comprehensive evaluation and ranking of all 169 genotypes (Tab. 1). CPI values ranged from -1.91 (KS-011) to 1.40 (KS-168), with a mean of 0.00 and standard deviation of 0.58. Based on CPI rankings, 19 superior genotypes representing the top 11.2% of the evaluated germplasm were systematically identified for advanced evaluation and potential cultivar development. Genotype KS-168 emerged as the highest-ranking selection (CPI = 1.40, Rank 1), distinguished by exceptional kernel yield (58.80%, z-score = 2.65), moderate nut size (1.56 cm), remarkably thin shell (0.53 mm, z-score = -1.65), and acceptable split ratio (87.50%). KS-103 ranked second (CPI = 1.32, Rank 2) with outstanding kernel yield (56.48%, z-score = 2.01), thin shell (0.53 mm), and impressive nut size (1.77 cm, z-score = 1.46), despite a moderately lower split ratio (77.40%). KS-132 held the third position (CPI = 1.19, Rank 3) with the highest kernel yield among top performers (57.00%, z-score = 2.15), exceptionally thin shell (0.50 mm, z-score = -2.01), and excellent split ratio (93.50%). The CPI scores theoretically range from negative to positive values depending on the distribution of standardized traits, with higher positive values indicating superior overall performance. In this study, observed CPI values ranged from –1.48 to 1.40, with the top-performing genotypes exceeding 0.88. The top ten geno-types consistently demonstrated CPI values exceeding 0.88, indicating substantial deviation above the population mean across weighted traits. KS-076 (Rank 4), KS-142 (Rank 5), and KS-143 (Rank 6) exhibited kernel yields above 54% with thin shells and excellent split ratios. KS-120 (CPI = 1.00, Rank 7) was particularly noteworthy for its extraordinarily large nut size (2.12 cm, z-score = 2.81), the highest among all 169 genotypes, combined with perfect split ratio (100%). KS-126 (Rank 8) demonstrated exceptional kernel yield (55.34%) paired with the thinnest shell among top performers (0.52 mm, z-score = -1.77). KS-157 (Rank 9) achieved the second-highest kernel yield overall (56.52%, z-score = 2.02), while KS-128 (Rank 10) displayed impressive kernel yield (52.57%) with outstanding split ratio (96.70%). Among genotypes ranked 11-19, KS-096, KS-141, KS-138, and KS-058 demonstrated high kernel yields exceeding 50% with favorable split characteristics. KS-058 (Rank 13) was particularly remarkable for achieving kernel yield of 55.45% despite compact size (1.46 cm), representing exceptional kernel-to-shell efficiency. KS-088 (Rank 17) possessed the thinnest shell among all genotypes (0.46 mm, z-score = -2.50), representing a critical genetic resource for shell thickness reduction. KS-133 (Rank 18) and KS-167 (Rank 19) offered balanced commercial attributes with large nut sizes (1.76 cm and 1.86 cm, respectively) and solid kernel yields above 51%. Standardized z-scores provided critical insights into trait performance. Shell thickness z-scores ranged from -2.50 to 2.71, with 16 of the 19 selected genotypes exhibiting negative values, confirming effective selection for thin-shelled characteristics. Kernel yield z-scores spanned from -3.13 (KS-110) to 2.65 (KS-168), representing a 20.92 percentage point difference. A total of 48 genotypes achieved perfect split ratios (100%).

Tab. 1:

Physical and chemical characteristics, Z-score values, and rankings of pistachio (Pistacia vera L.) genotypes

GenotypesSizeIn-shell WeightSplit RatioShell ThicknessKernel YieldCPIRankZ_SizeZ_In-shell WeightZ_Split RatioZ_Shell ThicknessZ_Kernel Yield
KS-0011.90.9100.00.844.7-0.2106.01.82.10.82.0-1.3
KS-0021.50.778.00.751.20.178.00.50.7-0.20.80.6
KS-0031.30.66.60.751.6-0.7144.0-0.5-0.4-3.30.30.7
KS-0041.00.553.30.747.7-0.8151.0-1.5-1.4-1.30.4-0.4
KS-0052.00.643.70.845.1-0.7145.02.30.0-1.71.0-1.1
KS-0060.90.4100.00.652.10.093.0-2.1-2.20.8-0.40.8
KS-0071.20.663.60.840.3-1.7168.0-0.9-0.9-0.82.1-2.5
KS-0081.30.667.70.651.10.085.0-0.5-0.5-0.6-0.60.5
KS-0091.10.5100.00.746.9-0.4135.0-1.3-1.20.80.3-0.7
KS-0101.40.682.30.751.60.275.00.1-0.20.00.00.7
KS-0111.20.657.50.939.0-1.9169.0-0.8-0.6-1.12.7-2.8
KS-0120.80.4100.00.644.5-0.7147.0-2.4-2.50.8-0.4-1.3
KS-0131.00.597.00.648.1-0.2107.0-1.4-1.50.6-0.4-0.3
KS-0141.00.596.80.647.9-0.3112.0-1.5-1.60.6-0.4-0.4
KS-0151.20.583.80.743.2-0.8148.0-1.0-1.00.10.0-1.7
KS-0161.20.697.00.746.6-0.3125.0-0.8-0.80.60.2-0.7
KS-0171.20.6100.00.649.20.176.0-0.7-0.80.8-0.70.0
KS-0181.00.496.80.647.5-0.3127.0-1.7-2.00.6-0.6-0.5
KS-0191.40.694.20.749.6-0.199.0-0.2-0.50.50.70.1
KS-0201.20.595.00.746.6-0.6140.0-0.9-1.10.60.8-0.7
KS-0211.00.4100.00.647.3-0.3116.0-1.6-2.00.8-0.7-0.5
KS-0221.30.693.70.746.6-0.4132.0-0.6-0.80.50.4-0.7
KS-0231.20.668.70.742.8-1.1158.0-0.8-0.8-0.60.8-1.8
KS-0241.00.575.00.741.0-1.3163.0-1.8-1.9-0.30.3-2.3
KS-0251.40.797.00.745.7-0.3120.00.20.10.60.7-1.0
KS-0261.40.696.00.652.90.540.00.0-0.30.6-0.61.0
KS-0271.40.678.00.651.90.361.00.0-0.2-0.2-0.60.8
KS-0281.30.693.70.843.8-0.9152.0-0.3-0.50.51.6-1.5
KS-0291.30.690.30.750.70.179.0-0.5-0.70.4-0.20.4
KS-0301.40.6100.00.747.7-0.2102.00.00.00.80.8-0.4
KS-0311.20.639.30.842.0-1.4165.0-0.8-0.4-1.91.3-2.0
KS-0320.90.481.00.649.9-0.4131.0-2.0-2.2-0.1-0.30.2
KS-0331.10.590.60.746.4-0.6143.0-1.0-1.20.40.5-0.8
KS-0341.40.690.90.747.5-0.3113.0-0.2-0.10.40.7-0.5
KS-0351.50.793.70.749.30.089.00.60.40.50.90.0
KS-0361.10.693.50.749.0-0.3124.0-1.3-0.90.50.5-0.1
KS-0371.40.665.60.749.3-0.3114.0-0.1-0.7-0.70.00.0
KS-0381.30.695.00.745.4-0.5139.0-0.6-0.30.60.9-1.1
KS-0391.30.6100.00.747.1-0.2105.0-0.5-0.40.80.3-0.6
KS-0401.20.697.00.841.0-1.2162.0-0.9-0.90.61.5-2.3
KS-0411.20.597.00.746.7-0.4134.0-1.0-1.30.60.2-0.7
KS-0421.30.693.50.750.10.090.0-0.3-0.70.50.20.3
KS-0431.00.5100.00.750.0-0.3119.0-1.7-1.90.80.30.2
KS-0441.00.588.20.843.2-1.1160.0-1.4-1.50.31.0-1.7
KS-0451.40.75.00.848.5-1.1161.00.10.6-3.42.0-0.2
KS-0461.10.590.30.744.6-0.8150.0-1.1-1.20.40.8-1.3
KS-0471.00.5100.00.653.50.448.0-1.7-1.90.8-1.41.2
KS-0481.00.5100.00.654.60.455.0-1.7-1.90.8-0.81.5
KS-0491.30.697.00.850.6-0.2108.0-0.3-0.50.61.50.4
KS-0500.80.495.00.749.8-0.4133.0-2.2-2.20.60.00.2
KS-0512.20.915.00.750.20.183.03.32.2-2.90.40.3
KS-0521.30.6100.00.751.40.268.0-0.5-0.40.80.00.6
KS-0531.50.697.00.746.6-0.3109.00.2-0.30.60.5-0.7
KS-0541.50.791.00.750.30.270.00.40.20.40.20.3
KS-0551.60.852.00.847.4-0.6141.00.81.0-1.31.5-0.5
KS-0561.30.643.70.745.3-0.9153.0-0.4-0.3-1.70.7-1.1
KS-0571.30.697.00.649.80.269.0-0.4-0.40.6-0.60.2
KS-0581.50.690.30.655.50.813.00.20.00.4-0.91.7
KS-0591.00.584.30.743.2-1.1159.0-1.7-1.90.10.5-1.7
KS-0601.30.681.80.748.4-0.4128.0-0.4-0.20.00.9-0.2
KS-0610.90.596.80.845.2-1.0156.0-1.9-1.80.61.3-1.1
KS-0621.80.893.50.845.2-0.3121.01.41.30.51.6-1.1
KS-0631.50.797.00.844.0-0.7146.00.40.40.62.1-1.5
KS-0641.30.795.00.747.1-0.2103.0-0.30.20.60.5-0.6
KS-0651.40.660.00.747.2-0.6142.00.2-0.2-1.00.8-0.6
KS-0661.00.596.80.751.9-0.196.0-1.8-1.60.6-0.10.8
KS-0671.00.5100.00.748.0-0.5137.0-1.7-1.40.80.5-0.3
KS-0681.30.650.00.745.6-0.9154.0-0.5-0.8-1.40.7-1.0
KS-0691.20.648.40.750.8-0.4130.0-0.7-0.3-1.50.40.4
KS-0701.20.670.00.842.2-1.4164.0-0.7-0.4-0.52.1-1.9
KS-0711.10.590.30.846.1-0.8149.0-1.3-1.00.41.1-0.9
KS-0721.40.694.00.744.7-0.5138.0-0.1-0.30.50.7-1.2
KS-0731.10.583.80.648.5-0.3123.0-1.3-1.50.1-0.3-0.2
KS-0741.20.640.00.746.4-0.9155.0-0.9-0.2-1.90.8-0.8
KS-0751.30.648.30.842.0-1.5167.0-0.4-0.1-1.52.1-2.0
KS-0761.60.793.50.554.91.14.00.60.40.5-1.51.6
KS-0771.40.693.70.750.30.274.00.1-0.10.50.00.3
KS-0781.70.716.00.546.4-0.2101.01.10.8-2.9-1.7-0.8
KS-0791.20.650.00.647.2-0.3126.0-0.9-0.6-1.4-1.2-0.6
KS-0801.60.775.00.654.40.720.00.60.4-0.3-1.01.4
KS-0811.50.797.00.650.50.453.00.30.10.6-0.30.4
KS-0821.60.77.00.751.9-0.3110.00.80.7-3.30.00.8
KS-0831.40.713.00.847.0-1.0157.0-0.20.4-3.01.0-0.6
KS-0841.50.791.70.650.90.446.00.20.40.4-0.40.5
KS-0851.60.784.00.746.6-0.197.00.60.70.10.3-0.7
KS-0861.80.897.00.854.30.722.01.71.70.61.01.4
KS-0871.50.759.00.749.20.095.00.40.6-1.0-0.10.0
KS-0881.20.595.00.552.60.817.0-0.9-1.00.6-2.50.9
KS-0891.60.790.60.646.70.358.00.80.70.4-1.2-0.7
KS-0901.40.796.90.652.30.721.00.10.20.6-1.30.9
KS-0911.40.7100.00.548.50.534.00.20.20.8-1.7-0.2
KS-0921.70.8100.00.744.20.088.01.31.40.80.3-1.4
KS-0931.80.884.30.648.20.449.01.51.60.1-0.3-0.3
KS-0941.50.781.20.746.50.094.00.20.8-0.1-0.2-0.8
KS-0951.50.790.90.647.60.451.00.20.60.4-1.3-0.5
KS-0961.70.893.70.550.60.911.01.01.20.5-1.50.4
KS-0971.80.868.70.747.60.184.01.51.3-0.60.2-0.5
KS-0981.60.793.30.549.00.724.00.60.30.5-2.0-0.1
KS-0991.50.775.00.652.90.723.00.30.4-0.3-1.41.0
KS-1001.40.796.00.743.5-0.3122.0-0.10.40.60.0-1.6
KS-1011.60.780.60.552.10.816.00.70.7-0.1-1.50.8
KS-1021.40.696.80.551.00.725.00.00.00.6-1.50.5
KS-1031.80.877.40.556.51.32.01.51.3-0.2-1.72.0
KS-1041.50.730.30.650.10.177.00.50.8-2.3-1.30.3
KS-1051.50.793.30.846.9-0.3118.00.50.40.51.3-0.6
KS-1061.90.973.30.749.80.535.01.92.2-0.4-0.20.2
KS-1071.60.751.50.752.20.272.00.60.9-1.30.00.8
KS-1081.70.710.00.653.30.181.01.20.6-3.2-0.61.1
KS-1091.30.680.60.749.80.091.0-0.3-0.2-0.1-0.10.2
KS-1101.30.6100.00.937.9-1.5166.0-0.4-0.40.82.3-3.1
KS-1111.50.756.60.751.10.087.00.20.3-1.1-0.10.5
KS-1121.60.793.50.752.60.628.00.80.50.5-0.21.0
KS-1131.40.796.80.651.20.538.00.10.20.6-0.80.6
KS-1141.40.7100.00.652.80.629.00.10.40.8-0.31.0
KS-1151.60.796.70.649.20.539.00.60.50.6-1.00.0
KS-1161.70.893.70.648.70.543.01.21.20.5-0.6-0.1
KS-1171.50.7100.00.749.50.364.00.50.50.80.20.1
KS-1181.70.878.00.648.60.362.01.01.1-0.2-0.4-0.2
KS-1191.60.797.00.648.20.454.00.70.60.6-0.6-0.3
KS-1202.10.9100.00.649.91.07.02.82.10.8-0.80.2
KS-1211.70.883.80.651.10.726.01.31.40.1-0.60.5
KS-1221.70.894.00.748.80.457.01.01.50.50.2-0.1
KS-1231.40.769.60.652.20.536.00.10.8-0.6-1.00.8
KS-1241.40.7100.00.749.70.365.00.10.20.8-0.10.1
KS-1251.30.6100.00.650.60.544.0-0.5-0.70.8-1.40.4
KS-1261.40.690.00.555.31.08.00.1-0.20.3-1.81.7
KS-1271.40.7100.00.750.10.271.00.00.20.80.30.3
KS-1281.60.796.70.652.60.910.00.60.70.6-1.40.9
KS-1291.60.741.90.650.70.273.00.90.9-1.8-0.70.4
KS-1301.30.6100.00.652.30.541.0-0.4-0.40.8-0.80.9
KS-1311.30.6100.00.651.40.537.0-0.4-0.30.8-1.20.6
KS-1321.40.693.50.557.01.23.0-0.2-0.20.5-2.02.2
KS-1331.80.8100.00.651.20.818.01.41.50.8-0.30.5
KS-1341.60.896.70.752.50.727.00.81.00.6-0.10.9
KS-1351.30.696.60.651.00.456.0-0.20.00.6-0.40.5
KS-1361.60.797.00.651.40.542.00.70.40.6-0.30.6
KS-1371.70.8100.00.850.10.267.01.01.40.81.30.2
KS-1381.60.791.10.653.90.814.00.80.40.4-0.91.3
KS-1391.70.8100.00.844.3-0.4136.01.21.20.82.0-1.4
KS-1401.40.793.90.751.10.363.0-0.10.30.50.00.5
KS-1411.40.794.00.653.80.912.00.00.20.5-1.41.3
KS-1421.70.894.00.654.61.05.01.01.10.5-0.81.5
KS-1431.60.791.60.654.41.06.00.60.70.4-1.31.4
KS-1441.80.915.60.748.3-0.3117.01.72.4-2.90.7-0.3
KS-1451.20.696.60.747.7-0.3111.0-0.7-0.30.60.5-0.4
KS-1461.40.762.50.649.60.086.00.00.4-0.9-0.40.1
KS-1471.40.778.00.748.50.092.00.10.2-0.2-0.2-0.2
KS-1481.60.7100.00.750.30.445.00.70.70.80.00.3
KS-1491.50.7100.00.748.60.266.00.30.70.80.0-0.2
KS-1501.90.7100.00.749.90.360.01.90.40.80.70.2
KS-1511.40.687.50.653.40.447.0-0.2-0.20.2-0.31.2
KS-1521.20.6100.00.652.10.633.0-0.7-0.40.8-1.20.8
KS-1531.40.697.00.651.60.632.0-0.20.00.6-1.00.7
KS-1541.30.693.30.650.40.452.0-0.50.00.5-0.90.3
KS-1551.40.693.70.651.40.630.0-0.2-0.20.5-1.40.6
KS-1561.30.683.80.649.10.182.0-0.5-0.40.1-0.70.0
KS-1571.50.794.00.656.51.09.00.40.30.5-0.72.0
KS-1581.80.762.50.751.40.180.01.60.4-0.90.70.6
KS-1591.40.790.60.849.0-0.2100.00.20.80.41.4-0.1
KS-1601.50.745.20.755.90.450.00.50.5-1.6-0.11.9
KS-1611.50.837.50.749.2-0.4129.00.41.3-2.00.90.0
KS-1621.20.6100.00.851.4-0.2104.0-0.7-0.80.81.40.6
KS-1631.60.7100.00.852.10.359.00.70.90.81.10.8
KS-1641.60.862.50.749.1-0.198.00.71.1-0.90.70.0
KS-1651.50.880.00.753.70.631.00.51.0-0.1-0.11.2
KS-1661.30.774.10.851.4-0.3115.0-0.20.3-0.41.60.6
KS-1671.90.874.10.652.50.819.01.81.8-0.4-0.40.9
KS-1681.60.787.50.558.81.41.00.60.50.2-1.72.7
KS-1691.60.796.80.652.40.815.00.80.80.6-0.80.9
Performance Characteristics of Selected Elite Genotypes

The 19 selected genotypes were classified into three distinct performance groups based on their phenotypic profiles and selection criteria (Tab. 2). The high-performing group comprised 11 genotypes (KS-076, KS-096, KS-103, KS-116, KS-121, KS-126, KS-133, KS-142, KS-143, KS-168, KS-169) characterized by superior kernel yield, balanced morphological traits, and excellent overall performance. These genotypes exhibited CPI values ranging from 0.50 to 1.40 and kernel yields from 48.68% to 58.80%. The efficient shell type group included four genotypes (KS-058, KS-086, KS-120, KS-128) distinguished by thin shells, compact dimensions, and high kernel-to-shell ratios, with CPI values between 0.73 and 1.06. The balanced type group contained four genotypes (KS-001, KS-128, KS-142, KS-157, KS-167) demonstrating uniform performance across multiple traits without extreme values in any single parameter, with CPI values ranging from -0.22 to 1.01. Notably, genotype KS-001, despite ranking 106th overall (CPI = -0.22), was selected for its exceptional size (1.86 cm, z-score = 1.82) and in-shell weight (0.86 g, z-score = 2.09), combined with perfect split ratio (100%), representing a valuable genetic resource for improving nut size in breeding programs. Conversely, genotype KS-093 (CPI = 0.40, Rank 49) was classified as an elite performer due to its consistently above-average performance across all evaluated traits, demonstrating balanced genetic potential without extreme deficiencies.

Tab. 2:

Pistachio (Pistacia vera L.) genotypes: Comprehesive evaluation of CPI, rankings, physical traits, selection reasons, and classification groups

GenotypesCPIRankSizeIn-shell WeightSplit RatioShell-ThicknessKernel YieldSelection ReasonGroup
KS-001-0.2106.01.90.9100.00.844.7High kernel yield and shell balanceBalanced type
KS-0580.813.01.50.690.30.655.5Compact size and efficient yieldEfficient shell type
KS-0761.14.01.60.793.50.554.9Excellent split ratioHigh-performing
KS-0860.722.01.80.897.00.854.3Thin shell and high kernel percentageEfficient shell type
KS-0930.449.01.80.884.30.648.2Top performer across traitsElite performer
KS-0960.911.01.70.893.70.550.6Excellent kernel yieldHigh-performing
KS-1031.32.01.80.877.40.556.5Thin shell and high kernel yieldHigh-performing
KS-1160.543.01.70.893.70.648.7High in-shell weight and kernel yieldHigh-performing
KS-1201.07.02.10.9100.00.649.9Thin shell and large kernelEfficient shell type
KS-1210.726.01.70.883.80.651.1High kernel yield with acceptable shellHigh-performing
KS-1261.08.01.40.690.00.555.3High split ratio and kernel yieldHigh-performing
KS-1280.910.01.60.796.70.652.6Thin shell with balanced kernel ratioBalanced type
KS-1330.818.01.80.8100.00.651.2High yield and kernel ratioHigh-performing
KS-1421.05.01.70.894.00.654.6Balanced performance across all traitsBalanced type
KS-1431.06.01.60.791.60.654.4High in-shell weight and kernel yieldHigh-performing
KS-1571.09.01.50.794.00.656.5Large nut size and good yieldBalanced type
KS-1670.819.01.90.874.10.652.5Large nut and uniform morphologyBalanced type
KS-1681.41.01.60.787.50.558.8High kernel yield and large sizeHigh-performing
KS-1690.815.01.60.796.80.652.4High split ratio and light shellHigh-performing
Detailed Pomological Characteristics of Superior Genotypes

Analysis of variance revealed highly significant differences (p < 0.05) among the 169 pistachio genotypes for all evaluated pomological traits, confirming substantial genetic diversity within the natural population (Tab. 3). Kernel yield exhibited the highest F-value (292.82), followed by shell thickness (80.72), split ratio (78.02), size (37.82), in-shell weight (29.12), and kernel weight (14.40). Duncan's Multiple Range Test identified statistically significant differences among the 19 selected genotypes for all evaluated traits (Tab. 4). Kernel yield ranged from 35.36% (KS-116) to 61.18% (KS-121). Genotypes KS-121 (61.18%) and KS-128 (57.14%) comprised the highest performing group, significantly superior to all others. The second group included 11 genotypes with kernel yields ranging from 41.31% to 49.66%, while groups and represented genotypes with progressively lower kernel percentages. Interestingly, KS-076 exhibited an anomalous kernel yield value (0.50%). Nut size measurements revealed three distinct groups. Genotype KS-001 (1.49 cm) formed a unique group, significantly larger than all other selections. The majority of genotypes (17 selections) clustered in group with sizes ranging from 0.99 to 1.37 cm, while KS-142 (0.85 cm) constituted the smallest group. In-shell weight showed limited variation, with 18 genotypes forming a homogeneous group ranging from 0.483 to 0.670 g, and only KS-001 (0.730 g) significantly heavier. Split ratio demonstrated remarkable uniformity among selected genotypes, with 18 selections forming a single homogeneous group ranging from 74.67 to 88.67%, while only KS-096 (54.67 %) exhibited significantly lower splitting percentage. Shell thickness measurements produced three distinct groups: KS-167 demonstrated the thickest shell (0.943 mm), 14 genotypes showed intermediate thickness (0.6470.840 mm), and KS-103 possessed the thinnest shell (0.593 mm). Kernel weight ranged from 0.203 to 0.337 g, with KS-168 (0.203_g) significantly lighter than the remaining 18 genotypes (0.2270.337 g).

Tab. 3

Analysis of variance (ANOVA) for physical and yield-related traits of pistachio (Pistacia vera L.) genotypes.

TraitSourceDFSSMSF-valueSignificance
SizeBetween Genotypes181.420.078737.82< 0.05**
Within Genotypes (Error)380.080.0021--
In-shell weightBetween Genotypes180.210.011629.12< 0.05**
Within Genotypes (Error)380.020.0004--
Split RatioBetween Genotypes182661.05147.836378.02< 0.05**
Within Genotypes (Error)3872.001.8947--
Shell thicknessBetween Genotypes180.420.023580.72< 0.05**
Within Genotypes (Error)380.010.0003--
Kernel yieldBetween Genotypes188054.66447.4812292.82< 0.05**
Within Genotypes (Error)3858.071.5282--
Kernel WeightBetween Genotypes180.060.003314.40< 0.05**
Within Genotypes (Error)380.010.0002--

DF: Degrees of freedom; SS: Sum of squares; MS: Mean square; Error: Residual variation within genotypes (experimental error representing unexplained variance after accounting for genotype effects);

**

Significant at p < 0.05.

Tab. 4

Physical and yield-related traits of pistachio (Pistacia vera L.) genotypes with mean ± standard error (SE)

GenotypeSize (cm)In-shell Weight (g)Split Ratio (%)Shell Thickness (mm)Kernel Yield (%)Kernel Weight (g)
KS-1211.27b± 0.000.527b± 0.00681.33a± 0.660.800b± 0.01561.18a± 0.600.300a± 0.011
KS-1281.25b± 0.010.547b± 0.00684.67a± 0.600.663b± 0.00657.14a± 0.590.320a± 0.015
KS-1201.14b± 0.010.483b± 0.01476.67a± 0.630.660b± 0.01149.66b± 0.870.273a± 0.006
KS-1671.35b± 0.000.513b± 0.00682.67a± 0.660.943a± 0.00649.33b± 0.660.257a± 0.003
KS-0931.37b± 0.010.640b± 0.01576.67a± 0.690.793b± 0.01748.83b± 0.440.303a± 0.006
KS-1031.37b± 0.010.600b± 0.01174.67a± 0.660.593c± 0.00648.51b± 0.510.303a± 0.003
KS-0961.28b± 0.010.613b± 0.00654.67b± 0.650.760b± 0.01047.40b± 0.700.293a± 0.006
KS-0860.99b± 0.000.563b± 0.00677.33a± 0.600.840b± 0.00546.79b± 0.910.283a± 0.006
KS-1571.23b± 0.010.547b± 0.00379.33a± 0.630.687b± 0.00646.76b± 0.900.257a± 0.003
KS-0011.49a± 0.080.730a± 0.02880.67a± 1.760.813b± 0.01744.67b± 0.880.337a± 0.008
KS-1691.13b± 0.010.573b± 0.00676.00a± 1.150.660b± 0.01143.67b± 0.660.260a± 0.000
KS-0581.30b± 0.010.647b± 0.00884.67a± 0.660.820b± 0.01143.31b± 0.650.293a± 0.003
KS-1420.85c± 0.010.633b± 0.00685.33a± 0.670.777b± 0.00842.40b± 0.870.277a± 0.012
KS-1431.13b± 0.010.647b± 0.00882.67a± 0.640.647b± 0.00642.00b± 1.150.250a± 0.010
KS-1331.10b± 0.050.613b± 0.00682.67a± 0.630.663b± 0.00641.31b± 0.650.250a± 0.010
KS-1261.19b± 0.010.670b± 0.01178.67a± 0.670.717b± 0.00637.10c± 0.580.237a± 0.006
KS-1681.07b± 0.000.543b± 0.00681.33a± 0.630.683b± 0.00835.79c± 0.400.203b± 0.003
KS-1161.04b± 0.020.613b± 0.02088.67a± 0.610.663b± 0.00635.36c± 0.670.227a± 0.017
KS-0760.97b± 0.010.533b± 0.00377.33a± 0.620.793b± 0.0060.50d± 0.000.273a± 0.008

Duncan's Multiple Range Test was performed at a = 0.05 significance level

  • Means followed by the same letter (superscript) are not significantly different

  • Different letters indicate statistically significant differences between genotypes

  • Each value represents the mean of 3 replications

Phenological Stages and Harvest Dates of Pistachio Genotypes in 2024

The phenological assessment of the 19 pistachio genotypes revealed considerable variation in the timing of bud burst, flowering stages, and harvest dates (Tab. 5). Bud burst occurred between 25 March (KS-103 and KS-126) and 3 April (KS-116 and KS-157), indicating early to late vegetative initiation among genotypes. The first bloom was observed between 31 March (KS-103) and 10 April (KS-157), while full bloom occurred from 3 April (KS-103) to 14 April (KS-157). The end of bloom ranged from 7 April (KS-103) to 20 April (KS-157), resulting in flowering durations of 6 to 10 days, reflecting differences in flowering intensity and synchronization. Harvest dates varied notably, spanning from 1 August (KS-093) to 15 August (KS-116 and KS-168), demonstrating genotype-dependent variation in fruit maturation. Overall, KS-103 consistently exhibited the earliest phenological stages, whereas KS-116 and KS-168 represented the latest genotypes in terms of harvest.

Tab. 5

Phenological stages and harvest dates of pistachio (Pistacia vera L.) genotypes in 2024

GenotypesBud BurstFirst BloomFull BloomEnd of BloomFlowering DurationHarvest Date
KS-00127-Mar-243-Apr-248-Apr-2414-Apr-249 days8-Aug-24
KS-05829-Mar-244-Apr-248-Apr-2412-Apr-248 days3-Aug-24
KS-07626-Mar-242-Apr-245-Apr-248-Apr-246 days14-Aug-24
KS-0862-Apr-248-Apr-2412-Apr-2416-Apr-248 days9-Aug-24
KS-09328-Mar-245-Apr-249-Apr-2414-Apr-249 days1-Aug-24
KS-09631-Mar-246-Apr-2410-Apr-2416-Apr-2410 days11-Aug-24
KS-10325-Mar-2431-Mar-243-Apr-247-Apr-247 days6-Aug-24
KS-1163-Apr-249-Apr-2413-Apr-2419-Apr-2410 days15-Aug-24
KS-12027-Mar-243-Apr-246-Apr-2410-Apr-247 days4-Aug-24
KS-12130-Mar-245-Apr-249-Apr-2415-Apr-2410 days13-Aug-24
KS-12625-Mar-241-Apr-244-Apr-249-Apr-248 days7-Aug-24
KS-1281-Apr-247-Apr-2411-Apr-2417-Apr-2410 days10-Aug-24
KS-13326-Mar-241-Apr-245-Apr-2410-Apr-249 days2-Aug-24
KS-14230-Mar-246-Apr-249-Apr-2415-Apr-249 days12-Aug-24
KS-14328-Mar-244-Apr-247-Apr-2411-Apr-247 days5-Aug-24
KS-1573-Apr-2410-Apr-2414-Apr-2420-Apr-2410 days11-Aug-24
KS-16729-Mar-245-Apr-248-Apr-2411-Apr-246 days8-Aug-24
KS-16827-Mar-244-Apr-247-Apr-2412-Apr-249 days15-Aug-24
KS-1691-Apr-248-Apr-2411-Apr-2416-Apr-248 days9-Aug-24
Hierarchical Clustering Analysis and Trait Correlations

Hierarchical cluster analysis based on standardized trait values grouped the 19 selected genotypes into distinct clusters according to their phenotypic similarity (Fig. 3). The heatmap visualization revealed clear patterns of trait performance across genotypes, with color intensity indicating standardized z-score values ranging from -4 (dark blue) to +4 (dark red). One cluster was characterized by high kernel yield and kernel weight (red coloration), including genotypes KS-126, KS-058, and KS-001. Another cluster exhibited strong performance in split ratio (yellow coloration) while maintaining moderate values for other traits. Shell thickness showed considerable variation across genotypes, with KS-076, KS-096, and KS-143 displaying notably thin shells (dark blue coloration). Pearson correlation analysis revealed significant relationships among evaluated traits (Fig. 4). Size exhibited a moderate positive correlation with kernel yield (r = 0.48), indicating that larger nuts tended to have higher kernel percentages. Size showed a weak positive correlation with in-shell weight (r = 0.27) and weak negative correlation with split ratio (r = -0.22). In-shell weight demonstrated a positive correlation with kernel weight (r = 0.26) but showed negligible correlations with split ratio (r = 0.04), shell thickness (r = 0.00), and kernel yield (r = -0.02). Shell thickness exhibited a weak negative correlation with kernel weight (r = -0.20). Kernel yield showed near-zero correlations with split ratio (r = -0.03), shell thickness (r = -0.01), and kernel weight (r = 0.03).

Fig. 3

Hierarchical clustering heatmap of kernel-related traits in pistachio genotypes. Traits include size, kernel yield, in-shell weight, shell thickness, split ratio, and kernel weight. The color scale represents standardized values (Z-scores), with red indicating higher and blue indicating lower trait values. Groups were classified based on similarity among genotypes.

Fig. 4

Correlation matrix of nut and kernel traits in pistachio genotypes. The color and size of the circles represent the strength and direction of Pearson’s correlation coefficients among traits, including size, in-shell weight, split ratio, shell thickness, kernel yield, and kernel weight. Positive correlations are shown in blue and negative correlations in red.

Discussion
Genetic Diversity and Selection Efficiency

The comprehensive evaluation of 169 pistachio genotypes revealed substantial genetic diversity across all pomological traits, with CPI values ranging from -1.91 to 1.40 (Tab. 1). This remarkable variation demonstrates the rich genetic reservoir present within the germplasm collection and confirms the effectiveness of multi-trait selection approaches in pistachio improvement programs. The highly significant differences (p < 0.05) detected through analysis of variance, particularly for kernel yield (F-value = 292.82) and shell thickness (F-value = 80.72), indicate the strong genetic control of these commercially critical traits. Similar levels of genetic diversity have been documented in wild and cultivated pistachio populations across Central Asia, where studies of the wild range of P. vera have led most scholars to conclude that the crop originated in this region (though the exact location remains unclear), with the only extant wild populations occurring in Turkmenistan, Uzbekistan, Tajikistan, Kyrgyzstan, southeastern Kazakhstan, northern Afghanistan, and northeastern Iran (Whitehouse, 1957; Zohary, 1952; Khanazarov et al., 2009). The substantial phenotypic variation observed in this study parallels findings from molecular diversity studies, where Talebi et al. (2012) linked a high level of genetic diversity among P. vera to the different geographic distribution zones of analyzed accessions. Although wild pistachio drupes are smaller than domesticated ones (Zeng et al., 2019), it is impossible to confidently differentiate individual specimens of the morphotypes, due to the wide reaction norm of developmental plasticity in the species. Even nuts collected from the same tree may vary considerably in shape and size, as one recent morphological study showed that domesticated pistachio endocarps range from 9.91 to 19.15 mm in length and between 6.5411.65 mm in width (Rousou et al., 2021). This genetic richness provides a robust foundation for selective breeding, as maintaining diverse genetic resources is crucial for developing cultivars adapted to varying environmental conditions and evolving market demands. The identification of 19 elite genotypes representing the top 11.2% of evaluated germplasm demonstrates the efficiency of the Composite Performance Index in capturing superior trait combinations while maintaining genetic diversity for long-term breeding objectives.

Morphological Characteristics and Commercial Quality

The elite genotypes identified in this study exhibited exceptional combinations of commercially valuable traits that directly address market preferences and processing efficiency requirements. Genotype KS-168 emerged as the highest-ranking selection with exceptional kernel yield (58.80%, z-score = 2.65) and remarkably thin shell (0.53 mm, z-score = -1.65), representing a significant advancement in combining productivity with processing efficiency. Similarly, KS-103 demonstrated outstanding performance with kernel yield of 56.48% and impressive nut size (1.77 cm, z-score = 1.46), illustrating that simultaneous improvement in multiple quality parameters is genetically feasible. The identification of KS 088 possessing the thinnest shell among all genotypes (0.46 mm, z-score = -2.50) represents a critical genetic resource, as shell thickness directly impacts kernel recovery efficiency and commercial value. This trait, combined with appropriate rootstock selection, is essential for optimizing pistachio production. P. vera rootstocks express sensitivity to specific climatic and environmental conditions and have varying levels of tolerance to fungus, which influences cultivation success and, consequently, the expression of desirable kernel characteristics (Ferguson, 2008; Tekin et al., 2020). The remarkable uniformity in split ratio among selected genotypes, with 48 genotypes achieving perfect splitting (100%), demonstrates that natural dehiscence is a well-established trait in this germplasm pool. Given the immense phenotypic diversity among specimens, clear evidence for morphological change through domestication remains challenging for archaeobotanical documentation, therefore further research with more focus on the dehiscence trait is necessary (Mir-Makhamad et al., 2022). The presence of genotypes like KS-120, combining extraordinarily large nut size (2.12 cm, z-score = 2.81) with perfect split ratio, provides valuable breeding material for developing cultivars that satisfy premium market segments demanding both size and processing quality.

Trait Relationships and Breeding Implications

The correlation analysis revealed complex trait relationships that have profound implications for breeding strategy formulation and parent selection. The moderate positive correlation between size and kernel yield (r = 0.48) indicates that larger nuts tended to produce higher kernel percentages, suggesting that simultaneous improvement of these commercially important traits is achievable without antagonistic genetic effects. However, the near-zero correlations between kernel yield and split ratio (r = -0.03), shell thickness (r = -0.01), and kernel weight (r = 0.03) indicate that these characteristics are largely independent and require separate selection attention. The genetic independence of these traits is consistent with observations that pistachios can be sexually propagated, and their wild forms are not significantly different from their domesticated forms, suggesting that the landraces of pistachios observed around the world today were recently developed (Mir-Makhamad et al., 2022). As a dioecious and wind-pollinated plant, it is necessary to have enough male trees in close proximity to a female for sufficient pollination, and suitable ratios of male and female trees with appropriate distribution in an orchard are important factors taken into consideration by farmers today (Abu-Zahra and Al-Abbadi, 2007; Hormaza and Wünsch, 2007). The complex pollination requirements may have been one of the reasons for the delayed dispersal of the crop historically (Vargas et al., 2008). Despite its tolerance to numerous environmental factors, productive pistachio cultivation requires special environmental conditions, with experts suggesting that one of the most important variables in pistachio production is a hot and long summer (25–35°C) with low humidity (Everest, 2021; Tekin et al., 2020). The weak correlation between in-shell weight and kernel weight (r = 0.26) suggests that increasing nut size does not necessarily result in proportional increases in kernel mass, highlighting the importance of shell thickness management. These findings underscore the need for multi-trait selection indices that can simultaneously optimize multiple quality parameters while accounting for trait independence.

Phenological Stages and Harvest Dates

Our findings reveal significant variability in the phenological stages and harvest dates among the pistachio genotypes evaluated. KS-103 exhibited the earliest bud burst and full bloom, whereas KS-116 and KS-168 were the latest to reach harvest, aligning with findings from Burgos-Miranda et al. (2025), who reported that 'Golden Hills' pistachios bloomed approximately one week earlier than 'Kerman' and were harvested four weeks earlier, indicating genotype-dependent differences in flowering and maturation timelines. Furthermore, the flowering duration varied from 6 to 10 days among the genotypes, which is consistent with the study by Neyshaburi et al. (2021), where the flowering period of male pistachio genotypes ranged from 4 to 12 days. This variation can influence pollination efficiency and nut set, suggesting that selecting a mix of early and late-flowering genotypes could optimize yield and extend the harvest window. These findings highlight the importance of genotype selection for optimized flowering management, pollination strategies, and staggered harvesting in pistachio orchards, providing opportunities to extend the harvest period and improve overall yield management.

Selection Strategies and Practical Applications

The classification of elite genotypes into three distinct performance groups, high-performing, efficient shell type, and balanced type, provides a practical framework for matching germplasm to specific breeding objectives and production systems. The high-performing group, characterized by superior kernel yields ranging from 48.68% to 58.80%, is particularly suitable for intensive commercial orchards where maximum productivity is prioritized. The global production of pistachios has increased dramatically over the past few decades, from around 50 thousand tons in 1970, to 500 thousand tons produced globally in 2000, to more than 1 million tons in 2020, with Iran, Turkey, and the USA continuing to maintain the highest pistachio production (FAOSTAT, 2025). The inclusion of genotype KS-001, despite its relatively low overall ranking (106th), exemplifies the importance of considering specific trait excellence rather than relying solely on composite indices. This selection philosophy aligns with historical domestication patterns, as domestication for many long-generation perennials is characterized by a shift from sexual reproduction to vegetative propagation (Zohary and Spiegel-Roy, 1975). In the case of the pistachio, the acquisition of grafting techniques would have been essential for domestication, as the tree does not lend itself to simple vegetative propagation (Zohary et al., 2012). Today, grafting is used to improve resistance to diseases, pests, and other biotic and abiotic stress, and can also lend to the creation of unusual growth forms and reduce overall plant size (Mudge et al., 2009; Wang et al., 2017). The 19 elite genotypes identified in this study represent immediate resources for cultivar release following regional trial evaluation, with their greatest value potentially lying in their use as parental material for breeding programs targeting specific trait combinations. Future research should focus on multi-location trials to assess genotype × environment interactions and stability of performance across diverse growing conditions, as P. vera grows in both semi-arid deserts and dry slopes of low mountains and foothills due to its tolerance of harsh environmental conditions (Vavilov, 1967; Whitehouse, 1957).

Conclusions

Based on our results, the use of a Composite Performance Index proved highly effective in distinguishing superior pistachio genotypes with desirable commercial and breeding attributes. 19 elite genotypes were identified, representing highly efficient phenotypes with contrasting strengths, including high kernel yield, thin shells, and excellent splitting ability. KS-168, KS-132, and KS-103 emerged as outstanding performers with consistently high CPI scores, while KS-058 and KS-088 demonstrated unique efficiency traits such as compact nut size combined with high kernel-to-shell ratios and ultra-thin shells. The clustering of superior genotypes into three phenotypic classes, high-yielding, efficient-shell, and balanced types, highlights the potential for targeted parental selection according to breeding objectives. Weak correlations among certain traits, such as shell thickness and kernel yield, indicate the possibility of simultaneous improvement without strong tradeoffs. The significant genetic variation observed across pomological traits confirms the richness of the evaluated germplasm and its suitability for structured improvement programs. The identified genotypes, especially those combining high kernel yield with reduced shell thickness, should be prioritized in advanced breeding and regional adaptability trials. Integrating CPI-based ranking with phenotypic clustering provides a strategic roadmap for accelerating development of pistachio cultivars with improved processing efficiency and market value. Overall, the identified elite genotypes, especially those combining high kernel yield with thin shells, should be prioritized in future hybridization, regional adaptation trials, and cultivar registration initiatives to accelerate pistachio improvement.

DOI: https://doi.org/10.2478/mittklbg-2025-0016 | Journal eISSN: 3061-063X | Journal ISSN: 3061-0621
Language: English
Page range: 228 - 248
Submitted on: Oct 23, 2025
|
Published on: Dec 31, 2025
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2025 Mederbek Zhaanbaev, Ozkan Kaya, published by High School and Federal Office of Viticulture and Pomology
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.