全 文 ：书玉米持绿相关犙犜犔整合图谱构建及一致性
犙犜犔区域内候选基因发掘
方永丰１，２，３，李永生３，白江平１，２，３，慕平３，孟亚雄１，２，３，
张金林４，王汉宁１，２，３，尚勋武３
（１．甘肃省作物遗传改良与种质创新重点实验室，甘肃 兰州７３００７０；２．甘肃省干旱生境作物学重点实验室，甘肃 兰州７３００７０；
３．甘肃农业大学农学院，甘肃 兰州７３００７０；４．兰州大学草地农业科技学院，甘肃 兰州７３００２０）
摘要：以玉米高密度遗传连锁图谱ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ为参考图谱，收集来自不同实验中的１７３个玉米持绿相关
数量性状位点（ｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｕｓ，ＱＴＬ）信息，利用ＢｉｏＭｅｒｃａｔｏｒ２．１软件，构建出玉米持绿相关 ＱＴＬ整合图
谱；采用元分析技术，在１，４，５，９号染色体上发掘出５个持绿“一致性”ＱＴＬ区间。根据“一致性”ＱＴＬ区间两端标
记在玉米物理图谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２上的位置，将“一致性”ＱＴＬ区间进行物理图谱定位，利用ＰｌａｎｔＧＤＢ（ｈｔｔｐ：／／
ｗｗｗ．ｐｌａｎｔｇｄｂ．ｏｒｇ／）在线区段批量下载工具（ｄｏｗｎｌｏａｄｒｅｇｉｏｎｄａｔａ）下载“一致性”区间的１４４５个预测基因序列并
进行生物信息学分析，发现预测基因主要参与具体的细胞过程，执行结合功能，催化、转移酶活性和氧化还原酶活
性等分子功能。根据“一致性”ＱＴＬ区间的基因位点名称，在ＮＣＢＩ中下载相关基因序列，与所在“一致性”ＱＴＬ区
间所有预测基因保守结构域进行比对，在５个“一致性”持绿 ＱＴＬ区间内初步确定８个持绿相关候选基因。利用
ＧＲＡＭＥＮＥ网站（ｈｔｔｐ：／／ｗｗｗ．ｇｒａｍｅｎｅ．ｏｒｇ／）的Ｃｍａｐ功能，将水稻持绿基因狊犵狉（ｓｔａｙｇｒｅｅｎ）转定位于玉米物理
图谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２上，找到与其同源的玉米候选基因ＧＲＭＺＭ２Ｇ０９１８３７＿Ｔ０１，其序列与已发表的玉米衰老诱导
叶绿体持绿蛋白基因狊犵狉１序列一致。
关键词：玉米；持绿性；数量性状位点（ＱＴＬ）；整合图谱；元分析
中图分类号：Ｓ５１３．０３；Ｑ９４３．２　　文献标识码：Ａ　　文章编号：１００４５７５９（２０１２）０４０１７５１１
　　玉米（犣犲犪犿犪狔狊）是重要的饲料、工业加工原料和粮食作物，其产量和品质对保障粮食安全和人民生活具有
举足轻重的战略作用［１］。持绿性是指在籽粒生理成熟期叶片因衰老进程延缓而保持绿色的特性，被认为是作物
的理想农艺性状［２］。玉米持绿型品种后期绿色叶片数和重量、叶面积指数、叶面积持续期等均优于常规品种，具
有较高的光合活性，较高的蛋白质、蔗糖及脂类含量，较低的纤维含量，其生物产量、籽粒产量和秸秆营养价值也
较高［３５］，并且对一些不良的条件如病害、旱灾、倒伏等具有较强的抗性［６］，是玉米品种的演进方向［７，８］。但由于玉
米的持绿性表现为一个典型的数量性状，而且只有在籽粒到达生理成熟期才能准确鉴定，使得用传统方法选育具
持绿性状的品种进展缓慢［９］。分子数量遗传学和分子标记技术的发展使得剖析数量性状的多基因遗传行为成为
可能，至今在玉米中已定位１７０多个持绿相关数量性状位点（ｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｕｓ，ＱＴＬ）［１０１６］，但由于各实验
所用的材料、作图群体、实验设计、统计方法、性状多少、标记多少、群体大小等互不相同，致使不同实验检测出的
ＱＴＬ数目和座位存在较大差异，虽然这些研究结果在探索玉米持绿性的遗传机制及其与环境的互作上提供了有
益的参考，但在玉米持绿性育种实践中应用仍有限［１７］。
元分析（ｍｅｔａａｎａｌｙｓｉｓ）是一种可以合并不同研究数据进行统计分析且可以对实际数据进行全面检验的方
法［１８］。它可以在整合ＱＴＬ的基础上，建立数学模型优化ＱＴＬ，缩小置信区间，提高ＱＴＬ定位的准确度和有效
性［１９］。随着生物信息学、比较基因组学理论和方法的发展以及水稻（犗狉狔狕犪狊犪狋犻狏犪）、玉米等重要作物全基因组测
序的完成，在全基因组序列水平，功能水平上进行比较基因组研究成为了可能［２０，２１］。２００４年，Ｃｈａｒｄｏｎ等［２２］利用
第２１卷　第４期
Ｖｏｌ．２１，Ｎｏ．４
草　业　学　报
ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ 　　
１７５－１８５
２０１２年８月
收稿日期：２０１２０２２０；改回日期：２０１２０４１８
基金项目：９７３计划前期研究专项（２０１２ＣＢ７２２９０２）和甘肃省科技攻关计划项目（２ＧＳ０６４Ａ４１００１０５）资助。
作者简介：方永丰（１９７６），男，甘肃临洮人，讲师，在读博士。Ｅｍａｉｌ：ｆａｎｇｙｆ＠ｇｓａｕ．ｅｄｕ．ｃｎ
通讯作者。Ｅｍａｉｌ：ｗａｎｇｈｎ＠ｇｓａｕ．ｅｄｕ．ｃｎ，ｓｈａｎｇｘｕｎｗｕ＠１６３．ｃｏｍ
元分析的方法对３１３个玉米花期相关ＱＴＬ信息进行整合分析，最后确定６２个“一致性”ＱＴＬ的准确位置，并通
过与水稻花期性状基因的共线性分析进行验证，发现玉米１９处“一致性”ＱＴＬ区域与水稻花期性状基因紧密相
关。目前，该方法已被广泛应用于各种作物不同性状的整合定位研究［２３，２４］。在玉米上，李雪华等［２５］、栗文娟
等［２６］分别对不同遗传背景的遗传图谱进行整合，构建了玉米抗旱性的一致性遗传图谱；王毅等［２７］、张耿等［２８］、王
帮太等［２９］、王晓丽等［３０］分别构建了玉米产量及构成因子一致性遗传图谱；吉海莲等［３１］构建了玉米抗丝黑穗病一
致性遗传图谱，但对玉米持绿性状ＱＴＬ的一致性遗传图谱的构建及元分析还未见报道。
本研究利用生物信息学方法，对前人定位的玉米持绿相关性状ＱＴＬ进行收集整理，利用ＢｉｏＭｅｒｃａｔｏｒ２．１
软件［１９］，以玉米高密度分子标记连锁图谱ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ为参考图谱，将原始图谱上不同来源的ＱＴＬ映
射到参考图谱上，构建玉米持绿性状ＱＴＬ整合图谱，通过元分析的方法进行“一致性”ＱＴＬ区间的鉴定，结合玉
米基因组序列信息以及生物信息学方法，发掘玉米持绿相关的基因组区域及其候选基因，为深入理解玉米持绿性
状的遗传控制机制和分子设计育种提供科学依据。
１　材料与方法
１．１　实验材料
以玉米基因组数据库（ｈｔｔｐ：／／ｗｗｗ．ｍａｉｚｅＧＤＢ．ｏｒｇ／）和已发表文献中的玉米持绿相关ＱＴＬ定位信息为研
究对象。
１．　实验方法
１．２．１　玉米持绿性相关ＱＴＬ数据的收集　从玉米基因组数据库（ｈｔｔｐ：／／ｗｗｗ．ｍａｉｚｅＧＤＢ．ｏｒｇ／）和已发表的文
献中对玉米成熟期绿色叶片数、成熟期绿叶面积、叶片叶绿素含量、光合速率等持绿相关性状ＱＴＬ定位信息进
行收集，将收集到的每个ＱＴＬ数据按ＢｉｏＭｅｒｃａｔｏｒ２．１软件的要求建立ＱＴＬ数据。ＱＴＬ的位置（最大可能性
的位置及其置信区间）、连锁系数（ｌｏｇａｒｉｔｈｍｏｆｏｄｄｓ，ＬＯＤ）（当ＬＯＤ值≥３．０时，表明标记和ＱＴＬ连锁）和贡献
率犚２（解释表型变异的比率）是ＱＴＬ的３个重要参数，如果某一ＱＴＬ的置信区间未知，则根据Ｄａｒｖａｓｉ和Ｓｏｌ
ｌｅｒ［３２］的公式推断其９５％的置信区间。
犆．犐．＝５３０／（犖×犚２） （１）
犆．犐．＝１６３／（犖×犚２） （２）
式中，犆．犐．指ＱＴＬ的置信区间（ｃｏｎｆｉｄｅｎｃｅｉｎｔｅｒｖａｌ），犖 代表作图群体的大小，犚２ 代表该 ＱＴＬ的贡献率，公式
（１）适用于回交和Ｆ２ 作图群体，公式（２）适用于重组自交系群体。
１．２．２　玉米持绿性相关ＱＴＬ信息的整合　根据收集到的玉米持绿ＱＴＬ信息确定其所涉及的染色体，以ＩＢＭ２
２００８Ｎｅｉｇｈｂｏｒｓ为参考图谱，将原始图谱和参考图谱上相关标记载入ＢｉｏＭｅｒｃａｔｏｒ２．１软件，构建图谱信息库，然
后将ＱＴＬ数据分不同性状载入对应图谱，利用齐序函数（ｈｏｍｏｔｈｅｔｉｃｆｕｎｃｔｉｏｎ）计算共有标记间距，将原始ＱＴＬ
的最大可能性位置及置信区间的两端标记按比例标注到参考图谱上，即映射（ｐｒｏｊｅｃｔｉｏｎ），如果某个ＱＴＬ的任一
位置标记不能映射，则去掉该ＱＴＬ［１９］。
１．２．３　玉米持绿性相关ＱＴＬ的元分析　利用ＢｉｏＭｅｒｃａｔｏｒ２．１软件中的元分析（ｍｅｔａａｎａｌｙｓｉｓ）程序估算“一致
性”ＱＴＬ存在的位置和置信区间。按照软件要求和综合图谱上ＱＴＬ置信区间信息，通过染色体步移法，对不同
性状ＱＴＬ分别进行元分析，对独立来源的同一染色体上相同性状且位于同一座位或有重叠座位的ＱＴＬ计算出
一个“一致性”ＱＴＬ，这个“一致性”ＱＴＬ包括１，２，３，４和 ＮＱＴＬ五个模型，其中可靠性检验值（ａｋａｉｋｅｔｙｐｅ
ｃｒｉｔｅｒｉａｖａｌｕｅｓ，ＡＩＣ）最小的模型为最优模型，即最优“一致性”ＱＴＬ，也是比较接近“真实”ＱＴＬ的模型，并通过高
斯定理最大似然比估算“一致性”ＱＴＬ存在的位置和置信区间［１９］。
１．２．４　玉米持绿性“一致性”ＱＴＬ区间的物理图谱定位及候选基因 ＧｅｎｅＯｎｔｏｌｏｇｙ（ＧＯ）分析　利用 Ｍａｉｚｅ
ＧＤＢ网站（ｈｔｔｐ：／／ｗｗｗ．ｍａｉｚｅＧＤＢ．ｏｒｇ／）中的ＧｅｎｏｍｅＢｒｏｗｓｅｒ工具将玉米物理图谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２序列信
息与ＩＢＭ２２００８连锁图谱的标记信息进行整合，将元分析所得“一致性”ＱＴＬ区间进行物理图谱定位，利用Ｐｌａｎｔ
ＧＤＢ（ｈｔｔｐ：／／ｗｗｗ．ｐｌａｎｔｇｄｂ．ｏｒｇ／）在线区段批量下载工具（ｄｏｗｎｌｏａｄｒｅｇｉｏｎｒａｔａ），下载“一致性”ＱＴＬ区间所
有候选基因序列，利用ＡｇＢａｓ（ｈｔｔｐ：／／ｗｗｗ．ａｇｂａｓｅ．ｍｓｓｔａｔｅ．ｅｄｕ）在线 ＧＯＳｌｉｍＶｉｅｗｅｒ功能，对候选基因进行
６７１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１２） Ｖｏｌ．２１，Ｎｏ．４
ＧｅｎｅＯｎｔｏｌｏｇｙ（ＧＯ）分析［３３］，以了解玉米持绿“一致性”ＱＴＬ区间候选基因所涉及的生物学功能。
１．２．５　玉米持绿性“一致性”ＱＴＬ范围内持绿相关基因位点候选基因发掘　在参考图谱ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ
上，对“一致性”ＱＴＬ区间及其邻近区域的持绿相关基因位点进行整理，根据“一致性”ＱＴＬ区间的基因位点名
称，在ＮＣＢＩ（ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／）中检索，即可得到该基因或相似功能基因的序列信息。利用
ＰｌａｎｔＧＤＢ在线区段批量下载工具，下载“一致性”区间所有蛋白序列，然后，利用Ｐｆａｍ网站（ｈｔｔｐ：／／ｐｆａｍ．ｓａｎ
ｇｅｒ．ａｃ．ｕｋ／）批量搜索工具（ｂａｔｃｈｓｅａｒｃｈ），对所有蛋白序列进行保守结构域分析。最后将这些分析的结果与
ＮＣＢＩ网站提供的基因注释信息进行对比与整合，发掘“一致性”ＱＴＬ区间及其邻近区域的持绿相关基因位点玉
米持绿相关候选基因。
１．２．６　基于基因比较定位的玉米持绿性候选基因发掘　从水稻基因组数据库 ＧＲＡＭＥＮＥ（ｈｔｔｐ：／／ｗｗｗ．
ｇｒａｍｅｎｅ．ｏｒｇ／）检索水稻持绿相关基因信息，利用ＧＲＡＭＥＮＥ网站的Ｃｍａｐ功能，将其转定位于玉米ＩＢＭ２００８
连锁图谱上；利用ＰｌａｎｔＧＤＢ网站在线区段批量下载工具分段下载该区间范围内已预测基因蛋白序列，利用
Ｐｆａｍ网站批量搜索工具，将下载的预测基因蛋白序列与水稻持绿基因的蛋白序列进行保守结构域分析，发掘与
水稻持绿基因同源的玉米持绿候选基因。
２　结果与分析
２．１　玉米持绿性相关ＱＴＬ信息的收集与整理
从玉米基因组数据库和已发表的相关文献中共收集了来源于Ｂ７３×Ｍｏ１７、Ｑ３１９×Ｍｏ１７、Ａ１５０３２×Ｍｏ１７、
黄Ｃ×１７８、４４４×Ｍｏ１７等７个作图群体、涉及４个玉米持绿相关性状的１７３个 ＱＴＬ。包含８１个叶绿素含量
ＱＴＬ，１６个成熟期绿叶数ＱＴＬ，４２个绿叶面积相关ＱＴＬ，３４个叶片光合特性ＱＴＬ（表１）。１７３个ＱＴＬ在１０
条玉米染色体上呈不均匀分布，其中，在第１，４，５染色体上较多，第３，６，１０染色体上较少；叶片叶绿素含量ＱＴＬ
主要分布在第１，４，５，７染色体上，成熟期绿叶面积和叶片光合特性ＱＴＬ主要集中在第１染色体上，而成熟期绿
叶数ＱＴＬ在各染色体上的分布差别不大。
表１　叶片持绿相关性状犙犜犔在各染色体上的分布
犜犪犫犾犲１　犇犻狊狋狉犻犫狌狋犻狅狀狅犳犙犜犔狊犳狅狉犿犪犻狕犲狊狋犪狔犵狉犲犲狀犪狀犱犻狋狊狉犲犾犪狋犲犱狋狉犪犻狋狊狅狀犮犺狉狅犿狅狊狅犿犲狊
性状
Ｔｒａｉｔ
ＱＴＬ类别
ＱＴＬｔｙｐｅ
ＱＴＬ数目
ＱＴＬｎｕｍｂｅｒ
染色体Ｃｈｒｏｍｏｓｏｍｅ
１ ２ ３ ４ ５ ６ ７ ８ ９ １０
叶绿素含量Ｃｈｌｏｒｏｐｈｙｃｏｎｔｅｎｔ ｑｃｈｌ ８１ １０ ５ ２ １６ ２１ ０ １０ ５ ８ ４
绿叶面积 Ｒｅｌａｔｉｖｅｇｒｅｅｎｌｅａｆａｒｅａ ｑｇｌ ４２ １５ ４ ２ ４ ７ ３ ０ ３ ３ １
成熟期绿叶数 Ｇｒｅｅｎｌｅａｆｎｕｍｂｅｒａｔｍａｔｕｒｅｔｉｍｅ ｑｓｔａｇｒ １６ ２ ２ １ １ ３ ２ ０ ２ ２ １
叶片光合特性Ｐｈｏｔｏｓｙｎｔｈｅｔｉｃｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｑｐｎｌ ３４ ８ １ ３ ４ ２ ３ ３ ３ ３ ４
总计 Ｔｏｔａｌ １７３ ３５ １２ ８ ２５ ３３ ８ １３ １３ １６ １０
整合的ＱＴＬＴｏｔａｌｎｕｍｂｅｒｏｆＱＴＬｏｎｒｅｆｅｒｅｎｃｅｍａｐ １５０ ３２ ８ ８ １９ ３３ ８ ４ １３ １６ ９
２．２　玉米持绿性相关ＱＴＬ的图谱整合
将原始图谱与参考图谱ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ上相关标记载入ＢｉｏＭｅｒｃａｔｏｒ２．１软件，利用其图谱映射程序
（Ｍａｐｓｐｒｏｊｅｃｔｉｏｎ）构建玉米持绿相关性状ＱＴＬ整合图谱。当ＬＯＤ值≥３．０时，１７３个ＱＴＬ中有１５０个被成功
整合到ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ图谱，被整合的ＱＴＬ在玉米１０条染色体上均有分布，但分布不均匀（表１），其中，
在第１，４，５，９染色体上集中了６６．７％的ＱＴＬ，其余６条染色体上分布的ＱＴＬ数较少，第７染色体上最少，仅整
合了４个ＱＴＬ。
玉米持绿相关性状ＱＴＬ在整合图谱上有明显成簇分布现象，存在 ＱＴＬ富集区域（图１），在第１染色体的
ｕｍｃ１５９８～ｕｍｃ２３９０、ｕｍｃ１３９５～ｂｎｌｇ４２１和ｕｍｃ１２４５～ｕｍｃ２０４７标记区间分别集中了８，１０和６个ＱＴＬ；第２染
色体的ｕｍｃ１８４５～ｕｍｃ１５４１标记区间集中了５个 ＱＴＬ；第３染色体的ｕｍｃ１３４７～ｕｍｃ１４４９和ｕｍｃ１２６６～
７７１第２１卷第４期 草业学报２０１２年
ｕｍｃ１２８６标记区间分别集中了３和２个ＱＴＬ；第４染色体的ｐｈｉ７２～ｕｍｃ１２８８和ｂｎｌｇ１２６～ｕｍｃ１７０２标记区间分
别集中了８和７个ＱＴＬ；第５染色体的ｕｍｃ１２４０～ｕｍｃ１４９６、ｕｍｃ１５８７～ｕｍｃ１８７０和ｐｈｉ８５～ｕｍｃ１７９２标记区间
分别集中了６，１２和７个ＱＴＬ；第６染色体的ｐｈｉ２６～ｕｍｃ２３１４和ｕｍｃ２３１７～ｕｍｃ１７９５标记区间分别集中了４和
３个ＱＴＬ；第７染色体的４个 ＱＴＬ主要分布在ｂｎｌｇ１６４２～ｕｍｃ１０１５标记区间；第８染色体的ｕｍｃ１１３９～
ｕｍｃ１４８３和ｐｈｉ１１５～ｃｓｕ３１ａ标记区间分别集中了５和６个ＱＴＬ；第９染色体的ｗｘ１～ｕｍｃ１１０７标记区间集中了
１２个ＱＴＬ；第１０染色体的ｂｎｌｇ１４５１～ｕｍｃ１０７７和ｕｍｃ１６７７～ｂｎｌ７．４９ａ标记区间各集中了４个 ＱＴＬ。一些
ＱＴＬ重叠区域与单个性状有关，如集中在第４染色体的ｐｈｉ７２～ｕｍｃ１２８８标记区段内的８个ＱＴＬ均为控制叶片
叶绿素含量的ＱＴＬ；还有一些ＱＴＬ重叠群由多个持绿相关性状ＱＴＬ组成，如第９染色体ｗｘ１～ｕｍｃ１１０７标记
区间集中了６个控制叶片叶绿素含量、２个成熟时绿叶面积、３个叶片光合特性和１个成熟时绿叶数的ＱＴＬ。
图１　第４和第５染色体持绿性相关犙犜犔在犐犅犕２犖犲犻犵犺犫狅狉狊２００８上的整合图谱
犉犻犵．１　犜犺犲犻狀狋犲犵狉犪犾犿犪狆狅犳狊狋犪狔犵狉犲犲狀犙犜犔狊狅犳犿犪犻狕犲狅狀狋犺犲犳狅狌狉狋犺犪狀犱犳犻犳狋犺犮犺狉狅犿狅狊狅犿犲犪狀犪犾狔狕犲犱犫狔犐犅犕２犖犲犻犵犺犫狅狉狊２００８
２．３　玉米持绿性相关ＱＴＬ的元分析
利用ＢｉｏＭｅｒｃａｔｏｒ２．１软件中元分析（ｍｅｔａａｎａｌｙｓｉｓ）程序对整合ＱＴＬ数目较多的第１，４，５，９染色体进行分
析，以每次 Ｍｅｔａ分析中ＡＩＣ值最小，确定１个“一致性”ＱＴＬ（图２），最终得到５个“一致性”ＱＴＬｓ（表２）。
在第１染色体上确定２个“一致性”ＱＴＬ，其中１个控制绿叶面积和光合特性的“一致性”ＱＴＬ（ＣＱＴＬ１），另
１个控制叶绿素含量和绿叶面积的“一致性”ＱＴＬ（ＣＱＴＬ２）（表２）；在第４染色体上确定１个控制叶片叶绿素含
８７１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１２） Ｖｏｌ．２１，Ｎｏ．４
量的“一致性”ＱＴＬ（ＣＱＴＬ３），介于标记ｕｍｃ１２７６～ｕｍｃ１７５７，置信区间为３３．９０～７０．２３ｃＭ，间距３６．３３ｃＭ（表
２）；在第５染色体上得到１个控制叶片叶绿素含量、成熟期绿叶面积及绿叶数的“一致性”ＱＴＬ（ＣＱＴＬ４），置信区
间为２０７．５７～２２２．２８ｃＭ，间距１４．７１ｃＭ，由标记ｕｍｃ１８５２和ｕｍｃ１４６８界定（表２）；在第９染色体上确定１个控
制叶片叶绿素含量和叶片光合特性的“一致性 ＱＴＬ”（ＣＱＴＬ５），置信区间为２２５．６３～２２６．２７ｃＭ，间距０．６４ｃＭ，
由标记ｐｈｉ０６５和ｂｎｌｇ１２７界定（表２）。与原来的ＱＴＬ相比，检测到的５个“一致性”ＱＴＬ的置信区间明显变窄。
图２　第４染色体持绿性相关犙犜犔元分析
犉犻犵．２　犜犺犲犿犲狋犪犪狀犪犾狔狊犻狊狅犳犙犜犔狉犲犾犪狋犲犱狋狅狊狋犪狔犵狉犲犲狀狋狉犪犻狋狊狅犳犿犪犻狕犲狅狀狋犺犲犳狅狌狉狋犺犮犺狉狅犿狅狊狅犿犲
表２　玉米持绿性相关犙犜犔的元分析
犜犪犫犾犲２　犜犺犲犿犲狋犪犪狀犪犾狔狊犻狊狅犳犙犜犔狉犲犾犪狋犲犱狋狅狊狋犪狔犵狉犲犲狀狋狉犪犻狋狊犻狀犿犪犻狕犲
“一致性”ＱＴＬ
ＣｏｎｓｅｎｓｕｓＱＴＬ
染色体
Ｃｈｒｏｍｏｓｏｍｅ
ＡＩＣ值
ＶａｌｕｅｏｆＡＩＣ
位置
Ｐｏｓｉｔｉｏｎ（ｃＭ）
置信区间
Ｃｏｎｆｉｄｅｎｃｅｉｎｔｅｒｖａｌ（ｃＭ）
间距
Ｉｎｔｅｒｖａｌ（ｃＭ）
置信区间标记
Ｍａｒｋｏｎｉｎｔｅｒｖａｌ
ＣＱＴＬ１ １．０４ １３９．６８ ３８６．４４ ３８４．５２～３８８．３７ ３．８５ ｕｍｃ２３９０～ｕｍｃ２２２８
ＣＱＴＬ２ １．０７ １３９．６８ ６６２．３７ ６４７．４７～６７７．２７ ２９．８０ ｕｍｃ１４８６～ｕｍｃ１２４５
ＣＱＴＬ３ ４．０１ ２４７．０２ ５２．０５ ３３．９０～７０．２３ ３６．３３ ｕｍｃ１２７６～ｕｍｃ１７５７
ＣＱＴＬ４ ５．０３ ２６２．１４ ２１４．９２ ２０７．５７～２２２．２８ １４．７１ ｕｍｃ１８５２～ｕｍｃ１４６８
ＣＱＴＬ５ ９．０３ １５３．３３ ２２５．９２ ２２５．６３～２２６．２７ ０．６４ ｐｈｉ０６５～ｂｎｌｇ１２７
２．４　玉米持绿性“一致性”ＱＴＬ区间的物理图谱定位及候选基因ＧｅｎｅＯｎｔｏｌｏｇｙ（ＧＯ）分析
根据“一致性”ＱＴＬ区间两端标记在玉米物理图谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２上的位置，将“一致性”ＱＴＬ进行物理图
谱定位，统计“一致性”ＱＴＬ区间内包含的基因个数（表３）。５个“一致性”ＱＴＬ区间总计包含１４４５个预测基因，
其中ＣＱＴＬ２覆盖范围最大（９．１１Ｍｂ），包含的基因个数最多（７０９个）；ＣＱＴＬ５覆盖范围最小，仅为０．１９Ｍｂ，其
包含的基因个数也最少（２８个）。单位长度基因个数最少的是ＣＱＴＬ１，为６１．３０个基因，单位长度基因个数最多
的是ＣＱＴＬ５，包含１４７．４个基因。
利用ＰｌａｎｔＧＤＢ在线区段批量下载工具（ｄｏｗｎｌｏａｄｒｅｇｉｏｎｄａｔａ），下载５个“一致性”ＱＴＬ区间所有基因序
列，利用ＡｇＢａｓ网站（ｈｔｔｐ：／／ｗｗｗ．ａｇｂａｓｅ．ｍｓｓｔａｔｅ．ｅｄｕ）的ＧＯＳｌｉｍＶｉｅｗｅｒ对候选基因进行ＧｅｎｅＯｎｔｏｌｏｇｙ分
析，发现１４４５个候选基因中，有５１３个基因参与植物体内的生物过程，１６５个基因是细胞组成成分，７６７个基因
具有分子功能，分别涉及２２种生物过程，２３种分子功能，１１种细胞成分。具有分子功能的基因中最主要涉及的
类是：结合功能，催化活性，转移酶活性以及氧化还原酶活性（图３Ａ）；从细胞组成分类来看，细胞及细胞组成成分
所占比例最大，其次为胞内物质、细胞质及细胞膜（图３Ｂ）；所参与的生物学过程按照比例从大到小依次包括：细
胞过程、代谢过程、大分子代谢过程、生物学过程、刺激的应答反应、核酸代谢过程、生物合成、细胞间通讯、运输、
分解代谢过程、细胞死亡、细胞内氨基酸代谢过程等（图３Ｃ）。
９７１第２１卷第４期 草业学报２０１２年
表３　“一致性”犙犜犔区间的物理图谱定位
犜犪犫犾犲３　犜犺犲犮狅狀狊犲狀狊狌狊犙犜犔狉犲犵犻狅狀狅狀犿犪犻狕犲狆犺狔狊犻犮犪犾犿犪狆狅犳犅７３
“一致性”ＱＴＬ
ＣｏｎｓｅｎｓｕｓＱＴＬ
染色体
Ｃｈｒｏｍｏｓｏｍｅ
物理图谱距离
Ｉｎｔｅｒｖａｌｏｆｐｈｙｓｉｃａｌｍａｐ（Ｍｂ）
间距
Ｉｎｔｅｒｖａｌ（Ｍｂ）
ＱＴＬ范围内基因个数
Ｇｅｎｅｓｎｕｍｂｅｒｉｎｉｎｔｅｒｖａｌ
ＣＱＴＬ１ １．０４ ７１．０２～７３．１９ ２．１７ １３３
ＣＱＴＬ２ １．０７ ２０５．７３～２１４．８５ ９．１１ ７０９
ＣＱＴＬ３ ４．０１ ２．８８～４．６８ １．８０ ２２３
ＣＱＴＬ４ ５．０３ １６．１５～２０．９０ ４．７５ ３５２
ＣＱＴＬ５ ９．０３ ２６．８２～２７．０１ ０．１９ ２８
总计Ｔｏｔａｌ １４４５
利用ＰｌａｎｔＧＤＢ在线区段批量下载工具（ｄｏｗｎｌｏａｄｒｅｇｉｏｎｄａｔａ），下载“一致性”ＱＴＬ区段的１４４５个候选基
因对应的蛋白序列，利用Ｐｆａｍ在线批量搜索工具（ｂａｔｃｈｓｅａｒｃｈ），对其蛋白序列进行保守结构域分析，发现在５
个“一致性”ＱＴＬ区段分别包含ＧＲＡＳ、ＣＢＦ、ＳＲＦ转录调控因子、叶绿体结合蛋白、植物激素诱导（响应）蛋白及
细胞分裂周期相关蛋白、细胞色素Ｐ４５０氧化酶、细胞分裂素脱氢酶、磷脂酶ＰＬＤ、ＲｕＢＰ羧化酶、糖基转移酶等已
报道的与持绿性相关的功能蛋白或转录因子家族的保守结构域。
图３　“一致性”犙犜犔区间候选基因犌犗分析
犉犻犵．３　犌犲狀犲犗狀狋狅犾狅犵狔犪狀犪犾狔狊犻狊狅犳犵犲狀犲狊狉犲犾犪狋犲犱狋狅犮狅狀狊犲狀狊狌狊犙犜犔
　Ａ：参与细胞组成的基因 Ｇｅｎｅｓｃｏｎｃｅｒｎｅｄｗｉｔｈｃｅｌｕｌａｒｃｏｍｐｏｎｅｎｔ；Ｂ：具有分子功能的基因 Ｍｏｌｅｃｕｌａｒｆｕｎｃｔｉｏｎｇｅｎｅ；Ｃ：参与生物过程的基因
Ｇｅｎｅｓｃｏｎｃｅｒｎｅｄｗｉｔｈｂｉｏｌｏｇｉｃａｌｐｒｏｃｅｓｓ．
２．５　玉米持绿性“一致性”ＱＴＬ范围内持绿相关基因位点候选基因发掘
在参考图谱ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ上，对“一致性”ＱＴＬ区间及其邻近区域的持绿相关基因位点进行整理，
发现７个与玉米持绿性有关的基因（表４），涉及到信号传递、胁迫应答、细胞组成、代谢调节等多种生理生化途
径。根据“一致性”ＱＴＬ区段的基因名称，在ＮＣＢＩ中搜索，即可得到该基因或相似功能基因的序列信息，下载后
在Ｐｆａｍ 网站搜索其保守结构域信息，与所在“一致性”ＱＴＬ区段的所有基因序列保守结构域进行比对发现，
０８１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１２） Ｖｏｌ．２１，Ｎｏ．４
ＣＱＴＬ１区间的基因 犿狊狏１与该区段内的 ＧＲＭＺＭ２Ｇ１７２３２２＿Ｔ０１基因序列保守结构域相同，序列同源比为
７４％；犵狊狋４２基因与该区段内的 ＧＲＭＺＭ２Ｇ３９８３５７＿Ｔ０４基因序列保守结构域相同，且序列相似度高达９５％。
ＣＱＴＬ２区间狅狆狉７基因与该区段的ＧＲＭＺＭ２Ｇ１４８２８１＿Ｔ０２基因序列一致，且保守结构域相同；犮狀犮狉１基因与该
区段的ＧＲＭＺＭ２Ｇ１３１２０５＿Ｔ０１基因序列保守结构域相同，且序列相似度达９０％。ＣＱＴＬ３区间犫狓３和犫狓４两个
基因分别与该区段的ＧＲＭＺＭ２Ｇ１６７５４８＿Ｔ０２和ＧＲＭＺＭ２Ｇ４７３４８０＿Ｔ０２基因序列一致，且保守结构域相同。
ＣＱＴＬ４区间狋犺犪９基因与该区段的ＧＲＭＺＭ２Ｇ８０２０２９＿Ｔ０１基因序列一致，且保守结构域相同。ＣＱＴＬ５区间
狆犲狆１基因与该区段的序列ＧＲＭＺＭ２Ｇ０８３８４１＿Ｔ０２保守结构域相同，序列同源比为８４％。
表４　“一致性”犙犜犔区间基因与玉米基因组预测基因序列的一致性
犜犪犫犾犲４　犆狅犺犲狉犲狀犮犲狅犳犵犲狀犲犾狅犮犻犪狀犱狆狉犲犱犻犮狋犲犱犵犲狀犲狊犳狉狅犿犿犪犻狕犲狊犲狇狌犲狀犮犲犱犪狋犪犫犪狊犲犻狀犮狅狀狊犲狀狊狌狊犙犜犔狉犲犵犻狅狀
“一致性”ＱＴＬ
ＣｏｎｓｅｎｓｕｓＱＴＬ
基因
Ｇｅｎｅ
预测基因号
ＧｅｎｅＩＤｆｒｏｍＧｒａｍｅｎｅ
保守结构域是否相同
Ｗｉｔｈ（Ｙ）ｏｒｗｉｔｈｏｕｔ（Ｎ）ｅｑｕａｌｃｏｎｓｅｒｖｅｄｄｏｍａｉｎ
序列同源比例
Ｃｏｎｓｅｒｖｅｄｐｒｏｐｏｒｔｉｏｎ（％）
ＣＱＴＬ１ 犿狊狏１ ＧＲＭＺＭ２Ｇ１７２３２２＿Ｔ０１ Ｙ ７４
犵狊狋４２ ＧＲＭＺＭ２Ｇ３９８３５７＿Ｔ０４ Ｙ ９５
ＣＱＴＬ２ 狅狆狉７ ＧＲＭＺＭ２Ｇ１４８２８１＿Ｔ０２ Ｙ １００
犮狀犮狉１ ＧＲＭＺＭ２Ｇ１３１２０５＿Ｔ０１ Ｙ ９０
ＣＱＴＬ３ 犫狓３ ＧＲＭＺＭ２Ｇ１６７５４８＿Ｔ０２ Ｙ １００
犫狓４ ＧＲＭＺＭ２Ｇ４７３４８０＿Ｔ０２ Ｙ １００
ＣＱＴＬ４ 狋犺犪９ ＧＲＭＺＭ２Ｇ８０２０２９＿Ｔ０１ Ｙ １００
ＣＱＴＬ５ 狆犲狆１ ＧＲＭＺＭ２Ｇ０８３８４１＿Ｔ０２ Ｙ ８４
　犿狊狏１：玉米条纹病抗性基因 Ｍａｉｚｅｓｔｒｅａｋｖｉｒｕｓｔｏｌｅｒａｎｃｅｇｅｎｅ１；犵狊狋４２：谷胱甘肽Ｓ转移酶基因ＧｌｕｔａｔｈｉｏｎｅＳｔｒａｎｓｆｅｒａｓｅｇｅｎｅ４２；狅狆狉７：１２羰植物
二烯酸还原酶基因１２ｏｘｏｐｈｙｔｏｄｉｅｎｏｉｃａｃｉｄｒｅｄｕｃｔａｓｅｇｅｎｅ７；犮狀犮狉１：肉桂酰辅酶 Ａ还原酶基因 ＣｉｎｎａｍｏｙｌＣｏＡｒｅｄｕｃｔａｓｅｇｅｎｅ；犫狓３，犫狓４：细胞色素
Ｐ４５０氧化酶基因Ｂｅｎｚｏｘａｚｉｎｏｎｅｓｙｎｔｈｅｓｉｓｇｅｎｅ３，ｇｅｎｅ４；狋犺犪９：叶绿体类囊体构成基因 Ｔｈｙｌａｋｏｉｄａｓｓｅｍｂｌｙｇｅｎｅ９；狆犲狆１：磷酸烯醇式丙酮酸羧化酶
基因Ｐｈｏｓｐｈｏｅｎｏｌｐｙｒｕｖａｔｅｃａｒｂｏｘｙｌａｓｅｇｅｎｅ１．
２．６　基于基因比较定位的玉米持绿性候选基因发掘
从ＧＲＡＭＥＮＥ网站检索到１个位于水稻遗传图ＩＧＣＮ１９９８第９染色体１３１．１～１３４．７ｃＭ区段的衰老诱导
叶绿体持绿蛋白基因（ｓｅｎｅｓｃｅｎｃｅｉｎｄｕｃｉｂｌｅｃｈｌｏｒｏｐｌａｓｔｓｔａｙｇｒｅｅｎｐｒｏｔｅｉｎ）狊犵狉（ＬＯＣ＿Ｏｓ０９ｇ３６２００．１）（ＧｅｎＢａｎｋ：
ＡＹ８５０１３４．１）。利用ＧＲＡＭＥＮＥ网站的Ｃｍａｐ功能，将该基因成功转定位于玉米ＩＢＭ２００８图谱的第７染色体
上，位于ｒｚ６１７～ｒｚ４０４标记区间，区间长度为１７９．４ｃＭ（图４Ａ）；利用 ＭａｉｚｅＧＤＢ网站中的ＧｅｎｏｍｅＢｒｏｗｓｅｒ工
具将Ｂ７３ＲｅｆＧｅｎ＿ｖ２序列信息与ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ标记信息进行整合，将水稻狊犵狉基因定位在玉米物理图
谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２的９４．４２～１４３．４１Ｍｂ范围，利用ＰｌａｎｔＧＤＢ网站（ｈｔｔｐ：／／ｗｗｗ．ｐｌａｎｔｇｄｂ．ｏｒｇ／）的在线区段
批量下载工具（ｄｏｗｎｌｏａｄｒｅｇｉｏｎｄａｔａ），分段下载该区间范围内已预测基因蛋白序列，与水稻持绿狊犵狉（ＬＯＣ＿
Ｏｓ０９ｇ３６２００．１）基因的蛋白序列进行保守结构域分析，最后在玉米第７染色体上（ｃｈｒ７：１４３０２０４５２Ｍｂ～
１４３０２２３７２Ｍｂ）找到与水稻狊犵狉基因同源的玉米基因ＧＲＭＺＭ２Ｇ０９１８３７＿Ｔ０１（图４Ｂ），其编码蛋白序列与水稻
狊犵狉基因一致（图４Ｃ）。
进一步将该预测基因序列与ＮＣＢＩ网站（ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／）提供的基因注释信息进行比对，
发现预测基因ＧＲＭＺＭ２Ｇ０９１８３７＿Ｔ０１与ＧｅｎＢａｎｋ中的玉米衰老诱导叶绿体持绿蛋白基因（ｓｅｎｅｓｃｅｎｃｅｉｎｄｕｃ
ｉｂｌｅｃｈｌｏｒｏｐｌａｓｔｓｔａｙｇｒｅｅｎｐｒｏｔｅｉｎ）狊犵狉１（ＡＹ８５０１３８．１）序列一致，表明ＧＲＭＺＭ２Ｇ０９１８３７＿Ｔ０１为玉米叶片持绿
的一个候选基因。
３　讨论
叶片衰老是一个受特定基因调控的主动的衰退过程，受不同生理生化途径的调控。对玉米持绿相关的４个
性状的１７３个ＱＴＬ进行整理，当ＬＯＤ值≥３．０时，１７３个ＱＴＬ中有１５０个被成功整合到ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ
１８１第２１卷第４期 草业学报２０１２年
图４　水稻持绿基因（狊犵狉）在玉米图谱的比较定位及其蛋白序列比对
犉犻犵．４　犆狅犾犻狀犲犪狉犻狋狔犮狅犿狆犪狉犲犱狋犺犲狊狋犪狔犵狉犲犲狀犵犲狀犲狊（狊犵狉）犻狀狏狅犾狏犻狀犵犻狀狉犻犮犲犪狀犱犿犪犻狕犲犪狀犱狋犺犲犻狉狆狉狅狋犲犻狀狊犲狇狌犲狀犮犲犪狀犪犾狔狊犻狊
　Ａ：在玉米ＩＢＭ２Ｎｅｉｇｈｂｏｒｓ２００８遗传图谱上水稻持绿基因（狊犵狉）的比较定位 Ｃｏｌｉｎｅａｒｉｔｙｃｏｍｐａｒｅｄｔｈｅｓｔａｙｇｒｅｅｎｇｅｎｅ（狊犵狉）ｉｎｖｏｌｖｉｎｇｉｎｒｉｃｅ
ａｎｄｍａｉｚｅｏｎＩＢＭ２Ｎｅｉｇｈｂｏｒｓ２００８ｍａｐ；Ｂ：水稻持绿基因（狊犵狉）在玉米物理图谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２上的比较定位 Ｍａｐｅｄｔｈｅｒｉｃｅｓｔａｙｇｒｅｅｎｇｅｎｅ
（狊犵狉）ｏｎＢ７３ＲｅｆＧｅｎ＿ｖ２；Ｃ：持绿基因（狊犵狉）蛋白序列同源比对结果 Ｔｈｅｂｌａｓｔｐｒｅｓｕｌｔｓｏｆｓｔａｙｇｒｅｅｎｇｅｎｅｐｒｏｔｅｉｎｓｅｑｕｅｎｃｅｏｆｒｉｃｅａｎｄｍａｉｚｅ．
高密度连锁图谱上，被整合的ＱＴＬ在１０条染色体上均有分布，这一结果显示，玉米的持绿性作为一个复杂的数
量性状，受众多基因的控制，作用机理复杂，涉及的生理调控途径也较多。Ｔｕｂｅｒｏｓａ等［３４］和严建兵等［３５］均发现
了ＱＴＬ聚集现象，本研究也发现不同材料涉及相同性状的ＱＴＬ在同一区域聚集，如集中在第４染色体的ｐｈｉ７２
～ｕｍｃ１２８８标记区段内的８个ＱＴＬ均为控制叶片叶绿素含量的ＱＴＬ，说明ＱＴＬ可能不是随机分布，而是位于
染色体的特定区域，预示着该区域的确存在“真实”的ＱＴＬ，这也是本研究仅对玉米持绿ＱＴＬ相对集中的第１，
４，５，９染色体进行了详细分析的原因。而且涉及不同性状的 ＱＴＬ也有聚集的现象，如第９染色体 ｗｘ１～
ｕｍｃ１１０７标记区间集中了６个控制叶片叶绿素含量、２个成熟时绿叶面积、３个叶片光合特性和１个成熟时绿叶
数的ＱＴＬ。这预示着在这一染色体区段可能密集着玉米持绿相关的ＱＴＬ簇，即基因簇，这种现象进一步说明了
不同性状之间的相互关系可归因于基因多效性或是基因之间的紧密连锁，为下一步基于一致性ＱＴＬ进行持绿
性相关的关键基因的定位和发掘提供了依据。
Ａｌｅｘａｎｄｒｏｖ等［３６］对玉米大规模ｃＤＮＡ测序结果分析以及Ｃａｒｏｌ等［３７］对玉米基因组２７４５５条全长ｃＤＮＡ的
ＧＯ注释结果分析表明，玉米基因组生物过程分类中，参与代谢过程的基因比例最大；细胞成分分类中，构成膜系
统基因所占比例最大；分子功能分类中，执行结合功能的基因所占比例最大。本研究中参与代谢过程、执行结合
功能的基因所占比例与整个基因组基因的ＧＯ分类结果一致，而细胞、细胞组成成分以及胞内物质所占比例远高
于膜系统所占比例，而且生物学过程中的细胞过程、大分子代谢过程、核酸代谢过程、细胞间通讯、运输、分解代谢
过程、细胞死亡、细胞内氨基酸代谢等过程所占比例远远高于基因组基因ＧＯ分布中的相应比例。由此可见，本
２８１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１２） Ｖｏｌ．２１，Ｎｏ．４
研究的１４４５个持绿候选基因主要参与了具体的细胞过程，执行结合功能，催化、转移酶活性和氧化还原酶活性
等分子功能。Ｗｏｎｇ等［３８］发现，控制类胡萝卜素４种组分合成的几个 ＱＴＬ与相关基因狔１和狏狆９位置重叠。
Ｚｈａｎｇ等［３９］发现，控制玉米花丝内可凝性球蛋白水平的一个ＱＴＬ与２个相关候选基因狆１和狆２位置重叠。本
研究的５个“一致性”ＱＴＬ区间在ＩＢＭ２２００８Ｎｅｉｇｈｂｏｒｓ图谱内都检测到持绿相关基因，在包含持绿相关基因的
ＭＱＴＬ区域，存在与这些基因序列相同或相似的预测基因序列，最终在５个玉米持绿“一致性”ＱＴＬ内，初步确
定８个玉米持绿性候选基因，下一步将验证其功能。
禾本科作物基因组存在高度保守性，控制重要性状的基因在不同作物中可能位于相同或相近的位置，并与相
同的分子标记相连锁，在某一作物中定位的基因或ＱＴＬ信息可以用于其他作物［２２］。水稻由于其基因组相对较
小，遗传转化体系成熟，且与禾本科粮食作物存在共线性，使其成为作物分子遗传学及基因组学研究的模式植物，
已经有一大批的控制水稻高产、优质、抗逆和营养高效等许多重要农艺性状的基因被成功分离［４０］。本研究通过
比较定位将水稻衰老诱导叶绿体持绿蛋白基因狊犵狉转定位于在玉米物理图谱Ｂ７３ＲｅｆＧｅｎ＿ｖ２上，采用染色体步
移的方法，找到与水稻狊犵狉基因同源的玉米基因ＧＲＭＺＭ２Ｇ０９１８３７＿Ｔ０１，将该基因序列与ＮＣＢＩ提供的基因注
释信息进行比对，发现与玉米衰老诱导叶绿体持绿蛋白基因狊犵狉１序列一致，为候选基因发掘提供了一条新途径。
４　结论
采用ＢｉｏＭｅｒｃａｔｏｒ２．１软件构建出含１５０个玉米持绿性相关ＱＴＬ的整合图谱，１５０个ＱＴＬ在１０条染色体
上均有分布，但存在明显富集和成簇分布的现象。对玉米持绿ＱＴＬ相对集中的第１，４，５，９染色体进行元分析，
得到５个玉米持绿“一致性”ＱＴＬｓ。持绿“一致性”ＱＴＬ区间共有１４４５个基因，其中，５１３个基因参与植物体内
的生物过程，１６５个基因是细胞组成成分，７６７个基因具有分子功能，分别涉及２２种生物过程，２３种分子功能，１１
种细胞成分。候选基因主要参与了具体的细胞过程，主要执行结合功能，催化、转移酶活性和氧化还原酶活性等
分子功能，包含ＧＲＡＳ、ＣＢＦ、ＳＲＦ转录调控因子、叶绿体结合蛋白、植物激素诱导／响应蛋白及细胞分裂周期相关
蛋白、细胞色素Ｐ４５０氧化酶、细胞分裂素脱氢酶、磷脂酶ＰＬＤ、ＲｕＢＰ羧化酶、糖基转移酶等已报道的与持绿性相
关的功能蛋白或转录因子家族的保守结构域。在包含持绿性相关基因位点的“一致性”ＱＴＬ区域，均存在与这些
基因序列相同或相似的预测基因，最终，在５个“一致性”持绿ＱＴＬ区间内初步确定了８个玉米持绿性相关候选
基因。利用比较基因组学方法，将水稻衰老诱导叶绿体持绿蛋白基因狊犵狉转定位于在玉米物理图谱Ｂ７３ＲｅｆＧｅｎ
＿ｖ２上，并找到与该基因同源的玉米候选基因ＧＲＭＺＭ２Ｇ０９１８３７＿Ｔ０１，为候选基因发掘提供了一条新途径。
参考文献：
［１］　乔岩，王汉宁，张成，等．玉米胚乳突变基因ａｅ连锁累赘的ＳＳＲ分析［Ｊ］．草业学报，２０１１，２０（１）：１４０１４７．
［２］　ＴｈｏｍａｓＨ，ＳｍａｒｔＣＭ．Ｃｒｏｐｓｔｈａｔｓｔａｙｇｒｅｅｎ［Ｊ］．ＡｎｎａｌｓｏｆＡｐｐｌｉｅｄＢｉｏｌｏｇｙ，１９９３，１２３：１９３２１９．
［３］　ＧｅｎｔｉｎｅｔｔａＥ，ＣｅｐｐｉＤ，ＬｅｐｏｒｉＣ，犲狋犪犾．Ａｍａｊｏｒｇｅｎｅｆｏｒｄｅｌａｙｅｄｓｅｎｅｓｃｅｎｃｅｉｎｍａｉｚｅ：Ｐａｔｔｅｒｎｏｆｐｈｏｔｏｓｙｎｔｈａｔｅｓａｃｃｕｍｕｌａｔｉｏｎ
ａｎｄｉｎｈｅｒｉｔａｎｃｅ［Ｊ］．ＰｌａｎｔＢｒｅｅｄ，１９８６，９７：１９３２０３．
［４］　ＭａＢＬ，ＤｗｙｅｒＬＭ．Ｎｉｔｒｏｇｅｎｕｐｔａｋｅａｎｄｕｓｅｏｆｔｗｏｃｏｎｔｒａｓｔｉｎｇｍａｉｚｅｈｙｂｒｉｄｓｄｉｆｆｅｒｉｎｇｉｎｌｅａｆｓｅｎｅｓｃｅｎｃｅ［Ｊ］．ＰｌａｎｔａｎｄＳｏｉｌ，
１９９８，１９９：２８３２９１．
［５］　ＰｏｍｍｅｌＢ，ＧａｌａｉｓＡ，ＣｏｑｕｅＭ，犲狋犪犾．Ｃａｒｂｏｎａｎｄｎｉｔｒｏｇｅｎａｌｏｃａｔｉｏｎａｎｄｇｒａｉｎｆｉｌｉｎｇｉｎｔｈｒｅｅｍａｉｚｅｈｙｂｒｉｄｓｄｉｆｆｅｒｉｎｇｉｎｌｅａｆ
ｓｅｎｅｓｃｅｎｃｅ［Ｊ］．ＥｕｒｏｐｅａｎＪｏｕｒｎａｌｏｆＡｇｒｏｎｏｍｙ，２００６，２４：２０３２１１．
［６］　ＢｅｋａｖａｃＧ，ＰｕｒａｒＢ，ＳｔｏｊａｋｏｖｉｈＭ，犲狋犪犾．Ｇｅｎｅｔｉｃａｎａｌｙｓｉｓｏｆｓｔａｙｇｒｅｅｎｔｒａｉｔｉｎｂｒｏａｄｂａｓｅｄｍａｉｚｅｐｏｐｕｌａｔｉｏｎｓ［Ｊ］．ＣｅｒｅａｌＲｅ
ｓｅａｒｃｈＣｏｍｍｕｎｉｃａｔｉｏｎｓ，２００７，３５（１）：１９３２０３．
［７］　董树亭，王空军，胡昌浩，等．玉米品种更替过程中群体光合特性的演变［Ｊ］．作物学报，２０００，２６（２）：２００２０４．
［８］　ＤｕｖｉｃｋＤＮ．Ｇｅｎｅｔｉｃｒａｔｅｓｏｆｇｒａｉｎｉｎｈｙｂｒｉｄｍａｉｚｅｙｉｅｌｄｓｄｕｒｉｎｇｔｈｅｐａｓｔ４０ｙｅａｒｓ［Ｊ］．Ｍａｙｄｉｃａ，１９９７，２３：１８７１９６．
［９］　郑洪建，沈雪芳，吴爱忠．玉米叶片保绿性遗传研究进展［Ｊ］．上海农业学报，２００７，２３（４）：９０９４．
［１０］　ＢｅａｖｉｓＷＤ，ＳｍｉｔｈＯＳ，ＧｒａｎｔＤ，犲狋犪犾．ＩｄｅｎｔｉｆｉｃａｔｉｏｎｏｆｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｉｕｓｉｎｇａｓｍａｌｓａｍｐｌｅｏｆｔｏｐｃｒｏｓｓｅｄａｎｄＦ４
ｐｒｏｇｅｎｙｆｒｏｍｍａｉｚｅ［Ｊ］．ＣｒｏｐＳｃｉｅｎｃｅ，１９９４，３４：８８２８９６．
［１１］　ＶｉｅｎｎｅＤ，ＬｅｏｎａｒｄｉＡ．ＧｅｎｅｔｉｃｓｏｆｐｒｏｔｅｏｍｅｖａｒｉａｔｉｏｎｆｏｒＱＴＬｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ：ａｐｐｌｉｃａｔｉｏｎｔｏｄｒｏｕｇｈｔｓｔｒｅｓｓｒｅｓｐｏｎｓｅｓｉｎ
３８１第２１卷第４期 草业学报２０１２年
ｍａｉｚｅ［Ｊ］．ＪｏｕｒｎａｌｏｆＥｘｐｅｒｉｍｅｎｔａｌＢｏｔａｎｙ，１９９９，５０：３０３３０９．
［１２］　ＺｈｅｎｇＨＪ，ＷｕＡＺ，ＣａｉＲ，犲狋犪犾．ＱＴＬｍａｐｐｉｎｇｏｆｍａｉｚｅ（犣犲犪犿犪狔狊Ｌ．）ｓｔａｙｇｒｅｅｎｔｒａｉｔｓａｎｄｔｈｅｉｒｒｅｌａｔｉｏｎｓｈｉｐｔｏｙｉｅｌｄ［Ｊ］．Ｐｌａｎｔ
Ｂｒｅｅｄｉｎｇ，２００９，１２８（１）：５４６２．
［１３］　王爱玉，张春庆．玉米叶绿素含量的ＱＴＬ定位［Ｊ］．遗传，２００８，３０（８）：１０８３１０９１．
［１４］　王春丽．玉米穗部性状及主要农艺性状发育动态的ＱＴＬ定位［Ｄ］．郑州：河南农业大学，２００５．
［１５］　刘宗华，谢惠玲，王春丽，等．氮胁迫和非胁迫条件下玉米不同时期叶绿素含量的 ＱＴＬ分析［Ｊ］．植物营养与肥料学报，
２００８，１４（５）：８４５８５１．
［１６］　于晓芳．玉米叶片水分利用效率及其相关性状的ＱＴＬ定位与分析［Ｄ］．包头：内蒙古农业大学，２０１０．
［１７］　王建康，李慧慧，张学才，等．中国作物分子设计育种［Ｊ］．作物学报，２０１１，３７（２）：１９１２０１．
［１８］　ＧｏｆｆｉｎｅｔＢ，ＧｅｒｂｅｒＳ．Ｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｉ：Ａｍｅｔａａｎａｌｙｓｉｓ［Ｊ］．Ｇｅｎｅｔｉｃｓ，２０００，１５５：４６３４７３．
［１９］　ＡｒｃａｄｅＡ，ＬａｂｏｕｒｄｅｔｔｅＡ，ＦａｌｑｕｅＭ，犲狋犪犾．ＢｉｏＭｅｒｃａｔｏｒ：ｉｎｔｅｇｒａｔｉｎｇｇｅｎｅｔｉｃｍａｐｓａｎｄＱＴＬｔｏｗａｒｄｓｄｉｓｃｏｖｅｒｙｏｆｃａｎｄｉｄａｔｅ
ｇｅｎｅｓ［Ｊ］．Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ，２００４，２０（１４）：２３２４２３２６．
［２０］　李立芹，黄玉碧，王西瑶．马铃薯 ＷＲＫＹ６基因的克隆、序列分析与原核表达研究［Ｊ］．草业学报，２０１１，２０（２）：１７７１８３．
［２１］　江腾，林勇祥，刘雪，等．苜蓿全基因组 ＷＲＫＹ转录因子基因的分析［Ｊ］．草业学报，２０１１，２０（３）：２１１２１８．
［２２］　ＣｈａｒｄｏｎＦ，ＶｉｒｌｏｎＢ，ＭｏｒｅａｕＬ，犲狋犪犾．Ｇｅｎｅｔｉｃａｒｃｈｉｔｅｃｔｕｒｅｏｆｆｌｏｗｅｒｉｎｇｔｉｍｅｉｎｍａｉｚｅａｓｉｎｆｅｒｒｅｄｆｒｏｍｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｉ
ｍｅｔａａｎａｌｙｓｉｓａｎｄｓｙｎｔｅｎｙｃｏｎｓｅｒｖａｔｉｏｎｗｉｔｈｔｈｅｒｉｃｅｇｅｎｏｍｅ［Ｊ］．Ｇｅｎｅｔｉｃｓ，２００４，１６８：２１６９２１８５．
［２３］　ＬｏｍｂａｒｄＶ，ＤｅｌｏｕｒｍｅＲ．Ａｃｏｎｓｅｎｓｕｓｌｉｎｋａｇｅｍａｐｆｏｒｒａｐｅｓｅｅｄ：Ｃｏｎｓｔｒｕｃｔｉｏｎａｎｄｉｎｔｅｇｒａｔｉｏｎｏｆｔｈｒｅｅｉｎｄｉｖｉｄｕａｌｍａｐｓｆｒｏｍ
ＤＨｐｏｐｕｌａｔｉｏｎｓ［Ｊ］．ＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，２００１，１０３：４９１５０７．
［２４］　ＭａｒｃｅｌＴＣ，ＶａｒｓｈｎｅｙＲＫ，ＢａｒｂｉｅｒｉＭ，犲狋犪犾．ＡｈｉｇｈｄｅｎｓｉｔｙｃｏｎｓｅｎｓｕｓｍａｐｏｆｂａｒｌｅｙｔｏｃｏｍｐａｒｅｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆＱＴＬｓ
ｆｏｒｐａｒｔｉａｌｒｅｓｉｓｔａｎｃｅｔｏＰｕｃｃｉｎｉａｈｏｒｄｅｉａｎｄｏｆｄｅｆｅｎｃｅｇｅｎｅｈｏｍｏｌｏｇｕｅｓ［Ｊ］．ＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，２００７，１１４：
４８７５００．
［２５］　李雪华，李新海，郝转芳，等．干旱条件下玉米耐旱相关性状的 ＱＴＬ一致性图谱构建［Ｊ］．中国农业科学，２００５，３８（５）：
８８２８９０．
［２６］　栗文娟，刘志斋，石云素，等．基于元分析和生物信息学分析的玉米抗旱相关性状 ＱＴＬ一致性区间定位［Ｊ］．作物学报，
２０１０，３６（９）：１４５７１４６７．
［２７］　王毅，姚骥，张征锋，等．基于玉米综合ＱＴＬ图谱的比较分析及株高ＱＴＬ的统合分析［Ｊ］．科学通报，２００６，５１（１５）：１７７６
１７８６．
［２８］　张耿，祝丽英，陈景堂，等．玉米产量相关性状 ＱＴＬ通用图谱的构建及功能候选基因的开发［Ｊ］．华北农学报，２００８，
２３（６）：２２２７．
［２９］　王帮太，吴建宇，丁俊强，等．玉米产量及产量相关性状ＱＴＬ的图谱整合［Ｊ］．作物学报，２００９，３５（１０）：１８３６１８４３．
［３０］　王晓丽，李新海，王振华，等．玉米产量因子ＱＴＬ整合图谱构建与“一致性”ＱＴＬ确定［Ｊ］．核农学报，２００８，２２（６）：７５６
７６１．
［３１］　吉海莲，李新海，谢传晓，等．基于元分析的抗玉米丝黑穗病 ＱＴＬ比较定位［Ｊ］．植物遗传资源学报，２００７，８（２）：１３２
１３９．
［３２］　ＤａｒｖａｓｉＡ，ＳｏｌｅｒＭ．ＡｓｉｍｐｌｅｍｅｔｈｏｄｔｏｃａｌｃｕｌａｔｅｒｅｓｏｌｖｉｎｇｐｏｗｅｒａｎｄｃｏｎｆｉｄｅｎｃｅｉｎｔｅｒｖａｌｏｆＱＴＬｍａｐｌｏｃａｔｉｏｎ［Ｊ］．Ｂｅｈａｖ
ｉｏｒＧｅｎｅｔｉｃｅｓ，１９９７，２７：１２５１３２．
［３３］　ＡｓｈｂｕｍｅｒＭ．Ｇｅｎｅｏｎｔｏｌｏｇｙ：ｔｏｏｌｆｏｒｔｈｅｕｎｉｆｉｅａｔｉｏｎｏｆｂｉｏｌｏｇｙ．ＴｈｅＧｅｎｅＯｎｔｏｌｏｇｙＣｏｎｓｏｒｔｉｕｍ［Ｊ］．ＮａｔｕｒｅＧｅｎｅｔ，２０００，
２５（１）：２５２９．
［３４］　ＴｕｂｅｒｏｓａＲ，ＳａｌｖｉＳ，ＳａｎｇｕｉｎｅｔｉＭ，犲狋犪犾．ＭａｐｐｉｎｇＱＴＬｒｅｇｕｌａｔｉｎｇｍｏｒｐｈｏｐｈｙｓｉｏｌｏｇｉｃａｌｔｒａｉｔｓａｎｄｙｉｅｌｄ：Ｃａｓｅｓｔｕｄｉｅｓ，ｓｈｏｒｔ
ｃｏｍｉｎｇｓａｎｄｐｅｒｓｐｅｃｔｉｖｅｓｉｎｄｒｏｕｇｈｔｓｔｒｅｓｓｅｄｍａｉｚｅ［Ｊ］．ＡｎｎａｌｅｓＢｏｔａｂｉｃｉＦｅｎｎｉｃｉ，２００２，８９（７）：９４１９６３．
［３５］　严建兵，汤华，黄益勤，等．玉米和水稻重要性状ＱＴＬ的比较研究［Ｊ］．遗传学报，２００４，３１：１４０１１４０７．
［３６］　ＡｌｅｘａｎｄｒｏｖＮＮ，ＢｒｏｖｅｒＷ，ＦｒｅｉｄｉｎＳ，犲狋犪犾．ＩｎｓｉｇｈｔｉｎｔｏｃｏｒｎｇｅｎｅｓｄｅｒｉｖｅｄｆｒｏｍｌａｒｇｅｓｃａｌｅｃＤＮＡｓｅｑｕｅｎｃｉｎｇ［Ｊ］．Ｐｌａｎｔ
ＭｏｌｅｃｕｌａｒＢｉｏｌｏｇｙ，２００９，６９：１７９１９４．
［３７］　ＣａｒｏｌＳ，ＡｎｎｅＤ，ＤａｖｅＫ，犲狋犪犾．Ｓｅｑｕｅｎｃｉｎｇ，ｍａｐｐｉｎｇ，ａｎｄａｎａｌｙｓｉｓｏｆ２７，４５５ｍａｉｚｅｆｕｌｌｅｎｇｔｈｃＤＮＡｓ［Ｊ］．ＰｌｏｓＧｅｎｅｔｉｃｓ，
２００９，５（１１）：１０００７１３１０００７４１．
４８１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１２） Ｖｏｌ．２１，Ｎｏ．４
［３８］　ＷｏｎｇＪＣ，ＬａｍｂｅｒｔＲＪ，ＷｕｒｔｚｅｌＥＴ，犲狋犪犾．ＱＴＬａｎｄｃａｎｄｉｄａｔｅｇｅｎｅｓｐｈｙｔｏｅｎｅｓｙｎｔｈａｓｅａｎｄｚｅｔａｃａｒｏｔｅｎｅｄｅｓａｔｕｒａｓｅａｓｓｏｃｉ
ａｔｅｄｗｉｔｈｔｈｅａｃｃｕｍｕｌａｔｉｏｎｏｆｃａｒｏｔｅｎｏｉｄｓｉｎｍａｉｚｅ［Ｊ］．ＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，２００４，１０８（２）：３４９３５９．
［３９］　ＺｈａｎｇＰ，ＷａｎｇＹ，ＺｈａｎｇＪ，犲狋犪犾．ＡｍａｉｚｅＱＴＬｆｏｒｓｉｌｋｍａｙｓｉｎｌｅｖｅｌｓｃｏｎｔａｉｎｓｄｕｐｌｉｃａｔｅｄＭｙｂｈｏｍｏｌｏｇｏｕｓｇｅｎｅｓｗｈｉｃｈ
ｊｏｉｎｔｌｙｒｅｇｕｌａｔｅｆｌａｖｏｎｅｂｉｏｓｙｎｔｈｅｓｉｓ［Ｊ］．ＰｌａｎｔＭｏｌｅｃｕｌａｒＢｉｏｌｏｇｙ，２００３，５２（１）：１１５．
［４０］　肖景华，吴昌银，韩斌，等．中国水稻功能基因组研究进展［Ｊ］．中国科学Ｃ辑：生命科学，２００９，３９（１０）：９０９９２４．
犆狅狀狊狋狉狌犮狋犻狅狀狅犳犻狀狋犲犵狉犪狋犻狅狀犙犜犔犿犪狆犪狀犱犻犱犲狀狋犻犳犻犮犪狋犻狅狀狅犳犮犪狀犱犻犱犪狋犲犵犲狀犲狊犳狅狉狊狋犪狔犵狉犲犲狀犻狀犿犪犻狕犲
ＦＡＮＧＹｏｎｇｆｅｎｇ１，２，３，ＬＩＹｏｎｇｓｈｅｎｇ３，ＢＡＩＪｉａｎｇｐｉｎｇ１，２，３，ＭＵＰｉｎｇ３，ＭＥＮＧＹａｘｉｏｎｇ１，２，３，
ＺＨＡＮＧＪｉｎｌｉｎ４，ＷＡＮＧＨａｎｎｉｎｇ１，２，３，ＳＨＡＮＧＸｕｎｗｕ３
（１．ＧａｎｓｕＫｅｙＬａｂｏｒａｔｏｒｙｏｆＣｒｏｐＩｍｐｒｏｖｅｍｅｎｔａｎｄＧｅｒｍｐｌａｓｍＥｎｈａｎｃｅｍｅｎｔ，Ｌａｎｚｈｏｕ７３００７０，Ｃｈｉｎａ；
２．ＧａｎｓｕＰｒｏｖｉｎｃｉａｌＫｅｙＬａｂｏｒａｔｏｒｙｏｆＡｒｉｄｌａｎｄＣｒｏｐＳｃｉｅｎｃｅ，Ｌａｎｚｈｏｕ７３００７０，Ｃｈｉｎａ；
３．ＣｏｌｅｇｅｏｆＡｇｒｏｎｏｍｙ，ＧａｎｓｕＡｇｒｉｃｕｌｔｕｒａｌＵｎｉｖｅｒｓｉｔｙ，Ｌａｎｚｈｏｕ７３００７０，Ｃｈｉｎａ；
４．ＣｏｌｅｇｅｏｆＰａｓｔｏｒａｌＡｇｒｉｃｕｌｔｕｒｅＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，
ＬａｎｚｈｏｕＵｎｉｖｅｒｓｉｔｙ，Ｌａｎｚｈｏｕ７３００２０，Ｃｈｉｎａ）
犃犫狊狋狉犪犮狋：Ｌｅａｆｓｅｎｅｓｃｅｎｃｅｉｓａｃｈａｒａｃｔｅｒｏｆｐｌａｎｔｄｅｖｅｌｏｐｍｅｎｔ．Ｄｅｌａｙｅｄｌｅａｆｓｅｎｅｓｃｅｎｃｅｏｒｓｔａｙｇｒｅｅｎｃｈａｒａｃｔｅｒ
ｈａｓｒｅｃｅｉｖｅｄｃｏｎｓｉｄｅｒａｂｌｅａｔｔｅｎｔｉｏｎｂｅｃａｕｓｅｏｆｉｔｓｎｅｇａｔｉｖｅｉｍｐａｃｔｏｎｐｈｏｔｏｓｙｎｔｈｅｓｉｓａｎｄｉｔｓｒｏｌｅｉｎｎｕｔｒｉｅｎｔｒｅ
ｄｉｓｔｒｉｂｕｔｉｏｎｗｉｔｈｉｎｔｈｅｐｌａｎｔ．Ｉｎｔｈｅｐａｓｔｆｅｗｄｅｃａｄｅｓ，ａｗｅａｌｔｈｏｆＱＴＬｓ（ｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｕｓ）ｍａｐｐｉｎｇｄａｔａ
ｆｏｒｓｔａｙｇｒｅｅｎａｎｄｉｔｓｒｅｌａｔｅｄｔｒａｉｔｓｉｎｍａｉｚｅｈａｓｂｅｅｎｐｒｏｄｕｃｅｄｕｓｉｎｇｍｏｌｅｃｕｌａｒｍａｒｋｅｒａｐｐｒｏａｃｈｅｓ．Ｉｎｏｒｄｅｒｔｏ
ｕｎｌｏｃｋｔｈｅｆｕｌｐｏｔｅｎｔｉａｌｏｆｔｈｅｉｎｆｏｒｍａｔｉｏｎｃｏｎｔａｉｎｅｄｉｎｔｈｅｓｅｉｎｄｅｐｅｎｄｅｎｔｅｘｐｅｒｉｍｅｎｔｓ，１７３ｓｔａｙｇｒｅｅｎａｎｄｒｅ
ｌａｔｅｄＱＴＬｓｉｎｍａｉｚｅｗｅｒｅｃｏｌｅｃｔｅｄｆｒｏｍｍａｉｚｅｇｅｎｏｍｉｃｄａｔａｂａｓｅ（ｍａｉｚｅＧＤＢ），ａｎｄｗｅｒｅｃｏｍｐｉｌｅｄｔｏｃｏｎｓｔｒｕｃｔ
ｔｈｅｉｎｔｅｇｒａｔｉｏｎＱＴＬｍａｐｆｏｒｓｔａｙｇｒｅｅｎｉｎｍａｉｚｅ（犣犲犪犿犪狔狊）．Ｔｈｅｈｉｇｈｓｃａｌｅｇｅｎｅｔｉｃｌｉｎｋａｇｅｍａｐｏｆｍａｉｚｅ
ＩＢＭ２２００８ＮｅｉｇｈｂｏｒｓｗａｓｕｓｅｄａｓｒｅｆｅｒｅｎｃｅａｎｄｔｈｅｍａｐｗｅｒｅｃｏｎｓｔｒｕｃｔｅｄｗｉｔｈＢｉｏＭｅｒｃａｔｏｒ２．１ｓｏｆｔｗａｒｅｂａｓｅｄ
ｏｎｔｈｅｍｅｔｈｏｄｏｆｍｅｔａａｎａｌｙｓｉｓ．ＦｉｖｅｃｏｎｓｅｎｓｕｓＱＴＬｆｏｒｍａｉｚｅｓｔａｙｇｒｅｅｎｗｅｒｅｉｄｅｎｔｉｆｉｅｄｏｎｍａｉｚｅｃｈｒｏｍｏ
ｓｏｍｅｓ１，４，５ａｎｄ９，ｒｅｓｐｅｃｔｉｖｅｌｙ．ＡｆｔｅｒｐｒｏｊｅｃｔｔｈｅｃｏｎｓｅｎｓｕｓＱＴＬｒｅｇｉｏｎｏｎｍａｉｚｅｐｈｙｓｉｃａｌｍａｐｏｆＢ７３Ｒｅｆ
Ｇｅｎ＿ｖ２，ａｔｏｔａｌｏｆ１４４５ｃａｎｄｉｄａｔｅｇｅｎｅｓｉｎｔｈｅｃｏｎｓｅｎｓｕｓＱＴＬｒｅｇｉｏｎｗｅｒｅｄｉｓｃｏｖｅｒｅｄａｎｄｄｏｗｎｌｏａｄｅｄｆｒｏｍ
ＰｌａｎｔＧＤＢ（ｈｔｔｐ：／／ｗｗｗ．ｐｌａｎｔｇｄｂ．ｏｒｇ／）．Ｍｅａｎｗｈｉｌｅ，ａｃｃｏｒｄｉｎｇｔｏｔｈｅｒｅｓｕｌｔｓｏｆｇｅｎｅｏｎｔｏｌｏｇｙａｎａｌｙｓｉｓ，ｍｏｓｔ
ｃａｎｄｉｄａｔｅｇｅｎｅｓｍｉｇｈｔｔａｋｅｐａｒｔｉｎｔｈｅｃｅｌｄｅｖｅｌｏｐｍｅｎｔ，ｃａｔａｌｙｔｉｃｐｒｏｃｅｓｓ，ｅｎｚｙｍａｔｉｃａｃｔｉｖｉｔｙａｎｄｏｔｈｅｒｂｉｏｌｏｇｉ
ｃａｌｐｒｏｇｒｅｓｓｅｓ．Ｔｈｅｓｅｑｕｅｎｃｅｓｏｆｅｉｇｈｔｓｔａｙｇｒｅｅｎｃａｎｄｉｄａｔｅｇｅｎｅｓｗｅｒｅｉｄｅｎｔｉｆｉｅｄｂｙｃｏｍｐａｒａｔｉｖｅｌｙａｎａｌｙｚｉｎｇ
ｔｈｅｓｅｑｕｅｎｃｅｓｏｆｓｔａｙｇｒｅｅｎｇｅｎｅｓｅｑｕｅｎｃｅｆｒｏｍｏｔｈｅｒｃｒｏｐｓ．Ｕｓｉｎｇａｓｙｎｔｏｎｉｃｃｏｎｓｅｒｖａｔｉｏｎａｐｐｒｏａｃｈ，ｗｈｉｃｈｉｓ
ｕｓｅｄｔｏｃｏｍｐａｒａｔｉｖｅｌｙａｎａｌｙｚｅｔｈｅｒｉｃｅｇｅｎｅｔｉｃｍａｐａｎｄｍａｉｚｅｐｈｙｓｉｃａｌｍａｐ，ｏｎｅｒｉｃｅｇｅｎｅ（狊犵狉）ｔｈａｔｉｓｃｏｄｉｎｇ
ｆｏｒｔｈｅｓｅｎｅｓｃｅｎｃｅｉｎｄｕｃｉｂｌｅｃｈｌｏｒｏｐｌａｓｔｓｔａｙｇｒｅｅｎｐｒｏｔｅｉｎｗａｓｐｒｏｊｅｃｔｅｄｏｎｔｏｔｈｅｐｈｙｓｉｃａｌｍａｐｏｆｍａｉｚｅ
Ｂ７３ＲｅｆＧｅｎ＿ｖ２．Ｔｈｅｒｅｆｏｒｅ，ｔｈｅｒｅｓｕｌｔｓｓｕｇｇｅｓｔｅｄｔｈａｔｔｈｉｓａｐｐｒｏａｃｈｐｒｏｖｉｄｅｄａｎｕｓｅｆｕｌｔｏｏｌｆｏｒｉｄｅｎｔｉｆｙｉｎｇＱＴＬ
ｃｏｎｆｅｒｒｉｎｇｓｔａｙｇｒｅｅｎａｎｄｏｔｈｅｒｃａｎｄｉｄａｔｅｇｅｎｅｓｉｎｍａｉｚｅ．
犓犲狔狑狅狉犱狊：犣犲犪犿犪狔狊；ｓｔａｙｇｒｅｅｎ；ｑｕａｎｔｉｔａｔｉｖｅｔｒａｉｔｌｏｃｕｓ（ＱＴＬ）；ｉｎｔｅｇｒａｔｉｏｎｍａｐ；ｍｅｔａａｎａｌｙｓｉｓ
５８１第２１卷第４期 草业学报２０１２年