摘要：為探究‘蒙農(nóng)三葉草1號’（Trifolium ambiguum Bieb. ‘Mengnong No.1’）與根瘤菌共生過程中相關(guān)基因的分子機制及代謝通路，本研究選取了接種根瘤菌的‘蒙農(nóng)三葉草1號’以及未接種的對照組為研究對象，通過轉(zhuǎn)錄組測序，比較了共生與非共生狀態(tài)下的轉(zhuǎn)錄組差異。結(jié)果表明，與未接種組相比，接種組共篩選出1105個差異表達(dá)基因（Differentially expressed genes，DEGs）。進一步分析發(fā)現(xiàn)，這些DEGs涉及多個代謝通路和生物學(xué)過程，包括苯丙烷生物合成、氨基糖和核苷酸糖代謝、半胱氨酸和甲硫氨酸代謝等。其中，677個基因上調(diào)，428個基因下調(diào)。通過實時熒光定量PCR（qRT-PCR）對部分基因進行驗證，證實了轉(zhuǎn)錄組測序的可靠性。本研究揭示了‘蒙農(nóng)三葉草1號’在根瘤菌共生狀態(tài)下的轉(zhuǎn)錄組特征，為深入理解植物與根瘤菌共生機制提供了重要參考。

關(guān)鍵詞：‘蒙農(nóng)三葉草1號’；接種根瘤菌；差異表達(dá)基因；轉(zhuǎn)錄組學(xué)

中圖分類號：S541.2""" 文獻標(biāo)識碼：A""""" 文章編號：1007-0435（2024）09-2759-10

收稿日期：2023-12-21；修回日期：2024-03-12

基金項目：國家自然科學(xué)基金（31960353）資助

作者簡介：

曹克璠（1998-），男，滿族，內(nèi)蒙古呼倫貝爾人，博士研究生，主要從事牧草種質(zhì)資源研究，E-mail：caokefan1003@163.com；*通信作者Author for correspondence，E-mail：wangmj_0540@163.com

doi：10.11733/j.issn.1007-0435.2024.09.009

引用格式：

曹克璠， 索榮臻， 張慧敏，等.接種根瘤菌對‘蒙農(nóng)三葉草1號’根部轉(zhuǎn)錄水平的影響［J］.草地學(xué)報，2024，32（9）：2759-2768

CAO Ke-fan， SUO Rong-zhen， ZHANG Hui-min，et al.The Impact of Rhizobium Inoculation on the Root Transcriptome of ‘Mengnong Clover No.1’［J］.Acta Agrestia Sinica，2024，32（9）：2759-2768

The Impact of Rhizobium Inoculation on the Root Transcriptome of

‘Mengnong Clover No.1’

CAO Ke-fan1， SUO Rong-zhen1， ZHANG Hui-min1， MA Yi-ming1，

WU Qian1， BAO Li-gao2， WANG Ming-jiu1*

（1. College of Grassland and Resources and Environment/Key Laboratory of Grassland Resources of Ministry of Education，

Inner Mongolia Agricultural University， Hohhot， Inner Mongolia 010018， China； 2. Inner Mongolia Agricultural and Animal

Husbandry Technology Extension Center， Hohhot， Inner Mongolia 010030， China）

Abstract：To explore the molecular mechanisms and metabolic pathways of genes related to the symbiotic relationship between ‘Mengnong Clover No.1’ （T. ambiguum Bieb. ‘Mengnong No.1’） and rhizobia，this study selected ‘Mengnong Clover No.1’ inoculated with rhizobia and an uninoculated control group as research subjects. Through transcriptomic sequencing，we compared the transcriptome differences between symbiotic and non-symbiotic states. The results showed that，compared with the uninoculated group，a total of 1105 differentially expressed genes （DEGs） were identified with significant changes of expression levels in all symbiotic samples. Among them，677 genes were upregulated and 428 genes were downregulated. The reliability of the transcriptomic sequencing was confirmed by real-time quantitative PCR （qRT-PCR） for some of these genes. Further analysis revealed that these DEGs were involved in multiple metabolic pathways and biological processes，including biosynthesis of phenylpropanoids，metabolism of amino sugars and nucleotide sugars，as well as cysteine and methionine metabolism. This study reveals the transcriptomic characteristics of ‘Mengnong Clover No.1’ under symbiotic conditions with rhizobia，and provides important references for a deeper understanding of the symbiotic mechanisms between plants and rhizobia.

Key words：‘Mengnong Clover No. 1’;Rhizobium inoculation;Differentially expressed genes （DEGs）;Transcriptomics

在全球人口持續(xù)增長的背景下，對高效且可持續(xù)的農(nóng)業(yè)生產(chǎn)方式的需求變得尤為迫切［1］。特別是在維持土壤肥力和健康方面，土壤中的氮素作為促進作物生長的關(guān)鍵養(yǎng)分，其高效利用成為現(xiàn)代農(nóng)業(yè)發(fā)展的重點［2］。在這種情況下，自然界中的生物固氮（Biological nitrogen fixation，BNF）過程提供了一個可持續(xù)的氮素獲取途徑［3］。特別是豆科植物與根瘤菌（Rhizobium）之間的共生關(guān)系，即共生固氮（Symbiotic nitrogen fixation，SNF），能夠?qū)⒋髿庵械牡獨廪D(zhuǎn)化為植物可直接利用的形式［4］，從而為農(nóng)業(yè)生產(chǎn)提供了一種不依賴化學(xué)肥料的氮源［5］。

高加索三葉草（Trifolium ambiguum M.Bieb.），也稱為小球三葉草或庫拉三葉草［6-7］。與其它豆科植物相比，高加索三葉草擁有更發(fā)達(dá)的主根和地下分蘗系統(tǒng)［8］，使其顯示出顯著的耐旱和耐寒能力，同時具備較高的營養(yǎng)價值，尤其是其較高的粗蛋白含量［9］。然而，種植高加索三葉草面臨的挑戰(zhàn)包括較低的種子產(chǎn)量和幼苗生長的緩慢，特別是在非原產(chǎn)地區(qū)，高加索三葉草往往難以與固氮菌形成有效的共生關(guān)系［10］。2012年，內(nèi)蒙古農(nóng)業(yè)大學(xué)培育出了‘蒙農(nóng)三葉草1號’，作為中國首個自主培育的高加索三葉草品種，其展現(xiàn)了良好的抗旱、耐寒和耐澇性能［11］。‘蒙農(nóng)三葉草1號’在內(nèi)蒙古地區(qū)表現(xiàn)出了結(jié)瘤固氮的能力，為探究其固氮能力的潛在變化，研究團隊將‘蒙農(nóng)三葉草1號’同樣移栽至青海省西寧市及青藏高原的高海拔地區(qū)。在這些新環(huán)境中，‘蒙農(nóng)三葉草1號’顯示出了自發(fā)結(jié)瘤固氮的現(xiàn)象。為了進一步探究這種現(xiàn)象的分子機制，研究團隊從青海地區(qū)的‘蒙農(nóng)三葉草1號’中收集根瘤菌，進行后續(xù)研究。

轉(zhuǎn)錄組測序（High-throughput RNA sequencing，RNA-Seq）技術(shù)在探索根瘤菌與植物共生過程中差異基因的表達(dá)上發(fā)揮了關(guān)鍵作用。例如Yuan等人對大豆接種了兩個不同屬的根瘤菌（日本慢生根瘤菌和費氏中華根瘤菌）后，大豆出現(xiàn)了較為明顯的轉(zhuǎn)錄水平的反應(yīng)［12］?？紤]到共生過程的復(fù)雜性，結(jié)瘤過程中的動態(tài)轉(zhuǎn)錄組學(xué)分析顯得尤為重要。Niyikiza等分析根瘤菌侵染12，22和36天的根瘤，鑒定出9669個差異表達(dá)基因，并發(fā)現(xiàn)根瘤菌通過調(diào)控植物次級代謝產(chǎn)物合成通路從而影響共生關(guān)系［13］。針對大豆根部，Hayashi等利用轉(zhuǎn)錄組學(xué)分析了野生根瘤菌和nodc-突變菌株侵染2天的大豆根部侵染區(qū)，鑒定出2915個受根瘤菌侵染調(diào)控的基因，其中包含一系列與赤霉素生物合成相關(guān)的基因［14］。Barros de Carvalho等通過轉(zhuǎn)錄組分析了根瘤菌侵染10天的大豆根部材料，鑒定出3210個受根瘤菌侵染顯著調(diào)控的基因，通過富集分析，發(fā)現(xiàn)受調(diào)控的基因主要參與了代謝過程以及刺激反應(yīng)調(diào)控［15］。

本研究集中于探討‘蒙農(nóng)三葉草1號’與根瘤菌共生關(guān)系的轉(zhuǎn)錄組學(xué)特征，并對接種根瘤菌的‘蒙農(nóng)三葉草1號’和未接種根瘤菌的對照組進行了細(xì)致的轉(zhuǎn)錄組學(xué)比較分析，以揭示在共生固氮過程中發(fā)生的分子層面的變化。本研究旨在深入理解根瘤菌如何影響植物的基因表達(dá)，以及這種共生關(guān)系如何促進植物在缺氮環(huán)境下的適應(yīng)和生長。該研究加深了對植物與微生物間相互作用的認(rèn)識，并為發(fā)展更可持續(xù)的農(nóng)業(yè)生產(chǎn)方法提供了關(guān)鍵的科學(xué)支持。

1" 材料與方法

1.1" 植物材料及處理

材料為2012年由內(nèi)蒙古農(nóng)業(yè)大學(xué)培育而成的內(nèi)蒙古牧草新品種‘蒙農(nóng)三葉草1號’（登記號：N010）。本試驗的供試菌株是在前期試驗中分離并鑒定出的一種根瘤菌菌株［9］。試驗于2021年4月在內(nèi)蒙古呼和浩特市內(nèi)蒙古農(nóng)業(yè)大學(xué)草原與資源環(huán)境學(xué)院草地資源教育部重點實驗室人工氣候室中進行。選取籽粒飽滿的‘蒙農(nóng)三葉草1號’種子數(shù)粒，進行表面消毒。將消毒后的種子置于鋪有濾紙的無菌培養(yǎng)皿中，26℃避光培養(yǎng)3～4 d。將發(fā)芽后的種子根部移栽至裝滿蛭石（已滅菌）的10 cm×10 cm的黑色營養(yǎng)缽中，置于28℃，16/8小時光照/黑暗循環(huán)、相對濕度70%的可控環(huán)境室中。每缽移栽3株，澆等量的Hoagland無氮營養(yǎng)液（NS1010-N-50L），培養(yǎng)1 d。在無菌條件下用移液槍將根瘤菌菌懸液加入幼苗根部，每盆1 mL，并以未接種根瘤菌的幼苗為對照組。每盆每天定時補充10 mL營養(yǎng)液。培養(yǎng)3 d后進行根部取樣，并保存在—80℃冰箱備用。

1.2" ‘蒙農(nóng)三葉草1號’RNA提取及文庫的制備

采用Trizol法提取總RNA［16-17］，并用Nanodrop GX方法，RNAse-free H2O檢測RNA樣本的純度和濃度，以確保使用合格的樣本進行轉(zhuǎn)錄組測序。RNA樣品檢測合格后，用帶有Oligo（dT）的磁珠富集真核生物mRNA，以短片段mRNA為模板，用六堿基隨機引物（Random Hexamers）合成第一條cDNA鏈，再利用緩沖液、dNTPs，RNase H和DNA Polymerase I合成第二條cDNA鏈，通過AMPure XP beads純化cDNA，純化的雙鏈cDNA再進行末端修復(fù)、加A尾并連接測序接頭，然后用AMPure XP beads進行片段大小選擇，最后通過PCR富集得到cDNA文庫，cDNA文庫制備完成。使用Agilent 2100生物分析儀驗證后，使用Illumina HiSeqTM 2000測序平臺進行測序。

1.3" 轉(zhuǎn)錄組De novo組裝和功能注釋

為了保證信息分析質(zhì)量，對得到的原始序列數(shù)據(jù)（Raw reads）進行過濾，去除里面含有帶接頭的（Adapter）、N比例大于10%以及低質(zhì)量的reads，得到過濾后序列數(shù)據(jù)（Clean reads）［18］。使用SOAPdenovo進行拼接得到contig，組裝后得到unigene，將unigene序列與蛋白質(zhì)數(shù)據(jù)庫NR，NT，Swiss-Prot，KEGG，COG和GO做blast比對，結(jié)果被用于進一步確定unigene的序列方向。如果不同數(shù)據(jù)庫的結(jié)果存在沖突，則按照NR，Swiss-Prot，Pfam，EggNOG，GO，KEGG優(yōu)先順序確定序列方向。與上述數(shù)據(jù)庫都比對不上的unigene，使用ESTScan軟件預(yù)測其編碼區(qū)域并確定其序列方向。利用Blasta2 go（http：//www.geneontology.org）獲得unigene的GO注釋，將unigene比對到KEGG數(shù)據(jù)庫中進行注釋，最后進行Pathway分析［19］。

1.4" 差異表達(dá)基因分析

采用DESeq2［20］進行樣品組間的差異表達(dá)分析，獲得兩個條件之間的差異表達(dá)基因集。在差異表達(dá)分析過程中采用了公認(rèn)有效的Benjamini-Hochberg方法對原有假設(shè)檢驗得到的顯著性P值進行校正，并最終采用校正后的P值，即假發(fā)現(xiàn)率（False discovery rate，F(xiàn)DR）作為差異表達(dá)基因篩選的關(guān)鍵指標(biāo)。利用GO和KEGG Pathway分析篩選出的基因，得到兩組樣本間差異表達(dá)基因的GO功能分類注釋、顯著性富集分析以及Pathway富集分析。

1.5" qRT-PCR驗證

為了驗證差異表達(dá)基因的結(jié)果，隨機選擇10個差異表達(dá)基因，通過qRT-PCR進行測量。利用Primer Premier 5.0 軟件設(shè)計引物（表1）。根據(jù)TB Green Premix Ex Taq TM Ⅱ試劑盒（Takara公司）說明書進行操作。反應(yīng)體系為TB Green Premix Ex Taq TM Ⅱ 10 μL，上下游引物各1 μL，cDNA 1 μL，RNase-Free H2O 7 μL。使用2-△△Ct法計算基因相對表達(dá)量，所有基因的表達(dá)分析均設(shè)置了3次生物學(xué)重復(fù)［21］。

2" 結(jié)果與分析

2.1" 短片段測序和組裝

本研究運用IlluminaHiSeqTM2000平臺對接種根瘤菌組與對照組的6個樣本進行了轉(zhuǎn)錄組測序，旨在深入探究植物根系在與根瘤菌互作過程中的基因表達(dá)變化。對于接種根瘤菌組，原始序列數(shù)據(jù)數(shù)量分別為43 469 000，42 533 232和42 563 226條。經(jīng)過去除低質(zhì)量的reads、帶接頭的reads以及重復(fù)冗余序列處理后，得到Clean reads數(shù)量分別為43 196 038，42 231 532和42 214 868條。對照組的Raw reads數(shù)量為42 821 958，42 543 248和42 467 922條，經(jīng)過相同處理步驟，得到Clean reads數(shù)量分別為42 459 818，42 038 840和42 020 092條。在對數(shù)據(jù)質(zhì)量的評估中，結(jié)果表明，在接種根瘤菌組和對照組中，Q20（序列質(zhì)量得分超過20的百分比）均達(dá)到94%以上，GC含量均超過42%（表2）。這充分說明本次轉(zhuǎn)錄組測序的質(zhì)量符合國際水平，為后續(xù)的生物信息學(xué)分析提供了有力支持。在進一步比較接種根瘤菌組和對照組的數(shù)據(jù)時，觀察到原始讀數(shù)（Raw reads）在接種根瘤菌組中略高于對照組。然而，經(jīng)過數(shù)據(jù)清洗后，兩組的清洗讀數(shù)（Clean reads）趨向一致。這提示植物在與根瘤菌互作時可能發(fā)生了基因表達(dá)的激增，為進一步探索植物與根瘤菌相互作用提供了引人深思的線索。

本研究生成的數(shù)據(jù)集可在NCBI序列讀取存檔（SRA）數(shù)據(jù)庫中找到，BioProject編號為PRJNA935687（https：//www.ncbi.nlm.nih.gov/bioproject/PRJNA935687）。

2.2" 基因功能注釋

將unigenes在NR，Swiss-Prot，Pfam，EggNOG，GO和KEGG數(shù)據(jù)庫進行blast比對，以預(yù)測unigene序列的功能。最終獲得96 844個unigenes注釋，EggNOG數(shù)據(jù)庫注釋的信息最多，共有36 519個unigenes獲得注釋，占所有unigenes的37.71%，其他數(shù)據(jù)庫注釋結(jié)果如表3所示。

GO（Geneontology）是一個基因功能分類體系，被劃分為生物學(xué)過程、細(xì)胞組分和分子功能三大類別。對‘蒙農(nóng)三葉草1號’轉(zhuǎn)錄組進行GO分類，可以確定基因行使的主要功能。利用GO共得到23 171個功能注釋，被劃歸到46個GO分類中。在分子功能的15個分類中，涉及催化活性和結(jié)合作用的最多；在細(xì)胞組分的12個分類中，大部分歸屬于膜部分、細(xì)胞部分與細(xì)胞器；在生物學(xué)過程的19個分類中，代謝過程與細(xì)胞過程涉及到的基因最多，種類與數(shù)量見圖1。

為鑒定三葉草中涉及的代謝通路，將測序所得的96 844條unigenes與KEGG數(shù)據(jù)庫進行比對分析，結(jié)果共有18 007個unigenes得到注釋，參與20個代謝通路。其中，主要涉及翻譯（1886個）、碳水化合物代謝（1382個）、折疊、排序和降解（1266個）、氨基酸代謝（916個）和運輸與分解代謝（895個），排名前20的代謝途徑見表4。

2.3" 差異表達(dá)基因鑒定

對獲取的轉(zhuǎn)錄組數(shù)據(jù)進行整理，篩選差異表達(dá)基因（DEGs），unigene表達(dá)采用FPKM法計算，篩選閾值為FDR（False discovery rate）lt;0.001，|log2Ratio|≥1，選擇至少有兩倍差異的轉(zhuǎn)錄本，最后鑒定出1105個差異表達(dá)基因，其中上調(diào)基因677個，下調(diào)基因428個（圖2）。

2.4" 差異表達(dá)基因功能注釋

為了探究差異表達(dá)基因的功能信息，通過BLAST軟件將篩選出來的1105條差異表達(dá)基因與GO數(shù)據(jù)庫進行比對，能夠得到注釋信息的差異表達(dá)基因一共有79個?；蚬δ苤饕譃?個方面：細(xì)胞組分、分子功能以及生物過程（圖3）。其中，注釋到細(xì)胞組分方面的差異表達(dá)基因最多，共有62個分布在7個亞類中，主要富集于細(xì)胞組分、膜部分和膜；其次是注釋到生物過程方面共8個亞類的51個基因，主要富集于代謝過程、細(xì)胞過程和對刺激的響應(yīng)；注釋到分子功能方面共4個亞類的基因有45個，主要富集在結(jié)合、催化活性和轉(zhuǎn)運活性等。

2.5" 差異表達(dá)基因KEGG富集分析

將篩選出來的差異表達(dá)基因與KEGG數(shù)據(jù)庫進行比對分析代謝途徑，結(jié)果顯示共有21個差異表達(dá)基因成功富集于22個代謝通路中。其中，6條（P≤0.05）通路被鑒定為顯著富集。苯丙烷類生物合成通路中注釋到了4個基因，氨基糖和核苷酸糖代謝通路中注釋到了3個基因，半胱氨酸和甲硫氨酸代謝通路中注釋到了2個基因。核黃素代謝、異黃酮類物質(zhì)生物合成和磺胺素代謝通路中均富集到了1個差異表達(dá)基因。

2.6" 實時定量PCR驗證

為了驗證RNA-seq的數(shù)據(jù)一致性，利用qRT-PCR的技術(shù)對隨機選擇的15個差異表達(dá)基因進行驗證。結(jié)果如圖5所示，qRT-PCR結(jié)果與RNA-seq的表達(dá)結(jié)果一致，說明轉(zhuǎn)錄組測序數(shù)據(jù)準(zhǔn)確。

3" 討論

近年來，RNA測序（RNA-seq）已被廣泛用于提升對基因調(diào)控和關(guān)鍵生物過程中信號網(wǎng)絡(luò)的理解，成為研究生物基因表達(dá)調(diào)控的重要工具［22］。在本研究采用Illumina HiSeq技術(shù)對接種根瘤菌后的‘蒙農(nóng)三葉草1號’進行全面的轉(zhuǎn)錄組測序，通過與未接種根瘤菌的對照組進行比較，結(jié)果獲得了21.21Gb的Clean reads，其中Q20和Q30均大于90%，表明測序數(shù)據(jù)質(zhì)量良好，為后續(xù)的深入分析提供了堅實的基礎(chǔ)。

通過基因的表達(dá)差異分析，共鑒定到677個上調(diào)基因，428個下調(diào)基因，利用BLAST軟件與GO數(shù)據(jù)庫進行比對后發(fā)現(xiàn)，差異表達(dá)基因主要富集在細(xì)胞組成、膜部分和膜上。這一結(jié)果表明，根瘤菌侵染過程中，植物細(xì)胞結(jié)構(gòu)和組織定位發(fā)生了顯著的變化。這種變化可能是植物為了適應(yīng)共生菌的侵入和根瘤的形成而調(diào)整細(xì)胞壁和膜組成的結(jié)果［23］。細(xì)胞膜作為感知外界信號的重要部位，其上的受體和通道蛋白的表達(dá)變化，可能是對根瘤菌信號分子的直接響應(yīng)［24］。

在生物過程方面，差異表達(dá)基因主要富集在代謝過程、細(xì)胞過程和對刺激的響應(yīng)中。這表明固氮共生過程要求植物調(diào)整其代謝機制，以支持根瘤菌的生存和發(fā)育，同時也涉及到植物細(xì)胞分裂、分化以及程序性細(xì)胞死亡等一系列復(fù)雜過程的調(diào)整［25］。

在分子功能方面，本研究所鑒定的差異表達(dá)基因顯示出在結(jié)合、催化活性和轉(zhuǎn)運活性等方面的顯著富集。這一發(fā)現(xiàn)揭示了根瘤菌感染過程中植物細(xì)胞在分子層面上的復(fù)雜調(diào)整和精細(xì)化反應(yīng)［26］。首先，結(jié)合活性的變化指出了植物細(xì)胞對根瘤菌感染作出反應(yīng)的分子機制。例如，一些蛋白質(zhì)可能增加了與特定信號分子的結(jié)合能力，這可能包括根瘤菌分泌的信號分子或植物自身產(chǎn)生的應(yīng)答分子。這種結(jié)合活性的改變可能是植物感知根瘤菌并啟動相應(yīng)防御機制的關(guān)鍵步驟［27］。其次，催化活性的變化反映了植物代謝路徑的調(diào)整［28］。在固氮共生過程中，植物需要合成新的化合物以支持共生菌，或者調(diào)整現(xiàn)有的代謝途徑以適應(yīng)根瘤菌的需求［29］。因此，某些酶的活性可能會增加，以催化新的或改變的化學(xué)反應(yīng)。最后，轉(zhuǎn)運活性的變化可能與植物細(xì)胞對營養(yǎng)物質(zhì)和信號分子的運輸需求增加有關(guān)。共生狀態(tài)要求植物和根瘤菌之間有有效的物質(zhì)交換，包括營養(yǎng)物質(zhì)、信號分子和其他小分子。因此，相關(guān)的轉(zhuǎn)運蛋白可能在這一過程中起到關(guān)鍵作用，以確保共生菌的正常生長和固氮活動的順利進行［30］。

在本研究中，對差異表達(dá)基因進行KEGG富集分析旨在深入理解根瘤菌感染對‘蒙農(nóng)三葉草1號’代謝途徑的影響。通過這一分析，可以識別出在根瘤菌感染過程中被激活或調(diào)節(jié)的關(guān)鍵生物代謝通路，從而提供了對共生機制更深入的理解。

分析結(jié)果顯示，共有6個顯著富集的代謝通路，反映了植物在根瘤菌感染后的多方面代謝調(diào)整。其中，苯丙烷生物合成通路的顯著富集表明植物可能通過增強這一途徑來增加其抗逆性，尤其是在對生物脅迫的防御反應(yīng)中。特別關(guān)注的是該通路中的過氧化物酶（酶編號1.11.1.7）［31］。過氧化物酶在植物體中扮演著多種功能，尤其是在參與木質(zhì)素合成過程中。木質(zhì)素作為植物細(xì)胞壁的重要組分，不僅增強了植物的結(jié)構(gòu)強度，而且在植物與微生物的互作中發(fā)揮著重要作用。根瘤菌感染初期，植物可能通過調(diào)節(jié)過氧化物酶的活性來管理與根瘤菌的相互作用，影響根瘤的形成和固氮能力［32］。這表明，苯丙烷代謝產(chǎn)物可能作為信號分子參與調(diào)控根瘤的形成，植物細(xì)胞壁的重構(gòu)過程中過氧化物酶可能起到關(guān)鍵作用。在植物與根瘤菌的互作中，1.11.1.7編碼的過氧化物酶的作用可能是雙重的：一方面，它被激活以調(diào)節(jié)根瘤的發(fā)育；另一方面，它可能被適度抑制以防止過度的防御反應(yīng)對共生關(guān)系造成損害［33］。

氨基糖和核苷酸糖的代謝通路也顯示出顯著變化，這些代謝途徑在細(xì)胞壁的構(gòu)建和能量轉(zhuǎn)換中發(fā)揮關(guān)鍵作用［34］。這一結(jié)果暗示，根瘤菌的接種可能激發(fā)了植物細(xì)胞壁重構(gòu)和能量代謝過程的調(diào)整，為根瘤的形成創(chuàng)造條件。此外，半胱氨酸和甲硫氨酸的代謝通路也被顯著富集，這兩種氨基酸在植物的抗氧化過程中具有重要作用。植物可能需要調(diào)整這些氨基酸的代謝來應(yīng)對由根瘤菌感染引起的氧化壓力，或為固氮過程中的蛋白質(zhì)合成提供必要的構(gòu)建塊［35］。

異黃酮生物合成通路中，CYP81E酶的表達(dá)變化也值得關(guān)注。CYP81E屬于細(xì)胞色素P450家族，在植物的次生代謝過程中起著重要作用，包括在植物抗性和植物-微生物相互作用中的功能［36］。這種酶在植物與微生物的互作過程中可能上調(diào)，以合成相關(guān)的防御分子，如異黃酮，它們可能參與抑制病原菌的侵染或調(diào)節(jié)共生微生物的相互作用［37］。CYP81E家族成員催化的化合物可能直接影響根瘤的形態(tài)發(fā)育和功能，調(diào)節(jié)植物的生長發(fā)育，對支持固氮活動至關(guān)重要。這些變化可能是植物響應(yīng)共生菌信號、激活防御機制以及調(diào)節(jié)共生關(guān)系平衡的表現(xiàn)。本研究中發(fā)現(xiàn)的其他顯著富集通路，如核黃素代謝和磺胺素代謝，可能反映了根瘤菌接種對植物次級代謝產(chǎn)物合成、能量代謝和細(xì)胞內(nèi)信號傳導(dǎo)等方面的影響。

總體而言，本研究揭示了植物與根瘤菌互作是一個涉及多層次的復(fù)雜相互調(diào)節(jié)過程。從根瘤的形成到固氮能力的激活，都涉及到細(xì)胞壁的重構(gòu)和分子功能的調(diào)整。植物與根瘤菌的共生關(guān)系在分子層面上呈現(xiàn)出復(fù)雜而微妙的相互作用，這些發(fā)現(xiàn)為未來深入理解植物與根瘤菌共生機制提供了寶貴的線索。

4" 結(jié)論

本研究通過轉(zhuǎn)錄組分析揭示了‘蒙農(nóng)三葉草1號’與根瘤菌共生關(guān)系中的關(guān)鍵分子變化，共篩選出1105個差異表達(dá)基因，其中677個基因上調(diào)，428個基因下調(diào)。這些基因的變化主要集中在細(xì)胞組成和膜部分，反映了植物細(xì)胞在共生狀態(tài)下的結(jié)構(gòu)和功能調(diào)整，特別是苯丙烷生物合成通路、氨基糖和核苷酸糖代謝通路、半胱氨酸和甲硫氨酸代謝通路在共生固氮過程中顯著調(diào)整。這些發(fā)現(xiàn)為理解植物與根瘤菌共生機制提供了新的視角，強調(diào)了植物與微生物相互作用的復(fù)雜性。

參考文獻

［1］" 張琰. ZmABA2參與玉米干旱調(diào)控的研究［D］. 合肥：安徽農(nóng)業(yè)大學(xué)，2022：19-20

［2］" 王瑞，張福耀，詹鵬杰，等. 基于RNA-Seq篩選高粱低氮脅迫相關(guān)候選基因［J］. 作物學(xué)報，2024，50（3）：669-685

［3］" 鮑根生，李媛，馮曉云，等. 高寒區(qū)氮添加和間作種植互作對燕麥和豌豆根系構(gòu)型影響的研究［J］. 草業(yè)學(xué)報，2024，33（3）：73-84

［4］" KOSTYTSCHEW S，SCHULGINA O. ber Alkoholgrung. XIX. Mitteilung. ber einige grung sfhige Mikroben im Macerationssaft［J］. Biological Chemistry，1929，182（50）：1-2

［5］" VANCE G. Legumes：importance and constraints to greater use［J］. Plant Physiology，2003，131（3）：872-877

［6］" BRYANT W G. Caucasian clover （Trifolium ambiguum Bieb.）：A review［J］. Journal of The A Ustralian Institute of Agricultural Science，1974，1（1）：11-19

［7］" TOWNSEND C E. Phenotypic diversity for agronomic characters and frequency of self-compatible plants in trifolium ambiguum［J］. Canadian Journal of Plant Science，1970，50（3）：331-338

［8］" ALLINSON D W，SPEER G S，TAYLOR R W，et al. Nutritional characteristics of kura clover （Trifolium ambiguum Bieb.） compared with other forage legumes［J］. Journal of Agricultural Science，1985，104（1）：227-229

［9］" 劉嘉偉，王明玖，曹克璠，等. 基于代謝組學(xué)的‘蒙農(nóng)三葉草1號’與結(jié)瘤相關(guān)的根系代謝物研究［J］. 草地學(xué)報，2023，31（9）：2703-2711

［10］曹克璠，劉嘉偉，索榮臻，等. 接種根瘤菌對‘蒙農(nóng)三葉草1號’結(jié)瘤固氮及生長的影響［J］. 草地學(xué)報，2023，31（12）：3876-3886

［11］劉嘉偉. 高加索三葉草根瘤菌多樣性及結(jié)瘤固氮機理研究［D］. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2023：10-12

［12］YUAN S，LI R，CHEN S，et al. RNA-Seq analysis of differential gene expression responding to different rhizobium strains in soybean （Glycine max） roots［J］. Frontiers in Plant Science，2016，7：721

［13］NIYIKIZA D，PIYA S，ROUTRAY P，et al. Interactions of gene expression，alternative splicing， and DNA methylation in determining nodule identity［J］. The Plant Journal，2020，103（5）：1744-1766

［14］HAYASHI S，REID D E，LORENC M T，et al. Transient nod factor-dependent gene expression in the nodulation-competent zone of soybean ［Glycine max（L.） Merr.］ roots：Soybean root responses during early nodulation［J］. Plant Biotechnology Journal，2012，10（8）：995-1010

［15］BARROS DE CARVALHO G A，BATISTA J S S，MARCELINO-GUIMARES F C，et al. Transcriptional analysis of genes involved in nodulation in soybean roots inoculated with Bradyrhizobium japonicumstrain CPAC 15［J］. BMC Genomics，2013，14（1）：153

［16］崔慧琳. 結(jié)球甘藍(lán)YL-1遺傳轉(zhuǎn)化體系的建立及花色基因BolC.cpc-1的功能驗證［D］. 合肥：安徽農(nóng)業(yè)大學(xué)，2018：8-10

［17］王麗娜. 楊樹PtPEPCK1調(diào)控碳氮平衡的分子機理研究［D］. 哈爾濱：東北林業(yè)大學(xué)，2018：18-20

［18］周玲芳，尚驍堯，張鐵軍，等. 蒺藜苜蓿（Medicago truncatula）葉片衰老的轉(zhuǎn)錄組分析［J］. 草地學(xué)報，2021，29（10）：2158-2168

［19］樊燕燕. 小針茅與大針茅對干旱脅迫的生理響應(yīng)及分子調(diào)控機制［D］. 呼和浩特：內(nèi)蒙古大學(xué)，2023，18-22

［20］NING L，DAWSON J A，THOMSON J A，et al. EBSeq：an empirical Bayes hierarchical model for inference in RNA-seq experiments［J］. Bioinformatics，2013（8）：2073-2073

［21］吳林，楊蘇渝，楊陽，等. 基于轉(zhuǎn)錄組測序挖掘馬鈴薯塊莖采后發(fā)芽關(guān)鍵基因［J］. 西南大學(xué)學(xué)報（自然科學(xué)版），2023，45（11）：70-80

［22］韓世明. 二代測序技術(shù)在銀杏轉(zhuǎn)錄組和擬南芥突變體分析中的應(yīng)用研究［J］. 武漢：華中師范大學(xué)，2015：12-14

［23］ZHANG S，WANG T，LIMA R M，et al. Widely conserved AHL transcription factors are essential for NCR gene expression and nodule development in Medicago［J］. Nature Plants，2023，9（2）：280-288

［24］劉艷茹. 大豆根瘤菌泌出蛋白NopT調(diào)控共生結(jié)瘤的分子機制分析［D］. 大慶：黑龍江八一農(nóng)墾大學(xué)，2023：18-22

［25］劉麗. 單細(xì)胞轉(zhuǎn)錄組鑒定大豆根瘤菌共生互作的關(guān)鍵基因及應(yīng)用［D］. 泰安：山東農(nóng)業(yè)大學(xué)，2023：20-26

［26］祁平，鄭凱杰，趙曉宇，等. 根瘤菌侵染早期大豆根系的轉(zhuǎn)錄組分析［J］. 大豆科學(xué)，2021，40（3）：289-298

［27］DINKINS R D，HANCOCK J A，BICKHART D M，et al. Expression and variation of the genes involved in rhizobium nodulation in red clover［J］. Plants，2022，11（21）：2888

［28］KE D，F(xiàn)ENG S，HU Y，et al. The soybean NADPH oxidase gene GmRbohK is crucial for nodule formation in symbiosis［J］. Canadian Journal of Plant Science，2023，103（5）：481-493

［29］劉淑霞. 根瘤菌SD101提高紫花苜蓿耐鹽生理和分子機制研究［D］. 大慶：黑龍江八一農(nóng)墾大學(xué)，2022：50-54

［30］刁瑞寧，楊莎，張佳蕾，等. 基于轉(zhuǎn)錄組學(xué)和代謝組學(xué)解析氮調(diào)控花生結(jié)瘤固氮的機理［J］. 山東農(nóng)業(yè)科學(xué)，2023，55（4）：1-10

［31］李宏宇. 磷素水平對大豆結(jié)瘤固氮及相關(guān)代謝途徑的調(diào)控研究［D］. 哈爾濱：東北農(nóng)業(yè)大學(xué)，2022，60-62

［32］陳惠昕. 羅漢松根瘤的發(fā)育解剖學(xué)及其轉(zhuǎn)錄組學(xué)研究［D］. 南寧：廣西大學(xué)，2021，58-61

［33］刁瑞寧，楊莎，張佳蕾，等. 轉(zhuǎn)錄組和代謝組分析鈣調(diào)控花生結(jié)瘤固氮的機理［J］. 花生學(xué)報，2022，51（4）：1-7，17

［34］CHAKRABORTY S，DRISCOLL H E，ABRAHANTE J E，et al. Salt stress enhances early symbi otic gene expression in Medicago truncatula and induces a stress-specific set of rhizobium-responsive genes［J］. Molecular Plant-Microbe Interactions，2021，34（8）：904-921

［35］CERVANTES-PREZ S A，THIBIVILLIERS S，LAFFONT C，et al. Cell-specific pathways recruite d for symbiotic nodulation in the Medicago truncatula legume［J］. Molecular Plant，2022，15（12）：1868-1888

［36］CHEN Y，LI W，JIA K，et al. Metabolomic and transcriptomice analyses of flavonoid biosynthesis in apricot fruits［J］. Frontiers in Plant Science，2023，14（12）：103-109

［37］CERVERA-TORRES C，ARTHIKALA M K，LARA M，et al. Comprehensive analysis of Phaseolus vulgaris SnRK gene family and their expression during rhizobial and mycorrhizal symbiosis［J］. Genes，2022，13（11）：2107

（責(zé)任編輯" 閔芝智）