高志勇 袁庭攀 趙亞蘭 尹三霞

摘要：介紹了玉米（Zea mays L.）基因組學(xué)的研究發(fā)展歷程，概述了玉米基因組學(xué)的研究進(jìn)展。重點(diǎn)論述了玉米基因組的結(jié)構(gòu)特征，玉米功能基因組學(xué)的發(fā)展現(xiàn)狀，玉米重要農(nóng)藝性狀的基因克隆，包括株高和株型相關(guān)基因的克隆，產(chǎn)量相關(guān)性狀基因的克隆以及品質(zhì)性狀相關(guān)基因的克隆，并展望了玉米基因組學(xué)的發(fā)展前景。

關(guān)鍵詞：玉米（Zea mays L.）；基因組學(xué)；農(nóng)藝性狀；基因克隆

中圖分類號：S513；Q785 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2016）05-1094-05

DOI：10.14088/j.cnki.issn0439-8114.2016.05.002

Research Progress on Maize Genomics and Gene Cloning of Agronomic Traits

GAO Zhi-yong1，2，YUAN Ting-pan1，ZHAO Ya-lan1，YIN San-xia1

（1. School of Chemistry and Life Science，Weinan Normal University，Weinan 714099，Shaanxi， China；

2. Key Laboratory of River Wetland Ecology and Environment in Shaanxi Province，Weinan 714099，Shaanxi， China）

Abstract： The paper introduces the development process of maize genomics， and summarizes the research progress of maize genomics. It focuses on the structural characteristics of the maize genome， functional genomics， and gene cloning of important agronomic traits， including plant height and plant type， yield， and quality traits. The future work was prospected on the development of maize genomics in the end.

Key words： maize（Zea mays L.）； genomics； agronomic traits； gene cloning

大約在1萬年前，在美洲中部墨西哥馴化了當(dāng)?shù)氐念愂袷?，形成了玉米（Zea mays L.），從此，玉米一直在種植和培育。在歷史的長河中，玉米基因組經(jīng)過了多輪的復(fù)制。大約在7 000萬年前，發(fā)生了古多倍體祖先的全基因組復(fù)制；5 000萬～12 000萬年前，進(jìn)行了額外全基因組復(fù)制等。通過這些事件，玉米和它的近親高粱（Sorghum bicolor）分離開來。玉米的基因組有10條染色體，在進(jìn)化中經(jīng)歷了動態(tài)變化，使得玉米基因組在結(jié)構(gòu)上顯示出復(fù)雜性。在過去的大約300萬年中，發(fā)生了長末端重復(fù)反轉(zhuǎn)錄轉(zhuǎn)座子（LTR retrotransposons）的增殖，使得玉米基因組發(fā)生了大規(guī)模的急劇擴(kuò)增，基因組的大小達(dá)到了2.3×109。目前，作為一種重要的單子葉模式植物，尤其是由于其功能基因組學(xué)的發(fā)展和進(jìn)步，玉米用在了許多重要的基礎(chǔ)研究上，主要的研究領(lǐng)域有基因和染色體的物理連鎖、基因的遺傳和功能、發(fā)育和細(xì)胞生物學(xué)、生物信息學(xué)的交叉、端粒的特征、基因的轉(zhuǎn)座、核仁的起源、數(shù)量遺傳和育種以及表觀遺傳學(xué)研究等。

1 玉米基因組的結(jié)構(gòu)特征

2005年美國國家科學(xué)基金會、能源部、農(nóng)業(yè)部資助玉米全基因組的測序工作。來自于美國的亞利桑那大學(xué)和圣路易斯華盛頓大學(xué)等科研機(jī)構(gòu)的150名科研人員承接了這項(xiàng)工作。研究小組選取了B73自交系為測序材料，工作人員經(jīng)過了4年多的艱苦努力，終于完成了測序工作，研究成果發(fā)表在了2009年11月份的《科學(xué)》期刊上[1]。測序工作中，科研人員利用細(xì)菌人工染色體（Bacterial artificial chromosome，BAC）（n=16 848）和來源于完整的遺傳和物理圖譜的Fosmid克?。╪=63），進(jìn)行了全基因組的測序，并且運(yùn)用了光學(xué)圖譜進(jìn)行了強(qiáng)度比。其中通過大約4～6倍的鳥槍法（Shotgun）對單克隆進(jìn)行了測序，對一些獨(dú)特的區(qū)域又做了人工與自動完善，獲得了首版的B73參考基因組（B73 RefGen_v1）。B73 RefGen_v1揭示了玉米基因組的大小達(dá)到了2.3×109，并預(yù)測了32 000多個基因，其中99.8%的基因定位在了染色體上。Messing等[2]早年通過對玉米的全基因組中的大量插入克隆片段末端的測序，測序了基因組中的6.2 kb的片段，共累積了大約1/8的全基因組序列，其中對大量的基因和轉(zhuǎn)座子進(jìn)行了描述?？偟睦鄯e的序列有307 Mb，其中重復(fù)序列占了58%，基因區(qū)域占了7.5%，預(yù)測了大約59 000個基因。

在玉米的基因組中，由于存在高密度重復(fù)序列和復(fù)雜的結(jié)構(gòu)，并且還有活性高的Mutator 轉(zhuǎn)座子，所以要對基因組進(jìn)行精確的預(yù)測，會存在很大的困難。這個基因組中有幾百個轉(zhuǎn)座子基因家族，總共占了接近85%，在基因組中非均勻地散布，它們的轉(zhuǎn)座，造成了大量的基因復(fù)制和捕獲，該基因組中有接近85%是由非均勻分散于整個基因組中的幾百個轉(zhuǎn)座子家族組成。這些轉(zhuǎn)座元件影響了數(shù)以千計(jì)的基因片段的捕獲和復(fù)制，改變了著絲點(diǎn)的組成、大小及位置。

在玉米的基因組中，轉(zhuǎn)座因子（Transposable element，TE）含量非常豐富，B73 RefGen_v1中顯示轉(zhuǎn)座子占了85%，其中有豐富的Ⅰ類RNA轉(zhuǎn)座子和活躍的Ⅱ類DNA 轉(zhuǎn)座子[3]。在玉米中，最先發(fā)現(xiàn)了轉(zhuǎn)座子，并且 Spm/En系統(tǒng)、PIF/Harbinger、Hat（Ac轉(zhuǎn)座子）、Mutator 超基因家族、MITE家族最初也都是在玉米中發(fā)現(xiàn)的。

在玉米中還存在著長末端重復(fù)反轉(zhuǎn)錄轉(zhuǎn)座子，并且也非常豐富。大約82%的轉(zhuǎn)座子已經(jīng)得到了鑒定。在各類轉(zhuǎn)座子家族中，最為復(fù)雜的是Mutator 超家族，其元件的大小和序列有著極大的差異，其中包含有262個MULEs（Mutator-like elements that contain gene fragments），MULEs 含有226個核基因的片段。通過對40 000個Mutator 活性區(qū)域的非冗余轉(zhuǎn)座子插入位點(diǎn)進(jìn)行擴(kuò)增和測序，對應(yīng)地繪制到了B73 RefGen_v1中，試驗(yàn)結(jié)果顯示，在減數(shù)分裂中，Mutator 插入位點(diǎn)重組率最高，許多轉(zhuǎn)座子都插入在含豐富基因及重組活性高的染色體的末端。

在植物、動物和真菌中，存在著一類稱作Helitrons 的轉(zhuǎn)座子。Helitrons 在玉米中含量豐富，轉(zhuǎn)座活性高，轉(zhuǎn)座方式采用滾環(huán)模式。在玉米中，有8個Helitrons 家族，大約總共有20 000個拷貝，在獲取基因片段時，其表現(xiàn)異?；钴S。在早先的動植物基因組研究中，顯示Helitrons大部分存在于基因富集區(qū)[4]。在B73 RefGen_v1中，長末端重復(fù)反轉(zhuǎn)錄轉(zhuǎn)座子的含量超過了全基因組的75%，它們在基因組中分散存在，它們的406個基因家族中，大部分的拷貝數(shù)小于10個。LTR retrotransposons 具有家族特異性，在染色體上分布也不均勻，像類Copia 因子（Copia-like elements）常常在基因豐富的常染色質(zhì)區(qū)域大量存在，而類Gypsy 因子（Gypsy-like elements）卻大量存在于基因稀疏的異染色質(zhì)區(qū)域，研究顯示長末端重復(fù)反轉(zhuǎn)錄轉(zhuǎn)座子含有的核基因片段大于180個。

隨著功能基因組學(xué)發(fā)展，產(chǎn)生了海量數(shù)據(jù)，也為研究提供了很好的生物信息學(xué)的平臺。玉米遺傳和基因組信息資源MaizeGDB（http：//www.maizegdb.org/）提供了玉米相關(guān)功能基因的注釋和信息。美國谷類作物數(shù)據(jù)庫（http：//harvest.ucr.edu/）收集了大麥、短柄草屬、柑橘、咖啡、豇豆、大豆、水稻、小麥等作物的表達(dá)譜數(shù)據(jù)及相關(guān)的一些分子信息。而另一谷類作物數(shù)據(jù)庫Gramene Database（http：//www.gramene.org/）是一個以各種作物基因組信息為主的數(shù)據(jù)庫，同時具備高級的各基因組間的分析功能。英國UK CROPNET（http：//ukcrop.net/）是英國農(nóng)作物生物信息學(xué)網(wǎng)絡(luò)數(shù)據(jù)庫，擁有很多其自己開發(fā)的數(shù)據(jù)庫和分析軟件，同時也收集了相關(guān)的一些文獻(xiàn)和美國該領(lǐng)域的一些信息。水稻（http：//ars-genome.cornell.edu/rice/）、小麥（http：//www.tigr.org/tdb/e2k1/tae1/）和Panzea（http：//panzea.org/database）等專一作物的數(shù)據(jù)庫，由美國TIGR institute維護(hù)的小麥基因組數(shù)據(jù)，提供小麥的基因組及基因注釋，同時，還提供其它谷類作物的同源基因數(shù)據(jù)，如玉米、大麥、高粱、水稻等主要農(nóng)作物，為玉米功能基因組學(xué)研究提供大量不同數(shù)據(jù)和信息。

2 玉米功能基因組學(xué)研究進(jìn)展

在玉米的研究史上，2009年的玉米自交系B73的測序完成可以說是一個大的飛躍，為玉米基因組學(xué)的研究提供了巨大的數(shù)據(jù)支撐。而后隨著分子生物學(xué)的迅猛發(fā)展，測序技術(shù)的進(jìn)一步完善，B73 RefGen_v2隨之問世，使對玉米基因組的認(rèn)識更加全面和精細(xì)，給玉米的分子生物學(xué)研究帶來了更大的便利。最近，在由冷泉港實(shí)驗(yàn)室和深圳華大基因研究院等單位主導(dǎo)完成的玉米遺傳變異研究論文發(fā)表在Nature Genetics[5]上。Chia等[5]對105個野生和栽培玉米品種進(jìn)行了測序和分析，成功構(gòu)建了第二代玉米單體型圖譜（簡稱Maize HapMap 2），并對玉米的遺傳多樣性進(jìn)行了全面分析。在研究中，他們建立了精密的群體遺傳學(xué)評分模型（Population-genetics scoring model），共鑒定了5 500多萬個遺傳分子標(biāo)記，同時發(fā)現(xiàn)染色體結(jié)構(gòu)變異（SV）在玉米野生種和栽培種的基因組中是普遍存在的，并且推測這些結(jié)構(gòu)變異與很多重要的農(nóng)藝性狀具有密切的關(guān)聯(lián)。

此外，在研究影響玉米基因組大小的主要因素時，研究人員發(fā)現(xiàn)染色體結(jié)（一種密集的DNA結(jié)構(gòu)，Chromosomal knobs）的存在或缺失造成了“種內(nèi)”玉米基因組大小出現(xiàn)很大的差異；而在“種間”進(jìn)行比較時，他們發(fā)現(xiàn)玉米基因組大小的變化則主要與大量的轉(zhuǎn)座子有關(guān)。

考古學(xué)及遺傳學(xué)證實(shí)玉米大約馴化自1萬多年前，在馴化過程中，玉米經(jīng)歷了一次特殊的表型轉(zhuǎn)變，這使其更加容易滿足人類的需求。通過對75個玉米株系（包括野生玉米、美洲各地的傳統(tǒng)品種和現(xiàn)代改良玉米品系）進(jìn)行全基因組重測序，并對玉米馴化進(jìn)行了全面的評估分析。研究發(fā)現(xiàn)玉米在經(jīng)過人工馴化之后又產(chǎn)生了新的遺傳多樣性，并推測這很有可能是由于野生近緣物種的基因滲入所導(dǎo)致的。研究人員發(fā)現(xiàn)了數(shù)百個具有強(qiáng)烈選擇信號的基因，并推測這些基因可能在玉米的馴化過程中發(fā)揮著重要的作用[6]。此外，研究數(shù)據(jù)還表明在玉米的馴化過程中，千百年前古代農(nóng)民應(yīng)用的人工馴化方法似乎比現(xiàn)代育種學(xué)家所使用的方法發(fā)揮了更大的作用。

通過馴化，人類已經(jīng)大大地改變了大芻草形態(tài)，成為今天可以辨認(rèn)的玉米。這個系統(tǒng)作為研究復(fù)雜性狀適應(yīng)、基因組進(jìn)化以及其遺傳學(xué)和進(jìn)化的模式。為了研究馴化如何重新形成玉米幼苗的轉(zhuǎn)錄組，分析38個不同的基因型玉米和24類蜀黍的基因型的18 242個基因表達(dá)譜，玉米與大芻草相比，超過600個基因具有顯著不同的表達(dá)水平。此外，超過1 100個基因表現(xiàn)出顯著共表達(dá)改變。這項(xiàng)研究不僅證明了隨著馴化玉米轉(zhuǎn)錄組發(fā)生改變，鑒定了幾個負(fù)責(zé)玉米進(jìn)化的基因，而且明確了基因表達(dá)結(jié)合群體遺傳學(xué)分析可以獲得基因互補(bǔ)的信息[7]。

3 玉米重要農(nóng)藝性狀基因的克隆

目前，許多模式植物已經(jīng)測序，研究也更加深入，基因功能的研究方法也更加完善，玉米的功能基因組學(xué)研究也成為了分子生物學(xué)研究的熱點(diǎn)領(lǐng)域之一。定位克隆（Positional cloning）在早期的克隆基因研究中應(yīng)用較多。相對于其他一些模式植物，玉米的基因組較大，重復(fù)序列多，重要基因的定位克隆令人望而卻步。隨著玉米全基因組序列的發(fā)表，DNA測序技術(shù)進(jìn)步，以及玉米全基因組關(guān)聯(lián)分析和比較基因組學(xué)技術(shù)發(fā)展等，應(yīng)用定位克隆技術(shù)也成功地克隆了玉米的重要關(guān)鍵基因和QTL。整合圖位克隆，比較基因組和關(guān)聯(lián)分析等多種方法來加速玉米株高主效QTL的克隆，即先使用高密度的遺傳圖譜迅速將玉米主效株高QTLs轉(zhuǎn)變?yōu)閱蝹€孟德爾因子，將其表型效應(yīng)放大，然后再逐一地精細(xì)定位；接著利用比較基因組將不同作物之間各種不同的功能基因組平臺串聯(lián)起來進(jìn)行基因候選；再使用定點(diǎn)的關(guān)聯(lián)分析、基因的跨物種轉(zhuǎn)化或者不同作物的株高突變體來驗(yàn)證基因的功能。

3.1 株高和株型相關(guān)基因的克隆

莖稈是玉米株型的重要組成部分，是影響產(chǎn)量的重要因素。在玉米栽培和馴化過程中，其株型發(fā)生了兩次革命性的變化：即由多分蘗大芻草到獨(dú)稈和由高稈玉米農(nóng)家種到半矮稈的玉米自交系的轉(zhuǎn)化。兩者變革使玉米的株型緊湊、密植不倒伏成為可能，為玉米的高產(chǎn)奠定了基礎(chǔ)。

目前，已知控制株高的基因在單子葉和雙子葉植物高度保守。玉米dwarf8（d8）突變體莖稈極端矮小，葉片縮短，花期變長。它編碼一個DELLA蛋白，是一個顯性對赤霉素不敏感的矮稈基因。dwarf3（d3）編碼一個細(xì)胞色素P450（CYP88A）蛋白，作用于赤霉素合成途徑[8]。玉米br2突變體莖稈節(jié)間縮短，二葉片、雄蕊的大小正常，br2基因控制生長素在莖稈中的極性運(yùn)輸。玉米na1編碼一個與擬南芥DET2同源的蛋白，調(diào)節(jié)油菜素內(nèi)酯生物合成，進(jìn)而影響植株的高低。

在當(dāng)前的研究中，玉米株高等重要農(nóng)藝性狀QTL的克隆進(jìn)展緩慢，主要原因表現(xiàn)在：大多數(shù)玉米株高等重要農(nóng)藝性狀受多個遺傳效應(yīng)低的QTL控制；表型自然變異的功能性位點(diǎn)多位于調(diào)控區(qū)，過于隱蔽，不易定位。在B73 X Mo17的群體中共定位了6個控制株高的QTL，解釋了73%的表型變異，加性效應(yīng)在6.1 cM到9.1 cM之間。玉米開花相關(guān)基因vgt1（Vegetative to generative transition 1）功能性位點(diǎn)在該基因上游70 kb啟動子區(qū)[9]。控制玉米分蘗數(shù)的tb1開放讀碼框上游60～70 kb位置上插入了Hopscotch轉(zhuǎn)座子，使得玉米tb1基因表達(dá)上調(diào)，從而玉米多分蘗變成獨(dú)稈。

株型緊湊的玉米品種莖葉夾角小，上部葉片趨于直立。透光性好，宜于密植。最近，有人利用GWAS全基因組關(guān)聯(lián)分析，從玉米基因組中發(fā)現(xiàn)了160萬個遺傳位點(diǎn)。在此基礎(chǔ)上鑒定了與莖葉夾角性狀相關(guān)的基因。玉米品種之間莖葉夾角差異最大可達(dá)80。但是，發(fā)現(xiàn)單基因的最大效應(yīng)僅能達(dá)到1.5。因此莖葉夾角的性狀是通過多個基因產(chǎn)生微小效應(yīng)疊加的結(jié)果[10]。

3.2 產(chǎn)量相關(guān)性狀的基因克隆

玉米產(chǎn)量性狀也是復(fù)雜的數(shù)量性狀，受多基因和環(huán)境變化的影響。例如，用3×87-1的441個永久F2群體，定位了27個包括穗行數(shù)、百粒重等產(chǎn)量性狀的QTL。對71個歐洲玉米自交系掃描了9個多態(tài)性位點(diǎn)，證明其中6個與玉米開花期相關(guān)[11]。最近，在玉米第四染色體上克隆fea2基因，控制穗行數(shù)，并且存在一個產(chǎn)生基因表達(dá)量差異的功能性變異，但其位置仍然未知[12]。

落粒性也是一個重要的產(chǎn)量性狀。利用玉米祖先大芻草和玉米高代回交，含886個體的大群體迅速將兩個玉米主效落粒QTL精細(xì)定位[13]。該基因克隆策略是利用圖位克隆、關(guān)聯(lián)分析和比較基因組學(xué)手段結(jié)合馴化綜合征在禾本科植物類似分子遺傳學(xué)基礎(chǔ)的原理，首先同時克隆了高粱一個質(zhì)量基因Sh1、兩個玉米落?；?、一個水稻落粒QTL。高粱Sh1落?；蛐蛄写嬖谌齻€功能變異，分別是調(diào)控Sh1表達(dá)的單倍型，Sh1失去第二和第三外顯子的大片段缺失，一個是調(diào)控Sh1第四外顯子被異常剪切的剪切位點(diǎn)變異。通過比較基因組學(xué)定位獲得四個同源區(qū)段，兩個玉米落粒QTL、一個水稻落粒QTL和一個谷子落粒QTL分別位于這4個區(qū)段內(nèi)。因而，實(shí)現(xiàn)了QTL跨物種的基因克隆。

同時，玉米抗倒伏性狀是玉米產(chǎn)量相關(guān)性狀的關(guān)鍵因素，與莖稈組織細(xì)胞壁結(jié)構(gòu)，纖維素組成和機(jī)械強(qiáng)度密切相關(guān)，作物倒伏或折斷是玉米高產(chǎn)穩(wěn)產(chǎn)和機(jī)械化收獲的重要生物學(xué)基礎(chǔ)。作物抗倒伏性是由作物倒伏基因或者稱作脆稈基因所控制的。水稻中的脆稈基因有較為詳細(xì)的報導(dǎo)，但是在玉米中的研究甚少。已知玉米中有6個脆稈基因，分別是BK1、BK2、BK3、btn1、bk-x、bk2，而這些報道只是簡單遺傳和表型分析。利用基因同源性，推測bk2可能編碼類似糖基錨定蛋白，調(diào)控次生細(xì)胞壁纖維素合成，進(jìn)而影響莖稈的機(jī)械強(qiáng)度[14]。因此，玉米的脆稈基因研究有待于功能基因組等相關(guān)技術(shù)的綜合運(yùn)用。

3.3 品質(zhì)性狀相關(guān)基因的克隆

玉米子粒的形成過程涉及一系列基因時空表達(dá)的改變和調(diào)節(jié)。利用4～6 d和7～23 d的胚乳構(gòu)建cDNA文庫。分析其中Cotigs和Singletons，排除貯藏蛋白基因外，至少有5 000個基因發(fā)生了變化[15]。DNA芯片分析表明，在玉米幼穗4個發(fā)育階段，有2 794個基因表達(dá)發(fā)生變化[16]。針對玉米子粒動態(tài)形成過程，分別對授粉后不同灌漿階段粒重、子粒灌漿速率、淀粉積累相關(guān)酶活性和可溶性碳水化合物濃度進(jìn)行QTL定位分析，表明粒重QTL定位分析結(jié)果受基礎(chǔ)材料、定位群體類型和大小、子粒灌漿時期、遺傳圖譜標(biāo)記飽和度、定位方法和環(huán)境條件等多種因素的單獨(dú)和綜合影響。通過構(gòu)建NILs精細(xì)定位和圖位克隆主要目標(biāo)QTL，具有更為重要的理論意義。

子粒的氨基酸組成是衡量蛋白質(zhì)品質(zhì)的重要指標(biāo)，尤其是賴氨酸和色氨酸的含量。Opaque-2基因突變體是高賴氨酸突變體，賴氨酸含量達(dá)4%左右[17]，這是由于Opaque-2子粒中zein蛋白含量大幅度下降，而非zein蛋白含量上升。但是，其子粒是軟質(zhì)胚乳，產(chǎn)量相對較低，因而限制了其應(yīng)用。后來由于Opaque-2修飾因子發(fā)現(xiàn)和應(yīng)用，可使玉米子粒賴氨酸含量增加1倍。玉米Opaque-2基因控制的子粒，在普通玉米中賴氨酸含量為0.20%，色氨酸為0.06%，而高賴氨酸玉米中分別達(dá)到0.48%和0.13%，是普通玉米的兩倍多。并且高賴氨酸玉米子粒中的組氨酸、精氨酸、甘氨酸、蛋氨酸、天門冬氨酸等的含量略高于普通玉米。但是，在實(shí)際應(yīng)用過程中如何將高賴氨酸的基因和軟質(zhì)胚乳修飾的多個數(shù)量性狀基因?qū)牍歉捎N自交系當(dāng)中仍然是巨大的挑戰(zhàn)。

4 展望

隨著玉米研究的更加深入，EMS誘變和Mutator 誘變突變體庫的建立，使得更多的玉米質(zhì)量性狀的突變體被發(fā)現(xiàn)，定位克隆技術(shù)隨之發(fā)展更為迅速。目前，新的定位克隆技術(shù)在克隆基因上更為有效，更加完善了玉米基因組序列，玉米的功能基因組學(xué)研究將進(jìn)入一個嶄新的發(fā)展時期。

參考文獻(xiàn)：

[1] SCHNABLE P S， WARE D， FULTON R S， et al. The B73 maize genome： Complexity， diversity， and dynamics[J]. Science， 2009， 326（5956）： 1112-1125.

[2] MESSING J， BHARTI A K， KARLOWSKI W M， et al. Sequence composition and genome organization of maize[J]. Proceedings of the National Academy of Sciences of the United States of America， 2004， 101（40）： 14349-14354.

[3] SANMIGUEL P， TIKHONOV A， JIN Y K， et al. Nested retrotransposons in the intergenic regions of the maize genome [J]. Science， 1996， 274（5288）： 765-768.

[4] LAI J， LI Y， MESSING J， et al. Gene movement by Helitron transposons contributes to the haplotype variability of maize [J]. Proceedings of the National Academy of Sciences of the United States of America， 2005， 102（25）： 9068-9073.

[5] CHIA J M， SONG C， BRADBURY P J， et al. Maize HapMap2 identifies extant variation from a genome in flux[J]. Nature Genetics， 2012， 44（7）： 803-807.

[6] HUFFORD M B， XU X， VAN HEERWAARDEN J， et al. Comparative population genomics of maize domestication and improvement [J]. Nature Genetics，2012，44（7）：808-811.

[7] SWANSON-WAGNER R， BRISKINE R， SCHAEFER R，et al. Reshaping of the maize transcriptome by domestication[J]. Proceedings of the National Academy of Sciences，2012，109（29）： 11878-11883.

[8] HELLIWELL C A， SULLIVAN J A， MOULD R M， et al. A plastid envelope location of Arabidopsis ent-kaurene oxidase links the plastid and endoplasmic reticulum steps of the gibberellin biosynthesis pathway[J]. The Plant Journal，2001，28（2）： 201-218.

[9] SALVI S， TUBEROSA R， CHIAPPARINO E， et al. Toward positional cloning of Vgt1， a QTL controlling the transition from the vegetative to the reproductive phase in maize[J]. Plant Molecular Biology， 2002， 48（5）： 601-613.

[10] TIAN F， BRADBURY P J， BROWN P J， et al. Genome-wide association study of leaf architecture in the maize nested association mapping population[J]. Nature Genetics，2011， 43（2）： 159-162.

[11] ANDERSEN J R， SCHRAG T， MELCHINGER A E， et al. Validation of Dwarf8 polymorphisms associated with flowering time in elite European inbred lines of maize （Zea mays L.） [J]. Theoretical and Applied Genetics， 2005， 111（2）： 206-217.

[12] BOMMERT P， NAGASAWA N S， JACKSON D. Quantitative variation in maize kernel row number is controlled by the FASCIATED EAR2 locus[J]. Nature Genetics， 2013， 45（3）： 334-337.

[13] LIN Z， LI X， SHANNON L M， et al. Parallel domestication of the Shattering1 genes in cereals[J]. Nature Genetics， 2012， 44（6）：720-724.

[14] CHING A， DHUGGA K S， APPENZELLER L， et al. Brittle stalk 2 encodes a putative glycosyl phosphatidylinositol-anchored protein that affects mechanical strength of maize tissues by altering the composition and structure of secondary cell walls[J]. Planta，2006，224（5）：1174-1184.

[15] LAI J， DEY N， KIM CS， et al. Characterization of the maize endosperm transcriptome and its comparison to the rice genome[J]. Genome Research，2004，14（10a）：1932-1937.

[16] ZHU Y， FU J， ZHANG J， et al. Genome-wide analysis of gene expression profiles during ear development of maize[J]. Plant Molecular biology，2009，70（1-2）：63-77.

[17] MERTZ E T， BATES L S， NELSON O E. Mutant gene that changes protein composition and increases lysine content of maize endosperm[J]. Science，1964，145（3629）：279-280.