玉米葉部性狀的QTL定位與候選基因分析

收稿日期：2024-01-17

基金項目：國家自然科學(xué)基金項目（32060460、32160451）；貴州省科技支撐計劃項目[黔科合支撐（2022）重點029]；貴州省山地農(nóng)業(yè)關(guān)鍵核心技術(shù)攻關(guān)項目（GZNYGJHX-2023012）；貴州省育種科研基礎(chǔ)平臺創(chuàng)新能力建設(shè)項目[黔科合服企（2022）014]；黔農(nóng)科院國基后補助項目[（2021）16號、（2022）02號、（2022） 09號]；貴州喀斯特山區(qū)重要作物生物育種平臺建設(shè)項目[黔科合中引地（2023）033]；貴州省育種科研基礎(chǔ)平臺創(chuàng)新能力建設(shè)項目[黔科合服企[（2022）014]；貴州省科研機構(gòu)創(chuàng)新能力建設(shè)項目[（2022）007]

作者簡介：郭 爽（1997-），女，貴州晴隆人，碩士研究生，研究方向為作物遺傳育種。（Tel）18985987996；（E-mail）18985987996@163.com

通訊作者：吳 迅， （Tel）18798004601;（E-mail）wuxunyong@126.com;陳澤輝，（Tel）13985441082;（E-mail）chenzh907@sina.com

摘要： 葉片在玉米生長過程中發(fā)揮著重要作用，它能夠有效地進行光合作用，為玉米提供營養(yǎng)物質(zhì)，通過影響耐密性等影響產(chǎn)量提升。本研究選用QR273和T32為親本，構(gòu)建150份F2、F2∶3家系材料，結(jié)合基因型和不同環(huán)境中葉部性狀的表型評價數(shù)據(jù)，利用完備區(qū)間作圖法進行數(shù)量性狀座位（QTL）定位。結(jié)果發(fā)現(xiàn)，2個環(huán)境下共檢測到85個葉部性狀相關(guān)QTL，其中有12個全株葉片數(shù)相關(guān)QTL、14個穗上葉片數(shù)相關(guān)QTL、22個葉長相關(guān)QTL、17個葉寬相關(guān)QTL、20個葉夾角相關(guān)QTL。結(jié)合公共數(shù)據(jù)庫和生物信息學(xué)分析方法共篩選出7個候選基因。其中Zm00001d013612編碼微管蛋白，參與調(diào)控細胞骨架結(jié)構(gòu)組成；Zm00001d053543參與油菜素甾醇介導(dǎo)的信號通路；Zm00001d031291編碼的蛋白質(zhì)具有組蛋白乙?；δ埽籞m00001d031292參與富含羥脯氨酸糖蛋白家族基因表達的調(diào)控；Zm00001d031296調(diào)控鉀離子跨膜轉(zhuǎn)運蛋白活性；Zm00001d031300、Zm00001d031303參與碳水化合物代謝過程。蛋白質(zhì)功能分析結(jié)果表明，這7個候選基因均參與細胞分化，與植物的生長發(fā)育息息相關(guān)。本研究結(jié)果將為深度揭示玉米葉部性狀變異的遺傳基礎(chǔ)提供更豐富的理論支持。

關(guān)鍵詞： 玉米；葉部性狀；數(shù)量性狀座位（QTL）；候選基因

中圖分類號： S513"" 文獻標識碼： A"" 文章編號： 1000-4440（2024）10-1777-10

QTL mapping and candidate gene analysis of leaf related traits in maize

GUO Shuang^1，2， NIE Lei^1，2， HE Yue^1，2， WANG Dong^1，2， TU Liang²， LIU Pengfei²， JIANG Yulin^2，3，GUO Xiangyang²， WANG Angui²， ZHU Yunfang²， WU Xun^2，3， CHEN Zehui²

（1.Agricultural College of Guizhou University， Guiyang 550025， China;2.Institute of Upland Food Crops， Guizhou Academy of Agricultural Sciences， Guiyang 550006， China;3.Key Laboratory of Crop Genetic Resources and Germplasm Innovation in Karst Region， Ministry of Agriculture and Rural Affairs， Guiyang 550006， China）

Abstract： Leaf plays an important role in the growth process of maize， and can effectively carry out photosynthesis， provide nutrients for maize， as well as affect yield increase by affecting density tolerance. In this paper， 150 F2 and F2∶3 family materials were constructed， using QR273 and T32 as the parents. Based on the genotype data and the phenotypic evaluation data of leaf related traits under different environments， quantitative trait locus （QTL） mapping was carried out by using inclusive composite interval mapping method. The results showed that， a total of 85 QTLs for leaf related traits were detected under two environments， including 12 QTLs relevant to total leaf number， 14 QTLs relevant to leaf number above the ear， 22 QTLs relevant to leaf length， 17 QTLs relevant to leaf width， and 20 QTLs relevant to leaf angle. Combined with the public database and bioinformatics analysis methods， seven candidate genes controlling leaf related traits were screened. Wherein Zm00001d013612 encoded tubulin and participated in the regulation of cytoskeleton structure composition， Zm00001d053543 involved in the brassinosteroid-mediated signaling pathway， and Zm00001d031291 encoded protein had the function of histone acetylation， Zm00001d031292 involved in the regulation of genes encoding hydroxyproline-rich glycoprotein family. Zm00001d031296 regulated the activity of potassium ion transmembrane transporter. Zm00001d031300 and Zm00001d031303 involved in carbohydrate metabolism. Functional analysis of proteins showed that， the seven candidate genes all involved in cell differentiation， which were closely related with plant growth and development. These results can provide relatively more theoretical support for deeply revealing the genetic basis of maize leaf related trait variation.

Key words： maize；leaf related traits；quantitative trait locus （QTL）；candidate gene

玉米（Zea mays L.）作為全球主要的糧飼兼用作物及主要的工業(yè)原料，在保障國家糧食安全及畜牧業(yè)發(fā)展中發(fā)揮著至關(guān)重要的作用^[1]。葉部性狀作為影響玉米種植密度的重要組成因素，對植物的光合作用、水分吸收、干物質(zhì)積累以及產(chǎn)量形成等都起著重要作用^[2]。玉米葉部性狀是受多個微效基因調(diào)控的數(shù)量性狀，包括全株葉片數(shù)、穗上葉片數(shù)、葉長、葉寬、葉夾角等，易受環(huán)境影響。目前，已有較多研究者對玉米葉部性狀進行了研究。奉杰等^[3]以用玉米自交系ZNC442和SCML0849為親本構(gòu)建的131份F2∶3家系為材料，在1號、2號、3號、5號、6號、8號、9號、10號染色體上定位到15個與穗上葉片數(shù)相關(guān)的數(shù)量性狀座位（QTL），貢獻率最大的為13.03%；定位到21個與葉夾角相關(guān)的QTL，可解釋2.47%～9.43%的表型變異。羅巧玲等^[4]以鄭58和B73構(gòu)建的雙親分離群體為材料，共定位到11個控制玉米葉長的QTL，單個QTL的表型貢獻率為3.69%～25.53%。王會濤等^[5]利用豫82和豫87-1構(gòu)建的重組自交系群體進行試驗，共定位到5個與葉長相關(guān)的QTL，單個QTL的表型貢獻率為8.72%～13.43%；與葉寬相關(guān)的QTL有7個，單個QTL的表型貢獻率為6.93%～11.68%；與葉夾角相關(guān)的QTL有7個，單個 QTL 的表型貢獻率為 7.08%～11.59%。這些研究結(jié)果為玉米葉部性狀遺傳機制解析和基于分子輔助的葉部性狀遺傳改良利用提供了較為豐富的分子遺傳學(xué)依據(jù)。但是，比較不同的研究成果，發(fā)現(xiàn)所用的材料、分子標記類型和數(shù)量都存在差異，所得到的QTL區(qū)段也存在差異。玉米葉部性狀不僅有著復(fù)雜的遺傳機制和基因調(diào)控網(wǎng)絡(luò)，還受到遺傳背景的影響。本研究擬以150份由玉米骨干自交系QR273和T32為親本構(gòu)建的F2、F2∶3家系為材料，結(jié)合簡化基因組測序（GBS）結(jié)果和多個環(huán)境下葉部性狀表型數(shù)據(jù)，進行QTL定位，利用生物信息學(xué)分析方法對定位結(jié)果進行候選基因的預(yù)測，研究結(jié)果可為玉米葉部性狀改良提供更為豐富的理論支持。

1 材料與方法

1.1 試驗材料

本研究選用的親本材料為QR273和T32，構(gòu)建150份F2、F2∶3家系。雙親均為熱帶玉米骨干自交系，并在葉部性狀表現(xiàn)上具有明顯差異，親本T32葉部性狀表現(xiàn)為葉夾角大、葉較長、葉較寬，親本QR273葉部性狀表現(xiàn)為葉夾角小、葉較短、葉較窄，全株葉片數(shù)較少。

1.2 田間試驗及性狀調(diào)查

2021年冬天，在海南種植親本及F2代群體材料，將套袋自交后得到的F2∶3群體種植于甘肅省張掖市、貴州省貴陽市2地，試驗采取隨機區(qū)組設(shè)計，設(shè)置行長為3.00 m，行距為0.65 m，重復(fù)2次。在授粉20 d后，參考石云素^[6]對葉部性狀的考察標準，調(diào)查全株葉片數(shù)、穗上葉片數(shù)、葉長、葉寬、葉夾角，調(diào)查標準和參數(shù)設(shè)置如下：

全株葉片數(shù)：玉米第一節(jié)間以上的葉片數(shù)，單位為片；穗上葉片數(shù)：玉米有效穗部葉及以上葉片數(shù)，單位為片；葉長：有效穗上第一葉葉舌至葉尖長度，單位為cm；葉寬：有效穗上第一葉葉中部的寬度，單位為cm；葉夾角：有效穗上第一葉與主莖間的角度，單位為°。

1.3 基因型鑒定

對田間的F2代植株進行單株幼苗的DNA提取，將提取液送至北京康普生生物有限公司進行基因組測序，詳細的過程參考Wu等^[7]的方法。

1.4 數(shù)據(jù)整理與分析

利用Excel2021和SPSS26.0軟件對調(diào)查數(shù)據(jù)進行整理和分析，計算2個環(huán)境下的材料在群體間的全株葉片數(shù)（LN）、穗上葉片數(shù)（LAN）、葉長（LL）、葉寬（LW）、葉夾角（LA）的平均值與變異系數(shù)、方差、峰度、偏度和相關(guān)性系數(shù)等，家系遺傳力的計算公式為：

H2per mean=δ2gδ2g+δ2gee+δ2εer

δ2g為基因型方差，δ2ge為基因型×環(huán)境互作方差，δ2ε為誤差方差，m為環(huán)境個數(shù)，ｒ為重復(fù)數(shù)。

1.5 QTL定位及候選基因預(yù)測

去掉等位基因頻率（MAF）lt;0.05以及缺失率大于20%的單核苷酸多態(tài)性（SNP）位點后得到68 994個高質(zhì)量SNP標記。QTL定位使用的軟件為IciMapping4.1，其中作圖步長為1 cM、逐步回歸標記進入的概率（PIN）為0.001、對數(shù)優(yōu)勢比（LOD值）為2.5。對獲得的QTL結(jié)果進行一致性分析，一致性QTL判定標準參照奉杰等^[3]的方法，即與前人定位到的同一性狀的QTL區(qū)段存在重疊區(qū)段。隨后，利用生物學(xué)信息法在maize GDB的GBrowse網(wǎng)站（https：// maizegdb.org/gbrowse）上對一致性QTL區(qū)段進行候選基因及其功能注釋的查找，再結(jié)合前人的相關(guān)研究報道，篩選出關(guān)鍵候選基因。

2 結(jié)果與分析

2.1 葉部性狀表型分析

從表1可看出，自交系QR273的穗上葉片數(shù)平均值高于T32，而T32的全株葉片數(shù)、葉長、葉夾角平均值高于QR273。對于F2∶3群體的分析結(jié)果表明，5個性狀的偏度和峰度的絕對值均小于3.00，符合正態(tài)分布，且各性狀在分離群體中差異顯著（P＜0.05），符合QTL定位分析的條件。在不同環(huán)境中遺傳力存在一定的差異，其中葉長和葉夾角間遺傳力差異較大，葉長的遺傳力在貴陽為0.90，在張掖為0.62；葉夾角的遺傳力在貴陽為0.60，在張掖為0.91。

2.2 相關(guān)性分析

葉部性狀間的相關(guān)性分析結(jié)果見表2。結(jié)果顯示，全株葉片數(shù)與穗上葉片數(shù)、葉寬，穗上葉片數(shù)與葉寬，在2個環(huán)境下均呈極顯著正相關(guān)，而葉長與葉夾角在貴陽環(huán)境下呈極顯著負相關(guān)，穗上葉片數(shù)與葉長在貴陽環(huán)境下呈顯著正相關(guān)。

2.3 高密度物理圖譜構(gòu)建

對獲得的SNP位點進行過濾后得到68 994個高質(zhì)量SNP標記。標記數(shù)量在1號染色體到10號染色體上依次為8 725個、8 824個、7 623個、8 901個、7 174個、5 432個、5 592個、6 407個、5 782個、4 534個，標記數(shù)量的大小情況表現(xiàn)為： 4號染色體＞2號染色體＞1號染色體＞3號染色體＞5號染色體＞8號染色體＞9號染色體＞7號染色體＞6號染色體＞10號染色體，染色體長度分別為308.42Mb、243.67 Mb、237.96 Mb、250.18 Mb、226.35 Mb、181.35 Mb、185.78 Mb、182.31 Mb、162.90 Mb、152.37 Mb，平均標記間距為30.89 kb。

2.4 控制葉部性狀QTL

如表3和圖1所示，共檢測到85個葉部相關(guān)性狀QTL。全株葉片數(shù)相關(guān)QTL有12個，分布于1號、2號、4號、5號、6號、7號染色體上，其中張掖有6個，貴陽有6個，單個QTL可解釋的表型貢獻率為1.07%～9.40%；穗上葉片數(shù)相關(guān)QTL有14個，分布于10條染色體上，其中張掖有8個，貴陽有6個，單個QTL可解釋的表型貢獻率為4.24%～11.98%；葉長相關(guān)QTL有22個，除6號染色體以外的9條染色體上均有分布，其中張掖有7個，貴陽有15個，單個QTL可解釋的表型貢獻率為2.70%～8.44%；葉寬相關(guān)QTL有17個，分布于1號、2號、3號、4號、5號、9號、10號染色體上，其中張掖有10個，貴陽有7個，單個QTL可解釋的表型貢獻率為1.58%～12.91%；葉夾角相關(guān)QTL有20個，10條染色體上均有分布，其中張掖有6個，貴陽有14個，單個QTL可解釋的表型貢獻率為2.43%～7.43%。由表3還可以看出，共檢測到有2個穗上葉片數(shù)相關(guān)QTL和1個葉寬相關(guān)QTL貢獻率大于10.00%，并且在10號染色體上發(fā)現(xiàn)1個QTL富集區(qū)域，即在10號染色體上的138 770 494～139 388 961 bp區(qū)間檢測到1個同時控制玉米葉寬和葉夾角的QTL，表現(xiàn)出一因多效的現(xiàn)象。

2.5 “一致性”QTL及候選基因預(yù)測

對比前人的研究結(jié)果^[8-12]，共找到7個“一致性”QTL，分別控制葉長、葉寬和葉夾角的變異。在這7個“一致性”QTL中查找相關(guān)基因的功能注釋，再根據(jù)前人的研究成果，初步預(yù)測出7個候選基因（表4）。蛋白質(zhì)功能分析結(jié)果顯示，Zm00001d013612編碼微管蛋白β鏈并且參與調(diào)控細胞骨架結(jié)構(gòu)組成、微管細胞骨架組織的組成；Zm00001d053543參與油菜素甾醇介導(dǎo)的信號通路；Zm00001d031291編碼組蛋白乙酰化并且參與DNA甲基化調(diào)控；Zm00001d031292 參與富含羥脯氨酸糖蛋白家族基因表達的調(diào)控；Zm00001d031296調(diào)控鉀離子跨膜轉(zhuǎn)運蛋白活性；Zm00001d031300、Zm00001d031303參與碳水化合物代謝過程。這些基因參與的蛋白質(zhì)編輯功能可能與植物葉片發(fā)育密切相關(guān)。

3 討論

3.1 本研究與前人研究結(jié)果的比較

本研究共定位到85個葉部性狀相關(guān)QTL，其中有7個QTL與前人定位到的區(qū)段存在重合，本研究在4號染色體232.9～233.6 Mb和10號染色體33.9～34.1 Mb定位到的葉長相關(guān)QTL，與趙文明^[8]利用簡單重復(fù)序列（SSR）標記對豫82和沈137構(gòu)建的F2∶3家系群體進行QTL定位，在第4染色的120.2～244.5 Mb（B73RefGen_v4）和第10染色體的13.5～85.2 Mb（B73RefGen_v4）定位到的葉長相關(guān)QTL存在重疊，但本研究鑒定出的QTL區(qū)間更小。此外，張曠野^[9]通過SSR標記對掖478和齊319構(gòu)建的300份重組自交系進行QTL定位，在1號染色體的184.8～193.6 Mb（B73RefGen_v4）定位到1個與葉寬相關(guān)QTL，大小為8.8 Mb。本研究在1號染色體185.3～186.1 Mb定位到1個葉寬相關(guān)QTL，大小為0.8 Mb，定位區(qū)間顯著小于前人定位結(jié)果。路明等^[10]利用SSR標記以掖 478和丹 340構(gòu)建的F2、F2∶3群體為材料，在164.2～194.9 Mb定位到1個控制葉夾角的QTL，與本研究在3號染色體的175.19～175.35 Mb定位到的葉夾角相關(guān)QTL存在重疊。常立國等^[11]以許178×K12構(gòu)建的150份重組自交系為材料，利用SSR標記對葉部性狀進行QTL定位，在第3染色體的172.8～190.3 Mb（B73RefGen_v4）和第10染色體的115.2～128.2 Mb（B73RefGen_v4）各定位到1個與葉夾角相關(guān)的QTL，與本研究定位結(jié)果相似，即本研究在3號染色體175.2～175.4 Mb（B73RefGen_v4）和10號染色體116.9～117.1 Mb（B73RefGen_v4）也各定位到1個與葉夾角相關(guān)的QTL，大小均為0.2 Mb，定位區(qū)間進一步縮小。劉鵬飛等^[12]在5號染色體的4.3～55.7 Mb（B73RefGen_v4）定位到1個控制玉米葉夾角的QTL，本研究在3號染色體175.2～175.4 Mb（B73RefGen_v4）和10號染色體116.9～117.1 Mb內(nèi)也各定位到1個與葉夾角相關(guān)的QTL，在這些“一致性”QTL表達較穩(wěn)定、可靠性較高的基礎(chǔ)上，本研究定位到的QTL較前人定位區(qū)間更小，更利于下一步精細定位和候選基因的圖位克隆。

3.2 候選基因預(yù)測分析

本研究共篩選出7個關(guān)鍵候選基因，這些基因參與的激素調(diào)控和蛋白質(zhì)表達，可能與植物葉片生長發(fā)育密切相關(guān)。其中，Zm00001d013612編碼微管蛋白β鏈并且參與調(diào)控細胞骨架結(jié)構(gòu)組成、微管細胞骨架組織的組成。根據(jù)前人的研究結(jié)果，細胞骨架包括微絲、微管和中間纖維^[13]，微管在植物中具有維持植物細胞形態(tài)^[14-15]、細胞運動^[16]、細胞內(nèi)物質(zhì)運輸^[17-18]、細胞分裂^[19]以及細胞壁形成等重要作用^[20]。Zm00001d053543參與油菜素甾醇介導(dǎo)的信號通路，在植物中參與莖的伸長、細胞分裂和分化的調(diào)節(jié)^[21-24]，參與了幾乎所有的植物生長發(fā)育過程，被證實在水稻中能夠調(diào)控株高、葉夾角、籽粒大小^[25]。Zm00001d031291編碼組蛋白乙?；⑶覅⑴cDNA甲基化調(diào)控，影響植物的分生組織形成^[26]、細胞的增殖^[27]等過程，與器官的生長過程密切相關(guān)，在細胞內(nèi)通過影響基因的轉(zhuǎn)錄和表達，從而影響細胞分化的過程^[28]。Zm00001d031292參與富含羥脯氨酸糖蛋白家族表達基因的調(diào)控，而此蛋白質(zhì)是高等植物細胞壁中特有的，具有強固細胞壁的作用^[29]。Zm00001d031296與鉀離子跨膜轉(zhuǎn)運蛋白活性密切相關(guān)，鉀營養(yǎng)能夠影響油菜葉片生長過程中葉肉細胞的分布和形態(tài)，對葉面積擴張、CO2和H2O的運輸過程具有重要作用^[30]。 Zm00001d031300、Zm00001d031303參與碳水化合物的代謝過程，碳水化合物能夠參與植物細胞壁的形成，是植物的主要能量來源，Blumstein等^[31]研究發(fā)現(xiàn)非結(jié)構(gòu)性碳水化合物能夠向細胞發(fā)出時間流逝的信號，觸發(fā)葉片發(fā)育，這些基因可能通過影響細胞的生長發(fā)育或參與激素的調(diào)控等方式影響植物葉部性狀的生長發(fā)育，可初步預(yù)測為關(guān)鍵候選基因。本研究結(jié)果可為后續(xù)的基因精細定位、克隆以及有效開展玉米葉部性狀分子改良研究提供支撐。

3.3 葉部性狀QTL在育種上的利用

目前也有較多研究者通過QTL定位發(fā)現(xiàn)控制玉米葉部性狀的關(guān)鍵基因，Gao等^[32]鑒定到葉寬相關(guān)QTL的候選基因ZmNL4在被敲除后玉米葉片寬度顯著降低， Xia等^[33]通過全基因組關(guān)聯(lián)分析發(fā)現(xiàn)了1個位于基因RHW1上的SNP位點與玉米殼葉寬度顯著相關(guān)，并證實RHW1-ZCN4調(diào)控途徑會影響玉米殼葉的寬度，Tian等^[34]研究證實ZmRAVL1能夠通過影響葉鞘和韌帶區(qū)域正向厚壁組織細胞數(shù)目從而正向調(diào)節(jié)玉米葉夾角。植物中還存在許多基因通過參與激素調(diào)控來控制葉部性狀。例如擬南芥中ROTUNDIFOLIA3/4編碼的CYP90C1能夠通過合成油菜素內(nèi)脂，影響細胞的極性擴展，從而對葉片長度起到調(diào)控作用^[35-37]。Uzair等^[38]研究發(fā)現(xiàn)窄葉突變體生長素響應(yīng)因子ARF的蛋白質(zhì)翻譯效率降低后會導(dǎo)致生長素調(diào)控的葉片細胞分裂和擴展受阻，從而影響葉片寬度。Best等^[39]研究發(fā)現(xiàn)影響油菜素甾醇合成的基因Nana plant2突變后會導(dǎo)致玉米矮小、葉片形態(tài)變化和葉夾角增大^[40]?？偟膩碚f，植物葉部性狀的生長發(fā)育是一個復(fù)雜的過程，具體的調(diào)控網(wǎng)絡(luò)還有待進一步的深入挖掘。QTL定位是數(shù)量遺傳學(xué)的主要研究課題之一，高分辨率的遺傳標記圖譜能夠定位到更精確的QTL區(qū)段，基于這些研究成果，后續(xù)可對所得到的關(guān)鍵遺傳區(qū)段進行精細定位，利用亞細胞定位、β-葡萄糖醛酸糖苷酶（GUS）染色、遺傳轉(zhuǎn)化（過表達、干擾、基因編輯等）對候選基因進行功能驗證，開展蛋白質(zhì)互作相關(guān)試驗，通過酵母雙雜、免疫共沉淀（CO-IP）、雙分子熒光互補（BiFC），或是蛋白質(zhì)與核酸互作研究等方式對目標基因的功能機制進行更深入的研究。綜上所述，本研究結(jié)果能夠?qū)罄m(xù)深度揭示玉米葉部性狀變異的遺傳機制、相關(guān)候選基因克隆以及基于分子輔助選擇育種提供參考。

4 結(jié)論

本研究共檢測到85個葉部性狀相關(guān)QTL，單個QTL的貢獻率為1.07%～12.91%，與前人研究結(jié)果進行比對，共發(fā)現(xiàn)7個“一致性”QTL，并且定位到的區(qū)間更小。共預(yù)測出7個候選基因，分別是Zm00001d013612、Zm00001d053543、Zm00001d031291、Zm00001d031292、Zm00001d031296、Zm00001d031300、Zm00001d031303，所參與的細胞生長發(fā)育、蛋白質(zhì)和激素調(diào)控在作物生長發(fā)育過程中起著重要作用，可能是玉米葉部性狀變異的關(guān)鍵候選基因，可以作為后續(xù)QTL的精細定位、候選基因圖位克隆和功能研究新的參照。

參考文獻：

[1] 張 強，趙曉陽，張 潔，等. 國審玉米品種寶景186及其栽培技術(shù)[J]. 中國種業(yè)，2023（7）：127-129.

[2] 高祥斌，張秀省，蔡連捷. 觀賞植物葉面積測定及相關(guān)分析[J]. 福建林業(yè)科技，2009，36 （2）：231-234，251.

[3] 奉 杰，劉鵬飛，涂 亮，等. 玉米株型相關(guān)性狀的QTL定位與候選基因分析[J]. 玉米科學(xué)，2023，31（4）：34-42.

[4] 羅巧玲，孫 強，張瀟譽，等. 玉米葉長QTL定位和全基因組選擇[J]. 新疆農(nóng)業(yè)大學(xué)學(xué)報，2021，44（1）：14-19.

[5] 王會濤，柳華峰，鄭耀剛，等. 玉米葉型相關(guān)性狀QTL定位及上位性效應(yīng)分析[J]. 分子植物育種，2018，16（15）：4955-4963.

[6] 石云素. 玉米種質(zhì)資源描述規(guī)范和數(shù)據(jù)標準[M]. 北京：中國農(nóng)業(yè)出版社，2006.

[7] WU X， GUO X Y， WANG A G， et al. Quantitative trait loci mapping of plant architecture-related traits using the high-throughput genotyping by sequencing method[J]. Euphytica，2019，215：1-13.

[8] 趙文明. 玉米株型相關(guān)性狀QTL定位與分析[D]. 鄭州：河南農(nóng)業(yè)大學(xué)，2008.

[9] 張曠野. 玉米葉夾角和葉寬的遺傳分析及QTL定位[D]. 沈陽：沈陽農(nóng)業(yè)大學(xué)，2018.

[10]路 明，周 芳，謝傳曉，等. 玉米雜交種掖單13號的SSR連鎖圖譜構(gòu)建與葉夾角和葉向值的QTL定位與分析[J]. 遺傳，2007，29（9）：1131-1138.

[11]常立國，何坤輝，劉建超，等. 不同環(huán)境條件下玉米葉夾角的QTL定位[J]. 玉米科學(xué)，2016，24（4）：49-55.

[12]劉鵬飛，蔣 鋒，王漢寧，等. 玉米葉夾角和葉向值的QTL定位[J]. 核農(nóng)學(xué)報，2012，26（2）：231-237.

[13]王潔茹，戴 銀，沈?qū)W懷，等. 冠狀病毒與宿主細胞骨架相互作用的研究進展[J]. 中國獸醫(yī)科學(xué)，2021，51（11）：1433-1437.

[14]GARDINER J. The evolution and diversification of plant microtubule-associated proteins[J]. The Plant Journal，2013，75（2）：219-229.

[15]SUN W， LI Y， ZHAO Y X， et al. The TsnsLTP4， a nonspecific lipid transfer protein involved in wax deposition and stress tolerance[J]. Plant Molecular Biology Reporter，2015，33：962-974.

[16]WANG X， ZHU L， LIU B Q， et al. Arabidopsis MICROTUBULE-ASSOCIATED PROTEIN18 functions in directional cell growth by destabilizing cortical microtubules[J]. The Plant Cell，2007，19（3）：877-889.

[17]ZHOU S， CHEN Q H， LI X Y， et al. MAP65-1 is required for the depolymerization and reorganization of cortical microtubules in the response to salt stress in Arabidopsis[J]. Plant Science，2017，264：112-121.

[18]LI S P， CHEN M， YU D L， et al. EXO70A1-mediated vesicle trafficking is critical for tracheary element development in Arabidopsis[J]. Plant Cell，2013，25（5）：1774-1786.

[19]WASTENEYS G O. Progress in understanding the role of microtubules in plant cells [J]. Current Opinion In Plant Biology，2004，7（6）：651-660.

[20]馬彩霞. 陸地棉GhTUBB1基因的功能研究[D]. 太原：山西農(nóng)業(yè)大學(xué)，2022.

[21]ASAMI T， NAKANO T， FUJIOKA S. Plant brassinosteroid hormones[J]. Vitamins amp; Hormones，2005，72：479-504.

[22]LI J H， LI Y H， CHEN S Y， et al. Involvement of brassinosteroid signals in the floral-induction network of Arabidopsis[J]. Journal of Experimental Botany，2010，61（15）：4221-4230.

[23]BASIT F， LIU J X， AN J Y， et al. Brassinosteroids as a multidimensional regulator of plant physiological and molecular responses under various environmental stresses[J]. Environmental Science and Pollution Research， 2021，28：44768-44779.

[24]張存家，劉小青，韓 婧，等. 基于轉(zhuǎn)錄組測序的紅樹秋茄葉片發(fā)育中差異表達基因分析[J]. 分子植物育種，2023，21（24）：8075-8083.

[25]楊彥釗，童紅寧. 水稻中油菜素甾醇功能機制解析與分子設(shè)計利用[J]. 生命的化學(xué)，2021，41（6）：1171-1180.

[26]HOLLENDER C， LIU Z C. Histone deacetylase genes in Arabidopsis development[J]. Journal of Integrative Plant Biology，2008，50（7）：875-885.

[27]LIU X C， YANG S G， ZHAO M L， et al. Transcriptional repression by histone deacetylases in plants[J]. Molecular Plant，2014，7（5）：764-772.

[28]李甜甜. 小茴香揮發(fā)油促進擬南芥生長的機制研究[D]. 昆明：云南中醫(yī)藥大學(xué)，2023.

[29]宋鳳鳴，鄭 重，葛起新. 富含羥脯氨酸糖蛋白在植物-病原物相互作用中的積累、作用及調(diào)控[J]. 植物生理學(xué)通訊，1992，（2）：141-145.

[30]HU W S， LU Z F， GU H H， et al. Potassium availability influences the mesophyll structure to coordinate the conductance of CO2 and H2O during leaf expansion[J]. Plant，Cell amp; Environment，2022， 45（10）：2987-3000.

[31]BLUMSTEIN M， OSEGUERA M， CASO-MCHUGH T， et al. Nonstructural carbohydrate dynamics’ relationship to leaf development under varying environments[J]. New Phytologist，2024，241（1）：102-113.

[32]GAO L L， YANG G H， LI Y F， et al. A kelch-repeat superfamily gene， ZmNL4， controls leaf width in maize （Zea mays L. ）[J]. The Plant Journal，2021，107（3）：817-830.

[33]XIA A A， ZHENG L M， WANG Z， et al. The RHW1-ZCN4 regulatory pathway confers natural variation of husk leaf width in maize[J]. New Phytologist，2023，239（6）：2367-2381.

[34]TIAN J G， WANG C L， XIA J L， et al. Teosinte ligule allele narrows plant architecture and enhances high-density maize yields[J]. Science，2019，365（6454）：658-664.

[35]KIM G T， TSUKAYA H， UCHIMIYA H. The ROTUNDIFOLIA3 gene of Arabidopsis thalianaencodes a new member of the cytochrome P-450 family that is required for the regulated polar elongation of leaf cells[J]. Genes amp; Development，1998，12（15）：2381-2391.

[36]KIM G T， FUJIOKA S， KOZUKA T， et al. CYP90C1 and CYP90D1 are involved in different steps in the brassinosteroid biosynthesis pathway in Arabidopsis thaliana[J]. The Plant Journal，2005，41（5）：710-721.

[37]崔曉峰，黃 海. 葉發(fā)育的遺傳調(diào)控機理研究進展[J]. 植物生理學(xué)報，2011，47（7）：631-640.

[38]UZAIR M， LONG H， ZAFAR S A， et al. Narrow Leaf21，encoding ribosomal protein RPS3A，controls leaf development in rice[J]. Plant Physiology，2021，186（1）：497-518.

[39]BEST N B， HARTWIG T， BUDKA J， et al. nana plant2 encodes a maize ortholog of the Arabidopsis brassinosteroid biosynthesis gene DWARF1，identifying developmental interactions between brassinosteroids and gibberellins[J]. Plant Physiology，2016，171（4）：2633-2647.

[40]張全艷，張培高，徐春霞，等. 玉米葉夾角的遺傳與分子調(diào)控研究進展[J]. 中國農(nóng)業(yè)科技導(dǎo)報，2021，23（10）：15-24.

（責(zé)任編輯：陳海霞）