Genotyping of Croatian Olive Germplasm with Consensus SSR Markers

: Leaf samples of 226 cultivated olive trees were collected from traditionally managed olive orchards and genotyped with eleven consensual SSR markers. The proportion of shared allele distance was used for the estimation of distances between olive genotypes. Cluster analyses were performed using a Fitch–Margoliash least-squares algorithm. The number of different genetic sub-groups of olive genotypes ( K ) was investigated using STRUCTURE analysis. The standardization of allele lengths was performed to enable the comparison SSR profiles of Croatian olive genotypes with olive profiles obtained with the same SSR primers in OleaDB and WOGB databases. Overall, 73 SSR profiles of known Croatian varieties and 53 profiles of unknown olive genotypes were differentiated. Synonyms were detected in 18 varieties, and we found intra-varietal differences in 15 varieties. Three genetic subgroups of olive genotypes were determined. Following allele length standardization achieved using nine referral samples, the genetic profiles of 126 cultivated olive genotypes were compared to OleaDB and WOGB databases, out of which 92 genotypes were found to be unique to Croatian olive germplasm. The results revealed the wide genetic diversity of olive germplasm beyond the known, registered varieties. The FAZ_oliveDB database containing the profiles of 126 Croatian olive genotypes was created and made available for public use.


Introduction
The olive tree crop (Olea europaea L.) is widely cultivated, generally being grown in the area of the Eastern region of Mediterranean basin along the European and African coasts [1].Olive was domesticated during the Chalcholitic Levant, about 6800-6300 years BP, and spread over the Mediterranean basin [2].However, our understanding of the heartland of primary olive domestication must be enlarged to the Levant and not only focus on the Jordan Valley [3].Shortly after these conclusions, Besnard et al. [4] pointed out that the northern Levant is the primary area of olive domestication.Since the period of olive domestication began, large numbers of cultivated olive varieties have spread until now thanks to the vegetative propagation of plants with the best trait combinations.This previously occurred as a result of random crosses between cultivated (ssp.sativa) and wild (ssp.oleaster) olive plants (i.e., feral forms of olive [5]) or as natural mutations.The "Olive Germplasm: Cultivars and World-Wide Collections" includes data on 1250 cultivars in 54 countries, conserved in more than hundred collections [6].
The olive is a species of great economic importance, not only in the countries of Mediterranean basin.Total world production in 2021 was 3,010,000 tons of olive oil [7].European production led the way, with the combined production of Spain, Italy, Greece, and Portugal standing at around 1,974,000 tones, and non-EU but IOC members such as Algeria, Argentina, Morocco, Egypt, Tunisia, and Turkey were estimated to have collectively produced 936,000 tons of olive oil.The olive tree has been cultivated in Croatia for over 2000 years [8].The early division of olive-growing regions in Croatia into southern, middle, and northern Dalmatia, the Kvarner islands, Istria, and inner Dalmatia [9] was largely supported by local microclimate regions.The production of olives in Croatia is proportional to the size of the Mediterranean part of the country, which is smaller comparing to the leading producing countries but still very significant in terms of germplasm diversity.Olive oil production was 32,026 hl [10] in 2021.In 2020, olives were grown on an area of 20,282 hectares, with an average yield of 1.6 t/ha and a total production of 33,230 tons of fruit.However, in the last 15 years, Croatian olive oil producers have won numerous silver and gold awards in international competitions for the organoleptic qualities of their olive oil.
Numerous studies of the local olive germplasm in Croatian territory were based on morphology [13][14][15][16][17][18][19][20].However, their morphological characterization was not strong enough for the precise differentiation of cultivars, especially of synonyms and homonyms, due to limitations caused by environmental influences.A great leap forward in the characterization of varieties was achieved using PCR-based DNA markers, which proved to be applicable for the investigation of genetic diversity, the identification of olive germplasm, and the investigation of intra-varietal variability.Among the others, the use of SSRs is one of the molecular techniques that may provide useful information on the level of polymorphism and diversity in the characterization of olive germplasm accessions [21] and may be applied for the discrimination of closely related cultivars and association mapping [22].SSR markers are far and away the ideal markers, being abundant and uniformly distributed over the genome, codominant, highly polymorphic and informative, highly reproducible, easily produced by PCR, and relatively simple to interpret.They produce data that are easily exchangeable among laboratories [23,24].
Olive genotypes from Croatia, either autochthonous or introduced, were genotyped only partly with different sets of SSR markers, comprising varieties from south Dalmatia [25] and Istria [26,27].Several studies focused only on the molecular characterization of distinct varieties [28,29], the relationships among olive varieties native to Croatia and Turkey [30], the clarification of possible cases of synonymy and homonymy [11], and the determination of intra-varietal diversity [31].In some reports, Croatian varieties were represented only by one or several of the most common varieties as part of other studies comprising genotypes directly sampled from the wider Mediterranean area [32][33][34][35], or by genotypes originating from Croatia that were sampled at olive collections in other countries [36][37][38].However, no comprehensive and systematic work on Croatian germplasm has been performed yet.Many other studies from the Mediterranean basin have been conducted recently in order to determine genetic variability and perform genotype identification.These studies cover different areas, from those locally focused on just one country, like Albania [39,40], Montenegro [35], Italy [41][42][43], Turkey [44][45][46][47], Greece [48], Iran [49], Lebanon [50], Tunisia [51][52][53], or Algeria [54], or even narrow producing areas like the West Bank [55], southern Anatolia [56], Sicily [57,58], southern Italy [59,60], parts of Spain [21,61], or Madeira [62], to those comprising several countries, but with fewer samples per country [34,[63][64][65][66]. Gomes et al. [67] analyzed the main Portuguese varieties, comparing them with several Italian and Spanish varieties using SSR and ISSR markers.This was carried out after a comparison had already been performed between the main Italian and Spanish varieties using SSR markers [68].With the objective of the standardization and the con-sensual application of SSR markers, and due to the need for more convenient methods of result exchange and data comparison, markers and protocols were proposed for the improvement of the application of SSRs in order to develop a robust method with which to track the origin of olive cultivars [69,70].These were overviewed [71] alongside the other molecular-marker techniques applied in olive germplasm studies.
The SSR markers are also widely used for the assessment of intra-varietal differences.There were studies that aimed to assess the on-farm variability of only two local cultivars of olives in the Alentejo region [72], to compare Iranian ecotypes and varieties with cultivars from the Mediterranean basin [73], and to compare olive clones from Morocco with varieties from other countries [74].The other studies analyzed the intra-varietal differences between three major Italian cultivars [75], the range of variability within the three major Iranian cultivars [76], and the intra-varietal differences and genetic relationships of some of the most widespread homonymic olive cultivars in the East Adriatic region [77].
Our hypothesis assumes the presence of huge genetic variability within the existing olive germplasm in Croatia.The objectives of the study are as follows: 1.
To genotype a large and representative number of olive samples using SSR markers, including olive traditional varieties and unnamed genotypes of economic importance; 2.
To identify synonyms and homonyms of Croatian varieties; 3.
To create a database of olive SSR profiles comparable to publicly available olive SSR databases and to detect genotypes which are unique only to the Croatian part of the Mediterranean.

Plant Material and DNA Isolation
In collaboration with the local community and olive producers, we localized, labelled, and sampled 226 cultivated trees of old varieties with uncertain status of synonymy, trees of the same variety with obvious in situ variability, and trees of unknown identity.All the trees were localized in old, extensive or traditionally managed olive orchards, covering five olive-growing zones in Croatia: 26 trees in Istria, 27 trees in Kvarner, 86 trees in northern Dalmatia, 21 trees in middle Dalmatia, and in 66 trees in southern Dalmatia.The age of the labelled trees ranged from at least 50 years to even more than 1000 years.The trees were selected and attributed according to the best practice and knowledge, respecting the descriptions published by previous authors [6,[13][14][15][16][17][18][19][20].Additionally, nine foreign varieties ('Ascolana tenera', 'Frantoio', 'Gentile di Chieti', 'Grossa di Spagna', 'Leccino', 'Moraiolo', 'Nostrana di Brisighella', 'Pendolino', and 'Picholine'), obtained from the nursery in Italy [Vivai Pietro Pacini, Via L. Galeotti, 1-51017 Pescia (PT)], were used as referral samples.Six of them ('Ascolana Tennera', 'Frantoio', 'Leccino', 'Pendolino', and 'Picholine') were introduced varieties that have been traditionally cultivated in Croatia for a long time.These referral varieties were chosen as the most appropriate options for covering the whole allele range of the primers applied in this research and enabling allele standardization among the databases.
Samples of fresh leaf tissues were collected in situ from the canopy and immediately put in zip bags with silica gel (in a ratio of 1 g of a fresh leaf sample to 5 g of silica gel) and dried for three weeks in the fridge at +4 • C. Dried leaves were ground into fine powder for 30 s at the vibrational frequency of 30 Hz using a mixer mill MM400 (Retsch GmbH, Haan, Germany).Genomic DNA was isolated using a DNA plant isolation kit (DNeasy Plant Mini Kit; Qiagene, Venlo, The Netherlands) following the protocol provided by manufacturer.The quality of the isolated DNA was checked in 0.8% agarose gel in a 0.5 × TBE buffer (45 mM Tris, 45 mM Boric acid, and 1 mM EDTA).The quantity of isolated DNA was measured in VersaFluorTM Fluorimeter (ver.170-2402EDU; Bio-Rad Laboratories, Hercules, CA, USA), and diluted with millipore water to a working concentration of 5 ng µL −1 for further molecular analysis.

Molecular Analysis
Among the markers previously designed by different authors [1,41,78,79], Baldoni et al. [70] suggested a consensus list of 11 the most effective markers for characterizing olive accessions and proved their use for comparison of the results obtained among different laboratories regarding different PCR conditions and techniques, cycling conditions, separation and visualization methods.Forward primers were labelled with fluorescent dyes (6FAM, VIC, and NED; Applied Biosystems, Foster City, CA, USA) (Table 1).The PCR reaction was performed in a VeritiTM 96 Well Thermal Cycler (Applied Biosystems, Foster City, CA, USA) in a reaction mix with a total volume of 10 µL (Table S1) at 72 • C; for 28 cycles of 30 s at 95 • C, 30 s at 55 • C, and 30 s at 72 • C; and a final extension step of 5 min at 72 • C. Amplified fragments were separated in the four capillary electrophoresis (Genetic Analyzer 3130; Applied Biosystems, Foster City, CA, USA) using GeneScan™ 500 LIZ™ dye of standard size (Applied Biosystems, Foster City, CA, USA).The PCR products of the 11 SSR primers of each olive DNA sample were mixed together two or three at a time, according to the dye and expected allele size, and put through capillary electrophoresis in the following combinations: standardization DCA3/DCA18/UDO43; DCA9/DCA14/DCA16; DCA5/GAPU71B/GAPU103A; and EMO90/GAPU101.The allele size of the amplified fragments was scored using the GeneMapper 4.0 software (Applied Biosystems, Foster City, CA, USA).

Allele Lenght Standardization
The standardization of allele lengths of Croatian olive genotypes based on the allele lengths of the 80 most important Mediterranean olive varieties from 11 Mediterranean countries, with data available in OleaDB [6], was accomplished using referral varieties of specific alleles for each locus.The standardization of allele lengths was performed to enable the comparison of the SSR profiles of Croatian olive genotypes obtained from the combinations of the eleven consensual SSR markers, with olive profiles obtained using the same SSR markers on the Olea database (Table S2).The standardized data of SSR allele lengths of Croatian olive genotypes were additionally standardized with the allele lengths of olive accessions of the Worldwide Olive Germplasm Banks of Córdoba (WOGBC) [37] and Worldwide Olive Germplasm Banks of Marrakech databases (WOGBM) [36].To standardize the length of alleles, we used the genotypes present in our research as referral samples, and those which are also present in the databases OleaDB, WOGBC, and WOGBM.

Data Analysis
Olive trees with identical genotypes were identified using the GenClone programme [80] and excluded from further analysis, with the exception of the analysis of variability within varieties.Information index (I) scores, the number of alleles (Na), the number of effective alleles (Ne), the number of rare alleles per locus (Nu), and the number of distinguished genotypes (NDG) were calculated for each SSR locus.The expected (He) and observed (Ho) heterozygosity for each SSR locus were estimated according to the work of Nei [81], while the polymorphic information content (PIC) was estimated according to the research of Botstain et al. [82].I, Na, Ne, Nu, NDG, He, Ho and PIC were calculated using GenAlEx 6.501 software [83].The probability of identity [84] for each SSR locus and the overall probability were estimated using software IDENTITY 1.0 [85].The power of discrimination-PD [86] was calculated for each SSR locus, where the allele frequency was replaced by the genotype frequency [1].
The proportion of shared allele distance (Dps, [87]) was used for the estimation of distances between olive genotypes using MICROSAT v.1.5software [88].Allele differences between pairs of Croatian olive genotypes were estimated using GenAlex 6.501 software [83].Cluster analyses were performed on SSR loci using a least-squares algorithm [89] with 1000 bootstraps using the FITCH programme of the Phylip software, ver.3.695 [90].The relationships among olive genotypes were visualized using Display Newick Trees constructed in MEGA7 [91].
For the analysis of intra-varietal variability, a matrix of allele differences among trees was used for the construction of a dendrogram using the UPGMA (unweighted pair group method using an arithmetic average) algorithm.Matrices of allele differences were calculated using GeneAlex program, ver.6.5 [83].UPGMA analysis was conducted using NTSYSpc software, ver.2.21 L [92].
A number of different genetic subgroups of olive genotype K (the modal value of ∆K) were investigated using STRUCTURE, ver.2.3.4 software [93].STRUCTURE analyses included a burn-in period (the initial stage of the sampling process) of 200,000 replicates, followed by 600,000 Markov chain Monte Carlo (MCMC) replicates for each run.Twenty repeat runs were carried out to quantify the amount of variation in the likelihood for each K (from K = 1 to K = 8) using an ADMIXTURE model and correlated allele frequencies.The posterior probability of the data lnP(K) for a given K can be used as an indication of the most likely number of real genetic groups [94].Therefore, the height of the modal value of the ∆K distribution was calculated in order to detect the number of real genetic subgroups K using Structure Harvester, v 0.6.92[95].The K that best described the data was chosen by examining the lnP(K) [96] and by calculating ∆K [94].The value of K with the highest mean log likelihood [lnP(K)] and ∆K statistic was selected.An olive genotype having a membership probability of Q > 75% indicated its belonging to a specific genetic group (K), while genotypes with membership probabilities of Q < 75% belonged to admixed groups of genotypes of mixed origin [97].

Results
We analyzed 226 olive trees grown in Croatia and assessed additional foreign referral samples of 9 trees (in total 235 trees) at eleven SSR loci.Duplicates (91 trees) were excluded from further statistical analyses, reducing the number of trees to 135 trees of different genotypes (Table S3).The analyses of these 135 genotypes resulted in 138 polymorphic and reproducible alleles, out of which 30% (42) had a frequency of less than 1%.Allele frequencies varied from 0.004 to a maximum of 0.567 for the allele 191 bp at the locus DCA14.The mean number of alleles per SSR locus was 12.5 (Table 2).The number of alleles obtained ranged from 5 (EMO90) to 19 (UDO43), while allele size ranged from 117 bp (GAPU71B) to 259 bp (DCA3).The number of effective alleles (Ne) per locus ranged from 1.36 (EMO90) to 7.20 (DCA9).The expected heterozygosity (He) ranged from 0.611 (DCA5) to 0.861 (DCA9), with a mean value of 0.766.The observed heterozygosity (Ho) was lowest at locus DCA14 (0.622), and highest at locus UDO43 (0.985).Loci DCA9 and DCA16 showed lower value Ho in relation to He.The mean value of Ho was higher (0.823) than the mean value of He (0.766), indicating high genetic variability within analyzed the germplasm.Most of SSR loci (9) have higher Ho then He, which is also reflected in the negative value of fixation index (F) (Table 2).The PIC value of the 11 SSR loci was above 0.5 and ranged from 0.56 (DCA5) to 0.86 (DCA9), with a mean value of 0.74.The average number of distinguished genotypes per SSR locus (NDG) was 27.The SSR locus DCA9 showed the highest number of identified genotypes (NDG = 42) while the loci EMO90 showed the lowest (NDG = 11).The probability of identity (PI) of genotypes varied from 0.019 (DCA14) to 0.210 (DCA5), while the cumulative PI value for all SSR loci was 5.86 × 10 −15 .The obtained PIC, NDG and PI values indicated that the SSR loci had very good discriminating power, successfully differentiating 73 SSR profiles of known Croatian varieties (including the cases of intra-varietal variability), 9 profiles of foreign varieties, and 53 profiles of unknown olive genotypes, from a total 135 SSR profiles (Table S3).The unknown (UNK) olive genotypes are the genotypes not documented or known by name.They also belong to the Olea europaea subsp.europaea species (not the Oleaster) and they are normally cultivated, but their genetic profiles are not identical to the profiles of any native Croatian or introduced variety.Within Croatian germplasm, 34 alleles were detected as rare or private (Table 3).These alleles were confirmed via repeated amplification.There were 34 private alleles within the set of 138 polymorphic alleles detected within Croatian germplasm, and 25 out of 135 genotypes possessed one to five private alleles.However, comparing SSR profiles of Croatian olive germplasm with the profiles of olive varieties in the Olea database and olive accessions in the WOGBC and WOGBM databases, we found that the number of rare alleles decreased from the initial 34 to 12 rare alleles, only represented within and unique to germplasm originating from the Croatian part of Mediterranean (Table 3).

Synonyms and Homonyms
Synonyms were detected in 18 cultivars, mostly in 'Oblica' (9 synonyms), 'Drobnica' (7 synonyms), and 'Plominka' (7 synonyms).Few synonyms were detected in the remaining 15 varieties of local significance (Table 4), the production of which is limited to small geographic areas.Homonyms were not detected in this research, but we did detect numerous cases of intra-varietal variability (Table 5).Intra-varietal differences were detected for a variable number of loci in 15 varieties and for one unknown genotype (Table 5).The number of variable loci ranged from one locus in seven varieties and one unknown genotype to two loci in 'Plominka', 'Oblica', 'Bjelica dubrovačka', and 'Aserija', three loci in 'Piculja' and 'Karbunčela', four loci in 'Krvavica mljetska', and five loci in 'Drobnica'.The loci DCA9 (27 allelic combinations), UDO43 (18 allelic combinations), and Gapu103A (8 allelic combinations) proved to be the most sensitive in terms of the detection of differences within varieties.Some of the alleles were differentiated in only two base pairs per locus.These small variations were proved by repeated PCR amplification.UPGMA clustering (Figure 1) based on allelic differences showed variation in olive trees within 15 varieties and one unknown genotype.

Genetic Distances and Relationships of Croatian Olive Germplasm
The proportion of shared allele distance (Dps) between 135 olive genotypes ranged from 0.045 (27 pairs of genotypes) to 0.955 (12 pairs of genotypes) with a mean value of 0.629.
Olive genotypes were divided into five clusters (C1-C5) usinge Fitch-Margoliash least-squares algorithm cluster analysis (Figure 2a).Varieties and unknown genotypes with intra-varietal variability are only represented by one representative in the cluster.The genotypes from the same region clustered mostly into the same cluster, but some of them were grouped into other clusters, mostly with the genotypes from another regions.The genotypes from Istria and Kvarner (46% of genotypes) clustered mostly within the cluster C1.Genotypes from the northern Dalmatia clustered mostly within the cluster C2 (50% of genotypes), and those from southern Dalmatia clustered within cluster C3 (45% of genotypes).The genotypes from middle Dalmatia mostly clustered into the clusters C3 (27%) and C4 (27%).The traditional varieties 'Drobnica', 'Oblica', and 'Plominka', cultivated in a wide area from Kvarner to southern Dalmatia, were grouped in separate clusters.'Drobnica' a variety mostly presented in Kvarner, middle and southern Dalmatia, was grouped within the cluster C2.Recently, the trees of 'Drobnica' were found even in Istria, in the north.'Oblica', the leading autochthonous variety in Dalmatia, and partially produced also in Kvarner, was grouped within cluster C5.The variety 'Plominka', mainly cultivated in Kvarner region, Istria, and southern Dalmatia, was also grouped within cluster C5.
Bayesian clustering analysis was used to test the genetic structures of Croatian olive germplasm, and one to eight clusters (K) were tested.The average log probability lnP(K) values increased up to the third testing cluster, after which the rate of change in the log probability value decreased (Figure 3).The relationship between K and ∆K indicated the clustering of olive genotypes into three distinctive genetic groups (Figure 2b).Genotypes of mixed origin (with membership probabilities of Q < 75%) are marked by black frames in Figure 2b.Genetic subgroups K1 (marked in red), K2 (blue), and K3 (green) consisted of 32, 59, and 44 genotypes, respectively.The varieties were mostly represented within the genetic subgroup K3 (43%), followed by genetic subgroup K2 (35%), and the least represented within genetic subgroup K1 (22%) (Figure 2b).The unknown genotypes were mostly distributed within the subgroup K2 (60%), while a slightly smaller part of unknown genotypes were distributed within genetic subgroups K1 (21%, and K3 (19%).The majority of olive genotypes originating from Istria and Kvarner were distributed within genetic subgroup K1, the majority of genotypes from the north and some genotypes from the middle Dalmatia were within the subgroup K3, while the majority of genotypes from the middle and all genotypes from southern Dalmatia were distributed within the subgroup K2.
The mean value of membership coefficients of Croatian olive genotypes, sampled in five olive-growing areas in Croatia, was calculated by each growing area on the basis of membership coefficients for K = 3 for olive genotypes (pie charts, Figure 4).The genetic basis of the plants originated from the three genetic subgroups.The most represented genetic subgroup in southern Dalmatia was K2 (blue) (52%).Moving northward, the share of that genetic subgroup dropped drastically to only 6% in Istria.The share of genetic subgroup K1 (red) was lower in southern Dalmatia (15%) than Istria (45%).The share of genetic subgroup K3 (green) was stable compared to remaining genetic subgroups, with the share between 33% in southern Dalmatia and 49% in Istria and middle Dalmatia (Figure 4).The mean value of membership coefficients of Croatian olive genotypes, sampled in five olive-growing areas in Croatia, was calculated by each growing area on the basis of membership coefficients for K = 3 for olive genotypes (pie charts, Figure 4).The genetic basis of the plants originated from the three genetic subgroups.The most represented genetic subgroup in southern Dalmatia was K2 (blue) (52%).Moving northward, the share of that genetic subgroup dropped drastically to only 6% in Istria.The share of genetic subgroup K1 (red) was lower in southern Dalmatia (15%) than Istria (45%).The share of genetic subgroup K3 (green) was stable compared to remaining genetic subgroups, with the share between 33% in southern Dalmatia and 49% in Istria and middle Dalmatia (Figure 4).

Discussion
In addition to the reports on global olive germplasm based on collections in situ [37,38], there are many published papers with a focus on olive germplasm ex situ and its diversity within certain geographic areas [46,48,50,53,54,58,60,72,98,99].Regarding Croatian olive germplasm, different sets of SSR loci were used to genotype traditional olive varieties cultivated in southern Dalmatia [25] and in Istria [27].Our report was based on the first comprehensive research on the Croatian olive germplasm using 11 SSR markers, which were suggested and recommended as consensual markers for the genetic characterization of olive germplasm [70].The variation in the allele size in the comparison of SSR profiles of domestic germplasm with other established international databases of SSR pro-

Discussion
In addition to the reports on global olive germplasm based on collections in situ [37,38], there are many published papers with a focus on olive germplasm ex situ and its diversity within certain geographic areas [46,48,50,53,54,58,60,72,98,99].Regarding Croatian olive germplasm, different sets of SSR loci were used to genotype traditional olive varieties cultivated in southern Dalmatia [25] and in Istria [27].Our report was based on the first comprehensive research on the Croatian olive germplasm using 11 SSR markers, which were suggested and recommended as consensual markers for the genetic characterization of olive germplasm [70].The variation in the allele size in the comparison of SSR profiles of domestic germplasm with other established international databases of SSR profiles was standardized, using referral genotypes that possess specific alleles for each locus as a standardization bridge.
The focus of this research was the determination of the genetic diversity of olives in Croatian part of the Mediterranean, with a special emphasis on the DNA identification of lesser-known varieties and unknown genotypes found in old and extensive olive orchards.Lesser-known varieties and unknown genotypes are mainly of local significance, but they are still important, not only as a pool of genetic variability, but also as an economic niche for the production of virgin olive oils, and as a basis for branding the final product with a mark of local origin authentic to the region (island, town) where it is produced.The results of this research represent an important contribution in terms of establishing a database of SSR profiles of domestic olive germplasm comparable with the databases of olive germplasm OleaDB [6], WOGBC [37], and WOGBM [36,38].
The results of SSR analyses revealed high genetic diversity within the olive germplasm in Croatia.Altogether, 126 genotypes originating from Croatian part of the Mediterranean were identified using 11 SSR loci.The major of SSR loci (nine out of 11) showed higher values of observed heterozygosity (Ho) compared to expected heterozygosity (He), as well as negative values of the inbreeding coefficient (F).These parameters indicate the existence of high genetic variability within the analyzed germplasm.The existence of high genetic variability within the germplasm is explained by the emergence and accumulation of spot mutations during the long period of the clonal reproduction of olive varieties [37,54,70,100].
The mean value of PIC (>0.7) for the set of consensual SSR loci was high, and as such was suitable for the assessment of molecular diversity and the identification of olive genotypes [27,64].These SSR loci proved to be enough powerful for the separation of genotypes within Croatian olive germplasm.The average number of distinguished genotypes per SSR locus was multiple times higher (NDG avg = 27) in comparison to the SSR loci applied in previous reports [54,58,59].The cumulative probability of the identity value (PI) of these loci was very low, and this showed that SSR loci used in this study had high discriminating power in genotyping of Croatian olive genotypes, which was in accordance with previous reports [36,60,101].
The grouping of the olive genotypes based on the results of the SSR analyses was not related to their geographic distribution along the Adriatic coast.Despite being grouped into three genetic subgroups according to their genetic origin, the results confirmed the high discriminatory ability of the consensus SSR markers.The high diversity of the genotypes in our research was a result of their different genetic origin and probably a lack of selection pressure.
Three genetic subgroups (K = 3) were determined via Bayesian clustering analysis, which was used to test the genetic structure of Croatian olive genotypes.All olive genotypes were distributed within these three genetic subgroups.Three fundamental genetic subgroups were also determined in previous studies of the olive genotypes from the whole Mediterranean [33,36,70].These three genetic subgroups are distributed in the eastern, middle, and western Mediterranean.The presence of the three genetic subgroups was determined in middle Mediterranean [36] and the values obtained corresponded to our results for the Croatian part of the middle Mediterranean.The main reason for that was probably the possibility of the spreading of generative and vegetative propagation material within the entire Mediterranean basin [36].
In our research, consensual SSR primers proved to be a tool of choice for distinguishing variable genotypes within a variety.In that sense, the loci DCA9, UDO43, and Gapu103A were the most effective, which corresponded with the identical conclusion [70] reached on the usefulness of this subset of markers in terms of segregating variable genotypes within the same variety.Considering the variability within several olive varieties in eastern Adriatic, the loci DCA9 and Gapu103A were the most effective [77].In our research, the intra-varietal variability was noticed in the very old and important varieties cultivated today, like 'Oblica', 'Drobnica', 'Karbunčela', and 'Piculja'.The variety 'Oblica' is the most represented olive variety along the whole Adriatic, counting for 60% of the total olive trees in Croatia [8].
The variety 'Drobnica' is mostly represented in northern Dalmatia [102], but is also present in southern Dalmatia [18].The name 'Drobnica' is mostly used in the northern part of Dalmatia, while the name 'Sitnica' is mostly used in southern part of Dalmatia [9].This research shows that the variety 'Drobnica' is also represented in Istria, that 'Karbunčela' is mostly represented in northern Dalmatia, in the islands and in the area around Zadar, while 'Piculja' is mostly represented in southern Dalmatia [13].The highest intra-varietal variability was detected for the variety 'Drobnica'.The clones of 'Drobnica' differed in as much as six out of 22 alleles.This kind of variability might be explained as polyclonal variations that arise as a consequence of the vegetative propagation of different genotypes [75,77].
Within 'Piculja' and 'Karbunčela', both varieties of local importance, variations were found in up to three out of 22 alleles.This variation might be explained as a result of the accumulation of somatic mutations [36].Variation was found within 'Oblica', the most represented variety in Croatia [8], in up to 1 out of the total of 22 alleles.'Oblica' is probably the variety that does not tend to accumulate somatic mutations.The variation in up to two alleles among clones is due to somatic mutation occurring in different branches of the ancient mother plants used for vegetative propagation [75].
A list of varieties, their names, and synonyms is provided in one of the first documented sources on olive varieties [13].Working based on descriptions and data collected in situ, Bulić [13] states that 'Dužica' and 'Duška' are synonyms.Also, synonymy was attributed to following three pairs of varieties: 'Mrčakinja' and 'Piculja', 'Paštrica' and 'Žutica', 'Puljizica' and 'Puljka' [13].In our research, we found unique genetic profiles for all of these putative synonyms, and hence established their varietal uniqueness, in spite of similar names and previous statements.Contrary to the attribution to varieties 'Murgulja' and 'Krvavica' of the characteristic of being synonyms [13,16], we found unique genetic profile for all of these varieties.However, there are two different varieties, known as 'Krvavica': 'Krvavica mljetska', found in southern Damlatia (island if Mljet), and 'Krvavica skradinska', found in middle Dalmatia (region of Skradin).With the addition of 'Murgulja', which was considered as a synonym [13], these genotypes are three different varieties, each with unique genetic profiles.The variety 'Velika lastovka' is cultivated in southern Dalmatia, but we discovered this variety far away in the Kvarner region and northern Dalmatia.
Within the scope of this research, for the first time, we created a comprehensive database of SSR profiles of Croatian olive germplasm FAZ_OliveDB (University of Zagreb, Faculty of Argiculture Olive Database) (Table S3) using the set of consensual SSR primers [70].The length of the alleles was standardized with the alleles in OleaDB [6], WOGBC [37], and WOGBM [36,38], enabling the comparison of SSR profiles of Croatian origin with SSR profiles of olives in these databases, contributing to the accurate determination of the identity of genotypes from the Adriatic coast of Croatia.
The matching of the SSR profiles of the varieties used as referral samples was performed via the comparison of genotypes on the basis of different SSR data, except for varieties 'Grossa di Spagna' and 'Picholine'.The profile of 'Grossa di Spagna' in the databases OleaDB [6] and WOGBC [37] is identical to the profile of variety 'Konservolia', while in the database WOGBM [36,38] it is identical to the variety 'Amphisis', which is considered as a local synonym of the variety 'Konservolia' [34].The SSR profile of variety 'Picholine' from OleaDB [6] is different from the profile in WOGBC [37] and WOGBM [36] by six alleles, but the difference with the results in our research is five alleles.
In the continuation of the discussion, when comparing the SSR profiles of olive genotypes obtained in our research with the SSR profiles of olive genotypes in the world databases, in addition to the OleaDB [6], we will use the WOGB database [38] because it combines the SSR profiles from both the WOGBC [37] and WOGBM [36] databases.
The SSR profile of the variety 'Karbuncela' (MAR00513) in the WOGB database differs from the SSR profile of 'Karbunčela' found in our studies.According to Bulić [13], the variety 'Karbunčela' is probably identical to the variety 'Carboncella', which was introduced from Southern Italy to the Zadar area, where it has spread in cultivation.According to our analyses, the SSR profiles of varieties 'Karbunčela' and 'Moraiolo' (syn.Carboncella [34]) match, confirming previous assumptions [13] about the origin of the variety and the origin of the name of the variety.By comparing the SSR profile of variety 'Karbuncela' (MAR00513) in the WOGB database with the SSR profile of 'Karbunčela' in our studies, we determined that the SSR profile of 'Karbuncela' in the WOGB database is actually identical to the SSR profile of variety 'Paštrica' found in our studies.Therefore, the accession named 'Karbuncela' (MAR00513) in WOGB database is actually the variety 'Paštrica'.
The SSR profiles of both variety 'Mezanica' (MAR00511) and variety 'Sitnica' (MAR00514) in the WOGB database are identical, but their profile differs from the SSR profile of variety 'Mezanica' tested in our research.The SSR profiles of variety 'Drobnica', together with its clonal variants identified in our research, are very similar to the SSR profiles of varieties 'Mezanica' and 'Sitnica' in the WOGB database.Also, the name of variety 'Sitnica' is a synonym of the variety 'Drobnica' [9].Taking into account these two facts, it is most likely that the WOGB database consists of clonal variants of 'Drobnica', but not 'Mezanica' (MAR00511), whose SSR profile is different in our studies.
The varieties 'Uljarica' (MAR00512) and 'Žabarka' (MAR00505), included in the WOGB database, also have the same SSR profile.However, in our research we determined that it is not Uljarica/Žabarka in WOGB, but 'Bijelica dubrovačka'.In our research, different SSR profiles were determined for 'Žabarka' and 'Uljarica' (Table S3).The SSR profile of 'Uljarica' was identical to the SSR profile of 'Gentile di Chieti', which was the reference variety in our research.The profile of our referral variety was also identical to the SSR profile of 'Gentile di Chieti' (MAR0015) in WOGB database.With the profiles of both accessions, MAR00512 and MAR00505 probably correspond to the variety 'Bjelica dubrovačka'.
The SSR profile of variety 'Oblica' determined in our study matches the SSR profile of 'Oblica' in OleaDB [6] and in WOGB database ('Oblica'; COR000706).The SSR profile of variety 'Lumbardeška' (MAR00500) in the WOGB database is the same as the SSR profile of 'Oblica' in the WOGB database and OleaDB [6].This profile is also identical to the profile of 'Oblica' in our research.According to previous statements [13], the variety 'Lumbardeška' is a synonym of 'Oblica', as confirmed here.
The SSR profiles of varieties 'Istarska bjelica', 'Oblica' and 'Simjaca' (syn. of 'Slimjača' in our research), taken from the databases Olea and WOGB, are identical to the SSR profiles in our research, while the profiles of varieties 'Drobnica', 'Lastovka', and 'Plominka' from these three databases differ in a range from one to a maximum of three alleles per locus compared to the profiles of the same varieties in our research.
The SSR profile of the genotype UNK_56 determined in our study differs in just one allele from both, the SSR profile of variety 'Frangivento' (syn.'Cipressino' [34]) in the OleaDB [6], and the profile of variety 'Cipressino' (COR000090) in WOGB database.Therefore, the genotype UNK_56 from our study is probably variety 'Cipressino', which is sporadically present in Istria, the origin of the genotype for our research.
By comparing the SSR profiles of 135 olive genotypes analyzed in our study with the SSR profiles of olive samples in the OleaDB [6] and WOGB [38] database, we found a total of 92 olive varieties and unknown genotypes with unique genetic profiles in the producing area of the Mediterranean part of Croatia.

Conclusions
This study is the first systematic genotyping of olive germplasm in the Croatian part of Mediterranean.We performed our genetic analyses of Croatian olive germplasm by applying a set of consensus SSR primers, and we conducted our study among a representative number of samples, fully covering the five olive-producing regions in Croatia, results in a comprehensive database of SSR profiles of Croatian olive germplasm "FAZ_OliveDB".The SSR profiles reveal high genetic diversity within the Croatian olive germplasm.Out of a total of 226 analyzed trees, a total of 132 genotypes were determined, including foreign introduced varieties.Wide genetic diversity was revealed among the unknown genotypes represented in production, but genetic uniqueness was confirmed for 73 local varieties and 53 unknown genotypes, latter being presented mainly in extensive production.Three different genetic subgroups were determined by STRUCTURE analysis.The varieties were mostly represented in the genetic subgroup K3, while unknown genotypes were mainly distributed in subgroup K2.Synonyms were found in 18 varieties, but the majority of them were synonymous with the varieties 'Oblica', 'Drobnica' and 'Plominka'.No homonyms were discovered.Small allele variations responsible for clonal variability were discovered in 15 varieties.The lengths of the alleles among databases were standardized in order to enable the comparison of SSR profiles of Croatian germplasm with the olive germplasm from OleaDB and WOGB databases.Out of 132 different genotypes, the profiles of 92 of them are shown for the first time in our database after not being determined in any of the previous reports.These genotypes are local varieties, or genotypes without name that are unique to producing regions in Croatia.Also, several cases of varietal mismatches were detected among the olive databases.Regarding the economic potential of many of Croatian olive genotypes in extensive production, this research is a contribution to their characterization, and is significant for claiming varietal uniqueness.It is an important step for further recognition of the olive oils specific to certain geographic regions or production zone.Knowing the exact genetic profile is important for nurseries in the production of authentic varieties, but also for the introduction and spread of interesting genotypes to wider production.Further work should be focused on the propagation of all the identified genotypes with unique SSR profiles and planting them into collection olive groves, where phenotypic data can be collected with due accuracy.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/horticulturae10040417/s1,Table S1: Final concentrations of components in PCR reactions.;Table S2: The alleles obtained on 11 SSR loci standardized among all three olive databases using referral genotypes.;Table S3: The FAZ_OliveDB (University of Zagreb, Faculty of Agriculture Olive Database) containing the genetic profiles of genotypes of Croatian olive germplasm, including the referral genotypes analyzed on 11 SSR loci from the consensus list of the SSR markers for olive genotyping [70].The table indicates the sites and geographic data for samples collected in Croatia.The last two columns indicate the presence of the same genotypes and their names as present in the databases WOGB and OleaDB, i.e., the synonyms.

Figure 1 .
Figure 1.Dendrogram of 128 olive trees of 16 Croatian olive varieties based on Squared Euclidean distances (i.e., allelic differences) using UPGMA method.Dendrogram represents variation among trees within varieties.

Figure 2 .
Figure 2. The dendrogram: (a) based on the proportion of shared allele distance (Dps) and cluster analysis using Fitch-Margoliash least-squares algorithm of 99 genotypes of Olea europaea based on 11 consensus SSR markers.Olive genotypes were divided into five clusters (C1-C5).Bootstrap values over 50% based on 1000 resamplings of the data set are indicated based on the branches of the dendrogram.(b) Bayesian cluster analysis of the 99 samples of Olea europaea.Genetic subgroups (K) are indicated at the right side of the bar plot.The strips of different colors represent membership coefficient for K = 3 (genetic subgroup K1 (red), K2 (blue), and K3 (green)), calculated using STRUC-TURE v.2.3.4 software.The length of colored strips represents the estimated membership proportion of genotype in corresponding genetic subgroup.Genotypes of mixed origin (with membership probabilities of Q < 75%) are graphically marked by black frame in the dendrogram.

Figure 2 .
Figure 2. The dendrogram: (a) on the proportion of shared allele distance (Dps) and cluster analysis using Fitch-Margoliash least-squares algorithm of 99 genotypes of Olea europaea based on 11 consensus SSR markers.Olive genotypes were divided into five clusters (C1-C5).Bootstrap values over 50% based on 1000 resamplings of the data set are indicated based on the branches of the dendrogram.(b) Bayesian cluster analysis of the 99 samples of Olea europaea.Genetic subgroups (K) are indicated at the right side of the bar plot.The strips of different colors represent membership coefficient for K = 3 (genetic subgroup K1 (red), K2 (blue), and K3 (green)), calculated using STRUC-TURE v.2.3.4 software.The length of colored strips represents the estimated membership proportion of genotype in corresponding genetic subgroup.Genotypes of mixed origin (with membership probabilities of Q < 75%) are graphically marked by black frame in the dendrogram.

Figure 3 .
Figure 3. Bayesian clustering of Croatian olive genotypes.The average log probability value LnP(K) (dashed line) and modal value of ΔK (solid line) or distinctive genetic subgroups K are shown.The ΔK values are based on 20 replicates for each value of K using the ADMIXTURE method.

Figure 3 .
Figure 3. Bayesian clustering of Croatian olive genotypes.The average log probability value LnP(K) (dashed line) and modal value of ∆K (solid line) or distinctive genetic subgroups K are shown.The ∆K values are based on 20 replicates for each value of K using the ADMIXTURE method.Horticulturae 2024, 10, x FOR PEER REVIEW 15 of 23

Figure 4 .
Figure 4. Map of Croatian Adriatic region with the five main olive-growing regions (Istria, Kvarner, northern, middle, and southern Dalmatia) marked.Pie charts representing membership coefficient for K = 3 (genetic subgroup K1 (red), K2 (blue), and K3 (green)), calculated using STRUCTURE averaged by areas designated by growing conditions.

Figure 4 .
Figure 4. Map of Croatian Adriatic region with the five main olive-growing regions (Istria, Kvarner, northern, middle, and southern Dalmatia) marked.Pie charts representing membership coefficient for K = 3 (genetic subgroup K1 (red), K2 (blue), and K3 (green)), calculated using STRUCTURE averaged by areas designated by growing conditions.

Table 2 .
Genetic diversity parameters estimated for 11 SSR primer pairs (loci) in 135 olive genotypes.The parameters reported for each SSR primer pair: the number of different (Na), effective (Ne) and rare or unique (Nu) alleles; the number of distinguished genotypes (NDG); the observed (Ho) and expected (He) heterozygosity; polymorphic information content (PIC); the fixation index or inbreeding coefficient (F); and the probability of identity (PI).

Table 3 .
List of genotypes with rare alleles within Croatian germplasm.The rare alleles unique to Croatian germplasm are marked with asterisk (*).

Table 4 .
Croatian olive varieties analyzed, with the number of sampled trees, common names of varieties, and its synonyms detected as a result of SSR analyses.

Table 4 .
Cont. taking into account the whole Adriatic olive production region, the synonyms are not limited to those stated herein.The synonyms in this table refer only to the location where samples were taken.For e.g., 'Oblica', a variety present at the whole Adriatic region, might have other synonyms in areas outside the ten areas where the samples were taken. 1

Table 5 .
Allele combinations (in bp) within genotype at loci with observed intra-varietal variation.The number of tree samples (n) with defined allelic combinations are presented in brackets.