全 文 :书高羊茅单株耐热性相关分子标记的
筛选及其与越夏性的关系研究
陈群1,袁晓君2,何亚丽1
(1.上海交通大学农业与生物学院,上海200240;2.上海大学生命科学技术学院,上海200444)
摘要:筛选与高羊茅单株耐热性和越夏性相关的分子标记有利于建立高效育种体系。以来自于耐热性不同的12
个高羊茅种质资源的216个单株为耐热性分子标记筛选的基础群体,在人工高温(40℃/35℃,昼/夜各12h)和夏
季田间自然条件下分别进行了耐热性和越夏性鉴定。运用简单重复序列区间(intersimplesequencerepeat,ISSR)
和随机扩增多态性DNA(randomamplifiedpolymorphismDNA,RAPD)分子标记技术,采用改良的集团分离分析
法(modifiedbulkedsegregantanalysis,BSA),从800个RAPD和100个ISSR分子标记中筛选与耐热性相关分子
标记,并研究其与单株越夏性的关系。结果表明,单株的耐热性与越夏性的相关程度较低(决定系数为3.4%);引
物351扩增的750bp条带(351T750)与耐热性显著相关,而引物138扩增的950bp条带(138T950)则与耐热性
和越夏性均显著相关。可以推论,耐热性强和/或越夏性强单株材料需要分别在人工高温和夏季田间自然条件下
加以鉴定筛选;351T750可应用于耐热性强单株的辅助选择,而138T950则可以用于耐热性强与越夏性强单株的
辅助选择。
关键词:高羊茅;耐热性;越夏性;简单重复序列区间分子标记(ISSR);随机扩增多态性DNA(RAPD);集团分离分
析法(BSA)
中图分类号:S816;S540.1 文献标识码:A 文章编号:10045759(2013)05008412
犇犗犐:10.11686/cyxb20130510
高羊茅(犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪)又称苇状羊茅,多年生冷季型和丛生型禾草[1],是温带地区广泛应用的草坪草
和牧草[2]。高羊茅为异花授粉植物,由3个染色体组组成的异源六倍体(PPG1G1G2G2,2狀=6狓=42),基因组包
含约50亿个碱基对,大小为(5.27~5.83)×106kb[1]。在多个冷季型草种中,虽然高羊茅的耐热性最强,周年绿
色期最长,但其适宜的地上部生长温度为15~24℃,根系为10~18℃[3],在我国温带和亚热带的夏季也常会遭受
高温胁迫,使生长迟缓或停滞,抗病虫害能力下降,从而造成牧草产量减产,草坪质量下降。因此,改良高羊茅耐
热性,培育耐热新品种,成为高羊茅育种中的当务之急。
在常规育种工作中,鉴于高羊茅是以种子繁殖为主的异花授粉植物,一般通过对不同种源材料杂交所得的后
代进行连续几代的随机授粉和选择达到群体改良的目的,因而其品种多为异质杂合群体[4]。品种群体内存在的
大量变异严重阻碍了基于表型性状的品种鉴定方法的应用。传统品种、品系(含杂交亲本、自交系及供无性繁殖
用的营养体)的鉴定方法是基于对农艺性状表型的识别,须耗费大量的时间和人力,而且常常会受到环境因素影
响,从而造成对多基因性状选择的干扰[4,5]。而耐热性作为一个耐逆性状,往往受到多基因的控制,其表现出的
复杂性使得对其进行形态选择变得非常困难[6]。分子标记辅助选择(molecularmarkerassistedselection,MAS)
被认为是目前改良植物抗逆性的一个有效方法,这一方法应用的关键在于首先需要找到与耐热性相关基因或数
量性状基因座(QTLs)紧密连锁且能稳定检测的理想分子标记。这些分子标记的获得不但有利于推动高羊茅耐
热性性状的育种进程,同时也将为揭示耐热性性状的分子遗传规律提供信息。近年来随着禾草转基因和表达研
究的不断深入[79],开发出更多的和耐逆性状相关的分子标记也有助于相关基因的定位和克隆。
84-95
2013年10月
草 业 学 报
ACTAPRATACULTURAESINICA
第22卷 第5期
Vol.22,No.5
收稿日期:20121203;改回日期:20130221
基金项目:国家自然科学基金面上项目(30972134)资助。
作者简介:陈群(1986),男,江西鹰潭人,在读硕士。Email:cq21914@sjtu.edu.cn
通讯作者。Email:heyali@sjtu.edu.cn
在众多的分子标记中,随机扩增多态性DNA(randomamplifiedpolymorphismDNA,RAPD)分子标记技术
操作简单、快速、省时、DNA用量少、设计引物无需预先知道序列信息,但是实验的稳定性和重复性较差。简单重
复序列区间(intersimplesequencerepeat,ISSR)分子标记技术实验操作简单、快速、高效、模板用量少,稳定性和
多态性高,标记为显性标记,符合孟德尔遗传规律,广泛地运用到分子标记工作当中[1014]。一般来说,在筛选与目
的基因紧密连锁的分子标记研究中,首先要通过多年回交和筛选建立近等基因系,工作量非常大[15]。Michel
more等[16]发展的集团分离分析法(bulkedsegregantanalysis,BSA)首先从一对具有目标基因的表型差异的亲
本所产生的任何一种分离群体中,根据目标基因的表型分别选取一定数量的植株,构成2个亚群或集团,将每群
的DNA等量混合,形成2个相对性状的“基因池”(genepool),然后用合适的分子标记对2个基因池进行分析,
在两群间表现多态性的分子标记遗传上与目标性状基因座位相连锁。由于BSA法构池过程中2个基因池中所
要标记的差异突出出来,而其他基因或位点的分布基本相同,使得2个DNA池总体表现出来的遗传背景近似相
等,这样大大提高了检测出与目的基因紧密连锁的分子标记的可能性[15]。
本研究在人工高温和田间夏季自然条件下,对来自于12个耐热性不同的种质资源的216个高羊茅单株进行
耐热性和越夏性鉴定,以极端耐热和不耐热单株分别建立耐热与不耐热基因池(拟BSA法),在800个RAPD和
100个ISSR引物中筛选与耐热性相关的分子标记,并研究其与越夏性的相关性,为提高高羊茅的耐热性和越夏
性育种提供有效的辅助选择工具。
1 材料与方法
1.1 供试材料
本试验的植物材料是由上海交通大学草坪草抗逆育种课题组经多年多点试验鉴定为耐热性强、中、差的3个
商业品种[17]和9个间接引自中国草种质资源中期库未经人工选择干预和未经耐热性鉴定的野外种质资源,其编
号与来源见表1。
1.2 植株准备
试验于2010年10月初,在上海交通大学农业与生物学院人工气候室内,将上述12个材料的种子播于盛满
园艺介质(上海孙桥农业技术有限公司生产,草炭∶珍珠岩=7∶3=v∶v)的10cm直径底部带出水口的育苗盆
中,浇水保持湿润使种子发芽与幼苗生长。播种21d后,从每个材料中随机抽取18株长势长相相似的种子萌发
苗,移栽至3个6×12穴的秧盘(又称“穴盘”)中,即每盘12行,每行6穴,每穴孔栽种1株。对于不同种质资源
材料的耐热性比较采用的是随机区组试验设计,重复3次,每材料在每重复(区组,穴盘)中6个单株。每穴容积
为2cm×2cm×3cm。栽了植株的秧盘置于底部平铺基质的托盘内,以保证秧盘平放和水分均衡供应。人工气
候室内温度设置为25℃/20℃(昼/夜各12h),光照强度为540μmol/(m
2·s)、相对湿度70%。每周每穴施1次
50mLHoagland营养液,隔天浇水1次,每7d修剪1次。育苗期间喷施甲基托布津/代森锰锌混合药液防治褐
斑病2次。热胁迫开始前2d将植株修剪到5cm高。每天目测检查介质中湿度,保持充足但不过量的水分供应。
1.3 热胁迫处理与单株耐热性排序和评分
2011年3月29日将人工气候室温度调至40℃/35℃(12h/12h),光强540μmol/(m
2·s),相对湿度70%,
对材料进行热胁迫。试验期间每天浇水1次,严格控制水量,使水分均衡供应。每周每穴施1次 Hoagland营养
液,50mL。为了降低不同部位受到的光照强度与温度的不同影响带来的试验误差,每天横向平移交换秧盘所在
位置,3d为1个平移周期,然后将秧盘反转180°后进行同样的平移处理。热胁迫期间每天考察植株在高温胁迫
下的表现,以目测法判定植株耐受热胁迫的天数,即自热胁迫当天开始至植株表现为叶色发黄和失水萎蔫严重的
天数 (下简称“耐受天数”)。在植株出现这种极度热伤害症状后,从热胁迫的人工气候室中取出,移栽到新的装
有园艺介质的秧盘中,置于常温下恢复生长。
同批取出的植株耐受热胁迫的天数相同,但不同单株间仍然能够目测区分出伤害程度不同,依热胁迫伤害最
重至相对最轻排序。例如热胁迫28d有首批11株出现萎蔫症状,其受害最重的为耐热性最差的植株,不耐热性
排名第1,其余依次排名至第11。第2批取出的植株中最不耐热的排第12,依次类推,总共216株植株的不耐热
性从1排至216位,第216位植株为最耐热的植株。
58第22卷第5期 草业学报2013年
表1 供试材料来源及特性
犜犪犫犾犲1 犜犺犲狊狅狌狉犮犲狊犪狀犱犮犺犪狉犪犮狋犲狉犻狊狋犻犮狊狅犳狋犺犲狋犲狊狋犲犱犿犪狋犲狉犻犪犾狊
材料编号
Materialsentry
产地或商品名称
Placeoforiginorcommercialname
类型和特点
Typeandmaincharacteristics
来源
Distributor
TF01 日本Japan 国外野生资源 Wildresourcesabroad 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF02 阿根廷Argentina 国外野生资源 Wildresourcesabroad 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF03 宁夏回族自治区 NingxiaAutonomous
RegionoftheHuiNationality
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF04 新疆伊犁昭苏 YiliZhaosuofXinjiang
AutonomousRegion
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF05 新疆博勒赛里木湖BolerSailimuLake
ofXinjiangAutonomousRegion
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF06 新疆谢家沟 XiejiagouofXinjiangAu
tonomousRegion
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF07 甘肃黄羊镇 HuangyangTownofGan
suProvince
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF08 新疆塔城TachengofXinjiangAutono
mousRegion
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF09 新疆乌市南戈壁SouthernGebiofUru
mqiCityofXinjiangAutonomousRegion
国内野生资源Domesticwildresources 国家草种质资源中期库StateMedium
termRepositoryforGrassGermplasms
TF10 家园CV.Plantation 耐热和越夏性中等品种Thecultivarsmedi
umtoleranttoheatandsummerstress
中国种子集团公司ChinaNationalSeed
GroupCO.,LTD.
TF11 强劲CV.Inferno 耐热和越夏强品种 Thecultivarsstrongly
toleranttoheatandsummerstress
北京 布 莱 特 草 业 有 限 公 司 Beijing
BrightGrassIndustryLTD.
TF12 辉煌BrandHuihuang 商品组合Commercialblend 北京 布 莱 特 草 业 有 限 公 司 Beijing
BrightGrassIndustryLTD.
采用改进的9分制评分方法对216株的耐热性进行评分。把首批取出的耐受天数为28d的11株(犖=11)
中最不耐热植株的评分定为1分,最后一批取出的耐受天数为52d的38株植株中最不耐热的植株评为7分。
由于首尾2批植株的耐受天数相隔24d,且2期最差植株的评分相差6分,则平均每延长4d取出的批次的最差
植株加1分[(52-28)/(7-1)=4],而同一批次内各植株的评分首尾相差0.25分,这样,不耐热性排序每滞后1
位加上0.25/(狀-1)分。根据上述评分方法,总结出单株的评分(K)公式为:
犓=1+犱-284 +0.25×
狀-1
犖-1
其中,犱为热胁迫耐受天数,狀为同一耐受天数中植株的不耐热性排序位次,犖 为同一耐受天数的植株数。例如,
第1批移出人工气候室的植株,为热胁迫耐受天数28d(犱=28)的单株,小计有11株(犖=11),其中最不耐热的
植株评分为1,排序第2的不耐热植株的评分则为1.025。又例如,第3批取出的,耐受天数为32d的植株有11
株,其排序第5的植株的评分则为2.100分。以此类推,各耐受天数植株数与其评分标准列于表2中。
1.4 单株越夏性鉴定
把热胁迫处理后植株置于室内常温下恢复生长1个月,于2011年秋季把216株中能够恢复生长的193株移
栽到大田。株行距35cm×45cm,每月修剪1次,剪草高度12cm。适时施肥和进行病、虫、草害管理。经过
68 ACTAPRATACULTURAESINICA(2013) Vol.22,No.5
2012年夏季自然高温和干旱胁迫后,于2012年8月22日按改进的9分制对田间单株进行评分,进行越夏性鉴
定。夏季高温干旱胁迫后枯死的植株评为1分,尚有1~2个成活茎蘖的单株给2分,6分为可接受的植株,8分
为分蘖多、长势依然旺盛的植株,9分为没有出现受害症状的植株。
表2 单株耐热性排序与评分方法
犜犪犫犾犲2 犚犪狀犽犻狀犵犪狀犱狉犪狋犻狀犵狊狅犳犺犲犪狋狋狅犾犲狉犪狀犮犲犳狅狉狊犻狀犵犾犲狆犾犪狀狋狊
植株停止热胁迫批次
Batchnumberof
plantsending
heatstress
热胁迫耐受天数
Heatstressdaysof
plantscould
survive(d)
批次内株数
Plantsnumber
withinthe
batch(N)
各批次内最不耐热植株评分
Ratinggradeforthemost
heatsensitiveplants
withinthebatch
批次内植株耐热性排序每滞后1位增加
的评分单位Increasingratingunitfor
plantswitheachonelaginranking
ofheattolerancewithinbatch
1 28 11 1.00 0.025
2 30 12 1.50 0.023
3 32 11 2.00 0.025
4 38 34 3.50 0.008
5 39 16 3.75 0.017
6 40 22 4.00 0.012
7 42 14 4.50 0.019
8 43 19 4.75 0.014
9 45 5 5.25 0.063
10 47 14 5.75 0.019
11 50 20 6.50 0.013
12 52 38 7.00 0.007
1.5 DNA提取和拟BSA法构建基因池
在耐热性鉴定后,分别取各单株叶片约0.2g,采用十六烷基三甲基溴化铵(CTAB)法[18]提取基因组DNA。
用0.8%的琼脂糖凝胶电泳检测,并结合紫外分光光度仪检测所提取DNA的完整性和浓度。再将所提DNA用
pH8.0的TrisEDTA缓冲液稀释成相同浓度(20ng/μL),-20℃保存备用。选取最耐热和最不耐热的各11个
单株的DNA等量混合,分别构建耐热基因池(HT,heattolerance)和不耐热基因池(HS,heatsensitive)。
1.6 PCR反应体系和产物的检测
采用的RAPD和ISSR引物,均根据加拿大大不列颠哥伦比亚大学(UniversityofBritishColumbia)公布的
编号和序列,由上海生工生物技术服务有限公司合成,其中RAPD引物编号1~800,ISSR引物编号801~900。
聚合酶链式反应(PCR)在96孔板上进行。采用10μL反应体系:5单位(U)/μLTaq聚合酶1μL,10×缓冲液
(Mg2+ Plus)1μL,10mmol/L引物0.5μL,2.5mmol/LdNTP混合液0.2μL,20ng/μLDNA模板1μL,50%
甘油2μL,双蒸水4.3μL。Taq酶、10×缓冲液、dNTP混合液均购自中国大连Takara公司。PCR扩增条件为
94℃下预变性3min;进入循环,94℃下变性10s,复性温度(T=Tm-2℃)下复性30s,72℃下延伸45s,共35
个循环;最后在72℃下延伸6min,结束后于4℃下保存。
扩增产物加入4μL1∶1混合的1×SYBR? GreenI(北京鼎国生物有限公司)和溴酚蓝(Takara公司)混合
溶液,静置10min,取6μL混合液,用2%的琼脂糖凝胶电泳分离。用自动凝胶成像系统(Tanon3500,上海天能
科技有限公司)进行紫外观察并拍照记录。上述所有反应均重复2次。
1.7 数据分析
数据用 MicrosoftOfficeExcel2007统计分析软件进行方差分析、平均数的多重比较(LSD法)、χ
2 独立性检
验和相关分析等处理。
78第22卷第5期 草业学报2013年
2 结果与分析
2.1 耐热性排序与评分和构池单株的确定
各批次取出(不同热胁迫耐受天数)的植株数目见
表2。最不耐热的植株可以耐受28d热胁迫,最耐热
的单株可以耐受52d及以上的热胁迫强度,而一半植
株耐受热胁迫的天数为42d。结合耐受天数和热伤害
的形态,按最耐热到最不耐热的顺序编号为1~216,1
为最不耐热单株,216为最耐热单株,从而得到216个
单株的耐热性排序。在恢复生长的过程中,有23个单
株死亡。在幸存的193株单株中,选取最耐热和最不
耐热的单株各11株构成不耐热基因池(HS)和耐热基
因池(HT)的植株(表3)。
对单株的耐热性评分结果(图1)表明最早受到热
伤害,且伤害程度最大的植株(即最不耐热植株)的评
分最低,为1分。形态上出现热伤害症状最晚,且伤害
程度最轻的植株(即最耐热植株)的得分最高,为7.25
分。全部单株的平均得分为4.63分,标准差1.81,峰
度系数为-0.83,偏度系数为-0.196,为略微左偏的
趋近正态分布(图1)。193株中,评分大于等于7分的
耐热单株有37株,占19%,而评分小于等于3分的不
耐热单株为21株,占10.9%。
表3 不同材料的耐热性与构池单株来源
犜犪犫犾犲3 犎犲犪狋狋狅犾犲狉犪狀犮犲狅犳狋犺犲犿犪狋犲狉犻犪犾狊犪狀犱狋犺犲狅狉犻犵犻狀
狅犳狊犻狀犵犾犲狆犾犪狀狋狊犮狅犿狆狅狊犲犱狋狅犇犖犃狆狅狅犾狊
材料编号
Materials
entry
质量评分
Visual
quality
耐热基因池植株数
Plantscomposed
toheattolerance
DNApool
不耐热基因池植株数
Plantscomposed
toheatsensitive
DNApool
TF09 5.87a 5
TF11 5.44ab 1 1
TF08 5.25abc 1
TF07 5.18abc 3
TF01 5.07abcd 2
TF12 4.77abcd 1 2
TF05 4.54bcde
TF02 4.34bcde
TF10 4.17cde 2
TF04 4.02de
TF03 3.51e 1
TF06 3.46e 3
不同字母表示显著差异(犘<0.05),LSD0.05=1.117。
Differentlettersindicatesignificantdifferenceat犘<0.05,LSD0.05=
1.117.
对不同材料的耐热性评分进行方差分析结果表
图1 单株耐热性评分分布
犉犻犵.1 犇犻狊狋狉犻犫狌狋犻狅狀狅犳狉犪狋犻狀犵犵狉犪犱犲狅犳
犺犲犪狋狋狅犾犲狉犪狀犮犲狅犳狊犻狀犵犾犲狆犾犪狀狋狊
明,不同材料间耐热性差异显著。对不同材料耐热性
评分的平均数多重比较结果(表3)表明,TF09(新疆
乌市南戈壁)的耐热性最强,单株平均分为5.87,显著
高于评分为4.54及以下的TF05等材料。与TF09的
耐热性没有显著差异的材料有TF11,TF08等5个,
均标有字母a。最不耐热的材料为 TF06(新疆谢家
沟),与其没有显著差异的有TF03(宁夏),TF04等5
个材料,它们的平均数均标记有字母e。
耐热基因池植株大多来自于标有字母a的较耐热
种质资源,其中5/11来自于最耐热的TF09。不耐热
基因池植株主要来源于标有字母e的材料。耐热的材
料中选到最耐热的单株较多,不耐热的材料中选到不
耐热单株较多,耐热中间型的材料中则能选到2种类型的植株。
2.2 越夏性评分
单株越夏性评分最低为1分(枯死植株),占42.4%。把越夏性评分小于等于2的植株(有1~3蘖存活)作
为极端越夏性差的植株有130株,占67.4%。越夏性评分大于等于7的植株有19株,占5.7%。平均2.63,标准
差2.15,峰度0.56,偏度1.34(图2)。
对不同材料的越夏性评分进行方差分析结果表明,不同材料间差异显著。对不同材料越夏性评分平均数的
多重比较结果(图3)表明,越夏性最强的材料为TF07(甘肃黄羊镇),与之没有显著差异的材料有TF09,TF11等
4个材料,均标有字母a。越夏性最弱的为TF02(阿根廷),与之没有显著差异的有TF04,TF06,TF12,它们均标
88 ACTAPRATACULTURAESINICA(2013) Vol.22,No.5
记有字母c。TF01等4个标记有bc字母的材料越夏
图2 单株越夏性评分分布
犉犻犵.2 犇犻狊狋狉犻犫狌狋犻狅狀狅犳狊狌犿犿犲狉狋狅犾犲狉犪狀犮犲
狉犪狋犻狀犵犵狉犪犱犲狊犳狅狉狊犻狀犵犾犲狆犾犪狀狋狊
性为中等。
图3 不同材料的越夏性评分的比较
犉犻犵.3 犕狌犾狋犻狆犾犲犮狅犿狆犪狉犻狊狅狀狊狅犳狊狌犿犿犲狉狋狅犾犲狉犪狀犮犲
狉犪狋犻狀犵犵狉犪犱犲犫犲狋狑犲犲狀犲犪犮犺犿犪狋犲狉犻犪犾狊
不同字母表示显著差异(犘<0.05),LSD0.05=1.435。Differentletters
indicatesignificantlydifferenceat犘<0.05,LSD0.05=1.435.
2.3 耐热性和越夏性的相关性
分别对各材料内单株的耐热性与越夏性的相关分
析结果(表4)表明,只有TF01和TF06的耐热性和越
夏性之间的相关系数显著不等于零(犘<0.05),相关
程度(决定系数狉2)分别为32%和29%。对各材料耐
热性评分平均数与越夏性评分平均数的相关分析结果
表明相关性显著(犘<0.05),相关程度达43%。虽然
对所有193株的耐热性与越夏性的相关性分析结果也
表明相关显著,但相关程度极低,决定系数仅为3.4%。
2.4 池间多态性分子标记
用800个 RAPD引物对2个基因池进行多态性
筛选的结果表明,能扩增出条带的引物有632个,占总
引物数的79%,平均每个引物能扩增出2.7条清晰的
条带,但最终获得池间具有多态性的引物仅有6个,为
54,68,78,138,351和511引物(图4)。用100个
ISSR对2个基因池进行多态性的筛选,能扩增出清晰
条带的引物有62个,占引物总数的62%,平均能扩增
出5.2条清晰的条带,最终获得池间具有多态性的引
物有2个,为843和856引物(图4)。这8个引物的碱
基序列及退火温度与其在耐热和不耐热池间的特异性
条带见表5。
2.5 池间多态性分子标记在构池单株中的验证
上述筛选出的多态性条带根据引物序号及其片段
大小,分别命名为54S550、68S1000、78T800、78T900、
138T950、351T750、351T500、511S700、511S900、
511S1100、843S650,存在亮度差异的2个条带命名为856S950、856S300。对从耐热、不耐热DNA池中筛选出
来的6个RAPD和2个ISSR引物,用建成耐热、不耐热DNA池的各11个单株进行验证。结果在池间扩增出的
多态性条带和构池单株扩增结果相符的引物只有138(图5)和351(图6)。将耐热性强弱和目标条带的有无进行
卡方独立性测验,结果表明,351T750和耐热性有显著相关(χ
2=4.91,犘<0.05),而138T950标记和耐热性虽
没有显著相关(χ
2=2.54,犘>0.05),但在不耐热基因池的11个单株中没有该条带的吻合度也达到82%。
2.6 池间多态性分子标记在基础群体单株中的验证
把在池间与构池单株间多态性表现一致的2个标记“138T950”和“351T750”在193个基础群体单株中进行
扩增。结果(表6)表明,在37株耐热性强(评分≥7)单株中,138T950标记扩增出条带的有17株单株,占
45.9%,而在不耐热的21株(评分≤2.25)中扩增出条带的植株数为3株,仅占14.3%,卡方独立性检验结果显
示该标记与耐热性显著相关(χ
2=4.63,犘<0.05)。同样,351T750标记与耐热性也显著相关(χ
2=11.79,犘<
0.05),在耐热单株中扩增出条带的单株占59.5%,显著高于在不耐热单株的9.5%。此外,138T950标记与群
体单株的越夏性具有显著相关(χ
2=4.94,犘<0.05),在越夏性强单株(评分≥7)中扩增出条带的单株比例为
57.9%,显著高于越夏性弱(评分≤2)的29.2%,但351T750标记与群体单株的越夏性未表现出显著相关性
(χ
2=0.02,犘>0.05)(表6)。
98第22卷第5期 草业学报2013年
表4 耐热性与越夏性的相关系数和决定系数
犜犪犫犾犲4 犆狅狉狉犲犾犪狋犻狅狀犪狀犱犱犲狋犲狉犿犻狀犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋狊犫犲狋狑犲犲狀犺犲犪狋犪狀犱狊狌犿犿犲狉狋狅犾犲狉犪狀犮犲
变异来源
Sourceofvariation
简单相关系数
Simplecorrelationcoefficient
决定系数
Determinationcoefficient(狉2)
自由度
Degreeoffreedom
最小显著相关系数
狉0.05
概率
犘
TF01 0.57 0.32 15 0.482 <0.05
TF02 0.31 0.10 15 0.482 >0.05
TF03 0.21 0.05 13 0.514 >0.05
TF04 0.01 0.00 14 0.497 >0.05
TF05 0.20 0.04 13 0.514 >0.05
TF06 0.54 0.29 13 0.514 <0.05
TF07 -0.14 0.02 14 0.497 >0.05
TF08 0.13 0.02 15 0.482 >0.05
TF09 0.02 0.00 13 0.514 >0.05
TF10 0.03 0.00 14 0.497 >0.05
TF11 0.11 0.01 16 0.468 >0.05
TF12 -0.34 0.12 14 0.497 >0.05
材料平均数 Meanofmaterials 0.66 0.43 10 0.576 <0.05
所有单株Alplants 0.18 0.03 191 0.138 <0.05
图4 由6个犚犃犘犇(54,68,78,138,351和511)和2个犐犛犛犚(843和856)引物在耐热植株犇犖犃组成的
耐热基因池(犜)和不耐热植株犇犖犃组成的不耐热基因池(犛)间犘犆犚扩增产物的多态性
犉犻犵.4 犘狅犾狔犿狅狉狆犺犻狊犿狊犳狉狅犿犘犆犚犪犿狆犾犻犮犪狋犻狅狀狑犻狋犺6犚犃犘犇(54,68,78,138,351犪狀犱511)犪狀犱2犐犛犛犚(843犪狀犱856)
狆狉犻犿犲狉狊犫犲狋狑犲犲狀犫狌犾犽犲犱犇犖犃狆狅狅犾狊犳狉狅犿犺犲犪狋狋狅犾犲狉犪狀狋犻狀犱犻狏犻犱狌犪犾狊(犜)犪狀犱犺犲犪狋狊犲狀狊犻狋犻狏犲犻狀犱犻狏犻犱狌犪犾狊(犛)
3 讨论和结论
3.1 单株耐热性鉴定的可靠性
改良植株耐热性的一个主要挑战是确定可靠的筛选方法和选择有效的耐热性评价指标来检测植株的耐热
性。一般而言,田间的环境变异大,不适合评价和选育耐热单株或群体。人工气候室中温度可控,理论上可以从
苗期鉴定到生殖生长期。在耐热性鉴定指标的选择上,人们往往更倾向于把热胁迫后植株的形态与多个生理指
标相结合的综合评价方法,比用单一指标的评价结果更可靠[19,20]。但这种多指标的评价方法往往需要进行破坏
性取样,需要较大的处理群体,而且耗时费力,一般只适用于对少数几个品种或处理的植株群体的耐热性评价。
本研究是要对216个单株进行耐热性评价,以便获得构建BSA法研究体系的耐热与不耐热基因池所需要的极端
09 ACTAPRATACULTURAESINICA(2013) Vol.22,No.5
表5 池间多态性引物与犘犆犚扩增条带
犜犪犫犾犲5 犅犪狀犱狊狅犳狆狅犾狔犿狅狉狆犺犻狊犿狊犪犿狆犾犻犳犻犲犱犫狔犘犆犚狑犻狋犺狆狉犻犿犲狉狊犫犲狋狑犲犲狀狋犺犲犇犖犃狆狅狅犾狊
引物编号
Primerentry
序列
Sequence
退火温度
Annealingtemperature(℃)
耐热基因池扩增条带
BandsinT(bp)
不耐热基因池扩增条带
BandsinS(bp)
54 5′GTCCCAGAGC3′ 41.00 - 550
68 5′GAGCTCGCGA5′ 41.00 - 1000
78 5′GAGCACTAGC5′ 36.90 800,900 -
138 5′GCTTCCCCTT5′ 36.90 950 -
351 5′CTCCCGGTGG5′ 45.10 750,500 -
511 5′GAATGGTGAG5′ 32.80 - 700,900,1100
843 5′CTCTCTCTCTCTCTCTRA5′ 53.88 - 650
856 5′ACACACACACACACACYA5′ 53.88 950(强Strong) 950(弱 Weak)
300(弱 Weak) 300(强Strong)
T,耐热基因池 Heattolerantgenepool;S,不耐热基因池 Heatsensitivegenepool.
图5 引物138在构池耐热(犎犜)和不耐热(犎犛)单株上的扩增电泳图谱
犉犻犵.5 犘狉狅犳犻犾犲狅犳犪犿狆犾犻犳犻犮犪狋犻狅狀狆狉狅犱狌犮狋狊犫狔狆狉犻犿犲狉138犻狀犺犲犪狋狋狅犾犲狉犪狀狋(犎犜)
犪狀犱犺犲犪狋狊犲狀狊犻狋犻狏犲(犎犛)狆犾犪狀狋狊犮狅犿狆狅狊犲犱狋狅犇犖犃狆狅狅犾狊
图6 随机多态性引物351在构池耐热(犎犜)和不耐热(犎犛)单株上的扩增电泳图谱
犉犻犵.6 犘狉狅犳犻犾犲狅犳犪犿狆犾犻犳犻犮犪狋犻狅狀狆狉狅犱狌犮狋狊犫狔狆狉犻犿犲狉351犻狀犺犲犪狋狋狅犾犲狉犪狀狋(犎犜)
犪狀犱犺犲犪狋狊犲狀狊犻狋犻狏犲(犎犛)狆犾犪狀狋狊犮狅犿狆狅狊犲犱狋狅犇犖犃狆狅狅犾狊
耐热和不耐热的单株。由于单株多和单株的茎叶量少,不能满足需要破坏性取样的多个生理指标的测定需求。
植物在热胁迫后最明显的伤害症状是叶片发黄和植株失水萎蔫。本研究根据各单株在热胁迫过程中的热伤害症
状出现的先后和不同程度确定各个植株耐受热胁迫的天数和进行9分制评分,得到了216株的耐热性排序,方法
简便,快速高效。但是,对于216个单株而言,试验中没有重复,试验误差无法估计,这种评价结果的可靠性值得
关注。为了提高单株耐热性比较结果的可靠性,唯有依赖于降低试验误差,尽量使各单株所处的温、光、水、肥等
的环境条件相似。主要采取了以下控制误差的措施:托盘内垫满园艺介质,使置于其上的秧盘内植株处于同样的
水平面上,避免穴盘的高低不平造成介质中水肥含量的差异;定量供应每棵单株的水和营养液;循环平移和180°
旋转穴盘,使秧盘在一个热胁迫周期内受到的温、光影响相似。这些控制试验误差的措施对提高试验结果的准确
性的作用可以从12个不同材料的耐热性比较结果得到反映。对于参试的12个种质资源材料而言,每材料具有
19第22卷第5期 草业学报2013年
18个单株,构成有重复的试验,可以进行统计分析,估计误差,其误差大小可以反映该试验的精确性。统计分析
结果表明,不同材料的耐热性差异显著,说明材料间的平均耐热性差异显著大于由材料内单株基因型和试验误差
引起耐热性差异,误差引起的耐热性差异相对较小。在12个参试材料中,“强劲(TF11)”在本次研究中耐热性位
于第2,越夏性位于第3,分别与位于第1的材料间没有显著差异,这与在田间的草坪群体的越夏性鉴定结果[17]
和温室中的耐热性鉴定结果(未发表)一致,均名列前茅。“家园”和“辉煌”的试验结果也与先前的研究结果[17]吻
合。本次试验3个材料单株鉴定结果与多年多点草坪鉴定结果趋势一致也表明了本研究中单株耐热性与越夏性
鉴定结果的可靠性。
表6 “138犜950”和“351犜750”在基础群体中耐热性和越夏性两极单株中的检测结果
犜犪犫犾犲6 犞犪犾犻犱犪狋犻狅狀狅犳“138犜950”和“351犜750”犿犪狉犽犲狉狊犻狀狆犾犪狀狋狊狅犳犺犲犪狋犪狀犱狊狌犿犿犲狉
狋狅犾犲狉犪狀狋犪狀犱狊犲狀狊犻狋犻狏犲犵狉狅狌狆狊犻狀狋犺犲犫犪狊犻犮狆狅狆狌犾犪狋犻狅狀
单株来源Plantorigin 耐热性 Heattolerance
单株(强,弱)
Plants(tolerant,
sensitive)
138T950
NT NS
351T750
NT NS
越夏性Summertolerance
单株(强,弱)
Plants(tolerant,
sensitive)
138T950
NT NS
351T750
NT NS
观察值Observation 37,21 17 3 22 2 19,130 11 38 7 46
具有标记植株百分率Percentage
ofplantswithmarkers(%)
45.9 14.3 59.5 9.5 57.9 29.2 36.8 35.4
χ2检验χ2test χ2=4.63
犘<0.05
χ2=11.79
犘<0.05
χ2=4.94
犘<0.05
χ2=0.02
犘>0.05
NT,NS表示检测到相关标记在耐热性强、弱和越夏性强、弱单株上出现的单株数。
NTandNSrepresentthenumberoftheplantswithbandofthemarkerfoundintheheattolerant(NT)orsensitive(NS)andsummertolerant(NT)
orsensitive(NS)groups.
从构建耐热与不耐热基因池的单株来源来看,耐热基因池单株主要来自耐热性强的材料,不耐热单株主要来
自于耐热性差的材料,而耐热性中等的材料中既可筛选出少数耐热性强的单株,也能筛选到耐热性差的单株。这
一现象,既反映了育种上 “优中选优”的单株选择的一般法则,也进一步反映出本研究单株耐热性鉴定结果的可
靠性。
3.2 耐热性与越夏性的关系
本试验中热胁迫温度处理为40℃/35℃(12h/12h),高羊茅的热胁迫耐受天数最短为28d,最长为52d及
以上,这远远超过了我国长江中下游地区夏季高温天数的总和。高羊茅在长江中下游地区会出现越夏困难,显然
不是单纯的高温胁迫所致。本研究中193个单株的耐热性评分与越夏性评分的相关性分析表明,因自由度大,显
著水平狉0.05降低,而测验出相关显著 (显著不为零),但并不能说明相关性好,描述两者相关程度的应为决定系
数。本研究中单株耐热性与越夏性的决定系数为3.4% (狉2),表明单株的越夏性仅有3.4%与耐热性有关,还有
96.6%的影响越夏性的因子与耐热性无关。单株耐热性与越夏性无关的实际表现是,耐热强的单株在田间的越
夏性不一定强,而耐热性弱的,在田间自然温度下,越夏性也不一定弱。分析其原因主要是越夏性与耐热性的内
涵不同,夏季胁迫是多种胁迫因素的共同作用,高羊茅越夏性是多种抗性组成的综合抗性,包括耐热性、抗病性、
耐涝性、耐强光照等。此外,越夏性的内涵随各年份夏季气候的不同而异。例如:某些干旱年份的夏季胁迫因子
就是干旱、暴晒和高温,该年植物的越夏性就是指对暴晒、干旱和高温,以及由此条件下带来的生物胁迫的抗性。
而对于夏季多雨的年份,越夏性包含的是植株对高温高湿和由此带来的病虫草害的生物逆境胁迫的抗性。所以
当把越夏性作为考察指标时应关注当年的气候条件,明确越夏性的具体内涵。由于单纯的耐热性鉴定结果不等
同于田间越夏性结果,耐热单株和品种的选择不一定能够提高新品种(系)的越夏性。对于以提高自然条件下综
合抗性为目标的高羊茅育种而言,提高高羊茅的越夏性比提高单纯的耐热性更有实际应用价值。而对于开展耐
29 ACTAPRATACULTURAESINICA(2013) Vol.22,No.5
热性机理研究的材料选择,应以人工热胁迫试验结果为依据,选择耐热性强与弱的单株为比较材料,而不应以田
间的越夏性表现为选择依据。
由参试材料的耐热性与越夏性鉴定结果可以看到,来自我国新疆乌市南戈壁的高羊茅野生资源(TF09)和来
自甘肃黄羊镇的野生资源(TF07)的耐热性和越夏性分别名列前茅,与多年多点试验中表现出色的国外引进品种
“强劲”没有显著差异,说明我国自然界也蕴藏着丰富的优质高羊茅种质资源,值得驯化和开发利用。
3.3 耐热性与越夏性分子标记辅助选择
分子标记辅助育种在我国作物育种上应用的研究相对较多,而在草坪草、牧草上的研究相对较少。对高羊茅
分子标记的研究主要体现在遗传图谱的构建、比较基因组作图、亲缘关系和系统分类、遗传多样性的检测和品种
指纹图谱的绘制及非整倍体和属间杂种鉴定等方面[21]。Xu等[22]利用限制性片段长度多态性(restrictionfrag
mentlengthpolymorphism,RFLP)遗传图谱上的数量性状基因座(quantitativetraitlocus,QTLs)和相关分子的
紧密连锁,通过 MAS技术可以提高高羊茅的重要经济性状。但迄今为止,尚未检索到有关高羊茅耐热性DNA
分子标记的研究报道。
本研究得到的与高羊茅单株耐热性和越夏性相关的分子标记分别为351T750和138T950,其在同一基础
群体内不同样本容量的耐热性两极单株组群中,条带的有无与耐热性的强弱的相关性不同。对耐热性处于两极
的22个构池单株进行扩增验证,结果仅有351T750与单株耐热性显著相关,而138T950条带虽在不耐热的11
株中有82%的吻合度,但与22株构池单株的耐热性没有显著相关。当增加验证单株数至58株(37株耐热+21
株不耐热)后,351T750和138T950标记均与耐热性显著相关。同样的扩增结果,若以耐热性排序的中位数为
分界点进行耐热性分组,即把排序前96位的单株归类为耐热性弱单株,排序后97的为耐热性强单株,相关性分
析结果又与22株基因池内单株的表现相似,以351T750与耐热性呈显著相关(χ
2=9.21,犘<0.05),而138
T950与耐热性分组无显著相关(χ
2=1.04,犘>0.05)。虽然138T950标记与单株的耐热性相关随验证群体的
大小不同而有变化,但在与田间193株的越夏性相关分析中表现出稳定的显著相关性。两分子标记与耐热性和
越夏性的相关性变化,可能的原因是高羊茅作为异缘六倍体异花授粉植物,比水稻(犗狉狔狕犪狊犪狋犻狏犪)等常规作物的
基因组更为复杂(约比水稻大14倍),分子标记在与基因组结合的过程中容易出现随机性的不稳定结合,在试验
过程中表现出池间和单株间的多态性不够稳定。前人在对小麦(犜狉犻狋犻犮狌犿犪犲狊狋犻狏狌犿)[15]和紫花苜蓿(犕犲犱犻犮犪犵狅
狊犪狋犻狏犪)[23]等其他多倍体植物的相关分子标记筛选结果也有类似的现象报道。鉴于该两标记在耐热性和越夏性
鉴定的两极组群单株中具有显著相关性,实际应用中可结合单株的热胁迫与夏季胁迫结果,在表型耐热和越夏强
单株中进一步进行该两分子标记的辅助选择,从而提高单株选择准确率,加快育种进程。
许占友等[24]认为,BSA法的成功运用受到2个因素的影响:1)分子标记与目的基因实际距离的远近影响检
测效率;2)构建每个池所选择的植株数量影响检测能力;索广力[15]指出在构建基因池的时候,还需要考虑所用材
料的基因组的大小,基因组太大如小麦,建池的株数不宜大,否则虚假多态性太高。与已报道的BSA法的研究群
体(F2)不同,本研究以来自于耐热性不同种质资源的单株为分子标记筛选的基础群体,通过人工耐热性鉴定,以
最不耐热和最耐热性的11个单株构建了耐热和不耐热基因池(拟BSA法),通过RAPD和ISSR筛选池间稳定
的多态性分子标记,也成功得到了351T750和138T950两个与极端耐热和不耐热单株群体显著相关的分子标
记,其中138T950与单株在田间的越夏性也存在显著相关,可用于高羊茅耐热和越夏强单株的辅助选择。研究
结果初步表明,这种改进的拟BSA法在性状相关分子标记的筛选上是可行的,有效的。为了提高标记的可重复
性和稳定性,进一步的工作是将特异性条带进行回收、测序,并根据测序结果开发特异性更强的测序扩增区段标
记(SCAR),并应用于单株辅助选择。
4 结论
单株或材料的耐热性与越夏性相关程度较低,耐热强单株越夏性不一定强,常绿草坪草品种的选育应以提高
越夏性为育种目标。138T950与耐热性和越夏性均显著相关,而351T750标记仅与耐热性显著相关,可分别用
于在热胁迫和夏季胁迫后抗性强单株中的进一步筛选,以提高单株选择效率。
39第22卷第5期 草业学报2013年
参考文献:
[1] SealA.DNAvariationin犉犲狊狋狌犮犪[J].Heredity,1983,50(3):225236.
[2] SleperD.Breedingtalfescue[J].PlantBreedingReviews,1985,(3):313342.
[3] 韩烈保.草坪管理学[M].北京:北京农业大学出版社,1994.
[4] HoC,McMaughS,WiltonA,犲狋犪犾.DNAamplificationvariationwithincultivarsofturftypecouchgrasses(犆狔狀狅犱狅狀spp.)[J].
PlantCelReports,1997,16(11):797801.
[5] BudakH,ShearmanR,ParmaksizI,犲狋犪犾.Comparativeanalysisofseededandvegetativebiotypebuffalograssesbasedon
phylogeneticrelationshipusingISSRs,SSRs,RAPDs,andSRAPs[J].TheoreticalandAppliedGenetics,2004,109(2):280
288.
[6] 刘刊,耿士均,商海燕,等.草坪草抗性研究进展[J].草业科学,2012,29(7):716.
[7] 霍松,陈慧,朱琼华,等.象草4CL基因片段的克隆及 RNAi表达载体构建[J].草业学报,2012,21(1):296301.
[8] 安惠惠,马晖玲,李坚,等.农杆菌介导的LyzGFP基因对匍匐翦股颖PennA1转化和表达的研究[J].草业学报,2012,
21(2):141148.
[9] 张妙青,张吉宇,刘志鹏,等.垂穗披碱草 MADSbox基因 WM8克隆及分析[J].草业学报,2012,21(4):141150.
[10] 桂枝,袁庆华,蔺海明.RAPD标记技术及其在牧草研究中的应用[J].草业科学,2001,18(2):5056.
[11] 胡盨,王永清,陶炼.三角紫叶酢浆草ISSR反应体系的建立与优化[J].草业学报,2011,20(5):142150.
[12] 解新明,卢小良.利用 RAPD标记分析狼尾草属牧草品种间的遗传关系[J].草业学报,2005,14(2):5256.
[13] 李惠英,娄燕宏,胡涛,等.中国高羊茅种质资源遗传多样性的 RAPD分析[J].草业学报,2010,19(6):208214.
[14] 田志宏,邱永福,严寒,等.用 RAPD标记分析高羊茅的遗传多样性[J].草业学报,2007,16(1):5863.
[15] 索广力.利用 RAPD—BSA技术筛选小麦耐盐突变位点的分子标记[J].植物学报(英文版),2001,43(6):598602.
[16] MichelmoreRW,ParanI,KesseliR.Identificationofmarkerslinkedtodiseaseresistancegenesbybulkedsegregantanaly
sis:arapidmethodtodetectmarkersinspecificgenomicregionsbyusingsegregatingpopulations[J].ProceedingsoftheNa
tionalAcademyofSciences,1991,88(21):98289832.
[17] 鲍真晶,闫军辉,何亚丽.高羊茅抗逆无性系的选择和其子一代在上海的草坪特性鉴定[J].上海交通大学学报(农业科学
版),2012,30(3):7281.
[18] ClarkMS.PlantMolecularBiologyALaboratoryManual[M].Berlin:SpringerVerlagPress,1997.
[19] HuangB,JiangY.Droughtandheatstressinjurytotwocoolseasonturfgrassesinrelationtoantioxidantmetabolismand
lipidperoxidation[J].CropScience,2001,41(2):436442.
[20] 王志勇.冷季型草坪草耐热性的研究进展[J].草业科学,2006,23(10):8491.
[21] 马啸,张新全,周永红,等.高羊茅的分子标记应用进展[J].草业学报,2006,15(2):18.
[22] XuW,SleperD,ChaoS.Genomemappingofpolyploidtalfescue(犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪Schreb.)withRFLPmarkers[J].
TheoreticalandAppliedGenetics,1995,91(6):947955.
[23] 孟芳,袁庆华,苏德荣,等.紫花苜蓿抗褐斑病ISSR分子标记研究[J].草地学报,2007,15(6):561565.
[24] 许占友,张爱民,张树臻.小麦恢复基因 RAPD分子标记的研究[A].第一届国际杂交小麦研讨会论文集:杂交小麦[C].
北京:中国农业大学出版社,1998.
49 ACTAPRATACULTURAESINICA(2013) Vol.22,No.5
犛犮狉犲犲狀犻狀犵犿狅犾犲犮狌犾犪狉犿犪狉犽犲狉狊犳狅狉犺犲犪狋狋狅犾犲狉犪狀犮犲犪狀犱犻狋狊狉犲犾犪狋犻狅狀
狋狅狊狌犿犿犲狉狋狅犾犲狉犪狀犮犲犻狀狋犪犾犳犲狊犮狌犲狊犻狀犵犾犲狆犾犪狀狋狊
CHENQun1,YUANXiaojun2,HEYali1
(1.SchoolofAgricultureandBiology,ShanghaiJiaoTongUniversity,Shanghai200240,China;
2.SchoolofLifeSciences,ShanghaiUniversity,Shanghai200444,China)
犃犫狊狋狉犪犮狋:Toestablishanefficientbreedingsystemintalfescue(犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪),molecularmarkersre
latedtoheatandsummertoleranceswerescreenedformarkerassistedselectionbreeding.Initialy,216single
plantsfrom12germplasmsdifferinginheattolerancewereusedasabasicpopulationtoidentifytheirheat
stress(40℃/35℃,day/night,12h/12h)toleranceinanartificialclimatechamberandtheirsummerstress
toleranceinexperimentalfieldsinShanghai.RandomamplifiedpolymorphismDNA(RAPD)andintersimple
sequencerepeat(ISSR)wereassessedwithmodifiedbulkedsegregantanalysis(BSA)toscreenformolecular
markersamongthe800RAPDand100ISSRprimers.Thecorrelationbetweenheatandsummertolerance
amongsingleplantswasquitelowwithonly3.4%ofthedeterminationcoefficient;Primer“351”amplifieda
750bpband(351T750)whichwassignificantlyrelatedtoheattolerance,while“138”amplifieda950bpband
relatedtobothheatandsummertolerance.Heatandsummertolerantsingleplantsshouldbescreenedunder
artificialheatandnaturalsummerstressesrespectively;Band351T750couldbeusedonlyforheattolerant
plantselection,while138T950couldbeusedforbothheatand/orsummertolerantplantselection.
犓犲狔狑狅狉犱狊:talfescue;heattolerance;summertolerance;intersimplesequencerepeat(ISSR);randomampli
fiedpolymorphismDNA(RAPD);modifiedbulkedsegregantanalysis(BSA
檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵檵
)
欢迎订阅2014年 《农产品质量与安全》
主管 中华人民共和国农业部 主办 中国农业科学院
支持单位 农业部农产品质量安全监管局
协办单位 农业部农产品质量安全中心 中国绿色食品发展中心
承办单位 中国农业科学院农业质量标准与检测技术研究所
主要栏目:本刊特稿、本刊专稿、政策法规、质量安全监管、无公害农产品、绿色食品、有机农产品、农产品地理
标志、农业标准化、检验检测、学科建设与发展、研究与探讨、安全生产技术、地方经验交流、海外博览、农业标准公
告、信息与动态等。
读者对象:与农产品质量安全、农业质量标准和检验检测有关的各级行政管理、科研教学、检验监测、技术推
广、生产企业等部门的相关人员。
本刊为双月刊,逢双月10日出版。大16开本,彩色四封,80页。全国各地邮局(所)均可订阅,也可直接到
本刊编辑部办理订阅手续。邮发代号:82-223。每册定价:10.00元,全年共60.00元。
通讯地址:北京市中关村南大街12号中国农科院质标所《农产品质量与安全》编辑部 邮政编码:100081
联系电话/传真:(010)82106521、82106522 Email:aqs@caas.net.cn
欢迎各界朋友赐教、赐稿、订阅和刊登广告
59第22卷第5期 草业学报2013年