姚琦馥，魏佳泰，李水琴，劉 航，馬 建

(1.銅仁學院農(nóng)林工程與規(guī)劃學院/貴州省梵凈山地區(qū)生物多樣性保護與利用重點實驗室，貴州銅仁 554300； 2.四川農(nóng)業(yè)大學小麥研究所，四川成都 611130)

大麥(HordeumvulgareL.,2n=14,HH)是世界上僅次于小麥、玉米和水稻的第四大重要谷類作物[1]。較小麥而言，大麥產(chǎn)量低、面筋含量低，且易倒伏[2]，但大麥卻含有種類豐富的營養(yǎng)物質(zhì)，如高蛋白、高膳食纖維、高維生素、低脂肪、低糖等[3-4]，因此大麥在保健方面具有很大的潛在價值。在種植面積有限的情況下，提高大麥單位面積產(chǎn)量意義重大。

植物葉片的形狀、大小可以調(diào)節(jié)植物的生長發(fā)育，如株高、粒長、粒寬、分蘗、產(chǎn)量等[5-6]，旗葉長和寬是旗葉大小的關鍵組成因素[7-9]，是谷類作物提高產(chǎn)量的兩個重要農(nóng)藝性狀，對產(chǎn)量具有重要作用[10-11]。旗葉長和寬均屬于數(shù)量性狀(quantitative trait locus,QTL)，通常受多個基因控制，受環(huán)境影響也較大[12-14]。與旗葉長和寬相關的QTL已經(jīng)在小麥、水稻等谷類作物中被廣泛報道[15-17]。Farooq等[18]利用水稻材料IR64及其基因導入系對旗葉長和寬進行了QTL定位，在2號和4號染色體共定位到2個控制旗葉長的QTL；在1號和4號染色體上共定位到2個控制旗葉寬的QTL。Bing等[19]利用水稻品種珍汕97和IRAT109構建的包含180個株系的F9代重組自交系群體，在2、3、4、5、6、7和10號染色體上共檢測到8個控制旗葉長的QTL，在4、5和6號染色體上共檢測到5個控制旗葉寬的QTL。Ma等[11]于11個不同生態(tài)環(huán)境下在小麥2D和5B染色體上共檢測到2個控制旗葉長的位點QFlw.sicau-2D和QFll.sicau-5B.3，在2D染色體上檢測到1個控制旗葉寬的位點QFlw.sicau-2D，可解釋4.09%～35.68%的表型變異。Liu等[6]對大麥旗葉形態(tài)相關的QTL進行了鑒定，但目前尚未見有人發(fā)現(xiàn)穩(wěn)定表達的主效QTL，嚴重影響了通過調(diào)節(jié)旗葉形態(tài)來優(yōu)化大麥株型結構，進而來實現(xiàn)大麥產(chǎn)量提升的目標。因此，本研究利用DArT標記對1個大麥重組自交系(RIL)群體進行分析，并基于多環(huán)境對該群體進行旗葉形態(tài)鑒定，以期獲得控制旗葉長和寬的穩(wěn)定表達的主效位點，為解析大麥旗葉形態(tài)遺傳機制和分子輔助育種奠定基礎。

1 材料與方法

1.1 試驗材料

供試材料為1個通過單粒傳方法獲得的F10代大麥RIL群體。該群體包含125個株系，其親本為經(jīng)過多年多點鑒定表現(xiàn)穩(wěn)定的J36528和BMJ89材料，均由四川農(nóng)業(yè)大學小麥研究所保存。母本J36528旗葉較父本BMJ89長，而父本BMJ89旗葉較母本J36528寬。

1.2 試驗方法

將RIL群體和親本分別于2016、2017、2018和2019年種植在四川農(nóng)業(yè)大學崇州試驗基地(分別用2016CZ、2017CZ、2018CZ和2019CZ表示)，于2018和2019年也種植在四川農(nóng)業(yè)大學雅安試驗田(分別用2018YA和2019YA表示)進行表型測定。試驗田按照隨機區(qū)組設計，包含2個重復。每個株系播種2行，行距0.3 m，行長1.5 m，株距0.1 m。試驗材料均按單粒播種法進行播種，田間管理方式按常規(guī)管理方法進行。

1.3 旗葉長和寬表型數(shù)據(jù)的測定

于小麥開花期，每個株系選取長勢良好且均一的單株用直尺進行旗葉長和寬的測定(不取每一行兩邊的單株，排除邊際效應)。每個株系選取6個單株進行測量，同時計算出每個株系旗葉長和寬的平均值，用于后續(xù)分析。

1.4 遺傳圖譜的構建和QTL定位

按照Wenzl等[20]的方法利用DArT標記(https://www.diversityarrays.com/)掃描J365283/BMJ89 RIL群體125個株系中的92個單株和親本，并利用已有SSR標記[21-23]對92個株系的群體進行分型[24]。對于存在共分離標記的位點，僅保留1個標記進行圖譜構建。利用JointMap 4.0構建連鎖圖譜。利用IciMapping 4.0軟件進行QTL檢測，找到LOD>2.5的位點，同時計算出每個QTL的表型貢獻率和加性效應。將在2個環(huán)境以上檢測獲得且表型貢獻率大于10%的QTL視為主效穩(wěn)定表達的位點。用Triticeae Multi-omics Center網(wǎng)站(http://202.194.139.32/blast/viroblast.php)，將本研究鑒定到的以及前人檢測到的QTL側翼標記序列與大麥基因型‘Morex’參考基因組[25]進行比對(默認參數(shù))，獲得QTL對應的物理位置。

1.5 數(shù)據(jù)分析

使用Excel 2019整理數(shù)據(jù)；使用SAS 9.4軟件計算旗葉長和寬的最佳線性無偏預測(best linear unbiased prediction，BLUP)值和遺傳力(H2)；使用SPSS 17.0軟件分析旗葉長和寬在不同環(huán)境及BLUP值的相關性；使用Origin 9.3軟件繪制旗葉長和寬的頻率分布圖。

2 結果與分析

2.1 親本及RIL群體旗葉長和寬的表型

對6個環(huán)境表型數(shù)據(jù)處理后發(fā)現(xiàn)，除了缺失數(shù)據(jù)外，J36528的旗葉長均極顯著高于BMJ89，J36528的旗葉寬均極顯著低于BMJ89(表1)。RIL群體在不同環(huán)境的旗葉長和旗葉寬均出現(xiàn)超親分離現(xiàn)象，且旗葉長和寬均具有較高的遺傳力，分別為0.93和0.82(表1)，表明這2個性狀主要受遺傳因子影響，受環(huán)境影響較小。旗葉長和寬在每個環(huán)境下均呈現(xiàn)近似正態(tài)分布，符合數(shù)量遺傳性狀的特點，適合進行QTL定位(圖1)。

圖1 不同環(huán)境下旗葉長和寬的表型頻率分布Fig.1 Frequency distribution of flag leaf length and width of the RIL population in different environments

表1 不同環(huán)境下親本及RIL群體旗葉長和寬的表型統(tǒng)計結果Table 1 Phenotypic statistics of flag leaf length and width of the parents and RIL population in different environments

2.2 旗葉長和寬在不同環(huán)境下的相關性

從表2可以看出，旗葉長和寬在不同環(huán)境及BLUP值之間均顯著相關，旗葉長的相關系數(shù)介于0.975～0.997之間，而旗葉寬的相關系數(shù)介于0.955～0.994之間，進一步說明旗葉長寬受環(huán)境影響較小。

表2 旗葉長和寬在不同環(huán)境下的相關性Table 2 Correlation coefficients of flag leaf length and width in different environments

2.3 遺傳圖譜的構建

基于447個DArT和8個SSR標記構建大麥群體的遺傳連鎖圖譜，分析結果(表3)顯示，染色體上標記數(shù)量介于26～99之間，平均每條染色體有65個標記；圖譜包含7個遺傳連鎖群，總長1 034.96 cM，每條染色體圖譜長度介于94.96 cM ～218.80 cM，平均每條染色體圖譜長度為147.85 cM。標記密度介于0.96～6.70 cM/標記之間，平均標記密度為2.27 cM/標記。

表3 構建的大麥遺傳圖譜中標記在染色體上的分布Table 3 Distribution of markers in the constructed barley genetic map

2.4 旗葉長和寬的QTL檢測結果

基于遺傳連鎖圖譜和測量獲得的4年6個環(huán)境的表型數(shù)據(jù)，通過IciMapping 4.0軟件共檢測到12個控制旗葉長和寬的QTL(表4)。其中，檢測到7個與旗葉長相關的QTL，分布在2H、5H、6H和7H染色體上。其中，5H染色體上的QFll.sicau-JB-5H.2在3個環(huán)境及BLUP值中被檢測到，可解釋12.02%～19.95%的表型變異，LOD值為3.14～3.61，為穩(wěn)定表達的主效QTL，加性效應為正值，表明加性效應來自母本J36528。6H染色體上的QFll.sicau-JB-6H.1在2個環(huán)境及BLUP值中被檢測到，可解釋 16.69%～16.99%的表型變異，LOD值為 2.89～4.49，為穩(wěn)定表達的主效QTL，加性效應為負值，表明其加性效應來自父本BMJ89。QFll.sicau-JB-2H.2和QFll.sicau-JB-7H在2個環(huán)境中被檢測到，加性效應分別來自親本J36528和BMJ89可分別解釋9.81%～10.31%和 14.28%～ 14.45%的表型變異。其余3個QTL只在單一的環(huán)境中被檢測到，可解釋13.79%～ 28.84%的表型變異。

表4 基于大麥RIL群體檢測到的與旗葉長和寬相關的QTLTable 4 QTLsfor flag leaf length and width detected in the J36528/BMJ89 (JM) population

檢測到5個與旗葉寬相關的QTL，分布在2H、4H、5H、6H和7H染色體上。其中QFlw.sicau-JB-4H在4個環(huán)境及BLUP值中被檢測到，可解釋15.03%～24.49%的表型變異，LOD值為2.59～6.91，為穩(wěn)定表達的主效QTL，加性效應為負值，表明其加性效應來自父本BMJ89。其余4個QTL均只在單一環(huán)境中被檢測到，可解釋了14.68%～32.66%的表型變異。

為了進一步揭示3個穩(wěn)定表達的主效位點(與旗葉長相關的QFll.sicau-JB-5H.2和QFll.sicau-JB-6H.1以及與旗葉寬相關的QFlw.sicau-JB-4H)的遺傳效應，根據(jù)以上3個位點的側翼標記，將RIL群體中的92個株系分為2種類型，一種是攜帶相應QTL增效位點的株系，另一種是攜帶相應QTL減效位點的株系。進一步對2種類型進行檢測，結果(圖2)表明，對于3個QTL而言，含有各自QTL增效位點株系的旗葉長或旗葉寬在各個環(huán)境下都顯著或者極顯著高于含有減效位點的株系(除了QFll.sicau-JB-5H.2位點的2019YA環(huán)境)。

*和**分別表示各環(huán)境下親本J36528和BMJ89之間在0.05和0.01水平上差異顯著。白色框代表帶有BMJ89等位基因的純合株系，灰色框代表帶有J36528等位基因的純合株系。＊代表異常值；□代表平均數(shù)。*and**indicate significant differences between J36528 and BMJ89 at the 0.05 and 0.01 levels,respectively.White box represents homozygous lines with BMJ89 allele,and gray box represents homozygous lines with J36528 allele.＊ represents outlier;□ represents mean value.圖2 旗葉長和寬主效QTL的遺傳效應Fig.2 Genetic effect of major QTLs for flag leaf length and width

2.5 與旗葉長相關的2個主效且穩(wěn)定表達QTL的聚合效應

本研究檢測到2個控制大麥旗葉長的主效且穩(wěn)定表達的QTL，即QFll.sicau-JB-5H.2和QFll.sicau-JB-6H.1。為了進一步分析這2個QTL的聚合效應，本研究根據(jù)QTL的側翼標記對含有這2個QTL不同效應位點的株系進行分析。結果(圖3)發(fā)現(xiàn)，含有2個QTL增效位點的株系(20個)的旗葉長顯著或極顯著高于不含任何位點株系(12個)以及僅含有QFll.sicau-5H.2或QFll.sicau-6H.1增效位點的株系，增幅分別為73.28%、14.33%和20.60%。僅含有QFll.sicau-5H.2或QFll.sicau-6H.1增效位點株系的旗葉長也極顯著高于不含任何位點的株系，增幅分別為51.58%和43.68%。但含有QFll.sicau-5H.2和QFll.sicau-6H.1增效位點的株系之間旗葉長無顯著差異。

+和-分別代表攜帶和不攜帶目標QTL增效位點的株系。圖中百分比數(shù)據(jù)表示兩個株系之間旗葉長的變幅。*和**分別表示株系之間在0.05和0.01水平上差異顯著。+ and - represent lines with and without the positive alleles of the target quantitative trait loci (QTL),respectively.The percentage data in the figure represent the variation of flag leaf length between two lines. * and ** indicate significant differences between lines at 0.05 and 0.01 levels,respectively.圖3 兩個旗葉長主效QTL的聚合效應Fig.3 Multiplying effect of the two major QTLs on flag leaf length of the RIL population

2.6 主效QTL在大麥參考基因組的物理位置

將主效位點QFlw.sicau-JB-4H、QFll.sicau-JB-5H.2和QFll.sicau-JB-6H.1的側翼標記序列錨定到大麥參考基因組上，發(fā)現(xiàn)QFlw.sicau-JB-4H位于大麥4H染色體長臂584.05～612.36 Mbp區(qū)間，QFll.sicau-JB-5H.2位于5H染色體長臂379.10～540.62 Mbp區(qū)間，QFll.sicau-JB-6H.1位于6H染色體短臂32.39～ 33.59 Mbp區(qū)間(圖4)。

圖4 大麥旗葉長和寬3個主效QTL的遺傳圖譜和物理圖譜Fig.4 Genetic maps and physical maps of the three major QTLs for flag leaf length and width

3 討 論

3.1 加快推進挖掘大麥旗葉形態(tài)調(diào)控位點

大麥作為第四大谷類作物，其產(chǎn)量增長對人們的生活有著重要意義。旗葉對于谷物產(chǎn)量的積極作用已有很多報道，例如Sakamoto等[9]發(fā)現(xiàn)，旗葉較直的水稻新品種可以增加光合作用所需的光捕捉量，進而增加籽粒灌漿所需的氮素儲存量，最終提高籽粒產(chǎn)量；Ma等[11]發(fā)現(xiàn)，旗葉大小和形態(tài)是小麥“理想型”的兩個重要因素，適當?shù)钠烊~長和寬可以增加小麥產(chǎn)量；Simon等[12]研究也發(fā)現(xiàn)，小麥的旗葉角度和旗葉面積是決定產(chǎn)量的重要性狀。除了與產(chǎn)量相關的報道外，在小麥、水稻和玉米中，還有很多關于旗葉形態(tài)調(diào)控位點鑒定的報道[13-15,26]。而在大麥中，前人對于大麥旗葉的研究多集中在光合特性等生理生化指標方面[27-29]，而鮮有關于大麥旗葉相關性狀位點的報道，到目前為止，對于大麥旗葉遺傳調(diào)控機制的認識還很有限。大麥作為重要糧食作物和飼料加工作物，提升其產(chǎn)量和品質(zhì)非常必要，本研究獲得的大麥旗葉相關性狀定位結果，為進一步推進大麥旗葉形態(tài)的遺傳研究和精細定位奠定了基礎，也為產(chǎn)量的提升提供了更多的可能。

3.2 大麥RIL群體旗葉受環(huán)境的影響較小

前人研究發(fā)現(xiàn)，旗葉形態(tài)受環(huán)境影響較大[11,30]，而本研究發(fā)現(xiàn)，大麥旗葉長和寬在不同環(huán)境下極顯著相關，即受到環(huán)境的影響較小，原因還有待進一步探究。本研究也檢測到3個主效且穩(wěn)定表達的QTL，說明這些位點在實際育種中具有潛在利用價值。

3.3 與旗葉形態(tài)相關的穩(wěn)定表達的主效QTL

本研究對RIL群體125個株系中的92個單株進行了標記掃描，構建的圖譜準確可信，基于該圖譜成功鑒定到3個主效且在多個環(huán)境穩(wěn)定表達的旗葉長和寬QTL。迄今為止，關于大麥旗葉長和寬位點鑒定的報道極少。Liu等[6]在4H染色體長臂上鑒定獲得1個控制大麥旗葉寬的位點qFLW4-18，經(jīng)過比對發(fā)現(xiàn)，其位于4H染色體601.76 Mbp處，與本文鑒定到的旗葉寬位點QFlw.sicau-JB-4H位于同一個區(qū)間，因此，推測兩者可能為等位基因。而Liu等[6]鑒定獲得的大麥旗葉長QTL主要分布于2H和3H染色體上，在其他染色體未見報道。因此，推測QFll.sicau-JB-5H.2和QFll.sicau-JB-6H.1可能為控制大麥旗葉長的新位點。值得注意的是，在小麥第5部分同源群染色體長臂已檢測到主效且穩(wěn)定表達的旗葉長位點，如QFll.sicau-5B.3[11]、QFll.sau-SY-5B[10]和QFll.cau-5B[31]。由于小麥和大麥染色體具有一定的線性相關性，因而，QFll.sicau-JB-5H.2有可能與這些位點具有同源性。此外，小麥春化基因Vrn-B1也位于5B染色體長臂上[32-33]。這些位于小麥和大麥第5部分同源群長臂的旗葉相關QTL與春化基因的相關性有待進一步揭示。

3.4 與旗葉長相關的主效且穩(wěn)定表達QTL的聚合效應

對控制同一性狀的不同位點進行聚合，通常會起到一定的疊加效應[34]。本研究對控制旗葉長的主效位點QFll.sicau-JB-5H.2和QFll.sicau-JB-6H.1進行聚合分析后發(fā)現(xiàn)，二者疊加后旗葉長極顯著增加。這為大麥育種中通過聚合不同位點來改良旗葉性狀提供了理論依據(jù)。