閆 偉 李 元 宋茂興 張曠野 孫銘澤 瞿 會(huì) 李鳳海鐘雪梅 朱 敏 杜萬(wàn)里 呂香玲

沈陽(yáng)農(nóng)業(yè)大學(xué)特種玉米研究所， 遼寧沈陽(yáng) 110866

?

玉米抗灰斑病QTL元分析及其驗(yàn)證

閆 偉 李 元 宋茂興 張曠野 孫銘澤 瞿 會(huì) 李鳳海鐘雪梅 朱 敏 杜萬(wàn)里 呂香玲*

沈陽(yáng)農(nóng)業(yè)大學(xué)特種玉米研究所， 遼寧沈陽(yáng) 110866

摘 要：玉米灰斑病是危害玉米生產(chǎn)的主要病害之一， 目前對(duì)抗灰斑病基因數(shù)目、位置及作用方式仍然不清楚， 這嚴(yán)重制約著玉米抗灰斑病育種進(jìn)展。本研究利用元分析方法分析并整理了14篇玉米抗灰斑病QTL文獻(xiàn)的信息， 共篩選確定了13個(gè)一致性QTL區(qū)間。利用以自交系81162為輪回親本、自交系CN165為非輪回親本構(gòu)建的回交導(dǎo)入群體根據(jù)連鎖不平衡原理對(duì)13個(gè)一致性QTL進(jìn)行驗(yàn)證， 在13個(gè)一致性QTL區(qū)間共獲得20多個(gè)偏分離位點(diǎn)。第1和第4染色體上偏分離最嚴(yán)重， 其他染色體上偏分離度較小。說明第1和第4染色體上存在著效應(yīng)較大的抗病QTL。第1染色體標(biāo)記umc2227、bnlg1832、umc1243、umc2025、umc1515、umc1297、umc1461處供體基因頻率均在50%以上， 可能存在幾個(gè)連鎖的抗病基因。第4染色體上基因位于標(biāo)記bnlg2291和umc1194之間。研究為精細(xì)定位供體CN165中第1和第4染色體上的抗灰斑病QTL奠定了基礎(chǔ)。

關(guān)鍵詞：灰斑??； 元分析； 一致性QTL； 回交導(dǎo)入系； 主效QTL

本研究由國(guó)家自然科學(xué)基金項(xiàng)目（31301322）和遼寧省科技特派項(xiàng)目（2014215031）資助。

This study was supported by the National Natural Science Foundation of China （31301322） and the Special Talent Appointment in Technology of Liaoning Province （2014215031）.

第一作者聯(lián)系方式∶ E-mail∶ yanwei1452@hotmail.com， Tel∶ 18512468548

URL∶ http∶//www.cnki.net/kcms/detail/11.1809.S.20160314.1444.006.html

玉米灰斑病是由尾孢菌侵染玉米葉部引起的真菌病害， 又被稱為玉米尾孢菌葉斑?。?-2］。玉蜀黍尾孢菌（Cecrospora zeae-maydis Tehon & Daniels）為玉米灰斑病的真菌侵染病原， 主要?dú)埓嬖谕寥乐械牟≈参餁報(bào)w上越冬， 其分生孢子通過風(fēng)和雨水飛濺散布到玉米植株上。溫和潮濕的條件下灰斑病會(huì)進(jìn)一步侵染， 導(dǎo)致嚴(yán)重的葉片衰老和莖稈腐爛。于20世紀(jì)20年代在美國(guó)首次發(fā)現(xiàn)灰斑病， 如今已成為危害玉米生產(chǎn)的主要病害［3］。長(zhǎng)期的生產(chǎn)實(shí)踐證明， 選育和推廣抗病品種是防治玉米灰斑病的最有效的途徑之一［4-5］。

研究表明， 玉米對(duì)灰斑病的抗性屬多基因控制的數(shù)量性狀， 其不同的基因作用方式和效應(yīng)不同， 抗性基因的累積可使品種產(chǎn)生較高的抗性。因此玉米抗灰斑病主效QTL的發(fā)掘?qū)τ衩卓共∮N具有重要意義［6］。自20世紀(jì)90年代以來(lái)， 國(guó)內(nèi)外學(xué)者開展了一些玉米抗灰斑病的基因定位研究， 發(fā)掘出數(shù)量較多的抗病QTL［7-20］。但由于不同研究間存在材料遺傳背景、發(fā)病環(huán)境條件、圖譜標(biāo)記密度等差異， 定位結(jié)果間一致性較差。目前對(duì)抗灰斑病基因數(shù)目、位置及作用方式仍然不清楚， 這嚴(yán)重制約著玉米抗灰斑病育種進(jìn)展。

元分析（Meta-analysis）統(tǒng)計(jì)方法是對(duì)眾多現(xiàn)有實(shí)證文獻(xiàn)的再次統(tǒng)計(jì)， 通過將不同文獻(xiàn)的QTL位點(diǎn)統(tǒng)一投射到參考圖譜上， 根據(jù)每個(gè)QTL的中心位置及其置信區(qū)間， 即可以根據(jù)最大概率的原理估計(jì)出真實(shí)QTL的位置、縮小QTL的置信區(qū)間。目前， 高密度的玉米IBM遺傳圖譜為不同試驗(yàn)QTL整合和一致性圖譜構(gòu)建奠定了基礎(chǔ)［21］。Chardon等［22］提出利用overview和元分析的方法對(duì)玉米花期相關(guān)的313個(gè)QTL信息進(jìn)行整合分析， 考慮QTL定位中作圖群體及其大小、QTL效應(yīng)等因素， 最終獲得62個(gè)一致性QTL。Truntzler等［23］利用元分析整合了11個(gè)不同材料的與青貯玉米品質(zhì)相關(guān)的QTL原始圖譜， 得到一致性QTL。這種基于性狀定位的比較研究不僅有助于發(fā)掘一致性主效基因位點(diǎn)及其連鎖標(biāo)記， 同時(shí)為基因/QTL的克隆和利用提供了新手段［24］。

雜交、回交與自交方法相結(jié)合， 構(gòu)建高世代回交導(dǎo)入系（或近等基因系）基因定位的理論， 成為目前開展性狀理論研究和育種實(shí)踐相結(jié)合的一個(gè)有效途徑。因?yàn)闃?gòu)建的回交導(dǎo)入系， 其每一導(dǎo)入系與導(dǎo)入系， 或與輪回親本僅在某幾個(gè)染色體片段上有差異， 這樣的植株適用于性狀分析，可作為目標(biāo)性狀遺傳與生理學(xué)基礎(chǔ)研究分析的材料； 同時(shí)導(dǎo)入系在遺傳背景上繼承了輪回親本的主要優(yōu)良基因，基本保持了原有輪回親本的一般配合力與特殊配合力，可以快速應(yīng)用于育種［25］。如徐建龍等［26］、Li等［27］、李芳等［28］、Xie等［29］分別利用回交導(dǎo)入系進(jìn)行了水稻抗旱、抗病及產(chǎn)量品質(zhì)相關(guān)的QTL研究； Ho等［30］、呂香玲等［31］、Li等［32］、Salvi等［33］、Pea等［34］、Teng等［35］分別利用導(dǎo)入系進(jìn)行了玉米抗旱、抗病和產(chǎn)量相關(guān)QTL的定位研究。

本文利用元分析方法對(duì)不同玉米群體的抗灰斑病QTL進(jìn)行整合分析， 并對(duì)分析結(jié)果進(jìn)行有條件的篩選，確定玉米抗灰斑病一致性 QTL。然后利用回交導(dǎo)入系群體進(jìn)行一致性QTL驗(yàn)證， 發(fā)掘玉米抗灰斑病主效QTL。

1 材料與方法

1.1 玉米抗灰斑病QTL文獻(xiàn)下載及格式文件整理

從NCBI、SpringerLink等網(wǎng)站下載已發(fā)表的玉米抗灰斑病QTL相關(guān)文獻(xiàn)， 根據(jù)元分析軟件格式要求按照QTL名稱、染色體位置、置信區(qū)間、臨近標(biāo)識(shí)、群體類型等信息整理數(shù)據(jù)。

1.2 玉米抗灰斑病QTL信息整合及元分析

通過閱讀文獻(xiàn)， 將整理好的每個(gè)QTL的相關(guān)信息按一定格式載入BioMercator V4.2.1軟件中， 以IBM2 2008 Neighbors作為參考圖譜對(duì)所有QTL整合并繪制成整合圖譜。根據(jù)整合結(jié)果， 分別對(duì)第1、第2、第3、……、第10染色體進(jìn)行元分析， 選取模型AIC值最小的為最適合模型，以包括2個(gè)及以上群體和3個(gè)及以上QTL為篩選條件， 確定一致性QTL的最可能位置及置信區(qū)間。

1.3 利用回交導(dǎo)入系進(jìn)行一致性QTL驗(yàn)證

1.3.1 供試材料 以CN165為供體和81162為受體通過2代回交和2代自交構(gòu)建出BC2F3群體， 自交系CN165為熱帶血緣材料， 高抗多種玉米病害， 植株屬于平展型， 較高大； 81162是吉林農(nóng)業(yè)科學(xué)院選育的優(yōu)良自交系， 株型緊湊， 葉片夾角小， 但對(duì)玉米灰斑病敏感。

1.3.2 灰斑病抗性鑒定 2012年在沈陽(yáng)農(nóng)業(yè)大學(xué)試驗(yàn)田種植親本自交系和BC2F3群體， 親本各種植2行， BC2F3群體種植60行， 行長(zhǎng)4 m， 每行17株， 利用人工接種技術(shù)在該群體內(nèi)篩選抗玉米灰斑病單株并自交留種， 每一果穗構(gòu)建為一個(gè)回交導(dǎo)入系。2013年種植親本自交系和回交導(dǎo)入系并人工接種鑒定。采用完全隨機(jī)區(qū)組排列， 單行區(qū)種植每個(gè)回交導(dǎo)入系， 行長(zhǎng)4 m， 每行17株， 2次重復(fù)。鑒定用的玉米灰斑病菌是從沈陽(yáng)地區(qū)采集的玉米感病葉片上的病斑， 按常規(guī)分離法分離培養(yǎng)的純培養(yǎng)物。采用玉米葉粉碳酸鈣瓊脂培養(yǎng)基（MLPCA）進(jìn)行分生孢子的培養(yǎng)。用無(wú)菌水配制成接種用分生孢子懸液， 將其濃度調(diào)至2.5×103個(gè) mL-1。在玉米植株喇叭口期（第9～第11葉期），用噴嘴處裝有20 mL注射器針頭的手提式注射器， 從植株喇叭口處平行插入， 將每株10 mL的病菌孢子懸液注入植株心葉， 接種所有單株。

按抗病性調(diào)查分級(jí)方法［36］， 于玉米乳熟期田間目測(cè)調(diào)查接種部位上、下葉片的發(fā)病情況， 記載病情級(jí)別。當(dāng)鑒定圃中的感病對(duì)照丹340達(dá)到其相應(yīng)的感病程度（7級(jí)）以上時(shí)， 鑒定則視為有效。1級(jí)為葉片上無(wú)病斑或僅有零星病斑， 病斑占葉面積少于或等于5%， 屬高抗（HR）； 3級(jí)為葉片上有少量病斑， 占葉面積6%～10%， 屬抗（R）； 5級(jí)為葉片上病斑較多， 占葉面積11%～30%， 屬中抗（MR）； 7級(jí)為葉片上有大量病斑， 病斑相連， 占葉面積31%～70%，屬感病（S）； 9級(jí)為葉片基本被病斑覆蓋， 葉片枯死， 屬高感（HS）。

1.3.3 回交導(dǎo)入系的基因型 采用CTAB法提取54個(gè)抗玉米灰斑病單株基因組DNA。在MaizeGDB數(shù)據(jù)庫(kù)中根據(jù)獲得的一致性QTL區(qū)間所在的染色體區(qū)域查找SSR標(biāo)記信息， 下載引物序列， 由生工生物工程（上海）股份有限公司合成引物。SSR標(biāo)記分析包括PCR擴(kuò)增、聚丙烯酰胺凝膠電泳、銀染等程序， 參考LABORATARY PROTOCOLS分析方法［37］， 并根據(jù)經(jīng)驗(yàn)優(yōu)化檢測(cè)技術(shù)。首先在親本間進(jìn)行多態(tài)性篩選， 獲得具有多態(tài)性的引物分析導(dǎo)入系基因型。在同等遷移速率下， 用“A”表示后代單株擴(kuò)增條帶與親本CN165條帶相一致， “B”表示與親本81162條帶相一致，“H”表示雜合條帶， “U”表示缺失條帶。

1.3.4 主效QTL分析方法 根據(jù)連鎖不平衡理論分析電泳檢測(cè)結(jié)果， 當(dāng)所選單株的某一標(biāo)記位點(diǎn)基因型頻率與理論期望值之間存在顯著差異時(shí)， 認(rèn)為存在該選擇位點(diǎn)。通常選擇P≤0.0001的顯著水平來(lái)減小誤差。其中供體位點(diǎn)的基因頻率高， 表明存在對(duì)供體基因的有利選擇，反之則不利； 同時(shí)可以推斷在高頻率基因型的標(biāo)記處存在性狀相關(guān)的QTL。當(dāng)這種偏分離相關(guān)性是由供體純合子引起時(shí)， 表明基因間存在著加性效應(yīng)； 如果偏分離是由供體雜合子引起， 表明基因間存在著超顯性效應(yīng)； 若偏分離是由供體純合子和雜合子共同引起， 則表明存在著部分或者完全顯性。

2 結(jié)果與分析

表1 玉米抗灰斑病QTL的文獻(xiàn)信息Table 1 Reference informations of resistance to gray leaf spot of corn

2.1 玉米抗灰斑病QTL文獻(xiàn)及整理

通過瀏覽NCBI等網(wǎng)站， 本文整理了14篇玉米抗灰斑病QTL文獻(xiàn)的信息， 共得到98個(gè)QTL區(qū)間（表1）。

2.2 玉米抗灰斑病QTL的整合與一致性QTL確定

整合已經(jīng)定位的98個(gè)QTL。表明， 在玉米的10條染色體上均有玉米抗灰斑病QTL存在， 每條染色體上分布的數(shù)量不同， 主要集中在第1、第2、第4、第5、第7、第8、第9染色體上。第6染色體上整合的QTL數(shù)最少（3個(gè)QTL），第4條染色體上整合的QTL數(shù)最多（21個(gè)QTL）。各QTL整合到IBM2 2008 Neighbors參考圖譜后的區(qū)間大小差異較大，從最小的1.7 cM到最大的277.83 cM。

對(duì)98個(gè)抗灰斑病QTL進(jìn)行元分析。在軟件給出的5個(gè)分析模型中， 選擇AIC值最小的模型， 共得到34個(gè)QTL區(qū)間（表2）。本文以群體數(shù)≥2、整合QTL數(shù)≥3為篩選條件，最終確定13個(gè)一致性QTL。這些QTL分別位于第1染色體3個(gè)， 置信區(qū)間分別為405.29～409.74 cM、525.28～538.44 cM 和927.07～973.79 cM； 第2染色體4個(gè)， 置信區(qū)間分別為280.51～288.84 cM、372.29～383.96 cM、468.84～486.22 cM 和624.34～635.58 cM； 第4染色體4個(gè)， 置信區(qū)間分別為175.09～187.61 cM、288.88～292.71 cM、433.84～444.70 cM和524.51～525.48 cM； 第7染色體1個(gè)， 置信區(qū)間為186.17～259.98 cM；第8染色體1個(gè)， 置信區(qū)間為378.64～390.84 cM （圖1）。

2.3 回交導(dǎo)入系抗病性表現(xiàn)

2012年對(duì)BC2F3群體單株接種鑒定， 共獲得54個(gè)抗病表現(xiàn)3級(jí)以上的抗病單株并自交留種獲得其回交導(dǎo)入系。2013年對(duì)54個(gè)回交導(dǎo)入系重復(fù)鑒定， 從每一回交導(dǎo)入系每一重復(fù)選取連續(xù)10株調(diào)查病斑覆蓋度， 并由10個(gè)單株病斑覆蓋度平均水平確定該導(dǎo)入系的抗感級(jí)別。對(duì)一個(gè)導(dǎo)入系， 選取2次重復(fù)鑒定結(jié)果中發(fā)病級(jí)別較高的結(jié)果作為導(dǎo)入系的最終抗感級(jí)別。54個(gè)回交導(dǎo)入系的抗病表現(xiàn)鑒定結(jié)果見表3。其中抗病性表現(xiàn)為1級(jí)的家系有21個(gè)， 3級(jí)的家系有19個(gè)， 5級(jí)的家系有14個(gè)。

2.4 標(biāo)記位點(diǎn)基因型分析

根據(jù)獲得的13個(gè)一致性QTL結(jié)果及玉米MaizeGDB信息選取一致性 QTL區(qū)間的標(biāo)記進(jìn)行親本多態(tài)性篩選，獲得74對(duì)在親本間表現(xiàn)多態(tài)性的引物。利用該74對(duì)親本多態(tài)性引物對(duì) 54個(gè)抗病單株進(jìn)行標(biāo)記基因型分析， 并計(jì)算供體基因頻率。供體基因頻率最高達(dá)89%， 其中達(dá)50%以上的標(biāo)記有15個(gè)， 分別位于第1和第4染色體上（表4）。

表2 玉米抗灰斑病QTL的元分析結(jié)果Table 2 QTL meta-analysis of resistance to gray leaf spot of corn

2.5 主效QTL分析

本試驗(yàn)54個(gè)導(dǎo)入系DNA采自BC2F3世代， 根據(jù)計(jì)算可知， 在非選擇的情況下， 群體中某一單株在單一標(biāo)記位點(diǎn)上基因型的理論概率P（AA）∶P（Aa）∶P（aa）=3∶2∶27， 分別代表了來(lái)源于供體親本的基因純合、基因雜合和來(lái)源于輪回親本的基因純合的理論概率。根據(jù)此理論值，計(jì)算3種基因型的χ2值， 利用χ2值的顯著性檢驗(yàn)判斷標(biāo)記位點(diǎn)附近是否存在抗病性QTL。

表4表明， 74個(gè)標(biāo)記位點(diǎn)中20多個(gè)位點(diǎn)出現(xiàn)了嚴(yán)重偏分離現(xiàn)象。其中第 1染色體標(biāo)記位點(diǎn) χ2值最大達(dá)到121.46， 第 2染色體達(dá)64.7， 第 4染色體達(dá)213.98， 第 7染色體達(dá)18.98， 第8染色體達(dá)24.76。由此說明第1和第4染色體上存在主效抗病 QTL， 第 2染色體上抗病 QTL效應(yīng)次之， 其他染色體上的抗病QTL效應(yīng)可能較小。

表3 54個(gè)導(dǎo)入系抗病性表現(xiàn)Table 3 Resistance performance of 54 backcross introgression lines

表4 標(biāo)記位點(diǎn)供體基因頻率及卡平方值Table 4 Frequency of donor parent's genotype and χ2

（續(xù)表4）

（續(xù)表4）

圖2 卡平方顯著性檢驗(yàn)Fig.2 Significance testing by χ2test

根據(jù)IBM2 2008 Neighbors標(biāo)記順序， 分析染色體上χ2值變化。由圖2可見， 第1染色體上χ2值沒有顯著的一個(gè)高峰點(diǎn)， 而是在標(biāo)記 umc2227、bnlg1832、umc1243、umc2025、umc1515、umc1297、umc1461處出現(xiàn)一個(gè)平緩的峰線。第2染色體上χ2高峰點(diǎn)位于標(biāo)記phi435417處。第4染色體上χ2高峰點(diǎn)位于標(biāo)記bnlg2291和umc1194處。第7染色體上χ2高峰點(diǎn)位于標(biāo)記umc2333處。第8染色體上χ2高峰點(diǎn)位于標(biāo)記phi080處。這說明在這些位點(diǎn)附近可能存在抗病QTL。

3 討論

迄今， 已有較多玉米抗灰斑病 QTL定位的研究， 然而仍然未能找到通用性標(biāo)記被育種家應(yīng)用于輔助育種。近年來(lái)， 元分析和生物信息學(xué)方法在整合不同 QTL信息并進(jìn)一步挖掘候選基因上得到了廣泛的應(yīng)用。這一基于性狀的比較定位途徑可以減少分子育種的基本投入， 為開發(fā)復(fù)雜數(shù)量性狀的主效基因及標(biāo)記， 進(jìn)而為克隆基因/QTL提供了技術(shù)支持。本研究通過對(duì)前人研究的總結(jié)及元分析處理， 共得到13個(gè)一致性QTL。利用回交導(dǎo)入群體根據(jù)連鎖不平衡原理對(duì)13個(gè)一致性 QTL進(jìn)行驗(yàn)證， 在13個(gè)一致性QTL區(qū)間共獲得20多個(gè)偏分離位點(diǎn)。第1和第4染色體上偏分離最嚴(yán)重， 其他染色體上偏分離度較少。說明第1和第4染色體上存在著效應(yīng)較大的抗病QTL。從χ2高峰點(diǎn)來(lái)看， 第 1染色體比較復(fù)雜， 標(biāo)記 umc2227、bnlg1832、umc1243、umc2025、umc1515、umc1297、umc1461 處 χ2值基本接近， 均較高， 供體基因頻率也均在 50%以上，這說明第1染色體抗病基因表現(xiàn)十分復(fù)雜， 可能存在幾個(gè)緊密連鎖的抗病基因。這也解釋了為什么多數(shù)研究在第1染色體上均定位到QTL， 但位置差異較大。第4染色體上χ2值最高， 而且供體基因頻率達(dá)到最大值89%， 這說明在本研究供體親本中第 4染色體上存在一個(gè)效應(yīng)較大的抗病基因。在前人的研究中， 第 4條染色體上獲得的 QTL數(shù)最多（21個(gè)）。如Bubeck等［19］發(fā)現(xiàn)在10條染色體上均存在抗病QTL； Saghai Maroof等［12］發(fā)現(xiàn)在第1、第2、第4、第8、第10染色體上有QTL； Gordon等［8，38］發(fā)現(xiàn)抗病QTL主要分布于第2和第4染色體的長(zhǎng)臂上； Zwonitzer等［16］定位的QTL分別位于染色體bins 1.05～1.06，1.08～1.09， 2.02～2.03， 3.04～3.05， 4.05， 7.02， 8.07～8.09，10.05； Benson等［39］發(fā)現(xiàn)在第1、第2、第4染色體上存在QTL。綜合玉米抗灰斑病基因（QTL）定位研究表明， 玉米10條染色體上均可能存在抗病 QTL， 而且不同材料控制灰斑病抗性的 QTL數(shù)目與位置不同， 多數(shù)研究都含有2～3個(gè)主效QTL。如Gordon等［8］分別在染色體bin2.09和bin 4.08定位到2個(gè)主效QTL； Lehmensiek等［9］研究發(fā)現(xiàn)染色體 bin1.04和 bin2.04各存在一個(gè)效應(yīng)值較高的主效QTL； Veiga等［13］得到位于bin4.03、bin4.04和bin4.05的3個(gè)主效QTL； Chung等［15］定位出2個(gè)主效QTL， 分別位于bin1.04和bin8.05； Zwonitzer等［16］檢測(cè)到3個(gè)主效QTL，分別位于bin1.05、bin1.08和bin2.02； Asea等［17］分別在染色體bin2.09和bin4.08檢測(cè)到主效QTL； Balint等［18］定位到3個(gè)與抗病有關(guān)的主效QTL， 分別位于bin1.05、bin2.04 和bin4.05。本研究中供體CN165中至少含有2個(gè)以上主效抗病QTL， 它們分別位于第1和第4染色體。

References

［1］王桂清， 陳捷.玉米灰斑病抗病性研究進(jìn)展.沈陽(yáng)農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2000， 31∶ 418-422 Wang G Q， Chen J.The research progress of resistance to gray leaf spot of maize.J Shenyang Agric Univ， 2000， 31∶ 418-422 （in Chinese with English abstract）

［2］Wang J， Levy M， Dunkle L D.Sibling species of Cercospora associated with gray leaf spot of maize.Phytopathology， 1998， 88∶1269-1275

［3］Ward J M J， Stromberg E L， Nowell D C， Nutter F W J.Gray leaf spot∶a disease of global importance in maize production.Plant Dis， 1999， 83∶ 884-895

［4］李世強(qiáng)， 陳威， 譚靜， 王坤， 番興明.玉米灰斑病的抗性機(jī)理研究進(jìn)展.生物技術(shù)進(jìn)展， 2011， （1）∶ 112-115 Li S Q， Chen W， Tan J， Wang K， Fan X M.The progress of resistance mechanism on maize gray leaf spot.Curr Biotechnol， 2011，（1）∶ 112-115 （in Chinese with English abstract）

［5］曹國(guó)輝.玉米灰斑病及抗性研究.玉米科學(xué)， 2009， 17∶152-155 Cao G H.The research advance on resistance to grey leaf spot in maize.J Maize Sci， 2009， 17∶ 152-155 （in Chinese with English abstract）

［6］呂香玲， 李新海， 陳陽(yáng)， 史振聲， 李鳳海， 傅俊范.玉米種質(zhì)抗灰斑病鑒定與評(píng)價(jià).玉米科學(xué)， 2011， 19（6）∶ 125-128 Lü X L， Li X H， Chen Y， Shi Z S， Li F H， Fu J F.Evaluation and identification of resistance to gray leaf spot （gls） in maize germplasm.J Maize Sci， 2011， 19（6）∶ 125-128 （in Chinese with English abstract）

［7］Zhang Y， Xu L， Fan X M， Tan J， Chen W.QTL mapping of resistance to gray leaf spot in maize.Theor Appl Genet， 2012， 125∶1797-1808

［8］Gordon S G， Bartsch M， Matthies I， Gevers H O， Lipps P E.Linkage of molecular markers to Cercosporazeae-maydis resistance in maize.Crop Sci， 2004， 44∶ 628-636

［9］Lehmensiek A， Esterhuizen A M， Staden D V， Nelson S W， Retief A E.Genetic mapping of gray leaf spot（GLS） resistance genes in maize.Theor Appl Genet， 2001， 103∶ 797-803

［10］Danson J， Lagat M， Kimani M， Kuria A.Quantitative trait loci （QTLs） for resistance to gray leaf spot and common rust diseases of maize.Afr J Biotechnol， 2008， 7∶ 3247-3254

［11］Clements M J， Dudley J W， White D G.Quantitative trait loci associated with resistance to gray leaf spot of corn.Phytopathology， 2000， 90∶ 1018-1025

［12］Saghai Maroof M A， Yue Y G， Xiang Z X， Stromberg E L，Rufener G K.Identification of quantitative trait loci controlling resistance to gray leaf spot disease in maize.Theor Appl Genet，1996， 93∶ 539-546

［13］Veiga A D， Pinho R G V， Resende L V.Quantitative trait loci associated with resistance to gray leaf spot and grain yield in corn.Ciênc Agrotec， 2012， 36∶ 31-38

［14］Pozar G， Butruille D， Silva H D， McCuddin Z P， Penna J C V.Mapping and validation of quantitative trait loci for resistance to Cercosporazeae-maydis infection in tropical maize.Theor Appl Genet， 2009， 118∶ 553-564

［15］Chung C L， Poland J， Kump K， Benson J， Longfellow J， Walsh E，Balint K P， Nelson R.Targeted discovery of quantitative trait loci for resistance to northern leaf blight and other diseases of maize.Theor Appl Genet， 2011， 123∶ 307-326

［16］Zwonitzer J C， Coles N D， Krakowsky M D， Arellano C， Holland J B， McMullen M D， Pratt R C， Balint-Kurti P J.Mapping resistance quantitative trait Loci for three foliar diseases in a maize recombinant inbred line population-evidence for multiple disease resistance？ Phytopathology， 2010， 100∶ 72-79

［17］Asea G， Vivek B S， Lipps P E， Pratt R C.Genetic gain and cost efficiency of marker-assisted selection of maize for improved resistance to multiple foliar pathogens.Mol Breed， 2012， 29∶515-527

［18］Balint K P J， Wisser R， Zwonitzer J C.Use of an advanced intercross line population for precise mapping of quantitative trait loci for gray leaf spot resistance in maize.Crop Sci， 2008， 48∶1696-1704

［19］Bubeck D M， Goodman M M， Beavis W D， Grant D.Quantitative trait loci controlling resistance to gray leaf spot in maize.Crop Sci， 1993， 33∶ 838-847

［20］Berger D K， Carstens M， Korsman J N， Middleton F， Kloppers F J， Tongoona P， Myburg A A.Mapping QTL conferring resistance in maize to gray leaf spot disease caused by Cercospora zeina.BMC Genet， 2014， 15∶ 60， DOI∶ 10.1186/1471-2156-15-60

［21］Lee M， Sharopova N， Beavis W D， Grant D， Katt M， Blair D，Hallauer A.Expanding the genetic map of maize with the intermated B73 × Mo17 （IBM） population.Plant Mol Biol， 2002， 48∶453-461

［22］Chardon F， Virlon B， Moreau L， Falque M， Joets J， Decousset L，Murigneux A， Charcosset A.Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome.Genetics， 2004， 168∶ 2169-2185

［23］Truntzler M， Barrière Y， Sawkins M C， Lespinasse D， Betran J，Charcosset A， Moreau L.Meta-analysis of QTL involved in silage quality of maize and comparison with the position of candidate genes.Theor Appl Genet， 2010， 121∶ 1465-1482

［24］Lü X L， Li X H， Xie C X， Hao Z F， Ji H L， Shi L Y， Zhang S H.Comparative QTL mapping of resistance to sugarcane mosaic virus in maize based on bioinformatics.Front Agric China， 2008，24∶ 365-371

［25］Bernacchi D， Beck-Bunn T， Eshed Y， Lopez J， Petiard V， Uhlig J，Zamir D， Tanksley S.Identification of QTLs for traits of agronomic importance from Lycopersicon hirsutum.Theor Appl Genet， 1998， 97∶ 381-397

［26］徐建龍， 薛慶中， 羅利軍， 黎志康.近等基因?qū)胂刀ㄎ凰究沟厩?shù)量性狀位點(diǎn)的研究初報(bào).浙江農(nóng)業(yè)學(xué)報(bào)， 2002，14（1）∶ 14-19 Xu J L， Xue Q Z， Luo L J， Li Z K.Preliminary report on quantitative trait loci mapping of false smut resistance using near-isogenic introgression lines in rice.Acta Agric Zhejiangensis，2002， 14（1）∶ 14-19 （in Chinese with English abstract）

［27］Li Z K， Fu B Y， Gao Y M， Xu J L， Ali J， Lafitte H R， Jiang Y Z，Domingo R J， Vijayakumar C H M， Maghirang R， Zheng T Q，Zhu L H.Genome-wide introgression lines and their use in genetic and molecular dissection of complex phenotypes in rice （Oryza sativa L.）.Plant Mol Biol， 2005， 59∶ 33-52

［28］李芳， 程立銳， 許美容， 徐建龍， 黎志康.利用品質(zhì)性狀的回交選擇導(dǎo)入系挖掘水稻抗紋枯病 QTL.作物學(xué)報(bào)， 2009， 35∶1729-1737 Li F， Cheng L R， Xu M R， Xu J L， Li Z K.QTL mining for sheath blight resistance using the backcross selected introgression lines for grain quality in rice.Acta Agron Sin， 2009， 35∶ 1729-1737 （in Chinese with English abstract）

［29］Xie X B， Jin F X， Song M H， Suh J P， Hwang H G， Kim Y G，McCouch S R， Ahn S N.Fine mapping of a yield-enhancing QTL cluster associated with transgressive variation in an Oryza sativa × O.rufipogon cross.Theor Appl Genet， 2008， 116∶ 613-622

［30］Ho J C， McCouch S R， Smith M E.Improvement of hybrid yield by advanced backcross QTL analysis in elite maize.Theor Appl Genet， 2002， 105∶ 440-448

［31］呂香玲， 李新海， 郝轉(zhuǎn)芳， 吉海蓮， 史利玉， 張世煌.基于近等基因?qū)胂蛋l(fā)掘玉米抗甘蔗花葉病毒主效基因.玉米科學(xué)，2007， 15（3）∶ 9-14 Lü X L， Li X H， Hao Z F， Ji H L， Shi L Y， Zhang S H.Identification of major QTL for resistance to sugarcane mosaic virus in maize based on introgression lines analysis.J Maize Sci， 2007，15（3）∶ 9-14 （in Chinese with English abstract）

［32］Li Y， Liu C， Shi Y， Song Y， Wang T， Li Y.Genetic location and evaluation of chromosomal segments for drought tolerance at flowering stage in maize using selected backcross populations.Plant Breed， 2009， 128∶ 342-347

［33］Salvi S， Corneti S， Bellotti M， Carraro N， Sanguineti M C， Castelletti S， Tuberosa R.Genetic dissection of maize phenology using an intraspecific introgression library.BMC Plant Biol， 2011，11∶ 4

［34］Pea G， Aung H H， Frascaroli E， Landi P， Pè M E.Extensive genomic characterization of a set of near-isogenic lines for heterotic QTL in maize （Zea mays L.）.BMC Genomics， 2013， 14∶ 61

［35］Teng F， Zhai L， Liu R， Bai W， Wang L， Huo D， Tao Y， Zheng Y，Zhang Z.ZmGA3ox2， a candidate gene for a major QTL， qPH3.1，for plant height in maize.Plant J， 2013， 73∶ 405-416

［36］董懷玉， 姜鈺， 王麗娟， 徐秀德.玉米種質(zhì)資源抗灰斑病鑒定與評(píng)價(jià).植物遺傳資源學(xué)報(bào)， 2005， 6∶ 441-443 Dong H Y， Jiang Y， Wang L J， Xu X D.Evaluation on maize germplasm resources for resistance to gray leaf spot.J Plant Genet Resour， 2005， 6∶ 441-443 （in Chinese with English abstract）

［37］CIMMYT.Applied Molecular Genetics Laboratory.Laboratory Protocol.Mexico∶ CIMMYT， 1998

［38］Gordon S G， Lipps P E， Pratt R C.Heritability and components of resistance to Cercospora zeae-maydis derived from maize inbred VO613Y.Phytopathology， 2006， 96∶ 593-598

［39］Benson J M， Poland J A， Benson B M， Stromberg E L， Nelson R J.Resistance to gray leaf spot of maize∶ genetic architecture and mechanisms elucidated through nested association mapping and near-isogenic line analysis.PLoS Genet， 2015， 11∶e1005045

DOI：10.3724/SP.J.1006.2016.00758

*通訊作者（

Corresponding author）∶ 呂香玲， E-mail∶ lvxiangling521@126.com， Tel∶ 13664105185

收稿日期Received（）∶ 2015-10-15； Accepted（接受日期）∶ 2016-03-02； Published online（網(wǎng)絡(luò)出版日期）∶ 2016-03-14.

Meta-analysis and Validation of QTL for Resistance to Gray Leaf Spot in Maize

YAN Wei， LI Yuan， SONG Mao-Xing， ZHANG Kuang-Ye， SUN Ming-Ze， QU Hui， LI Feng-Hai， ZHONG Xue-Mei， ZHU Min， DU Wan-Li， and Lü Xiang-Ling*
Special Maize Institute， Shenyang Agricultural University， Shenyang 110866， China

Abstract：Gray leaf spot （GLS） is one of the most severe leaf diseases of maize worldwide.The breeding progress for resistance to GLS has been seriously hindered by less knowledge about QTL number， QTL intervals and mechanism of GLS.Based on meta-analysis， we conformed 13 consensus QTL regions from 14 articles on resistance to gray leaf spot of maize.One backcross with inbred line 81162 as recurrent parent and inbred line CN165 as donor parent， was used to test the consensus of those 13 QTL regions on the basic of linkage disequilibrium getting more than 20 partial separated loci.High level of partial separation indicated that there were QTLs with high effects of resistance to GLS on chromosome 1 and chromosome 4.On chromosome 1， the donor genes closed to markers of umc2227， bnlg1832， umc1243， umc2025， umc1515， umc1297， and umc1461 showed the frequency over 50%.Therefore， we inferred that there were a few of highly linked QTLs on chromosome 1.The resistant QTL on chromosome 4 was located between markers bnlg2291 and umc1194.Consequently， this study could lay a foundation for the QTL fine mapping on chromosome 1 and chromosome 4 in the donor parent CN165.

Keywords：Gray leaf spot； Meta-analysis； Consensus QTLs； Backcross introgression lines； Major QTL