Strong Genetic Differentiation of Submerged Plant Populations across Mountain Ranges: Evidence from Potamogeton pectinatus in Iran

Biogeographic barriers for freshwater biota can be effective at various spatial scales. At the largest spatial scale, freshwater organisms can become genetically isolated by their high mountain ranges, vast deserts, and inability to cross oceans. Isolation by distance of aquatic plants is expected to be stronger across than alongside mountain ridges whereas the heterogeneity of habitats among populations and temporary droughts may influence connectivity and hamper dispersal. Suitable aquatic plant habitats became reduced, even for the widespread submerged Potamogeton pectinatus L. (also named Stuckenia pectinata) giving structure to various aquatic habitats. We compared the level of genetic diversity in a heterogeneous series of aquatic habitats across Iran and tested their differentiation over distances and across mountain ranges (Alborz and Zagros) and desert zones (Kavir), with values obtained from temperate region populations. The diversity of aquatic ecosystems across and along large geographic barriers provided a unique ecological situation within Iran. P. pectinatus were considered from thirty-six sites across Iran at direct flight distances ranging from 20 to 1,200 km. Nine microsatellite loci revealed a very high number of alleles over all sites. A PCoA, NJT clustering and STRUCTURE analysis revealed a separate grouping of individuals of southeastern Iranian sites and was confirmed by their different nuclear ITS and cpDNA haplotypes thereby indicating an evolutionary significant unit (ESU). At the level of populations, a positive correlation between allelic differentiation Dest with geographic distance was found. Individual-based STRUCTURE analysis over 36 sites showed 7 genetic clusters. FST and RST values for ten populations reached 0.343 and 0.521, respectively thereby indicating that allele length differences are more important and contain evolutionary information. Overall, higher levels of diversity and a stronger differentiation was revealed among Iranian P. pectinatus than previously observed for temperate European regions, due to regional differences across mountain ranges over long distances.


Introduction
Biogeographic barriers for freshwater biota can be considered at three spatial scales, namely inter and within basin and continental. At the largest spatial scale, freshwater organisms are isolated by their high mountain ranges, inability to cross oceans and vast deserts [1]. Gene flow between populations of aquatic plants is restricted by the discontinuous nature of their habitat embedded in another terrestrial landscape and as such, lakes or rivers could be considered as islands of suitable habitats. Under an Island Migration Model, increasing geographical distance between populations is expected to lead to enhanced genetic isolation, which essentially corresponds to a stepping-stone population structure and requires that dispersal is primarily local rather than long-distance [2,3,4,5,6].
Besides merely distance as a limit to connectivity, especially large deserts and mountain ranges act as additional natural barriers. The flora of the mountains provides many examples of subspecies with populations and very closely related species that show restricted gene flow due to these inherent geographic barriers [7,8]. Isolation by distance is stronger across than alongside mountain ridges [9,10,11].
Drought as a climatic event is one of the most threatening challenges that occur due to climatic change and global warming [12]. The heterogeneity of habitats and environment among populations probably influence gene flow in disrupting or dispersal, resulting in isolation by environment (IBE) model in which genetic differentiation is positively correlated with environmental differentiation [13,14]. Drought has been a problematic phenomenon for aquatic ecosystems in Mediterranean temporary pools [15] and in low rainfall regions e.g. Central and SW Asia [16] during recent years because it declines water quality and decreases aquatic plant occurrences.
The importance of submerged vegetation in aquatic ecosystems management and biodiversity has been highlighted [17,18,19]. Many aquatic plants have mixed modes of reproduction and different models of dispersal [20,21]. In submerged water plants (macrophytes), sexual reproduction by seeds produces new genotypes [21]. In spite of the low clonal diversity between and within populations [20], several studies showed that aquatic plant populations often had high clonal diversity [22,23,24,25] nearly similar to non-clonal plant species [26]. Among aquatic plant taxa, fennel pondweed, Potamogeton pectinatus L. is one of most diverse submerged aquatic species [27,28]. This species maintains many useful physiological traits such as tolerance to a wide range of nutrients, ability to grow in oligotrophic to eutrophic waters [29], and utility to progress quality of water by removing nutrients [30]. These useful traits made the species also a potentially important biotic tool for cleaning up polluted waters by heavy metals absorption [31,32]. Plant parts are an important source of food for many waterfowl [23].
P. pectinatus reproduces both sexually and vegetatively through tubers emerging from the rhizomes and acting as propagules. Asexual reproduction is thought to be responsible for short distance dispersal while sexual reproduction with seeds is more important to ensure long distance dispersal and long-term survival [21,29]. There are several studies about genetic variation of P. pectinatus. Hettiarachchi & Triest [33] used isozyme markers to detect genotypes within different habitats such as freshwater and brackish water, and showed a geographic trend in genotypes among freshwater populations. Using isozyme polymorphism, Hollingsworth et al. [2] confirmed the importance of local clonal growth in P. pectinatus. Hangelbroek et al [23] used RAPD markers and showed that sexual reproduction has a very important role within a lake population. Hollingsworth et al. [2] also showed that isolation by distance was determined as a regulatory factor for gene flow at distances of more than 1000 km along the Atlantic and Baltic Sea coastal habitats, proposing that dispersal between populations is mainly by seeds. The effect of migrating waterfowls on plant distribution pattern and seed dispersal over large distances have been reported [22,34,35]. Also Mader et al. [22] and King et al. [34] who used RAPD and ISSR markers, respectively concluded that genetic distance between populations of P. pectinatus increased with geographic distance. Following the development of nine polymorphic microsatellite loci for P. pectinatus [36] an evolutionary divergence and less gene flow between two types of habitat on a regional scale [37] were shown, namely for the brackish water populations of the Baltic Sea lagoons versus freshwater inland populations of Central Europe. Genetic differentiation of P. pectinatus but also other submerged plant taxa was found between coastal brackish water and freshwater inland habitats in lowland regiosn of W. Europe by Triest et al. [38]. They also detected that upstream forest ponds can be detected as source populations and defections for clonal diversity to recolonize the stressful downstream river habitat. A more detailed study of the clonal diversity and fine-scaled spatial structure at individual level of P. pectinatus in river populations and managed pond [39] revealed that the ponds populations had a higher amount of gene and clonal diversity than those of rivers. Han et al. [40] using AFLP showed that P. pectinatus maintained a high level of variation within and between two contrasting and distant lakes in China. Those previous studies over the last decades confirmed that P. pectinatus has a mixed reproduction system, an overall high level of gene diversity and an increased differentiation over long distances of hundreds of kilometers. All previous studies were done in temperate and temperate-cold areas having many permanent water bodies and often at close vicinity. Therefore, it can be hypothesized that P. pectinatus populations might exhibit a stronger genetic differentiation over similar distances across mountain ranges and in a drier climate zone with often temporary habitats of very isolated wetlands areas and rivers habitats that show strong fluctuations due to high evaporation levels and use for irrigation in agriculture.
The specific biogeographic location of Iran [41] with its several large mountain ranges and deserts present a totally different suite of habitat types at a presumed low connectivity. The Alborz mountain ridge is located in the north whereas the Zagros mountain region is extending from northwest to southeast of Iran. The mountains of Kerman are located at the eastern edge of Zagros. Dashte Kavir is a desert region in the center of Iran. These geographic characteristics are assumed to create a differentiation and isolation between north and south populations and thus can pose an additional barrier to dispersal besides the distance to overcome. Along the same side of mountain ranges, these biogeographic barriers promote homogenization of aquatic biota [42]. Additionally, the mountains of Kerman located in the southeast of Iran are a potential barrier between north and south populations. Historically, these barriers identified as dominant factor in determining the composition of regional organisms [43,44]. The world's second largest saline lake, Urmieh lake, is located in the north west region. Iran has a great diversity of aquatic ecosystems many of which have been recorded in the Ramsar Convention [45]. The largest portion of freshwater bodies in Iran is composed of many large rivers and lagoons. Different challenges such as the effects of human-based threats disrupt their ecological characteristics [46].
Our aim was to test, by considering sites with P. pectinatus populations over a long gradient of regions across mountain ranges as potential barriers, whether in such habitats 1. P. pectinatus can maintain its usually high levels of allele and gene diversity.
2. Differentiation over distances and across mountain ranges and around deserts would be more pronounced and structured than in temperate regions.
Therefore we considered a sampling design with a large number of sites across the various biogeographic regions and used the same microsatellite loci as a previous European studies, following a preliminary taxon verification with nuclear ITS and maternal chloroplast marker genes.

The species
Fennel pondweed, P. pectinatus L. (= Stuckenia pectinata, see Lindqvist et al. [47] for taxonomic revision) is a submerged aquatic plant with a nearly cosmopolitan distribution ranging from sub-arctic to tropic regions, of different types of waterbodies and adapted to a large variation in water depth, flow regime, trophic status and salinity [48,37]. Potamogeton pectinatus reproduces both sexually and vegetatively through propagules emerging from the rhizomes [29]. It is an important source of food for many waterfowls [23]. This species harbors many useful physiological traits such as tolerance to a wide range of nutrients, capability to grow in oligotrophic to eutrophic waters [29], and ability to improve quality of water by removing nutrients [30]. These useful traits made the species a potentially important candidate for cleaning up polluted waters by heavy metals absorption [31,32].
Study area and sampling design P. pectinatus plants were collected in 2015-2016 in thirty six sites from wetlands, lakes and rivers in Iran (Table 1, Fig 1). No specific permission was required for these locations/activities and field studies did not involved endangered or protected species. The geographic distance between pairs of populations ranged from 20 to 1200 km. Because the populations of this species are still declining due to human activities and recent drought in Iran; in each site we collected 1-20 individual shoots (ramets) at 2-3 m intervals. Design was to cover the entire distribution range and diversity of aquatic habitats within Iran. Our samples usually were in the main channel in downstream parts of rivers with low speed.
Samples were dried on silica gel prior to extraction and a voucher specimen for each population was deposited at HUI (Herbarium of the University of Isfahan). One subset (one individual from one sampling site) of the collection was used for cpDNA haplotypes and ITS (accession number: DQ840279.1) characterization (results not shown in this paper). Identification of P. pectinatus populations was based on cpDNA sequence variability in trnH-psbA referring to P. pectinatus haplotypes and on two nuclear spacers (ITS1 and ITS2)-as in Triest and Sierens [49]. A more detailed study of genetic diversity and structure at the population level was done on a subset of 10 populations where a sufficient large number of different genets could be considered.
For each specimen we measured 29 morphological characters. Among these, the leaf apex and leaf width (36 measurements per site) could differentiate morphologically the specimens characterized for a nuclear ITS-A or B. B. Mean rank for Krustal-Wallis was between 1.12-3.11 mm for the group of specimens sharing ITS-B. Morphological measurements were tested for significant differences with one-way ANOVA (Kruskal-Wallis) and pairwise Mann-Whitney U test.

DNA extraction and microsatellite amplification
Genomic DNA was extracted from leaf tissue using CTAB by Gawel and Jarret [50] with some modifications. All individuals were genotyped at nine nuclear microsatellite loci: Potpect24, Potpect26, Potpect28, Potpect32, Potpect34, Potpect37, Potpect39, Potpect40, and Potpect42 [36]. Amplification of these nine loci was not performed in three multiplex PCR as described in Nies & Reusch [36], but into a unique multiplex PCR. This single amplification   [36,37]. To explain this phenomenon, these authors support the idea that the time that has elapsed after the polyploidisation event was sufficient to cause genetic divergence of the microsatellites on the different sets of homologous chromosomes [37]. After excluding forty eight repeated multilocus genotypes (MLGs) from the data set, we calculated genetic diversity measures on 133 individuals using FSTAT 2.9.3 [51] and GenAlex 6.5 [52]. At locus level, we estimated the number of alleles (A), effective number of alleles (Ae), observed heterozygosity (Ho), expected and unbiased expected heterozygosity (He and uHe), Nei's heterozygosity at subpopulation level (Hs) and total (Ht), Wright's F statistics [53] with inbreeding coefficient (F IS ), total inbreeding (F IT ), subpopulation differentiation (F ST ) and gene flow (Nm) using GenAlex. Weir & Cockerham [54] estimation of F IT (CapF), F IS (small f) and F ST (Theta) following Rousset [55] were calculated with FSTAT. Inbreeding values for the total populations are given by small f and jackknifed over loci (FSTAT).  We tested for recent bottlenecks in each population under the two-phase model (TPM) with 95% single-step mutations and 5% multiple-step mutations (Wilcoxon's test 1-tailed) using bottleneck 1.2.02 [56]. Genetic structure was assessed with a three-level analysis of molecular variance (AMOVA) [57] using hierarchical and standardized fixation indices for 5 regions, 10 populations and 68 individuals (GenAlex). Pairwise FST-values were calculated between all pairs of populations and tested for significant differentiation using 999 permutations. Isolation-by-distance between pairs of populations and their geographical distances was tested with ϴ/1-ϴ considering straight flight distances, log transformed [55] between populations in a Mantel test using 1000 randomizations [58]. An estimator of actual differentiation D est [59] was calculated between all pairs of populations using SMOGD [60] and used in a Mantel test with log transformed straight flight distances. To infer the population structure on basis of assignment of individual genotypes into groups, a Bayesian clustering method [61] was carried out using STRUCTURE version 2.3. We tested for the number of clusters (K) in ten independent runs from 1 to 14 (10,000 burn-in, 10,000 Markov chain Monte Carlo replicates in each run), without using sampling site as a prior to assess convergence of the estimated ln probability of data, ln (PD). Runs were carried out assuming admixture and an independent model of allele frequencies. The number of clusters was determined from the K with the highest posterior probability and using the second-order rate of change of the likelihood function ΔK, as suggested by Evanno et al. [62]. For testing the evidence of scoring error, Evidence of large allele dropout and Evidence of null allele we used Micro-checker 2.2 software. Also we used POPTREE2 [63] software for cluster analysis. Isolation by distance and R ST was tested by calculating multilocus estimates of kinship coefficient (Fij) between all pairs of individuals implemented in the software SPAGeDi v1.4 [64].

Microsatellite loci properties
As the primers were previously developed and used on European Potamogeton pectinatus populations, we tested these cross-amplified loci with microchecker, but there was no evidence for null alleles, large allele dropout or scoring error due to stuttering. Over all sites and loci, the nine microsatellite loci revealed a total of 130 alleles. The mean number of alleles was highest for potpec 39 (3.3) and lowest for potpec24 (1.9) with an overall observed heterozygosity (Ho = 0.2-0.7) nearly similar to the overall expected heterozygosity (He = 0.3-0.7) ( Table 2).
From the 10 considered populations, the percent of individuals assigned to the same population was higher (76%) than those assigned to other populations (24%). The populations of Amirkelaieh (2), Azbaran (3) and Kaniborazan (21) did not share any allele.
A PCoA at the level of all individuals revealed a gradient along a first axis that explained 28.9% of all variation but where several individuals were clearly separated as a group (Figs 3  and 4). Most of these individuals were also tested for ITS and showed the ITS-B sequence. A PCoA confined to ten populations that explained 28.8% of all variation revealed a very clear separation of most populations from each other. A gradient can be observed from individuals of N. Alborz, around lake Urumieh, S.E. Zagros, N. Zagros and SE Zagros to Kerman mountain ranges.
Percentage of variation explained by the first three axes in PCoA at 10 populations is 28.80% and 17.94% for first and second axes, respectively.
The range of F ST was between 0.156-0.595 among the ten populations. The lowest F ST was between Izeh (19) and Amirkelaieh (2) and the highest was between Azbaran (3) and Yaschaman (35) ( Table 4).
The range of D est values were between from 0.001 to 0.925. There was no relationship at all between Fst and D est (Pearson correlation = 0.007) indicating that allelic differentiation can be more important as an explanatory variable of IBD. There was a positive relationship between F ST and geographic distance (R 2 = 0.35 and P = 0.01) but there was an improved positive correlation between D est and Ln (1+Km) (P = 0.01, R 2 = 0.56) (Fig 5). An AMOVA showed that the percentage of molecular variance among regions and populations were 15% and 18%, respectively ( Table 5).
The genetic structure of 10 populations of Potamogeton pectinatus in Iran showed 6 genetic clusters (Fig 6). The samples from Yaschaman (35) are placed separately. Genetic analyses for 10 populations at 5 distance classes showed that F ST and R ST values were 0.343 and 0.485, respectively. Table 6 indicated special genetic structure analysis for these classes.
The higher values for R ST indicated that allele length differences are more important over the entire study area and evolutionary significant through low historical connectivity. In the second distance class, F ST and R ST had a similar value (0.32-0.33). In the largest distance class, the R ST value (0.698) and F ST value (0.496) were highest. Also in all distance classes, except second distance class, the values for R ST were significantly higher than for F ST . F ST for all loci

Allele and gene diversity
In spite of preliminary studies on the vegetative mode of reproduction of P. pectinatus and previous genetic studies on submerged aquatic plants based on dominant markers [29,33,2], as well as in other studies on co-dominant markers [36,39] a high level of allelic and gene diversity was shown. Several studies showed that aquatic plant populations often had high clonal diversity [22,23,24,25] nearly similar to non-clonal plant species [26]. Among aquatic plant taxa, fennel pondweed, P. pectinatus L. is one of most diverse submerged aquatic species [27,  28]. Across the wide study area we observed 130 different alleles for 133 individuals whereas Triest et al [39] and Nies & Reusch [36] observed 56 and 65 alleles with microsatellites, respectively in a much larger amount of individuals. This can be explained from the geographic distances that reached only up to 250 km whereas our study covered sites at distance up to1200 km, potentially leading to detection of much more allelic variants and a greater isolation by distance between populations. The lower number of samples per site clearly allowed to obtain many more alleles than in temperate regions and to detect relevant genetic structures in agreement with biogeography and mountain ranges, the sample size did not prevent us from observing the strong differentiation (with more samples/site, alleles would remain different across sites and regions). Other researchers also used the low sample size for example Rodriguez-Bonilla et al [65] in genetic diversity evaluation of sweet potato in Puerto Rico used 137 landraces, Population Genetics study of the Sao Tome Caecilian (Gymnophiona: Dermophiidae: Schistometopum thomense) revealed strong geographic structuring with 138 specimens by  Stoeltin et al [66], Jackrel and Wootton [67] used low sample size in assessment variation of riparian plants within and among species shapes river communities. Under an Island Migration Model, increasing geographical distance among populations is expected to lead to enhanced genetic isolation, which essentially corresponds to a steppingstone population structure and requires that dispersal is primarily local rather than long-distance [2,3,4,5,6]. Maximum geographic distance can have a strong effect on among population diversity [68]. With regard to genetic differentiation, observed F ST value in this study was 0.5  that was relatively high due to spatial differences of studied populations and stochastic founder events as in Triest et al [39] and Nies & Reusch [36]. However, the allelic differentiation (D est ) was much more pronounced because of the large amount of alleles detected across the different areas. Also the He value in our study is higher than in central or western Europe [36,39]. Genotypic variation and recombination in P. pectinatus due to sexual reproduction is most likely explaining the observed high genetic diversity [22,23]. The mixed mode of reproduction in P. pectinatus has an important effect on its genetic structure [20] and we also observed some clonal repeats among ramets within a site. P. pectinatus is a polymorphic species with high phenotypic plasticity [28]. As a phenotypic plasticity phenomenon, the plants may produce more seeds that can lead to high genetic diversity and heterozygosity [69,70].

Geographic differentiation/IBD?/Connectivity
With regard to our results, we found strong genetic differentiation between Iranian sites of Potamogeton pectinatus. Mountains act as strong genetic barriers but river valleys recognized as corridors for gene flow [10]. The distribution of Iranian Potamogeton seems to be linked to mountain regions as indicated by the STRUCTURE results (for 10 populations) (Fig 6) and by the PCoA at individual level (for all 36 sites) (Fig 3). We obtained clear groups according to their 'Mountain positions': NE Alborz (3); NW Alborz (2); alongside of Urmieh lake (21,24); SW Zagros (19,27,30), SE Zagros (22), N Zagros (34) and Kerman (35). Also allelic differentiation (D est values) showed very low or absent connectivity across mountains north (Albroz) and south (Zagros/Kerman). A strong correlation between genetic and geographic distances revealed a pattern of IBD across the distribution range of P. pectinatus in Iran. We found a strong positive correlation between D est and geographic distances (Fig 5). The differentiation was larger for populations separated by the Alborz mountains than those separated by the Zagros mountains.
There is however some connectivity between both sides of Zagros mountains namely between SW Zagros (19,25,27,30) and N Zagros (4, 5, 6, 7, 8, 9, 18 etc.). This area in Iran was reported as a route for bird migration [71]. It thus can be assumed that e.g. between Izeh and Isfahan this route forms a ridge between high Zagros and low Zagros. Equally important could be the wider pass with upstream rivers reaching geographically rather close towards the watershed. A connectivity over longer distance exists between Vimcheh in Isfahan (N. Zagros) and Yaschaman (Kerman) with the corridor on the same side of the mountain area as a probable explanation (= ITS-B and cpDNA-380), together with the bird migration maintaining this connectivity. The connectivity along the same side of Zagros Mts is higher than for all other pairwise comparisons. The populations South of Zagros Mts, showed low differentiation between SW Zagros and SE Zagros and this was confirmed by the STRUCTURE results. Recent studies also indicate a plant species migration along altitudinal gradients [11]. A study of the population genetic diversity of Rhodiola dumulosa (Crassulaceae) between rivers in Mountains ridges also showed strong differentiation [5]. From that PCoA of individuals, there is a clear geographical gradient along the first axis and corresponding to their position alongside or across mountains.

Microsatellites with dispersal evidence from maternal cpDNA of trnH-psbA
The most common haplotype (cp 378) can be observed in most regions but cp 379 in only two regions (NW Alborz;up of Lake Urmieh and N.Zagros) despite their large allelic differentiation of microsats. Thus the mutation of the chloroplast mononucleotide repeat (1 base only) is most likely independent as event in both regions and not an indication of long distance dispersal. Another shorter variant (cp377) also occurs in only two regions (NE of Alborz and N.Zagros but remain very different for their microsats. So again, it most likely is an independent mutation event in both regions and not necessarily an indication of long distance dispersal. However, the haplotype variant cp380 and combined with a unique nuclear ITS-B variant could be an indication of historical dispersal over longer distances. These slow evolving genes are similar but the faster evolving microsats show large differences between pop 1, 5 and 15/35. Therefore, we assume that the isolation of 15/35 is not a recent one. Potamogeton represents the largest genus in Potamogetonaceae, including about 100 species and 50 interspecific hybrids worldwide. Hybridization therefore can be another complexity that is thought to be relatively frequent, because complex series of polyploidy and aneuploidy are present within the genus [28]. The hexaploid P. pectinatus having 78 chromosomes could form hybrids with other linear leafed taxa such as P. filiformis and P. amblyphyllus [72]. This issue might have an effect on the genetic diversity of Potamogeton populations. In our study we have detected only P. pectinatus haplotypes in correlation with specific microsatellite alleles for each site. No chloroplast capture was found. Consequently, we assume the absence of hybrid specimen in our samples.

Habitat
It is proven that habitat diversity often makes ecological barriers against gene flow that leads to interpopulation genetic separation [73,74]. Different populations produce different selection pressures in the presence of environmental conditions [75,73]. Iran has different ecological conditions that lead to different plant vegetations and biogeographic regions. In recent years, changing environmental conditions such as global warming, drought or in some cases grazing by fishes could be interesting for knowing how genetic diversity is retained in such situations [74]. Iran lies in a region that confronts droughts. This phenomenon could have led to the bottleneck effect in Vimcheh and Izeh because of temporary habitats at low level of river water. In Azbaran, there is a high level of fish grazing pressure and also probably water bird feed on seeds and tubers that might have led to repeated bottlenecks ( Table 3). The highest allelic and heterozygosity level was observed in N.Alborz and Kaniborazan population in West of Alborz-Azarbayejan Province. These provinces is documented to harbor the largest ecological diversity of Iran. Assignment of individuals to the own Kaniborazan population was highest because this population had the highest number of private alleles. In aquatic ecosystems, flow regimes and seasonal changing of water level lead to spread the seed to different habitats, inducing genetic diversity of aquatic populations [40]. Our results showed few shared alleles giving a stronger population genetic structure. The flora of the mountains provides many examples of closely related species or subspecies with populations that show limited gene flow due to these inherent geographic barriers [7,8]. In this study the individuals corresponding to populations of Aligoodarz, Borujen, Yaschaman and Googhar were separated from other P. pectinatus populations with 100% bootstrap. Additional evidence from ITS and cpDNA was provided and indicated these as an evolutionary significant unit (ESU). The morphological characters such as leaf sheath, width of leaf and tip of leaf were different from all other populations. These populations clearly belong to an evolutionary lineage within P. pectinatus and do not represent another species (e.g. P. amblyphyllus) because their ITS in a NJtree was separated (97% bootstrap) from P. amblyphyllus (S. amblyophylla) of NCBI (next paper under submission).
As a conclusion, when considering P. pectinatus populations over a long gradient of regions across mountain ranges and around deserts in Iran, their genetic differentiation over distances was more pronounced and structured than for temperate regions. However, across such habitats P. pectinatus maintained an even higher level of allele and gene diversity than in temperate regions.