芥菜型油菜苗期氮效率相關(guān)性狀全基因組關(guān)聯(lián)分析

曾健，王茹夢，龔盼，楊瀟，尹曉奇，李江鶴，陳世龍，姚蕾，宋海星，康雷，張振華

芥菜型油菜苗期氮效率相關(guān)性狀全基因組關(guān)聯(lián)分析

曾健1，王茹夢1，龔盼1，楊瀟1，尹曉奇1，李江鶴1，陳世龍1，姚蕾1，宋海星1，康雷2，張振華1

1湖南農(nóng)業(yè)大學資源學院，長沙 410128；2湖南農(nóng)業(yè)大學農(nóng)學院，長沙 410128

【目的】通過對芥菜型油菜苗期氮效率相關(guān)性狀進行全基因組關(guān)聯(lián)分析，挖掘與油菜氮效率相關(guān)性狀顯著關(guān)聯(lián)的SNP位點，預測相關(guān)候選基因，為揭示油菜氮高效分子機制和創(chuàng)制氮高效種質(zhì)提供理論依據(jù)。【方法】以153份芥菜型油菜品系作為關(guān)聯(lián)群體，設置低氮和正常氮2個處理，每個處理設置3個重復，并在2年（2021和2022年）進行2次完全重復的營養(yǎng)液培養(yǎng)試驗。計算植株根冠比和地上部氮含量的相對值（低氮/正常氮），以此作為氮效率相關(guān)性狀進行全基因組關(guān)聯(lián)分析（genome-wide association study，GWAS），挖掘氮效率相關(guān)候選基因。【結(jié)果】芥菜型油菜氮效率相關(guān)性狀表現(xiàn)出豐富的變異，變異幅度為0.21—2.44，變異系數(shù)為22.92%—26.19%。全基因組關(guān)聯(lián)分析檢測到45個顯著關(guān)聯(lián)的SNP位點，其中，第一次根冠比相對值（relative root-shoot ratio 1，RRSR1）和第二次根冠比相對值（relative root-shoot ratio 2，RRSR2）共同關(guān)聯(lián)到的顯著SNP位點有16個，這些位點可解釋的表型變異率為10.69%—15.39%。第一次地上部相對氮含量（relative shoot nitrogen concentration 1，RSNC1）和第二次地上部相對氮含量（relative shoot nitrogen concentration 2，RSNC2）共同關(guān)聯(lián)到的顯著SNP位點共有29個，可解釋的表型變異率為13.22%—23.96%。在顯著SNP位點上下游200 kb范圍內(nèi)共篩選出15個與氮效率相關(guān)候選基因，包括5個與硝酸鹽轉(zhuǎn)運相關(guān)的基因（、、、和）、3個與氮代謝相關(guān)的基因（、和）、4個與生長發(fā)育相關(guān)的基因（BjuCOBL8、、和BjuUP3）和3個與逆境響應相關(guān)的基因（BjuNTP7、BjuJUB1和BjuPYL6）。【結(jié)論】共檢測到45個與氮效率相關(guān)性狀顯著關(guān)聯(lián)的SNP位點，篩選出15個可能與芥菜型油菜氮效率相關(guān)的候選基因。

芥菜型油菜；氮效率相關(guān)性狀；全基因組關(guān)聯(lián)分析；候選基因

0 引言

【研究意義】氮（N）是植物生長所必需的“生命元素”。以氮肥為代表的化學肥料的施用是促進作物增產(chǎn)的關(guān)鍵因素之一。我國由于不合理或過量施用氮肥，氮肥平均利用率僅為30%左右，顯著低于發(fā)達國家[1-2]，也導致一系列的生態(tài)環(huán)境問題[3]，提高作物的氮素營養(yǎng)效率已成為我國實現(xiàn)農(nóng)業(yè)綠色高效和可持續(xù)發(fā)展的迫切任務。芥菜型油菜具有耐貧瘠、抗（耐）病蟲害、抗裂角、抗倒伏等特點，它不僅在干旱地區(qū)具有種植優(yōu)勢，還是進行油菜遺傳改良的寶貴資源[4-5]，探究其氮高效的分子遺傳機制，有利于為探索油菜氮肥利用率品種改良提供新的技術(shù)方法。【前人研究進展】植物吸收利用的兩大無機氮源為硝態(tài)氮和銨態(tài)氮，其中，硝態(tài)氮是旱地作物油菜根系吸收的主要氮素形態(tài)[6-7]。為了更有效地利用土壤中的硝態(tài)氮，植物在長期進化過程中形成了一套高效的硝酸鹽吸收轉(zhuǎn)運系統(tǒng)，分為高親和轉(zhuǎn)運系統(tǒng)（high-affinity transport system，HATS）和低親和轉(zhuǎn)運系統(tǒng)（low-affinity transport system，LATS）。其轉(zhuǎn)運能力受外界硝酸鹽濃度影響[8]，當硝酸鹽濃度低于250 mmol·L-1時，HATS在轉(zhuǎn)運系統(tǒng)中發(fā)揮主要作用，硝酸鹽濃度高于250 mmol·L-1時，LATS在轉(zhuǎn)運系統(tǒng)中發(fā)揮主要作用[9]。植物通過木質(zhì)部將根中的NO3-轉(zhuǎn)運到地上部或直接在根部將其還原為NO2-，NO2-繼續(xù)被還原為NH4+，再同化為植物體所需的有機物質(zhì)[10]。參與該同化過程的主要有硝酸還原酶（nitrate reductase，NR）、亞硝酸還原酶（nitrite reductase，NIR）、谷氨酰胺合成酶（glutamine synthetase，GS）和谷氨酸合酶（glutamate synthase，GOGAT）。NR和NiR主要參與硝酸鹽的兩步還原過程，NR將NO3-還原為NO2-，NiR將NO2-還原為NH4+，而NR的活性遠低于NiR，是硝酸鹽還原過程中的限速酶[11]。GS與GOGAT則將NH4+同化為谷氨酰胺和谷氨酸。對于油菜而言，其生物量被證實與氮素利用效率（nitrogen use efficiency，NUE）存在關(guān)聯(lián)。顧熾明等[12]在2個氮水平條件下對162份油菜種質(zhì)進行了氮效率評價和基因型差異分析，結(jié)果顯示，油菜生物量、主根長、側(cè)根長、根冠比、氮累積量及氮吸收和利用效率在不同氮水平處理間差異極顯著，因此，該研究認為油菜的苗期以生物量作為評價標準來衡量油菜NUE是一種行之有效的途徑。隨著測序技術(shù)的快速發(fā)展和成本的不斷下降，全基因組關(guān)聯(lián)分析成為一種獲取候選基因的有效途徑。Ahmad等[13]通過GWAS、加權(quán)基因共表達網(wǎng)絡和差異表達基因（differentially expressed genes，DEGs）的綜合分析，發(fā)現(xiàn)了16個基因參與了低氮脅迫下油菜的根系發(fā)育，揭示了氮脅迫的遺傳基礎，為油菜氮素利用率的遺傳改良提供了有價值的SNP位點。Li等[14]研究發(fā)現(xiàn)甘藍型油菜氮素吸收和利用效率、光合作用及碳氮代謝是影響其氮利用效率的主要原因，并確定了可以提高甘藍型油菜氮利用效率的候選基因。He等[15]對甘藍型油菜的氮同化基因（nitrogen assimilation genes，NAG）進行了全基因組分析，共鑒定出67個NAG編碼與氮同化有關(guān)的主要酶，還發(fā)現(xiàn)了12個響應氮饑餓的樞紐基因，這些基因可能在氮素利用中發(fā)揮重要作用，為進一步研究NAG在油菜氮素同化中的功能提供了依據(jù)，同時也為NAG在不同養(yǎng)分脅迫下的響應提供了新的思路。【本研究切入點】在長期高氮肥的施用和選擇育種過程中，導致氮高效等重要基因資源的丟失，以致目前主栽品種氮肥利用率普遍偏低[16-18]。而相對于我國廣泛種植的甘藍型油菜，芥菜型油菜具有極豐富的遺傳變異資源，且更耐低氮、耐貧瘠，是挖掘油菜氮高效優(yōu)異基因的寶貴材料。【擬解決的關(guān)鍵問題】本研究通過對153份芥菜型油菜種質(zhì)在不同氮水平下進行苗期氮效率相關(guān)性狀的調(diào)查，以此開展油菜氮效率相關(guān)性狀的全基因組關(guān)聯(lián)分析，鑒定油菜氮效率候選基因，為揭示油菜氮高效分子機制和創(chuàng)制氮高效油菜新品種提供理論基礎和優(yōu)異基因資源。

1 材料與方法

1.1 試驗材料

以153份芥菜型油菜品種（系）組成的自然群體為試驗材料，包括國內(nèi)品種84個、國外品種58個及地理來源未知品種11個。

1.2 試驗設計

試驗于2021和2022年在湖南農(nóng)業(yè)大學進行，參考Hoagland營養(yǎng)液配方配制試驗所用營養(yǎng)液，并根據(jù)芥菜型油菜的營養(yǎng)特性進行適當修改，設置低氮（0.25 mmol·L-1N）和正常氮（4 mmol·L-1N）2個處理，各處理設置3個生物學重復。晝夜溫度設為22℃，光照/黑暗時間為16 h/8 h，光照強度為300— 320 μmol·m-2·s-1，相對濕度為60%—75%。選擇圓潤飽滿種子，將種子用去離子水浸泡24 h，期間換水5—6次。將浸泡后的種子均勻播種在塑料育苗盤表面固定的紗布上，育苗盤中加適量的超純水。育苗7 d后，每個材料選擇長勢一致的3株幼苗移栽到盛有40 L Hoagland營養(yǎng)液的黑色塑料水產(chǎn)盒中。根據(jù)芥菜型油菜幼苗營養(yǎng)需求，出苗后的植株先移到1/4濃度的全量營養(yǎng)液，然后，低氮處理的植株用0.25 mmol·L-1N的營養(yǎng)液培養(yǎng)，正常氮處理的植株依次更換為1/2濃度營養(yǎng)液和全量營養(yǎng)液。培養(yǎng)過程中，每隔3 d更換1次營養(yǎng)液，24 h不間斷通氣，培養(yǎng)18 d收獲。整個試驗重復2次。

1.3 測定指標、方法及計算方法

將樣品收獲后105 ℃殺青30 min，然后65 ℃烘干至恒重，用萬分之一電子天平稱量地上部生物量（shoot dry weight，SDW）與根系生物量（root dry weight，RDW），計算植株根冠比（root-shoot ratio，RSR）=RDW/SDW。采用H2SO4-H2O2消煮法和K9840自動凱式定氮儀測定地上部氮含量（shoot nitrogen concentration，SNC）[19]。通過Excal 2019軟件對表型數(shù)據(jù)進行描述性統(tǒng)計并作圖，以SPSS 19.0軟件對數(shù)據(jù)進行相關(guān)性分析和方差分析。

1.4 SNP、親緣關(guān)系與群體結(jié)構(gòu)分析

根據(jù)湖南農(nóng)業(yè)大學劉忠松課題組480份芥菜型油菜種質(zhì)，重測序后鑒定出約453萬個SNP[20]。利用VCFTOOLS軟件提取出本試驗所需的153份種質(zhì)SNP多態(tài)性信息，共獲得4 529 618個SNP。使用PLINK軟件過濾掉次要等位基因頻率小于5%的位點及缺失率大于80%的位點，最終獲得3 944 772個高質(zhì)量的SNP。親緣關(guān)系分析利用TASSEL V5.0軟件進行，用ADMIXTURE V1.30軟件對=2—10進行聚類運算，交叉驗證錯誤（cross-validation errors，CV）值最低的為最優(yōu)群體結(jié)構(gòu)[21]。

1.5 全基因組關(guān)聯(lián)分析

通過對2次重復試驗（R1和R2）測定的根冠比和地上部氮含量的相對值（低氮/正常氮）進行關(guān)聯(lián)分析。利用TASSEL V5.0軟件中的一般線性模型GLM（群體結(jié)構(gòu)矩陣作為固定效應）和混合線性模型MLM（群體結(jié)構(gòu)和親緣關(guān)系矩陣同時作為固定效應）分別對關(guān)聯(lián)群體進行全基因組關(guān)聯(lián)分析，通過2種模型QQ圖的比較，選擇最優(yōu)模型。Bonferroni閾值為-log10()，顯著關(guān)聯(lián)SNP寬松閾值設為4.0，嚴謹閾值設為4.5。關(guān)聯(lián)分析中所用到的曼哈頓圖和QQ圖均由R軟件中的CM plot包完成。

1.6 候選基因篩選和單倍型分析

基于已公布的芥菜型油菜基因組序列[20]，在全基因組關(guān)聯(lián)分析中檢測到的顯著關(guān)聯(lián)的SNP位點上下游200 kb范圍內(nèi)篩選出與氮效率相關(guān)的基因。進一步利用目前芥菜型油菜基因組的注釋信息及其與擬南芥基因的共線性關(guān)系，根據(jù)擬南芥中已經(jīng)明確功能的氮效率相關(guān)基因，篩選出候選區(qū)段內(nèi)與氮效率相關(guān)同源基因。利用Haploview軟件進行候選基因單倍型分析，并基于候選基因中MAF＞0.05的所有SNP，對單倍型進行分類。

2 結(jié)果

2.1 氮效率相關(guān)性狀的表型變異

通過對153份芥菜型油菜進行氮效率相關(guān)性狀分析發(fā)現(xiàn)，在低氮/正常氮條件下，這些性狀的相對值均表現(xiàn)出廣泛的表型變異。相關(guān)性分析表明，2次重復試驗的根冠比和地上部氮含量的相對值均顯著相關(guān)，相關(guān)系數(shù)超過0.85（表1）。在第一次試驗中，根冠比相對值（RRSR1）變異幅度為0.59—2.44，平均值為1.46，變異系數(shù)為23.97%；地上部氮含量相對值（RSNC1）變異幅度為0.29—0.83，平均值為0.48，變異系數(shù)為22.92%。第二次試驗中，根冠比相對值（RRSR2）變異幅度為0.66—1.83，平均值為1.16，變異系數(shù)為23.28%；地上部氮含量相對值（RSNC2）變異幅度為0.21—0.78，平均值為0.42，變異系數(shù)為26.19%。2次重復試驗的根冠比和地上部氮含量的相對值均呈正態(tài)分布（圖1）。表明這些氮效率相關(guān)性狀受多基因控制，屬于數(shù)量性狀遺傳。

表1 氮效率相關(guān)性狀的表型變異

RRSR：根冠比相對值；RSNC：地上部氮含量相對值；1和2表示同一性狀的第1和第2次試驗；**表示在＜0.01水平下差異顯著。下同

RRSR: relative values of root-shoot ratio; RSNC: relative values of shoot nitrogen concentration; 1 and 2 represent experiment replicate 1 and replicate 2 of the same trait; **Significant at＜0.01. The same as below

圖1 氮效率相關(guān)性狀的頻次分布圖

2.2 親緣關(guān)系分析和群體結(jié)構(gòu)分析

利用3 944 772個高質(zhì)量SNP位點對153份芥菜型油菜種質(zhì)進行親緣關(guān)系和群體結(jié)構(gòu)分析。親緣關(guān)系頻次分布圖（圖2-a）及熱圖顯示（圖2-b），約60%供試材料親緣關(guān)系值等于0，表明整個群體材料之間的親緣關(guān)系較遠，具有足夠的遺傳差異性，符合進行全基因組關(guān)聯(lián)分析的要求。

群體結(jié)構(gòu)分析表明，在=6時，交叉驗證錯誤最小，即=6時為最佳亞群數(shù)目（圖3-a）。而PCA分析結(jié)果表明，前2種成分解釋的遺傳方差分別為12.7%和7.57%，這153份芥菜型油菜群體被劃分為6個亞群（圖3-b）。因此，將=6時生成的群體結(jié)構(gòu)矩陣，用于性狀與標記的關(guān)聯(lián)分析。

2.3 關(guān)聯(lián)分析模型的選擇

采用TASSEL V5.0軟件中的一般線性模型GLM（群體結(jié)構(gòu)矩陣作為固定效應）和混合線性模型MLM（群體結(jié)構(gòu)和親緣關(guān)系矩陣同時作為固定效應）分別對153份芥菜型油菜進行氮效率相關(guān)性狀的全基因組關(guān)聯(lián)分析。通過比較2種模型下QQ圖的分布（圖4），發(fā)現(xiàn)在這些性狀中GLM模型始終優(yōu)于MLM模型。因此，選擇GLM模型作為最終分析模型來尋找關(guān)聯(lián)位點。

a：親緣關(guān)系頻率圖；b：親緣關(guān)系熱圖 a: Frequency plot of relatives; b: Kinship heatmap

a：交叉驗證錯誤值；b：主成分PC1、PC2的散點分布 a: Cross validation error values; b: Plots of the first two principal components (PC1 and PC2)

圖4 氮效率相關(guān)性狀在GLM模型和MLM模型中的QQ圖

2.4 氮效率相關(guān)性狀的全基因組關(guān)聯(lián)分析

選取GLM模型對153份芥菜型油菜材料氮效率相關(guān)性狀進行GWAS分析。由于RRSR1和RRSR2關(guān)聯(lián)到的共同位點較少，采用寬松閾值-log10()＞4（表2、圖5-a和圖5-c），而RSNC1和RSNC2關(guān)聯(lián)到的共同位點較多，所以篩選顯著關(guān)聯(lián)SNP時采用嚴謹閾值-log10()＞4.5（表2、圖5-b和圖5-d）。結(jié)果表明，RRSR1和RRSR2共同關(guān)聯(lián)到的顯著SNP位點有16個，分別位于A01、A02、A04和A05染色體上。這些位點可解釋的表型變異率為10.69%— 15.39%。RSNC1和RSNC2共同關(guān)聯(lián)到的顯著SNP位點共有29個，分別位于A01、A02、A06、A07和B01染色體上。這些位點可解釋的表型變異率為13.22%—23.96%。

2.5 候選基因預測及單倍型分析

共篩選出15個與氮效率相關(guān)候選基因（表3）。其中，RSNC關(guān)聯(lián)到29個顯著SNP位點，篩選出7個氮效率相關(guān)候選基因，包括、、、、、和。RRSR關(guān)聯(lián)到16個顯著SNP位點，篩選出8個氮效率相關(guān)候選基因，分別為、、、、、、和。通過對其進行單倍型分析（圖6），結(jié)果表明，候選基因、、、、、、、、和的單倍型間的RRSR或RSNC均存在顯著差異，單倍型Hap1為它們的優(yōu)勢單倍型。其中，、、和在關(guān)聯(lián)群體中可分為3種單倍型；、和可分為4種單倍型；和可分為5種單倍型；可分為6種單倍型。

圖5 氮效率相關(guān)性狀的全基因組關(guān)聯(lián)分析

表2 關(guān)聯(lián)分析中同一性狀檢測到的共同顯著SNP位點

RRSR和RSNC：兩批數(shù)據(jù)根冠比和地上部氮含量相對值的平均。箱線圖頂部的不同字母表示在P＜0.05水平差異顯著

表3 顯著性SNP位點側(cè)翼序列200 kb內(nèi)的氮效率相關(guān)候選基因

3 討論

3.1 芥菜型油菜氮效率GWAS分析

全基因組關(guān)聯(lián)分析是一種用于探究多種性狀與基因組中遺傳變異間關(guān)系的分析方法，也是解析數(shù)量性狀遺傳機制的有力工具。隨著甘藍型油菜基因組序列的發(fā)布以及SNP芯片的開發(fā)，加快了油菜重要數(shù)量性狀的定位研究[22]。目前，全基因組關(guān)聯(lián)分析已經(jīng)在油菜種子品質(zhì)[23]、干旱脅迫[24]、鹽脅迫[25]、株高[26]、千粒重[27]和開花時間[28]等許多方面實現(xiàn)了應用研究，鑒定了大批控制油菜重要性狀的位點。而關(guān)于芥菜型油菜氮效率相關(guān)性狀的全基因組關(guān)聯(lián)分析及其氮效率相關(guān)基因挖掘的研究卻相對較少。芥菜型油菜是一種重要的農(nóng)業(yè)作物，其氮效率是影響產(chǎn)量和質(zhì)量的關(guān)鍵性狀之一。本研究采用GLM模型對153份芥菜型油菜的氮效率相關(guān)性狀進行全基因組關(guān)聯(lián)分析，共檢測到45個與氮效率相關(guān)性狀顯著關(guān)聯(lián)的SNP位點，分布在A01、A02、A04、A05、A06、A07和B01染色體上。由于GWAS已被用于芥菜型油菜氮效率相關(guān)性狀研究[29-31]，因此，將本研究中的關(guān)聯(lián)結(jié)果與以前的工作進行了比較。遺憾的是，在它們之間沒有發(fā)現(xiàn)相同的SNP。但令人鼓舞的是，AKHATAR等[29]檢測到的9個SNP在本研究檢測到的SNP位點（A02：31293392、A05：22237279和A07：29692636—29693369）的3—9 mb范圍內(nèi)。而前人研究檢測到的一些位點在本研究中未被發(fā)現(xiàn)，可能是受到試驗條件、供試材料及用于關(guān)聯(lián)分析的表型性狀等因素的影響。

3.2 氮效率相關(guān)候選基因預測

本研究挖掘出15個可能與芥菜型油菜氮效率有關(guān)的候選基因，位于A02、A05上的5個基因、、、和編碼硝酸鹽轉(zhuǎn)運蛋白，為硝酸鹽轉(zhuǎn)運蛋白家族基因，介導擬南芥硝酸鹽的吸收和轉(zhuǎn)運[32]。為NRT2家族基因，擬南芥在種子硝酸鹽積累中起著特定的作用[33]，另外還有研究表明，花生能夠提高植物氮代謝過程的氮素利用率[34]。擬南芥在鹽脅迫下促進硝酸鹽向地上部的轉(zhuǎn)運[35]。的擬南芥同源基因，編碼液泡硝酸鹽轉(zhuǎn)運蛋白，介導擬南芥液泡內(nèi)NO3?/H+交換[36]。同時，Li等[14]通過對不同氮素利用效率的甘藍型油菜的基因組和轉(zhuǎn)錄組分析篩選到液泡硝酸鹽轉(zhuǎn)運蛋白。突變增加了其對多胺的抗性，而多胺的運輸或代謝與植株地上部薄壁組織中硝酸鹽的運輸有關(guān)[37]。位于A01、A04上的4個與植物生長發(fā)育相關(guān)的候選基因BjuCOBL8、、和BjuUP3，COBL8參與調(diào)節(jié)細胞壁擴張，與CesA基因家族的一些成員共同表達[38-39]。是一個脫落酸受體基因，是水稻生長發(fā)育和抗旱性的關(guān)鍵調(diào)控因子[40]。是SAUR41亞家族基因中的一員，擬南芥SAUR41亞家族基因是受脫落酸誘導的，可調(diào)節(jié)細胞擴張、離子穩(wěn)態(tài)和耐鹽性[41]。UP3為過氧化物酶體蛋白，植物過氧化物酶體是高度動態(tài)的細胞器，介導一系列對發(fā)育至關(guān)重要的代謝過程[42]。位于A01、A04和A07上的3個基因、和參與植物體內(nèi)氮代謝的過程。編碼天冬酰胺合成酶，參與植物氮代謝的過程[43-44]。陳敬東[45]通過對232份甘藍型油菜種質(zhì)進行全基因組關(guān)聯(lián)分析并結(jié)合轉(zhuǎn)錄組分析，篩選到1個與編碼擬南芥天冬酰胺合成酶家族蛋白的ASNS同源基因，這與本研究中篩選到的候選基因位于同一蛋白家族，可能具有相似的功能。可能在根系初級氮同化中起主要作用[46]。對葉酸的生物合成起著至關(guān)重要的作用而葉酸在氮代謝中具有調(diào)節(jié)作用，因此，可能間接地調(diào)節(jié)氮代謝[47-48]。位于A01、A02和A04上的3個基因、和參與植物對非生物脅迫的反應，JUB1是一個NAC轉(zhuǎn)錄因子，有研究表明通過直接（或部分間接）調(diào)控赤霉素和油菜素內(nèi)酯的生物合成和信號通路，負向調(diào)控整合植物生長和非生物脅迫反應[49]。為BHLH家族轉(zhuǎn)錄因子，可能非生物脅迫反應中發(fā)揮功能[50]。可能參與各種脅迫應答，并在脅迫響應中發(fā)揮功能[51]。HU等[52]通過對418份不同品種甘藍型油菜全基因組關(guān)聯(lián)研究發(fā)現(xiàn)了56個重要農(nóng)藝性狀的628個相關(guān)位點候選基因，包括植株構(gòu)型和產(chǎn)量性狀。將其研究結(jié)果與本研究篩選到的候選基因進行比較，沒有發(fā)現(xiàn)相同的候選基因，這可能是由于試驗條件、供試材料及用于關(guān)聯(lián)分析的表型性狀等因素所導致的。候選基因單倍型分析結(jié)果表明，、、、、、、、、和的單倍型之間存在顯著差異，這意味著這些基因在自然群體中可能存在優(yōu)異的自然等位變異。因此，在后續(xù)的功能驗證中可進一步研究這些基因的差異位點，將有助于揭示芥菜型油菜氮高效的遺傳機制，為提高油菜氮效率和創(chuàng)制氮高效品種提供遺傳理論指導。

4 結(jié)論

獲得45個與氮效率相關(guān)性狀顯著關(guān)聯(lián)的SNP位點。最終篩選出15個可能與芥菜型油菜氮效率直接或間接相關(guān)的候選基因，分別編碼不同的轉(zhuǎn)運蛋白、生長調(diào)節(jié)因子、跨膜蛋白、激素受體及代謝相關(guān)蛋白，可能通過參與跨膜運輸、生長發(fā)育、逆境響應及代謝等途徑來調(diào)控油菜氮效率。

[1] JIN J Y. Changes in the efficiency of fertiliser use in China. Journal of the Science of Food and Agriculture, 2012, 92(5): 1006-1009.

[2] 袁力行, 申建波, 崔振嶺, 張福鎖. 植物營養(yǎng)學科發(fā)展報告. 農(nóng)學學報, 2018, 8(1): 39-43.

YUAN L X, SHEN J B, CUI Z L, ZHANG F S. Advances and prospects of plant nutrition. Journal of Agriculture, 2018, 8(1): 39-43. (in Chinese)

[3] GALLOWAY J N, TOWNSEND A R, ERISMAN J W, BEKUNDA M, CAI Z C, FRENEY J R, MARTINELLI L A, SEITZINGER S P, SUTTON M A. Transformation of the nitrogen cycle: Recent trends, questions, and potential solutions. Science, 2008, 320(5878): 889-892.

[4] 徐愛遐. 芥菜型油菜（）遺傳多樣性及其黃籽與芥酸性狀的分子標記[D]. 楊凌: 西北農(nóng)林科技大學, 2008.

XU A X. Genetic diversity and molecular markers of yellow seed and erucic acid in[D]. Yangling: Northwest Agriculture and Forestry University, 2008. (in Chinese)

[5] 劉淑艷, 劉忠松, 官春云. 芥菜型油菜種質(zhì)資源研究進展. 植物遺傳資源學報, 2007, 8(3): 351-358.

LIU S Y, LIU Z S, GUAN C Y. Advances in germplasm of oilseed. Journal of Plant Genetic Resources, 2007, 8(3): 351-358.(in Chinese)

[6] TANG Y F, SUN X C, HU C X, TAN Q L, ZHAO X H. Genotypic differences in nitrate uptake, translocation and assimilation of two Chinese cabbage cultivars [L. ssp.(L.)]. Plant Physiology and Biochemistry, 2013, 70: 14-20.

[7] 張振華. 作物硝態(tài)氮轉(zhuǎn)運利用與氮素利用效率的關(guān)系. 植物營養(yǎng)與肥料學報, 2017, 23(1): 217-223.

ZHANG Z H. The relationship between nitrate transport and utilization in crop and nitrogen utilization efficiency. Journal of Plant Nutrition and Fertilizer, 2017, 23(1): 217-223.(in Chinese)

[8] SINGH R P, TRIPATHI R D, SINHA S K, MAHESHWARI R, SRIVASTAVA H S. Response of higher plants to lead contaminated environment. Chemosphere, 1997, 34(11): 2467-2493.

[9] CRAWFORD N M, FORDE B G. Molecular and developmental biology of inorganic nitrogen nutrition. The Arabidopsis Book, 2002, 1: e0011.

[10] SCHEIBLE W R, LAUERER M, SCHULZE E D, CABOCHE M, STITT M. Accumulation of nitrate in the shoot acts as a signal to regulate shoot‐root allocation in tobacco. The Plant Journal, 1997, 11(4): 671-691.

[11] CAMPBELL W H. Nitrate reductase structure, function and regulation: bridging the gap between biochemistry and physiology. Annual Review of Plant Physiology and Plant Molecular Biology, 1999, 50: 277-303.

[12] 顧熾明, 韓配配, 胡瓊, 李銀水, 廖祥生, 張志華, 謝立華, 胡小加, 秦璐, 廖星. 甘藍型油菜苗期氮效率評價. 中國油料作物學報, 2018, 40(6): 851-860.

GU C M, HAN P P, HU Q, LI Y S, LIAO X S, ZHANG Z H, XIE L H, HU X J, QIN L, LIAO X. Nitrogen efficiency evaluation in rapeseed (L.) at seedling stage. Chinese Journal of Oil Crop Sciences, 2018, 40(6): 851-860. (in Chinese)

[13] AHMAD N, SU B, IBRAHIM S, KUANG L Q, TIAN Z, WANG X F, WANG H Z, DUN X L. Deciphering the genetic basis of root and biomass traits in rapeseed (L.) through the integration of GWAS and RNA-Seq under nitrogen stress. International Journal of Molecular Sciences, 2022, 23(14): 7958.

[14] LI Q, DING G D, YANG N M, WHITE P J, YE X S, CAI H M, LU J W, SHI L, XU F S. Comparative genome and transcriptome analysis unravels key factors of nitrogen use efficiency inL.. Plant, Cell & Environment, 2020, 43(3): 712-731.

[15] HE X Y, ZHANG H, YE X S, HONG J, DING G D. Nitrogen assimilation related genes in: Systematic characterization and expression analysis identified hub genes in multiple nutrient stress responses. Plants, 2021, 10(10): 2160.

[16] HUANG Y, WANG H, ZHU Y, HUANG X, LI S, WU X, ZHAO Y, BAO Z, QIN L, JIN Y, CUI Y, MA G, XIAO Q, WANG Q, WANG J, YANG X, LIU H, LU X, LARKINS B A, WANG W, WU Y. THP9 enhances seed protein content and nitrogen-use efficiency in maize. Nature, 2022, 612(7939): 292-300.

[17] LIU Y, WANG H, JIANG Z, WANG W, XU R, WANG Q, ZHANG Z, LI A, LIANG Y, OU S, LIU X, CAO S, TONG H, WANG Y, ZHOU F, LIAO H, HU B, CHU C. Genomic basis of geographical adaptation to soil nitrogen in rice. Nature, 2021, 590(7847): 600-605.

[18] TANG W, YE J, YAO X, ZHAO P, XUAN W, TIAN Y, ZHANG Y, XU S, AN H, CHEN G, YU J, WU W, GE Y, LIU X, LI J, ZHANG H, ZHAO Y, YANG B, JIANG X, PENG C, ZHOU C, TERZAGHI W, WANG C, WAN J. Genome-wide associated study identifies NAC42-activated nitrate transporter conferring high nitrogen use efficiency in rice. Nature Communications, 2019, 10(1): 5279.

[19] SIDDIQI M Y, GLASS A D M. Utilization index: A modified approach to the estimation and comparison of nutrient utilization efficiency in plants. Journal of Plant Nutrition, 1981, 4(3): 289-302.

[20] KANG L, QIAN L W, ZHENG M, CHEN L Y, CHEN H, YANG L, YOU L, YANG B, YAN M L, GU Y G, WANG T Y, SCHIESSL S V, AN H, BLISCHAK P, LIU X J, LU H F, ZHANG D W, RAO Y, JIA D H, ZHOU D G, XIAO H G, WANG Y G, XIONG X H, MASON A S, CHRIS PIRES J, SNOWDON R J, HUA W, LIU Z S. Genomic insights into the origin, domestication and diversification of. Nature Genetics, 2021, 53(9): 1392-1402.

[21] ALEXANDER D H, NOVEMBRE J, LANGE K. Fast model-based estimation of ancestry in unrelated individuals. Genome Research, 2009, 19(9): 1655-1664.

[22] CHALHOUB B, DENOEUD F, LIU S Y, PARKIN I A P, TANG H B, WANG X Y, CHIQUET J, BELCRAM H, TONG C B, SAMANS B, CORRéA M, DA SILVA C, JUST J, FALENTIN C, KOH C S, LE CLAINCHE I, BERNARD M, BENTO P, NOEL B, LABADIE K, ALBERTI A, CHARLES M, ARNAUD D, GUO H, DAVIAUD C, ALAMERY S, JABBARI K, ZHAO M X, EDGER P P, CHELAIFA H, TACK D, LASSALLE G, MESTIRI I, SCHNEL N, LE PASLIER M C, FAN G Y, RENAULT V, BAYER P E, GOLICZ A A, MANOLI S, LEE T H, THI V H D, CHALABI S, HU Q, FAN C C, TOLLENAERE R, LU Y H, BATTAIL C, SHEN J X, SIDEBOTTOM C H D, WANG X F, CANAGUIER A, CHAUVEAU A, BéRARD A, DENIOT G, GUAN M, LIU Z S, SUN F M, LIM Y P, LYONS E, TOWN C D, BANCROFT I, WANG X W, MENG J L, MA J X, PIRES J C, KING G J, BRUNEL D, DELOURME R, RENARD M, AURY J M, ADAMS K L, BATLEY J, SNOWDON R J, TOST J, EDWARDS D, ZHOU Y M, HUA W, SHARPE A G, PATERSON A H, GUAN C Y, WINCKER P. Plant genetics. Early allopolyploid evolution in the post-Neolithicoilseed genome. Science, 2014, 345(6199): 950-953.

[23] LI F, CHEN B Y, XU K, WU J F, SONG W L, BANCROFT I, HARPER A L, TRICK M, LIU S Y, GAO G Z, WANG N, YAN G X, QIAO J W, LI J, LI H, XIAO X, ZHANG T Y, WU X M. Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (L.). DNA Research, 2014, 21(4): 355-367.

[24] ZHANG J, MASON A S, WU J, LIU S, ZHANG X C, LUO T, REDDEN R, BATLEY J, HU L Y, YAN G J. Identification of putative candidate genes for water stress tolerance in canola (). Frontiers in Plant Science, 2015, 6: 1058.

[25] 賀亞軍, 吳道明, 游婧璨, 錢偉. 油菜耐鹽相關(guān)性狀的全基因組關(guān)聯(lián)分析及其候選基因預測. 中國農(nóng)業(yè)科學, 2017, 50(7): 1189-1201.

HE Y J, WU D M, YOU J C, QIAN W. Genome-wide association analysis of salt tolerance related traits in Brassica napus and candidate gene prediction. Scientia Agricultura Sinica, 2017, 50(7): 1189-1201. (in Chinese)

[26] LI F, CHEN B Y, XU K, GAO G Z, YAN G X, QIAO J W, LI J, LI H, LI L X, XIAO X, ZHANG T Y, NISHIO T, WU X M. A genome-wide association study of plant height and primary branch number in rapeseed (). Plant Science, 2016, 242: 169-177.

[27] 張春, 趙小珍, 龐承珂, 彭門路, 王曉東, 陳鋒, 張維, 陳松, 彭琦, 易斌, 孫程明, 張潔夫, 傅廷棟. 甘藍型油菜千粒重全基因組關(guān)聯(lián)分析. 作物學報, 2021, 47(4): 650-659.

ZHANG C, ZHAO X Z, PANG C K, PENG M L, WANG X D, CHEN F, ZHANG W, CHEN S, PENG Q, YI B, SUN C M, ZHANG J F, FU T D. Genome-wide association study of 1000-seed weight in rapeseed (L.). Acta Agronomica Sinica, 2021, 47(4): 650-659. (in Chinese)

[28] LI L, LONG Y, ZHANG L B, DALTON-MORGAN J, BATLEY J, YU L J, MENG J L, LI M T. Genome wide analysis of flowering time trait in multiple environments via high-throughput genotyping technique inL.. PLoS ONE, 2015, 10(3): e0119425.

[29] AKHATAR J, GOYAL A, KAUR N, ATRI C, MITTAL M, SINGH M P, KAUR R, RIALCH I, BANGA S S. Genome wide association analyses to understand genetic basis of flowering and plant height under three levels of nitrogen application in(L.) Czern & Coss. Scientific Reports, 2021, 11: 4278.

[30] GUPTA S, AKHATAR J, KAUR P, SHARMA A, SHARMA P, MITTAL M, BHARTI B, BANGA S S. Genetic analyses of nitrogen assimilation enzymes in(L.) Czern & Coss. Molecular Biology Reports, 2019, 46(4): 4235-4244.

[31] GUPTA N, GUPTA M, AKHATAR J, GOYAL A, KAUR R, SHARMA S, GOYAL P, MUKTA A, KAUR N, MITTAL M, SINGH M P, BHARTI B, SARDANA V K, BANGA S S. Association genetics of the parameters related to nitrogen use efficiency inL.. Plant Molecular Biology, 2021, 105(1): 161-175.

[32] CHEN S Y, GU T Y, QI Z A, YAN J, FANG Z J, LU Y T, LI H, GONG J M. Two NPF transporters mediate iron long-distance transport and homeostasis in. Plant Communications, 2021, 2(6): 100244.

[33] CHOPIN F, ORSEL M, DORBE M F, CHARDON F, TRUONG H N, MILLER A J, KRAPP A, DANIEL-VEDELE F. TheATNRT2.7 nitrate transporter controls nitrate content in seeds. The Plant Cell, 2007, 19(5): 1590-1602.

[34] 王娟, 陳皓寧, 石大川, 于天一, 閆彩霞, 孫全喜, 苑翠玲, 趙小波, 牟藝菲, 王奇, 李春娟, 單世華. 花生高親和硝酸鹽轉(zhuǎn)運蛋白基因AhNRT2.7a響應低氮脅迫的功能研究. 中國農(nóng)業(yè)科學, 2022, 55(22): 4356-4372.

WANG J, CHEN H N, SHI D C, YU T Y, YAN C X, SUN Q X, YUAN C L, ZHAO X B, MOU Y F, WANG Q, LI C J, SHAN S H. Functional analysis of AhNRT2.7a in response to low-nitrogen in peanut. Scientia Agricultura Sinica, 2022, 55(22): 4356-4372. (in Chinese)

[35] TAOCHY C, GAILLARD I, IPOTESI E, OOMEN R, LEONHARDT N, ZIMMERMANN S, PELTIER J B, SZPONARSKI W, SIMONNEAU T, SENTENAC H, GIBRAT R, BOYER J C. The Arabidopsis root stele transporter NPF2.3 contributes to nitrate translocation to shoots under salt stress. The Plant Journal, 2015, 83(3): 466-479.

[36] VON DER FECHT-BARTENBACH J, BOGNER M, DYNOWSKI M, LUDEWIG U. CLC-b-mediated NO3-/H+exchange across the tonoplast ofvacuoles. Plant and Cell Physiology, 2010, 51(6): 960-968.

[37] TONG W, IMAI A, TABATA R, SHIGENOBU S, YAMAGUCHI K, YAMADA M, HASEBE M, SAWA S, MOTOSE H, TAKAHASHI T. Polyamine resistance is increased by mutations in a nitrate transporter gene NRT1.3 (AtNPF6.4) in. Frontiers in Plant Science, 2016, 7: 834.

[38] PERSSON S, WEI H R, MILNE J, PAGE G P, SOMERVILLE C R. Identification of genes required for cellulose synthesis by regression analysis of public microarray data sets. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102(24): 8633-8638.

[39] SCHINDELMAN G, MORIKAMI A, JUNG J, BASKIN T I, CARPITA N C, DERBYSHIRE P, MCCANN M C, BENFEY P N. COBRA encodes a putative GPI-anchored protein, which is polarly localized and necessary for oriented cell expansion in. Genes & Development, 2001, 15(9): 1115-1127.

[40] SANTOSH KUMAR V V, YADAV S K, VERMA R K, SHRIVASTAVA S, GHIMIRE O, PUSHKAR S, RAO M V, SENTHIL KUMAR T, CHINNUSAMY V. The abscisic acid receptor OsPYL6 confers drought tolerance to indica rice through dehydration avoidance and tolerance mechanisms. Journal of Experimental Botany, 2021, 72(4): 1411-1431.

[41] QIU T, QI M Y, DING X H, ZHENG Y Y, ZHOU T J, CHEN Y, HAN N, ZHU M Y, BIAN H W, WANG J H. The saur41 subfamily of small auxin up RNA genes is abscisic acid inducible to modulate cell expansion and salt tolerance inseedlings. Annals of Botany, 2020, 125(5): 805-819.

[42] QUAN S, YANG P F, CASSIN-ROSS G, KAUR N, SWITZENBERG R, AUNG K, LI J Y, HU J P. Proteome analysis of peroxisomes from etiolatedseedlings identifies a peroxisomal protease involved in β-oxidation and development. Plant Physiology, 2013, 163(4): 1518-1538.

[43] GAUFICHON L, MARMAGNE A, YONEYAMA T, HASE T, CLéMENT G, TRASSAERT M, XU X L, SHAKIBAEI M, NAJIHI A, SUZUKI A. Impact of the disruption of ASN3-encoding asparagine synthetase ondevelopment. Agronomy, 2016, 6(1): 12.

[44] SIECIECHOWICZ K A, JOY K W, IRELAND R J. The metabolism of asparagine in plants. Phytochemistry, 1988, 27(3): 663-671.

[45] 陳敬東. 基于GWAS與RNA-SEQ技術(shù)挖掘甘藍型油菜氮高效相關(guān)QTLs與候選基因[D]. 武漢: 江漢大學,2022.

CHEN J D. QTLs and candidate genes related to nitrogen efficiency in Brassica napus were identified using GWAS and RNA-seq[D]. Wuhan: Jianghan University, 2022. (in Chinese)

[46] COSCHIGANO K T, MELO-OLIVEIRA R, LIM J, CORUZZI G M.gls mutants and distinct Fd-GOGAT genes. Implications for photorespiration and primary nitrogen assimilation. The Plant Cell, 1998, 10(5): 741-752.

[47] 闞靜, 李莉, 許激揚. 葉酸的生物合成及其代謝工程研究進展. 中國生化藥物雜志, 2009, 30(4): 284-286.

KAN J, LI L, XU J Y. Research progress in folic acid biosynthesis and its metabolic engineering. Chinese Journal of Biochemical Pharmaceutics, 2009, 30(4): 284-286. (in Chinese)

[48] LI W C, LIANG Q J, MISHRA R C, SANCHEZ-MU?OZ R, WANG H, CHEN X, VAN DER STRAETEN D, ZHANG C Y, XIAO Y L. The 5-formyl-tetrahydrofolate proteome links folates with C/N metabolism and reveals feedback regulation of folate biosynthesis. The Plant Cell, 2021, 33(10): 3367-3385.

[49] SHAHNEJAT-BUSHEHRI S, TARKOWSKA D, SAKURABA Y, BALAZADEH S.NAC transcription factor JUB1 regulates GA/BR metabolism and signalling. Nature Plants, 2016, 2: 16013.

[50] NOWAK K, GAJ M D. Stress-related function of bHLH109 in somatic embryo induction in. Journal of Plant Physiology, 2016, 193: 119-126.

[51] 岳路明. 核苷酸轉(zhuǎn)移酶在擬南芥非生物脅迫應答中的功能研究[D]. 深圳: 深圳大學, 2017.

YUE L M. study on the functional of nucleotidyl transferase in abiotic stress response of[D]. Shenzhen: Shenzhen University, 2017. (in Chinese)

[52] HU J H, CHEN B Y, ZHAO J, ZHANG F G, XIE T, XU K, GAO G Z, YAN G X, LI H G, LI L X, JI G X, AN H, LI H, HUANG Q, ZHANG M L, WU J F, SONG W L, ZHANG X J, LUO Y J, PIRES C, BATLEY J, TIAN S L, WU X M. Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding. Nature Genetics, 2022, 54(5): 694-704.

Genome-Wide Association Study of Nitrogen Efficiency Related Traits at Seedling Stage in

ZENG Jian1, WANG RuMeng1, GONG Pan1, YANG Xiao1, YIN XiaoQi1, LI JiangHe1, CHEN ShiLong1, YAO Lei1, SONG HaiXing1, KANG Lei2, ZHANG ZhenHua

1College of Resources, Hunan Agricultural University, Changsha 410128;2College of Agronomy, Hunan Agricultural University, Changsha 410128

【Objective】The genome-wide association analysis was performed to identify SNP loci significantly associated with nitrogen use efficiency (NUE) traits inat seedling stage and to predict the relevant candidate genes, providing a theoretical basis for revealing the molecular mechanism of nitrogen use efficiency in rapeseed and creating nitrogen-efficient germplasm. 【Method】The population of 153resources was used as the analysis population. Two treatments, low N and normal N, were established using three replicates for each treatment, and two replicated nutrient culture trials were conducted over a two-year period (2021 and 2022). The relative values of root-shoot ratio and shoot nitrogen concentration (low/normal N) were calculated and utilized as NUE traits for a genome-wide association study (GWAS) aimed at exploring candidate genes for NUE. 【Result】NUE traits ofresources exhibited abundant variation, ranging from 0.21-2.44 with coefficients of variation of 22.92%-26.19%. The GWAS identified 45 significant SNP loci, among which 16 overlapped between the first relative root-shoot ratio (RRSR1) and the second relative root-shoot ratio (RRSR2), accounting for a phenotype variance range of 10.69%-15.39%. Additionally, 29 significant SNP loci were shared between the first relative shoot nitrogen concentration (RSNC1) and the second relative shoot nitrogen concentration (RSNC2), explaining a phenotype variance range of 13.22%-23.96%. 15 candidate genes for NUE were identified within 200 kb upstream and downstream regions of significant SNP loci, including 5 genes related to nitrate transport (,,,and), 3 genes associated with nitrogen metabolism (,and), 4 genes involved in plant growth and development (,,and) and 3 genes participated in stress response (,and). 【Conclusion】45 SNP loci were detected significantly associated with NUE traits and 15 candidate genes for NUE were identified in this study.

Brassica juncea; nitrogen use efficiency traits; GWAS; candidate genes

10.3864/j.issn.0578-1752.2023.20.003

2023-04-10；

2023-06-09

國家自然科學基金聯(lián)合基金重點項目（U21A20236）、國家自然科學基金面上項目（32072664）、湖南省科技創(chuàng)新領(lǐng)軍人才項目（2022RC3053）、湖南省自然科學基金杰出青年基金（2021JJ0004）

曾健，E-mail：1214223485@qq.com。王茹夢，E-mail：15526461183@163.com。曾健和王茹夢為同等貢獻作者。通信作者康雷，E-mail：kanglei@hunau.edu.cn。通信作者張振華，E-mail：zhzh1468@163.com

（責任編輯 李莉）