Morphometric characterization
Agronomically interesting parameters, including fruit weight, thickness of pomegranate peel, and weight of 100 Arils are key traits for the fresh market and breeding programs (Arlotta et al., 2022). Data in Table 3 shows that phenotypical parameters varied, ranging from 92.27 g (A6) to 538.18 g (A9), 1.45 mm (A20) to 6.52 mm (A4), and 18.26 g (A14) to 55.48 g (A19), with an average of about 299.21 g, 3.47 mm, and 38.59 g for the fruit weight, thickness of the pomegranate peel, and weight of 100 Arils, respectively. Morphometric characterization of different accessions, particularly fruit weight, could depend widely on genetic and ecological condition (Melgarejo-Sánchez et al., 2015). Therefore, our results are comparable with those of other researchers, including Ejjilani et al., (2022) confirmed that fruit weight varied for 14 pomegranate genotypes from 168.0 to 558.0 g, with an average of about 339.5 g, and the thickness of the pomegranate peel ranged from 1.2 to 3.4 mm. Lahouel and Belhadj (2022) reported that the average fruit weight varies between a maximum value of 317.39 g for genotype (Messaad) and 135.85 g for genotype (Amourah) as a minimum value. Fruit weight ranged from 103.38 to 505.00 g (Akbarpour et al. 2009). Arlotta et al. (2022) reported that minimum and maximum differences among pomegranate genotypes ranged from 255.87 g to 825.82 g and were observed with an average fruit weight of 413.96 g. Adiletta et al. (2018) observed that Granato’, ‘Roce’ and ‘Wonderful’ showed a thin thickness of skin (3-3.9 mm), while ‘Dente di Cavallo’, ‘San Pietro’ and Mondrone Dolce’ were ranked among pomegranate genotypes possessing thick skin (5.1-6 mm). A previous researcher disputed that for 100 Arils (g) weighted among differing genotypes, the mean value for the weight of 100 arils was 34.7 g (Chater et al. 2018) and also Wetzstein et al. (2011) reported that the weight of 100 Arils for ‘Wonderful’ was reported to be 35.7 g, which is very close to and in agreement with our results. It is clearly noted that the pomological variations recorded on the 14 cultivars herein evaluated fall within the variation intervals for each single trait, with evidently a remarkable phenotypic effect with regards to the environment effect as well. Most of the fruit traits are greatly affected by the orchard management and environmental conditions, although in pomegranate these traits are also controlled by multiple genes (Harel-Beja et al. 2015).
Table 3
Shows a comparison of the genotypes for the morphological properties of the pomegranate accession.
Accessions
|
Fruit weight (g)
|
Thickness of pomegranate peel (mm)
|
weight of 100 Arils (g)
|
A1
|
339.14 c-f
|
4.18 de
|
42.98 cde
|
A2
|
352.74 b-e
|
4.45 cde
|
44.63 cd
|
A3
|
185.74 gh
|
2.19 jk
|
34.80 f-i
|
A4
|
382.05 bc
|
6.52 a
|
41.08 c-g
|
A5
|
367.97 bcd
|
3.46 e-i
|
41.91 c-f
|
A6
|
92.27 h
|
2.87 f-j
|
27.76 i
|
A7
|
450.82 ab
|
2.37 ijk
|
44.36 cd
|
A8
|
237.39 fg
|
4.12 def
|
28.51 i
|
A9
|
538.18 a
|
5.52 abc
|
47.90 bc
|
A10
|
368.32 bcd
|
2.83 g-j
|
39.72 d-h
|
A11
|
256.35 efg
|
2.54 h-k
|
33.88 ghi
|
A12
|
275.10 d-g
|
4.81 bcd
|
35.49 e-i
|
A13
|
119.58 h
|
3.63 d-i
|
30.09 i
|
A14
|
118.11 h
|
5.64 ab
|
18.26 j
|
A15
|
245.08 fg
|
3.84 d-g
|
33.54 ghi
|
A16
|
305.55 c-f
|
3.81 d-h
|
43.12 cde
|
A17
|
112.27 h
|
3.72 d-h
|
32.89 hi
|
A18
|
390.93 bc
|
2.74 g-j
|
53.78 ab
|
A19
|
338.73 c-f
|
2.88 f-j
|
55.48 a
|
A20
|
389.93 bc
|
1.45 k
|
40.59 c-h
|
A21
|
318.83 c-f
|
2.14 jk
|
40.86 c-g
|
A22
|
247.03 fg
|
2.58 g-k
|
39.04 d-h
|
A23
|
349.83 b-e
|
2.88 f-j
|
35.42 e-i
|
A24
|
399.00 bc
|
2.19 jk
|
40.10 c-h
|
Minimum
|
92.27
|
1.45
|
18.26
|
Maximum
|
538.18
|
6.52
|
55.48
|
Mean
|
299.21
|
3.47
|
38.59
|
Standard deviation
|
121.23
|
1.35
|
8.78
|
Influence of genotypes on biochemical traits
Results regarding biochemical characterization of the genotypes are summarized in Table 4, which shows that significant differences have been detected among 24 pomegranate accessions. The maximum TFC and TPC content were recorded in accession 3 with 60.07 µg/ml and accession 1 with 63.26 µg/ml, the minimum level was observed in accession 10 by 49.03 µg/ml and accession 19 with 38.91 µg/ml, respectively. The mean values of TFC and TPC were 54.50 and 47.97 µg/ml, respectively. In addition, antioxidant variability also was observed: accession 10 had the highest level of antioxidant, which was 31.16 µg/ml while accession 12 had lower content with 15.79 µg/ml and the mean value of 21.08 µg/ml was also observed. The total soluble solids content (TSS) showed significant variability with an interval of variation of 8.97–18.70 Brix with mean value of 14.59 Brix. The highest TSS content was recorded in accession 9, and the lowest levels were observed in accession 19. Our results are comparable to those found by Meena et al. (2021), who demonstrated that the flavonoid content of juice extracted from five different cultivars of pomegranates varied from 15.72 to 20.3 (CE mg 100 g-1) whereas the phenol content varied from 53.63 to 104.33 (GAE mg 100 g-1). In addition, total phenol content among genotypes was investigated that ‘Granato’ and ‘Roce’ showed a similar (1100 mg GAE L− 1) and significantly higher than ‘San Pietro’, ‘Mondrone Dolce’ and ‘Dente di Cavallo’ (750 mg L− 1); the highest concentration was found in cv. ‘Wonderful’ (Adiletta et al. 2018). Recent studies by El Moujahed et al. (2022) demonstrate that the total flavonoid content of juice extracted from different Moroccan local cultivars of pomegranates ranged from 62.47 to 140.44 mg GAE/100 ml, whereas the total phenol content for local varieties ranged from 885 mg to 1348 mg GAE/L. Nikdel et al. (2016) confirmed that ‘Rubab’ had the highest content of total flavonoids (4.93 mg GAE/100 ml) and 'Shahvar’ had the minimum of total flavonoids (2.24 mg GAE/100 ml), Rubab’ had the highest content of total phenolics (6.36 mg GAE/100 ml), and the lowest content of total phenolics was found in ‘Ghand’ as (1.78 mg GAE/100 ml). Reading antioxidant activity, the highest and the lowest antioxidant activity were detected in ‘Rubab’ (86.77%) and ‘Shahvar’ (79.54%), respectively. Additionally, levels of antioxidant activity in other research were 10.37–67.46% for seven cultivars of pomegranate juice in Turkey (Tezcan et al. 2009) and 18.6–42.8% for eight pomegranate juices from Iran (Mousavinejad et al. 2009).
Regarding total soluble solids (TSS), which estimate the level of dissolved sugars but also the presence of other soluble compounds such as acids, salts, water-soluble vitamins, and other chemical compounds, our results are similar to those found by (Adiletta et al. 2018). In the investigated pomegranate genotypes, the TSS range was from 17.0 to 18.5°Brix; in particular, ‘Granato’ and ‘Roce’ had the lowest total soluble solid content, corresponding to the lowest reducing sugars. Ejjilani et al. (2022) confirmed that the total soluble solids content (TSS) showed a significant variability ranging from 13.00 to 17.00 Brix. The TSS values ranged from 16.167 (Salakhani) to 17.467 ºBrix (Wonderful) in fresh juice (Abdulrahman et al. 2021). Hasnaoui et al. (2011) also reported similar TSS values in wild and cultivated pomegranate fruits, with respective ranges of 17.57–19.99 ºBrix and from 13.13 to 16.55 ºBrix. Another researcher confirmed that TSS was reported between 14.31 and 15. ºBrix (Melgarejo et al. 2011), 13.73–17.60 ºBrix (Mena et al. 2011) and 12.36–16.32 ºBrix (Martínez et al. 2012). 11.4–16.2 °Brix in Iran (Sarkhosh et al. 2009).
It is clear that the differences in chemical composition among the twenty-four accessions confirmed that varietal distinction impacts the biosynthesis of the proximate constitution and phytochemical components in pomegranate accessions. Therefore, these variabilities indicated the presence of an interesting variant among the pomegranate accessions for the ability to accumulate these compounds and suggested differences in genetic potential between them. Many researchers confirmed that phytochemical parameters may be influenced by many agents, including the growth year, growing region, cultivar, and maturity degree of the fruit (Abdulrahman et al. 2021). In addition, it is important to emphasize that several factors related to genotype, climate, and agricultural practices may explain differences in terms of juice (Ejjilani et al. 2022). Furthermore, researchers verified that the polyphenol and antioxidant potential levels of pomegranate juice could be significantly affected by the genotype (anthocyanin composition, and phenolic content) of some pomegranate cultivars with temperature during the maturity period and latitude of the growing region (Li et al. 2015). Moreover, the desirable trait in pomegranate juice is TSS, which is definitely important for fresh consumption, the food industry, and processing as it combines sweetness and flavor (Khadivi et al. 2020). These results showed that the levels of sugar and other physicochemical properties were different among various cultivars of pomegranate, which could be due to the existence of high genetic heterogeneity within the cultivars (Tehranifar et al. 2010). Eventually, morphological characteristics were used for descriptive purposes and are now commonly utilized to differentiate plant varieties. Furthermore, the ecosystem has a strong influence on morphological traits, and there are time and cost considerations. Therefore, morphological criteria alone are not sufficient to distinguish some different varieties that are morphologically similar. Therefore, the molecular fingerprinting of a plant variety is extremely important for protecting plant breeders’ rights (Pasquali et al. 2022).
Table 4
Phytochemical and antioxidant analysis of 24 pomegranate juices accessions.
Accessions
|
TFC (µg/ml)
|
TPC (µg/ml)
|
Antioxidant (µg/ml)
|
TSS (Brix)
|
1
|
58.01 abc
|
63.26 a
|
28.47 ab
|
16.90 cd
|
2
|
58.95 ab
|
60.64 ab
|
19.33 d-i
|
17.70 bc
|
3
|
60.07 a
|
58.20 bc
|
18.91 e-i
|
17.80 b
|
4
|
56.70 abc
|
59.51 abc
|
19.48 d-i
|
17.40 bcd
|
5
|
54.64 abc
|
43.78 fg
|
16.93 ghi
|
14.23 g
|
6
|
52.40 abc
|
56.14 cd
|
17.07 ghi
|
15.77 f
|
7
|
58.76 ab
|
45.09 fg
|
19.05 e-i
|
15.93 ef
|
8
|
53.52 abc
|
61.95 ab
|
20.96 d-g
|
15.60 f
|
9
|
57.64 abc
|
53.90 de
|
23.02 cde
|
18.70 a
|
10
|
49.03 c
|
40.79 gh
|
31.16 a
|
16.60 de
|
11
|
53.52 abc
|
41.72 gh
|
26.06 bc
|
16.77 d
|
12
|
52.02 abc
|
51.84 e
|
15.79 i
|
15.70 f
|
13
|
52.96 abc
|
43.03 fgh
|
20.96 d-g
|
15.43 f
|
14
|
49.59 bc
|
42.10 gh
|
21.95 c-f
|
17.30 bcd
|
15
|
52.21 abc
|
41.54 gh
|
19.05 e-i
|
9.53 l
|
16
|
52.58 abc
|
43.22 fgh
|
22.88 c-f
|
12.27 ij
|
17
|
52.96 abc
|
43.41 fgh
|
16.01 hi
|
13.40 h
|
18
|
53.71 abc
|
45.28 fg
|
22.10 c-f
|
12.17 j
|
19
|
55.21 abc
|
38.91 h
|
22.03 c-f
|
8.97 l
|
20
|
54.64 abc
|
42.28 fgh
|
20.54 d-h
|
11.23 k
|
21
|
56.33 abc
|
44.34 fg
|
19.97 d-i
|
13.00 hi
|
22
|
53.15 abc
|
41.72 gh
|
23.94 cd
|
13.17 h
|
23
|
55.77 abc
|
41.91 gh
|
18.20 f-i
|
13.80 gh
|
24
|
53.52 abc
|
46.78 f
|
22.17 c-f
|
10.83 k
|
Maximum
|
60.07
|
63.26
|
31.16
|
18.70
|
Minimum
|
49.03
|
38.91
|
15.79
|
8.97
|
Mean
|
54.50
|
47.97
|
21.08
|
14.59
|
Standard deviation
|
3.98
|
7.89
|
3.92
|
2.73
|
ISSR markers with polymorphism parameters
To analyze the genetic diversity of plants, a molecular marker study was used. In the current study, all ISSR primers gave scorable, good amplification products and presented high polymorphisms among the 24 pomegranate genotypes examined. The twelve random ISSR primers were utilized and generated 78 scorable polymorphic bands (Table 5 and Fig. 1). The number of amplified bands in our study ranged from 3 to 12 for ISSR markers (ISSR25, ISSR845) and ISSR847, respectively. In total 83 bands were scored among which 78 bands were polymorphic and five bands were monomorphic. The mean values, highest and lowest of polymorphic bands were (6.5,11 and 3) respectively. Our results showed that the polymorphism information content (PIC) values were recorded ranged from 0.58 to 0.90 with an average value of 0.79 per ISSR marker. The major allele frequency ranged from 0.21 to 0.58, with 0.31 as the average allele per marker. (58%) were observed for the ISSR25 marker, which has major alleles with the highest frequency. Values of the diversity gene were detected in the cultivated pomegranates at 0.61–0.91, with an average value of 0.81 per ISSR marker. In the genetic diversity research, the number of polymorphism bands with each marker technique employed may be extra-straightly linked to the ability to detect genetic differences among different accessions.
Therefore, our consequences are comparable to those of Talib et al. (2011), who found that the total number of polymorphic markers and percentage of polymorphism were 64% and 36.99%, respectively. The mean PIC value was 0.163, and the lowest and highest PIC values were 0.099 (ISSR5 and ISSR6) and 0.257 (ISSR11), respectively. Researchers demonstrated that nine primers were detected for the diversity and the genetic differences among different pomegranate varieties; total bands gave 88 bands, of which 72 bands had a polymorphism of 80.66% (Almiahy and Jum'a, 2017). Al-Mousa et al. (2019) evaluated the genetic variation among five pomegranate genotypes using twenty ISSR markers. Twelve ISSR primers were successfully used as fingerprinting tools and amplified 137 DNA fragments, of which 78 were polymorphic (56.93%). Primers 17, NLSSR3, and 16, which showed the highest polymorphism percentages of 90.41, 80.00, and 76.47%, respectively, with the highest number of unique bands of 6, 6, and 9, respectively. Genotype C amplified the highest number of DNA fragments (115) and unique bands (13). (PIC) ranged from 0.449 to 0.768, while genetic diversity ranged from 0.56 to 0.80. These results indicated that the ISSR technique was sufficiently informative and powerful to assess genetic variability in pomegranate. It is clear that a good technique to identify the pomegranate genotype is ISSR, these results are conducted with the results of Ajal et al. (2014). In addition, researchers confirmed the analysis of genetic diversity among the different cultivated pomegranate genotypes from Azerbaijan using the ISSR marker (Hajiyeva et al. 2018).
Informative molecular markers are commonly exploited for genetic diversity measurement improvement, genotypes that identify genotypes, populations for breeding, and conservation applications. As the genomic regions are highly polymorphic, different approaches to markers will produce varied consequences. The accuracy of the information should improve as the number of markers and genome coverage expand. The ISSR markers can be utilized to detect samples of plants that are closely linked in different populations. The ISSR markers with the most polymorphisms demonstrate how each locus can be used to estimate genetic diversity. Because primers with more alleles cover more of the genome, they are better at detecting variations in genetic make-up. It's possible that the variation in primer polymorphism is due to the fact that different primers have different nucleotide sequences. The results of this study show that ISSR markers are active in analyzing pomegranate genetic diversity and that the accessions studied have important genetic variation. An increased PIC value designates a higher degree of ISSR marker polymorphism in the genetic divergence analysis.
Table 5
Summary of PCR-ISSR amplified products including marker names, number of amplified bands, number of polymorphic bands, number of monomorphic bands, major allele frequencies, gene diversity, and PIC value that were obtained from 24 pomegranate accessions using twelve primers.
Marker
name
|
Number of amplifed bands
|
Number of polymorphic bands
|
Number of monomorphic bands
|
Major allele frequency
|
Gene diversity
|
PIC
|
ISSR6
|
10
|
10
|
0
|
0.38
|
0.72
|
0.67
|
ISSR9
|
9
|
9
|
0
|
0.33
|
0.83
|
0.81
|
ISSR10
|
7
|
7
|
0
|
0.29
|
0.85
|
0.84
|
ISSR11
|
8
|
7
|
1
|
0.25
|
0.86
|
0.85
|
ISSR12
|
4
|
4
|
0
|
0.29
|
0.80
|
0.77
|
ISSR25
|
3
|
3
|
0
|
0.58
|
0.61
|
0.58
|
ISSR811
|
7
|
6
|
1
|
0.25
|
0.85
|
0.83
|
ISSR812
|
6
|
5
|
1
|
0.29
|
0.83
|
0.81
|
ISSR841
|
10
|
9
|
1
|
0.21
|
0.90
|
0.90
|
ISSR845
|
3
|
3
|
0
|
0.29
|
0.81
|
0.79
|
ISSR847
|
12
|
11
|
1
|
0.21
|
0.91
|
0.90
|
ISSR849
|
4
|
4
|
0
|
0.33
|
0.78
|
0.74
|
Total
|
83
|
78
|
5
|
|
|
|
Mean
|
6.91
|
6.5
|
0.41
|
0.31
|
0.81
|
0.79
|
Clustering and population structure analysis of pomegranate genotypes.
Multivariate statistical methods are important in studying genetic diversity; one of them, cluster analysis, divides individuals into graphs based on intervals. The unweighted pair-group method (UPGMA) was used based on the Jaccard similarity coefficients to analyze the data, assessing and clustering for connections among pomegranate accessions (Rasul et al. 2022). The dissimilarity coefficients ranged from 0.23 (G4 vs. G5) to 0.63 (G13 vs. G14), and all 24 pomegranate accessions were grouped into 5 groups (A, B, C, D, and E) with a mean dissimilarity of 0.49 for twelve ISSR markers (Fig. 3). Cluster A includes only (G13 and G14), cluster B (G8 and G10), cluster C (G7, G9, G11, G19, G23, and G24), cluster D (G16, G6, G4, G5, G15, G2, G12, G3, G1, G18, and G22), and cluster E (G17, G20, and G21). Many researchers were recorded regarding pomegranate genetic diversity and their relationships, including Al-Mousa et al. (2019). They evaluated the genetic variation among five pomegranate genotypes using twenty ISSR markers and found that the genetic distance ranged from 0.26 to 0.37. Some genotypes showed wide divergence (C and A; C and D), while genotypes A and B were closely related. The unweighted pair-group method with arithmetic average (UPGMA) dendrogram grouped the genotypes into two clusters. In addition, Ajal et al. (2014) demonstrated the genetic relationship and diversity of Moroccan pomegranate germplasm among cultivars, and two main groups were designed, 1 and 2, based on the UPGMA dendrogram. In addition, studies indicate that the characterization of genetic markers for cultivar is also very promising by ISSR, which indicates an appreciably higher level of polymorphism in pomegranate plants. Almiahy and Jum'a (2017) reported that they evaluated the genetic diversity and the relationship among ten genotypes of pomegranates from different geographical regions of Iraq. According to the results of ISSR genetic distance using the UPGMA method, they showed that the 10 genotypes were distributed into two main groups, A and B.
Molecular variance (AMOVA)
The analysis of molecular variance (AMOVA) for the 24 pomegranate accessions in the ISSR investigation confirmed 87% of the total difference within the populations and 13% of the differences among populations (Table 6). Clustering and structure analyses were confirmed based on AMOVA results.
Table 6
Illustrated the molecular variance (AMOVA) of the five populations for 24 pomegranate accessions.
Source
|
df
|
SS
|
MS
|
Est. Var.
|
%
|
P-Value
|
Among Population
|
3
|
89.125
|
29.708
|
2.369
|
13%**
|
0.001
|
Within Population
|
20
|
315.125
|
15.756
|
15.756
|
87%**
|
0.001
|
Total
|
23
|
404.250
|
|
18.125
|
100%
|
|
Structure analysis among 24 accessions by ISSR Markers
Allele frequencies using STRUCTURE analysis was utilized to determine for 24 pomegranate accessions by Evanno et al. (2005). Results showed that four groups or sub-populations were represented by colors, including group 1 (red line), group 2 (green line), group 3 (blue line), and group 4 (yellow line) according to Delta K genotypes (Fig. 4B). In addition, many genotype combinations were observed with more than one background, including (2, 5, 6, 7, 8, 10, 11, 12, 13, 15, 16, 17, 18, 23, and 24), and the rest shown is another background that may have had a complex history from the gene flow among taxa that caused connecting intercrossing or practicable consequences. In addition, altering weather conditions within the locations may be contributing to the high variability between accessions. The right number of clusters (K) was determined and detected in a sample of individuals that the peak was started at 3 and ISSR markers number 12 showed the real K value with the highest value of K = 4 (Fig. 4A). Unlike the dendrogram plot, our analysis shows that the genotypes can be divided into four distinct populations. The accessions distribution patterns in dendrogram plots are also the most similar to the distribution patterns in STRUCTURE analysis. The clusters' accessions contain the majority of the diversity that may be imply the presence of admixed among accessions.