Genetic diversity and differentiation of Olea europaea subsp. cuspidata (Wall. & G.Don) Cif. in the Hajar Mountains of Oman

Olea europaea subsp. cuspidata (Wall. & G. Don) Cif. is one of the six subspecies important for domestication of olive described as having valuable breeding traits. It is distributed from South Africa to the Middle East and the Mediterranean region to China mainly at mid to high altitudes with adequate precipitation. The currently rather fragmented Olea europaea subsp. cuspidata populations likely represent outposts of a formerly larger and more connected area and are hence vulnerable to additional stresses such as erratic weather extremes and anthropogenic influences. This may result in further fragmentation of stands, shifts in demography, increased clonal growth, lower genetic diversity, and higher differentiation between olive stands. In this study, we present data on Olea europaea subsp. cuspidata growing in the Hajar Mountains, Sultanate of Oman, a large area subjected to many anthropogenic disturbances. We used dendrological traits and genetic parameters to infer its stand and population genetic structure, levels of genetic diversity, and differentiation. To accomplish this, 366 individual trees in total were inventoried and 188 of them genotyped using 12 microsatellite markers. Dendrological measures indicated disturbances mainly due to human activities. Height and crown area of individuals averaged 3.7 m (ranging from 0.2 to 12 m) and 14 m (ranging from 1 to 18 m), respectively. The stand conditional index was mainly 50%. The identified genetic variation within the sampled areas was relatively moderate. The expected and observed heterozygosity equalled 0.55–0.58 and 0.62–0.64, respectively. All individual genotypes were unique indicating the absence of clones. Though significant genetic differentiation was found, there was no evidence of population substructure suggesting effective pollen and long-distance seed dispersal in the past. However, changes in gene flow patterns may be reflected in the next (sibling) generation. Further studies are necessary to detect any changes in the future and to develop sustainable conservation strategies for this species.


Introduction
Plant communities in tropical and subtropical regions are commonly exposed to extreme weather and climatic conditions including high-temperature fluctuations, variable precipitation patterns, and intense radiation. Along with recent climate changes, anthropogenic activities, such as overexploitation, habitat destruction and fragmentation, and introduction of new species, have accelerated ecosystem alterations (Goudie and Viles 2013) resulting in often massive ecological disturbances and fragile biological conditions with unpredictable consequences. Becklin et al. (2016) reported that alteration of precipitation regimes and temperature may disturb the regular patterns of plant morphology, physiology, and life history, which can lead to low plant viability and germination success of seeds (Anderson 2016). Furthermore, disturbances may cause shifts in demography which in turn cause alterations in genetic settings (genetic diversity and differentiation, genetic bottlenecks, gene flow patterns) worldwide (Kassa et al. 2017;MacLaren 2016;Tardieu 2013;Wesche et al. 2005). Responses to changing and harsher environmental conditions can result in a shift in the distribution of mountain plant communities (Grabherr et al. 1994;Lenoir et al. 2008) and a switch from sexual to asexual reproduction (Braatne et al. 1996), respectively, among others.
Such adverse conditions and their effects are also assumed for the open forest community of high altitudinal regions in the Hajar Mountains (the highest area of the Sultanate of Oman with a maximum altitude of 2976 m a.s.l.) dominated by Olea europaea subsp. cuspidata (Wall. & G.Don) Cif. and Juniperus seravschanica Kom. (Zeravschan or Persian juniper). Studies by Al-Kalbani and Price (2015;data range 1979-2012 and Al Sarmi and Washington (2011;data range 1980 indicated that in one decade, mean winter and annual temperatures have risen by 0.85°C and 0.27°C, respectively, while the annual precipitation dropped by 67 mm and 9 mm, respectively. Therefore, in addition to recent and accelerated anthropogenic disturbances such as urban encroachments (housing area expansion, road constructions) and touristic activities (campfires), this mountain range is considered highly vulnerable (Al Charaabi and Al-Yahyai 2013;Al-Kalbani et al. 2014;Patzelt 2015a, b;Victor 2008). Consequently, these changes are likely to affect the stands of Olea europaea subsp. cuspidata, nonetheless clear-cut causes are not yet understood.
Olea europaea subsp. cuspidata is one of six subspecies distributed in Asia, Africa, the Mediterranean region, and Macaronesia with less oil yield and smaller fruits than the economically important and widely cultivated ancient cultigen O. europaea subsp. europaea L. (European olive) in the Mediterranean Basin (Díaz et al. 2006;Rallo et al. 2000). Olea europaea subsp. cuspidata is the most widespread subspecies distributed naturally from the southern to north-eastern part of Africa and from the Arabian Peninsula to China and is restricted to regions with adequate precipitation (± 600 mm) including higher elevation areas. It was introduced into Australia, Hawaii, New Zealand, and Norfolk Island (Green M  . This diploid (2n = 46) and monoecious tree species is mostly outcrossing and wind pollinated, though insect pollination has been also observed (Canale and Loni 2010;Kassa et al. 2017). Seeds are typically dispersed by small mammals and birds (Abiyu et al. 2015;Cuneo and Leishman 2006). So far, clonal reproduction for this subspecies has not yet been proven, though related subspecies demonstrate clonal growth (Baali-Cherif and Besnard 2005). The tree plays an important role as forage and shade tree for agro-pastoral activities, and fruits are eaten by local people (Ghazanfar 2018).
This study was set up to test the hypotheses that (1) remoteness from human activity positively affects stand structure and viability of O. europaea subsp. cuspidata, (2) clonal growth of the species may exist in drier (sloped) conditions promoted by lack of water that suppresses successful germination of seeds, and (3) genetic structure shows bottleneck effects and signs of inbreeding.
Knowledge about this subspecies is insufficient, and the overall goal of this study therefore was to evaluate the vulnerability of the subspecies based on genetic data that can be potentially used to develop conservation strategies for this wild valuable genetic resource.

Study sites and plant sampling
The field study was conducted in Al Jabal Al Akhdar (hereafter named as Akhdar) and Jabal Shams (herein named Shams) area in the Hajar Mountains range of 2110-2507 m a.s.l. (Fig. 1) in October 2016 and approved by the Ministry of Environmental and Climate Affairs (MECA, permit number 15/2016) of Sultanate of Oman. Three plateaus and their adjacent slopes (* smallest sampling units with distances of 0.5-1.0 km between them) were selected per each of the two mountains ((3 plateaus ? 3 slopes) 9 2 mountains = 12 sites in total). The field design furthermore included six nested locations (3 locations (plateaus ? slopes) 9 2 mountains = 6 locations) (Table 1). In addition, the two mountains were analysed as two separate populations. All slope sites selected were north-exposed to reduce bias and to avoid exposure effects on plant growth and stand density. Per site, one reference point was assigned randomly. Around this predefined reference point, a group of the closest 30-32 individuals was recorded ( Fig. 2a, b). Finally, 366 individual trees were sampled. For each individual, spatial coordinates were recorded with a hand-held GPS receiver (Vista HCx eTrex, accuracy ± 2 m, Garmin Ò Ltd., Southampton, Ireland). The following individual dendrological parameters were recorded: height and crown diameter (both in m), number of alive and dead trunks, diameter at breast height (DBH), and diameter at base (both in cm) for individuals smaller than 1.3 m. Tree condition index was used as proxy for tree vitality and estimated based on foliage cover, partitioned into four classes (100%, 75%, 50% and 25%, Fisher and Gardner 1995), knowing that this parameter is rather subjective. In addition, the presence of human disturbance (pruning and cutting) was recorded. Five healthy leaves from each individual were collected and air dried. Out of the 366 individuals, a subset of 188 individuals (15-16 individuals nearest to the reference

DNA extraction and microsatellite marker analysis
Total genomic DNA was extracted from 188 individuals using one air-dried leaf per tree and the DNeasy TM 96 Plant Kit (Qiagen GmbH, Hilden, Germany) according to the manufacturers' protocol.
Data analysis GenAlEx v6.5 (Peakall and Smouse 2012) was used to calculate percentage of polymorphic loci (PPL), the mean number of all (N a ), effective (N e ) and private alleles (N PA ), and observed (H o ) and expected (H e ) heterozygosity (Nei 1973), including the standard error of the mean each. For comparative analysis, the samples were grouped based either on their location (six groups, each including plateau and adjacent slope), plateau versus slope site (two groups), or Akhdar vs. Shams (two groups). Genetic parameters were calculated for all geographic groupings. To check for generation-dependent effects on genetic parameters, age class (assuming increasing DBH as proxy for increasing tree age) quartiles (B 7.2 cm, B 10.7 cm, B 16.1 cm and [ 16.1 cm; n = 160) were created from trees higher than 1.3 m. The GenAlEx software was also used to check for the probability of identity (probability that two unrelated individuals will have the same genotype) and to detect identical multi-locus genotypes (match function). The MICRO-CHECKER software (van Oosterhout et al. 2004) was used to test for the presence of null alleles and to estimate their frequency. To check for usability of markers, Polymorphic Information Content (PIC) was calculated as PIC = 1 -P (P i ) 2 , where P i is the proportion of samples carrying the ith allele of a particular locus. GenAlEx was applied to check for deviations from Hardy-Weinberg equilibrium (HWE). The Genepop 4.3 software (Rousset 2008) was used to calculate F IS values (Weir and Cockerham 1984) for all markers. Significance values were corrected for multiple testing using the false discovery rate (FDR) method in the p.adjust function with R 3.4.3 software (R Core Team 2017). The BOTTLE-NECK v. 1.2.02 software (Piry et al. 1999) was used to test whether observed heterozygosity and allele frequency distribution were affected by bottleneck events in the recent past. Non-excess of heterozygosity assumes L-shaped allele frequency distributions,  (Sefc et al. 2000) b nSSRs (Carriero et al. 2002) c nSSRs (Cipriani et al. 2002) d EST-SSRs (Arbeiter et al. 2017) which are expected under mutation-drift equilibrium and lack of bottlenecks. The program was run for the two-phased model (TPM) of mutation choosing a 0.7 proportion of the stepwise mutation model (SMM) in the TPM and a 0.3 variance of the geometric distribution for TPM with 10,000 iterations. The Wilcoxon rank test was applied to test for significance. AMOVA, fixation indices (F ST , F IS and F IT ), the standardized F 0 ST (Meirmans 2006), Jost's estimate of differentiation D est (Jost 2008), pairwise genetic differentiation with 999 permutations, and Nei's genetic distance between groups were calculated using GenAlEx.
A Principal Coordinate Analysis (PCoA) based on Neis unbiased genetic distance (Nei 1978) was also performed with GenAlEx.
The STRUCTURE v. 2.3.4 software (Pritchard et al. 2010) was used to examine the entire sample for population structure by inferring the number of potential subpopulations or genetic clusters (K) in a sample based on the Bayesian model-based clustering method and considering correlated allele frequencies for an admixture model. The admixture model potentially improves the clustering of closely related species or populations (Lind and Gailing 2013). The program was set to 10 independent runs for each value of potential clusters (K). Each K from 1 to 10 included a burn in period of 50,000 MCMC followed by 10 6 iterations. One of the 12 SSR markers, OeUP-16, likely contained null alleles at high frequency based on the MICRO-CHECKER analysis (see results below) and was hence excluded from STRUCTURE analysis since the presence of null alleles might adversely affect the efficiency of assigning individuals correctly in STRUCTURE (Carlsson 2008).
The STRUCTURE HARVESTER 0.6.94 software (Earl and VonHoldt 2012) was used to extract the most likely number of clusters (K) based on the rate of change in the log probability of data between successive K values as suggested by Evanno et al. (2005). We used the CLUMPAK software (Kopelman et al. 2015) to generate the summary and the graphical presentation of STRUCTURE results.
To check for large-scale spatial genetic patterns, isolation-by-distance (Mantel test with 999 permutations) was assessed with GenAlEx. Fine-scale spatial genetic structure (FSGS) was analysed for each site by spatial autocorrelation analysis using SPAGeDi v. 1.5 (Hardy and Vekemans 2002). To quantify FSGS, the Sp statistic (Vekemans and Hardy 2004) was calculated as Sp=-b F /(1 -F 1 ), where b F is the regression slope of F(d) over the natural logarithm of distance, and F 1 is the average kinship coefficient referred to as a co-ancestry estimate (Loiselle et al. 1995) based on all loci between pairs of individuals of the first distance class. The kinship coefficients were plotted against 10 m distance classes. The distance classes were defined so that 50% of all trees were represented at least once in the interval, and the coefficient of variation (CV) of the number of times when each individual was represented was B 1. Upper and lower 95% confidence intervals of kinship coefficients were determined through 10,000 permutations of individuals among distance classes.

Dendrological parameters
The height and crown width of individuals averaged 3.7 m (ranged from 0.2 to 12.0 m) and 14 m (ranged from 1 to 18 m) (Table 3), respectively. The number of alive and dead trunks varied only marginally from three to four (Table 3). DBH size classes showed a bell-shaped distribution, which was slightly rightskewed due to a few very old trees (Fig. 3). Trees smaller than 1.3 m showed an L-shaped distribution and were heavily browsed, resulting in a dwarfed habitus. DBH and diameter at base varied strongly among sites and locations, although less between mountains (data not shown). The plateau sites exhibited larger DBH values then their adjacent slopes, a trend that was not observed for diameter at base. Only one individual had fruits. The condition index averaged in total 50%, while trees from plateaus of Shams showed the lowest vitalities. In terms of cuttings (branches and stems), human influence was pronounced on Akhdar, while almost absent on Shams (Table 3).

Genetic richness, diversity, and differentiation
In total, 147 different microsatellite alleles were found for all 12 microsatellite markers. Probability of observing identical genotypes by chance was extremely low (4.3 -13 ), and no identical multi-locus genotypes were found in the entire sample. The mean number of alleles per locus was 12.3 with the lowest number of alleles (N a = 2) at both EST-SSR loci OeUP-42 and -45 and the highest number of alleles (N a = 31) at locus ssrOeUA-DCA9 (Table 4). The effective number of alleles per locus averaged 5.6. The observed heterozygosity (H o ) ranged from 0.03 (OeUP-45) to 0.89 (ssrOeUA-DCA18) with a mean value of 0.57. The average expected heterozygosity was higher than the observed one (H e = 0.64 vs. H o = 0.57). Fixation index (F IS ) and genetic differentiation (F ST ) for the 12 markers averaged 0.09 (ranging from -0.41 to 0.15) and 0.02 (ranging from 0.02 to 0.04), respectively (Table 4). The standardized differentiation measures F' ST and Josts D est averaged 0.03 (ranging from 0.003 to 0.18) and 0.02 (ranging from 0.0001 to 0.17), respectively (Table 4). PIC ranged between 0.28 and 0.82 (Table 4).
The HWE tests per locus and in the six locations showed overall low indication of homozygote excess (about 25% of the markers demonstrated homozygote excess, mostly statistically insignificant and likely due to null allele presence) ( Table 5). The OeUP-16 marker indicated a pronounced excess of homozygotes with significant deviation from HWE in all six locations, except for Akhdar 2 ( Table 5). The ssrOeUA-DCA4, ssrOeUA-DCA9, GAPU71, UDO99-011, and OeUP-41 markers were only occasionally indicating deviations from HWE and null allele excess. Deviations from HWE and excess of homozygotes were observed for the same markers also for the mountain and site groups. F IS ranged from -0.064 to 0.396.
Microsatellite loci were polymorphic for any grouping ( Table 6). The average number of alleles (N a ) per grouping was nearly twice as high as the number of effective alleles (N e ; Table 6). The number of private alleles (N PA ) ranged between 3 and 12 (location), 36-51 (mountain), and 24-41 (site)  (Table 6). Genetic diversity parameters declined with increasing age, while number of private alleles did not follow this trend (Table 6). An L-shaped distribution of allele frequencies was observed for any grouping, indicating the absence of pronounced bottleneck effects (data no shown). Pairwise Neis genetic distance ranged from 0.03 to 0.09, being highest for Shams 3 against all others ( Table 7). The F-index values were low, but differed significantly from zero for any grouping (Table 8). While F IS and F IT between groupings differed only marginally, the F ST value between the two mountains (Akhdar and Shams) and among the six locations was about two and four times higher, respectively, as compared to plateau and slope sites (Table 8). However, the F ST values were very small in general.

Population structure
The PCoA did not reveal any clear substructure among the sampled individuals (Fig. 4) demonstrated by the low discriminatory power of the first two coordinates that explained 12.9% of the total genetic variation observed. Together with the third coordinate they still explained only 17.8% of the total variation. These data were in agreement with the STRUCTURE analysis results, indicating little or no genetic structure. The most likely number of clusters (K) was eight, although less likely it could also equal four or six (Fig. 5). However, most individuals expressed some degree of balanced share of cluster assignments irrespective of the actual number of clusters and grouping. Therefore, neither PCoA and STRUCTURE nor F ST values revealed any obvious population substructure or strong differentiation among O. europaea subsp. cuspidata trees.  Large-scale spatial analyses based on a Mantel test revealed a low but significant correlation between genetic and geographic distances (r = 0.139, p = 0.04). Fine-scale spatial analyses revealed a positive significant spatial autocorrelation at the Akhdar 3 slope and Shams 2 plateau sites in the first (0-10 m) and second (10-20 m) distance intervals, respectively (Fig. 6). A negative significant autocorrelation was observed for the Shams 2 plateau site at the third (20-30 m), for the Akhdar 3 slope site at the seventh (70-80 m), and for the Shams 3 slope site at the second (10-20 m) distance intervals. Overall, the observed patterns demonstrated no or weak autocorrelation.

Dendrology
The obtained dendrological data indicated clearly location and site-specific responses. Signs of human activities such as pruning of branches and logging of trees on Shams seemed to be less visible, which can be historically explained: until recently, only few and small settlements exist on Shams. This could be attributed to a scarce water availability and rather steep terrain (higher water run-off) in comparison to the neighbouring Jabal Al Akhdar. The generally sinking groundwater levels in the region (Al-Kalbani and Price 2015) and the overall more severe water approachability on Jabal Shams might thus also affect the existence of O. europaea subsp. cuspidata, visible in pronounced die-backs on Shams. Additionally, livestock browsing was apparently higher on Shams considering the number of dwarf shrubs (diameter at base) as proxy for high livestock intensities. This observation might be related to the slightly different livelihoods on both peaks; Shams serves predominantly as a rangeland area for the piedmont villages, while mountain villages on Akhdar traditionally practice terrace gardening with nowadays increasing off-farm activities such as tourism and administrative jobs.
The bell-shaped distribution of size class[ 1.30 m clearly indicates an ongoing over-aging process of the  (Bekele 2000) with typical L-shaped age class distributions indicative for relatively high proportions of young trees. In our study, trees with heights smaller than\ 1.30 m indicate high recruitment rates (L-shape distribution). The older aged individuals with heights smaller than \ 1.30 m were also present, but have been kept short due to browsing.
Given the different ecological conditions on both mountains, hypothesis 1 cannot be satisfactorily answered, as human influence seems higher on Akhdar, while livestock significantly affects stands on Shams.

Genetic richness and diversity
Although clonal propagation would be potentially possible and the sampling design should allow us to detect it, we did not find any evidence of vegetative propagation. Similarly, Kassa et al. (2017) did not find clonal colonies (genets) of O. europaea subsp. cuspidata in Ethiopia. However, Baali-Cherif and Besnard (2005) identified genets of O. europaea subsp. laperrinei in Algeria based on nine nSSRs, of which ssrOeUA-DCA3 and -DCA 9 have been also used in our study. The rate of clonal reproduction in Algeria was relatively high -16% of genets were represented by two or more ramets and the area covered was more than 80 m 2 in one of the locations. The absence of genets in our sample indicates an exclusive sexual reproduction mode in this population, probably also for the entire subspecies. This observation would merit further studies in additional stands of the Hajar Mountains as well as other populations, for instance  in Dohfar (S-Oman) and neighbouring countries. Therefore, the clonal propagation proposed in hypothesis 2 must be rejected.
Genetic diversity was rather moderate (H o-= 0.55-0.58, H e = 0.62-0.64), although slightly lower as compared to the same species in an Ethiopian forest studied by Kassa et al. (2017) using nSSR    Loiselle et al. (1995) with 95% confidence interval (grey area) along spatial distance (m) for each site of 188 Olea europaea subsp. cuspidata trees in the Hajar Mountains. Circles indicate significant autocorrelation (see geographic data on the sites in Table 1) (Kassa et al. 2017) (H o = 0.76, ranging from 0.40 to 0.92 and H e = 0.79, ranging from 0.52 to 0.93) for the same subset of EST-SSR markers suggesting possibly a lower adaptive variation for the populations under study. Nevertheless, our data on genetic diversity are in agreement with the SSR data found in other published studies of related subspecies. The observation of decreasing diversity with increasing size classes, however, was unexpected and has been to our knowledge only rarely observed as for instance for Quercus rubra in southern New England (USA) by Gerwein and Kesseli (2006). Authors suggested some recovery of alleles when forest cover increases, which might be true also for our study indicating a formerly stronger fragmented population with reduced gene flow subjected to loss of genetic diversity. Alternatively, selected genotypes with augmented growth could have resulted in low genetic diversity with increasing size. Our results, however, do not provide evidence for selected genotypes with comparable genetic composition. Moreover, since we have used DBH as a proxy for age, but have not determined tree age by annual ring counts, these observations remain vague and, if they are not accidental, could have been caused by other factors.
Genetic differentiation and population structure Genetic differentiation among geographic groups was very low in general, but significant in all cases. The low differentiation among the six sampling locations (F ST = 0.012, F 0 ST = 0.036) and even lower differentiation between sites (F ST = 0.004, F 0 ST = 0.009) and mountains (F ST = 0.007, F 0 ST = 0.020) suggest effective gene flow (at least in the recent past) and a still effectively large population size, supported by an absent or at least a weak fine-scale genetic structure. As not all markers demonstrated significant excess of homozygotes, the fixation index (F IS ) cannot serve as a proxy for inbreeding.
The low genetic differentiation measures found for the samples can be explained by long-distance dispersal. Extensive wind pollination and long distance seed dispersal by mammals and birds may ensure comparatively high genetic diversity even in isolated mountain populations of O. europaea subsp. cuspidata as likewise suggested for related O. europaea subspecies (Baali-Cherif and Besnard 2005;Besnard et al. 2007Besnard et al. , 2009Besnard et al. , 2013Breton et al. 2006;Kassa et al. 2017). Potential mammal vectors for O. europaea subsp. cuspidata seeds (fruits) dispersal in the Hajar Mountains are livestock (goats and sheep) and feral donkeys nowadays. Livestock species very likely contributed to the dispersal activity in the last decades and still do. However, other taxa such as birds of the genera Blackcaps (Sylvia sp.) and Thrushes (Turdus sp.) are known to be likely effective seed dispersal agents for O. europaea (Rey et al. 1997). Although a reduced gene flow can be expected due to a rough terrain, absence of sufficient flowering trees, and the fragmentation of stands (Abiyu 2012;Kassa et al. 2017), we did not find any environmental or demographic effect on subpopulation structure. Even a * 20 km distance between the mountain peaks seems to have no effect on the substructure. Perhaps, it is due to the connection via a mountain crest that hardly goes beyond 1800 m. Wadis (deep valley occasionally filled with water) between mountains are narrow and are unlikely to be major barriers for anemochory and zoochory. This finding was furthermore substantiated by the low correlation between genetic and geographic distances. However, this correlation was significant demonstrating that gene flow is or was relatively limited and that isolation by distance might have played (and certainly still does) an important role in differentiation of olive populations.
In comparison to other population genetic studies of olive subspecies that mostly employed random sampling approaches, an effect of our sampling design on the obtained results should also be considered. Although no closely related or clonal trees have been found, neighbouring trees were slightly more related than more distant trees based on the spatial autocorrelation analysis. Finally, the STRUCTURE analysis did find an ambiguous number of clusters, which can be interpreted as a lack of a strong population substructure. A large individual genetic admixture and rather arbitrary assignment of individuals to different clusters mean no obvious substructure except for the Shams 3 slope site that behaved slightly different. However, it is known that population structure algorithms in BAPS, GENELAND and STRUCTURE computer programs tend to overestimate K, especially when isolation-by-distance is strong (Frantz et al. 2009;Guillot et al. 2005;Pritchard et al. 2000). Therefore, the low, but nonetheless significant isolation-by-distance found in our study can partly explain the rather large number of clusters inferred in our sample by STRUCTURE analysis, although those numbers could be rather illusive. Nevertheless, simulations by Guillot and Santos (2009) found that such interferences are stronger particularly in populations with small dispersal distances and spatially irregular patterns (clumped sampling as ours). Moreover, Pritchard et al. (2010) even stated that they permit allele frequencies to be correlated across populations in STRUCTURE to allow accurate assignments of individuals in very closely related individuals, though possibly at the cost of making it likely to overestimate the number of clusters. Hence, we should not overemphasize the number of clusters, especially when population structure and overall differentiation measures are low. However, effects of more recent fragmentation events could be reflected in the next offspring generation suggesting the need for additional studies, particular with respect to a very low recruitment and die-backs that seems to be chronic for the population in the Hajar Mountains. Finally, hypothesis 3 must be rejected as no evidence on genetic bottlenecks and signatures of inbreeding were found.

Conclusions
We could not identify the driving forces of low regeneration and die-back of trees, but suggest cumulative abiotic and biotic effects as the most plausible cause for demographic change. Remoteness alone appears to be not an explanatory factor for the viability of stands. Therefore, we failed to answer hypothesis 1. The absence of identical multi-locus genotypes led us conclude that clonal growth is non-existent at least for the studied individuals and reject hypothesis 2. The studied populations contained relatively high levels of genetic variation and low differentiation based on both nSSR markers that are likely selectively neutral and EST-SSR markers that could be under selection. Overall, bottleneck effects and inbreeding expected under the prevalent conditions have not been found prompting us to reject also hypothesis 3. However, comparatively high levels of genetic variation do not preclude that this population could become vulnerable, especially considering their fragmentation, the lack of regeneration, over-aging, and ongoing climate variations. The unexpectedly higher diversity found in younger size classes merits further studies of sibling generations to better address the questions regarding sustainability of this population, and whether they need special protection.
Funding Open Access funding enabled and organized by Projekt DEAL. There was no specific funding source.
Data availability Data can be shared upon request.

Compliance with ethical standards
Conflict of interest The authors declare there are no competing interests.
Ethics approval Not applicable.
Informed consent All authors agreed on the publication of the present work.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.