高山松的二倍体杂种起源-来自三个基因组的证据

利用松科植物特殊的遗传体系（叶绿体基因组一父系遗传、线粒体基因组—母系遗传、核基因组一双亲遗传），我们对高山松及其两个亲本种进行了广泛的群体取样，通过线粒体基因nadl、叶绿体基因rbcL和trnL-F基因间区以及低拷贝核基因4CL的序列分析或PCR-RFLP分析，为高山松同倍体杂种起源假说提供了翔实的遗传学证据，同时在个体水平上探讨了高山松不同群体的遗传组成、群体遗传结构、基因交流方向、群体建立过程以及杂种基因组的进化。具体结果如下： 1．细胞质基因组分析 1)线粒体基因nudl分析 本研究对油松、高山松和云南松的19个群体、295个个体的线粒体基因nadl的一个内含子进行了序列分析或PCR-RFLP分析，共检测到3种线粒体DNA单倍型-A、B和C。油松所有的取样群体仅含单倍型A;除BX群体外，所有的云南松群体仅含单倍型B; 10个高山松群体中，5个群体固定单倍型A，4个群体固定单倍型B，1个群体(ZD)分布有A和B两种单倍型。2)叶绿体rbcL基因分析 对同一组群体的rbcL基因进行序列分析或PCR-RFLP分析，共检测到两个变异位点和三种叶绿体单倍型（TT、TC和GC）。TT和GC分别是油松和云南松种特异性叶绿体单倍型，而在高山松群体里则三种单倍型均有分布，而且TC单倍型广泛地分布在7个杂种群体中，该单倍型很可能来源于点突变或第三个已灭绝的亲本。rbcL基因检测到的高山松群体分化系数很高(Gst=0.533)。 3)叶绿体trn L-F区序列分析 叶绿体trnL-F分子标记检测到的不同单倍型的差异主要是由引物“e”下游120碱基处一个多聚T结构的长度变异所致（叶绿体SSR位点）。10个高山松群体中共检测到5种叶绿体单倍型，其中两种主要的单倍型(9T和11T)分别为油松和云南松的种特异性单倍型，其他单倍型均为非典型单倍型。群体遗传结构分析表明：杂种群体表现最高的遗传多样性，而且trnL-F分析得到的高山松群体的分化系数也很高( Gst=0.443)。 总之，对高山松、油松和云南松的同一组群体取样进行的细胞质基因组分析表明：高山松群体分布有油松和云南松种特异性的线粒体和叶绿体单倍型，该细胞质DNA单倍型的地理分布为假说“高山松为油松和云南松的的二倍体杂种”提供了翔实的遗传学证据。油松和云南松在不同的杂种群体中分别做父本和母本，即两亲本在杂交过程中发生了双向基因交流。群体遗传结构分析发现高山松群体表现最高的遗传多样性，而且群体间的分化系数很高。不同的杂种群体在遗传组成上的差异表明他们经历过不同的建立和进化历史。从线粒体和叶绿体单倍型的地理分布可以看出杂种群体的建立曾经历强烈的奠基者效应和回交。青藏高原的隆升对高山松的起源、杂种群体的适应辐射以及保持产生了重要的影响。川西南和滇西北作为青藏高原的东边边界，很可能是当初云南松和油松分布的重叠区及杂交地带，即高山松的起源地。 2．核基因4CL分析 对高山松、油松和云南松的19个群体、32个个体的低拷贝核基因4CL进行了克隆及序列分析，获得的78条序列可分为两种类型（类型A和类型B）。这两种类型明显的差别是类型A相对于类型B在内含子区有- 20bp的缺失。以华山松的3条序列为外类群，对得到的78条序列进行基因谱系分析，发现所有的序列分成明显的两支，分别对应于类型A和类型B，而且每一支均包含三个种的部分序列，表明4CL基因在这三个种分化之前就已发生重复。另一个明显的特点是某个种的一条序列与另一个种的序列比其与同种的其他序列关系更近，可能因基因交流（杂交和渐渗）、非共祖、致同进化和重组等进化事件所致。三种松树中共检测到4CL基因序列的两种类型和六个亚类型，高山松群体中没有发现杂种独特的类型或亚类型。高山松和云南松共享三种序列亚类型以及最多的序列多态性，表明这两个种之间曾存在广泛的基因交流。

独花兰和扇脉杓兰的传粉生态学

兰科植物精巧的花部结构及其独特的传粉机制为自然选择理论和异花受粉优势学说提供了强有力的证据。开展兰科植物的传粉生态学研究对进一步探讨植物进化中的一些关键问题（如生殖隔离、物种形成、适应和繁育系统进化等）具有重要意义。本文通过对独花兰 (Changnienia amoena) 和扇脉杓兰 (Cypripedium japonicum) 两种兰科植物进行植物与传粉者之间关系、种内形态变异、花粉散布等方面的研究，探讨其适应进化方式、传粉系统的进化趋势，同时为开发适用的分子标记，在独花兰中初步摸索和探讨了微卫星标记的分离。主要研究结果如下： 1. 濒危植物独花兰的传粉生态学 在神农架2个地点进行了连续2年的野外观测和实验。结果表明，独花兰是一种自交亲和、需要昆虫传粉的欺骗性植物。在2个地点独花兰的传粉者种类不同，在龙门河三条熊蜂 (Bombus (Diversobombus) trifasciatus) 是主要传粉者，在关门山仿熊蜂 (B. (Tricornibombus) imitator) 是唯一的传粉者。尽管它们出现的丰度和觅食行为不同，但其传粉行为和携粉部位非常相同，因此，可以看作是一个功能群（functional group）。从种群水平上看，开花个体的分布式样不是显著的偏斜曲线图，传粉者的有效访问主要集中在开花的前、中期。由于没有花蜜，传粉者在花内停留的时间很短，花粉块输出和传粉发生在熊蜂劳而无获的访问、退出花时。自然条件下，独花兰的结实率很低，只有6%-12%，并呈现明显的年份和地点变化。传粉者限制是结实率低的主要原因。与同亚族近缘种布袋兰相比，分布于中国中部的独花兰受冰期的影响很小，其传粉环境相对稳定，在进化历史中，形成了一种稳定的传粉系统。本研究结果纠正了前人对独花兰繁殖方式的错误认识，并为其保护策略的制订提供了重要资料。 2. 独花兰种群大小、传粉者访问和花粉流 为了检验独花兰种群大小、空间格局与传粉者访问之间的关系，连续2年对10个独花兰种群的花粉块输出和输入进行了统计分析，同时对花粉块进行标记研究其散布式样。结果显示，花粉块输出比例与种群大小之间关系在2003年呈显著负相关，但在2004年这种负相关关系不显著。尽管独花兰的空间分布格局在不同种群不同，但传粉者的访问除了在有蜂巢的种群为聚集式访问外，其余种群内均为随机访问式样。花粉块的散布式样呈尖峰态分布，多数花粉块散布很短的距离，少数散布较远。花粉块的平均散布距离在龙门河为7.3 m，在关门山为10.6 m。由于各种群之间相隔很远，种群之间花粉块交流受到限制，很少有种群外的花粉块输入。花粉块传递只在种群内近缘个体中进行，有可能导致种内遗传或形态分化。 3. 独花兰种内形态变异及其适应意义 对庐山、新宁和神农架3个地点15个独花兰自然种群的形态变异进行了研究，探讨了形态多样性水平和地理变异式样及其可能的适应机制。结果表明，在物种水平上独花兰形态性状存在丰富的变异。单因素方差分析显示3个地区间多个形态性状存在极显著差异（P<0.01），UPGMA聚类分析也表明这3个地区分别形成明显不同的分支，说明3个地区种群植物形态已经出现分化。在神农架地区，龙门河和关门山两个地点间出现明显的形态分化，而这两个地点的传粉者在形体大小表现出显著差异。适合度与形态性状之间的相关性分析显示，独花兰种内形态分化是以传粉者为媒介的自然选择作用的结果。移栽实验显示，本地传粉者对本地和移栽种群花的访问表现出一定的选择性。 4. 扇脉杓兰的传粉生态学研究 与独花兰同域分布、花期不遇的扇脉杓兰是杓兰属的一个特殊类群，具有典型“Japonicum”型唇瓣类型。对神农架6个种群连续两年野外观测和实验结果表明，扇脉杓兰是一种自交亲和、需要昆虫传粉的欺骗性植物。在种群水平上开花个体的分布式样呈显著的偏斜曲线图，在短时间内迅速到达盛花期。扇脉杓兰的传粉者是三种熊蜂，其访问频次很低，访问时间很短，有效访问主要集中在开花的前、中期。自然结实率只有4.3％－8.5％，人工授粉实验证明，传粉者限制是结实率低的主要原因，且传粉者限制程度在花期不同阶段和不同地点存在差异。对欺骗性植物聚集生长、开花有利于吸引传粉者这一假说进行检验，结果显示，克隆群丛大小与传粉者有效访问（花粉块输出比例）呈弱负相关关系。聚集生长的结实率与分散生长的没有显著差异，说明繁殖成功与开花个体的密度无关。扇脉杓兰的花粉集结为未蜡质化的花粉块，是杓兰属中粘质花粉和蜡质花粉块的中间进化类型。传粉生态学和形态学证据显示C. acaule和扇脉杓兰是东亚-北美间断分布的姐妹种对。 5. 独花兰微卫星位点的分离 为了深入研究独花兰种群的基因流和交配系统，用Dynabeads和选择性杂交法分离微卫星位点、筛选引物。首先将基因组DNA用合适的内切酶消化为400-1000 bp大小的片段，然后用生物素标记的简单重复寡核苷酸序列做探针与其杂交，杂交复合物用抗生链霉素蛋白包裹的磁珠吸附，经过一系列洗脱、沉淀，得到富含重复片段的DNA，然后在进行克隆、测序，利用SSR两侧翼区设计引物，经过多态性分析可得到微卫星标记。结果显示，35个克隆序列中18个（51.4％）含有重复片段（SSRs），说明用Dynabeads富集效率明显高于传统方法。到目前为止所得引物共4对。这一实验方法的探索为以后分离杓兰属等兰科植物微卫星位点、进而进行群体遗传学研究奠定了一定的基础。

松科nrDNA 的进化研究Evolution of nuclear ribosomal DNA in Pinaceae

核核糖体DNA（nrDNA）已被作为一个重要的标记，用于推断很多分类等级上的系统发育关系。相对于在被子植物中的快速致同进化，nrDNA在裸子植物中的致同进化速率低，且ITS和5S-NTS区有着较大的长度变异，这种现象在松科植物中尤为明显。在本研究中，我们克隆并测定了银杉属的5S rDNA以及冷杉属、银杉属、雪松属、油杉属、长苞铁杉属、金钱松属与铁杉属的ITS序列。基于获得的新数据，再结合前人报导的其它属的数据，我们探讨了如下四个问题： （1）松科 nrDNA ITS1 亚重复单位的组成、分布及进化;（2）ITS1区的长度变异与亚重复单位数目的关系以及它们的系统学意义;（3）松科ITS1的二级结构特征;（4）银杉5S rDNA编码区及非转录间隔区的结构特征。主要研究结果如下： 1. ITS区的序列分析ITS区的克隆及序列分析发现：（1） 松科ITS1的长度变异范围为 944-3271 bp, 这是目前已报导的真核生物中属间ITS变异最大的类群之一;（2） 所有松科植物的ITS区域都包含亚重复单位，亚重复单位的数目从2到9，并且这些亚重复单位可分为两种类型，即不含保守核心序列（5’-GGCCACCCTAGTC ） 的长亚重复单位（LSR）和含上述保守核心序列的短亚重复单位（SSR）;（3） ITS1区的巨大长度变异主要归因于亚重复单位的数量变异; （4） ITS1区的GC含量与 它的序列长度和亚重复单位的数目有一定关系，并能够提供一些系统发育信息，特别是支持云杉属、松属和银杉属三者具有很近的亲缘关系。 2. ITS1亚重复单位的系统发育分析为了研究亚重复单位的进化关系，我们用最大似然法和最大简约法构建了松科ITS1亚重复单位的系统发育树。结果表明：（1）在ML和MP树中可发现有共同的五个分支; （2） 银杉比松科其它属拥有更多的SSR，且该属的所有9个SSR在系统树中构成一个单系支，表明它们是在银杉属内发生重复的;（3）一些SSR在属间和种间具有同源性，可为nrDNA ITS 的进化历史以及松科的系统发育 研究提供重要信息;（4）亚重复单位的多次重复以及伴随的重组可能是导致LSR 和SSR在松科不同属中分布式样不同的原因。 3. 松科ITS1的二级结构 用 Mfold 3.2 软件对松科所有11个属的ITS1区进行了二级结构预测，共获得了563个最低自由能折叠。结合以前关于松科二级结构的报导，我们分析的结果表明：（1） 松科ITS1的二级结构主要由几个延展的发夹结构组成;（2） 构象的复杂性与亚重复的数目呈正相关;（3）配对的亚重复单位通常在保守核心区（5’-GGCCACCCTAGTC ） 处有部分重叠，并且构成一个长茎，而其它的亚重复单位通常会自身折叠，且保守核心区的部分出现在发夹结构的环中。 4. 银杉5S rDNA 序列分析 我们对来自银杉不同群体的3个个体的5S rDNA进行了克隆，共获得 45 条序列，分析结果表明：（1） 绝大多数银杉5S rDNA编码区长度为120 bp, 以GGG 开头，以CTC结尾，编码区出现的碱基替代主要为转换;（2） 银杉与其它裸子植物相比，5S rDNA基因编码区具很高的相似性（90-99%）; （3）间隔区含有一个poly-C和一个poly-T结构、两个TC丰富区以及五个GC丰富区。根据长度和序列特征，银杉的5S rDNA间隔区可分为三种类型：Type A 长751-764 bp，Type B 长770-807 bp （含一个32 bp的插入），Type C 长581-594 bp; （5）长间隔区（Type A，Type B ）中含有两个148-175 bp的串联亚重复单位，该亚重复单位与短间隔区（Type C ）中的一段143 bp的序列具有较高的相似性（56.0-66.8%）。 5. 银杉5S rRNA的二级结构 Mfold 3.2 预测结果表明：（1）银杉5S rRNA二级结构包括5个双螺旋区（干区）（Ⅰ-Ⅴ）、2个发夹结构环区（C和D）、3个中间环区（B1、B2 和 E）和1个铰链区（A）, 铰链区为三个双螺旋的结合处;（2） 二级结构中的环区通常比双螺旋区更加保守;（3）在5个双螺旋中，I 和 IV 区有较高的碱基替代率。

Isolation and characterization of microsatellite markers for the endangered Taxus yunnanensis

Eleven polymorphic microsatellite loci were characterized for the endangered conifer Taxus yunnanensis. Eight loci were isolated through SSR-anchored PCR, one locus was developed by cross-species amplification tests, while the last two loci were obtained

Genetic diversity study of Caspian brown trout (Salmo trutta caspius Kessler, 1877) population in five rivers in the southern part of Caspian Sea, Iran

For study the genetic diversity of Caspian brown trout population in five rivers in the southern part of Caspian Sea in Iran 182 number generators in the fall and winter of 1390 were collected in Chalus, Sardab Rud, Cheshmeh Kileh, Kargan Rud and Astara rivers. Then about 3-5 g of soft and fresh tissue from the bottom fin fish removed and were fixed in ethanol 96°. Genomic DNA was extracted by using ammonium acetate, then quantity and quality of the extracted DNA were determined by using spectrophotometry and horizontal electrophoresis in 1% agarose gel. The polymerase chain reaction was performed by using 16 SSR primers and sequencing primers (D-Loop) and the quality of PCR products amplified by SSR method were performed by using horizontal electrophoresis in 2% agarose gel. Alleles and their sizes were determined by using vertical electrophoresis in 6% polyacrylamide gel and silver nitrate staining method. Gel images were recorded by gel documentarian, the bands were scored by using Photo- Capt software and statistical analysis was performed by using Gene Alex and Pop Gene software. Also the PCR sequencing products after quality assessment by usinghorizontal electrophoresis in 1.5% agarose gel were purified and sent to South Korea Bioneer Corporation for sequencing. Sequencing was performed by chain termination method and the statistical analysis was performed by using Bio- Edit, Mega, Arlequin and DNA SP software. The SSR method, 5 pairs of primers produced polymorphic bands and the average real and effective number of alleles were calculated 5.60±1.83 and 3.87±1.46 in the Cheshmeh Kileh river and 7.60±1.75 and 5.48±1.32 in the Karganrud river and the mean observed and expected heterozygosity were calculated 0.44 ±0.15 and 0.52 ±0.16 in the Cheshmeh Kileh river and 0.50 ±0.11 and 0.70±0.13 in the Karganrud river. Analysis of Molecular Variance results showed that significant differences in genetic diversity between and within populations and between and within individuals in the studied rivers (P<0.01). The sequencing method identified 35 different haplotype, the highest number of polymorphic position (251) and haplotype (14) were observed in the Chalus river. The highest mean observed number of alleles (2.24±0.48) was calculated in the Sardabrud river, the highest mean observed heterozygosity (1.00±0.03) was calculated in the Chalus river and the highest mean nucleotide diversity (0.13±0.07) was observed in the Sardabrud river and mean haplotype diversity was obtained (1) in three studied rivers. The overall results show that there are no same population of this fish in the studied rivers and Karganrud and Chalus rivers in the SSR and sequencing methods had the highest levels of genetic diversity.

Population genetic structure of Rutilus frisii kutum in Caspian Sea using biometry and microsatellite molecular technique

In order to carry out Biometric studies, 75 samples were caught from 3 locations ( Tajan river, Sefidrud and Shirud) using Salic and the length (±1 mm) and weights (± 5 gr) of samples were determined. Using One-way ANOVA by SPPSS software, there wasn’t significant difference between locations in length and fecondity (P ≥0.01(, but there was significant difference between Shirud and tajan samples with sefidrud in weight ) P≤0.01(. In order to carry out genetic variation studies, 210 fish were caught from 3 different regions of the Iranian coastline (Khoshkrud, Tonekabon, Gorganrud) and 1 region in Azerbaijan (Waters of the Caspian Sea close to Kura River mouth) during 2008-2009 . Genomic DNA was extracted of fin using the phenol-chloroform. The quantity and quality of DNA from samples were assessed by spectrophptometer and 1% agarose gel electro-phoresis. PCR was carried out using 15 paired microsatellite primers. PCR products were separated on 8% polyacrylamide gels that were stained using silver nitrate. Molecular weight calculate using UVTech software. The recorded microsatellite genotypes were used as input data for the GENALEX software version 6 package in order to calculate allele and genotype frequencies, observed (Ho) and (He) expected heterozygosities and to test for deviations from Hardy-Weinberg equilibrium. Genetic distance between two populations was estimated from Nei standard genetic distance and genetic similarity index (Nei, 1972). Genetic differentiation between populations was also evaluated by the calculation of pairwise estimates of Fst and Rst values. From 15 SSR markers were used in this investigation, 9 of them were polymorph. Average of expected and observed heterozygosity was 0.54 and 0.49 respectively. Significant deviations from Hardy-Weinberg expectations were observed in all of location except Anzali lagoon- autumn in AF277576 and EF144125, Khoshkrud in EF144125 and Gorganrud and Kura in AF277576. Using Fst and Rst there was significant difference between locations ) P≤0.01(. According to Fst , the highest population differentiation (Fst= 0.217) was between Gorganrud and Khoshkrud that have the lowest Nm and the lowest (Fst= 0.086) was between Gorganrud and Tonekabon that have the highest Nm. Using Rst the highest population differentiation (Rst= 0.271) was between Tonekabon and spring Anzali lagoon and the lowest (Rst= 0.026) was between Tonekabon and Autumn Anzali 159 lagoon. Also the difference between Spring Anzali lagoon and Autumn Anzali lagoon was noticeable (Fst=0.15). AMOVA analysis with consideration of 2 sampling regions (Iran and Azerbaijan) and 7 sampling locations (Iran: Khoshkrud, Tonekabon, Gorganrud, Spring Anzali lagoon and Autumn Anzali lagoon ; Azerbaijan: the Kura mouth) revealed that almost all of the variance in data namely 83% )P≤0.01( was within locations, Genetic variances among locations was 14% )P≤0.01( and among regions was 3% )P≤0.01(. The genetic distance was the highest (0.646) between Gorganrud and Autumn Anzali lagoon populations, whereas the lowest distance (0.237) was between Gorganrud and Tonekabon River. Result obtained from the present study show that at least 2 different population of Rutilus frissi kutum are found in the Caspian sea,which are including the kura river population and the southern Caspian sea samples and it appears that there is more than one population in southern Caspian sea that should be attantioned in artifical reproduction Center and stoke rebilding.

Mouse Toxicity of Anabaena flos-aquae from Lake Dianchi, China

Some species of the genera Anabaena can produce various kinds of cyanotoxins, which may pose risks to environment and human health. Anabaena has frequently been observed in eutrophic freshwater of China in recent years, but its toxicity has been reported only in a few studies. In the present study, the toxicity of an Anabaena flos-aquae strain isolated from Lake Dianchi was investigated. Acute toxicity testing was performed by mouse bioassay using crude extracts from the lyophilized cultures. The mice exposed to crude extracts showed visible symptoms of toxicity and died within 10-24 h of the injection. Serum biochemical parameters were evaluated by the use of commercial diagnostic kits. Significant alterations were found in the serum biochemical parameters: alkaline phosphatase (AKP), gamma-glutamyl transpepticlase (gamma-GT), aspartate amino transferase (AST), alanine amino transferase (ALT), AST/ALT ratio, total protein content, albumin content, albumin/globulin (A/G) ratio, blood urea nitrogen (BUN), serum creatinine (Ssr), and total antioxidative capacity (T-AOC). Histopathological observations were carried out with hematoxylin and eosin (HE) stain under light microscope. Severe lesions were seen in the livers, kidneys, and lungs of the mice injected with crude extracts. The alterations of biochemical parameters were in a dose-dependent manner, and the severities of histological lesions were in the same manner. Based on biochemical and histological studies, this research firstly shows the presence of toxin-producing Anabaena species in Lake Dianchi and the toxic effects of its crude extracts on mammals. (C) 2008 Wiley Periodicals, Inc. Environ Toxicol 24: 10-18, 2009.

Development of EST-SSRs by an efficient FIASCO-based strategy: A case study in rare minnow (Gobiocyrpis rarus)

Type I markers are useful for comparative mapping and other genetic analyses, but relatively difficult to develop. In the present study a microsatellite (SSR)-enriched cDNA library was constructed for the first time using the fast isolation by AFLP of sequences containing repeats (FIASCO) method in a small fish, Chinese rare minnow (Gobiocypris rarus). A total of 97.4% of the expressed sequence tags (ESTs) contained targeted CA-repeats, in which 29 unique EST-SSRs were identified. Ten out of the 28 loci for which primer pairs were designed were polymorphic with alleles ranging from three to seven (mean 4.50). Some of these EST-SSRs can be amplified in other species. These results proved that cDNA-FIASCO is an efficient way to isolate novel EST-SSRs in a fish.

Construction and characterization of a Lipotes vexillifer genomic DNA BAC library

We constructed a genomic DNA library for Lipotes vexillifer (L. vexillifer), the Baiji or Yangtze River dolphin, one of the most endangered mammals in the world. The library consists of 149,000 BAC clones, with an average insert size of 83 kb, representing approximately 3.4 haploid genome equivalents. PCR amplification of four known L. vexillifer genes yielded two to four positive clones each. To demonstrate the utility of this library, we isolated and sequenced the L. vexillifer alpha lactalbumin gene, which is a gene specific to mammals and one which has been widely used as molecular tool in phylogenetic analysis. We also end-sequenced 20 randomly selected clones, resulting in the identification of at least five new L. vexilliter genes, five SSR loci, and one SINE locus. These results suggest that this library is a valuable resource for candidate gene cloning, physical mapping, and genome sequencing of this important and threatened species.

A genetic map of large yellow croaker Pseudosciaena crocea

Genetic linkage maps were constructed for large yellow croaker Pseudosciaena crocea (Richardson, 1846) using AFLP and microsatellite markers in an F-1 family. Five hundred and twenty-three AFLP markers and 36 microsatellites were genotyped in the parents and 94 F-1 progeny. Among these, 362 AFLP markers and 13 SSR markers followed the 1:1 Mendelian segregation ratio (P > 0.05). The female genetic map contained 181 AFLP and 7 microsatellite markers forming 24 linkage groups spanning 2959.1 cM, while the male map consisted of 153 AFLP and 8 microsatellite markers in 23 linkage groups covering 2205.7 cM. One sex linked marker was mapped to the male map and co-segregated with the AFLP marker agacta355, suggesting an XY-male determination mechanism and this may be useful in the breeding of monosex populations. (c) 2007 Elsevier B.V. All rights reserved.

基于安全状态域的网络评估模型

将基于攻击图的评估与依赖标准的评估相结合,提出了一种基于安全状态域(security state region,简称SSR)的网络安全评估模型(security-state-region-based evaluation model,简称SSREM).该模型将攻击的影响分为攻击能力改变和环境改变,通过两者之间的因果关系建立数学模型,提出了安全状态域趋向指数的概念,借助Matlab进行攻击趋势的曲面拟合,进而进行安全状态域的划分和网络的安全性评估.实验结果表明,依据SSREM进行的评估能够通过安全状态城和安全状态域趋向指数反映网络进入不同状态的难易程度,对网络安全性量化评估具有借鉴意义.

A linkage map of common carp (Cyprinus carpio) based on AFLP and microsatellite markers

P>Common carp (Cyprinus carpio) is an important fish for aquaculture, but genomics of this species is still in its infancy. In this study, a linkage map of common carp based on Amplified Fragment Length Polymorphism (AFLP) and microsatellite (SSR) markers has been generated using gynogenetic haploids. Of 926 markers genotyped, 151 (149 AFLPs, two SSRs) were distorted and eliminated from the linkage analyses. A total of 699 AFLP and 20 microsatellite (SSR) markers were assigned to the map, which comprised 64 linkage groups and covered 5506.9 cM Kosambi, with an average interval distance of 7.66 cM Kosambi. The normality tests on interval map distances showed a non-normal marker distribution. Visual inspection of the map distance distribution histogram showed a cluster of interval map distances on the left side of the chart, which suggested the occurrence of AFLP marker clusters. On the other hand, the lack of an obvious cluster on the right side showed that there were a few big gaps which need more markers to bridge. The correlation analysis showed a highly significant relatedness between the length of linkage group and the number of markers, indicating that the AFLP markers in this map were randomly distributed among different linkage groups. This study is helpful for research into the common carp genome and for further studies of genetics and marker-assisted breeding in this species.

中国卧龙自然保护区不同海拔川滇高山栎(Quercus aquifolioides)群体的遗传变异

海拔梯度造成的环境异质性，如崎岖的地形、复杂的植被结构以及花期延迟等可能会极大地影响到物种的形态和遗传变异格局。理解物种形态和遗传变异的海拔格局对于物种多样性的管理和保护是非常重要的。尽管植物群体遗传学是一个飞速发展的研究领域，然而与海拔相关的形态变异、遗传变异及群体间遗传差异的研究却很少。到目前为止，还不清楚遗传变异与海拔之间是否必然的相关性。 川滇高山栎是一种重要的生态和经济型树种，广泛分布于中国西南的四川、西藏、贵州和云南省的高海拔地区，在保持水土、调节气候方面起着十分重要的作用。尽管主要受阳光限制而仅分布于阳坡，但其海拔梯度范围较大，表明川滇高山栎对不同的环境具有很强的适应性。本文通过叶型及生理响应、微卫星分子标记和扩增性片段长度多态性方法，试图探索川滇高山栎叶沿海拔梯度的形态和生理响应及其沿海拔梯度的遗传变异格局，为川滇高山栎的保护和利用提供进一步的遗传学理论依据和技术指导。 对叶形、含氮量及碳同位素的试验结果表明，平均比叶面积、气孔密度、气孔长度和气孔指数等气孔参数随海拔的升高呈非线性变化。在海拔大于2800 m时，川滇高山栎的比叶面积、气孔长度和气孔指数都随海拔升高而降低，但是在海拔小于2800 m时，这些指标都随海拔的升高而增大。相对而言，单位叶面积的含氮量和碳同位素则表现出相反的变化模式。另外，比叶面积是决定碳同位素沿海拔梯度变化的最重要参数。本研究结果表明，海拔2800 m附近是川滇高山栎生长和发育的最适地带，在这里生长的植物叶片厚度更薄、气孔更大、叶碳同位素值更小。 利用六对微卫星引物对五个不同海拔川滇高山栎群体遗传多样性进行研究，结果表明，群体内表现出较高的遗传多样性，平均每位点等位基因数11.33个，平均期望杂合度达0.820。群体间差异较小，分化仅为6.6％。聚类分析也并没有显示出明显的海拔格局。然而低频率等位基因却与海拔呈显著性正相关(R2=0.97, P < 0.01)，表明在高海拔处，川滇高山栎以更多的稀有基因来适应恶劣的环境条件。本试验结果表明由海拔梯度形成的选择性压力对川滇高山栎群体的遗传变异影响并不明显。 为了进一步探讨川滇高山栎群体遗传变异与海拔之间的相互关系，我们还对其进行了扩增性片段长度多态性分析。结果表明：(1)随海拔的升高(从群体WL2到群体WL5)，群体内遗传变异降低，而群体间遗传差异增加；(2)低海拔群体WL1表现出最低的遗传变异性(HE = 0.181)，同时与其余四个群体间呈现出最大的遗传差异性(平均FST = 0.0596)；(3)在除去低海拔群体WL1后，Mantel检测表明群体间遗传距离与海拔距离之间表现出正相关性。另外，研究结果还表明，遗传变异受生境条件(过度的湿热环境)及人为干扰(火烧、砍伐和放牧)的影响，这一点至少在低海拔群体WL1上发生了作用。 通过叶形态、生理及DNA分子水平的研究，结果表明叶形态特征和碳同位素与海拔紧密相关，与海拔之间呈非线性变化，海拔2,800 m附近是川滇高山栎生长和发育的最适地带。海拔梯度在一定程度上会影响到川滇高山栎群体的遗传变异结构，但在这样一个狭窄的地理分布区域里，这种影响并不足以导致群体间较大的遗传分化。同时生境条件及人为干扰也是影响遗传变异的限制性因子，不容忽视。 Altitudinal gradients impose heterogeneous environmental conditions, such as rugged topography, a complex pattern of vegetation and flowering delay, and they likely furthermore markedly affect the morphological and genetic variation pattern of a species. Understanding altitudinal pattern of morphological and genetic variation at a species is important for the management and conservation of species diversity. Although plant population genetics is a fast growing field of research, there are only few recent investigations, which analyzed the genetic differentiation and changes of intra-population variation along altitudinal gradients. At present, it is still unclear whether there are some common patterns of morphological and genetic variation with altitude. Quercus aquifolioides Rehder & E.H. Wilson, which is an important ecological and economical endemic woody plant species, is widely distributed in the Yunnan and Sichuan provinces, Southwest China. Its large range of habitat across different altitudes implies strong adaptation to different environments, although it is mainly restricted to sunny, south facing slopes. It plays a very important role in preventing soil erosion, soil water loss and regulating climate, as well as in retaining ecological stability. In this paper, we tried to understand the altitudinal pattern of morphological and genetic variation along altitudinal gradients through the experiments of leaf morphological and physiological responses, microsatellite analysis and AFLP markers. In leaf morphological and physiological responses experiment, we measured leaf morphology, nitrogen content and carbon isotope composition (as an indicator of water use efficiency) of Q. aquifolioides along an altitudinal gradient. We found that these leaf morphological and physiological responses to altitudinal gradients were non-linear with increasing altitude. Specific leaf area, stomatal length and index increased with increasing altitude below 2,800 m, but decreased with increasing altitude above 2,800 m. In contrast, leaf nitrogen content per unit area and carbon isotope composition showed opposite change patterns. Specific leaf area seemed to be the most important parameter that determined the carbon isotope composition along the altitudinal gradient. Our results suggest that near 2,800 m in altitude could be the optimum zone for growth and development of Q. aquifolioides, and highlight the importance of the influence of altitude in research on plant physiological ecology. Genetic variation and differentiation were investigated among five natural populations of Q. aquifolioides occurring along an altitudinal gradient that varied from 2,000 to 3,600 m above sea level in the Wolong Natural Reserve of China, by analyzing variation at six microsatellite loci. The results showed that the populations were characterized by relatively high intra-population variation with the average number of alleles equaling 11.33 per locus and the average expected heterozygosity (HE) being 0.779. The amount of genetic variation varied only little among populations, which suggests that the influence of altitude factors on microsatellite variation is limited. However, there is a significantly positive correlation between altitude and the number of low-frequency alleles (R2=0.97, P < 0.01), which indicates that Q. aquifolioides from high altitudes has more unique variation, possibly enabling adaptation to severe conditions. F statistics showed the presence of a slight deficiency of heterozygosity (FIS=0.136) and a low level of differentiation among populations (FST=0.066). The result of the cluster analysis demonstrates that the grouping of populations does not correspond to the altitude of the populations. Based on the available data, it is likely that the selective forces related to altitude are not strong enough to significantly differentiate the populations of Q. aquifolioides in terms of microsatellite variation. To further elucidate genetic variation pattern of Q. aquifolioides populations under sub-alpine environments, genetic variation and differentiation were investigated along altitudinal gradients using AFLP markers. The altitudinal populations with an average altitude interval of 400 m, i.e. WL1, WL2, WL3, WL4 and WL5, correspond to the altitudes 2,000, 2,400, 2,800, 3,200 and 3,600 m, respectively. Our results were as follows: (i) decreasing genetic variation (ranging from 0.253 to 0.210) and increasing genetic differentiation with altitude were obtained from the WL2 to the WL5 population; (ii) the WL1 population showed the lowest genetic variation (HE = 0.181) and the highest genetic differentiation (average FST = 0.0596) with the other four populations; (iii) the positive correlation was obtained using Mantel tests between genetic and altitude distances except for the WL1 population. Our results suggest that altitudinal gradients may have influenced the genetic variation pattern of Q. aquifolioides populations to some extent. In addition, habitat environments (unfavorable wet and hot conditions) and human disturbances (burning, grazing and felling) were possible influencing factors, especially to the low-altitude WL1 population. The present study shows that there were close correlations between morphological features and carbon isotope composition in our data. This indicates that a coordinated plant response modified these parameters simultaneously across different altitudes. Around 2,800 m altitude there seems to be an optimum zone for growth and development of Q. aquifolioides, as indicated by thinner leaves, larger stomata and more negative d13C values. All available evidence indicates altitudinal gradients may have influenced the genetic variation pattern of Q. aquifolioides to some extent. Decreasing genetic variation and increasing genetic differentiation with altitude was obtained except for the WL1 population. And the environment of habitats and human disturbances were also contributing factors, which impact genetic variation pattern, especially to the low-altitude WL1 population.

青藏高原东缘青杨（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.

粗枝云杉( Picea asperata Mast.)天然群体的遗传变异

西南地区在我国的经济发展和生态环境建设中占重要地位，但也是我国生态环境最脆弱的地区之一，生态系统退化，生态功能减弱，严重制约着西南林业的可持续经营与发展。本项目采用DNA 分子标记SSR 研究不同生境条件下粗枝云杉群体的遗传变异及其时空分布格局，考察遗传变异与复杂的山地生态环境间的潜在联系，系统地揭示粗枝云杉天然群体与环境系统相互作用的生态适应与分子进化机制。粗枝云杉适应性强，生长迅速，在植树造林和工业用材方面占有重要地位，研究成果可为中国西南部亚高山天然林的可持续经营及退化生态系统的恢复与重建提供理论依据和科学指导。主要研究结果如下： 1. SSR 位点变异丰富，等位基因频率的分布格局多样。7 个SSR 标记全是多态位点，每位点的等位基因数变化范围为13~24，平均为19.9 个。SSR 位点的等位基因片段长度范围变化较大。73.1%的等位基因变异遵循逐步突变模型(SSM)而发生1 个重复基元的变化，22.3%和4.6%的变异分别按两阶段突变模型(TMP)发生1 个重复基元以上的变化和在SSR 位点侧翼区发生1 个碱基变化的插入-删除事件。 2. 粗枝云杉拥有中等偏高水平的遗传多样性和相对大的群体间遗传分化。通过分析代表10 个群体的250 个个体在7 个SSR位点的变化，调查了源自中国西南山区的粗枝云杉的微卫星变异。相当高的遗传多样性和强烈的群体分化发生在粗枝云杉中， 其群体平均Nei's 期望杂合度为0.707 ， 群体间遗传距离为0.121~0.224(FST)和0.100~0.537(RST)。然而，群体间遗传距离与地理距离之间无相关性，从而排除了简单的距离分离模式并暗示迁移不是影响粗枝云杉遗传变异格局的主要因素。事实上，使用私有等位基因估算的基因流数量非常低，仅等于0.753。等位基因置换检验(Allele permutation tests)揭示逐步突变及遗传漂变都对群体间分化有贡献。另外，在多数位点检测到显著的群体间遗传差异，这个结果说明自然选择，假设通过环境压力，是引起粗枝云杉微地理分化的主要因素之一。根据SSR基因型，250 个粗枝云杉个体的70%被正确地归类入其各自的来源群体，结果表明微卫星(SSR)对区分来自中国不同生态地理位点的粗枝云杉基因型是有效的。 3. 在SSR、RAPD 和AFLP 位点，显著的群体间遗传结构被发现的，但三种标记间遗传分化程度和群体遗传关系有差异。利用来自10 个群体的247 个个体，我们报告了关于样本粗枝云杉群体间遗传关系的总体看法。根据各自对评价遗传关系的信息能力和适用性，SSR、RAPD 和AFLP 标记被选用，三种技术非常有效地区别这些基因型。使用的SSR、RAPD 和AFLP 标记分别估计平均Dice 相似性系数。Mantel 检验产生显著但相对低的共表型适合度(RAPD = 0.63£AFLP = 0.60和SSR = 0.75)。比较三种标记系统，RAPD 和AFLP 共表型指数相对高地相关(r =0.59)，而RAPD 和SSR 及SSR 和AFLP 之间的相关系数分别是0.53 和0.35。所有系统树，包括不同标记资料结合获得的系统树，反映了多数群体依据它们的地理条件而成某种特定关系。结果暗示单个或结合标记系统能用来深入洞察粗枝云杉遗传研究，并且不同标记系统合并资料能提供更可靠的信息。 Southwestern region plays an important role in economic developmentand ecological construction in China. Yet, it is also one of the weak regionsof ecological environment in China with degraded ecosystem and imperfectfunction, which restricts the sustaining management and development ofsouthwestern forestry. The genetic variation and spatial distribution patternof P. asperata populations originating from different habitats wereinvestigated using SSR molecular markers in this study. The correlationsbetween genetic variation and ecological and environmental conditionswere detected, and the interaction between P. asperata populations andenvironmental system and the mechanism of ecological adaption -molecular evolution were revealed. Given the significant ecological andeconomic roles of the fast-growing and wide-adaptive species in reforestation and production of pulp wood and timber, the study couldprovide a strong theoretical evidence and scientific direction for thesustaining management of subalpine natural forest, and the afforestationand rehabilitation of degraded ecosystem. The results are as follows: 1. The genetic variation at SSR loci was abundant and the distributionof allelic frequencies was uneven. All seven loci were polymorphic, and thenumber of alleles per locus varied from 13 to 24 with a mean valueequaling 19.9. The allele sizes at SSR loci were found to vary widely.73.1% of allelic variation followed stepwise mutation model (SSM) whichresults increase or decrease by one repeat type, and 22.3% and 4.6% wereresulted from two-phase mutation model (TMP) with allele size varying bymore than one repeat type and from insertion-deletion events in theflanking regions at SSR loci with a single basepair changing, respectively. 2. P. asperata possessed a moderate to high level of genetic diversityand considerable genetic differentiation. Microsatellite variation of P.asperata. originating from the mountains of southwestern China wasinvestigated by analyzing variation at seven SSR loci in 250 individualsrepresenting ten populations. A fair degree of genetic diversity and strongpopulation subdivision occurred with the mean gene diversity (H) of 0.707,and genetic distances among populations varying between 0.121 and 0.224(FST) and between 0.100 and 0.537 (RST). However, inter-populationgenetic distances showed no correlation with geographic distances between the population sites. This ruled out a simple isolation by distance modeland suggested that migration does not have a great impact. In fact, theamount of gene flow, detected using private alleles, was very low, equalingonly 0.753. Allele permutation tests revealed that stepwise-like mutations,coupled with genetic drift, could contribute to population differentiation.Moreover, significant genetic differences between populations weredetected at most loci. The results indicate that natural selection, presumablythrough environmental stress, may be one of the main factors causingmicro-geographical differentiation in the genetic structure of P. asperata.Based on SSR genotypes, 70% of the 250 individuals were correctlyclassified into their sites of origin. This suggests that microsatellites (SSRs) are effective in distinguishing genotypes of P. asperata originating fromdiverse eco-geographical sites in China. 3. Using a set of 247 individuals from ten P. asperata populations wereport an overview on the genetic relationship among the sampled P.asperata populations. RAPD, AFLP and SSR were used in terms of theirinformativeness and applicability for evaluate relationship and all threetechniques discriminated the genotypes very effectively. Mean Dicesimilarities coefficient were estimated using RAPD, AFLP and SSR,respectively. The Mantel test resulted in a significant but relatively low fit(RAPD = 0.63, AFLP = 0.60 and SSR = 0.75) of cophenetic values.Comparing the three marker systems to each other, RAPD and AFLP cophenetic indices were highly correlated (r = 0.59), while correlationcoefficient between RAPD and SSR was r = 0.53 and between SSR andAFLP was r = 0.35. For all markers a relatively high similarity indendrogram topologies was obtained although some differences wereobserved. All the dendrograms, including that obtained by the combineduse of all the marker data, reflect some relationships for most of thepopulations according to their geographic conditions. The results indicatethat single or combined marker system could be used to insight into geneticstudy in P. asperata and the combined data of different marker systems canprovide more reliable information.