全 文 ：书能源草芦竹遗传多样性的犐犛犛犚分析
曾汉元１，２，魏麟２，刘鹏２，刘选明３
（１．湖南农业大学生物科学技术学院，湖南 长沙４１０１２８；２．怀化学院生命科学系，湖南 怀化４１８００８；
３．湖南大学生命科学研究院，湖南 长沙４１００８２）
摘要：采用ＩＳＳＲ分子标记技术，对我国６省１１个芦竹种群共８９个个体的遗传多样性水平和遗传结构进行了研
究。从ＵＢＣ８０１～９００中筛选出７个有效引物，共扩增出４０个位点，其中多态性位点为３８个。ＰＣＲ产物的分子量
为１００～２０００ｂｐ。在物种水平上，Ｓｈａｎｎｏｎ多样性指数为０．４９４，Ｎｅｉ指数为０．３２９，遗传分化系数Ｇｓｔ为０．８２６，平
均基因流Ｎｍ 为０．０９７。种群内的平均遗传多样度Ｈｓ为０．０４０。表明芦竹在物种水平上具有较高的遗传多样性而
种群内的分化极小。ＵＰＧＭＡ法聚类结果表明，芦竹１１个种群分为３个类群：浙江和江苏的芦竹聚为一簇，湖南、
贵州和广西的芦竹聚为一簇，云南种群单独为一支。不同种群芦竹亲缘关系的远近与它们的地理距离有一定的关
系，但不完全一致。
关键词：芦竹；ＩＳＳＲ；遗传多样性；种群
中图分类号：Ｑ９４３ 文献标识码：Ａ 文章编号：１００４５７５９（２０１３）０３０２６６０８
犇犗犐：１０．１１６８６／ｃｙｘｂ２０１３０３３５
芦竹 （犃狉狌狀犱狅犱狅狀犪狓）属于禾本科芦竹属，系多年生高大丛生草本，具有发达的根状茎，产于我国广东、海
南、广西、贵州、云南、四川、湖南、江西、福建、台湾、浙江和江苏等省，生于河岸、道旁、沙质壤土上，亚洲、非洲和大
洋洲热带地区广布［１］。它的用途广泛，根状茎入药，具有清热泻火、止呕生津的功效［２］；嫩枝叶是牲畜的良好青饲
料［１］，秆是优质造纸原料［３］和良好的管乐簧片材料，并可编织用具和制造人造丝，还是优良护堤植物［１］；某些生态
型的芦竹对镉、镍、铅和汞等重金属有良好的富集作用，可用于对重金属污染土壤的修复［４，５］；芦竹生物质产量
高，生态适应性强［６］，富含纤维素［７］，是生产燃料酒精的好原料，是欧洲国家筛选出的四大能源牧草之一，在欧洲
还建立了“芦竹网络”，旨在将芦竹作为生物质能源及造纸原料加以推广利用［８］。目前，木质纤维素燃烧转化为
热、气和电的技术正在普及［６］，有关芦竹木质纤维素转化为燃料酒精的物理和化学法预处理工艺以及纤维素乙醇
发酵技术研究成为热点［９１６］，预处理的方法有稀酸水解法、汽爆法和超临界液化提取法等，目的是为了分离和去
除木质纤维素原料中的木质素，从而提高纤维素和半纤维素的转化率。代莹等［１３］对分别采用１％Ｈ２ＳＯ４ 和１０％
ＮａＯＨ预处理的黑麦草（犔狅犾犻狌犿狆犲狉犲狀狀犲）酶解制取燃料酒精的最佳条件进行了研究，叶艳平等［１４］探讨了青霉纤
维素酶的酶解温度、酶解时间、ｐＨ 值和酶载量等因素对芦竹酶解效果的影响，当酶解时间为７２ｈ、酶解温度
４５℃、ｐＨ值４．８、酶载量３０ＦＰＵ／ｇＤＭ 时，葡萄糖产率为（６１．７±１．２）％、木糖产率为（４０±６．９）％，当增加β－
葡萄糖苷酶酶活至４０ＣＢＵ／ｇＤＭ 时，葡萄糖产率为６９％、木糖产率达９８％。Ｇａｌｅｔｔｉ等［１５］报道了芦竹酶解生产
乙酰丙酸和γ戊内酯的工艺。此外，对芦竹燃料酒精生产中的副产物半纤维素和木质素的利用研究也有报道，
可用于生产活性炭［１７］、抗氧化剂［１８］和木糖［１９］等高附加值产品，从而实现原料充分利用和产品价值最大化。
ＩＳＳＲ分子标记技术比ＲＦＬＰ和ＳＳＲ技术操作要简单，比 ＲＡＰＤ的稳定性和重复性要好，而且可揭示比
ＲＦＬＰ、ＲＡＰＤ和ＳＳＲ更多的多态性［２０］，因此，目前ＩＳＳＲ已在基因定位、遗传作图［２１］、系统与进化、遗传多样性
及品种鉴定［２２２４］等研究中被广泛应用，李飞飞等［２３］采用ＩＳＳＲ分子标记成功区分了黄花苜蓿（犕犲犱犻犮犪犵狅犳犪犾犮犪
狋犪）、多变苜蓿（Ｍｅｄｉｃａｇｏｖａｒｉａ）及紫花苜蓿（犕犲犱犻犮犪犵狅狊犪狋犻狏犪），刘欢等［２４］采用ＩＳＳＲ分子标记技术对来自中国、
加拿大、美国和欧洲的４２份饲用燕麦（犃狏犲狀犪狊犪狋犻狏犪）品种及６份裸燕麦（犃狏犲狀犪狊犪狋犻狏犪）品种进行了分析，所得结
２６６－２７３
２０１３年６月
草 业 学 报
ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ
第２２卷 第３期
Ｖｏｌ．２２，Ｎｏ．３
收稿日期：２０１２１１１２；改回日期：２０１２１２２０
基金项目：湖南省自然科学基金项目（０９ＪＪ６０４５）资助。
作者简介：曾汉元（１９６６），男，湖南绥宁人，教授，在读博士。Ｅｍａｉｌ：ｚｅｎｇｈａｎｙｕａｎ＠１６３．ｃｏｍ
通讯作者。Ｅｍａｉｌ：ｘｍｌ０５＠１２６．ｃｏｍ
果与形态学及其他标记分类结果一致，说明ＩＳＳＲ
分子标记应用于种内甚至近缘种间的分析都是可行
的。本研究采用ＩＳＳＲ分子标记技术对采自我国６
省共１１个种群的芦竹进行分析，旨在探讨它们的遗
传变异和亲缘关系，为芦竹种质资源的筛选和利用
提供理论依据。
１ 材料与方法
１．１ 材料
本试验材料芦竹于２００８年１月－２００９年５月
采集于野外，包括１１个种群共８９个样品，同一种群
不同个体的采样直线距离在２００ｍ以上。采样时
选取嫩叶，用硅胶干燥保存。材料来源见表１。
１．２ 方法
１．２．１ 总ＤＮＡ提取与检测 采用改良的ＣＴＡＢ
法［２５］进行提取。在提取的ＤＮＡ样品中加入５０μＬ
ＴＥ溶液（含有１０μｇ／ｍＬＲＮａｓｅＡ）。得到的ＤＮＡ
用０．８％的琼脂糖凝胶电泳和紫外分光光度计检测
浓度和纯度。保存于－２０℃冰箱中。
１．２．２ ＩＳＳＲ有效引物的筛选 根据加拿大哥伦比
亚大学公布的第九套序列，选择编号为ＵＢＣ８０１～
９００的引物，由上海生工公司合成。在８９份样品
中，随机选取４个样品，分别用１００个引物进行
ＰＣＲ扩增。结果筛选出８０７，８０８，８１０，８１２，８１７，８３４
和８６８共７个引物为有效引物，扩增出具较高多态
性的清晰条带。然后用这７个引物对所有样品进行
ＰＣＲ扩增。有效引物的编号及其序列见表２。
１．２．３ ＩＳＳＲ试验条件的优化 考虑到 Ｍｇ２＋ 浓
度、ｄＮＴＰ浓度、模板ＤＮＡ量及退火温度等均会影
响ＩＳＳＲ扩增结果，因此，需要对试验条件进行优化
研究。
ＩＳＳＲ－ＰＣＲ反应体系的筛选：扩增前，将各样
品ＤＮＡ稀释成２０ｎｇ／μＬ，分别采用３个不同体系
进行对照试验（表３）。
结果以反应体系２的扩增结果最好，因此，后面
的试验全部采用反应体系２进行。
引物退火温度：参照引物的Ｔｍ值，本研究设置
了４８，４９，５０，５１和５２℃共５个退火温度进行扩增。
结果表明，引物８０７和８１７退火温度以５２℃扩增效
果最好，而其他引物则以退火温度５０℃时扩增效果
最好，因此，除引物８０７和８１７的退火温度确定为
５２℃外，其他的引物退火温度均为５０℃。
表１ 不同种群芦竹材料说明
犜犪犫犾犲１ 犜犺犲犱犲狊犮狉犻狆狋犻狅狀狅犳犃．犱狅狀犪狓犻狀犱犻犳犳犲狉犲狀狋狆狅狆狌犾犪狋犻狅狀狊
种群号
Ｎｏ．ｏｆ
ｐｏｐｕｌａｔｉｏｎ
采集地
Ｌｏｃａｌｉｔｙ
所采集的样品份数
Ｓａｍｐｌｅｎｕｍｂｅｒｓ
ｏｆｃｏｌｅｃｔｅｄ
样品编号
Ｎｏ．ｏｆ
ｓａｍｐｌｅ
１ 江苏阜宁Ｆｕｎｉｎｇ，Ｊｉａｎｇｓｕ ８ １～８
２ 江苏南京Ｎａｎｊｉｎｇ，Ｊｉａｎｇｓｕ ９ ９～１７
３ 浙江嵊州Ｓｈｅｎｇｚｈｏｕ，Ｚｈｅｊｉａｎｇ １０ １８～２７
４ 浙江慈溪Ｃｉｘｉ，Ｚｈｅｊｉａｎｇ １０ ２８～３７
５ 湖南衡阳 Ｈａｎｙａｎｇ，Ｈｕｎａｎ ８ ３８～４５
６ 湖南凤凰Ｆｅｎｇｈｕａｎｇ，Ｈｕｎａｎ １０ ４６～５５
７ 贵州铜仁Ｔｏｎｇｒｅｎ，Ｇｕｉｚｈｏｕ １０ ５６～６５
８ 浙江海盐 Ｈａｉｙａｎ，Ｚｈｅｊｉａｎｇ ４ ６６～６９
９ 湖南通道Ｔｏｎｇｄａｏ，Ｈｕｎａｎ ８ ７０～７７
１０ 广西桂林Ｇｕｉｌｉｎ，Ｇｕａｎｇｘｉ ８ ７８～８５
１１ 云南洱源Ｅｒｙａｎ，Ｙｕｎｎａｎ ４ ８６～８９
表２ 试验筛选出的犐犛犛犚有效引物及其序列
犜犪犫犾犲２ 犐犛犛犚犲犳犳犻犮犻犲狀狋狆狉犻犿犲狉狊犪狀犱狋犺犲犻狉
狊犲狇狌犲狀犮犲狊狌狊犲犱犻狀狋犺犻狊犲狓狆犲狉犻犿犲狀狋
引物编号ＰｒｉｍｅｒＮｏ． 序列Ｓｅｑｕｅｎｃｅ（５′－３′）
ＵＢＣ８０７ ＡＧＡＧＡＧＡＧＡＧＡＧＡＧＡＧＴ
ＵＢＣ８０８ ＡＧＡＧＡＧＡＧＡＧＡＧＡＧＡＧＣ
ＵＢＣ８１０ ＧＡＧＡＧＡＧＡＧＡＧＡＧＡＧＡＴ
ＵＢＣ８１２ ＧＡＧＡＧＡＧＡＧＡＧＡＧＡＧＡＡ
ＵＢＣ８１７ ＣＡＣＡＣＡＣＡＣＡＣＡＣＡＣＡＡ
ＵＢＣ８３４ ＡＧＡＧＡＧＡＧＡＧＡＧＡＧＡＧＹＴ
ＵＢＣ８６８ ＧＡＡＧＡＡＧＡＡＧＡＡＧＡＡＧＡＡ
表３ 犐犛犛犚－犘犆犚反应体系
犜犪犫犾犲３ 犐犛犛犚－犘犆犚狉犲犪犮狋犻狅狀狊狔狊狋犲犿 μＬ
项目
Ｉｔｅｍ
体系１
Ｓｙｓｔｅｍ１
体系２
Ｓｙｓｔｅｍ２
体系３
Ｓｙｓｔｅｍ３
ｄｄＨ２Ｏ １１ １１ １１
１０×Ｂｕｆｆｅｒ ２．４ ２．２ ２．２
２５ｍｍｏｌ／ＬＭｇＣｌ２ １．２ １．０ ０．８
１０μｍｏｌ／Ｌ引物Ｐｒｉｍｅｒ ０．７ ０．６ ０．７
ｄＮＴＰｓ ０．５ ０．５ ０．４
１Ｕ／μＬＴａｑＤＮＡ聚合酶
ＴａｑＤＮＡｐｏｌｙｍｅｒａｓｅ
０．２ ０．２ ０．２
总体积Ｔｏｔａｌｖｏｌｕｍｅ（μＬ） １６ １５．５ １５．３
１０×Ｂｕｆｆｅｒ：１０ｍｍｏｌ／ＬＴｒｉｓ－ＨＣｌ，ｐＨ９．０，５０ｍｍｏｌ／ＬＫＣｌ，０．１％
ＴｒｉｔｏｎＸ－１００．
７６２第２２卷第３期 草业学报２０１３年
１．２．４ ＩＳＳＲ扩增及产物检测 扩增反应在ＰＴＣ２００ＰＣＲ扩增仪（美国）上进行，扩增步骤及参数如下：预变性
９４℃３ｍｉｎ，变性９４℃４５ｓ，退火Ｔｍ４５ｓ（退火温度依据引物变化，并根据所得结果进行调整），延伸７２℃１ｍｉｎ，
３５个循环，后延伸７２℃５ｍｉｎ。扩增完毕，将样本取出，置于４℃冰箱保存，等待电泳。产物检测：１％琼脂糖凝胶
在１×ＴＡＥ电泳缓冲液中电泳，电压：１５０Ｖ，电泳时间为３ｈ。采用ＤＬ２０００ＤＮＡ Ｍａｒｋｅｒ标记：２０００，１０００，
７５０，５００，２５０和１００ｂｐ。电泳后，将凝胶ＤＮＡ在含溴化乙锭的水溶液中浸泡染色１５ｍｉｎ，在凝胶成像分析系统
中观察并照相。
１．２．５ 数据分析方法 按照同一位置上扩增产生条带的有无进行统计，有带的（含弱带）记为“１”，无带的记为
“０”，缺失的记为－，记录电泳带谱。利用Ｐｏｐｇｅｎ３２软件计算多态位点数、多态位点比率、Ｓｈａｎｎｏｎ多样性指数
（Ｉ）、Ｎｅｉ遗传多样性指数（Ｈ）、种群间基因分化系数（Ｇｓｔ）、基因流（Ｎｍ）、种群间的遗传一致度和遗传距离。根据
Ｎｅｉ’ｓ遗传距离按不加权成对群算术平均法 （ＵＰＧＭＡ）进行聚类分析。
２ 结果与分析
２．１ ＩＳＳＲ扩增的多态位点分析
７个有效引物对１１个芦竹种群ＤＮＡ进行扩增，共扩增出４０个位点，其中多态位点有３８个，多态位点比率
为９５％（表４）。ＰＣＲ产物的分子量在１００～２０００ｂｐ，不同引物扩增出的多态性位点数从４个到７个不等，形成
了带形丰富、片断大小及其组合不同的电泳图谱（图１～图４），这说明了芦竹在物种水平上具有较高的遗传多样
性。
图１ 引物犝犅犆８０７在种群１～５中的扩增结果
犉犻犵．１ 犐犛犛犚犪犿狆犾犻犳犻犮犪狋犻狅狀狉犲狊狌犾狋狊狅犳狆狉犻犿犲狉狊犝犅犆８０７犻狀狆狅狆狌犾犪狋犻狅狀狊狅犳１－５
Ｍ：标记 Ｍａｒｋｅｒ；１～４５为样品编号，与表１中的相同。１－４５ａｒｅｔｈｅｓａｍｅａｓＴａｂｌｅ１．
图２ 引物犝犅犆８０７在种群６－１１中的扩增结果
犉犻犵．２ 犐犛犛犚犪犿狆犾犻犳犻犮犪狋犻狅狀狉犲狊狌犾狋狊狅犳狆狉犻犿犲狉狊犝犅犆８０７犻狀狆狅狆狌犾犪狋犻狅狀狊狅犳６－１１
Ｍ：标记 Ｍａｒｋｅｒ；４６～８９为样品编号，与表１中的相同。４６－８９ａｒｅｔｈｅｓａｍｅａｓＴａｂｌｅ１．
８６２ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．３
图３ 引物犝犅犆８１０在种群１～５中的扩增结果
犉犻犵．３ 犐犛犛犚犪犿狆犾犻犳犻犮犪狋犻狅狀狉犲狊狌犾狋狊狅犳狆狉犻犿犲狉狊犝犅犆８１０犻狀狆狅狆狌犾犪狋犻狅狀狊狅犳１－５
Ｍ：标记 Ｍａｒｋｅｒ；１～４５为样品编号，与表１中的相同。１－４５ａｒｅｔｈｅｓａｍｅａｓＴａｂｌｅ１．
图４ 引物犝犅犆８１０在种群６～１１中的扩增结果
犉犻犵．４ 犐犛犛犚犪犿狆犾犻犳犻犮犪狋犻狅狀狉犲狊狌犾狋狊狅犳狆狉犻犿犲狉狊犝犅犆８１０犻狀狆狅狆狌犾犪狋犻狅狀狊狅犳６－１１
Ｍ：标记 Ｍａｒｋｅｒ；４６～８９为样品编号，与表１中的相同。４６－８９ａｒｅｔｈｅｓａｍｅａｓＴａｂｌｅ１．
２．２ 芦竹的遗传变异和遗传分化
芦竹在物种水平上的Ｓｈａｎｎｏｎ指数为０．４９４，Ｎｅｉ指数为０．３２９。各种群内的Ｓｈａｎｎｏｎ指数（Ｉ）的多态性为
０～０．０９６，Ｎｅｉ指数的多态性为０～０．０６６（表４）。因此，芦竹在物种水平上的遗传多样性指数相对较高，但种群
内的遗传多样性指数都相当低。
芦竹种群的遗传分化情况见表５。１１个种群总的遗传多样度 Ｈｔ为０．３４１，种群间遗传分化系数 Ｇｓｔ为
０．８２６，即８２．６％的遗传变异存在于种群间，１７．４４％的遗传变异存在于种群内。种群内遗传多样度Ｈｓ为０．０４０，
平均基因流为０．０９７。因此，芦竹种群间存在较大的差异，但种群内个体的差异极小。
引物不同，其遗传分化系数所占的比率不同，其中引物ＵＢＣ８０８的最高，占到９２．３％，引物ＵＢＣ８３４的最低，
为６６．０％（表５）。
２．３ 芦竹不同种群的遗传距离与聚类分析
芦竹不同种群间遗传一致度和遗传距离的计算结果见表６。
浙江嵊州芦竹与浙江慈溪芦竹的遗传距离最小（０．０２６），遗传一致度最大（０．９７４）；湖南通道芦竹与浙江海盐
芦竹间的遗传距离最大（０．７４３），遗传一致度最小（０．４７６）。
根据种群间的遗传距离，采用ＵＰＧＭＡ法对１１个种群进行聚类绘图（图５），结果大体分成三类群，浙江和江
苏的芦竹种群为一个类群，湖南、贵州和广西的芦竹种群为一个类群，云南种群单独为一支。因此，它们的亲缘关
系的远近与地理距离远近有关，但不完全相同。
９６２第２２卷第３期 草业学报２０１３年
表４ 芦竹种群的遗传多样性
犜犪犫犾犲４ 犌犲狀犲狋犻犮犱犻狏犲狉狊犻狋狔犪狀犪犾狔狊犻狊狅犳１１狆狅狆狌犾犪狋犻狅狀狊狅犳犃．犱狅狀犪狓
种群
Ｐｏｐｕｌａｔｉｏｎｓ
样本数
Ｓｉｍｐｌｅｎｕｍｂｅｒ
位点数
Ｔｏｔａｌｌｏｃｉ
多态位点数
Ｐｏｌｙｍｏｒｐｈｉｃｌｏｃｉ
多态位点比率Ｐｒｏｐｏｒｔｉｏｎｏｆ
ｐｏｌｙｍｏｒｐｈｉｃｌｏｃｉ（％）
Ｓｈａｎｎｏｎ指数
Ｓｈａｎｎｏｎ’ｉｎｄｅｘ（Ｉ）
Ｎｅｉ指数
Ｎｅｉ’ｉｎｄｅｘ（Ｈ）
１ ８ ４０ ３ ７．５ ０．０５０ ０．０３５
２ ９ ４０ ３ ７．５ ０．０４８ ０．０３３
３ １０ ４０ ３ ７．５ ０．０４７ ０．０３２
４ １０ ４０ ５ １２．５ ０．０８３ ０．０５９
５ ８ ４０ ４ １０．０ ０．０６４ ０．０４５
６ １０ ４０ ６ １５．０ ０．０９６ ０．０６６
７ １０ ４０ ４ １０．０ ０．０９５ ０．０６６
８ ４ ４０ ０ ０ ０ ０
９ ８ ４０ ３ ７．５ ０．０５１ ０．０３６
１０ ８ ４０ ４ １０．０ ０．０６４ ０．０４５
１１ ４ ４０ ４ １０．０ ０．０６０ ０．０４１
物种水平Ｓｐｅｃｉｅｓｌｅｖｅｌ ８９ ４０ ３８ ９５．０ ０．４９４ ０．３２９
表５ 芦竹１１个种群的遗传分化
犜犪犫犾犲５ 犌犲狀犲狋犻犮犱犻犳犳犲狉犲狀狋犻犪狋犻狅狀狊犪犿狅狀犵１１狆狅狆狌犾犪狋犻狅狀狊狅犳犃．犱狅狀犪狓
引物
Ｐｒｉｍｅｒｓ
样点数
Ｎｕｍｂｅｒｓｏｆ
ｓａｍｐｌｅｓｉｚｅ
总遗传多样度
Ｔｏｔａｌｇｅｎｅｄｉｖｅｒｓｉｔｙ
（Ｈｔ）
种群内遗传多样度
Ｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙｗｉｔｈｉｎ
ｐｏｐｕｌａｔｉｏｎ（Ｈｓ）
种群间遗传分化系数
Ｇｅｎｅｔｉｃｄｉｆｆｅｒｅｎｔｉａｔｉｏｎａｍｏｎｇ
ｐｏｐｕｌａｔｉｏｎ（Ｇｓｔ）
基因流
Ｇｅｎｅｆｌｏｗ
（Ｎｍ）
ＵＢＣ８０７ ８９ ０．４３１ ０．０６２ ０．８６８ ０．０８３
ＵＢＣ８０８ ８９ ０．４０５ ０．０３４ ０．９２３ ０．０５２
ＵＢＣ８１０ ８９ ０．２８４ ０．０２８ ０．７３９ ０．０６８
ＵＢＣ８１２ ８９ ０．２９８ ０．０２８ ０．８９３ ０．０７９
ＵＢＣ８１７ ８９ ０．４１０ ０．０５０ ０．８６５ ０．１２５
ＵＢＣ８３４ ８９ ０．２９７ ０．０４３ ０．６６０ ０．１３８
ＵＢＣ８６８ ８９ ０．２６３ ０．０３６ ０．８３２ ０．１３６
平均 Ｍｅａｎ ８９ ０．３４１ ０．０４０ ０．８２６ ０．０９７
标准差Ｓｔａｎｄａｒｄｄｅｖｉａｔｉｏｎ ０．０７１ ０．０１３ ０．０９４ ０．０３５
表６ 芦竹不同种群间的遗传一致度和遗传距离
犜犪犫犾犲６ 犌犲狀犲狋犻犮犻犱犲狀狋犻狋狔犪狀犱犱犻狊狋犪狀犮犲犫犲狋狑犲犲狀犱犻犳犳犲狉犲狀狋狆狅狆狌犾犪狋犻狅狀狊狅犳犃．犱狅狀犪狓
种群编号ＰｏｐｕｌａｔｉｏｎＮｏ． １ ２ ３ ４ ５ ６ ７ ８ ９ １０ １１
１ ０．７６４ ０．６７１ ０．６６５ ０．５７１ ０．６６０ ０．６６８ ０．５７３ ０．６９８ ０．７７９ ０．５７６
２ ０．２６９ ０．７７４ ０．７７２ ０．５６１ ０．５８７ ０．６０７ ０．７２１ ０．５６８ ０．７１８ ０．５９６
３ ０．３９８ ０．２５６ ０．９７４ ０．６４９ ０．６１２ ０．６４８ ０．７９７ ０．５６２ ０．６７９ ０．６０３
４ ０．４０８ ０．２５９ ０．０２６ ０．６３３ ０．５９４ ０．６３２ ０．７９９ ０．５２６ ０．６４９ ０．５８０
５ ０．５６１ ０．５７８ ０．４３３ ０．４５８ ０．６９３ ０．８１４ ０．７９３ ０．６２６ ０．６６６ ０．５１１
６ ０．４１６ ０．５３３ ０．４９２ ０．５２０ ０．３６７ ０．８４７ ０．６５４ ０．５６５ ０．７０２ ０．６２２
７ ０．４０３ ０．４９９ ０．４３４ ０．４５９ ０．２０６ ０．１６６ ０．７９２ ０．６６３ ０．７２５ ０．６３２
８ ０．５５６ ０．３２７ ０．２２７ ０．２２４ ０．２３３ ０．４２４ ０．２３３ ０．４７６ ０．５２１ ０．５８５
９ ０．３５９ ０．５６５ ０．５７８ ０．６４３ ０．４６８ ０．５７０ ０．４１１ ０．７４３ ０．８４６ ０．５５４
１０ ０．２５０ ０．３３２ ０．３８７ ０．４３２ ０．４０６ ０．３５４ ０．３２１ ０．６５３ ０．１６８ ０．５９７
１１ ０．５５２ ０．５１８ ０．５０５ ０．５４６ ０．６７１ ０．４７５ ０．４５９ ０．５３６ ０．５９１ ０．５１６
注：上三角是遗传一致度，下三角是遗传距离。
Ｎｏｔｅ：Ｎｅｉ’ｓｇｅｎｅｔｉｃｉｄｅｎｔｉｔｙ（ａｂｏｖｅｄｉａｇｏｎａｌ）ａｎｄｇｅｎｅｔｉｃｄｉｓｔａｎｃｅ（ｂｅｌｏｗｄｉａｇｏｎａｌ）．
０７２ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．３
３ 讨论
图５ 芦竹１１个种群的犝犘犌犕犃遗传聚类图
犉犻犵．５ 犝犘犌犕犃犱犲狀犱狉狅犵狉犪犿狅犳１１狆狅狆狌犾犪狋犻狅狀狊狅犳犃．犱狅狀犪狓
３．１ 芦竹种群内遗传多样性极低的原因
Ｂａｌｏｇｈ等［２６］在研究芦竹雌、雄配子发育中观察到
大约９０％的子房不发育，只有１０％左右的子房能够发
育，形成颖果。在这些能够发育的子房壁内观察到１
个倒生的胚珠，具双珠被，内含大量的细胞壁厚的珠心
细胞，大小均匀，偶见孢原细胞，但这些孢原细胞不能
发育出大孢子母细胞，因此不能产生卵细胞；在花粉发
育中，不能形成生殖核，因此，这样发育出来的颖果实
质上只有果皮，没有种子，是不育的。将这样的颖果放
在适宜的条件下催芽，未见萌发，也证明了这一点。
Ｈａｄｄａｄｃｈｉ等［２７］的研究表明，芦竹的体细胞染色体为
２ｎ＝７Ｘ＝８４，为不平衡的倍性水平，也是导致不育的
原因。因此，芦竹为克隆植物。芦竹采用根状茎和茎
秆进行无性繁殖，因而几乎没有基因交流，导致芦竹种群内的遗传多样性指数极低。再者，在芦竹采样时，同一地
区的采样间距不是很大，可能会导致所采集的部分个体来源于相同的基株，这也是导致所测定的同一采样地区的
芦竹遗传多样性指数低的原因。Ａｈｍａｄ等［２８］采用序列相关扩增多态性（ｓｅｑｕｅｎｃｅｒｅｌａｔｅｄａｍｐｌｉｆｉｃａｔｉｏｎｐｏｌｙ
ｍｏｒｐｈｉｓｍ，ＳＲＡＰ）和转座因子ＴＥ（ｔｒａｎｓｐｏｓａｂｌｅｅｌｅｍｅｎｔ）分子标志，对美国和法国的外来物种芦竹进行了研
究，发现法国的４个芦竹样品有一个多态性位点，而美国的所有芦竹样品的ＤＮＡ指纹图谱是一致的，属于一个
芦竹基株的克隆。这一研究结果也说明，同一个基株经过克隆繁殖产生的众多个体在较短的历史时期内几乎没
有遗传差异，因而种群内的遗传多样性极低。
３．２ 芦竹种群间遗传多样性较高的原因
杨允菲等［２９］和Ｂａｒｓｏｕｍ［３０］对部分克隆植物繁殖方式的研究表明，它们的有性繁殖比率极低，主要靠克隆繁
殖来更新种群。克隆繁殖的缺点是子代扩展距离有限，不利于种群间基因传递，减少了生物多样性出现的可能
性［３１３３］。然而，克隆植物芦竹却为亚洲、非洲和大洋洲热带地区广布种［１］，并且在物种水平上具有较高的遗传多
样性，Ｈａｄｄａｄｃｈｉ等［２７］对澳大利亚５８个芦竹种群的ＩＳＳＲ分析得出的Ｓｉｍｐｓｏｎ指数为０．８８１～１．０００，种群间的
遗传分化系数为０．４８５～１．０００；笔者对１１个芦竹种群的ＩＳＳＲ分析得出芦竹在物种水平上的平均Ｓｈａｎｎｏｎ指数
为０．４９４，种群间的遗传分化系数为０．６６０～０．９２３，都证明了这一点。其可能原因是：１）芦竹是起源于白垩纪的
古老物种，在起源与进化早期，它的有性生殖能力强，从而把它的后代扩展到世界上多个热带地区。由于进行有
性生殖会引起种群间的基因交流，因而使得今天分布在不同地区的芦竹具有较高的遗传多样性。但后来因为选
择了根状茎和茎秆的克隆繁殖策略，有性生殖能力逐渐丧失。２）芦竹有多种生态类型［３４３６］，可以在多种不同生境
中进化并积累遗传变异，从而使物种水平具有较高的遗传多样性。Ｈａｄｄａｄｃｈｉ等［２７］对澳大利亚５８个芦竹种群
的研究中发现３８种遗传型，认为产生如此丰富的遗传型的原因是：澳大利亚的芦竹可能来自多个基株或者来自
体细胞的突变或（和）染色体的多倍化。中国的芦竹遗传多样性较高可能也存在类似的原因。
３．３ 芦竹核心种质构建的取样方法的思考
笔者本次只采集了我国６省１１个地区的芦竹样品进行遗传多样性分析，而芦竹为亚洲、非洲和大洋洲热带
地区广布种，因此，今后有必要补充其他地区的芦竹样品。另外，研究表明，不同地区的芦竹种群间的遗传差异较
大，而地理距离较近的样品间的遗传差异极小，同一基株发展起来的芦竹个体间无遗传差异。因此，今后采样时
要扩大地理间距。詹世雄和曾宪威［３７］总结了构建核心种质的取样策略，有恒量取样、比例取样、对数取样和遗传
多样性取样法等。笔者认为，在芦竹取样时，要注意采集地理间距远的和各种生态型的芦竹样品，要综合考虑各
样品的植物学性状、农艺学性状和遗传多样性，通过聚类分析和主成分分析，从而构建芦竹的核心种质资源库。
１７２第２２卷第３期 草业学报２０１３年
参考文献：
［１］ 中国植物志编委会．中国植物志（第九卷第二分册）［Ｍ］．北京：科学出版社，２００２：２０２１．
［２］ 国家中医药管理局．中华本草［Ｍ］．上海：上海科学技术出版社，１９９９：６４５９６４５０．
［３］ 吴武汉．芦竹，一种高产优质的造纸原料［Ｊ］．天津造纸，１９９３，（４）：２８２９．
［４］ ＰａｐａｚｏｇｌｏｕＥＧ．犃狉狌狀犱狅犱狅狀犪狓Ｌ．ｓｔｒｅｓｓｔｏｌｅｒａｎｃｅｕｎｄｅｒｉｒｒｉｇａｔｉｏｎｗｉｔｈｈｅａｖｙｍｅｔａｌａｑｕｅｏｕｓｓｏｌｕｔｉｏｎｓ［Ｊ］．Ｄｅｓａｌｉｎａｔｉｏｎ，
２００７，２１１：３０４３１３．
［５］ 韩志萍，胡正海．芦竹对不同重金属耐性的研究［Ｊ］．应用生态学报，２００６，１６（１）：１６１１６５．
［６］ 余醉，李建龙，李高扬．芦竹作为清洁生物质能源牧草开发的潜力分析［Ｊ］．草业科学，２００９，２６（６）：６２６９．
［７］ 曾汉元，杨洋，姚元枝，等．不同居群芦竹纤维素和木质素含量的比较研究［Ｊ］．中国农学通报，２０１２，２８（１９）：２２５２２８．
［８］ ＬｅｗａｎｄｏｗｓｋｉＩ，ＳｃｕｒｌｏｃｋＪＭＯ，ＬｉｎｄｖａｌＥ，犲狋犪犾．Ｔｈｅｄｅｖｅｌｏｐｍｅｎｔａｎｄｃｕｒｒｅｎｔｓｔａｔｕｓｏｆｐｅｒｅｎｎｉａｌｒｈｉｚｏｍａｔｏｕｓｇｒａｓｓｅｓａｓ
ｅｎｅｒｇｙｃｒｏｐｓｉｎｔｈｅＵＳａｎｄＥｕｒｏｐｅ［Ｊ］．ＢｉｏｍａｓｓａｎｄＢｉｏｅｎｅｒｇｙ，２００３，２５（４）：３３５３６１．
［９］ 商成祥．国内外纤维素燃料生物汽油发展综述［Ｊ］．化学工程与装备，２００９，（１１）：１３１１３２．
［１０］ Ｊｅｇｕｉｒｉｍ Ｍ，ＴｒｏｕｖéＧ．Ｐｙｒｏｌｙｓｉｓｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎｄｋｉｎｅｔｉｃｓｏｆ犃狉狌狀犱狅犱狅狀犪狓ｕｓｉｎｇｔｈｅｒｍｏｇｒａｖｉｍｅｔｒｉｃａｎａｌｙｓｉｓ［Ｊ］．Ｂｉｏｒｅ
ｓｏｕｒｃｅＴｅｃｈｎｏｌｏｇｙ，２００９，１００（１７）：４０２６４０３１．
［１１］ ＳｃｏｒｄｉａＤ，ＣｏｓｅｎｔｉｎｏＳＬ，ＬｅｅＪＷ，犲狋犪犾．Ｄｉｌｕｔｅｏｘａｌｉｃａｃｉｄｐｒｅｔｒｅａｔｍｅｎｔｆｏｒｂｉｏｒｅｆｉｎｉｎｇｇｉａｎｔｒｅｅｄ（犃狉狌狀犱狅犱狅狀犪狓Ｌ．）［Ｊ］．
ＢｉｏｍａｓｓａｎｄＢｉｏｅｎｅｒｇｙ，２０１１，３５（７）：３０１８３０２４．
［１２］ ＡｙｓｕＴ，ＫüüｋＭ Ｍ．Ｌｉｑｕｅｆａｃｔｉｏｎｏｆｇｉａｎｔｒｅｅｄ（犃狉狌狀犱狅犱狅狀犪狓Ｌ．）ｂｙｓｕｐｅｒｃｒｉｔｉｃａｌｆｌｕｉｄｅｘｔｒａｃｔｉｏｎ［Ｊ］．Ｆｕｅｌ，２０１３，１０３：
７５８７６３．
［１３］ 代莹，郭和蓉，张兴龙，等．一年生黑麦草制取燃料酒精发酵工艺的初步研究［Ｊ］．草业学报，２０１２，２１（１）：２２０２２８．
［１４］ 叶艳平，周玉杰，丁一刚，等．青霉纤维素酶酶解芦竹的研究［Ｊ］．化学与生物工程，２０１２，２９（５）：２２２５．
［１５］ ＧａｌｅｔｔｉＡＭＲ，ＡｎｔｏｎｅｔｔｉＣ，ＲｉｂｅｃｈｉｎｉＥ，犲狋犪犾．Ｆｒｏｍｇｉａｎｔｒｅｅｄｔｏｌｅｖｕｌｉｎｉｃａｃｉｄａｎｄｇａｍｍａｖａｌｅｒｏｌａｃｔｏｎｅ：Ａｈｉｇｈｙｉｅｌｄ
ｃａｔａｌｙｔｉｃｒｏｕｔｅｔｏｖａｌｅｒｉｃｂｉｏｆｕｅｌｓ［Ｊ］．ＡｐｐｌｉｅｄＥｎｅｒｇｙ，２０１３，１０２：１５７１６２．
［１６］ ＳｃｏｒｄｉａＤ，ＣｏｓｅｎｔｉｎｏＳＬ，ＬｅｅＪＷ，犲狋犪犾．Ｂｉｏｃｏｎｖｅｒｓｉｏｎｏｆｇｉａｎｔｒｅｅｄ（犃狉狌狀犱狅犱狅狀犪狓Ｌ．）ｈｅｍｉｃｅｌｕｌｏｓｅｈｙｄｒｏｌｙｓａｔｅｔｏｅｔｈａ
ｎｏｌｂｙ犛犮犺犲犳犳犲狉狊狅犿狔犮犲狊狊狋犻狆犻狋犻狊ＣＢＳ６０５４［Ｊ］．ＢｉｏｍａｓｓａｎｄＢｉｏｅｎｅｒｇｙ，２０１２，３９：２９６３０５．
［１７］ ＳｕｎＹＹ，ＹｕｅＱＹ，ＧａｏＢＹ，犲狋犪犾．Ｃｏｍｐａｒｉｓｏｎｏｆａｃｔｉｖａｔｅｄｃａｒｂｏｎｓｆｒｏｍ犃狉狌狀犱狅犱狅狀犪狓ＬｉｎｎｗｉｔｈＨ４Ｐ２Ｏ７ａｃｔｉｖａｔｉｏｎｂｙ
ｃｏｎｖｅｎｔｉｏｎａｌａｎｄｍｉｃｒｏｗａｖｅｈｅａｔｉｎｇｍｅｔｈｏｄｓ［Ｊ］．ＣｈｅｍｉｃａｌＥｎｇｉｎｅｅｒｉｎｇＪｏｕｒｎａｌ，２０１２，１９２：３０８３１４．
［１８］ ＲｉｂｅｃｈｉｎｉＥ，ＺａｎａｂｏｎｉＭ，ＧａｌｅｔｔｉＡＭＲ，犲狋犪犾．ＰｙＧＣ／ＭＳｃｈａｒａｃｔｅｒｉｚａｔｉｏｎｏｆａｗｉｌｄａｎｄａｓｅｌｅｃｔｅｄｃｌｏｎｅｏｆ犃狉狌狀犱狅犱狅狀犪狓，
ａｎｄｏｆｉｔｓｒｅｓｉｄｕｅｓａｆｔｅｒｃａｔａｌｙｔｉｃｈｙｄｒｏｔｈｅｒｍａｌｃｏｎｖｅｒｓｉｏｎｔｏｈｉｇｈａｄｄｅｄｖａｌｕｅｐｒｏｄｕｃｔｓ［Ｊ］．ＪｏｕｒｎａｌｏｆＡｎａｌｙｔｉｃａｌａｎｄＡｐｐｌｉｅｄ
Ｐｙｒｏｌｙｓｉｓ，２０１２，９４：２２３２２９．
［１９］ ＳｈａｔａｌｏｖＡＡ，ＰｅｒｅｉｒａＨ．Ｘｙｌｏｓｅｐｒｏｄｕｃｔｉｏｎｆｒｏｍｇｉａｎｔｒｅｅｄ（犃狉狌狀犱狅犱狅狀犪狓Ｌ．）：Ｍｏｄｅｌｉｎｇａｎｄｏｐｔｉｍｉｚａｔｉｏｎ０ｆｄｉｌｕｔｅａｃｉｄ
ｈｙｄｒｏｌｙｓｉｓ［Ｊ］．ＣａｒｂｏｈｙｄｒａｔｅＰｏｌｙｍｅｒｓ，２０１２，８７：２１０２１７．
［２０］ ＡｒｃａｄｅＡ，ＡｎｓｅｌｉｎＦ，ＲａｍＰａｎｔＰＦ，犲狋犪犾．ＡｐｐｌｉｃａｔｉｏｎｏｆＡＦＬＰ，ＲＡＰＤａｎｄＩＳＳＲｍａｒｋｅｒｓｔｏｇｅｎｅｔｉｃｍａｐｐｉｎｇｏｆＥｕｒｏｐｅａｎ
ａｎｄＪａｐａｎｅｓｅｌａｒｃｈ［Ｊ］．ＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，２０００，１００：２９９３０７．
［２１］ ＤｉｒｌｅｗａｎｇｅｒＥ，ＰｒｏｎｉｅｒＶ，ＰａｒｖｅｒｙＣ．ＧｅｎｅｔｉｃｌｉｎｋａｇｅｍａｐｏｆＰｅａｃｈ（犘狉狌狀狌狊狆犲狉狊犻狕犲犪Ｌ．）ｕｓｉｎｇｍｏｒｐｈｏｌｏｇｉｃａｌａｎｄｍｏｌｅｃｕ
ｌａｒｍａｒｋｅｒｓ［Ｊ］．ＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，１９９８，９７：８８８８８９．
［２２］ ＦｅｒｎａｎｄｅｚＭＥ，ＦｉｇｕｅｉｒａｓＡＭ，ＢｅｎｉｔｏＣ．ＴｈｅｕｓｅｏｆＩＳＳＲａｎｄＲＡＰＤｍａｒｋｅｒｓｆｏｒｄｅｔｅｃｔｉｎｇＤＮＡｐｏｌｙｍｏｒｐｈｉｓｍ，ｇｅｎｏｔｙｐｅ
ｉｄｅｎｔｉｆｉｃａｔｉｏｎａｎｄｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙａｍｏｎｇｂａｒｌｅｙｃｕｌｔｉｖａｒｓｗｉｔｈｋｎｏｗｎｏｒｉｇｉｎ［Ｊ］．ＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，２００２，
１０４：８４５８５１．
［２３］ 李飞飞，崔大方，羊海军，等．中国新疆紫花苜蓿复合体３个种的遗传多样性及亲缘关系研究［Ｊ］．草业学报，２０１２，２１（１）：
１９０１９８．
［２４］ 刘欢，慕平，赵桂琴．燕麦种质资源遗传多样性ＩＳＳＲ研究［Ｊ］．草业学报，２０１２，２１（４）：１１６１２４．
［２５］ 王关林，方宏筠．植物基因工程（第二版）［Ｍ］．北京：科学出版社，２００２：７４２７４４．
［２６］ ＢａｌｏｇｈＥ，ＪｏｈｎＭ ＨＪ，ＣｚａｋóＭ，犲狋犪犾．Ｄｅｆｅｃｔｉｖｅｄｅｖｅｌｏｐｍｅｎｔｏｆｍａｌｅａｎｄｆｅｍａｌｅｇａｍｅｔｏｐｈｙｔｅｓｉｎ犃狉狌狀犱狅犱狅狀犪狓 Ｌ．
（Ｐｏａｃｅａｅ）［Ｊ］．ＢｉｏｍａｓｓａｎｄＢｉｏｅｎｅｒｇｙ，２０１２，４５：２６５２６９．
［２７］ ＨａｄｄａｄｃｈｉＡ，ＧｒｏｓｓＣＬ，ＦａｔｅｍｉＭ．Ｔｈｅｅｘｐａｎｓｉｏｎｏｆｓｔｅｒｉｌｅ犃狉狌狀犱狅犱狅狀犪狓 （Ｐｏａｃｅａｅ）ｉｎｓｏｕｔｈｅａｓｔｅｒｎＡｕｓｔｒａｌｉａｉｓａｃｃｏｍ
２７２ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．３
ｐａｎｉｅｄｂｙｇｅｎｏｔｙｐｉｃｖａｒｉａｔｉｏｎ［Ｊ］．ＡｑｕａｔｉｃＢｏｔａｎｙ，２０１３，１０４：１５３１６１．
［２８］ ＡｈｍａｄＲ，ＬｉｏｗＰＳ，ＳｐｅｎｃｅｒＤＦ，犲狋犪犾．Ｍｏｌｅｃｕｌａｒｅｖｉｄｅｎｃｅｆｏｒａｓｉｎｇｌｅｇｅｎｅｔｉｃｃｌｏｎｅｏｆｉｎｖａｓｉｖｅ犃狉狌狀犱狅犱狅狀犪狓ｉｎｔｈｅＵｎｉｔ
ｅｄＳｔａｔｅｓ［Ｊ］．ＡｑｕａｔｉｃＢｏｔａｎｙ，２００８，８８（２）：１１３１２０．
［２９］ 杨允菲，刘庚长，张宝田．羊草种群年龄结构及无性繁殖对策的分析［Ｊ］．植物学报，１９９５，３７（２）：１４７１５３．
［３０］ ＢａｒｓｏｕｍＮ．Ｒｅｌａｔｉｖｅｃｏｎｔｒｉｂｕｔｉｏｎｓｏｆｓｅｘｕａｌａｎｄａｓｅｘｕａｌｒｅｇｅｎｅｒａｔｉｏｎｓｔｒａｔｅｇｉｅｓｉｎ犘狅狆狌犾狌狊狀犻犵狉犪ａｎｄ犛犪犾犻狓犪犾犫犪ｄｕｒｉｎｇｔｈｅ
ｆｉｒｓｔｙｅａｒｓｏｆｅｓｔａｂｌｉｓｈｍｅｎｔｏｎａｂｒａｉｄｅｄｇｒａｖｅｌｂｅｄｒｉｖｅｒ［Ｊ］．ＥｖｏｌｕｔｉｏｎａｒｙＥｃｏｌｏｇｙ，２００２，１５：２５５２７９．
［３１］ 夏立群，李建强，李伟．论克隆植物的遗传多样性［Ｊ］．植物学通报，２００２，１９（４）：４２５４３１．
［３２］ 汤俊兵，肖燕，安树青．根茎克隆植物生态学研究进展［Ｊ］．生态学报，２０１０，３０（１１）：３０２８３０３６．
［３３］ 张玉芬，张大勇．克隆植物的无性与有性繁殖对策［Ｊ］．植物生态学报，２００６，３０（１）：１７４１８３．
［３４］ 赵丽萍，许卉．芦竹在滨海盐碱地的开发应用及栽培技术［Ｊ］．北方园艺，２００７，（７）：１６４１６５．
［３５］ 韩志萍，胡晓斌，胡正海．芦竹修复镉汞污染湿地的研究［Ｊ］．应用生态学报，２００５，１６（５）：９４５９５０．
［３６］ 文仕知，李铁华，陈亮明．湘中紫色页岩区芦竹生物量及养分分布规律研究［Ｊ］．中南林业科技大学学报，２００８，２８（５）：
１５．
［３７］ 詹世雄，曾宪威．植物核心种质构建方法研究进展［Ｊ］．中国农学通报，２０１０，２６（３）：２７９２８２．
犌犲狀犲狋犻犮犱犻狏犲狉狊犻狋狔狅犳犲狀犲狉犵狔犵狉犪狊狊犃狉狌狀犱狅犱狅狀犪狓狉犲狏犲犪犾犲犱犫狔犐犛犛犚犿犪狉犽犲狉狊
ＺＥＮＧＨａｎｙｕａｎ１，２，ＷＥＩＬｉｎ２，ＬＩＵＰｅｎｇ２，ＬＩＵＸｕａｎｍｉｎｇ３
（１．ＣｏｌｅｇｅｏｆＢｉｏｓｃｉｅｎｃｅａｎｄＢｉｏｔｅｃｈｎｏｌｏｇｙ，ＨｕｎａｎＡｇｒｉｃｕｌｔｕｒａｌＵｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ４１０１２８，Ｃｈｉｎａ；
２．ＤｅｐａｒｔｍｅｎｔｏｆＬｉｆｅＳｃｉｅｎｃｅ，ＨｕａｉｈｕａＣｏｌｅｇｅ，Ｈｕａｉｈｕａ４１８００８，Ｃｈｉｎａ；３．Ｉｎｓｔｉｔｕｔｅ
ｏｆＬｉｆｅＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，ＨｕｎａｎＵｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ４１００８２，Ｃｈｉｎａ）
犃犫狊狋狉犪犮狋：Ｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙａｎｄｓｔｒｕｃｔｕｒｅｏｆ１１犃狉狌狀犱狅犱狅狀犪狓ｐｏｐｕｌａｔｉｏｎｓ（ａｂｏｕｔ８９ｉｎｄｉｖｉｄｕａｌｓｉｎｔｏｔａｌ）ｆｒｏｍ６
ｐｒｏｖｉｎｃｅｓｗｅｒｅａｎａｌｙｓｅｄｗｉｔｈＩＳＳＲ（ｉｎｔｅｒｓｉｍｐｌｅｓｅｑｕｅｎｃｅｒｅｐｅａｔ）ｍａｒｋｅｒｓ．ＦｒｏｍＵＢＣ８０１－９００，ｓｅｖｅｎｅｆｆｉ
ｃｉｅｎｔＩＳＳＲｐｒｉｍｅｒｓｗｅｒｅｓｃｒｅｅｎｅｄａｎｄａｔｏｔａｌｏｆ４０ｌｏｃｉｗｅｒｅｄｅｔｅｃｔｅｄｏｆｗｈｉｃｈ３８ｗｅｒｅｐｏｌｙｍｏｒｐｈｉｃ．Ｔｈｅｌｅｎｇｔｈ
ｏｆＰＣＲｐｒｏｄｕｃｔｓｒａｎｇｅｄｆｒｏｍ１００ｔｏ２０００ｂｐ．Ａｔｔｈｅｓｐｅｃｉｅｓｌｅｖｅｌ，ｔｈｅＳｈａｎｎｏｎｉｎｄｅｘ（Ｉ）ｗａｓ０．４９４，Ｎｅｉ’ｓ
ｇｅｎｅｄｉｖｅｒｓｉｔｙｉｎｄｅｘ（Ｈ）ｗａｓ０．３２９，ｔｈｅｃｏｅｆｆｉｃｉｅｎｔｏｆｇｅｎｅｔｉｃｄｉｆｆｅｒｅｎｔｉａｔｉｏｎ（Ｇｓｔ）ｗａｓ０．８２６，ａｎｄｔｈｅｇｅｎｅ
ｆｌｏｗ（Ｎｍ）ｗａｓ０．０９７．ｂｕｔｔｈｅａｖｅｒａｇｅｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙｗｉｔｈｉｎｐｏｐｕｌａｔｉｏｎｓ（Ｈｓ）ｗａｓｏｎｌｙ０．０４０．Ｔｈｉｓｓｈｏｗｅｄ
ｔｈａｔｔｈｅｇｅｎｅｄｉｖｅｒｓｉｔｙｏｆ犃．犱狅狀犪狓ｗａｓｈｉｇｈａｔｔｈｅｓｐｅｃｉｅｓｌｅｖｅｌｂｕｔｖｅｒｙｌｏｗｗｉｔｈｉｎｐｏｐｕｌａｔｉｏｎｓ．ＵＰＧＭＡ
ｃｌｕｓｔｅｒｉｎｇｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｅ１１ｐｏｐｕｌａｔｉｏｎｓｏｆ犃．犱狅狀犪狓ｃｏｕｌｄｂｅｄｉｖｉｄｅｄｉｎｔｏｔｈｒｅｅｇｒｏｕｐｓ：ｔｈｅｐｏｐｕｌａ
ｔｉｏｎｓｆｒｏｍＺｈｅｊｉａｎｇａｎｄＪｉａｎｇｓｕｃｌｕｓｔｅｒｅｄｉｎｔｏｏｎｅｇｒｏｕｐ，ｔｈｅｐｏｐｕｌａｔｉｏｎｓｆｒｏｍＨｕｎａｎ，ＧｕｉｚｈｏｕａｎｄＧｕａｎｇｘｉ
ｃｌｕｓｔｅｒｅｄｉｎｔｏａｓｅｃｏｎｄｇｒｏｕｐ，ａｎｄｔｈｅＹｕｎｎａｎｐｏｐｕｌａｔｉｏｎｉｎｔｏａｓｅｐａｒａｔｅ，ｔｈｉｒｄ，ｄｉｓｔｉｎｃｔｇｒｏｕｐ．Ｔｈｅｒｅｗａｓ
ｔｈｅｒｅｆｏｒｅｓｏｍｅｃｏｒｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｇｅｎｅｔｉｃｄｉｓｔａｎｃｅａｎｄｇｅｏｇｒａｐｈｉｃｄｉｓｔａｎｃｅ，ｂｕｔｉｔｗａｓｎｏｔａｂｓｏｌｕｔｅ．
犓犲狔狑狅狉犱狊：ＩＳＳＲ；ｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙ；犃狉狌狀犱狅犱狅狀犪狓；ｐｏｐｕｌａｔｉｏｎ
３７２第２２卷第３期 草业学报２０１３年