基于高通量测序的思茅松微卫星位点的特征分析

以思茅松针叶为材料,采用Illumina Hiseq2000平台测序,得到121 882条无冗余的序列,对这些序列进行SSR位点搜索,共获得3534个SSR位点,出现频率为2.9%,分布平均距离25kb。所搜索的SSR以单核苷酸重复类型最多,其次为三核苷酸和二核苷酸,而四、五、六核苷酸重复类型较少(1%)。单核苷酸重复类型中以A/T基元较丰富(43.58%);二核苷酸中出现频率最高的是AT/TA基元(14.04%),AG/CT次之,AC/GT和CG/CG较低;三核苷酸重复类型中AAG/CTT最多,AGC/CTG和AGG/CCT次之,以CCG/CGG、ACG/CGT和ACT/AGT较少(1%);四、五、六核苷酸类型中各重复基元相差不大,均较少。SSR数量随对应的重复类型重复次数的增加而降低,也随重复区段碱基长度的增加而降低。     

云南金花茶转录组SSR的分布及其序列特征

分析云南金花茶转录组中的SSR位点信息及其序列分布特征,为SSR引物的大量开发奠定前期基础。利用MISA软件分析了云南金花茶转录组数据155 011条Unigenes的SSR位点信息,在云南金花茶转录组中共检测到30 435个SSR位点,出现频率为19.63%,平均3.46 kb就出现1个SSR位点。而重复基元类型从二核苷酸到六核苷酸均有出现,其中二核苷酸重复基元为主要类型,占SSR总数的71.44%;其次是三核苷酸,占SSR总数的25.48%;四核苷酸到六核苷酸仅占3.08%。AG/TC与CT/GA是二核苷酸的优势重复基元,分别占总基元的26.75%与25.17%;CTT/GAA与AGA/TCT是三核苷酸的优势重复基元,分别占总数的2.64%与2.23%。根据云南金花茶转录组测序结果,云南金花茶转录组的SSR位点出现的频率较高,具有丰富的基元类型,并且有较高的重复次数,具有较高的多态性潜力。随机选出50对设计的SSR引物用于检测,其中有14对引物能扩增出清晰且单一的谱带。对云南金花茶转录组数据中SSR位点信息进行分析,可为其SSR引物的设计与筛选、SSR遗传多样性分析以及种质资源的保护提供理论依据。     

银杏转录组SSR位点分布及其序列特征分析

【目的】为了探究银杏叶片转录组中简单重复序列(Simple sequence repeats,SSR)的位点信息,以丰富银杏分子标记库,为日后银杏SSR分子标记的筛选和开发、SSR遗传多样性分析、遗传图谱构建和品种鉴定及种质资源保护奠定理论基础。【方法】利用Illumina HiSeq~(TM)2500高通量测序平台获得银杏叶片转录组数据,然后采用MicroSAtellite(MISA)软件对测序获得的299 373条Unigenes序列SSR位点的分布特征进行分析;再利用转录组数据开发SSR分子标记,对6份不同品种(单株)的银杏叶片材料进行扩增和筛选。【结果】通过组装得到了299 373条Unigenes序列,其总长为232 983.457 kb,其中,G+C的占比为42.82%。通过检测,共搜索到SSR位点17 821个,SSR位点的出现频率为5.95%,即平均长为13.07 kb的序列中就含有1个SSR位点,其中,有16 163条Unigenes序列含有1个以上的SSR位点;SSR位点重复类型中最主要的类型是二核苷酸和三核苷酸,分别占SSR位点总数的74.46%和23.37%,其中的AG/CT、AAG/TTC分别是二核苷酸、三核苷酸的优势重复单元,分别占SSR位点总数的31.19%和4.65%;对5种不同SSR位点重复类型的重复次数的统计结果表明,重复次数为5～11次的不同SSR位点类型占SSR位点总数的90.21%;基序长度为10～20 bp的SSR位点数占SSR位点总数的83.17%。利用6个不同品种(单株)的银杏叶片材料对所有引物进行筛选,得到了条带清晰、稳定且多态性好的引物。【结论】根据转录组测序数据中不同SSR位点的分析结果可知,转录组中SSR位点的出现频率相对较高、平均距离较短,说明转录组中SSR位点的数量和种类均较多。随着重复次数的增加,SSR位点的出现频率随之降低。重复次数和基序长度也都说明了转录组测序所得的SSR位点大部分具有多态性的潜能。试验筛选出了条带清晰、稳定且多态性好的引物,其不仅适用于后续的银杏遗传多样性分析,还可以用于银杏的遗传图谱构建和品种鉴定等方面的研究之中。     

杨树抗逆转录因子基因内SSR标记的开发

利用直接测序法开发小叶杨抗逆转录因子基因内一套新的核基因组SSR标记。通过对3个抗逆转录因子共21个成员在36个小叶杨基因型个体中的序列比较分析后共检测到31个SSR多态性位点,SSR出现的频率为1/1916bp。在小叶杨自然群体中,SSR多态性位点的碱基重复呈现出2～5碱基形式,基元重复次数变异范围为3～20次,其中以二碱基重复的位点较多,占总数的51.6%。在此基础上,依据SSR位点两侧的保守序列,设计31对SSR位点PCR扩增引物对。利用设计的引物检测所开发的SSR位点在杨属内22个基因型个体中PCR扩增的有效性及SSR位点的保守性。PCR扩增结果显示,93.5%的SSR位点能够在杨属内至少4个派内有效扩增,每对引物组合可检测到SSR多样性位点数3～12个,平均6个。基于抗逆转录因子基因内开发的SSR标记位点为分子标记辅助小叶杨抗逆性状育种提供工具。     

基于锦绣杜鹃花蕾转录组的SSR标记开发及应用

[目的]探究锦绣杜鹃(Rhododendron pulchrum Planch.) EST-SSR的类型、分布频率和分布特征,开发有效的SSR标记,并验证其在遗传多样性研究和跨物种转移中的应用潜力。[方法]采用RNA-seq技术对锦绣杜鹃‘紫鹤’品种的花蕾进行转录组测序,MISA软件对组装的unigenes内部的SSR位点进行检索,分析其类型、分布频率。利用Primer 3.0软件设计引物,并对锦绣杜鹃群体和其近缘种映山红群体进行多样性检测,利用POPGENE-PC 2.2软件计算等位基因数、有效等位基因数、多态性位点百分率、Shannon’s信息指数、Nei氏多样性指数、观察杂合度、期望杂合度等参数。[结果]锦绣杜鹃‘紫鹤’花蕾转录组共有49 527个unigenes(43 766 249 bp),从中筛选出16 120个SSR位点(24.46%),发生频率为1·(2.7 kb)~(-1)。微卫星的重复次数主要集中在5～24次之间,二核苷酸发生频率最高(9 624个,59.70%),其次是单核苷酸(3 738个,23.19%),频率最低的为五核苷酸(42个,0.26%)。频率最高的重复基序有A/T、AG/CT、AAG/CTT、AGG/CCT、ACC/GGT、AGC/GCT、AAAG/CTTT、AAGAG/CTCTT、AGAGGG/CCCTCT等。选用13个多态性SSR标记对锦绣杜鹃群体进行PCR扩增,共得到71个等位基因,平均每个位点3～9个;有效等位基因数在各SSR位点之间变化范围为1.684～5.930;观察杂合度H_O和期望杂合度H_E的变化范围分别为0.000～1.000和0.433～0.848,平均值分别为0.696±0.426和0.705±0.129;Shannon’s信息指数I和Nei氏多样性指数h的变化范围分别为0.736～1.961和0.406～0.831,平均值分别为1.376±0.339和0.683±0.131。此13个标记在映山红群体中的跨物种扩增成功率为100%,并反映出丰富的遗传多样性。2个群体间的遗传差异93.4%存在于群体内部。[结论]基于AG/CT重复基序开发的13个SSR标记具有高度的多态性,为后续杜鹃花属植物的遗传多样性研究、遗传图谱构建、基因定位、分子标记辅助育种等研究奠定了基础。大别山野生映山红资源具有较高的遗传多样性,且杂合度过剩,遗传变异主要发生在群体内部。     

穗花杉的转录组测序及其转录组特性分析

[目的]通过对红豆杉科穗花杉属的穗花杉进行转录组测序,为穗花杉的萜类合成途径和分类学研究提供支持。[方法]利用Hi Seq2500技术对穗花杉的茎叶进行转录组测序。[结果]穗花杉转录组测序共得到8.14 Gb的clean data。从头组装共获得82 884条unigene,其中有27 495条unigene被注释。此外,在82 884条unigene中共搜索到2 827个SSR位点,单核苷酸重复SSR出现频率最高(60.1%),其次为三核苷酸重复SSR(25.4%)。在穗花杉unigene中挖掘出了1个牦牛儿基焦磷酸合成酶(GPS)、1个法尼烯基焦磷酸合成酶(FPS)和5个牦牛儿基牦牛儿基焦磷酸合成酶(GGPS)的同源基因,同时得到紫杉二烯合成酶(TS)同源基因13个。[结论]发现了20个与萜类合成有关的unigene和2 827个SSR位点,为今后穗花杉萜类生物合成研究,特别是紫杉二烯合成酶编码基因的研究打下了基础,也为该物种系统分类地位及遗传多样性研究提供了新的遗传信息。     

基于转录组序列的楠木SSR分子标记开发

【目的】通过对楠木转录组数据的分析,开发楠木SSR分子标记技术,为楠木遗传多样性分析和资源保护提供保障。【方法】将楠木高通量转录组测序获得的67 331条Unigene进行简单重复序列(SSR)位点挖掘和分析,并结合引物验证,初步证明其可行性。【结果】通过软件分析,共获得9 405个SSR位点,出现频率为13.97%,涉及序列数量为6 667条,发生频率为9.90%。SSR序列中包括166种重复基元类型,其中,单核苷酸、二核苷酸和三核苷酸为优势重复类型,SSR位点数分别为3 890(41.36%),2 885(30.68%),2 489(26.46%)。SSR位点重复次数差别较大,其中重复10次比例最高,SSR位点数为1 621(17.24%),其次是5次(1 549,16.47%)和6次(1 495,15.90%)。主导重复基元类型为A/T(40.67%),AG/CT(27.59%),AAG/CTT(11.29%)。运用多态性分析的方式初步验证了SSR位点在楠木标记中的可行性。同时,利用Primer 3.0进行引物设计,随机筛选18对进行验证,9对引物可以扩增出预期大小的条带。【结论】通过对楠木高通量转录组序列的SSR信息的研究,获得6 667条SSR序列,主导重复类型为A/T,AG/CT和AAG/CTT。同时随机挑选18对引物对SSR分子标记方法进行验证。本文对楠木SSR标记的研究将有助于楠木基因挖掘、分子标记育种和资源保护。     

杜仲转录组微卫星特征分析

为给杜仲SSR分子标记的选择提供参考依据,利用MISA软件对杜仲转录组序列进行了微卫星查找,共获得重复单元长度为1~6碱基的29 067个微卫星序列。分析这些微卫星序列后发现,在杜仲转录组中,长度为二核苷酸的微卫星重复单元最为丰富,占43.7%;在单碱基重复和二碱基重复这两种类型中,最主要的优势重复单元分别是A/T和AG/CT;在三碱基、四碱基、五碱基重复类型中,(AAN)n、(ANAT)n、(AANAN)n分别为其对应的优势重复单元。杜仲转录组中,微卫星的频率随着微卫星总长度的增加而降低。     

基于桤木属转录组测序的SSR分子标记的开发

[目的]基于转录组数据开发适用于桤木属树种的SSR标记,揭示其在转录组序列中的分布类型及特征,为桤木属树种分子标记辅助育种提供有利工具。[方法]利用Micro SAtellite(MISA)软件对所有转录组序列进行SSRs搜索,并对SSR位点的数量、分布特征进行统计分析。设计100对SSR引物,采用琼脂糖凝胶电泳和毛细管电泳分离检测方法对3种不同倍性桤木属植物(12份材料)进行遗传多样性检测,确定引物多态性及通用性。[结果]85 769条Unigenes序列中发现8 678个SSR位点,分布在8 298条Unigenes中发生频率为9.67%,转录组序列中平均每14.04 kb长度就有一个SSR位点分布。其中,二核苷酸重复类型数量最多,占65.87%。根据转录组Unigenes序列,利用Primer 3软件共设计出4 531对符合要求的引物,挑选出的100对SSR引物中,获得18对多态性高、稳定性好的SSR引物。[结论]本研究可扩增出多态性位点的引物重复单元以二、三核苷酸重复为主。基于桤木属转录组序列的SSR标记开发是可行的,开发的引物为桤木属遗传多样性分析、分子育种、遗传图谱构建和功能基因的挖掘提供了丰富的候选分子标记。     

基于苦楝转录组测序的SSR分子标记开发

【目的】基于转录组测序数据开发苦楝SSR标记,为苦楝种质资源的选育、评价和遗传改良提供理论基础和科学依据。【方法】利用Illumina Hi Seq~(TM)2500平台对苦楝叶片进行高通量测序,分别利用Trinity、MISA软件对转录组数据进行拼接组装、序列检索及SSR分布类型和特征分析。采用Primer 3软件设计合成100对SSR引物,分别用2%琼脂糖凝胶电泳和PAGE凝胶电泳进行初步筛选和多态性筛选。【结果】共获得Unigenes 20 077条,平均长度为1 431.82 bp,N50为1 955 bp;对Unigenes进行SSR位点的发掘和分析,发现4 116条Unigenes序列中含有5 469个SSR位点,SSR的发生频率为20.50%,平均分布距离为8.74 kb;单、二、三核苷酸重复单元分别占总SSR的50.23%、24.43%和22.11%; SSR重复单元以6～10次重复的最多,占总数的51.27%。筛选出16对多态性引物并进行PCR扩增,共检测到56个等位基因片段,平均每个位点含3.5个,平均有效等位基因数为2.31个,平均Shannon多样性指数I为0.861,平均多态性信息含量PIC为0.507。【结论】苦楝转录组SSR位点具有较高的出现频率和分布密度,基元类型和重复次数相对较高,具有高多态性的潜能,可以进行有目的的引物设计和开发。开发的16对SSR引物在苦楝上可检测到较高的多态信息含量,进一步丰富了楝科现有SSR标记数据库资源,可为楝科植物在基因组水平上的遗传多样性分析和分子辅助育种等提供参考。     

基于马尾松干旱转录组的抗旱功能SSR位点分析

[目的]了解马尾松干旱胁迫转录组序列的功能分布和SSR位点分布特征,并探索与干旱关联的SSR位点。[方法]对马尾松幼苗进行持续干旱胁迫,选取干旱10、15、25 d及正常供水对照的马尾松针叶样品,通过总RNA提取、Illumina测序、raw reads去冗、Trinity拼接获得unigenes。利用Blast比对,对Unigene序列进行GO、KOG、KEGG功能注释及分类;利用Misa软件进行SSR位点批量搜寻,primer 3. 0软件进行规模化SSR引物设计;利用GOSeq(1. 10. 0)、KOBAS(v2. 0. 12)软件对差异表达的含SSR位点Unigene进行GO、KEGG显著性富集分析。[结果]马尾松转录组194 821个Unigene中,101 806个Unigene获得注释,包含64 943个GO功能注释、35 880个KOG功能注释及30 882个KEGG注释。搜寻到6 728个SSR位点,分布于6 367个Unigene中,SSR出现频率为3. 45%;重复类型以单、三、二核苷酸为主,分别占总SSR的35. 82%、33. 03%和25. 22%;重复基序以A/T、AT/AT、AG/CT、AGC/CTG、AAG/CTT为主;基序长度以10 20 bp的短序列SSR为主;基序重复次数以5～10次重复占优势;批量设计13 338对SSR引物。422个含SSR位点Unigene具有差异表达; KEGG富集分析发现,有11个含SSR位点的差异Unigene参与了光合作用、植物激素信号传导及类胡萝卜素合成等3个与干旱响应相关的代谢途径。[结论]从较高质量的马尾松转录组中获得101 806个具有注释的Unigene;从6 367个Unigene中挖掘出6 728个SSR位点; 422个含SSR位点的差异表达Unigene中筛选出11个与干旱关联的SSR功能位点,为抗旱功能基因定位及马尾松抗旱分子机制研究奠定基础。     

蕨类植物芒萁幼孢子体转录组高通量测序及特征分析

[目的]采用高通量测序技术Illumina Mi Seq250获得蕨类植物芒萁的孢子体转录组数据,以期为芒萁的生长、发育、代谢调控、微进化机制分析等提供重要的分子信息。[方法]应用生物信息学方法对测序获得的大量单基因簇(Unigene)进行基因功能注释、代谢途径及微卫星分析等。[结果]本研究共获得18 463 296条序列读取片段(reads),总碱基数为4.62Gbp序列信息,经序列组装最终得到63 169个Unigene,平均单条Unigene长度为863 bp,N50为1 587 bp,其中分布在200 500 bp长度区间Unigene占总数的55.4%。数据库中的序列同源性比较表明,26 826个Unigene与其他物种的已知基因具有不同程度的同源性。芒萁转录组中的Unigene根据GO功能大致可分为细胞组成、分子功能和生物学过程3大类47个分支,其中有大量的Unigene与细胞进程、绑定活性、代谢过程和催化活性相关。将Unigene与COG数据库进行比对,根据其功能大致可分为26类。以KEGG数据库作为参考,依据代谢途径可将Unigene定位到257个代谢途径分支。此外,利用MISA软件检索2 6碱基微卫星,共找到13 286个SSR。在不同长度微卫星中,三核苷重复数量最多,占总数的40.41%。在各重复基序类型中出现频率最高的为AG/CT(14.45%)与AAG/CTT(12.39%)。利用重复基序开发的多态性SSR标记,可应用于芒萁不同个体的基因型分型鉴定。[结论]本研究获得了较高质量的芒萁转录组数据库,揭示了芒萁孢子体生长发育过程中表达基因的功能总体特征,可为芒萁进一步的功能基因挖掘和分子标记规模化开发奠定基础。     

白蜡虫雄虫真蛹转录组分析

[目的]蛹期对于白蜡虫雄虫发育和交配具有重要意义,但目前缺乏对于白蜡虫真蛹基因转录情况的分析,本研究旨在获得白蜡虫真蛹转录组数据,了解其转录组特征。[方法]采用Illumina Hi Seq 2000进行高通量测序,并进行生物信息学分析;采用荧光定量PCR(qRT-PCR)检测3个热激蛋白基因(hsp)在不同虫态中的表达动态。[结果]共获得63 957条unigene、63 272个开放阅读框(ORF)、9 561个SSR分子标记,unigene平均长度674 bp,共有14 327条unigene获得了注释,Gene Ontology(GO)注释和Kyoto Encyclopedia of Genes and Genomes(KEGG)注释分析表明,很多基因参与的功能与白蜡虫蛹期的生理特征吻合。[结论]该转录组数据为白蜡虫基因功能鉴定和蛋白组学分析研究奠定了数据基础。     

白蜡虫热激蛋白序列分析

[目的]了解白蜡虫泌蜡高峰期所表达的热激蛋白及其序列特征。[方法]从白蜡虫转录组数据库筛选获得hsp10、hsp20、hsp40、hsp60、hsp70、hsp90基因序列,进行序列联配和进化树分析,并对hsp40、hsp60、hsp90基因在不同温度下mRNA相对表达量情况进行荧光定量PCR检测。[结果]从白蜡虫转录组数据库中共筛选出32条hsp基因序列,序列分析发现它们与目前已报道的序列具有很高一致性,大部分含有保守区域和特征序列,能够与同目昆虫基因聚为一支。hsp40、hsp60、hsp90基因均能被温度诱导。[结论]在白蜡虫泌蜡高峰期共鉴定到32个热激蛋白非重复序列基因(unigene):hsp10 unigene 3个,hsp20 unigene 1个,hsp40 unigene 11个,hsp60 unigene 1个,hsp70 unigene 13个,hsp90 unigene 3个。     

泡泡刺高通量转录组鉴定及其黄酮类代谢途径初步分析

[目的]为了更好的了解泡泡刺黄酮类生物合成途径及其潜在的抗逆作用,[方法]应用高通量RNA-seq测序技术对泡泡刺当年生新叶进行了转录组测序及相关生物信息学分析。[结果]表明:泡泡刺当年生叶中获得了13 013 444 total reads,总核苷酸数为1 171 209 960 bp(1.09 Gb),组装拼接后得到48 921条unigene序列;Nr、Swiss-Prot、COG、GO、KEGG等数据库分析显示,有30 407个Nr注释,19 671个Swissprot注释,9 273个COG功能注释,46 153个GO功能注释,13 654个KEGG注释,并从KEGG通路中找到参与黄酮类化合物合成途径的关键基因片段186个。[结论]本研究获得了较好的泡泡刺转录组序列信息,并且其中富含黄酮类化合物代谢途径的相关基因,这为其抗逆机制、药用价值的研究及开发利用提供宝贵资源。     

基于cpDNA rps16序列分析兰考泡桐与白花泡桐和毛泡桐的遗传关系

【目的】通过测序法分析兰考泡桐与白花泡桐和毛泡桐在叶绿体rps16序列上的遗传差异,旨在分析三者之间在叶绿体基因上的变化特点和规律,探讨其种间的遗传关系。【方法】选取兰考泡桐、白花泡桐和毛泡桐各15个样本,对其提取的DNA用PCR扩增获得特异片段,并将其纯化与测序。利用软件Clustal X 2.0对所得序列进行排序;运行MEGA 4软件,进行多序列比对,分析其序列特征,并计算出K2P遗传距离。【结果】(1)对获得的rps16序列进行测定分析,得兰考泡桐序列长度分别为932 933 bp;白花泡桐序列长度为932 bp;毛泡桐序列长度分别为916918 bp。对所得rps16序列进行排序后的长度为938 bp,平均GC含量为34.31%。3个种所代表的个体之间共有10个变异位点,占整个序列长度的1.07%。其中有9个变异位点属于碱基插入或缺失类型,占变异位点总数的90%,占整个序列长度的0.96%。有1个变异位点属于碱基替换类型,占整个变异位点总数的10%,占整个序列长度的0.11%。(2)整个rps16片段的序列共有10个变异位点,其中兰考泡桐与白花泡桐在总的变异位点上,具有一致的碱基位点9个,占总变异的90%。而兰考泡桐与毛泡桐相比,没有相同的碱基。【结论】根据三种泡桐的rps16序列的序列特征和变异位点的分析,表明在叶绿体遗传方面,兰考泡桐具有与白花泡桐更多相似的遗传物质,其亲缘关系较近。综上所述,推测兰考泡桐与白花泡桐可能来自同一母系遗传。     

红豆树优树子代遗传多样性及与生长相关性分析

[目的]研究红豆树优树自由授粉子代遗传多样性及其对生长的影响,比较天然居群子代和孤立木子代遗传多样性差异,揭示子代遗传多样性的变化规律及天然居群在子代遗传多样性维持中的作用,为红豆树遗传保育和优异种质挖掘提供科学依据。[方法]以来自浙、闽、赣、川等26个红豆树优树自由授粉家系为研究对象,利用11对SSR引物对765个子代群体进行遗传多样性评价,同时分析子代遗传多样性参数与种子、生长性状的相关性。[结果](1)红豆树优树子代群体具有较高的遗传多样性,有效等位基因数为7.766个,观测杂合度(H_O)和期望杂合度(H_E)分别为0.469和0.865。(2)除SSR8外,其余位点的观测杂合度均小于期望杂合度,表明子代群体绝大多数位点处于杂合子缺失状态。(3)红豆树不同家系的遗传多样性存在明显差异,12号家系的遗传多样性水平最高,8号家系则最低。(4)比较发现,天然居群子代的遗传多样性显著或极显著地高于孤立木子代。(5)F统计量和分子方差分析(AMOVA)均表明,红豆树优树子代群体的遗传变异主要存在于家系内,家系间的遗传分化相对较小。(6)相关性分析发现,子代遗传多样性参数与种子性状、子代年高生长量呈显著正相关(r=0.378~0.527)。[结论]较大的红豆树天然居群在维持其子代较高遗传多样性中发挥了重要作用,子代遗传多样性显著影响苗木生长,这为红豆树遗传保育和优良家系选择提供了理论依据。     

日本落叶松EST-SSR标记挖掘及特征分析

基于日本落叶松转录组测序(RNA-seq)结果,利用SSRIT工具查找SSR位点,检测到1 788个EST-SSR位点,平均长度是17.5 bp。进一步统计分析发现在落叶松中重复类型SSR位点出现次数最多的是六、三核苷酸,其余依次是五、二、四核苷酸重复类型;同时检测到656种重复基元类型,种类丰富。用Primer Premier 5.0设计500对引物,随机合成102对,以红豆杉科、柏科、松科材料为模板进行检测,结果显示在红豆杉科和柏科里的扩增率分别为95%、50%;在松科不同属的通用性为落叶松属99.01%、黄杉属67.4%、雪松属63.8%、冷杉属67.4%、云杉属70.2%和松属75.7%。利用101对SSR引物对7种不同落叶松品种亲缘关系进行了分析,获得79个多态性位点,在遗传相似系数0.69处将其区分为4个类群。