| Table 1. The geographical characteristics of four populations of C. sativa in Hyrcanian forest and the sample size in each population. | ||||
| Population name | Latitude, Longitude | Altitude (m a.s.l.) | Sample size (No.) | Population size (Hectare) |
| Veysroud (V) | 37°16´N, 49°15´E | 211–711 | 8 | 60 |
| Shafaroud (SH) | 37°30´N, 49°02´E | 200–360 | 8 | 50 |
| Siyahmazgi (S) | 40°97´N, 35°49´E | 290–350 | 8 | 40 |
| QaleRoudkhan (R) | 37°05´N, 49°14´E | 200–400 | 8 | 50 |

Fig. 1. The geographical positions of Hyrcanian populations of C. sativa in North of Iran. R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region.
| Table 2. Sequences of Castanea SSR primer pairs, PCR annealing temperature and PCR expected product size. | |||||||
| Marker | Number | SSR locus | Primer sequences (5’ → 3’) | Ta (°C) | Size of alleles (bp) | Reference | |
| NU SSR | 1 | CsCAT 1 | F | GAGAATGCCCACTTTTGCA | 54 | 177–212 | Marinoni et al. 2003 |
| R | GCTCCCTTATGGTCTCG | ||||||
| 2 | CsCAT 3 | F | CACTATTTTATCATGGACGG | 58 | 179–198 | ||
| R | CGAATTGAGAGTTCATACTC | ||||||
| 3 | CsCAT 6 | F | AGTGCTCGTGGTCAGTGAG | 54 | 199–238 | ||
| R | CAACTCTGCATGATAAC | ||||||
| 4 | CsCAT 7 | F | GAACATGATGATTGGCCTC | 54 | 179–223 | ||
| R | CCAAACATGACATATGTCCC | ||||||
| 5 | CsCAT 14 | F | CGAGGTTGTTGTTCATCATTAC | 54 | 133–151 | ||
| R | GATCTCAAGTCAAAAGGTGTC | ||||||
| 6 | CsCAT 15 | F | TTCTGCGACCTCGAAACCGA | 56 | 123–155 | ||
| R | GCTAGGGTTTTCATTTCTAG | ||||||
| 7 | CsCAT 17 | F | TTGGCTATACTTGTTCTGCAAG | 56 | 132–151 | ||
| R | GCCCCATGTTTTCTTCCATGG | ||||||
| 8 | CsCAT 41 | F | AAGTCAGCAACACCATATGC | 58 | 199–238 | ||
| R | CCCACTGTTCATGAGTTTCT | ||||||
| 9 | EMC 38 | F | TTTCCCTATTTCTAGTTTGTGATG | 58 | 227–261 | Buck et al. 2003 | |
| R | ATGGCGCTTTGGATGAAC | ||||||
| 10 | QpZAG110 | F | GGAGGCTTCCTTCAACCTACT | 54 | 207–223 | Steinkellner et al. 1997 | |
| R | GATCTCTTGTGTGCTGTATTT | ||||||
| CP SSR | 1 | CCMP2 | F | GATCCCGGACGTAATCCTG | 55 | 234 | Weising and Gardner 1999 |
| R | ATCGTACCGAGGGTTCGAAT | ||||||
| 2 | CCMP4 | F | AATGCTGAATCGAYGACCTA | 55 | 114 | ||
| R | CCAAAATATTBGGAGGACTCT | ||||||
| 3 | CCMP10 | F | TTTTTTTTTAGTGAACGTGTCA | 55 | 113 | ||
| R | TTCGTCGDCGTAGTAAATAG | ||||||
| 4 | CMCS2 | F | GAGCCATTCCCTTTTAGAAT | 55 | 141 | Sebastiani et al. 2004 | |
| R | TTGAAAACCGGTATAGTTCG | ||||||
| 5 | CMCS7 | F | AAGCGAGATGAATGAGTTTT | 55 | 205 | ||
| R | AAAATTGGATTGATTATTGACT | ||||||
| 6 | CMCS8 | F | GGTCTATTTTTCCACTCACAA | 55 | 179 | ||
| R | AGAAATAAACACCCCCATTA | ||||||
| 7 | CMCS12 | F | ATATTGGTAAAACGGCAACT | 55 | 216 | ||
| R | TTTATGGCATGAAACAACTC | ||||||
| 8 | CMCS14 | F | GGATTGTAACAAATTTTTCAGG | 55 | 176 | ||
| R | GTGCAAGGAATGTCGAACTA | ||||||
| F = forward primer; R = reverse primer; SSR = Simple sequence repeat; CP SSR = Chloroplast SSR; NU SSR = Nuclear SSR; Ta = PCR annealing temperature | |||||||
| Table 3. Genetic variability within four C. sativa populations based on SSR markers. | ||||||||||
| Marker | Pop | N | Na | Ne | I | Ho | He | UHe | F | |
| NU SSR | R | Mean | 8.000 | 3.000 | 2.142 | 0.807 | 0.538 | 0.463 | 0.494 | –0.147 |
| SE | 0.000 | 0.365 | 0.243 | 0.133 | 0.110 | 0.074 | 0.078 | 0.125 | ||
| S | Mean | 8.000 | 2.400 | 1.903 | 0.679 | 0.563 | 0.434 | 0.463 | –0.293 | |
| SE | 0.000 | 0.221 | 0.148 | 0.093 | 0.112 | 0.059 | 0.063 | 0.166 | ||
| SH | Mean | 8.000 | 3.300 | 2.435 | 0.962 | 0.675 | 0.559 | 0.597 | –0.230 | |
| SE | 0.000 | 0.213 | 0.197 | 0.081 | 0.075 | 0.044 | 0.046 | 0.130 | ||
| V | Mean | 8.000 | 2.700 | 2.174 | 0.781 | 0.475 | 0.477 | 0.508 | –0.007 | |
| SE | 0.000 | 0.260 | 0.223 | 0.122 | 0.096 | 0.073 | 0.078 | 0.115 | ||
| Total | Mean | 8.000 | 2.85 | 2.163 | 0.630 | 0.563 | 0.512 | 0.516 | –0.169 | |
| SE | 0.000 | 0.265 | 0.203 | 0.107 | 0.098 | 0.063 | 0.071 | 0.134 | ||
| NU SSR = Nuclear SSR; R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region; N = number of individuals analyzed; Na = number of alleles; Ne = effective number of alleles; I = Shannon index; Ho = observed heterozygosity; He = expected heterozygosity; UHe = unbiased expected heterozygosity; F = F-values; SE = Standard errors | ||||||||||
| Table 4. Analyses of molecular variance (AMOVA) for four populations of C. sativa from Hyrcanian forest by Nuclear SSR (NU SSR). Statistics include sums of squared deviations (SS); mean squared deviations (MS), variance component estimates (Est. Var.), the percentage of the total variance contributed by each component; estimator of relative genetic differentiation based on fraction of total variance of allele size between two subpopulations (Rst) and the probability of obtaining a more extreme component estimate by chance alone. | ||||||
| Marker | S.O.V. | D.f | SS | MS | Est. Var. | % of Variance |
| NU SSR | Between groups | 3 | 8.305.094 | 2.768.365 | 130.923 | 16% |
| Within groups | 60 | 40.415.375 | 673.590 | 673.590 | 84% | |
| Total Sum | 63 | 48.720.469 | 804.513 | 100% | ||
| Stat | Value(Φ) | P(rand > = data) | ||||
| Rst | 0.163 | 0.010 | ||||
| Table 5. Pairwise estimated of Nei’s genetic distance and the calculated Fst based on 10 nuclear SSR markers among four populations of C. sativa in Hyrcanian forest. | ||||
| Population | R | S | SH | |
| Nei distance | R | 0.000 | ||
| S | 0.215 | |||
| SH | 0.237 | 0.175 | 0.000 | |
| V | 0.012 | 0.095 | 0.208 | |
| Fst | R | 0.000 | ||
| S | 0.132 | 0.000 | ||
| SH | 0.185 | 0.148 | 0.000 | |
| V | 0.034 | 0.114 | 0.163 | |
| NU SSR = Nuclear SSR; R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region; Fst = Fixation index | ||||

Fig. 2. Population structure of C. sativa was estimated based on Structure analysis. In the figure, the individuals of Iranian Castanea were sorted by Bayesian clustering approaches based on 10 nuclear microsatellite loci. Each bar represents a single individual analyzed. Legend of number in brackets: 1 = R1; 2 = R2; 3 = R3; 4 = R4; 5 = R5; 6 = R6; 7 = R7; 8 = R8; 9 = S1; 10 = S2; 11 = S3; 12 = S4; 13 = S5; 14 = S6; 15 = S7; 16 = S8; 17 = SH1; 18 = SH2; 19 = SH3; 20 = SH4; 21 = SH5; 22 = SH7; 23 = SH10; 24 = SH12; 25 = V1; 26 = V2; 27 = V3; 28 = V4; 29 = V5; 30 = V6; 31 = V7; 32 = V8.
| Table 6. The proportion membership of each individual of the four populations in each of the 3 gene pools identified by Structure analysis. | ||||||||||
| Sample | Structuregene pools (% inferred ancestry) | Sample | Structuregene pools (% inferred ancestry) | |||||||
| red | green | blue | red | green | blue | |||||
| 1 | R1 | 0.972 | 17 | SH1 | 0.950 | |||||
| 2 | R2 | 0.750 | 18 | SH2 | 0.930 | |||||
| 3 | R3 | 0.921 | 19 | SH3 | 0.972 | |||||
| 4 | R4 | 0.927 | 20 | SH4 | 0.976 | |||||
| 5 | R5 | 0.951 | 21 | SH5 | 0.974 | |||||
| 6 | R6 | 0.977 | 22 | SH7 | 0.967 | |||||
| 7 | R7 | 0.896 | 23 | SH10 | 0.943 | |||||
| 8 | R8 | 0.973 | 24 | SH12 | 0.803 | |||||
| 9 | S1 | 0.965 | 25 | V1 | 0.747 | |||||
| 10 | S2 | 0.784 | 26 | V2 | 0.649 | |||||
| 11 | S3 | 0.781 | 27 | V3 | 0.911 | |||||
| 12 | S4 | 0.978 | 28 | V4 | 0.956 | |||||
| 13 | S5 | 0.912 | 28 | V5 | 0.910 | |||||
| 14 | S6 | 0.973 | 30 | V6 | 0.727 | |||||
| 15 | S7 | 0.952 | 31 | V7 | 0.820 | |||||
| 16 | S8 | 0.744 | 32 | V8 | 0.964 | |||||
| R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region | ||||||||||

Fig. 3. Assignment percentage of individuals based on structure analysis per each population (K = 3). R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region.

Fig. 4. Dendrogram depicting the distribution of genotypes of the four populations of C. sativa in Hyrcanian forest based on UMGMA method. R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region.
| Table 7. Results of BOTTLENECK tests elaborated on the four C. sativa populations analysed at 10 polymorphicnuclear microsatellite loci. | ||||||
| Test | Mutation Model | Population | ||||
| V | SH | S | R | |||
| Sign | P-value | IAM | 0.08 | 0.225 | 0.272 | 0.347 |
| SMM | 0.034* | 0.807 | 0.513 | 0.359 | ||
| TPM | 0.072 | 0.359 | 0.335 | 0.369 | ||
| Wilcoxon | P-value | IAM | 0.359 | 0.921 | 0.496 | 1.000 |
| SMM | 0.203 | 0.625 | 0.734 | 0.82 | ||
| TPM | 0.250 | 0.625 | 0.57 | 0.91 | ||
| Standardized differences | P-value | IAM | <0.001* | 0.468 | 0.003* | 0.396 |
| SMM | <0.001* | 0.272 | 0.04* | 0.057 | ||
| TPM | <0.001* | 0.114 | 0.012 | 0.208 | ||
| * Significant at level 0.05; P = Probability; IAM = Infinite allele model, TPM = Two phase model, SMM = Stepwise mutation model; R = QalehRoudkhan region; S = SiyahMazgi region; SH = Shafaroud region; V = Veysroud region | ||||||