趙艷艷，袁亞培，梁 雪，宮曉平，吳春紅，周秀文，郭 營，趙 巖，李斯深，孔凡美*

（1 作物生物學國家重點實驗室，土肥資源高效利用國家工程實驗室，山東農(nóng)業(yè)大學，山東泰安 271018；2 山東省菏澤市曹縣農(nóng)業(yè)局，山東曹縣 274400）

不同磷、鉀處理小麥苗期氮營養(yǎng)性狀的QTL分析

趙艷艷1，袁亞培1，梁 雪2，宮曉平1，吳春紅1，周秀文1，郭 營1，趙 巖1，李斯深1，孔凡美1*

（1 作物生物學國家重點實驗室，土肥資源高效利用國家工程實驗室，山東農(nóng)業(yè)大學，山東泰安 271018；2 山東省菏澤市曹縣農(nóng)業(yè)局，山東曹縣 274400）

【目的】對不同濃度磷、鉀處理下小麥苗期氮養(yǎng)分效率相關性狀進行 QTL 分析，以深入理解磷、鉀與氮養(yǎng)分效率的相互關系，為氮營養(yǎng)相關性狀的圖位克隆及分子標記輔助選擇育種奠定基礎。 【方法】采用苗期液培試驗，以“川 35050 × 山農(nóng) 483”組合衍生的小麥重組自交系群體 (131 個株系) 為研究材料，設置了中磷中鉀(MPMK)、高磷 (HP)、低磷 1 (LP1)、低磷 2 (LP2)、低磷 3 (LP3)，高鉀 (HK)、低鉀 1 (LK1)、低鉀 2 (LK2)、低鉀 3 (LK3) 共 9 個處理，對不同磷、鉀處理下的氮養(yǎng)分效率相關性狀進行研究，并結(jié)合分子標記遺傳圖譜，從整個基因組水平對與小麥苗期氮養(yǎng)分效率相關的 10 個性狀進行 QTL 定位及遺傳分析。 【結(jié)果】不同處理下的10 個性狀共檢測到 137 個 QTL，位于除 3D 外的 20 條染色體上，大部分 QTL (89.05%) 僅在單一處理下被定位到，有 3 個 QTL (QRnue-1A.2、QSnue-1A.1 和 QTnue-1A.1) 可在至少 4 個處理中被檢測到，有 5 個 QTL (QRnue-1A.1、QTnue-1A.1、QSnc-4A、QRnc-6A.3 和 QSnue-6B) 可同時在低磷和低鉀環(huán)境中被檢測到。本研究還檢測到至少包含 3 個以上 QTL 的 QTL 簇 17 個，分別位于 1A、1B、2B、2D、3A、3B、4A、4B、5D、6A、6B、6D 和 7A 染色體上，共涉及 66 個 QTL，占 QTL 總數(shù)的 48.18%。其中，有 5 個 QTL簇僅與特定磷、鉀處理有關，大多數(shù) QTL 簇均同時定位了不同磷、鉀處理的不同性狀，許多 QTL 簇位點還與前人定位的生物量、產(chǎn)量及其他養(yǎng)分有關。 【結(jié)論】磷、鉀的供應能夠顯著影響小麥苗期對氮素的吸收利用及其相關 QTL 的表達。影響苗期小麥氮養(yǎng)分效率相關性狀的 QTL 大多數(shù)僅在特定處理下被檢測到，但大多數(shù) QTL會形成 QTL 簇，構(gòu)成了控制氮養(yǎng)分效率的 QTL 熱點，許多熱點區(qū)域也與前人定位的許多成株期性狀如生物量、產(chǎn)量及其他養(yǎng)分效率有關，這些 QTL/基因密集區(qū)域及其特點的發(fā)現(xiàn)，為我們深入理解小麥氮養(yǎng)分效率的遺傳控制特點及其與磷、鉀養(yǎng)分供應的關系提供了新的視角，也為這些重要位點的克隆及其應用提供數(shù)據(jù)支持。

小麥；重組自交系群體；QTL；氮養(yǎng)分效率

氮素是作物生長發(fā)育必需的大量元素之一。我國部分高產(chǎn)攻關田小麥氮肥用量達 551 kg/hm2，而氮肥利用率僅為 23%～28%[1]。隨著全球人口的增加，為增加作物的產(chǎn)量，大量的氮肥資源投入農(nóng)田，不僅降低了肥料利用率[2]，造成了資源的浪費，還引發(fā)了嚴重的環(huán)境問題[3–4]。因此，培育高產(chǎn)、養(yǎng)分高效利用的新品種成為“第二次綠色革命”的重要目標。深入探討氮效率相關性狀的遺傳基礎是小麥氮養(yǎng)分效率遺傳改良的前提[5]。

氮營養(yǎng)相關性狀是典型的數(shù)量性狀，數(shù)量性狀基因座位 (quantitative trait locus, QTL) 分析可以把復雜的數(shù)量性狀分解成不同的 QTL，為復雜性狀遺傳控制的研究提供了一種有效的新方法[6]。到目前為止，與小麥氮素營養(yǎng)有關的 QTL 的分析主要集中在苗期或成熟期，涉及低氮或高氮處理。Quarrie 等[7]在低氮脅迫的環(huán)境下檢測到與產(chǎn)量性狀相關的QTL，主要分布在 4A、4B、6A、7A 和 7B 染色體上。An 等[8]在大田環(huán)境下共檢測到 17 個與氮吸收效率有關的 QTL，高氮和低氮環(huán)境下分別為 8 個和 9 個。2006 年，Laperche 等[9]在高氮和低氮水平下，檢測到與籽 粒產(chǎn)量、籽粒蛋白和氮含量相關的 QTL，在 4B和 5A 染色體上檢測到 2 個與根干重相關的 QTL，此外，在 2B 染色體上定位了 1 個與總干重有關的QTL。2007 年，Laperche 等[10]對產(chǎn)量、千粒重、收獲指數(shù)等性狀進行 QTL 分析，在高氮和低氮環(huán)境下檢測到的 QTL 數(shù)分別為 67 個和 51 個，同時在兩個環(huán)境中檢測到的 QTL 為 85 個。由于 QTL 受試驗環(huán)境的影響很大，這些 QTL 的穩(wěn)定性需要進一步驗證。特別是磷和鉀與氮素的吸收利用效率存在復雜的交互效應[11]，氮、磷、鉀的吸收和利用可能受共同的 QTL 控制[12]。因此，研究不同磷、鉀處理環(huán)境下的氮效率相關 QTL 對于深入理解氮、磷、鉀相互關系及氮營養(yǎng)性狀的 QTL 環(huán)境穩(wěn)定性具有重要意義，也為氮營養(yǎng)性狀的圖位克隆及分子標記輔助選擇育種 (marker assistant selection, MAS) 奠定基礎。

本研究以小麥重組自交系群體 (recombinant inbred lines, RIL) 為研究材料，研究了不同濃度磷、鉀處理對小麥苗期氮養(yǎng)分效率相關性狀的影響并對所有氮養(yǎng)分效率相關性狀進行了 QTL 分析。

1 材料與方法

1.1 試驗材料

本研究以“川 35050 × 山農(nóng) 483”組合衍生的小麥重組自交系群體 (recombinant inbred lines, RIL)(131個株系) 為研究材料。川 35050 是中國西南冬麥區(qū)的小麥品系，穗大、強筋；山農(nóng) 483 是黃淮麥區(qū)的小麥品系，中弱筋，豐產(chǎn)性好，該品系來自小麥種質(zhì)和小麥育種的骨干親本“矮孟?！?，“矮孟?！痹谛←溣N中得到了廣泛使用，以其為親本已培育出十多個小麥新品種，累計推廣 3 億多畝，增產(chǎn)小麥 100 多億公斤。

1.2 試驗設計

營養(yǎng)液培養(yǎng)試驗在山東農(nóng)業(yè)大學溫室進行，采用隨機區(qū)組設計，營養(yǎng)液以 Hoagland 營養(yǎng)液[13]為基礎，根據(jù)小麥的生理特點對其進行了優(yōu)化調(diào)整[12]。試驗共設置 9 個處理，分別為中磷中鉀 (MPMK)、高磷 (HP)、低磷 1 (LP1)、低磷 2 (LP2)、低磷 3 (LP3)，高鉀 (HK)、低鉀 1 (LK1)、低鉀 2 (LK2)、低鉀 3 (LK3)，各處理重復 3 次，具體方案見表 1。

表1 不同試驗磷、鉀濃度 (mmol/L)Table 1 P and K concentrations in the treatments

挑選大小一致的小麥重組自交系 (RIL) 和親本種子，經(jīng) 10% H2O2滅菌 5 min 后用蒸餾水沖洗，充分吸脹后取出置于發(fā)芽網(wǎng)上發(fā)芽。7 天后，每個株系選取長勢一致的幼苗移栽到育苗盤上，每穴 2 株，重復 3 次，并用海綿固定好。之后將育苗盤固定在盛有 20 L 營養(yǎng)液、外側(cè)不透光 (保證根系健康生長并防止藻類的生長) 的黑色長方形塑料盒中。移苗后用蒸餾水緩苗 1 天，之后每隔 3 天更換一次營養(yǎng)液，用 0.5% 的 NaOH 調(diào)節(jié) pH 在 6.0～6.2 之間，并用充氣泵不斷向營養(yǎng)液充氣，培養(yǎng) 30 天收獲。

試驗在日光溫室進行，氣溫 7.5℃～33.9℃，平均 18.5℃；相對濕度 10.1%～91.3%，平均 51.4%；光照強度 0～69.0 Klux，平均 5.9 Klux。數(shù)據(jù)每隔10 min 由 ZDR 系列記錄儀記錄。

1.3 測定項目及方法

收獲時，先將根系在蒸餾水中沖洗 10 min 后取出，用吸水紙吸干水分。將地上部和根系分別裝入紙袋，置于 60℃ 烘箱中烘干至恒重。用萬分之一天平對單株根干重和單株地上部干重進行稱重。單株總重為單株根和地上部的干重之和。

將烘干的小麥幼苗磨碎，稱取適量加入濃 H2SO4和 H2O2消解至無色透明。用凱氏定氮法[14]測定根系氮濃度和地上部氮濃度。主要計算公式[15]：

總的氮濃度=(單株地上部干重×地上部氮濃度+單株根干重×根系氮濃度)/(單株根干重+地上部干重)

根系(地上部)氮累積量=單株根(地上部)干重×根系(地上部)氮濃度

總的氮累積量=根系氮累積量+地上部氮累積量

累積量根冠比=根系氮累積量/地上部氮累積量

根系(地上部)氮利用率=單株根(地上部)干重/根系(地上部)氮濃度

總氮利用率=單株總重/總的氮濃度

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

方差分析 (ANOVA)、最小顯著性差異 (LSD) 分析和性狀間的相關系數(shù)分析均采用 SAS 軟件。廣義遺傳率的相關數(shù)據(jù)由 SAS 軟件中的 GLM 程序得到，計算公式為：

本研究利用已構(gòu)建的遺傳圖譜[17]進行 QTL 分析，該圖譜包含 21 條染色體上的 719 個標記，圖譜全長為 4008.4 cM，標記密度為 7.15 cM/ 位點。標記類型主要為 DArT 標記 (多樣性微陣列技術)、SSR 標記 (簡單序列重復多態(tài)性)、EST-SSR 標記 (基于表達序列標簽開發(fā)的重復多態(tài)性)。采用 Windows QTL Cartographer 2.5 軟件[18]對目標性狀進行 QTL 復合區(qū)間作圖 (composite-interval mapping)。檢測步長為 1 cM，每個 QTL 的認定峰值臨界值為 LOD ≥ 3.0，并在 P ≤ 0.05[19]的情況下進行 1000 次置換檢測，每個QTL 區(qū)間入選臨界值為 LOD ≥ 2.5[12]。Stoll 等將染色體同一位置同時檢測到 2 個或 2 個以上性狀的QTL 位點定義為一個 QTL 簇[20]。本研究將染色體同一位點定位到 3 個及以上不同性狀的 QTL，定義為一個 QTL 簇[12]。

2 結(jié)果與分析

2.1 苗期小麥部分性狀的表型變異及相關分析

由表 2 可以看出，10 個性狀的基因型和處理間變異都達到顯著水平 (P ≤ 0.01)，表明遺傳因素和環(huán)境因素對氮效率相關性狀的表型變異均有顯著影響。LSD 檢測結(jié)果 (表 3) 表明，調(diào)查性狀的平均值在 9 個處理下有比較大的差異。10 個性狀中除RSNC 外，遺傳力都在 50% 以上，其中遺傳力最大為 82.5% (TNCE)，最小為 32.8% (RSNC)。所有氮效率相關性狀均表現(xiàn)出連續(xù)變異，變異系數(shù)大部分在15% 以上，最大值為 31.54% (RSNC, HP)，最小值為4.41% (TNCE, HK)。除大部分營養(yǎng)含量性狀與累積量性狀呈顯著負相關或不相關，SNC 與 RSNC 負相關或不相關外，其他各性狀表現(xiàn)出正相關或極顯著正相關 (表 4)。

表2 小麥苗期氮營養(yǎng)相關性狀的方差分析Table 2 Variance analysis (ANOVA) for the N related traits of wheat at seedling stage

2.2 苗期小麥氮效率相關性狀的QTL分析

試驗中小麥苗期 9 個處理下的 10 個性狀共檢測到 160 個 QTL (單一 QTL 單一環(huán)境)，將不同處理條件下檢測到的相同 QTL 合并后共檢測到 137 個 QTL位點，位于小麥除 3D 外的 20 條染色體上。在 7 個氮營養(yǎng)吸收相關性狀 (RNCE、SNCE、TNCE、RNC、SNC、TNC、RSNC) 和 3 個氮利用效率性狀(RNUE、SNUE、TNUE) 上分別定位到 107 個 QTL和 30 個 QTL。單個 QTL 可解釋表型變異的 6.4% ～61.5%。有 15 個 (10.95%) QTL 可以在至少兩個處理中被檢測到，貢獻率為 7.4% ～ 37%。其中，QRnue-1A.2, QSnue-1A.1 和 QTnue-1A.1 這 3 個 QTL 可在至少 4 個處理中被檢測到 (表 5)。

2.3 重要的 QTL 簇

QTL 具有成簇分布的特征，染色體同一位置同時檢測到 3 個及以上不同性狀的 QTL 形成 QTL 簇。本研究共定位到 17 個 QTL 簇 (C1～C17)，分別位于1A、 1B、 2B、 2D、 3A、 3B、 4A、 4B、 5D、6A、 6B、 6D 和 7A 染色體上，共涉及 66 個 QTL，占 QTL 總數(shù)的 48.18% (表 5)。

C1 和 C3 同時定位的 QTL 數(shù)目最多。C1 位于1A 染色體 wPt-8455-wPt-730618 區(qū)間，定位了 5 個QTL (QSnc.1、QTnc、QRnue.1、QSnue.1、QTnue.1)(表 5)。單個 QTL 能夠解釋表型變異的10.1%～18.0%。有 3 個 QTL 可以在至少 3 個處理中被定位到。該 QTL 簇中各位點加性效應均來自川35050，表現(xiàn)為正連鎖關系。

C3 位于 1A 染色體 wPt-730213-wPt-731476 區(qū)間，定位了 6 個 QTL (QRnc.2、QSnc.3、QRsnc、QRnue.2、QSnue.2、QTnue.2)(表 5)。單個 QTL 能夠

解釋表型變異的 9.1%～23.2%。其中有 3 個 QTL 可在多個處理中被定位到。該 QTL 簇中僅有 QRsnc 的加性效應來自川 35050，其他各位點加性效應均來自山農(nóng) 483。

表3 小麥苗期 RIL 群體及父母本性狀測定結(jié)果Table 3 Phenotypic performance for seedling traits of the RILs and their parents

續(xù)表3Table 3 continued

表4 不同磷、鉀處理下小麥苗期氮養(yǎng)分效率各相關性狀間的相關系數(shù)Table 4 Correlation coefficients (r) between N efficiency related traits under different P and K treatments

表5 重要的 QTL 簇Table 5 Important QTL clusters

3 討論

3.1 不同磷、鉀處理對氮效率性狀的影響

本研究中，與對照處理相比，在高磷處理下，地上部氮濃度 (SNCE)、總氮濃度 (TNCE)、地上部氮累積量 (SNC)、總氮累積量 (TNC) 均增加，表明磷供應量增加可以促進小麥對氮的吸收。而低磷處理下，根系氮濃度 (RNCE)、根系氮累積量 (RNC) 均顯著增加。有研究表明缺磷使小麥根部生長素濃度提高，促進根分生組織細胞分裂[21]，這會導致根系生長相對增加。而不同磷濃度下，地上部與根系含氮量的差異，一方面可能是因為低磷減少了幼苗對氮的吸收，另一方面，為緩解磷供應不足，植物會增加同化產(chǎn)物向根部的分配，以促進根系生長增加磷的吸收，這將導致氮被更多的用于根系生長，從而減少地上部的供應，隨著磷濃度的升高，地上部氮含量和累積量會呈現(xiàn)增加的趨勢[22]。

高鉀處理下，地上部氮濃度 (SNCE)、總氮濃度(TNCE)、地上部氮累積量 (SNC)、總氮累積量 (TNC)均增加。而低鉀處理下，根系氮濃度 (RNCE)、根系氮累積量 (RNC) 均顯著增加。譚金芳等[23]的研究也表明，隨著施鉀量的增加，冬小麥的地上部氮素累積量呈逐漸增加趨勢。其生理原因可能有兩方面：一是光反應轉(zhuǎn)化的能量有很大一部分用于氮素的同化，施鉀提高了葉片含鉀量從而提高了凈光合速率，間接促進了氮素的吸收和積累；二是鉀素提高了葉片蛋白酶的活性，從而提高了對氮素的吸收[24]。此外，鉀可以促進物質(zhì)的轉(zhuǎn)運和高分子化合物的合成[25]，缺鉀造成植株體內(nèi)氮運輸受阻，在根系內(nèi)積累較多的氮素。

低磷處理下，根系氮利用率 (RNUE) 顯著增加，這主要是由于低磷會誘導根系生長，從而間接增加了氮素的利用效率。在低磷處理下，地上部氮利用率 (SNUE)、總的氮利用率 (TNUE) 顯著下降，這與Guo 等[12]所得到的結(jié)果一致。這可能是因為低磷抑制地上部生長，增加根冠比，地上部氮素利用效率隨之下降。在低鉀處理下，根系氮利用率 (RNUE)、地上部氮利用率 (SNUE)、植株總氮利用率 (TNCE) 都有所下降，這與 Guo 等[12]在低鉀處理下得到的結(jié)果一致?？梢?，鉀素對于氮素利用效率的影響很大。以上結(jié)果表明，磷和鉀的供應能夠顯著影響小麥的氮營養(yǎng)效率及其相關性狀。

3.2 不同濃度磷、鉀處理下，小麥苗期氮效率相關性狀的QTL

近年來，關于小麥氮營養(yǎng)效率相關性狀 QTL 的研究已有一些報道。這些報道涉及不同的生長環(huán)境，如在高氮和低氮處理下，有液培試驗[8–9]、盆栽試驗[26]和大田試驗[7–8,10,27–28]。但由于養(yǎng)分之間存在復雜的相互作用，且 QTL 受環(huán)境影響較大，而探討不同養(yǎng)分處理對 QTL 影響的報道很少。本研究對不同P、K 濃度，共 9 個處理條件下的 10 個苗期性狀進行了 QTL 分析，共定位到 137 個 QTL。其中有 122個 (89.05%) 僅在一個處理中被檢測到。不同處理條件下檢測到的 QTL 數(shù)量存在很大差異，9 個處理(MPMK、HP、LP1、LP2、LP3、HK、LK1、LK2、LK3) 檢測到的 QTL (QTL-處理) 數(shù)目分別為 25、19、16、16、22、17、13、19 和 13 個。其中對照處理檢測到的 QTL 數(shù)目最多。表明磷、鉀脅迫會顯著影響 N 養(yǎng)分效率相關 QTL 的表達。許多研究也表明，不同營養(yǎng)環(huán)境下同一營養(yǎng)性狀所檢測到的QTL的數(shù)量和位置存在很大差異[8,12,26–27]。這表明控制小麥氮養(yǎng)分效率的 QTL/ 基因系統(tǒng)在不同養(yǎng)分環(huán)境條件下有很大調(diào)整，需要從整體的概念去理解這種系統(tǒng)上的調(diào)整工程。

同時，我們也發(fā)現(xiàn)，有 15 個 QTL 至少 在 2～4個處理環(huán)境中被定位到，其中，QRnue-1A.2、QSnue-1A.1 和 QTnue-1A.1 這 3 個 QTL 可在至少 4個處理環(huán)境中被定位到 (圖 1)，表明這 3 個 QTL 在不同磷、鉀濃度環(huán)境中的表達頻率較高。另外，在這 15 個相對穩(wěn)定的 QTL 中，有 5 個可在低磷和低鉀環(huán)境中同時被檢測到，分別為 QRnue-1A.1、QTnue-1A.1、QSnc-4A、QRnc-6A.3 和 QSnue-6B。說明這幾個 QTL 在低磷和低鉀脅迫的環(huán)境下均可以發(fā)揮作用?？梢?，這些 QTL 對植物在低磷和低鉀條件下調(diào)控氮的利用效率均是有效的。其調(diào)控機理還需要更深入研究。這些表達相對穩(wěn)定的 QTL 對作物育種具有重要意義，值得深入了解。

本試驗檢測到的許多 QTL 前人的試驗中也在其他營養(yǎng)環(huán)境中檢測到 (表 6)[8,12,29-33]，這些 QTL 的環(huán)境穩(wěn)定性相對較好，也是值得關注的熱點區(qū)域。

3.3 與氮效率相關的一些重要QTL簇

控制不同性狀的 QTL 會定位到染色體相同的位點上，形成 QTL 簇。本研究共檢測到 76 個染色體位點，其中有 39 個位點只定位到 1 個 QTL，其余位點可在一個或幾個環(huán)境中定位到幾個不同性狀的 QTL，形成 QTL 簇。可見，QTL 簇的形成是一個普遍現(xiàn)象，反映了一個 QTL/ 基因在一個或幾個環(huán)境中的不同作用。本試驗共檢測到與至少 3 個性狀相關的 17個 QTL 簇 (C1～C17)，由 66 個 QTL 組成，占 QTL定位總數(shù)的 48.18% (圖 1，表 5)。這表明接近一半與氮效率性狀相關的 QTL 主要集中在少數(shù)幾個區(qū)域。

圖1 小麥苗期 QTL 位點在遺傳圖譜上的位置Fig. 1 Locations of QTLs for wheat seedling traits under different P and K treatments

續(xù)圖1Figure 1 continued

續(xù)圖1Figure 1 continued

表6 本研究與前人研究所定位 QTL 位點與性狀比較Table 6 Comparison of the QTL location and trait between previous and this study

不同磷、鉀處理條件對 QTL 簇也有很大影響。有些 QTL 簇僅在特定處理下存在，例如 C5、C9、C13 只在低磷下被檢測到，C2 只在高磷下被檢測到，C16 只在不同的鉀處理中被檢測到。然而，大多數(shù)的 QTL 簇在不同磷和鉀處理中均能夠被檢測到?？梢?，有些 QTL 簇僅在植物適應低磷或者低鉀處理時控制氮的養(yǎng)分效率，而有些位點則在低磷、低鉀條件下均發(fā)揮作用，這些位點對進一步研究小麥氮營養(yǎng)均有一定的參考價值。

續(xù)表6Table 6 continued

3.4 與其他性狀有關的 QTL 簇熱點區(qū)域

通過與前人的研究結(jié)果進行對比，發(fā)現(xiàn)許多QTL 簇位點也在不同環(huán)境中檢測到與其他的一些性狀相關 (表 7)。例如，在 1D 染色體的 wmc432 標記附近，我們定位到了根系氮含量、總氮累積量性狀的 QTL；Guo 等[12]定位到了主根數(shù)、株高、單株根干重、單株地上部干重、單株總重、地上部氮累積量、總氮累積量、地上部磷累積量、總磷累積量、根鉀累積量、地上部鉀累積量、總鉀累積量、地上部氮利用率、總氮利用率、根磷利用率、地上部磷利用率、總磷利用率、根鉀利用率、地上部鉀利用率、總鉀利用率性狀的 QTL；Sun 等[29]定位到了單株總干重、單株總鮮重、單株地上部干重、單株根鮮重、單株地上部鮮重性狀的 QTL；Kong 等[30]定位到了根系鉀利用率、地上部鉀利用率、千粒重、籽粒鉀利用率性狀的 QTL；Xu 等[31]定位到了株穗數(shù)性狀的 QTL。此外，位于 2B 染色體標記 wmc154b 附近的 QTL 簇 C5，檢測到與地上部氮累積量、總氮累積量、地上部氮利用率性狀有關，而 Guo 等[12]檢測到與單株總重、根氮累積量、地上部氮累積量、總氮累積量、地上部磷累積量、總磷累積量、地上部鉀累積量、總鉀累積量、地上部氮利用率、總氮利用率、根鉀利用率性狀有關；Sun 等[29]檢測到與單株總干重、單株地上部干重、單株地上部鮮重、根冠比性狀有關；Kong 等[30]檢測到與地上部鉀累積量、總鉀累積量和千粒重性狀有關?？梢?，這些熱點區(qū)域不僅與小麥苗期氮素養(yǎng)分效率的控制有關，還同時與眾多的其他元素營養(yǎng)性狀及產(chǎn)量性狀有關。這些位點可能包含在不同生育階段均起作用的不同或相同基因。這些 QTL 熱點區(qū)域，很可能是基因富集的區(qū)域，值得進一步研究。

4 結(jié)論

磷、鉀的供應能夠顯著影響小麥苗期的氮素吸收利用及其相關 QTL 的表達。影響苗期小麥氮養(yǎng)分效率相關性狀的 QTL 大多數(shù)僅在特定處理下被檢測到，但是大多數(shù) QTL 會被集中定位在特定染色體位點上形成 QTL 簇，構(gòu)成了氮養(yǎng)分效率控制的 QTL 熱點，這些熱點區(qū)域也與前人定位的許多成株期性狀有關，值得深入探討。

表7 本研究與前人研究所定位 QTL 簇的比較Table 7 Comparison of the QTL clusters between previous studies and this study

[1]張福鎖, 崔振嶺, 王激清. 中國土壤和植物養(yǎng)分管理現(xiàn)狀與改進策略[J]. 植物學通報, 2007, 24(6): 687–694. Zhang F S, Cui Z L, Wang J Q. Current status of soil and plant nutrient management in China and improvement strategies [J]. Chinese Bulletin of Botany, 2007, 24(6): 687–694.

[2]Schachtman D P , Shin R. Nutrient sensing and signaling: NPKS [J]. Annual Review of Plant Biology, 2007, 58: 47–69.

[3]Giles J. Nitrogen study fertilizes fears of pollution [J]. Nature, 2005, 433(7028): 791.

[4]Davidson E A. The contribution of manure and fertilizer nitrogen to atmospheric nitrous oxide since 1860 [J]. Nature Geoscience, 2009, 2(9): 659–662.

[5]Yan X, Wu P, Ling H, et al. Plant nutriomics in China: an overview [J]. Annals of Botany, 2006, 98(3): 473–482.

[6]Doerge R W. Mapping and analysis of quantitative trait loci in experimental populations [J]. Nature Reviews Genetics, 2002, 3(1): 43–52.

[7]Quarrie S, Steed A, Calestani C, et al. A high-density genetic map of hexaploid wheat (Triticum aestivum L.) from the cross ChineseSpring × SQ1 and its use to compare QTLs for grain yield across a range of environments [J]. Theoretical and Applied Genetics, 2005, 110(5): 865–880.

[8]An D, Su J, Liu Q, et al. Mapping QTLs for nitrogen uptake in relation to the early growth of wheat (Triticum aestivum L.) [J]. Plant and Soil, 2006, 284(1–2): 73–84.

[9]Laperche A, Devienne-Barret F, Maury O, et al. A simplified conceptual model of carbon/nitrogen functioning for QTL analysis of winter wheat adaptation to nitrogen deficiency [J]. Theoretical and Applied Genetics, 2006, 113(6): 1131–1146.

[10]Laperche A, Brancourt-Hulmel M, Heumez E, et al. Using genotype × nitrogen interaction variables to evaluate the QTL involved in wheat tolerance to nitrogen constraints [J]. Theoretical and Applied Genetics, 2007, 115(3): 399–415.

[11]Fageria N K, The use of nutrients in crop plants [M]. Boca Raton: CRC Press, FL, 2009.

[12]Guo Y, Kong F M, Xu Y F, et al. QTL mapping for seedling traits in wheat grown under varying concentrations of N, P and K nutrients [J]. Theoretical and Applied Genetics, 2012, 124(5): 851–865.

[13]Hoagland D R, Arnon D I. The water-culture method for growing plants without soil [J]. Circular California Agricultural Experiment Station, 1950, 347: 32.

[14]Kjeldahl J. A new method for the determination of nitrogen in organic matter [J]. Analytical and Bioanalytical Chemistry, 1883, 22(1): 366–382.

[15]Siddiqi M Y, Glass A D. Utilization index: a modified approach to the estimation and comparison of nutrient utilization efficiency in plants [J]. Journal of Plant Nutrition, 1981, 4(3): 289–302.

[16]Knapp S, Stroup W, Ross W. Exact confidence intervals for heritability on a progeny mean basis [J]. Crop Science, 1985, 25(1): 192–194.

[17]Wang Y Y, Sun X Y, Zhao Y, et al. Enrichment of a common wheat genetic map and QTL mapping for fatty acid content in grain [J]. Plant Science, 2011, 181(1): 65–75.

[18]Wang S, Basten C , Zeng Z. Windows QTL cartographer 2.5[M]. Raleigh, NC: Department of Statistics, North Carolina State University, 2007.

[19]Churchill G A, Doerge R W. Empirical threshold values for quantitative trait mapping [J]. Genetics, 1994, 138(3): 963–971.

[20]Stoll M, Kwitek-Black A E, Cowley A W, et al. New target regions for human hypertension via comparative genomics [J]. Genome Research, 2000, 10(4): 473–482.

[21]孫海國, 張福鎖. 小麥根系生長對缺磷脅迫的反應[J]. 植物學報, 2000, 42(9): 913–919. Sun H G, Zhang F S. Growth response of wheat roots to phosphorus deficiency [J]. Acta Botanica Sinica, 2000, 42(9): 913–919.

[22]Yuan H, Liu D. Signaling components involved in plant responses to phosphate starvation [J]. Journal of Integrative Plant Biology, 2008, 50(7): 849–859.

[23]譚金芳, 苗玉紅, 張許, 等. 鉀肥不同基追比對冬小麥產(chǎn)量、養(yǎng)分累積量及鉀肥利用效率的影響[J]. 江西農(nóng)業(yè)學報, 2009, 21(12): 46–48. Tan J F, Miao Y H, Zhang X, et al. Effects of ratio of basal to dressing application of potassium fertilizer on yield, nutrient accumulation of winter wheat and use efficiency of potassium fertilizer [J]. Acta Agriculturae Jiangxi, 2009, 21(12): 46–48.

[24]鄒鐵祥, 戴廷波, 姜東, 等. 鉀素水平對小麥氮素積累和運轉(zhuǎn)及籽粒蛋白形成的影響[J]. 中國農(nóng)業(yè)科學, 2006, 39(4): 686–692. Zou T X, Dai T B, Jiang D, et al. Potassium supply affected plant nitrogen accumulation and translocation and grain protein formation in winter wheat [J]. Scientia Agricultura Sinica, 2006, 39(4): 686–692.

[25]Servaites J C. Energy development loading of amino acids and sucrose into the phloem of soybeam [J]. Plant Physiology, 1979, 64(4), 546–550.

[26]Habash D, Bernard S, Schondelmaier J, et al. The genetics of nitrogen use in hexaploid wheat: N utilisation, development and yield [J]. Theoretical and Applied Genetics, 2007, 114(3): 403–429.

[27]Fontaine J X, Ravel C, Pageau K, et al. A quantitative genetic study for elucidating the contribution of glutamine synthetase, glutamate dehydrogenase and other nitrogen-related physiological traits to the agronomic performance of common wheat [J]. Theoretical and Applied Genetics, 2009, 119(4): 645–662.

[28]Laperche A, Le Gouis J, Hanocq E and Brancourt-Hulmel M. Modelling nitrogen stress with probe genotypes to assess genetic parameters and genetic determinism of winter wheat tolerance to nitrogen constraint [J]. Euphytica, 2008, 161(1–2): 259–271.

[29]Sun J J, Guo Y, Zhang G Z, et al. QTL mapping for seedling traits under different nitrogen forms in wheat [J]. Euphytica, 2013, 191(3): 317–331.

[30]Kong F M, Guo Y, Liang X, et al. Potassium (K) effects and QTL mapping for K efficiency traits at seedling and adult stages in wheat [J]. Plant and Soil, 2013, 373(1–2): 877–892.

[31]Xu Y, Wang R, Tong Y, et al. Mapping QTLs for yield and nitrogenrelated traits in wheat: influence of nitrogen and phosphorus fertilization on QTL expression [J]. Theoretical and Applied Genetics, 2014, 127(1): 59–72.

[32]Li S S, Jia J Z, Wei X Y, et al. A intervarietal genetic map and QTL analysis for yield traits in wheat [J]. Molecular Breeding, 2007, 20(2): 167–178.

[33]Su J Y, Zheng Q, Li H W, et al. Detection of QTLs for phosphorus use efficiency in relation to agronomic performance of wheat grown under phosphorus sufficient and limited conditions [J]. Plant Science, 2009, 176(6): 824–836.

QTL mapping of nitrogen related traits under different concentrations of P and K at wheat seedling stage

ZHAO Yan-yan1, YUAN Ya-pei1, LIANG Xue2, GONG Xiao-ping1, WU Chun-hong1, ZHOU Xiu-wen1, GUO Ying1, ZHAO Yan1, LI Si-shen1, KONG Fan-mei1*
( 1 State Key Laboratory of Crop Biology/National Engineering Laboratory for Efficient Utilization of Soil and Fertilizer Resources, Shandong Agricultural University, Tai’an, Shandong 271018, China; 2 Caoxian Bureau of Agriculture, Heze, Shandong 274400, China )

【Objectives】The objective of this study was to detect quantitative trait locus (QTL) for N use efficiency related traits at the seedling stage of wheat under different levels of phosphorus (P) and potassium (K) treatments. The results will deepen our understanding of the relationship between N and P/K, and this knowledge could be applied for map-based cloning and marker assisted selection (MAS) in wheat breeding. 【Methods】Population employed for QTL analysis was a set of 131 RILs derived from a Chuan 35050 × Shannong 483 cross. The 131 RILs and their parental lines were grown under hydroponic culture in greenhouse. Nine treatments were designed. The concentrations of the applied P and K treatments were as follows: moderate P and K (MPMK), high P (HP), low P1 (LP1), low P2 (LP2), low P3 (LP3), high K (HK), low K1 (LK1), low K2 (LK2) and low K3 (LK3). Combined with an enriched genetic map, we conducted research on the N use efficiency traits under different P and K treatments in order to map the QTLs (quantitative trait loci) related to N use efficiency of 10 traits at the wheat seedling stage and to conduct genetic analysis at the whole genome level. 【Results】 A total of 137 QTLs for the 10 seedling traits were detected on 20 chromosomes, with the exception of 3D. Of these QTLs, 122 (89.05%) were detected in only one treatment. Three relatively high-frequency QTLs (QRnue-1A.2, QSnue-1A.1 and QTnue-1A.1) were detected in four treatments at least. Five QTLs (QRnue-1A.1, QTnue-1A.1, QSnc-4A, QRnc-6A.3 and QSnue-6B) were detected in both the low P and low K treatments. Seventeen important QTL clusters (C1-C17) containing at least three or more traits which involved 66 out of 137 QTLs (48.18%) were mapped on chromosomes 1A, 1B, 2B, 2D, 3A, 3B, 4A, 4B, 5D, 6A, 6B, 6D and 7A. Most QTL clusters were detected in both P and K treatments, while five QTL clusters only associated with specific phosphorus or potassium treatments. Many QTL cluster sites in this study were also related to the biomass, yield and other nutrients in the predecessors’ investigation. 【Conclusions】 The supply of P and K significantly affects N efficiency of wheat at the seedling stage and the expression of related QTLs. Most N efficiency related QTLs were only detected in certain treatments, but many QTLs were located at the same hot sites and composed clusters (hot sites). Many of these sites were also co-located QTLs for biomass, yield and other nutrient traits in the previous study. The detection of these QTL/gene hot sites is valuable for further understanding about the genetic control of N efficiency related traits of wheat and their relationship with P and K environment.

wheat; recombinant inbred line (RIL); quantitative trait loci (QTL); N efficiency

S512.1.01

A

1008–505X(2016)06–1523–15

2016–01–18 接受日期：2016–05–18

國家自然科學基金（31201671）；山東省自然科學基金（ZR2010CZ003）項目資助。

趙艷艷（1989—），女，山東濱州人，碩士，主要從事植物營養(yǎng)遺傳研究。E-mail：snzhaoyanyan@163.com

* 通信作者 E-mail: fmkong@sdau.edu.cn