全 文 :书中国西南区鸭茅种质遗传变异的犛犛犚分析
谢文刚1,张新全1,马啸1,黄琳凯1,曾兵2
(1.四川农业大学草业科学系,四川 雅安625014;2.西南大学动物科学系,重庆402460)
摘要:用SSR标记对来自中国西南地区(四川、云南、贵州、重庆)的11份鸭茅种质的110个单株进行遗传多样性研
究。结果表明,1)21对引物共扩增出多态性条带143条,平均每对引物扩增出6.8条,多态性条带比率 (犘犘犅)达
90.7%;多态性信息含量(犘犐犆)范围为0.17(A01F24)~0.45(A01E02),平均为0.33;Shannon多样性指数(犐)范围
在0.0463~0.3477,表明供试鸭茅种质具有丰富的遗传多样性。2)AMOVA分析表明,63.93%的遗传变异存在
于种质内部,种质间变异仅为21.93%,二倍体与四倍体鸭茅群体间遗传差异不显著。3)对各地理类群基于Nei氏
无偏估计的遗传一致度可将11份鸭茅分为2大类:来自重庆和贵州的3份二倍体材料 YA20063、YA20064和
YA02101聚为I类,其余8份材料聚为II类,来自相同或相似生态地理环境、相同倍型的材料基本聚为一类,呈现
出一定的相关性。
关键词:鸭茅;种质资源;SSR;遗传多样性
中图分类号:S543+.303;Q943 文献标识码:A 文章编号:10045759(2009)04013809
鸭茅(犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪)又名鸡脚草、果园草,属于禾本科,鸭茅属(犇犪犮狋狔犾犻狊)全属共约5个种,主要分布于
欧亚温带和北非[1]。鸭茅是世界著名的优良牧草,在北美、欧洲和大洋洲等地广泛种植。自然界的鸭茅存在四倍
体(2狀=4x=28)和二倍体(2狀=2x=14)2种类型[2,3]。我国野生鸭茅资源丰富,生态地理类型多样,蕴藏着丰富
的基因资源,全国发现野生鸭茅生长地26个[4~7]。近年来,鸭茅在我国的四川、重庆、贵州、山西、甘肃等地被大
量应用于草地畜牧业及生态环境治理,取得了良好的经济和生态效益。
国内外学者已通过形态标记、细胞学标记和同工酶生化标记等,对鸭茅居群遗传结构与生态型遗传分化及种
质间遗传多样性进行研究[8~10]。分子标记技术因其简单、快捷、受环境影响小和可重复性强等优点,被广泛应用
在牧草遗传研究中。目前,国内已采用AFLP(amplificationfragmentlengthpolymorphism,扩增片段长度多态
性)、RAPD(randomamplifiedpolymorphicDNA,随机扩增多态性DNA)、ISSR(intersimplesequencerepeat,
简单重复序列间多态性)、SRAP(sequencerelatedamplifiedpolymorphism,相关序列扩增多态性)等分子标记对
鸭茅进行了种质间遗传多样性研究,并揭示出自然界鸭茅资源存在着丰富的遗传变异[11~14]。但有关中国西南地
区野生鸭茅种质内遗传分化和变异的研究未见报道。
简单序列重复(simplesequencerepeat,简称SSR)又称微卫星DNA,是建立在PCR(polymerasechainreac
tion,聚合酶链式反应)基础上的新型DNA指纹标记,以1~6个碱基为基本单元的串联重复序列,在植物基因组
中广泛分布[15]。由于SSR标记具有共显性、遗传方式简单、多态性好、操作简便、重复性好、稳定可靠等特点,广
泛应用于遗传图谱构建,种质鉴定和分子标记辅助育种等研究领域。目前已开展SSR标记研究的作物有水稻
(犗狉狔狕犪狊犪狋犻狏犪)、小麦(犜狉犻狋犻犮狌犿犪犲狊狋犻狏狌犿)、玉米(犣犲犪犿犪狔狊)等。而在草业上SSR应用不多,仅见于披碱草(犈犾狔
犿狌狊犱犪犺狌狉犻犮狌狊)、柱花草(犛狋狔犾狅狊犪狀狋犺犲狊spp.)、紫花苜蓿(犕犲犱犻犮犪犵狅狊犪狋犻狏犪)、结缕草(犣狅狔狊犻犪犼犪狆狅狀犻犮犪)等少数草
种[16~19]。
本研究以我国西南区11份鸭茅种质的110个单株为材料。采用SSR分子标记技术研究鸭茅种质遗传多样
性,为科学保护和合理利用鸭茅种质资源及新品种的培育提供理论依据。
138-146
2009年8月
草 业 学 报
ACTAPRATACULTURAESINICA
第18卷 第4期
Vol.18,No.4
收稿日期:20080909;改回日期:20081020
基金项目:国家“973”基础研究项目(2007CB108907)和国家自然基金项目(30871820)资助。
作者简介:谢文刚(1982),男,四川资中人,在读博士。Email:xwgorchardgrass@hotmail.com
通讯作者。Email:zhangxq@sicau.edu.cn
1 材料与方法
1.1 供试材料
研究所用的11份野生鸭茅种质为四川农业大学草业科学系从我国西南地区收集而得,其中5份来自四川,2
份来自云南,2份来自贵州,2份来自重庆,经染色体检测,其中7份为二倍体,4份为四倍体(表1)。每份种质选
取种子在温室发芽,2周后随机选取10个单株,11份种质共110个单株,并于2008年3月将其种于四川农业大
学教学试验基地。所用的100对鸭茅SSR引物序列由日本草地畜产种子协会饲料作物研究所才宏伟教授提供,
由上海生工生物技术有限公司合成。PCR所用的 Mg2+,Taq酶,dNTP,Buffer均购自博瑞克生物技术公司。
表1 供试材料
犜犪犫犾犲1 犇.犵犾狅犿犲狉犪狋犪犪犮犮犲狊狊犻狅狀狊狋犲狊狋犲犱犻狀狋犺犻狊狊狋狌犱狔
编号Accessioncode 来源地Origin 地理类群Geographicalgroups 倍性Chromosomeno.(2狀) 取样数Samplesize
YA90130 四川,茂县 Maocounty,Sichuan SC 14 10
YA912 四川,宝兴Baoxing,Sichuan SC 28 10
YA9070 四川,康定 Kangding,Sichuan SC 14 10
YA917 四川,汉源 Hanyuan,Sichuan SC 14 10
YA02106 四川,宝兴Baoxing,Sichuan SC 28 10
YA02111 云南,中甸Zhongdian,Yunnan YN 14 10
YA02116 云南,昆明 Kunming,Yunnan YN 28 10
YA20063 重庆,巫溪 Wuxi,Chongqing CQ 14 10
YA20064 重庆,巫山 Wushan,Chongqing CQ 14 10
YA01101 贵州,毕节Bijie,Guizhou GZ 28 10
YA02101 贵州,毕节Bijie,Guizhou GZ 14 10
注:SC,四川;YN,云南;CQ,重庆;GZ,贵洲。
Note:SC,Sichuan;YN,Yunnan;CQ,Chongqing;GZ,Guizhou.
1.2 基因组DNA提取
取生长良好的鸭茅幼嫩植物的叶片,采用SaghaiMaroof等[15]的CTAB(cetyltrimethylammoniumbro
mide,十六烷基三甲基溴化铵)方法提取DNA。通过0.8%琼脂糖凝胶电泳和紫外分光光度计检测DNA的浓
度和纯度,合格的样品于4℃冰箱保存备用。
1.3 SSR反应体系的建立和优化
利用正交设计L16(44)与单因素相结合的方法[20],对 Mg2+浓度、dNTP、引物浓度、Taq酶浓度进行4因素4
水平的正交试验,并在ThermoHybaidPCR仪上进行引物最佳退火温度的筛选。
1.4 SSR引物的筛选
选取经电泳检测质量较高且田间试验性状差异较大的材料DNA (分别为SAU00850,SAU02123,“宝兴”,
“安巴”)作模板,从100对引物中筛选出多态性高的引物,再对全部材料进行PCR扩增。
1.5 PCR扩增及电泳
PCR扩增在ThermoHybaidPCR仪上进行。扩增产物在6%聚丙烯酰胺凝胶上(丙烯酰胺∶甲叉=19∶1,
7.5mol/L尿素,1×TBE缓冲液)进行检测。样孔中每样品点样10μL,以博瑞克公司的D1500Marker为对
照,先在北京六一电泳仪厂的DYY6C电泳仪上进行200V恒定电压下30min的预电泳,然后在400V恒定电
压下电泳2h,再用0.1%的AgNO3进行银染色和在NaOH 液中显色,凝胶在灯光下用数码相机照相保存以供
分析。
1.6 数据统计及分析
对获得的清晰可重复的DNA条带进行统计,在相同迁移位置上有带的记为1,无带的记为0,构成遗传相似
矩阵。
931第18卷第4期 草业学报2009年
根据表征矩阵,统计SSR扩增产物的条带总数和多态性条带数量,计算多态性条带所占的比率(犘犘犅)和引
物多态性信息含量(犘犐犆)[21],犘犐犆犻=2犳犻(1-犳犻),式中,犘犐犆犻表示引物犻的多态性信息含量,犳犻表示有带所占的
频率,1-犳犻表示无带所占的频率。在哈-温(Hardy-Weinberg)平衡的基础上,计算各种质材料Nei基因多样
性指数(犎、犐)[22]。AMOVA分析计算种质内和种质间及地理内和地理间的变异分布[23]。
先采用POPGENE软件计算各材料Nei氏无偏估计的遗传一致度和遗传距离,再用NTSYSpc绘制 UPG
MA聚类图[24]。
以上参数计算分别在POPGENE1.31、AMOVA-PREP1.01、DCFA1.1中进行。
2 结果与分析
2.1 筛选引物及SSR扩增
从100对引物中筛选出21对条带清晰、多态性好的引物(表2),用于110份单株材料的PCR扩增(图1)。
正交优化试验所得适合鸭茅SSR-PCR反应的扩增程序为:94℃预变性4min,94℃变性30s,48~52℃变
性30s,72℃延伸1min,共35个循环,72℃10min,4℃保存。所得15μL优化反应体系为:模板 DNA10
ng/μL,dNTP240μmol/L,引物浓度0.4μmol/L,Taq酶1.0U,Mg
2+2.5mmol/L。
2.2 SSR多态性分析
21对引物共扩增出156条清晰可辨条带,平均每对引物扩增的条带为7.4条,最多能得到13条清晰带
(A03K22),最少有3条(A01C20);多态性条带总数为143条,平均每对引物扩增6.8条,引物的平均多态性比率
为90.7%;多态性信息含量范围为0.17(A01F24)~0.45(A01E02),平均为0.33(表3),表明SSR能检测到较多
的遗传位点,获得多态性较好的PCR结果。
表2 用于鸭茅犛犛犚分析的引物序列
犜犪犫犾犲2 犘狉犻犿犲狉狊犲狇狌犲狀犮犲狊狌狊犲犱犻狀犛犛犚犪狀犪犾狔狊犲狊狅犳犇.犵犾狅犿犲狉犪狋犪
引物编号Primercode 引物序列(5′-3′)Primersequence 引物编号Primercode 引物序列(5′-3′)Primersequence
A01B10
TCTTCCTTGGAAAACATCAA
ACTTGCTTACACGGTATCATG
A03C05
TAAGAATCGATCCTCCCG
ACCTTCTTCCACTCCGTC
A01C20
AACCAGTTGTTGATCTGCTT
TACGATATGTGCGTGTTGAT
A03K22
AGACTCTAGGGTGGCACAC
GTAGCACGCTAACGAGAGAT
A01E02
AGCTGTGGAGAAAAAAATGA
GATGCCATTAAGTTCAAAATG
A03N16
AACGCAACGTCTATGGTTAT
ACACGCACCCACACATAC
A01E14
ACCCGTTTTCTATCTCCAG
GTTCTAGCGTCGTGAGGG
A04C24
AGCAACATATCTTACTGCAATG
ATCAAACTCGAAAAGTTGTCA
A01F24
AAAATGTTTTATTCTCAGCCC
TGCAAGATGGAATGCTCT
A04O08
AGAGGTTAGATGGATGTAGGC
ATAGACCCATAGCATGTTGG
A01I13
ATGCTGTTTGATCACAGTCA
GTTGGACTGCCATTACTAGC
B01A02
TTCTCCATTAAGCCTCCAG
CAGCGAGTCCTTGTCGAC
A01K14
AAGGATGGCCTGATCTTC
GCAGAGGTCTTTCTCTTGG
B01A05
GAGAGCGGCAGAGTTATTC
AAAGGTCGATATCTCTATTCCA
A02A10
AGGTTACCGATAGTAAGTGGG
AGGGGATGGTTGGTTAGTAT
B01B19
AGAAGTTGGCCTGTCTCC
CTCTTCCTTCCTCCTTGG
A02G09
TACACGAGAGGGCAGATACT
CGTAACTTGAATCTTCCAGG
B01C15
GTCGATTGATGGTGTACGTA
TCTAGTGCTACTTGTATGCACC
A02I05
GCAAATGTCCACACCATT
CTACCACAGCGACATCAAG
B01E09
ACAACTCACAAACTCAAGAACA
GTGGACTCGGAGGAGAAG
A03B16
TCTGGAATCTCTCTGAAATCA
ATCTTGACCCTGATGTTCTG
041 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4
图1 引物犃03犆05对犢犃90130和犢犃917的犛犛犚扩增结果
犉犻犵.1 犛犛犚犪犿狆犾犻犳犻犲犱狉犲狊狌犾狋狊狑犻狋犺狆狉犻犿犲狉狆犪犻狉犖狅.犃03犆05狅狀犢犃90130犪狀犱犢犃917
编号1~10为YA90130,11~20为YA917,M为 MarkerCode1-10forYA90130,11-20forYA917,M:Marker
表3 鸭茅犛犛犚引物的多态性
犜犪犫犾犲3 犃犿狆犾犻犳犻犮犪狋犻狅狀狉犲狊狌犾狋狊犳狉狅犿21狆狉犻犿犲狉狆犪犻狉狊
引物编号
Primercode
条带总数
Totalnumbers
ofbands
多态性条带数
Numbersofpolymorphism
bands
多态性比率
Percentageofpolymorphism
bands(犘犘犅,%)
多态信息含量
Polymorphisminformation
content(犘犐犆)
A01B10 5 4 80.0 0.38
A01C20 3 3 100.0 0.22
A01E02 7 6 85.7 0.45
A01E14 9 8 88.9 0.37
A01F24 5 4 80.0 0.17
A01I13 6 6 100.0 0.42
A01K14 12 12 100.0 0.32
A02A10 10 10 100.0 0.38
A02G09 6 6 100.0 0.28
A02I05 10 9 90.0 0.33
A03B16 5 4 80.0 0.39
A03C05 8 7 87.5 0.34
A03K22 13 13 100.0 0.29
A03N16 9 8 88.9 0.39
A04C24 4 4 100.0 0.26
A04O08 5 4 80.0 0.34
B01A02 8 7 87.5 0.31
B01A05 5 4 80.0 0.35
B01B19 8 7 87.5 0.25
B01C15 9 9 100.0 0.33
B01E09 9 8 88.9 0.38
总共 Total 156 143
平均 Mean 7.4 6.8 90.7 0.33
2.3 鸭茅种质遗传变异分析
2.3.1 11份鸭茅种质遗传距离分析 对扩增结果采用NeiLi遗传一致度和遗传距离的计算方法,得到供试材
料相似性矩阵。11份种质的遗传一致度变异范围在0.7320~0.9382,平均为0.8351;遗传距离范围在0.0473
141第18卷第4期 草业学报2009年
~0.3120,平均为0.1796。由遗传相似矩阵可以看出,来自四川宝兴的YA02106和YA912遗传一致度最大,
遗传距离最小,表明他们的亲缘关系较近,它和重庆巫溪的二倍体材料YA20063遗传距离最大,达到0.3120,
两者之间亲缘关系最远。
2.3.2 鸭茅种质遗传多样性分析 从Nei氏遗传多样性指数来看,在倍型水平上,二倍体(犐,0.3153;犘犘犅,
91.25%)和四倍体类群(犐,0.4177;犘犘犅,94.23%)的遗传多样性相似(表4和5),四倍体种质间遗传多样性水
平相似性高(犘犘犅在68.59%~71.79%),二倍体种质
间存在较大遗传差异(犘犘犅 在10.25%~67.95%)。
从犘犘犅、犎、犐3个遗传多样性指数来看,除重庆和贵
州的3份二倍体种质(YA20063,YA20064,YA02
101)遗传多样性较低外,其余8份材料种质内遗传多
样性水平较高,多态性百分率均大于60%。3个遗传
多样性参数犘犘犅、犎 和犐均表现为种质YA01101>
YA912>YA02106>YA02116>YA02111>YA
917>YA9070>YA90130(表4),以上结果表明,
供试材料间差异明显,具有较为丰富的遗传多样性,四
倍体遗传多样性水平总体上高于二倍体。
AMOVA分析结果表明,63.93%的遗传变异存
在于种质内部,种质间变异为21.93%,二倍体与四倍
体鸭茅群体间遗传差异不显著,变异为14.13%。从
材料来源的地理分布看,地理间的变异为24.42%,地
理内的变异为75.58%。以上结果表明,鸭茅种质的
遗传变异主要存在于种质内部,这可能是由于鸭茅是
一种异花授粉植物,自交很少发生,同时他们广泛分布
于不同的生态地理区域,长期的自然选择使其保持了
较大的变异度和丰富的多样性(表6)。
表4 鸭茅种质遗传多样性指数
犜犪犫犾犲4 犌犲狀犲狋犻犮狆狅犾狔犿狅狉狆犺犻狊犿犻狀犱犲狓犲狊犳狅狉11
犵犲狉犿狆犾犪狊犿犪犮犮犲狊狊犻狅狀狊狅犳犇.犵犾狅犿犲狉犪狋犪
编号Accessioncode 犎 犐 犘犅 犘犘犅(%)
YA90130 0.1913 0.2965 100 64.10
YA912 0.2189 0.3375 110 70.51
YA9070 0.2084 0.3198 101 64.74
YA917 0.2102 0.3210 103 66.03
YA02106 0.2171 0.3359 108 69.23
YA02111 0.2119 0.3218 106 67.95
YA02116 0.2123 0.3332 107 68.59
YA20063 0.0315 0.0463 16 10.25
YA20064 0.0577 0.0921 27 17.30
YA01101 0.2289 0.3477 112 71.79
YA02101 0.0548 0.0822 25 16.03
犎:Nei氏遗传多样性;犐:香农信息多样性指数;犘犅:多态性条带数;
犘犘犅:多态性条带比率。
犎:Nei’sgenediversity;犐:Shannon’sinformationindexofdiversity;
犘犅:Numberofthepolymorphismbands;犘犘犅:Percentageofpolymor
phismbands。
表5 鸭茅不同倍性和地理类群遗传多样性指数
犜犪犫犾犲5 犌犲狀犲狋犻犮犱犻狏犲狉狊犻狋狔狑犻狋犺犻狀狆犾狅犻犱狔犾犲狏犲犾犪狀犱犵犲狅犵狉犪狆犺犻犮犪犾狉犲犵犻狅狀狊狅犳犇.犵犾狅犿犲狉犪狋犪
类群
Groups
多态性条带数Numbersof
polymorphismbands
香农信息多样性指数Shannon’s
informationindexofdiversity(犐)
多态性比率Percentageof
polymorphismbands(犘犘犅,%)
四倍体Tetraploid 146 0.4177 94.23
二倍体Diploid 142 0.3153 91.25
四川Sichuan(SC) 148 0.3927 94.87
云南Yunnan(YN) 129 0.3846 83.33
贵州Guizhou(GZ) 105 0.2679 67.31
重庆Chongqing(CQ) 31 0.2679 19.87
2.4 供试材料基于NeiLi遗传相似系数的聚类分析
利用POPGENE软件计算出各类群间的Nei氏遗传一致度和遗传距离的无偏估计值,然后用NTSYSpc绘
制UPGMA聚类图。在相似系数0.86水平上可将11份材料分为2大类(图2),来自重庆的 YA20063,
YA20064和YA02101三份二倍体与其余材料差异明显,花序下垂、前期生长速度缓慢、分蘖能力弱,牧草产量
低,完成生育期所需时间较长、种子体积小、千粒重轻、生态适应能力较弱,植株低矮,自成I类;其余来自四川,云
241 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4
南,贵州的材料聚为II类。该类较混杂,包括来自四川的5份材料,云南的2份材料和贵州的1份材料,在形态特
征方面该类大部分种质形态及农艺性状表现为花序直立、前期生长速度快、分蘖能力强、牧草产量高、完成生育期
所需时间较短、种子较大、千粒重较重、生态适应能力较强。该大类中相同倍型的材料(如YA9070、YA90130、
YA02111)优先聚为一类,然后相同或相近地理来源的材料基本能聚为一类。图2揭示了供试材料的聚类与其
形态特征、地理分布和倍型存在一定的相关性。
表6 11份鸭茅种质遗传分化犃犕犗犞犃分析
犜犪犫犾犲6 犃狀犪犾狔狊犻狊狅犳犿狅犾犲犮狌犾犪狉狏犪狉犻犪狀犮犲(犃犕犗犞犃)犳狅狉11犪犮犮犲狊狊犻狅狀狊狅犳犇.犵犾狅犿犲狉犪狋犪
变异来源Sourceofvariation 自由度DF 总体平方差SSD 平均平方方差 MSD 方差分量Variancecomponent 总变异Totalvariance(%)
倍性水平间Betweenploidylevel 1 310.81 310.81 3.93 14.13
种质间Amongaccessions 9 867.08 96.34 8.05 21.93
种质内 Withinaccession 99 1568.20 15.84 15.84 63.93
地区间Amongregions 7 988.18 141.17 6.77 24.42
地区内 Withinregion 142 2977.31 20.96 20.96 75.58
按材料地理来源分可将4个地理类群大致聚为3
类。4个地理类群间的遗传一致度较高,均大于0.824
(表7)。其中云南类群和四川类群遗传距离最小
(0.0350),表现出较高的遗传相似性,最先聚为一类,
再和贵州类群聚为一类,重庆类群与其他类群遗传距
离较远,单独聚为一类(图3)。
3 讨论与结论
3.1 鉴定分析材料的遗传多样性是研究物种起源进
化、发现新的基因资源、改良现有育种材料的基础工
作。获取育种资源的遗传多样性信息方法很多,如形
态鉴定、生化测试等。从前人研究结果看,这些方法存
在着数量有限、位点不够等不足之处,难以获得完整的
数据资料,从而难以准确反映材料间实质性遗传差异。
表7 各地理类群间犖犲犻氏遗传一致度和
遗传距离的无偏估计值
犜犪犫犾犲7 犖犲犻’狊狌狀犫犻犪狊犲犱犿犲犪狊狌狉犲狊狅犳犵犲狀犲狋犻犮犻犱犲狀狋犻狋狔犪狀犱
犵犲狀犲狋犻犮犱犻狊狋犪狀犮犲犫犲狋狑犲犲狀犳狅狌狉犵犲狅犵狉犪狆犺犻犮犵狉狅狌狆狊
地理类群
Groups
四川
Sichuan
云南
Yunnan
重庆
Chongqing
贵州
Guizhou
四川Sichuan - 0.9656 0.8597 0.9584
云南 Yunnan 0.0350 - 0.8243 0.9268
重庆Chongqing 0.1512 0.1932 - 0.9088
贵州 Guizhou 0.0425 0.0760 0.0956 -
注:Nei氏遗传一致度(对角线上方)和遗传距离(对角线下方)。
Note:Nei’sgeneticidentity(abovediagonal)andgeneticdistance(be
lowdiagonal).
分子标记具有信息量大、效率高、不受环境限制和影响等特点,它一般覆盖整个基因组,可以很好的揭示供试群体
的遗传多样性。利用SSR标记技术检测了中国西南地区11份鸭茅种质的遗传多样性。SSR标记结果表明,21
对引物可在110份鸭茅基因组DNA的156个位点扩增出条带,多态性条带比率为90.7%,在分子水平上证明了
中国西南区鸭茅种质资源遗传多样性丰富。本试验中,鸭茅的SSR标记表现出高的多态性(90.7%),高于彭燕
等[11]用AFLP技术对37份鸭茅种质研究的多态性(84.0%),也高于曾兵等[12,14]利用ISSR和SRAP技术对鸭
茅的遗传多样性研究所报道的结果(2种标记多态性分别为86.3%和84.38%),由于鸭茅资源在地理上分布广
泛,不同地区土壤的酸碱度、积温和日照时数等生态环境因子可能存在明显的差异,这些生态因子对鸭茅形成长
期的自然选择,加上自然突变、人为选择等其他因素,使得鸭茅在形态特征、生态类型、生理适应性和基因等方面
保持了较大的变异度和丰富的多样性。由此可见,SSR技术可以作为研究鸭茅遗传多样性的简单、有效的手段。
3.2 分子变异分析(AMOVA)进一步研究了鸭茅种质内遗传变异的分布,表明63.93%的遗传变异存在于种质
内部,而种质间变异较小,贡献率仅为21.93%。该结果与Koliker等[13]用RAPD方法研究发现的鸭茅种质内
具有较高的遗传变异(85.1%)的结论吻合。同时本研究所得到的鸭茅种质内的遗传变异模式与那些多年生且具
有异交繁育系统的植物具有相似的遗传变异模式[25]。这可能是由于鸭茅是一种异交植物,自交很少发生,而异
交种的遗传变异常发生在居群内或种质内,居群或种质间遗传分化较少[26]。Lumaret和Barrientos[27]用同工酶
341第18卷第4期 草业学报2009年
图2 11份鸭茅种质材料基于犖犲犻’狊遗传一致度的犝犘犌犕犃聚类
犉犻犵.2 犝犘犌犕犃犱犲狀犱狉狅犵狉犪犿犳狅狉犇.犵犾狅犿犲狉犪狋犪犫犪狊犲犱狅狀犖犲犻犔犻’狊犵犲狀犲狋犻犮犻犱犲狀狋犻狋狔
图3 4个鸭茅地理类群基于犖犲犻氏遗传一致度的聚类
犉犻犵.3 犇犲狀犱狉狅犵狉犪犿狅犳犳狅狌狉犵犲狅犵狉犪狆犺犻犮犪犾犵狉狅狌狆狊犫犪狊犲犱狅狀犖犲犻’狊狌狀犫犻犪狊犲犱犿犲犪狊狌狉犲狊狅犳犵犲狀犲狋犻犮犻犱犲狀狋犻狋狔
标记方法研究了形态上无明显差异的二倍体和四倍体鸭茅。结果显示,二倍体和四倍体的基因之间有着极高的
相似性,并认为自然界的鸭茅可能起源于各种比率杂交二倍体的同源化。本研究发现二倍体与四倍体群体间的
遗传分化较小,仅为14.13%,表明二者亲缘关系较近,基因之间有较高的相似性。Xu等[28]研究表明多倍体植物
存在更高的遗传变异和多样性,本研究发现四倍体群体比二倍体群体遗传多样性更为丰富,印证了这一结论。本
研究还发现四倍体种质遗传多样性相似性高(犘犘犅在68.59%~71.79%),但二倍体种质间存在较大遗传差异
(犘犘犅在10.25%~67.95%),与帅素容等[29]利用同工酶发现二倍体、四倍体种质内同源性较高,遗传差异较小
的结论有一定的差异。其原因可能有3个方面:一是试验材料和研究手段不同,本研究材料来源地分布广泛,其
地理和气候特征存在明显差异;另外同工酶是基因的产物,利用同工酶进行多样性研究时可能会有诸多不稳定因
素存在。比如同一材料不同器官酶带数和酶活性在同工酶类型上存在较大的变化。在不同的生长发育期,其过
氧化物酶谱表现出相当大的差异等,SSR分子标记则是针对DNA所作的分析,DNA的稳定性保证了SSR分析
结果的可靠性和稳定性。二是鸭茅是异花授粉植物,不同种质间花粉的传播是造成变异的重要因素之一,供试的
3份二倍体材料YA20064、YA20063和YA02101前期生长缓慢,生育期晚于其余材料,这可能影响了花粉的
传播,减少了种质杂交所引起的基因混杂,从而其独特的基因资源得以保留。由于内在的遗传变异最终会表现为
441 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4
材料栽培学形态和数量上的差异,所以更实用的分析应将分子标记显示的遗传差异与材料的形态学和生化特征
结合进行。
3.3 聚类分析表明,供试野生鸭茅材料的聚类与倍型特征有较强的相关性,且呈现出一定地域性分布规律。倍
型相同的材料优先聚为一类,然后相同地理分布及地形特征相近的材料聚在一起。如几份二倍体种质基本能优
先聚在一起,来自四川宝兴和重庆材料分别聚为一类。但这也并非绝对,如来自贵州的材料分别聚为I类和II
类,二倍体与四倍体在聚类时也不能完全分开,如四川汉源的二倍体YA917与四川宝兴和云南昆明的3个四倍
体聚为一类。这种聚类不完全符合地域性分布规律和倍性水平的原因可能有3个方面:一是花粉串粉导致了天
然杂交,鸭茅是异花授粉植物,在材料收集与保存过程中,忽略了种质隔离,不同地理来源的材料间发生了天然杂
交;二是基因突变,物种在适应自然选择的过程中,有个别的基因发生了突变,且该突变能较好的适应环境,这个
突变就能保留下来;三是在资源采集、保存、整理及实验室工作中(包括种子的收集、DNA提取、保存、PCR扩增
等方面)出现样本人为的混淆,试验中如果保持严谨的工作态度,人为误差将能降低到最小程度。
3.4 自然资源保护的一个主要目标是保护基因多样性,认识和了解物种的遗传变异,对于合理开发利用物种资
源具有重要意义。本研究结果表明,中国西南区鸭茅资源遗传多样性丰富,其遗传变异主要存在于种质内和地理
区域内。这就要求在鸭茅资源收集和保护中要注意广泛性和特殊性,即所有的生态类型分布区域要进行保护,同
时在资源收集时增加同一区域内的样本采集数量;对于一些特殊的遗传资源要设立有效的自然保护区如重庆和
贵州,从而使一些特殊的有利遗传性状得到保护。鉴于鸭茅是异花授粉植物,在广泛收集种质资源的同时,应注
意隔离防止串粉,妥善保存种质资源以利于利用。
参考文献:
[1] 贾慎修.中国饲用植物志(第一卷)[M].北京:农业出版社,1985.6973.
[2] LindnerR,GarciaA.Geneticdifferencesbetweennaturalpopulationsofdiploidandtetraploid犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪spp.izcoi[J].Grass
andForageScience,1997,52(3):291297.
[3] VilharB,VidicT,JoganN,犲狋犪犾.Genomesizeandthenucleoarnumberasestimatorofploidylevelin犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪in
theSolvenianAlps[J].PlantSystematicsandEvaluation,2002,234:113.
[4] 彭燕,张新全.鸭茅种质资源多样性研究进展[J].植物遗传资源学报,2003,4(2):179183.
[5] 张新全,杜逸,郑德诚.鸭茅二倍体和四倍体PMC减数分裂,花粉育性及结实性的研究[J].中国草地,1996,(6):3840.
[6] 张新全.几种禾草染色体核型和开花习性的研究[D].雅安:四川农业大学,1992.
[7] Reeves,DennisFM,StuartD,犲狋犪犾.Genomesizeisnegativelycorrelatedwithaltitudeinnaturalpopulationsof犇犪犮狋狔犾犻狊
犵犾狅犿犲狉犪狋犪[J].AnnalsofBotany,1998,82(supplementA):99105.
[8] 彭燕,张新全,曾兵.野生鸭茅植物学形态特征变异研究[J].草业学报,2007,16(2):6975.
[9] SugiyamaS.Geographicaldistributionandphenotypicdifferentiationinpopulationof犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪L.inJapan[J].Plant
Ecology,2003,169(2):295305.
[10] TunaM,SabancCO.Aevaluationofsomeagronomiccharactersofnaturalorchardgrass(犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪L.)popula
tionscolectedfromThraceregionofTurkey[J].BulgarianJournalofAgriculturalScience,2004,10(1):5764.
[11] 彭燕,张新全,刘金平,等.野生鸭茅种质遗传多样性的AFLP分子标记[J].遗传,2006,28(7):845850.
[12] 曾兵,张新全,范彦,等.鸭茅种质资源遗传多样性的ISSR研究[J].遗传,2006,28(9):10931100.
[13] KolikerR,StadelmannFJ,ReidyB,犲狋犪犾.Geneticvariabilityofforagegrasscultivars:Acomparisonof犉犲狋狌犮犪狆狉犪犲狀狊犻狊
Huds.,犔狅犾犻狌犿狆犲狉犲狀狀犲L.,and犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪L.[J].Euphytica,1999,106:261270.
[14] ZengB,ZhangXQ,LanY,犲狋犪犾.Evaluationofgeneticdiversityandrelationshipsinorchardgrass(犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪L.)
germplasmbasedonSRAPmarkers[J].CanadianJournalofPlantScience,2008,88:5360.
[15] SaghaiMaroofMA,SolimanKM,JorgensenRA.RibosomalDNAspacerlengthpolymorphismsinbarely:Ribosomalin
heritance,chromosomallocation,andpopulationdynamics[J].ProceedingNationalAcademyofSciencesUSA,1984,81:801.
[16] 蒋昌顺,马欣荣,邹冬梅,等.应用微卫星标记分析柱花草的遗传多样性[J].高技术通讯,2004,14(4):2530.
[17] 李永祥,李斯深,李立会,等.披碱草属12个物种遗传多样性的ISSR和SSR比较分析[J].中国农业科学,2005,38(8):
541第18卷第4期 草业学报2009年
15221527.
[18] 张阿英,胡宝忠,姜述君,等.影响紫花苜蓿SSR分析因素的研究[J].黑龙江农业科学,2002,(1):1517.
[19] 郭海林,刘建秀,高鹤,等.结缕草属优良品系SSR指纹图谱的构建[J].草业学报,2007,16(2):5359.
[20] 隋晓青,王.克氏针茅ISSR-PCR反应体系的建立与优化[J].草业学报,2008,17(3):7178.
[21] RoldanRuizI,DendauwJ,VanBockstaeleE,犲狋犪犾.AFLPmarkersrevealhighpolymorphicratesinryegrasses(犔狅犾犻狌犿
spp.)[J].MolecularBreeding,2000,6:125134.
[22] 车永和,李立会,何蓓如.冰草属植物遗传多样性取样策略基于醇溶蛋白的研究[J].植物遗传资源学报,2004,(5):216
221.
[23] ExcoffierL,SmousePE,QuattroJM.AnalysisofmolecularvarianceinferredfrommetricdistanceamongDNAhaplotypes:
ApplicationtohumanmitochondriaDNArestrictionsites[J].Genetics,1992,131:479491.
[24] RohlfFJ.NTSYSpc:NumerialTaxonomyandMultivariateAnalysisSystem (Version2.0)[M].NewYork:ExterSoft
were,1994.
[25] NybomH.ComparisonofdifferentnuclearDNAmarkersforestimatingintraspecificgeneticdiversityinplants[J].Molecular
Ecology,2004,13:11431155.
[26] JonesK,CarrolCP,Borril MA.Chromosomeatlasofthegenus犇犪犮狋狔犾犻狊L.[J].Cytologia,1961,26:333343.
[27] LumaretR,BarrientosE.Phylogeneticrelationshipsandgeneflowbetweensympatricdiploidandtetraploidplantsof犇犪犮狋狔犾
犻狊犵犾狅犿犲狉犪狋犪(Gramineae)[J].PlantSystematicsandEvolution,1990,169:12,8196.
[28] XuW W,SleperDA,KrauseGF.GeneticdiversityoftalfescuegermplasmbasedonRFLPs[J].CropScience,1991,34:
246252.
[29] 帅素容,张新全,白史且.不同倍性鸭茅同工酶比较研究[J].草业科学,1998,6:1116.
犌犲狀犲狋犻犮狏犪狉犻犪狋犻狅狀狅犳犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪犵犲狉犿狆犾犪狊犿犳狉狅犿犛狅狌狋犺狑犲狊狋犆犺犻狀犪犱犲狋犲犮狋犲犱犫狔犛犛犚犿犪狉犽犲狉狊
XIEWengang1,ZHANGXinquan1,MAXiao1,HUANGLinkai1,ZENGBing2
(1.DepartmentofGrasslandScience,SichuanAgriculturalUniversity,Ya’an625014,China;
2.DepartmentofAnimalScience,SouthwestUniversity,Chongqing402460,China)
犃犫狊狋狉犪犮狋:Simplesequencerepeat(SSR)molecularmarkerswereappliedtodetectthegeneticvariationof110
individualsof11犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪accessionsfromSouthwestChina(Sichuanprovince,Yunnanprovince,
GuizhouprovinceandChongqingcity).1)21primerpairswereselectedforSSRanalysisand143polymorphic
bandswerefoundwithanaverageof6.8bandsperprimerpair.Thepercentageofpolymorphismbands(犘犘犅)
was90.7%,thepolymorphisminformationcontent(犘犐犆)rangedfrom0.17(A01F24)to0.45(A01E02)with
anaverageof0.33,andtheShannon’sinformationindexofdiversity(犐)rangedfrom0.0463to0.3477,
showingconsiderablegeneticdiversitybetweenaltheaccessions.2)TheAMOVAanalysisshowedconsidera
blegeneticvariationwithinaccessions(63.39%)ratherthanbetweenthem (21.93%).Therewasnosignifi
cantgeneticvariationbetweendiploidandtetraploidgroups.3)UPGMAclustersbasedonaverageNei’sgenet
icidentitydivided11accessionsintotwogroups:3diploidaccessions(YA20063,YA20064,YA02101)col
lectedfromChongqingandGuizhouwereclassifiedintogroupIandtheremaining8accessionswereclassified
intogroupII.Inmostcases,materialsfromthesameregionandploidylevelwereinthesamegroupthus
showingageographicaldistributionofgeneticdiversityofthetestedmaterials.
犓犲狔狑狅狉犱狊:犇犪犮狋狔犾犻狊犵犾狅犿犲狉犪狋犪;germplasmresources;SSR;geneticdiversity
641 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4