青藏高原东缘青杨（Populus cathayana Rehd.）遗传多样性研究

青杨（Populus cathayana Rehd.）是青杨派杨树的主要树种之一，为我国特有乡土树种，其主要分布区之一是我国的青藏高原，集中分布地带在甘肃省中部及青海省东部，四川省西北部岷江上游和松潘等地区。本研究以青藏高原东缘青杨天然分布区的6个群体143个个体为材料，用AFLP、SSR和叶绿体SSR分子标记分析青杨天然群体的遗传多样性，分析其遗传结构和分化，比较6个群体间遗传多样性的高低和群体间的遗传关系。旨在为青杨基因资源评价、保护与保存、遗传改良策略制定等提供科学理论依据。通过以上研究，得出如下主要研究结果： 1 AFLP分子标记研究结果 采用4对选择性引物对6个青杨天然群体143个个体进行分析，扩增谱带分析共检测到175个位点，其中173个位点表现为多态，多态位点百分率高达98.9%。从整体上表现出较高的遗传多样性，Nei’s基因多样度(h)水平为0.306。从青杨天然群体位点分布来看，有高达20%的位点（32位点）为群体所特有，仅有9.14%的位点（16位点）在所有群体中存在。群体间的遗传分化极大，所有遗传变异中，有48.9%的遗传变异存在于群体间。在个体群丛（Individuals cluster）和主坐标（PCO analysis）分析中，青杨各群体未呈现任何地理模式，Mantel检测也显示各群体间遗传距离与地理距离无明显相关。研究认为，由于地理和空间上大尺度的隔离和地形地貌复杂使得群体间无法进行基因交流，导致群体间遗传分化极大，另外各群体在不同的选择压力下，经历各自独立的进化历程，这些都可能导致群体间遗传距离与地理距离的不相关。 2 SSR分子标记研究结果 在SSR分析中，7个位点在6个青杨天然群体143个个体中共检测到79个等位基因，每位点检测到的等位基因数在5-16之间，平均11.3个，总体上多态位点百分率达100%。平均观察杂合度和期望杂合度分别为0.792和0.802。Hardy-Weinberg平衡检验表明青杨大部分群体都处于非平衡状态，群体大部分位点都是偏离哈迪-温伯格平衡（76.3%），只有23.7%的测验满足哈迪-温伯格平衡。分析青杨天然群体内和群体间的遗传变异，基因分化系数（GST）为0.373，即有62.7%的遗传变异存在群体内，37.3%的遗传变异存在群体间。群体内的遗传变异高于群体间水平。根据各群体遗传距离UPGMA聚类分析，有来自相临分布区、近似气候类型的群体聚在一起的趋势，但Mantel检测反映遗传距离与地理距离间并无明显相关性。 3 cpSSR分子标记研究结果 分析来自青藏高原东缘6个青杨天然群体，所用cpSSR引物中有5对cpSSR引物（CCMP2、CCMP5、SCUO01、SCU03、SCU07）都表现较高的多态性，单个引物检测的片段数都在4以上。5对cpSSR引物共检测片段数26个，组成了12种叶绿体DNA单倍型。各群体的单倍型分布和频率有较大差异，群体单倍型多样性范围为0-0.4926，TS、JZ、PW和SHY群体单倍型多样性高于QHY和LED群体水平。本研究发现，分布在青藏高原东缘的青杨天然群体，群体间不存在共享的单倍型，各群体间存在极大的遗传分化（GST=0.9223）。从青藏高原东缘地区经历的地质历史事件来看，第四纪的冰期气候变迁可能是造成青杨现今遗传结构模式的主要因素之一。根据单倍型在各群体的分布情况，进行青杨群体聚类分析结果，各群体无明显的分组现象，青杨各群体也未呈现任何清晰地理模式。 由于不同分子标记在对群体遗传多样性检测能力与效率上存在差异，所以三种标记检测的青藏高原东缘青杨天然群体遗传多性水平也不尽一致，但在与用同种方法检测其它物种或同一物种不同种源群体比较，三种分子标记方法都揭示了青藏高原东缘青杨天然群体具有中等偏上的遗传多样性水平。结果分析表明，群体间遗传分化极大，这是由于青杨天然群体分布于青藏高原东缘，既有高原又有高山峡谷，由于地理和空间上大尺度的隔离和地形地貌复杂导致了基因流物理上的阻隔。三种分子标记研究结果经Mantel分析检测，遗传距离与地理距离之间都无明显相关性。较为一致的解释是，青杨分布区域地理和空间上大尺度的隔离和和地形地貌复杂导致群体之间不存在均匀扩散现象，另外各群体在不同的选择压力下，经历各自独立的进化历程，这些都可能导致群体间遗传距离与地理距离的不相关。 The wide geographical and climatic distribution of P. cathayana Rehd. indicates that there is a large amount of genetic diversity available, which can be exploited for conservation, breeding programs and afforestation schemes. The results are as follows: 1 Research results of AFLP genetic diversity In present study, genetic diversity was evaluated in the natural populations of P. cathayana originating from southern and eastern edge of the Qinghai-Tibetan Plateau of China by means of AFLP markers. For four primer combinations, a total of 175 bands were obtained, of which 173 (98.9%) were polymorphic. Six natural populations of P. cathayana possessed different levels of genetic diversity, high level of genetic differentiation existed among populations (GST=0.489) of P. cathayana. Individuals cluster and PCO analysis based on Jaccard’s similarity coefficient also showed evident population genetic structure with high level population genetic differentiation. The long evolutionary process coupled with genetic drift within populations, rather than contemporary gene flow, are the major forces shaping genetic structure of P. cathayana populations. Moreover, there is no correspondence between geographical and genetic distances in the populations of P. cathayana, seldom gene exchange among populations and different selection pressures may be the causes. Our finding of different levels of genetic diversity within population and high level of genetic differentiation among populations provided promising condition for further breeding or conservation programs. 2 Research results of SSR genetic diversity In this study, the genetic diversity of P. cathayana was investigated using microsatellite markers. In a total of 150 individuals collected from six natural populations in the southeastern part of the Qinghai-Tibetan Plateau in China, a high level of microsatellite polymorphism was detected. At the seven investigated microsatellite loci, the number of alleles per locus ranged from 5 to 16, with a mean of 11.3, the observed heterozygosities across populations ranged from 0.408 to 0.986, with a mean of 0.792, and the expected heterozygosities across populations ranged from 0.511 to 0.891, with a mean of 0.802. The proportion of genetic differentiation among populations accounted for 37.3% of the whole genetic diversity. The presence of such a high level of genetic diversity could be attributed to the features of the species and the habitats where the sampled populations occur: The southeastern part of the Qinghai-Tibetan Plateau is regarded as the natural distribution and variation center of the genus Populus in China. Variation in environmental conditions and selection pressures in different populations, and topographic dispersal barriers could be factors associated with the high level of genetic differentiation found among populations. The populations possessed significant heterozygosity excesses, which may be due to extensive population mixing at the local scale. The cluster analysis showed that the populations are not strictly grouped according to their geographic distances but the habitat characteristics also influence the divergence pattern. In addition, we suggest that population SHY should be regarded as an ecologically divergent species of P. cathayana. 3 Research results of cpSSR genetic diversity Genetic diversity of six natural populations of P. cathayana originating from the southeastern part of the Qinghai-Tibetan Plateau in China was studied by use of cpSSR markers. Based on 5 pairs of polymorphic primers screened from 12 pairs of primers, twenty-six different length fragments and twelve different kinds of haplotypes were reduced in 143 samples. There were significant variant haplotypes among the populations．There were no shared haplotypes found among populations, analysis of molecular variance indicated that a high proportion of the total genetic variance was attributable to variations among populations (92.23%). The pattern of genetic structure which is associated with spatial separation, variation in environmental conditions and selection pressures in different populations, is also the result of geological historical factor. A molecular phylogenetic tree based on the 12 haplotypes showed that the populations are not strictly grouped according to their geographic distances.

川西青杨组(Section Tacamahaca Spach)不同种的亲缘关系分析

青杨组（Section Tacamahaca Spach）杨树是我国重要的乡土经济树种，目前对其分子遗传变异和系统进化的研究还很少，尤其是在青杨组杨树遗传资源极为丰富的川西地区，杨树的分子进化及亲缘关系的研究极为缺乏，非常不利于该树种遗传资源的开发和利用。本研究从川西地区收集了青杨（Populus cathayana）、青海杨（P. prezewalskii）、滇杨（P. yunnanensis）、康定杨（P. kangdingensis）、西南杨（P. schneideri）、小叶杨（P.simonii）和三脉青杨（P. trinervis）这7 个青杨组树种的10 个群体，利用多种分子标记手段对其种间的亲缘关系进行比较，并结合形态和地史资料进行了全面的研究和评价，得到了如下的主要研究结果： 1. SSR 和ISSR 位点变异丰富。通过10 对引物对50 个杨树个体的DNA 样品进行了SSR 分析，所有位点展现了丰富的群体间和种间的多态性，多态位点率达到了100%，每位点的等位基因数变化范围为5 ~ 17，平均为11.9 个；通过11 条ISSR 随机引物对供试的混合DNA 样品进行分析，共检测到130 个标记，其中多态性标记为119 个，多态百分率为91.5%。研究认为，SSR 单个标记能展现高水平信息，而ISSR 单个引物能探测更多数量多态性。通过两个标记的遗传距离、聚类图和PCA 分析，表明：同一种内不同群体间的同源性最高；康定杨和西南杨有较近的亲缘关系；小叶杨和三脉青杨聚合在一起，显示了其相互较近的亲缘关系；滇杨与其它杨树种可能存在着较远的亲缘关系。 2. 采用4 对选择性引物对7 个青杨组杨树种10 个群体进行AFLP 分析，总共扩增出284 个位点，其中200 个位点显示出了多态性，多态位点百分比为70.4%，平均多态带为50 条。TE-AFLP 的分析总共扩增出192 个位点，其中139 个位点显示出了多态性，多态位点百分比为72.4%，平均多态带为34.7 条。比较的结果表明AFLP、TE-AFLP 的遗传信息含量比较接近，略小于ISSR，大约仅为SSR 的1/3；但这两个基于AFLP 的标记系统的信息探察能力也远大于ISSR 和SSR 标记系统。这两个分子标记的聚类结果，显示小叶杨、三脉青杨和滇杨三个种聚为一组，其中小叶杨与三脉青杨的亲缘关系更近；其它几个杨树种聚为一类，西南杨与青杨表现出较近的亲缘关系。 3. 所有7 对cpSSR 引物中，仅有4 个叶绿体位点在种间具有多态性，而在种内群体中并不具有多态性，共检测出13 个条带，组合成了4 种不同的单倍型；对于cpDNA的5 对引物，共检测出了73 条酶切片段，其中52 条是多态带，组合成了9 种不同的单倍型；而5 对mtDNA 通用引物未能检测出多态性的条带，表现出线粒体的保守性。叶绿体的聚类分析认为，小叶杨、三脉青杨和滇杨有较近的母性起源，且依次聚合；其余四种杨树聚为一类，并且康定杨与西南杨表现出最近的亲缘关系，并依次与青杨和青海杨聚合。 4. 根据本文的分子数据，结合形态和生境分布资料分析认为：青杨组杨树种内群体间的遗传变异程度是小于种间的遗传差异，显示了与传统分类一致的结果；三脉青杨和小叶杨有很近的亲缘关系，可能拥有相同的祖先类群；滇杨与小叶杨和三脉青杨之间具有一定的亲缘关系，特别是在其母性祖先的起源上有着一定的同源性；西南杨与青杨和康定杨均保持着较近的亲缘关系，且有可能是这两个种原始祖先杂交后所形成的。 Although western Sichuan is regarded as a natural distribution and variation center forthe Section Tacahamaca of the Populus species in China, little is currently known about themajority of poplar species occurring in this region. In the present study, molecular data wereutilized to determine the genetic relationships among Populus species in Section Tacamahacain western Sichuan including P. cathayana, P. prezewalskii, P. yunnanensis, P. kangdingensis,P. schneideri, P. simonii and P.trinervis. The results are as fellows: 1. The genetic variation at SSR and ISSR loci was abundant. All the 10 SSR loci werepolymorphic, and the number of alleles per locus varied from 5 to 17 with a mean valueequaling 11.9. Based on the 11 ISSR primers, 130 clear and reproducible DNA fragmentswere generated, of which 119 (91.5%) were polymorphic. Our results reveal that single SSRlocus can present more genetic information, while more polymorphic bands can be detectedby single ISSR primer. Moreover, the genetic distance, cluster and PCA analysisdemonstrated that: a close relationship among accessions of the same species and suggestedmonophyly in P. przewalskii and P. cathayana; P. schneideri is genetically highly similar to P.kangdingensis; P. trinervis and P. simonii have a close genetic affinity; P. yunnanensis isdistinct from the other species. 2. Genetic relationships of poplar species in Section Tacamahaca from western Sichuanwere evaluated by means of AFLP and TE-AFLP. For four AFLP primer combinations, atotal of 284 bands were obtained of which 200 (70.4%) were polymorphic with the average of50 polymorphic bands. For four TE-AFLP primer combinations, a total of 192 band wereobtained of which 139 (72.4%) were polymorphic with the average of 34.7 polymorphicbands. Our results indicate that the genetic information of AFLP is similar to that ofTE-AFLP, and little less than that of ISSR, but only about 1/3 of that of SSR. However, theability of information detection of the two AFLP-based markers is much higher than that ofISSR and SSR. In addition, the cluster analysis of AFLP, TE-AFLP and combined data revealthat: P. yunnanensis, P. trinervis and P. simonii clustered together, and P. trinervis and P.simonii showed more closed affinity; the other four poplar species clustered together, P.cathayana and P. schneideri showed more closed origin especially. 3. The cpSSR analysis for seven Populus species belonging to the Section Tacamahaca.Four out of the seven analyzed chloroplast loci were polymorphic, whereas none of the lociwere polymorphic across the accessions within a species. 13 bands and 4 different kinds ofhaplotypes were reduced. Based on 5 pairs of cpDNA primers, 73 fragments (52 polymorphic)and 9 kinds of haplotypes were produced. However, none of the polymorphic was detected bythe 5 mtDNA primer pairs, revealing conservation of mtDNA region. The cluster analysis ofcpDNA revealed that: similar maternal phylogeny among P. yunnanensis, P. trinervis and P.simonii; the other four species clustered together, P. schneideri and P. kangdingensis showedmore closed maternal lineage especially. 4. Our molecular data, morphological characters and nature habitat revealed that: sameto the traditional taxonomy assignment, genetic variation within a same Populus species islower than that among Populus species in Section Tacamahaca; P. yunnanensis may share itschloroplast ancestor with P. trinervis and P. simonii; moreover, sister genetic relationship of P.trinervis and P. simonii indicated their similar origin; P. schneideri clustered with P.kangdingensis and P. cathayana, respectively, and may have derived from an ancienthybridization event involving the ancestors of the two species.

中国沙棘(Hippophae rhamnoides subsp. sinensis)天然群体的遗传变异分析

近十年，植物群体遗传学的研究飞速发展，然而与海拔相关的植物群体遗传结构和遗传变异研究却相对较少。到目前为止，还不清楚遗传变异与海拔之间是否有一个通用的格局。在山区，各种生态因子，如温度、降水、降雪、紫外线辐射强度以及土壤成分都随海拔梯度急剧变化，造成了即使在一个小的空间区域，植被类型变化显著，这种高山环境的异质性和复杂性为我们研究植物群体遗传结构和分化提供了方便。沙棘(Hippophea)属于胡颓子科(Elaeagnaceae)为多年生落叶灌木或乔木，雌雄异株，天然种群分布极为广泛。中国沙棘(H. rhamnoides subsp. sinensis)是沙棘属植物中分布较广的一个亚种，种内形态变异非常丰富，加之其具有独特的繁育系统和广泛的生态地理分布，是研究沙棘属植物遗传变异和系统分化的理想材料。本文从1,800 m 到3,400 m 分5 个海拔梯度进行取样，用RAPD 和cpSSR 分子标记研究了卧龙自然保护区中国沙棘天然群体的遗传结构和遗传变异。5 个取样群体依次标记为A、B、C、D 和E，它们分别代表分布在海拔1,800，2,200，2,600，3,000 和3,400 m 的5 个天然群体。RAPD实验用11 条寡核苷酸引物，扩增得到151 个重复性好的位点，其中143 个多态位点，多态率达94.7%。在5 个沙棘群体中，总遗传多样性值(HT)为0.289，B群体内的遗传多样性值为0.315，这完全符合沙棘这种多年生、远交的木本植物具有高遗传变异的特性。5 个群体内遗传多样性随海拔升高呈低－高－低变异趋势，在2,200 m海拔处的B群体遗传多样性达最大值0.315，3,400 m海拔处的E群体则表现最小仅0.098。5 个群体间的遗传分化值GST＝0.406，也即是说有40.6%的遗传变异存在于群体间，1,800 m海拔处的A群体与其它群体的明显分离是造成群体间遗传分化大的原因。UPGMA聚类图和PCoA散点图进一步确证了5 个群体间的关系和所有个体间的关系。最后，经过Mantel检测，遗传距离与海拔表现了明显的相关性(r = 0.646, P = 0.011)。cpSSR 实验中，经过对24 对cpSSR 通用引物筛选，11 对引物能扩增出特异性条带，只有2 对引物(ccmp2 和ARCP4)呈现多态性。4 个等位基因共组合出4 种单倍型，单倍型Ⅰ出现在A 群体的所有个体和B 群体的8 个个体中，C、D、E 三个群体均不含有，而单倍型Ⅱ出现在C、D、E 三个群体的所有个体及B 群体的18 个个体中，A 群体不含有。另外两种单倍型Ⅲ和Ⅳ为稀有类型，仅B 群体中的4 个个体拥有。这种单倍型分布模式和TFPGA 群体聚类图揭示了，C、D、E 群体可能来源于同一祖先种，而A 群体却是由另一祖先种发展起来的，B 群体则兼具了这两种起源种的信息，这可能是因为在历史上的某一时期，在中国沙棘群体高山分化的过程中，B 群体处某个或者某些个体发生了基因突变，具备了适应高海拔环境的能力，产生了高海拔沙棘群体的祖先种。 In recent ten years, studies about population genetics of plants developed rapidly,whereas their genetic structure and genetic variation along altitudinal gradients have beenstudied relatively little. So far, it is uncleared whether there is a common pattern betweengenetic variation and altitudinal gradients. In the mountain environments, importantecological factors, e.g., temperature, rainfall, snowfall, ultraviolet radiation and soil substratesetc., change rapidly with altitudes, which cause the vegetation distribution varying typically,even on a small spatial scale. The mountain environments, which are heterogeneous andcomplex, facilitate and offer a good opportunity to characterize population genetic structureand population differentiation.The species of the genus Hippophae L. (Elaeagnaceae) are perennial deciduous shrubs ortrees, which are dioecious, wind-pollinated pioneer plants. The natural genus has a widedistribution extending from Northern Europe through Central Europe and Central Asia toChina. According to the latest taxonomy, the genus Hippophae is divided into six species and12 subspecies. The subspecies H. rhamnoides ssp. sinensis shows significant morphologicalvariations, large geographic range and dominantly outcrossing mating system. Thesecharacteristics of the subspecies are favourable to elucidate genetic variation and systemevolution. To estimate genetic variation and genetic structure of H. rhamnoides ssp. sinensisat different altitudes, we surveyed five natural populations in the Wolong Natural Reserve at altitudes ranging from 1,800 to 3,400 m above sea level (a.s.l.) using random amplifiedpolymorphic DNA markers (RAPDs) and cpSSR molecular methods. The five populations A,B, C, D, and E correspond to the altitudes 1,800, 2,200, 2,600, 3,000 and 3,400 m,respectively.Based on 11 decamer primers, a total of 151 reproducible DNA loci were yielded, ofwhich 143 were polymorphic and the percentage of polymorphic loci equaled 94.7%. Amongthe five populations investigated, the total gene diversity (HT) and gene diversity within population B equaled 0.289 and 0.315, respectively, which are modest for a subspecies of H.rhamnoides, which is an outcrossing, long-lived, woody plant. The amount of geneticvariation within populations varied from 0.098 within population E (3,400 m a.s.l.) to 0.315within population B (2,200 m a.s.l.). The coefficient of gene differentiation (GST) amongpopulations equaled 0.406 and revealed that 40.6% of the genetic variance existed amongpopulations and 59.4% within populations. The population A (1,800 m a.s.l.) differed greatlyfrom the other four populations, which contributes to high genetic differentiation. A UPGMAcluster analysis and principal coordinate analyses based on Nei's genetic distances furthercorroborated the relationships among the five populations and all the sampling individuals,respectively. Mantel tests detected a significant correlation between genetic distances andaltitudinal gradients (r = 0.646, P = 0.011).Eleven of the original 24 cpSSR primer pairs tested produced good PCR products, onlytwo (ccmp2 and ARCP4) of which were polymorphic. Four total length variants (alleles) werecombined resulting in 4 haplotypes. The haplotype was present in all individuals of Ⅰpopulation A and 8 individuals of populations B, the other three populations (C, D and Epopulations) did not share. The haplotype was present in all individuals of populations C, D Ⅱand E and 18 individuals of populations B, population A did not share. The other twohaplotypes and were rare haplotypes, which were only shared in 4 individuals of Ⅲ Ⅳpopulation B. The distribution of haplotypes and TFPGA population clustering map showedthat the populations C, D and E might be origined from one ancestor seed and population Amight be from another, whereas population B owned information of the two ancestor seeds. Itwas because that gene mutation within some individual or seed in the location of population Bwas likely to happen in the history of H. rhamnoides, which was the original ancestor of thehigh-altitude populations.

Chloroplast microsatellite markers for the Neotropical orchid genus Epidendrum, and cross-amplification in other Laeliinae species (Orchidaceae)

One of the most significant challenges confronting orchid researchers is the lack of specific molecular markers, mainly for species in the Neotropics. Here we report the first set of specific chloroplast microsatellite primers (cpSSR) developed for Neotropical orchids. In total, nine polymorphic cpSSR loci were isolated and characterized in four species occurring in the Brazilian Atlantic Rainforest: Epidendrum cinnabarinum, E. denticulatum, E. fulgens and E. puniceoluteum. Levels of intraspecific polymorphism were characterized using two populations for each species, with 13-20 individuals each. Allele numbers varied from two to three per locus, while the number of haplotypes ranged from three to six per species. Extensive differentiation among the taxa was detected. All markers were successfully cross-amplified in eight other different genera. These cpSSRs markers will enable novel insights into the evolution of this important Neotropical genus.

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines.

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.

Current range characteristics of Swiss stone pine (Pinus cembra L.) along the Carpathians revealed by chloroplast SSR markers.

We investigated the diversity pattern of nine Swiss stone pine (Pinus cembra L.) populations along the Carpathian range including the High Tatras, by using six chloroplast DNA microsatellites (cpSSR). Our aim was to detect genetically distinct regions by clustering of populations, and to tackle possible historical colonization routes. Our analysis referred to an investigated geographical range with the two most distant populations situated at about 500 air km. We found that the most diverse populations are situated at the two edges of the investigated part, in the Retezat Mts. (South Carpathians) and the High Tatras, and diversity decreases towards the populations of the Eastern Carpathians. Hierarchical clustering and NMDS revealed that the populations of the South Carpathians with the Tatras form a distinct cluster, significantly separated from those of the Eastern Carpathians. Moreover, based on the most variable chloroplast microsatellites, the four populations of the two range edges are not significantly different. Our results, supported also by palynological and late glacial macrofossil evidences, indicate refugial territories within the Retezat Mts. that conserved rich haplotype composition. From this refugial territory Pinus cembra might have colonized the Eastern Carpathians, and this was accompanied by a gradual decrease in population diversity. Populations of the High Tatras might have had the same role in the colonizing events of the Carpathians, as positive correlation was detected among populations lying from each other at a distance of 280 km, the maximum distance between neighbouring populations.