全 文 :书我国9种披碱草属植物的系统学关系
严学兵1,周禾2,王2,王成章1,郭玉霞1
(1.河南农业大学牧医工程学院,河南 郑州450002;2.中国农业大学草地研究所,北京100094)
摘要:采用单一遗传标记研究不同分类群之间的亲缘关系,由此可产生不同系统发育关系,对正确解释不同类群之
间的亲缘关系带来很多困难。本试验采用形态学、等位酶和微卫星3种遗传标记对我国9种披碱草属植物的亲缘
关系进行分析,结果表明,形态学标记把9种披碱草分为披碱草组(Sect.犜狌狉犮狕犪狀犻狀狅狏犻犪(Nevski)Tzvel.)和老芒麦
组(Sect.犈犾狔犿狌狊)两大类;等位酶标记虽然没有严格按照穗子的形状分为2组(披碱草组和老芒麦组),但同一组内
或相同的基因组构成具有较近的亲缘关系;微卫星标记没有把9种披碱草分成几个明显的组,且各种间表现出较
远的亲缘关系。3种遗传标记之间均表现出显著的相关关系,不同标记反映的亲缘关系与种本身之间,形态学标记
的遗传关系支持基于形态特征分类系统而与种本身的亲缘关系更加相符(犚=0.3276,犘<0.01),其次是等位酶标
记(犚=0.3006,犘<0.01),微卫星偏离相对较远(犚=0.1350,犘<0.05)。因此认为形态学是最直接、最能反映披
碱草属植物本身特点的标记方法,但要考虑指标之间的线性化,并非越多越好;等位酶分析是一种研究种群内部遗
传结构和变异比较好的方法,也是种、种群间遗传关系有价值的预测者;尽管SSR在属或种间分析不是最适用的,
但由于SSR是最多态的标记,对种内分析非常适合。
关键词:披碱草属;形态学;等位酶;微卫星;亲缘关系
中图分类号:S543+.901;Q944 文献标识码:A 文章编号:10045759(2009)03007412
披碱草属(犈犾狔犿狌狊)是禾本科(Poaceae)小麦族(Triticeae)重要的多年生属,主要分布在欧亚大陆,南、北美洲
北部也有少量分布。其垂直分布从海拔为几米的海滩一直到海拔5200m以上的喜马拉雅山区[1]。披碱草属植
物的多数物种为草原和草甸的重要组成成分,许多种类是优良的牧草,饲用价值极高。同时,有些种类还具有麦
类作物所缺乏的高产、优质、抗病、抗虫和抗逆等基因,这些优良基因可以通过远缘杂交、染色体工程技术、基因工
程技术等现代遗传和生物技术的方法从野生种类中转移到栽培小麦(犜狉犻狋犻犮狌犿犪犲狊狋犻狏狌犿)和大麦(犎狅狉犱犲狌犿狏狌犾
犵犪狉犲)的遗传背景中去。因此,作为麦类作物的野生近缘种遗传资源,披碱草属植物具有重要的经济价值。
关于披碱草属的分类系统,目前争议是将三者合并为1个属,即披碱草属(广义),还是将其划分为3个独立
的属。Bentham[2]和Nevski[3]等的北美和欧洲的一些学者,特别是Dewey[4]和Lve[5]以染色体组为分类依据,
认为形态特征在披碱草属的界定上不太重要,他们不考虑每一穗节的小穗数量而支持广义的披碱草属分类方法,
把鹅观草属、猬草属归于披碱草属之中。有不少学者,特别是我国的禾草学家则主张将鹅观草属、猬草属从披碱
草属中独立出来[6~8],这种分类系统主要依据穗部的形态特征。按照广义的披碱草属概念,全世界包括150多种
植物,我国80多种,是目前国际上最为通用的概念。而根据狭义的披碱草属概念,我国仅有12个物种[7]。
关于此分类系统的争议,国内学者也曾进行了研究和论述[7,9~13],但更多证据来自于形态和细胞学方面的研
究。随着分子生物技术的发展,人们期待从分子水平寻找第三方面的证据。周永红等[13]利用随机引物扩增多态
性技术(randomlyamplifiedpolymorphicDNA,RAPD)分析了10种披碱草属植物的系统学关系,结果表明披碱
草属中四倍体和六倍体物种各自聚为一支,根据老芒麦(犈.狊犻犫犻狉犻犮狌狊)与犬草(犈.犮犪狀犻狀狌狊)的亲缘关系较近,支持
细胞学分类的学者将犬草放入披碱草属[4,5]。遗传标记常用于研究其他不同分类群之间的亲缘关系,由此产生
许多不同的系统发育关系,但是由于大多采用单一分子标记,所以对于正确解释不同类群之间的亲缘关系带来很
多困难。本研究选用9种国内披碱草属植物 (含StH、StYH染色体组),利用3种分子标记分析,旨在检验3种
分子标记在反映披碱草属植物亲缘关系之间的差异性,为利用多方面的证据进行分类和种间关系探讨时提供新
74-85
2009年6月
草 业 学 报
ACTAPRATACULTURAESINICA
第18卷 第3期
Vol.18,No.3
收稿日期:20080806;改回日期:20081006
基金项目:国家自然科学基金(30871531)和青海省重大科技攻关项目(21031423)资助。
作者简介:严学兵(1974),男,河南周口人,副教授,博士(后)。Email:yxbbjzz@163.com
通讯作者。
思路。
1 材料与方法
1.1 试验材料
试验材料共有9种披碱草属植物(狭义的披碱草属概念):垂穗披碱草(犈.狀狌狋犪狀狊)、老芒麦、短芒披碱草(犈.
犫狉犲狏犻犪狉犻狊狋犪狋狌狊)、圆柱披碱草(犈.犮狔犾犻狀犱狉犻犮狌狊)、黑紫披碱草(犈.犪狋狉犪狋狌狊)、紫芒披碱草(犈.狆狌狉狆狌狉犪狉犻狊狋犪狋狌狊)、披碱
草(犈.犱犪犺狌狉犻犮狌狊)、青海披碱草(犈.犵犲犿犻狀犪狋狌狊)和无芒披碱草(犈.狊狌犫犿狌狋犻犮狌狊)的40个野生种群(表1)。种子采
收后在中国农业大学校园进行盆栽试验,每个种群播种3盆,每盆定苗10株。盆的大小为25cm(上口直径)×
17cm(下口直径)×29cm(高)。土壤的理化性质为pH值7.13,有效磷(OsenP)0.07mg/kg,有效钾0.28mg/
kg,全氮0.74%,全磷2.94%,有机质4.05%。
1.2 形态特征观测项目及测定方法
选取成熟期植株的表型性状(在同一种群内调查):叶片数、穗长、穗宽、穗节数、小穗长、小穗宽、生殖枝长、旗
叶与穗基部长度、外稃、内稃、芒长、颖长、颖宽、小穗柄、单株重、单穗重、旗叶宽、旗叶长、倒2叶片宽、倒2叶长、
株高均值、分蘖数是在成熟期随机测量30个穗子,求其平均值。生长速度测定分别在幼苗期(2003年12月)和
幼苗-分蘖期(2004年2月)测定株高进行计算而得。穗长、穗宽、穗节数、小穗长、小穗宽、生殖枝长、旗叶与穗
基部长度、外稃、内稃、芒长、颖长、颖宽、小穗柄等指标,均测10个穗子,测量部位为穗的中部,求其平均值。
1.3 等位酶分析
1.3.1 酶系统和酶样品制备 共对9种酶进行检测,5种结果稳定、可靠的酶系统用于遗传分析,分别为乙醇脱
氢酶 (ADH,EC1.1.1.1),天冬氨酸转氨酶(AAT,EC2.6.1.1),酯酶 (EST,EC3.1.1.2),过氧化物酶 (PRX,
EC1.II.1.7)和苹果酸脱氢酶 (MDH,EC1.1.1.37)。
每个种群取样5株,取二叶一心期幼苗单株,分3次加入150μL的TrisHCl提取缓冲液。冰镇研磨,用吸
管吸入1.5mL离心管中,保存于-70℃的低温冰箱中备用。点样前需在12000r/min的离心机中离心10min,
取30μL的上清液进行电泳。
1.3.2 聚丙烯酰胺凝胶电泳(SDS-PAGE,sodiumdodecylsulfatepolyacrylamidegelelectrophoresis) 采用
SDS-PAGE技术对酶系统进行检测。凝胶板规格为160mm×160mm×1.0mm,每板25个宽5mm的样品
槽。丙烯酰胺浓度7.0%的分离胶,丙烯酰胺浓度5.0%的浓缩胶。分离胶和浓缩胶缓冲液为TrisHCl(pH值
分别为8.9和6.8),电极缓冲液为Tris甘氨酸(pH值8.3),先在120V电压下电泳20~30min,将电压升为
200V继续电泳4~5h,直到溴酚蓝前沿移至距分离胶起点7.5cm处,停止电泳,剥胶染色。显色后凝胶片在
7%的冰醋酸中固定24h。酶的组织化学染色法按王中仁[14]的方法进行。
表1 披碱草属供试植物来源,海拔,经纬度和生境条件
犜犪犫犾犲1 犗狉犻犵犻狀,犪犾狋犻狋狌犱犲,犾狅狀犵犻狋狌犱犲,犾犪狋犻狋狌犱犲犪狀犱犺犪犫犻狋犪狋狅犳犈犾狔犿狌狊狊狆犲犮犻犲狊狌狊犲犱犻狀狋犺犻狊狊狋狌犱狔
种名
Species
染色体数
Chromosome
number
染色体组
Chromosome
来源
Origin
编号
Identity
海拔
Altitude
(m)
纬度
Latitude
(N)
经度
Longitude
(E)
生境
Habitat
垂穗披碱草
犈.狀狌狋犪狀狊
42 StYH
甘肃天祝Tianzhu,Gansu GJC 2915 37°11′38″102°48′14″ 高寒草甸Alpinemeadow
甘肃合作 Hezuo,Gansu GHC 3040 35°02′18″102°55′31″ 封育 的 天 然 草 地 Enclosed
rangeland
青海玛沁 Maqin,Qinghai QMC 3289 35°16′32″100°39′22″ 高寒草甸Alpinemeasow
青海同德Tongde,Qinghai QHHC 3487 34°42′36″100°42′37″ 高寒草甸Alpinemeadow
青海玛沁 Maqin,Qinghai QMDC1 3838 34°29′28″100°09′32″ 高寒草甸Alpinemeadow
青海玛沁 Maqin,Qinghai QMD(W1)3908 34°29′42″100°09′47″ 灌丛草地Shrub
青海玛沁 Maqin,Qinghai QMDC2 3916 34°40′27″100°24′52″ 高寒草甸Alpinemeadow
青海玛沁 Maqin,Qinghai QMXC 3724 35°02′05″099°12′38″ 河漫滩Swampalongtheriver
青海玛沁 Maqin,Qinghai QMHC 3867 34°32′06″100°25′16″ 高寒草甸Alpinemeadow
57第18卷第3期 草业学报2009年
续表1 Continued
种名
Species
染色体数
Chromosome
number
染色体组
Chromosome
来源
Origin
编号
Identity
海拔
Altitude
(m)
纬度
Latitude
(N)
经度
Longitude
(E)
生境
Habitat
短芒披碱草
犈.犫狉犲狏犻犪
狉犻狊狋犪狋狌狊
42 StYH
甘肃天祝Tianzhu,Gansu GJD 2883 37°12′15″102°47′42″ 路边Roadside
甘肃合作 Hezuo,Gansu GHD 3040 35°02′18″102°55′31″ 封育 的 天 然 草 地 Enclosed
rangeland
青海玛沁 Maqin,Qinghai QMDP 3280 35°16′08″100°39′29″ 高寒草甸Alpinemeadow
青海玛沁 Maqin,Qinghai QMHD 3628 34°32′06″100°25′14″ 山坡草地Sloperangeland
青海兴海Xinghai,Qinghai QXED 4186 35°27′41″ 99°29′13″ 水蚀沟内Ditcherodedbywater
青海玛沁 Maqin,Qinghai QMDL 3282 35°16′08″100°39′29″ 疏 林 草 地 Rangeland with
sparseforest
黑紫披碱草
犈.犪狋狉犪狋狌狊
42 StYH 甘肃天祝Tianzhu,Gansu GJH 2878 37°11′53″102°47′13″ 河漫滩草地Swampmeadow
老芒麦
犈.狊犻犫犻狉犻犮狌狊
28 StH
河北沽源Guyuan,Hebei HYL 1458 41°45′29″116°09′32″ 人工草地Pasture
河北沽源Guyuan,Hebei HGL(N)1350 41°53′26115°52′24″ 封育 的 典 型 草 原 Enclosed
steppe
甘肃天祝Tianzhu,Gansu GJL(Y) 2881 37°11′05″102°47′38″ 高寒草甸Alpinemeadow
甘肃天祝Tianzhu,Gansu GJLZ 2887 37°12′18″102°47′11″ 人工草地Pasture
甘肃合作 Hezuo,Gansu GHL(1) 2532 35°13′25″102°47′03″ 封育 的 天 然 草 地 Enclosed
rangeland
甘肃合作 Hezuo,Gansu GHL(2) 2693 35°28′03″102°55′31″ 封育 的 人 工 草 地 Enclosed
rangeland
青海玛沁 Maqin,Qinghai QML(Y)3282 35°16′08″100°39′29″ 疏 林 草 地 Rangeland with
sparseforest
青海同德Tongde,Qinghai QM3 3280 35°16′22″100°39′28″ 试验地Trialland
青海同德Tongde,Qinghai QHYF 3786 34°54′29″100°53′25″ 高寒草甸Alpinemeadow
披碱草
犈.犱犪犺狌狉犻犮狌狊
42 StYH
甘肃天祝Tianzhu,Gansu GJM 2874 37°10′18″102°47′31″ 排水沟边Ditch
甘肃永靖Yongjing,Gansu GYP 1624 35°09′31″103°04′02″ 路边Roadside
青海同德Tongde,Qinghai QHHM 3470 34°43′19″100°42′49″ 路边Roadside
青海披碱草
犈.犵犲犿犻狀犪狋狌狊
/ /
青海玛沁 Maqin,Qinghai QMX(W)4010 35°02′05″ 99°12′38″ 灌丛草地Shrubs
青海兴海Xinghai,Qinghai QXE(W)4034 35°27′44″ 99°29′12″ 河漫滩Swampalongtheriver
无芒披碱草
犈.狊狌犫犿狌狋犻犮狌狊
42 StYH 甘肃天祝Tianzhu,Gansu GJW 2883 37°12′15″102°47′42″ 路边Roadside
圆柱披碱草
犈.犮狔犾犻狀犱狉犻犮狌狊
42 StYH
河北沽源Guyuan,Hebei HGY(N)1353 41°52′45″115°51′57″ 封育 的 典 型 草 原 Enclosed
steppe
河北沽源Guyuan,Hebei HGY(K)1363 41°52′31″115°48′20″ 拓荒地 Wilderness
甘肃天祝Tianzhu,Gansu GJY 2918 37°11′33″102°48′18″ 高寒草甸Alpinemeadow
紫芒披碱草
犈.狆狌狉狆狌
狉犪狉犻狊狋犪狋狌狊
42 StYH
河北沽源Guyuan,Hebei HGZ(N)1357 41°52′01″115°52′27″ 路边Roadside
河北沽源Guyuan,Hebei HYZK 1443 41°45′28″116°10′18″ 人工草地Pasture
甘肃合作 Hezuo,Gansu GHY 2107 36°07′52″102°18′02″ 干河滩Dryriverway
甘肃临夏Linxia,Gansu GL(W) 1973 35°06′24″103°01′16″ 沟旁Ditch
青海玛沁 Maqin,Qinghai QMZ 3282 35°16′08″100°39′29″ 疏 林 草 地 Rangeland with
sparseforest
河北沽源Guyuan,Hebei HYZH 1447 41°45′39″116°09′04″ 路边Roadside
67 ACTAPRATACULTURAESINICA(2009) Vol.18,No.3
1.4 微卫星(simplesequencerepeats,SSR)分析
1.4.1 植物总DNA的提取方法 取每个种群30个单株幼苗期叶片,等量均匀混合,迅速保存于-80℃冰箱长
期备用。根据Doyle[15]的方法稍加改动提取种群DNA混合样。
1.4.2 SSR扩增 a)SSR引物:试验选择2类引物,一类是来自于披碱草属植物基因组的14条引物,其中7条
来自于犬草 (缩写为ECGA)[16],另外7条来自于阿拉斯加披碱草(犈.犪犾犪狊犽犪狀狌狊)基因组(缩写为EAGA)[17];另
一类是来自于小麦基因组的8条引物(缩写为 WMS)[18,19],具体见表2。
b)聚合酶链式反应(PCR,polymerasechainreaction)体系和程序:所有PCR反应在型号为PTC200的
PCR仪上进行,来自于犬草基因组引物的反应体系和PCR程序采用Sun等[16]的方法,来自于阿拉斯加披碱草基
因组的引物采用Sun等[17]的方法;另一类来自于小麦基因组的8条引物(缩写为 WMS)采用Rder等[18]和Pl
aschke等[19]的方法。PCR反应物在冷却后加1/5体积的上样缓冲液在95℃变性5min,然后迅速在-80℃冷却
备用。
表2 试验所选用犛犛犚引物
犜犪犫犾犲2 犈犾狔犿狌狊犪狀犱犜狉犻狋犻犮狌犿犿犻犮狉狅狊犪狋犲犾犻狋犲狆狉犻犿犲狉狊狌狊犲犱犻狀犛犛犚狊狋狌犱狔
引物Primers 序列Sequence(5′→3′) 引物Primers 序列Sequence(5′→3′)
ECGA22 GAAGGTGACTAGGTCCAAC
ATAGTCTCGGTCAGGCTC
EAGA102 CACAACCGTGGCAGTGGT
TCAACATTCAAATTTGGCAGAC
ECGA114 CTTATATCTTGTGGGTTATCAT
GATCTGATACGTGACATCTCA
EAGA103 GACTGCTAAACTAACAACAAAACT
ATCTCACACTGGCTCATGTC
ECGA4 TGATCAAGAAGAGGAACAT
GATAAGATCGTGACTCTCCT
EAGA104 AGCCAGACACAAGTGACTATG
GACTACGGTCCATAGATAGTATCT
ECGA89 TTAGCTCTTTACTTATTCAAAC
TCCTATGATCAAGCACAAG
WMS37 ACTTCATTGTTGATCTTGCAT
CGACGAATTCCCAGCTAAAC
ECGA210 CGACAACTAGTGGATCAAA
GAAGTACTCTCGAGAAGCTT
WMS2 CTGCAAGCCTGTGATCAACT
CATTCTCAAATGATCGAACA
ECGT36 TCACAGAGTGATTACATGCG
ACATGTAACTGGAGTCGAGC
WMS43 CACCGACGGTTTCCCTAGAGT
GGTGAGTGCAAATGTCATGTG
ECGA125 TGCTTCCAACTTGCTCA
TGCATCTGTGTGTCCACA
WMS44 GTTGAGCTTTTCAGTTCGGC
ACTGGCATCCACTGAGCTG
EAGA13 ACACACTGAGACTGCTTTGTCA
TGCGAAGGTTGAGTTATGCT
WMS46 GCACGTGAATGGATTGGAC
TGACCCAATAGTGGTGGTCA
EAGA51 TCCACTGTGCTACTTCTAGCT
AACTGAATGAATGTGTGCAGT
WMS47 TTGCTACCATGCATGACCAT
TTCACCTCGATTGAGGTCCT
EAGA101 GATGATGTCAATATGCTGCAAT
AGGCCAGACTTATGAACACTAT
WMS106 CTGTTCTTGCGTGGCATTAA
AATAAGGACACAATTGGGATGG
EAGA74 AGGTTGTATCTGCTGATACGATC
TGCAAGTCGAGCTGTCACTAG
WMS159 GGGCCAACACTGGAACAC
GCAGAAGCTTGTTGGTAGGC
1.4.3 变性聚丙烯酰胺凝胶电泳和染色 制备6%变性聚丙烯酰胺凝胶→电泳(电压300V,1×TAE电泳4~
5h)→0.2%硝酸银染色→凝胶成像系统扫描成像。
1.5 数据计算和统计处理
1.5.1 形态数据的分析与处理 以上述各形态学指标(数据太多且篇幅有限没有显示)作为原始数据,对每种披
碱草属植物形态指标计算其平均值。不同种间的差异用遗传相似性(geneticsimilarity,GS)[20]表示,不同种的
77第18卷第3期 草业学报2009年
披碱草材料间相似性的距离矩阵采用不加权成对算术平均法 (unweightedpairgroupmethodwitharithmetic
average,UPGMA)[21]进行聚类分析,建立亲缘关系树图。分别采用Excel和NTSYSPC(版本2.1)软件进行计
算和聚类分析。
1.5.2 等位酶数据的分析与处理 酶谱的遗传分析参考王中仁[14]的工作,结合谱带在天然小群体中的分离式
样和酶分子结构推断,并通过电泳分析获得不同个体基因型频率进而计算遗传参数。采用Popgen1.31软件计
算遗传相似性系数。不同种间采用遗传距离和UPGMA法进行聚类分析,建立亲缘关系树状图。每个披碱草属
植物的所有种群构成了本物种的遗传数据,进行遗传距离和不同数据间的相关性分析,计算运用NTSYSPC(版
本2.1)程序。
1.5.3 SSR数据的分析与处理 每个样品的电泳条带按有或无记录,电泳条带存在时赋值1,无则赋值0,缺失
用“.”表示。采用Popgen1.31软件对不同种群的披碱草材料间,以简单匹配系数为基础(simplematchcoeffi
cient),计算遗传相似性系数,用UPGMA法对披碱草属种进行聚类分析,建立亲缘关系树图,聚类分析运用NT
SYSPC(版本2.1)程序。
1.5.4 3种遗传数据间的 Mantel检验 3种标记数据之间的比较用NTSYSPC(版本2.1)程序,将基于3种遗
传标记的不同材料间的数据矩阵,选择Graphics的 Matrixcomparisonplot子程序,采用 Mantel检验[22]标准Z
和相关系数(犚)进行相关性分析。原材料本身的遗传距离矩阵:同种材料遗传距离为0,不同种材料的为1而建
立。
2 结果与分析
2.1 形态学标记的聚类和亲缘关系分析
遗传距离或相似性系数常用来判断群体之间的亲缘关系,当遗传距离为1(相似性系数为0)时,表明2个群
体完全不一样,无亲缘关系;当遗传距离为0(相似性系数为1)时,表明2个群体完全一样。把一个种的所有种群
作为一个样本,种作为一个组进行种间关系比较,对9种披碱草28个形态指标的平均值进行聚类分析,不同披碱
草种间相似性系数变化范围为0.8601~0.9821(图1)。我国9种披碱草属植物可以明显分为两大类:披碱草组
(Sect.犜狌狉犮狕犪狀犻狀狅狏犻犪(Nevski)Tzvel.)的3种披碱草为一类,老芒麦组(Sect.犈犾狔犿狌狊)的6种为第2类。老芒麦
组的短芒披碱草、无芒披碱草亲缘关系最近,首先聚在一起,后垂穗披碱草、老芒麦、青海披碱草,后来是黑紫披碱
草。披碱草组的3种披碱草中,披碱草和圆柱披碱草亲缘关系较近,然后和紫芒披碱草聚为一类。
2.2 等位酶标记的聚类和亲缘关系分析
在试验中,除了老芒麦是四倍体之外,其他均为六倍体植物。通过对披碱草属40个野生种群的200份样品
的5种酶进行遗传分析,9种披碱草属植物的亲缘关系,遗传相似性系数为0.6649~0.9627,垂穗披碱草和青海
披碱草最为相近,圆柱披碱草和黑紫披碱草亲缘关系最远。根据聚类的总体来看(图2),虽然没有严格按照穗子
的形状(披碱草组和老芒麦组)分为2组,但披碱草组的紫芒披碱草、圆柱披碱草、披碱草亲缘关系较近。在老芒
麦组内,垂穗披碱草和青海披碱草,短芒披碱和黑紫披碱草亲缘关系较近。另外,老芒麦组的垂穗披碱草和披碱
草组具有相同染色体组(StYH)的披碱草、圆柱披碱草和紫芒披碱草亲缘关系也较近(图3)。
2.3 SSR标记的聚类和亲缘关系分析
基于对披碱草属40个野生种群的微卫星标记的遗传分析(以引物 WMS43为例,图4),我国不同披碱草种
间相似性系数变化范围为0.4801~0.9071。根据遗传距离,用UPGMA法对9种披碱草进行聚类分析并绘制
树状图(图5)。从种间亲缘关系系统聚类图(图5)来看,披碱草属的9个种并没有分成几个明显的组,且各种间
表现出较远的亲缘关系。其中,老芒麦和青海披碱草亲缘关系最近(0.9071),其次是紫芒披碱草和短芒披碱草,
老芒麦和黑紫披碱草亲缘关系最远(0.4801)。
2.4 不同标记之间的比较
对3种标记方法(形态、等位酶和SSR)之间以及与种本身的相关关系见图6~8。对形态、等位酶和SSR三
组数据矩阵进行 Mantel检测,得出相关系数和显著性检验结果见表3。形态和等位酶、形态和SSR、等位酶和
SSR之间的变化均表现有一定相似性,经过进一步 Mantel检验发现(表3),形态和等位酶、形态和SSR、等位酶
87 ACTAPRATACULTURAESINICA(2009) Vol.18,No.3
和SSR之间的相关系数分别为0.3613,0.2993和0.1489,3种方法之间有显著的相关关系,其中形态和等位
图1 形态学标记披碱草属种间亲缘关系
犉犻犵.1 犇犲狀犱狉狅犵狉犪犿狅犳犿狅狉狆犺狅犾狅犵犻犮犪犾狋狉犪犻狋狊犫犲狋狑犲犲狀狊狆犲犮犻犲狊
图2 等位酶标记披碱草属种间亲缘关系
犉犻犵.2 犚犲犾犪狋犻狅狀狊犺犻狆犪犿狅狀犵犈犾狔犿狌狊狊狆犲犮犻犲狊犫狔犪犾狅狕狔犿犲犱犪狋犪
表3 披碱草属种群3种类型遗传标记数据关系
犜犪犫犾犲3 犚犲犾犪狋犻狅狀狊犺犻狆狊狅犳狋犺狉犲犲狋狔狆犲狊狅犳犵犲狀犲狋犻犮犱犪狋犪犪犿狅狀犵狆狅狆狌犾犪狋犻狅狀狊狅犳犈犾狔犿狌狊狊狆犲犮犻犲狊
项目
Item
自由度
N
X轴平均值
MeanX
Y轴平均值
MeanY
X轴标准差
SSx
Y轴标准差
SSy
检验结果 Manteltest
相关系数Relationcoefficient 狋值狋value 可信度Prob.randomZ<obs.Z
1 780 1.3754 0.2736 84.4311 16.8016 0.3613 3.9877 0.9986
2 780 1.3754 0.5832 84.4311 61.4440 0.2993 3.4132 0.9921
3 780 0.2736 0.5832 16.8016 61.4440 0.1489 1.6177 0.9616
4 780 1.3754 0.8808 84.4311 81.9115 0.3276 6.3988 1.0000
5 780 0.2736 0.8808 16.8016 81.9115 0.3006 5.8259 1.0000
6 780 0.5832 0.8808 61.4440 81.9115 0.1350 1.8553 0.9562
注:1代表形态(X轴)和等位酶(Y轴);2代表形态(X轴)和SSR(Y轴);3代表等位酶(X轴)和SSR(Y轴);4代表形态(X轴)和种本身(Y轴);
5代表等位酶(X轴)和种本身(Y轴);6代表SSR(X轴)和种本身(Y轴)。
Note:Number1denotesmorphologicaldata(Xaxis)andalozymedata(Yaxis);2morphologicaldata(Xaxis)andmicrosatelitedata(Yaxis);3
alozymedata(Xaxis)andmicrosatelitedata(Yaxis);4morphologicaldata(Xaxis)andspeciesdata(Yaxis);5alozymedata(Xaxis)andspeciesda
ta(Yaxis);6microsatelitedata(Xaxis)andspeciesdata(Yaxis).
97第18卷第3期 草业学报2009年
图3 过氧化物酶的电泳图谱
犉犻犵.3 犈犾犲犮狋狉狅狆犺狅狉犲狊犻狊狅犳狆犲狉狅狓犻狊狅犿犲犻狀狊犪犿狆犾犲狊狅犳犈犾狔犿狌狊狆狅狆狌犾犪狋犻狅狀狊
泳道的样品分别是Samplesinlanesare:QMZ,HYZK,QMHD,GHD,QMDL,HGZ(N),QMDP,QXED,GYP,GHY;
每个种群连续5个个体Fiveconsecutivesamplesforeachpopulation
图4 引物 犠犕犛43的扩增结果
犉犻犵.4 犈犾犲犮狋狉狅狆犺狅狉犲狊犻狊狅犳犪犿狆犾犻犳犻犲犱犫犪狀犱狊犫狔狆狉犻犿犲狉狊犠犕犛43
泳道的样品分别是Samplesinlanesfrom1to23are:1.GJY;2.GJD;3.QMDL;4.QMDC1;5.QHYF;6.HGY(N);7.GYP;8.QMZ;
9.GHD;10.QML(Y);11.QMHC;12.QMC;13.QMXC;15.QXED;18.GL(W);19.QHHC;20.GJH;21.GHY;23.GJW
酶、形态和SSR呈极显著(犘<0.01)的相关关系,等位酶和SSR显著(犘<0.05)相关。不同标记反映的亲缘关系
与种本身之间,形态学标记的遗传关系支持基于形态特征分类系统而与种本身的亲缘关系更加相符(犚=
0.3276,犘=1.0000),其次是等位酶标记(犚=0.3006,犘=1.0000),SSR标记偏离相对较远(犚=0.1350,犘=
0.9562)。
3 讨论
本试验采用3种遗传标记对我国披碱草属植物进行系统聚类分析,不同标记显示的亲缘关系有很大的差异。
形态学标记把我国9种披碱草属植物可以明显分为两大类:披碱草组和老芒麦组。以形态学特征为依据,披碱草
属共有2个组,披碱草组为颖宽长类群,属返祖进化的表现;另一个老芒麦组为颖狭小类群,这个类群很可能就是
由鹅观草属进化到猬草属的过渡类群。尤其是老芒麦,其颖特别细狭,先端渐尖成芒,而且花序也特别疏松、冗
长、蜿蜒下垂,外观近似鹅观草属植物,因而成为近代一些学者将鹅观草属归并到披碱草属中的一个重要原
因[4,5,23,24]。本研究的形态学特征中也可以找到支持广义披碱草属概念的证据,然而国内学者并不支持广义的披
碱草属概念[25]。本试验中的材料是依据我国的分类系统,即狭义的披碱草属概念,在我国仅包括12个物种,不
同的分类系统会缩小或放大不同地域分类群的亲缘关系。同时,仅披碱草属植物的穗部形态就存在很大的多样
08 ACTAPRATACULTURAESINICA(2009) Vol.18,No.3
图5 犛犛犚标记披碱草属种间亲缘关系
犉犻犵.5 犚犲犾犪狋犻狅狀狊犺犻狆犪犿狅狀犵犈犾狔犿狌狊狊狆犲犮犻犲狊犫狔犛犛犚犱犪狋犪
图6 披碱草属种群形态和等位酶标记之间的关系
犉犻犵.6 犚犲犾犪狋犻狅狀狊犺犻狆狅犳犿狅狉狆犺狅犾狅犵犻犮犪犾犪狀犱犪犾狅狕狔犿犲犱犪狋犪犪犿狅狀犵狆狅狆狌犾犪狋犻狅狀狊狅犳犈犾狔犿狌狊狊狆犲犮犻犲狊
性[26],因此采用外部形态特征对不同披碱草的亲缘关系进行研究会受到经典分类学的干扰。
虽然等位酶标记没有严格按照穗子的形状(披碱草组和老芒麦组)分为2组,但同一组内的披碱草亲缘关系
较近,同时具有相同染色体组的披碱草亲缘关系也较近。SSR标记没有把参试材料分成几个明显的组,且各种
间表现出较远的亲缘关系。就披碱草属种间关系,周永红等[27]采用RAPD技术对10种披碱草进行聚类分析,结
果是四倍体和六倍体分别聚为一类,基本与其他形态学、细胞学研究结果一致。由此可见,披碱草属植物之间的
亲缘关系密切,采取不同的遗传分子标记方法可能会取得不同的结果。在本试验中也证明了这一点,更确切地
说,形态学标记的遗传关系支持基于形态特征分类系统而与种本身的亲缘关系更加相符(犚=0.3276,犘=
1.0000),其次是等位酶标记(犚=0.3006,犘=1.0000),SS犚标记偏离相对较远(犚=0.1350,犘=0.9562)。
关于形态学与其他遗传标记之间,多数研究表明形态学性状与DNA标记在研究遗传多样性中具有显著的
相关性[28,29]。本试验中形态学距离和等位酶遗传距离之间相关性极显著,类似与Cox等[30]报道大豆(犌犾狔犮犻狀犲
犿犪狓)形态学距离和等位酶遗传距离间存在有较高的相关性(犚=0.60)。余汉勇等[31]应用形态、等位酶和SSR
18第18卷第3期 草业学报2009年
图7 披碱草属种群形态和犛犛犚标记之间的关系
犉犻犵.7 犚犲犾犪狋犻狅狀狊犺犻狆狅犳犿狅狉狆犺狅犾狅犵犻犮犪犾犪狀犱犛犛犚犱犪狋犪犪犿狅狀犵狆狅狆狌犾犪狋犻狅狀狊狅犳犈犾狔犿狌狊狊狆犲犮犻犲狊
图8 披碱草属种群等位酶和犛犛犚标记之间的关系
犉犻犵.8 犚犲犾犪狋犻狅狀狊犺犻狆狅犳犪犾狅狕狔犿犲犪狀犱犛犛犚犱犪狋犪犪犿狅狀犵狆狅狆狌犾犪狋犻狅狀狊狅犳犈犾狔犿狌狊狊狆犲犮犻犲狊
标记研究水稻(犗狉狔狕犪狊犪狋犻狏犪)矮仔占衍生品种的遗传差异,3类标记的遗传距离间虽存在有一定的相关性,但相
关系数较小(0.217~0.494)。其中,形态学性状与SSR标记间相关极显著(犚=0.494,犘=0.002),等位酶与其
他2类标记的相关性均不显著,表明3类标记性状的分析结果存在较大的差异。在等位酶与DNA标记间的相
关性研究中,不同分子标记有很大的不同。Buso等[32]在研究南美展颖野生稻4个种群的遗传结构中发现等位
酶和RAPD有较好的一致性,但RAPD在解决种间系统关系方面显得有点微弱[33]。
在不同的标记之间有着不同的结论,Sun等[34]报道在披碱草属不同种的20个种群中RAPD与SSR的犚=
28 ACTAPRATACULTURAESINICA(2009) Vol.18,No.3
0.80。在本试验中发现,形态、等位酶和SSR标记得到的数据之间虽然存在着明显的相关关系,但相关系数很
低,形态和等位酶之间相关系数最高,为0.3613(犘<0.01),等位酶和SSR之间最低,为0.1489(犘<0.05),形
态和SSR之间为0.2993(犘<0.01)居中。这与Thormann等[35]和Sun等[34]得出的不同标记之间的相关系数出
入较大。Sun等[36]用 Mantel的检验对其关系进行分析,发现存在弱的相关关系,RAPD与SSR仅为0.267,
RAPD与同工酶为0.204,SSR与同工酶为0.164的结果与本试验相似。Roberto等[37]在对小麦和其近缘
种———犬草(1种四倍体披碱草属植物)进行分子标记的比较研究发现,RAPD和同工酶之间相关性最高,为
0.93,RAPD和SSR之间为0.74,SSR和同工酶之间相关系数最低为0.68。至于和种本身之间的关系,也是
SSR最低,为0.72。虽然和本试验相关系数有一定的差距,但 Mantel检验结果是相同的,均处于显著相关水平,
相关性的强弱次序SSR也是最弱的。
对在这些树状图中发现的差异最科学的解释就是基于不同遗传标记所提供的遗传信息的不同,形态学特征
是植物最直接的外部体现,同工酶(等位酶)变异仅反映编码蛋白的基因之间差异。编码顺序是在巨大的选择压
力下维持其功能序列的。而重复序列通过扩增和漂移的方式使进化速度快于单烤贝片断。考虑到双向比较间表
现高水平多样性,因此SSR检测结果与其他检测方法之间缺乏相关性是有可能的[38]。由于形态学性状具有选
择性和易受环境影响,等位酶和SSR标记的选择中性和不受环境影响特征,造成3种标记在聚类分析中的差异。
另外,遗传标记技术中影响遗传关系的是分析中所用到的标记物或探针数量和数据统计的方法[39,40]。通常,当
分子标记物或探针数量增加时,精确度也随之增加。但形态学指标却不尽相同,不同形态指标有时可能存在线性
化掩盖或加强了某个方面的特征。不同的数据计算和统计方法也会对结果有一定的影响,表示亲缘关系的有:遗
传相似性系数、Jaccard相似性系数,还有遗传距离等,统计方法有采用线性化的,有采用标准化的,还有按照质量
性状(qualitativetraits)进行计算相似性的等。目前研究中所使用的数据计算和统计方法以及标记方法,采用了
不同数量的标记物,并产生不同数量的片断,也许是同工酶、SSR数据间相关系数不同的部分原因。鉴于以上因
素,由这些方法产生的系统发育树中的不一致并不令人感到惊奇。
不同分子标记在处理不同类群之间的系统关系适宜程度有所不同。以剪股颖属为例,RAPD、SSR及ISSR
(intersimplesequencerepeat)标记虽然都产生了谱带清晰、多态性程度高的现象,可以利用这3种分子标记进行
种间的系统学研究完全可行。当研究种下关系时,AFLP(amplifiedfragmentlengthpolymorphism)标记可提供
更多和精确的多态性状,当涉及到种特异性分析时,SCAR(sequencecharacterizedamplifiedregion)标记是首
选[41]。郭海林等[42]认为SSR更适用于属内不同品系的指纹图谱的构建,刘伟等[43]认为ISSR标记比较适合于
植物材料的来源分析。用哪种遗传标记是评价披碱草属(狭义)系统发育关系最好的方法,形态是反应披碱草属
植物变异最稳定也是比较可靠的特征,形态、等位酶和微卫星分析技术都是应用于种、种群遗传变异研究常用方
法,且2种分子标记均是共显性标记,以孟德尔方式遗传,还能够定量分析研究材料中的遗传差异,等位酶不容易
受环境的饰变,而微卫星却反应了环境的饰变作用。二者反应基因组不同的遗传变异,结合形态学,基本能够全
面反应基因组的遗传变异信息。结合本研究,可认为形态学是最直接、最能反映披碱草属植物本身特点的标记方
法,但要考虑指标之间的线性化,并非越多越好;等位酶分析是一种研究种群内部遗传结构和变异比较好的方法,
也是种、种群间遗传关系有价值的预测者;尽管SSR在属或种间分析它并不是最适用的,由于SSR是最多态的
标记,对于种内分析非常适合。
参考文献:
[1] 卢宝荣.披碱草属(犈犾狔犿狌狊L.)两个物种之订正[J].植物分类学报,2002,40(6):539545.
[2] BenthamG.NotesonGramineae[J].JournaloftheLinneanSociety.Botany.London,1881,(1):14134.
[3] NevskiSA.HordeaeBenth[A].In:KomarovVL,RoshevitsRY,ShishkinBK.FloraURSSII[M].Leningrad:FloraUS
RR,1934.
[4] DeweyDR.ThegenomicsystemofclassificationasaguidetointergenerichybridizationwiththeperennialTriticeae[A].In:
GustafsenJP.GeneManipulationinPlantImprovement[M].NewYork:PlenumPress,1984.209279.
38第18卷第3期 草业学报2009年
[5] Lve?.ConspectusoftheTriticeae[J].FeddesRepert,1984,95(78):425521.
[6] 蔡联炳.中国鹅观草属的分类研究[J].植物分类学报,1997,35(2):148177.
[7] 郭本兆.中国植物志[M].北京:科学出版社,1987.7104.
[8] BaumBR,YenC,YangJL.Roegneria:Itsgeneticlimitsandjustificationforitsrecognition[J].CanadianJournalofBotany,
1991,69:282294.
[9] 卢宝荣.披碱草属与大麦属系统关系的研究[J].植物分类学报,1997,35(3):193207.
[10] 严学兵,郭玉霞,周禾,等.我国披碱草属植物分类和遗传多样性研究现状[J].草业科学,2005,22(7):17.
[11] 张颖,周永红,张利,等.鹅观草属、披碱草属、猬草属和仲彬草属物种的RAMP分析及其系统学意义[J].西北植物学报,
2005,25(2):368375.
[12] 张颖,周永红,张利,等.小麦族鹅观草属、披碱草属、猬草属和仲彬草属细胞质基因组PCR-RFLP分析[J].遗传,2006,
28(4):449457.
[13] 周永红,郑有良,杨俊良,等.鹅观草属、披碱草属、猬草属和仲彬草属RAPD分析及其系统学意义[J].四川农业大学学报,
2001,19(1):1421.
[14] 王中仁.植物等位酶分析[M].北京:科学出版社,1996.1011.
[15] DoyleJ.IsolationofplantDNAfromfreshtissue[J].FocusLifeTechnologies,1989,12:1.
[16] SunGL,SalomonB,BothmerRV.Characterizationandanalysisofmicrosatelitelociin犈犾狔犿狌狊犮犪狀犻狀狌狊(犜狉犻狋犻犮犲犪犲:犘狅犪犮犲
犪犲)[J].TheoreticalandAppliedGenetics,1998,96:676682.
[17] SunGL,SalomonB,BothmerVR.Characterizationofmicrosatelitelocifrom犈犾狔犿狌狊犪犾犪狊犽犪狀狌狊andlengthpolymorphism
inseveral犈犾狔犿狌狊species(犜狉犻狋犻犮犲犪犲:犘狅犪犮犲犪犲)[J].Genome,1998,41:455463.
[18] RderMS,KorzunV,WendehakeK,犲狋犪犾.Amicrosatelitemapofwheat[J].Genetics,1998,149:20072023.
[19] PlaschkeJ,GanalM W,RderMS.Detectionofgeneticdiversityincloselyrelatedbreadwheatusingmicrosatelitemarkers
[J].TheoreticalandAppliedGenetics,1995,91:10011007.
[20] NeiM.Estimationofaverageheterozygosityandgeneticdistancefromasmalnumberofindividuals[J].Genetics,1978,89:
583590.
[21] RohlfFJ.NumericalTaxonomyandMultivariateAnalysisSystem[M].Version1.80.ExeterSoftware,NY:Setauket,
1993.
[22] MantelN.Thedetectionofdiseaseclusteringandageneralizedregressionapproach[J].CancerResearch,1967,27:209
220.
[23] TzvelevNN.Tribe3TriticeaeDum[A].PoaceaeUSSR[M].Leningrad:NaukaPubHouse,1976.105206.
[24] 蔡联炳.根据外部形态特征试论鹅观草属的亲缘演化关系[J].西北植物学报,1998,18(4):606612.
[25] 肖海峻,翟利剑,卢宏双,等.小麦族鹅观草属植物研究进展[J].草业科学,2007,24(4):4146.
[26] 鄢家俊,白史且,马啸,等.川西北高原野生老芒麦居群穗部形态多样性研究[J].草业学报,2007,16(6):99106.
[27] 周永红,郑有良,杨俊良,等.10种披碱草属植物的RAPD分析及其系统学意义[J].植物分类学报,1999,37(5):425
432.
[28] BeerSC,GoffredaJ,PhilipsTD,犲狋犪犾.Assessmentofgeneticvariationin犃狏犲狀犪狊狋犲狉犻犾犻狊usingmorphologicaltraits,
isozymes,andRFLPs[J].CropScience,1993,33:13861393.
[29] NbrilaJ,CrnokrakP.Comparisonofgeneticdifferentiationatmarkerlociandquantitativetraits[J].JournalofEvolutionary
Biology,2001,14:892903.
[30] CoxTS,MangYT,GormanMB,犲狋犪犾.Relationshipbetweencoefficientofparentageandgeneticsimilarityindicesinthe
soybean[J].CropScience,1985,25:529532.
[31] 余汉勇,魏兴华,王一平,等.应用形态、等位酶和SSR标记研究水稻矮仔占衍生品种的遗传差异[J].中国水稻科学,
2004,18(6):477482.
[32] BusoGSC,RangelPH,FerreiraME.AnalysisofgeneticvariabilityofSouthAmericanwildricepopulations(犗狉狔狕犪
犵犾狌狉狀狀犲狆犪狋狌犾犪)withisozymesandRAPDmarkers[J].MolecularEcology,1998,7:107117.
[33] 李杰勤,王丽华,詹秋文,等.应用RAPD标记对高粱属两物种之间遗传差异的研究[J].草业学报,2007,16(5):140144.
48 ACTAPRATACULTURAESINICA(2009) Vol.18,No.3
[34] SunGL,SalomonB,BothmerRV.Analysisoftetraploid犈犾狔犿狌狊speciesusingwheatmicrosatelitemarkersandRAPD
makers[J].Genome,1997,40:806814.
[35] ThormannCE,FerreiraME,CamagoLE.ComparisonofRFLPandRAPDmarkerstoestimatinggeneticrelationships
withinandamongcruciferousspecies[J].TheoreticalandAppliedGenetics,1994,88:973980.
[36] SunGL,DíazO,SalomonB,犲狋犪犾.Geneticdiversityin犈犾狔犿狌狊犮犪狀犻狀狌狊asrevealedbyisozyme,RAPDandmicrosatelite
markers[J].Genome,1999,42:420431.
[37] RobertoG,DessislavaSB,FrancoisF.Specificgeneticmarkersforwheat,spelt,andfourwildrelatives:Comparisonof
isozymes,RAPDs,andwheatmicrosatelites[J].Genome,2001,44(4):610621.
[38] RusselJR,FulerJD,MacaulayM,犲狋犪犾.Directcomparisonoflevelsofgeneticvariationamongbarleyaccessionsdetected
byRFLPs,AFLPs,SSRsandRAPDs[J].TheoreticalandAppliedGenetics,1997,95:714722.
[39] MessmerM M,MelchingerAE,WoodmanWL,犲狋犪犾.GeneticdiversityamongprogenitorsandelitelinesfromtheIowa
StiffStalkSynthetic(BSSS)maizepopulation:ComparisonofalozymeandRFLPdata[J].TheoreticalandAppliedGenetics,
1991,83:97107.
[40] SmithOS,SmithJSC,BowenSL,犲狋犪犾.NumbersofRFLPprobesnecessarytoshowassociationsbetweenlines[J].Maize
Genet.Coop.Newslett,1992,66:66.
[41] 张道远,TakamizoT.剪股颖属植物遗传分化及系统关系的分子标记研究[J].草业学报,2006,15(3):100106.
[42] 郭海林,刘建秀,高鹤,等.结缕草属优良品系SSR指纹图谱的构建[J].草业学报,2007,16(2):5359.
[43] 刘伟,张新全,李芳,等.西南区野生狗牙根遗传多样性的ISSR标记与地理来源分析[J].草业学报,2007,16(3):5561.
犜犺犲狆犺狔犾狅犵犲狀犲狋犻犮狉犲犾犪狋犻狅狀狊犺犻狆狊犫犲狋狑犲犲狀狀犻狀犲犈犾狔犿狌狊狋犪狓犪犳狉狅犿犆犺犻狀犪犫犪狊犲犱狅狀狋犺狉犲犲犵犲狀犲狋犻犮犿犪狉犽犲狉狊
YANXuebing1,ZHOUHe2,WANGKun2,WANGChengzhang1,GUOYuxia1
(1.SchoolofAnimalScienceandVetinaryofHenanAgriculturalUniversity,Zhengzhou450002,China;
2.GrasslandInstituteofChinaAgriculturalUniversity,Beijing100094,China)
犃犫狊狋狉犪犮狋:Geneticmarkersarecommonlyusedtoconstructsystematictreesinstudiesofphylogeneticrelation
shipsamongplanttaxa.Inthisstudy,ninenative犈犾狔犿狌狊speciesweresampledtoanalyzetheirphylogenetic
relationshipsusingthreegeneticmarkers.Nine犈犾狔犿狌狊specieswereseparatedintotwogroups:Sect.犜狌狉犮狕犪
狀犻狀狅狏犻犪(Nevski)Tzvel.andSect.犈犾狔犿狌狊.Alozymemarkersdidnotshowcleargroupings,but犈犾狔犿狌狊spe
ciesbelongingtothesamesectionorwiththesamegenomicconstitutionindicatedacloserelationship.Microsa
telitemarkersdidnotindicatecleargroupingsorrelationshipsamongthenine犈犾狔犿狌狊taxa.Asignificantrela
tionshipwasfoundbetweenthreegeneticmarkers.Amorphologicalmarkerwasthemostsignificantandgave
similarresultstotheactualrelationshipofindividualspecies.Thesecondwasthealozymemarkerandthethird
amicrosatelitemarker.Therefore,wesuggestthatmorphologicaltraitsarethemostdirectandefficient
geneticmarkersof犈犾狔犿狌狊taxa,whilealozymemarkersareavaluablepredictoroftherelationshipamongpop
ulationsorspecies,andmicrosatelitesaremoresuitableforstudyinginterpopulationratherthaninterspecies
relationshipof犈犾狔犿狌狊.
犓犲狔狑狅狉犱狊:犈犾狔犿狌狊;morphology;alozyme;microsatelite;relationship
58第18卷第3期 草业学报2009年