花生根部結(jié)瘤性狀QTL定位

黃 莉 陳偉剛 李威濤 喻博倫 郭建斌 周小靜 羅懷勇 劉 念 雷 永 廖伯壽 姜慧芳

花生根部結(jié)瘤性狀QTL定位

黃 莉 陳偉剛 李威濤 喻博倫 郭建斌 周小靜 羅懷勇 劉 念 雷 永 廖伯壽 姜慧芳*

中國農(nóng)業(yè)科學(xué)院油料作物研究所 / 農(nóng)業(yè)農(nóng)村部油料作物生物學(xué)與遺傳育種重點(diǎn)實驗室, 湖北武漢 430062

花生是我國重要的豆科油料作物和經(jīng)濟(jì)作物。根瘤是花生共生固氮的重要場所, 研究花生根瘤形成的遺傳基礎(chǔ), 有助于更好地研究花生根瘤固氮能力和固氮特性。然而, 關(guān)于花生根瘤形成的研究較少, 調(diào)控花生根瘤形成的遺傳機(jī)制不清楚。本研究通過對一個高世代RIL群體的根部結(jié)瘤性狀進(jìn)行調(diào)查, 鑒定到7份根部不結(jié)瘤家系, 根部不結(jié)瘤家系的葉綠素含量, 以及株高、單株鮮重、單株干重均顯著低于雙親。利用前期構(gòu)建的SSR標(biāo)記遺傳連鎖圖, 在A08和B07染色體上各鑒定到1個主效QTL和。通過InDel標(biāo)記加密, 將QTL的區(qū)間由4.7 Mb縮小至1.6 Mb, 遺傳變異解釋率由9.1%增加至16.4%;的區(qū)間由9.9 Mb縮小至1.8 Mb, 遺傳變異解釋率由7.1%增加至9.9%。根據(jù)基因功能注釋, 2個QTL區(qū)間分別鑒定到4個和2個結(jié)瘤素基因存在變異位點(diǎn)。本研究將為解析花生根瘤形成發(fā)育以及共生固氮提供理論依據(jù)。

花生; 結(jié)瘤; QTL定位; 連鎖分析

花生(L.)是我國重要的油料與經(jīng)濟(jì)作物之一, 是食用植物油和蛋白質(zhì)的重要來源, 在國民經(jīng)濟(jì)和社會發(fā)展占有重要地位。我國是世界上花生總產(chǎn)量和消費(fèi)量最大的國家, 而且花生的總產(chǎn)、單產(chǎn)、單位面積產(chǎn)油量以及種植業(yè)產(chǎn)值均居國內(nèi)油料作物首位[1]?；ㄉ鳛槎箍谱魑? 根瘤固氮對氮素積累的貢獻(xiàn)率高達(dá)50%, 是花生氮素供應(yīng)的主要方式[2]。根瘤固氮能力直接影響著植株地上部的生長和產(chǎn)量的形成, 根瘤固氮能力越強(qiáng), 植株固氮積累量及產(chǎn)量越高。對花生根瘤固氮相關(guān)指標(biāo)及其與產(chǎn)量的相關(guān)性分析表明根瘤數(shù)和根瘤鮮重與根瘤固氮積累量、供氮比例和莢果產(chǎn)量呈正相關(guān), 高產(chǎn)品種具有良好的固氮能力[3]。因此, 研究花生根瘤形成的遺傳基礎(chǔ), 有助于更好地研究花生根瘤固氮能力和固氮特性, 為花生根瘤高效固氮提供理論依據(jù)。

根瘤的形成和發(fā)育過程受結(jié)瘤因子信號通路(Nod Factor Signaling Pathway)和結(jié)瘤自調(diào)控信號通路(Autoregulation of Nodulation Signaling Pathway)調(diào)控。豆科植物在根際釋放代謝產(chǎn)物類黃酮物質(zhì), 根瘤菌感知類黃酮信號附著于根毛上, 激活結(jié)瘤基因的表達(dá), 釋放結(jié)瘤因子(Nod Factor, NF)到植物根際。根毛表皮細(xì)胞質(zhì)膜上的NFR1 (Nod Factor Receptor 1)和NFR5 (Nod Factor Receptor 5)受體識別結(jié)瘤因子后, 二者形成共受體將信號向下傳遞, 誘導(dǎo)產(chǎn)生鈣震蕩, 進(jìn)而通過依賴于Ca2+、鈣調(diào)素的激酶CCaMK與其磷酸化底物CYCLOPS的相互作用, 激活下游的轉(zhuǎn)錄因子, 如()、()和()等, 最后促進(jìn)下游早期結(jié)瘤相關(guān)基因表達(dá), 如()、(), 從而形成根瘤原基, 即根瘤[4-7]。由于根瘤菌的侵染, 根被誘導(dǎo)生成一種CLE類小肽, 通過植物的傳導(dǎo)結(jié)構(gòu)向地上運(yùn)輸, 并作為富亮氨酸類受體激酶(LRR-RLK)的配體在莖部與LRR-RLK互作, 并使之激活[8]。類受體蛋白激酶可能介導(dǎo)合成與結(jié)瘤有關(guān)的信號分子, 通過傳導(dǎo)結(jié)構(gòu)被運(yùn)送至根部, 調(diào)控根部結(jié)瘤的數(shù)量[9-10]。此外, microRNA也參與根瘤的形成與發(fā)育,能夠靶向切割百脈根的mRNA, 阻礙侵染線的形成并抑制根瘤發(fā)育[11]。光信號能夠通過誘導(dǎo)大豆地上組織合成TGACG-motif結(jié)合因子3/4 (GmSTF3/4)和 FLOWERING LOCUS T (GmFTs)蛋白, 并將這2個蛋白長距離轉(zhuǎn)運(yùn)至根中從而調(diào)控大豆根瘤形成[12]。綜上, 豆科植物根瘤的形成發(fā)育是一個復(fù)雜的過程, 受到多個基因的協(xié)同調(diào)控。

國內(nèi)外在百脈根、大豆等根瘤形成上的研究較多, 但花生根瘤已有的研究主要集中在根瘤菌與花生共生互作過程的生理特征, 關(guān)于根瘤形成發(fā)育的遺傳基礎(chǔ)研究較少?；ㄉ龅男纬砂l(fā)育過程可以分為侵染(1 d)、根瘤原基形成(4 d)、根瘤相似結(jié)構(gòu)形成(8 d)、未成熟根瘤形成(12 d)、成熟根瘤(21 d) 5個階段[13]。模式豆科植物中, 結(jié)瘤因子被受體識別后, 誘導(dǎo)根毛變形彎曲, 形成侵染線, 以“侵染線(Infection Thread)”方式進(jìn)入寄主植物表皮層細(xì)胞。根瘤菌在花生中是以“裂隙侵染(Crack-Entry)”方式通過寄主植物表皮傷口, 從細(xì)胞間隙侵染進(jìn)入, 不同于模式豆科植物[13-14]。研究人員認(rèn)為, 花生根瘤的形成機(jī)制與模式豆科植物根瘤形成機(jī)制存在不同,花生“裂隙侵染”形成根瘤的方式更為古早[15-16]。早期研究發(fā)現(xiàn)通過雜交2個根部結(jié)瘤的花生材料, 其重組自交系后代中會出現(xiàn)不結(jié)瘤家系材料, 遺傳學(xué)分析表明花生根瘤數(shù)受2個獨(dú)立基因調(diào)控, 當(dāng)2個隱性等位基因重組在一起時, 出現(xiàn)根部不結(jié)瘤的現(xiàn)象[17-18]。對侵染根瘤菌5 d后的根進(jìn)行轉(zhuǎn)錄組分析表明, 結(jié)瘤與不結(jié)瘤材料之間表達(dá)差異顯著的基因與一些已知的結(jié)瘤因子信號通路基因同源, 如、、、等[19]。

本研究利用一個高世代RIL群體, 通過對其根瘤多年鑒定, 發(fā)掘根部不結(jié)瘤種質(zhì), 對根部不結(jié)瘤種質(zhì)材料的葉綠素含量以及農(nóng)藝性狀進(jìn)行鑒定, 并通過QTL定位分析, 挖掘花生根部結(jié)瘤性狀的主效QTL, 為花生根瘤形成發(fā)育以及共生固氮提供理論依據(jù)。

1 材料與方法

1.1 試驗材料

以油料所培育的品種中花10號為母本、資源材料ICG12625為父本進(jìn)行雜交, 通過單粒傳法連續(xù)多代自交, 構(gòu)建了一個高世代的重組自交系群體(recombinant inbred lines, RIL), 共包含140個家系。其中, 中花10號為珍珠豆類型花生, ICG12625為赤道類型花生。

1.2 田間試驗及根部結(jié)瘤性狀考察

親本及群體材料分別于2017年和2018年種植于湖北武漢中國農(nóng)業(yè)科學(xué)院油料作物研究所試驗基地, 采用完全隨機(jī)區(qū)組設(shè)計, 2次重復(fù), 每個小區(qū)單行種植, 每行15株, 行間距30 cm, 每行株距10 cm, 常規(guī)田間管理。收獲時, 考慮邊際效應(yīng), 剔除兩邊的植株, 從每行材料中間選取10株, 用清水洗凈根部泥土, 對根部結(jié)瘤情況進(jìn)行觀察鑒定。

1.3 根部不結(jié)瘤材料葉綠素含量及農(nóng)藝性狀鑒定

將RIL群體中鑒定到的穩(wěn)定不結(jié)瘤的家系分別于2019年和2020年種植于武昌湖北武漢中國農(nóng)業(yè)科學(xué)院油料作物研究所試驗基地, 采用完全隨機(jī)區(qū)組設(shè)計, 2次重復(fù), 每個小區(qū)單行種植, 每行15株, 行間距30 cm, 每行株距10 cm, 常規(guī)田間管理。親本和不結(jié)瘤家系每份材料各取10株, 利用SPAD- 502型手持式葉綠素儀在飽果期測量主莖上倒數(shù)第3完全展開葉。收獲時, 調(diào)查選取的每株材料的主莖高和單株植株鮮重, 并調(diào)查單株莢果干重。

1.4 根部結(jié)瘤性狀QTL定位分析

根據(jù)田間調(diào)查的結(jié)果, 將結(jié)瘤性狀分為“結(jié)瘤”和“不結(jié)瘤”, 視為質(zhì)量性狀, 結(jié)合前期構(gòu)建的該群體SSR分子標(biāo)記遺傳連鎖圖[20], 利用Windows QTL Cartographer V. 2.5軟件復(fù)合區(qū)間作圖法(composite interval mapping, CIM)進(jìn)行QTL定位分析。LOD值大于3.0的QTL被認(rèn)為主效QTL。將QTL兩側(cè)SSR標(biāo)記的引物序列在PeanutBase數(shù)據(jù)庫(www. peanutbase.com)中與已發(fā)表的栽培種花生基因組序列進(jìn)行比對, 獲得該標(biāo)記在基因組上的物理位置。利用實驗室已有的親本中花10號和ICG12625的基因組重測序數(shù)據(jù), 使用GATK軟件包檢測QTL區(qū)間內(nèi)兩親本間的SNP和InDel差異位點(diǎn), 將差異片段>3 bp的InDel位點(diǎn)作為目標(biāo)InDel, 將其上下游各500 bp的DNA序列從PeanutBase數(shù)據(jù)庫中下載, 并利用Primer 5.0軟件設(shè)計正向和反向引物。設(shè)計完成的引物序列在PeanutBase數(shù)據(jù)庫中與栽培種花生基因組序列比對, 確認(rèn)該引物是否為目標(biāo)InDel位點(diǎn)的有效擴(kuò)增。利用這些開發(fā)的InDel標(biāo)記對RIL群體進(jìn)行擴(kuò)增, 對主效QTL進(jìn)行加密分析。

1.5 QTL區(qū)間多態(tài)性位點(diǎn)分析

使用ANNOVAR軟件對加密后的QTL區(qū)間內(nèi)的變異位點(diǎn)進(jìn)行注釋, 并預(yù)測變異對基因功能的影響, 根據(jù)變異位點(diǎn)在參考基因組上的位置信息(UTR區(qū)、基因間區(qū)、基因區(qū)或CDS區(qū)等), 判斷變異對基因的影響(同義突變、非同義突變等)。

1.6 QTL區(qū)間基因表達(dá)分析

利用R軟件包“biomaRt” (https://bioconductor. org/packages/release/bioc/html/biomaRt.html)將已發(fā)表的花生多組織轉(zhuǎn)錄組測序獲得的read counts數(shù)據(jù)[21]轉(zhuǎn)換成FPKM值, 分析QTL區(qū)間內(nèi)基因在花生根與根瘤中的表達(dá)水平。

2 結(jié)果與分析

2.1 花生根部結(jié)瘤性狀的鑒定

連續(xù)2年對140份RIL群體及其親本材料的根部結(jié)瘤情況進(jìn)行調(diào)查發(fā)現(xiàn), 親本中花10號和ICG12625根部均有根瘤, RIL群體中出現(xiàn)不結(jié)瘤家系。母本中花10號平均根瘤數(shù)為187個, 父本ICG12625平均根瘤數(shù)為92個。RIL群體中鑒定到7份家系在2年的環(huán)境中穩(wěn)定地表現(xiàn)為植株根部不結(jié)瘤(圖1)。

2.2 花生根部不結(jié)瘤家系性狀鑒定

將這7份不結(jié)瘤家系與結(jié)瘤家系及親本進(jìn)行對比發(fā)現(xiàn), 不結(jié)瘤家系葉片偏黃。對這7份不結(jié)瘤家系和親本材料的葉綠素含量進(jìn)行測定, 不結(jié)瘤家系2019年葉片的葉綠素含量變異范圍在20.41～29.22 SPAD, 2020年葉片的葉綠素含量變異范圍在24.37～ 29.97 SPAD, 而母本材料中花10號2年的葉綠素含量分別為43.91 SAPD和45.80 SPAD, 父本材料ICG12625兩年的葉綠素含量分別為40.92 SAPD和46.18 SPAD。不結(jié)瘤家系的葉片葉綠素含量均顯著低于2個親本材料的葉綠素含量。株高方面, 除L099家系跟母本中花10號株高無差異, 其他6份不結(jié)瘤家系的株高2年均顯著低于雙親材料, 2019年株高變異范圍為30.33～38.52 cm, 2020年株高變異范圍為28.23～39.63 cm。不結(jié)瘤家系的單株植株鮮重變異范圍為119.0～208.8 g, 單株莢果干重變異范圍為16.1～29.6 g, 均顯著低于雙親材料(表1)。

圖1 花生植株根部結(jié)瘤

A: 結(jié)瘤的根; B: 不結(jié)瘤的根。A: root with nodules; B: root without nodules.

2.3 花生根部結(jié)瘤性狀QTL定位

將該RIL群體根部結(jié)瘤性狀分為“結(jié)瘤”與“不結(jié)瘤”, 通過WinQTLCart2.5進(jìn)行QTL定位分析, 鑒定到2個主效QTL和, 分別位于連鎖群A08和B07。QTL的LOD值為3.7, 遺傳變異解釋率為9.1%, 區(qū)間兩側(cè)標(biāo)記分別為AHGS1880和AGGS1495。QTL的LOD值為3.1, 遺傳變異解釋率為7.1%, 區(qū)間兩側(cè)標(biāo)記分別為Ai07B24086和Ai07B26973。將QTL區(qū)間兩側(cè)SSR分子標(biāo)記的引物序列與四倍體栽培種花生Tifrunner基因組序列進(jìn)行比對, QTL定位到染色體A08上9.5～14.2 Mb區(qū)間,定位到染色體B07上118.7～128.6 Mb區(qū)間(表2)。

表1 不結(jié)瘤家系和親本的葉片葉綠素含量和農(nóng)藝性狀鑒定

表2 花生根部結(jié)瘤性狀QTL定位

QTL: quantitative trait locus; PVE: phenotypic variation explained.

2.4 花生根部結(jié)瘤性狀QTL加密分析

利用親本中花10號和ICG12625的重測序數(shù)據(jù),在染色體A08上9.5～14.2 Mb區(qū)間和B07上118.7～128.6 Mb區(qū)間鑒定親本間堿基差異。2個親本在QTL區(qū)間內(nèi)有647個SNP差異和117個InDel差異, 在區(qū)間內(nèi)有1609個SNP差異和564個InDel差異?；贗nDel差異, 在2個QTL區(qū)間分別開發(fā)了10個和19個InDel標(biāo)記, 并在RIL群體中進(jìn)行擴(kuò)增。利用這些InDel標(biāo)記在RIL群體中基因型, 在染色體A08上9.5～14.2 Mb區(qū)間加密構(gòu)建了一個包含10個標(biāo)記、遺傳距離17.6 cM的連鎖群, 在染色體B07上118.7～128.6 Mb區(qū)間加密構(gòu)建了一張包含19個標(biāo)記、遺傳距離18.8 cM的連鎖群。

利用這2個加密的連鎖群, 對不結(jié)瘤性狀進(jìn)行QTL分析, 在這2個連鎖群上仍然鑒定到2個主效QTL (表3和圖2)。區(qū)間為6.9～13.5 cM, LOD值由3.7增加到5.9, 遺傳變異解釋率由9.1%增加到16.4%, 比對基因組序列后, 該QTL區(qū)間物理位置由原來的4.7 Mb縮小至1.6 Mb。區(qū)間為9.9～15.2 cM, LOD值由3.1增加到3.7, 遺傳變異解釋率由7.1%增加到9.9%, 同時該QTL區(qū)間物理位置由原來的9.9 Mb縮小至1.8 Mb。通過標(biāo)記加密, 2個QTL區(qū)間物理位置均縮小至2 Mb以內(nèi)。

表3 標(biāo)記加密后的根部結(jié)瘤性狀主效QTL信息

QTL: quantitative trait locus; PVE: phenotypic variation explained.

圖2 標(biāo)記加密后的根部結(jié)瘤性狀主效QTL

2.5 根部結(jié)瘤性狀QTL加密區(qū)間變異位點(diǎn)分析

QTL加密后,區(qū)間1.6 Mb內(nèi)有237個多態(tài)性位點(diǎn), 其中包括198個SNP差異和39個InDel差異。對這些變異位點(diǎn)進(jìn)行了分析, 接近80%的變異位點(diǎn)位于基因間區(qū), 而在基因內(nèi)的變異位點(diǎn)中, 存在8個外顯子變異位點(diǎn), 其中有5個變異位點(diǎn)屬于非同義突變。區(qū)間1.8 Mb內(nèi)有269個多態(tài)性位點(diǎn), 其中包括199個SNP差異和70個InDel差異。約85%的變異位點(diǎn)位于基因間區(qū), 基因內(nèi)僅有3個變異位點(diǎn)位于外顯子區(qū)域, 分別為1個非同義突變、1個翻譯提前終止和1個非移碼缺失(表4)。

表4 主效QTL候選區(qū)間內(nèi)SNP和InDel注釋結(jié)果

此外, 加密后的區(qū)間和區(qū)間各包含88個和80個基因?；趨⒖蓟蚪M基因功能注釋, 這2個QTL區(qū)間中分別包含6個結(jié)瘤素基因, 其中基因、、和在兩親本之間沒有差異位點(diǎn), 其他8個結(jié)瘤素基因共有42個變異位點(diǎn), 包含26個基因間區(qū)變異位點(diǎn)和16個基因區(qū)變異位點(diǎn), 重要的變異位點(diǎn)信息見表5。

表5 基因區(qū)存在變異位點(diǎn)的結(jié)瘤素基因

2.6 根部結(jié)瘤性狀QTL區(qū)間內(nèi)基因表達(dá)分析

利用已發(fā)表的花生多組織的轉(zhuǎn)錄組數(shù)據(jù), 對QTL和區(qū)間內(nèi)的基因進(jìn)行表達(dá)分析。QTL區(qū)間內(nèi)有19個基因在根和根瘤中均不表達(dá), 有8個基因在根和根瘤中基因表達(dá)FPKM值均小于2。結(jié)瘤素基因和在根中的表達(dá)量顯著高于根瘤的表達(dá), 而和在根瘤中的表達(dá)顯著高于根的表達(dá), 并且基因在根中幾乎不表達(dá)。QTL區(qū)間內(nèi)有22個基因在根和根瘤中均不表達(dá), 有9個基因在根和根瘤中基因表達(dá)FPKM值均小于2。結(jié)瘤素基因和在根中的表達(dá)量顯著高于根瘤的表達(dá), 而和在根瘤中的表達(dá)顯著高于根的表達(dá)。同時在根中幾乎不表達(dá),在根瘤中幾乎不表達(dá)(表6)。

3 討論

花生與根瘤菌互作形成根瘤, 根瘤是花生進(jìn)行生物固氮的重要場所。本研究連續(xù)2年對140份RIL群體及其親本材料的根部結(jié)瘤情況進(jìn)行調(diào)查發(fā)現(xiàn), 親本材料中花10號和ICG12625根部均有根瘤, RIL群體中出現(xiàn)不結(jié)瘤的材料。前人研究分析發(fā)現(xiàn)花生不結(jié)瘤受2對隱性基因共同調(diào)控[17-18]。通過QTL定位分析, 鑒定到2個主效QTL和, 分別位于連鎖群A08和B07, 根據(jù)加性效應(yīng)值推測2個主效QTL分別來源于2個親本。分析不結(jié)瘤家系的基因型發(fā)現(xiàn), 父本ICG12625攜帶的隱性等位基因, 而母本中花10號攜帶的隱性等位基因。只有2個QTL的隱性等位基因聚合在一起時, 才會抑制根瘤的形成, 導(dǎo)致根部不結(jié)瘤。對花生核心種質(zhì)資源的多年田間觀察, 沒有觀察到不結(jié)瘤的材料, 說明雙隱性等位基因的組合在自然界中幾乎不存在。

表6 主效QTL區(qū)間內(nèi)結(jié)瘤素基因在根和根瘤中的表達(dá)

基于分子標(biāo)記和RIL群體的連鎖分析是QTL定位研究的重要方法之一。本研究利用前期已構(gòu)建的SSR分子標(biāo)記的遺傳連鎖圖[20]和RIL群體的根部結(jié)瘤情況, 在A08染色體和B07染色體上分別各定位到一個結(jié)瘤性狀的主效QTL。但是由于SSR標(biāo)記的遺傳連鎖圖定位精度較低, 2個QTL在物理圖譜上區(qū)間還比較大。近年來, 隨著高通量測序技術(shù)的發(fā)展, 基于基因組重測序技術(shù)開發(fā)SNP標(biāo)記和InDel標(biāo)記成為可能。但是由于SNP標(biāo)記需要通過測序技術(shù)檢測結(jié)果, 而InDel標(biāo)記的擴(kuò)增只需通過PAGE膠或者瓊脂糖膠檢測, 因此InDel標(biāo)記是成為加密QTL的最快捷有效的標(biāo)記。本研究利用親本基因組重測序的數(shù)據(jù), 在基于SSR標(biāo)記檢測到的2個主效QTL區(qū)間內(nèi)設(shè)計InDel標(biāo)記進(jìn)行QTL加密。加密后, 2個主效QTL的LOD值和遺傳變異解釋率均有所增加, 同時,區(qū)間由4.7 Mb縮小至1.6 Mb,區(qū)間由9.9 Mb縮小至1.8 Mb。通過InDel標(biāo)記加密, 將QTL區(qū)間均縮小至2 Mb, 這為后續(xù)挖掘候選基因提供了理論基礎(chǔ)。

根瘤的形成和發(fā)育離不開植物基因的調(diào)控和參與。結(jié)瘤素基因是指一些在根瘤形成發(fā)育過程中特異表達(dá)并參與共生固氮作用的植物基因, 對根瘤結(jié)構(gòu)的形成及發(fā)育、固氮產(chǎn)物的代謝都有重要作用[22]。根據(jù)結(jié)瘤素基因表達(dá)時間的早晚分為早期結(jié)瘤素基因()和晚期結(jié)瘤素基因()。在根瘤固氮之前表達(dá)的植物基因被稱為早期結(jié)瘤素基因, 一般在根瘤菌侵染幾個小時后到根瘤形成、固氮作用之前出現(xiàn), 主要參與根瘤的形成和發(fā)育。而晚期結(jié)瘤素基因則在發(fā)育完全的、成熟的根瘤中表達(dá), 參與維持根瘤的結(jié)構(gòu)以及共生固氮[23]。的表達(dá)參與了皮層細(xì)胞分裂, 導(dǎo)致根瘤原基的形成[24-25]; 敲除不影響侵染線的形成, 但會抑制根瘤的形成[26]; 大豆中miR393j-3p可以通過靶向調(diào)控而負(fù)向調(diào)控根瘤的形成[27]。此外, 蒺藜苜蓿中的和[28]、蠶豆中的和[29]、大豆中的[30]均被報道可能涉及豆科植物的根瘤形成和發(fā)育。本研究在2個主效QTL區(qū)間內(nèi)均鑒定到8個結(jié)瘤素基因在基因區(qū)有變異位點(diǎn), 同時部分結(jié)瘤素基因在根和根瘤中的表達(dá)有顯著差異, 但現(xiàn)有的結(jié)果還不能充分確定這些結(jié)瘤素基因是否與調(diào)控花生根瘤形成的關(guān)聯(lián)性, 還需要后續(xù)的試驗驗證。

4 結(jié)論

本研究鑒定到7份根部不結(jié)瘤的花生材料, 根部不結(jié)瘤家系的葉綠素含量以及株高、單株鮮重、單株干重均顯著低于雙親。利用連鎖分析, 在A08和B07染色體上各鑒定到1個主效QTL和。

[1] 廖伯壽. 我國花生生產(chǎn)發(fā)展現(xiàn)狀與潛力分析. 中國油料作物學(xué)報, 2020, 42: 161–166. Liao B S. A review on progress and prospects of peanut industry in China., 2020, 42: 161–166 (in Chinese with English abstract).

[2] 吳正鋒, 陳殿緒, 鄭永美, 王才斌, 孫學(xué)武, 李向東, 王興祥, 石程仁, 馮昊, 于天一. 花生不同氮源供氮特性及氮肥利用率研究. 中國油料作物學(xué)報, 2016, 38: 207–213. Wu Z F, Chen D X, Zheng Y M, Wang C B, Sun X W, Li X D, Wang X X, Shi C R, Feng H, Yu T Y. Supply characteristics of different nitrogen sources and nitrogen use efficiency of peanut., 2016, 38: 207–213 (in Chinese with English abstract).

[3] 鄭永美, 杜連濤, 王春曉, 吳正鋒, 孫學(xué)武, 于天一, 沈浦, 王才斌. 不同花生品種根瘤固氮特點(diǎn)及其與產(chǎn)量的關(guān)系. 應(yīng)用生態(tài)學(xué)報, 2019, 30: 961–968. Zheng Y M, Du L T, Wang C X, Wu Z F, Sun X W, Yu T Y, Shen P, Wang C B. Nitrogen fixation characteristics of root nodules in different peanut varieties and their relationship with yield., 2019, 30: 961–968 (in Chinese with English abstract).

[4] Oldroyd G E, Downie J A. Calcium, kinases and nodulation signalling in legumes., 2004, 5: 566–576.

[5] Smit P, Raedts J, Portyanko V, Debellé F, Gough C, Bisseling T, Geurts R. NSP1 of the GRAS protein family is essential for rhizobia Nod factor-induced transcription., 2005, 308: 1789–1791.

[6] Hirsch S, Kim J, Mu?oz A, Heckmann A B, Downie J A, Oldroyd G E. GRAS proteins form a DNA binding complex to induce gene expression during nodulation signaling in., 2009, 21: 545–557.

[7] Ferguson B J, Indrasumunar A, Hayashi S, Lin M H, Lin Y H, Reid D E, Gresshoff P M. Molecular analysis of legume nodule development and autoregulation., 2010, 52: 61–76.

[8] Oldroyd G E, Downie J A. Coordinating nodule morphogenesis with rhizobial infection in legumes., 2008, 59: 519–546.

[9] Krusell L, Madsen L H, Sato S, Aubert G, Genua A, Szczyglowski K, Duc G, Kaneko T, Tabata S, de Bruijn F, Pajuelo E, Sandal N, Stougaard J. Shoot control of root development and nodulation is mediated by a receptor-like kinase., 2002, 420: 422–426.

[10] Searle I R, Men A E, Laniya T S, Buzas D M, Iturbe-Ormaetxe I, Carroll B J, Gresshoff P M. Long-distance signaling in nodulation directed by a CLAVATA1-like receptor kinase., 2003, 299: 109–112.

[11] De Smet I. Lateral root initiation: one step at a time., 2012, 193: 867–873.

[12] Wang T, Guo J, Peng Y, Lyu X, Liu B, Sun S, Wang X. Light- induced mobile factors from shoots regulate rhizobium-triggered soybean root nodulation., 2021, 374: 65–71.

[13] Kanchan K, Anindya K, Zaigam R A, Emeric D, Dany S, Pierre C, Fabienne C, Maitrayee D. Transcriptomic analysis with the progress of symbiosis in ‘crack-entry’ legumehighlights its contrast with ‘infection thread’ adapted legumes., 2019, 32: 271–285.

[14] Sharma V, Bhattacharyya S, Kumar R, Kumar A, Ibanez F, Wang J, Guo B, Sudini H K, Gopalakrishnan S, DasGupta M, Varshney R K, Pandey M K. Molecular basis of root nodule symbiosis between Bradyrhizobium and ‘crack-entry’ legume groundnut (L.)., 2020, 9: 276.

[15] Kinkema M, Scott P T, Gresshoff P M. Legume nodulation: successful symbiosis through short- and long-distance signaling., 2006, 33: 707–721.

[16] Madsen L H, Tirichine L, Jurkiewicz A, Sullivan J T, Heckmann A B, Bek A S, Ronson C W, James E K, Stougaard J. The molecular network governing nodule organogenesis and infection in the model legume., 2010, 1: 10.

[17] Nigam S N, Nambiar P T C, Dwivedi S L, Gibbons R W, Dart P J. Genetics of nonnodulation in groundnut (L.) studies with single and mixed Rhizobium strains., 1982, 31: 691–693.

[18] Essomba N B, Coffelt T A, Branch W D, Van Scoyoc S W. Inheritance of stem color and non-nodulation in peanut., 1991, 18: 126–131.

[19] Peng Z, Liu F, Wang L, Zhou H, Paudel D, Tan L, Maku J, Gallo M, Wang J. Transcriptome profiles reveal gene regulation of peanut (L.) nodulation., 2017, 7: 40066.

[20] Liu N, Chen H, Huai D, Xia F, Huang L, Chen W, Wu B, Ren X, Luo H, Zhou X, Chen Y, Lei Y, Liao B, Jiang H. Four QTL clusters containing major and stable QTLs for saturated fatty acid contents in a dense genetic map of cultivated peanut (L.)., 2019, 39: 23.

[21] Clevenger J, Chu Y, Scheffler B, Ozias-Akins P. A developmental transcriptome map for allotetraploid., 2016, 7: 1446.

[22] van Kammen A. Suggested nomenclature for plant genes involved in nodulation and symbiosis., 1984, 2: 43–45.

[23] Nap J P, Bisseling T. Developmental biology of a plant- prokaryote symbiosis, the legume root nodule., 1990, 250: 948–954.

[24] Gleason C, Chaudhuri S, Yang T, Mu?oz A, Poovaiah B W, Oldroyd G E. Nodulation independent of rhizobia induced by a calcium-activated kinase lacking autoinhibition., 2006, 441: 1149–1152.

[25] Tirichine L, Imaizumi-Anraku H, Yoshida S, Murakami Y, Madsen LH, Miwa H, Nakagawa T, Sandal N, Albrektsen A S, Kawaguchi M, Downie A, Sato S, Tabata S, Kouchi H, Parniske M, Kawasaki S, Stougaard J. Deregulation of a Ca2+/calmodulin- dependent kinase leads to spontaneous nodule development., 2006. 441: 1153–1156.

[26] Wan X, Hontelez J, Lillo A, Guarnerio C, van de Peut D, Fedorova E, Bisseling T, Franssen H.andare both involved in nodule initiation and bacteroid development., 2007, 58:2033–2041.

[27] Yan Z, Hossain M S, Arikit S, Valdés-López O, Zhai J, Wang J, Libault M, Ji T, Qiu L, Meyers B C, Stacey G. Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation., 2015, 207: 748–759.

[28] Greene E A, Erard M, Dedieu A, Barker D G. MtENOD16 and 20 are members of a family of phytocyanin-related early nodulins., 1998, 36: 775–783.

[29] Vijn I, Yang W C, Pallisg?rd N, Ostergaard Jensen E, van Kammen A, Bisseling T. VsENOD5, VsENOD12 and VsENOD40 expression during Rhizobium-induced nodule formation onroots., 1995, 28: 1111–1119.

[30] de Blank C, Mylona P, Yang W C, Katinakis P, Bisseling T, Franssen H. Characterization of the soybean early nodulin cDNA clone GmENOD55., 1993, 22: 1167–1171.

Identification of major QTLs for nodule formation in peanut

HUANG Li, CHEN Wei-Gang, LI Wei-Tao, YU Bo-Lun, GUO Jian-Bin, ZHOU Xiao-Jing, LUO Huai-Yong, LIU Nian, LEI Yong, LIAO Bo-Shou, and JIANG Hui-Fang*

Oil Crops Research Institute, Chinese Academy of Agricultural Sciences / Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture and Rural Affairs, Wuhan 430062, Hubei, China

Peanut is an important oil and economic crop of legumes in China. Nodules are important for symbiotic nitrogen fixation in peanut. Dissecting the genetic basis of nodule formation could promote the understanding of nitrogen fixation ability and characteristics of nodules in peanut. However, there are few studies on nodule formation and the genetic basis of nodule formation remained unclear in peanut. In this study, the nodules were investigated in a high generation RIL population and seven lines were identified without nodule. The chlorophyll content, plant height, fresh weight per plant, and dry weight per plant of the lines without nodule were significantly lower than those of parents. Using the previous genetic linkage map based on SSR markers, major QTLsandwere identified on chromosome A08 and B07, respectively. The InDel markers in the two major QTLs were developed by the genome resequencing data of parents and performed amplification in the RILs. The interval of QTLsandwas reduced to 1.6 Mb and 1.8 Mb, respectively. Meanwhile, the phenotypic variation explained of these two QTLs was increased to 16.4% and 9.9%, respectively. There were four and two nodulin genes variation loci in the two target intervals, respectively. This study provided the understanding for unveiling the genetic basis of formation and development of nodule and symbiotic nitrogen fixation in peanut.

peanut; nodule formation; QTLs mapping; linkage analysis

2023-02-10;

2023-02-22.

10.3724/SP.J.1006.2023.24184

通信作者(Corresponding author):姜慧芳, E-mail: peanutlab@oilcrops.cn

E-mail: huangli5100@126.com

2022-08-09;

本研究由財政部和農(nóng)業(yè)農(nóng)村部國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(CARS-13)和中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程項目(2022-2060299- 089-031)資助。

This study was supported by the China Agriculture Research System of MOF and MARA (CARS-13) and the Agricultural Science and Technology Innovation Program of Chinese Academy of Agricultural Sciences (2022-2060299-089-031).

URL: https://kns.cnki.net/kcms/detail//11.1809.S.20230221.1719.014.html

This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).