Abstract：To identify the members of the LysM-RLKs family of Stylosanthes guianensis （Stylo） in response to Coletotrichum gloeosporioides and low phosphorus stress，we used transcriptomedata inoculated with Colletot richum gloeosporioides and treated with low phosphorus to identify and name LysM-RLKs family members， and predict their subcellular location and protein sequence by bioinformatics methods. The specificity of tissue expression and the expression pattrns in response to chitin，anthrax and low phosphorus treatment were analyzed through the related tests of Stylosanthes plants.The results showed that 1l genes，members of LysMRLKs family，were characterized by using the corresponding transcriptome data， namely S g L Y K I ， S g L Y K 2 a S g L Y K 2 b ， S g L Y K 3 a ， S g L Y K 3 b ， S g L Y K 4 a ， S g L Y K 4 b ， S g L Y K 5 ， S g L Y K 7 ， S g L Y K 8 and S g N F R 5 ： The prediction of protein localization showed that most of LysM-RLK proteins were localized on the plasma membrane，contained typical extracellular LysM domains，transmembrane domains and intracellular kinase domains.The prediction of protein structure showed that alpha helix and random coil were the major secondary structure of the proteins，and the LysM-RLK proteins were diferent in terms of the three-dimensional structure.The gene expression pattern analysis showed that all the identified LysM-RLKs members had tissue-specific expression pattern. The S g L Y K I ， S g L Y K 4 a ， S g L Y K 5 ， S g L Y K 7 and S g L Y K 8 may play important roles in responding to C ：gloeosporioides infection and chitin-induced immune responses，while S g L Y K 3 a may be involved in the responses to low phosphorus stress.This study could provide candidate genes for the breeding of anthracnose- and low phosphorus-resistant cultivars of stylo.

Key words： Stylosanthes guianensis； LysM-RLK family； Anthracnose； Low phosphorus stress； Expression patterns

植物模式識(shí)別受體（Pattern recognition receptor，PRR）主要有兩類蛋白：類受體樣激酶（Receptor-likekinases，RLKs）和類受體蛋白（Receptor-likeproteins，RLPs）。其中RLKs由胞外結(jié)構(gòu)域、跨膜結(jié)構(gòu)域和激酶結(jié)構(gòu)域組成。常見的胞外結(jié)構(gòu)域有富含亮氨酸重復(fù)序列（Leucine-richrepeats，LRR）溶菌素基序（Lysinmotifs，LysMs）和凝集素型基序（Lectin-typemotifs）[1]。胞外結(jié)構(gòu)域是LysMs的RLKs，含有許多參與識(shí)別病原菌侵染和進(jìn)行免疫防御信號(hào)傳導(dǎo)的防御蛋白，這些蛋白主要通過感知幾丁質(zhì)（Chitin）以誘導(dǎo)植物的免疫反應(yīng)[2]。幾丁質(zhì)是一種保守的病原體相關(guān)分子模式（Pathogen-associatedmolecularpattern，PAMP），是真菌細(xì)胞壁的主要成分[3]。最早被鑒定到能感知幾丁質(zhì)誘導(dǎo)免疫反應(yīng)的LysM-RLKs是擬南芥（Arabidopsisthaliana）中的CERK1 （Chitin elicitor receptor kinase 。

LysM-RLKs基因家族在植物感知病原侵入方面起到關(guān)鍵作用[1]。通過生物信息學(xué)分析和多組學(xué)測序技術(shù)，目前許多植物的LysM-RLKs家族成員已被鑒定出來，例如在擬南芥、百脈根（Lotusjaponicus）大豆（Glycinemax）中分別鑒定到5個(gè)、20個(gè)和27個(gè)LysM-RLKs成員[5-7]。研究表明，LysM-RLKs在幾丁質(zhì)觸發(fā)的免疫反應(yīng)中起著保守和核心作用，同時(shí)在介導(dǎo)叢枝菌根（AM）共生中起到重要作用[1]。例如，在水稻（Oryza satioa）中，OsCERK1與OsMYR1/LYK2形成復(fù)合體，OsLYK2與短鏈殼寡糖（CO4）結(jié)合，介導(dǎo)AM共生8；在百脈根中，LjCERK6是幾丁質(zhì)識(shí)別受體；在蒺藜苜蓿（Medi-cagotruncatula）中，MtLYK9/CERK1在免疫和共生中發(fā)揮雙重作用[9]。

柱花草（Stylosanthesguianensis，Stylo）是重要的熱帶豆科植物，它不僅是重要的熱帶牧草[10]，而且可在恢復(fù)退化土地[1、根系固氮12和保護(hù)生態(tài)環(huán)境13中發(fā)揮重要作用。由膠孢炭疽菌（Colletotri-chumgloeosporioides）引起的柱花草炭疽病（Styloanthracnose）是柱花草的主要病害[14-15]。利用抗病基因培育抗病品種是防控柱花草炭疽病的重要方法[16]。目前，已知的柱花草抗炭疽病基因很少，亟需發(fā)掘用于柱花草抗病育種的相關(guān)基因。另外，柱花草有很好的低磷耐受性[17]，發(fā)掘其抗低磷脅迫基因?qū)ε嘤缕贩N具有重要應(yīng)用價(jià)值。

LysM-RLKs在植物免疫和共生中具有重要作用，但在柱花草中尚未得鑒定和研究。本文利用柱花草接種炭疽菌和低磷處理轉(zhuǎn)錄組數(shù)據(jù)，通過生物信息學(xué)方法鑒定和命名柱花草的LysM-RLKs家族成員，并進(jìn)行幾丁質(zhì)處理、接種炭疽菌以及低磷處理的表達(dá)模式分析。研究成果可為柱花草LysM-RLKs基因家族的功能研究和柱花草抗逆品種的選育提供參考。

1 材料與方法

1. 1 植物材料

試驗(yàn)使用的植物材料是‘熱研5號(hào)'柱花草（ losanthesguianensis‘ReyanNo.5'），其種子由中國熱帶農(nóng)業(yè)科學(xué)院熱帶作物品種資源研究所草業(yè)研究中心提供。柱花草種子在熱水（ 中浸泡3分鐘后，在黑暗條件下發(fā)芽 。一批幼苗移栽到裝有蛭石和營養(yǎng)土 的育苗盒中，在溫室大棚中生長30d后，分別進(jìn)行組織采集、接種炭疽菌[18]；另一批幼苗在改良的Magnavaca營養(yǎng)液中水培14d[19]后，進(jìn)行磷處理。

1.2 生物信息學(xué)分析

1.2.1 柱花草LysM-RLKs家族成員的鑒定利用保存于NCBI數(shù)據(jù)庫（http：//www.ncbi.nlm.nih.gov）的柱花草接種炭疽菌[20]和磷處理[21]轉(zhuǎn)錄組數(shù)據(jù)，其中接種炭疽菌轉(zhuǎn)錄組數(shù)據(jù)編號(hào)為SUB8868028，磷處理轉(zhuǎn)錄組數(shù)據(jù)編號(hào)為PRJNA660932，使用TBtools的“BatchTranslateCDStoProtein\"工具將編碼序列（Codingsequence，CDS）翻譯成蛋白質(zhì)序列。從Pfam數(shù)據(jù)庫（http：//pfam.xfam.org/）下載LysM結(jié)構(gòu)域（PF01476）和Pkinase結(jié)構(gòu)域（PFOOO69）的隱馬爾可夫模型（HMM），使用TBtools的“SimpleHMMSearch\"工具初步篩選出柱花草LysM-RLKs家族成員，然后通過SMART（http：//smart.embl-heidelberg.de/）網(wǎng)站進(jìn)行結(jié)構(gòu)域驗(yàn)證，去掉不含有LysM和Pkinase結(jié)構(gòu)域的序列，最后將上述兩種方法獲得的家族成員用BioEdit軟件進(jìn)行BLAST比對(duì)，去除重復(fù)序列，從而獲得最終的柱花草LysM-RLKs家族成員。

1.2.2系統(tǒng)發(fā)育分析在擬南芥數(shù)據(jù)庫（https：//www.arabidopsis.org/browse/gene_family）下載擬南芥的LysM-RLKs基因家族成員，在Phytozome數(shù)據(jù)庫（https：//phytozome-next. jgi.doe.gov/）下載大豆和百脈根的基因家族成員。通過MAFFT軟件將擬南芥（At）、大豆（Gm）、百脈根（Lj）、柱花草（Sg）的LysM-RLKs基因家族成員進(jìn)行多序列比對(duì)，默認(rèn)參數(shù)。使用MEGA11軟件選擇鄰接法（Neighbor-Joining，NJ）構(gòu)建進(jìn)化樹，設(shè)置Boot-strap值為1000來檢測進(jìn)化樹的可靠性，選擇pdistance和pairwisedeletion，其他為默認(rèn)值。使用ITOL網(wǎng)站（https：//itol.embl.de/）進(jìn)行進(jìn)化樹美化。

1.2.3理化性質(zhì)及亞細(xì)胞定位預(yù)測分析使用WoLFPSORT網(wǎng)站（https：//wolfpsort.hgc. jp/）預(yù)測亞細(xì)胞定位。使用ExPASy網(wǎng)站（https：//web.expasy.org/protparam/）預(yù)測蛋白質(zhì)的長度、分子量、理論等電點(diǎn)、不穩(wěn)定指數(shù)和親水性。

1.2.4保守基序（Motif）和結(jié)構(gòu)域（Domain）分析利用保守基序在線預(yù)測工具M(jìn)EME（https：//meme-suite.org/meme）對(duì)SgLysM-RLKs蛋白序列進(jìn)行預(yù)測，保守基序數(shù)設(shè)置為10，其他參數(shù)設(shè)置為默認(rèn)。使用TBtools的“BatchSMART\"工具預(yù)測SgLysM-RLKs蛋白的保守結(jié)構(gòu)域，以及進(jìn)行可視化。1.2.5蛋白質(zhì)二、三級(jí)結(jié)構(gòu)預(yù)測在SPOMA網(wǎng)站（https：//npsa-prabi. ibcp. fr/cgi-bin/npsa_automat.pl？ page °eq /NPSA/npsa_sopma_f.html）預(yù)測蛋白質(zhì)二級(jí)結(jié)構(gòu)占比情況。在ALPHA-FOLD2網(wǎng)站（https：//colab.research. google.com/github/sokrypton/ColabFold/blob/main/AlphaFold2.ipynb）預(yù)測蛋白3D結(jié)構(gòu)，使用PyMOL軟件進(jìn)行可視化。

1.3SgLysM-RLKs家族成員響應(yīng)炭疽菌及低磷脅迫的轉(zhuǎn)錄組學(xué)分析

使用前面提到的兩組柱花草轉(zhuǎn)錄組數(shù)據(jù)，以 （FPKM + 0 . 1 ） 值表述SgLysM-RLKs基因在炭疽菌侵染以及低磷處理?xiàng)l件下的表達(dá)模式，使用TBtools的\"HeatMap\"工具制作熱圖。

1.4SgLysM-RLKs家族成員的表達(dá)模式分析

1.4.1組織特異性表達(dá)模式分析采集4周齡柱花草植株的根、葉、莖組織，檢測SgLysM-RLK家族成員組織特異性的表達(dá)水平，設(shè)置3個(gè)生物學(xué)重復(fù)。1.4.2幾丁質(zhì)處理根據(jù)高靜等22]的方法進(jìn)行柱花草無菌苗的培育。柱花草生長20d左右時(shí)，使用 幾丁質(zhì)噴灑在柱花草葉片上，分別在噴灑后 0 ， 0 . 2 5 ， 0 . 5 ， 1 ， 2 ， 4 和6h收集柱花草的葉片，每個(gè)時(shí)間點(diǎn)3個(gè)生物學(xué)重復(fù)。

1.4.3炭疽菌處理使用課題組保存的膠孢炭疽菌菌株WC-O2和DZ-19[23]，其中WC-02為弱致病菌，DZ-19為強(qiáng)致病菌。炭疽菌在馬鈴薯葡萄糖瓊脂（PDA）平板上生長7d后，轉(zhuǎn)移到CM培養(yǎng)液中培養(yǎng)3d，用雙層紗布過濾菌絲，收集孢子，并用無菌水稀釋至每毫升 個(gè)分生孢子；用含有 0 . 0 2 % sil-weetL-77的孢子懸浮液噴霧接種柱花草，在0，0.5，1，4，8，12和 收集柱花草的葉片，每個(gè)時(shí)間點(diǎn)3個(gè)生物學(xué)重復(fù)。

1.4.4低磷處理正常水培柱花草14d，挑選長勢良好、均一的植株分別在缺磷（LP， ）以及正常磷供給（HP，300 的改良Magnavaca營養(yǎng)液中繼續(xù)培養(yǎng)，每3d換1次營養(yǎng)液，10d后分別收取地上部和根部樣品，設(shè)置3個(gè)生物學(xué)重復(fù)。

1.4.5LysM-RLKs家族成員的表達(dá)模式分析所有材料均通過液氮速凍收集，并儲(chǔ)存在一 下直

至提取RNA。

使用TRNzolUniversalTotalRNA提取試劑盒（TIANGEN，貨號(hào)DP424）提取植物中的總RNA，使用HiScriptIII RT SuperMix forqPCR（Vazyme，貨號(hào)R323）試劑盒反轉(zhuǎn)為cDNA，使用ChamQUniversalSYBRqPCRMasterMix（ToloBio，貨號(hào)22204）試劑進(jìn)行實(shí)時(shí)定量聚合酶鏈?zhǔn)椒磻?yīng)（Real-TimeQuantitative Polymerase Chain Reaction，RT-qPCR），以 作為內(nèi)部參考基因[24]，表達(dá)水平通過 方法確定。

引物使用primer5軟件設(shè)計(jì)，使用NCBI網(wǎng)站的“Primer-BLAST\"工具（https：//www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi）檢測引物特異性，所有引物如表1所示。

使用Excel和GraphPadPrism軟件進(jìn)行數(shù)據(jù)整理分析。

2 結(jié)果與分析

2.1SgLysM-RLKs家族成員的系統(tǒng)發(fā)育分析與命名

通過分析柱花草轉(zhuǎn)錄組數(shù)據(jù)[20-21]，共鑒定到11個(gè)SgLysM-RLKs家族成員，其中在接種炭疽菌的轉(zhuǎn)錄組中鑒定到5個(gè)成員，在低磷處理的轉(zhuǎn)錄組中鑒定到了11個(gè)成員（其中5個(gè)與接種炭疽菌轉(zhuǎn)錄組重復(fù)）。這表明，SgLysM-RLKs在柱花草中產(chǎn)生了功能分化，不僅有抵抗病原菌的作用，而且還進(jìn)化出了應(yīng)對(duì)低磷脅迫的作用。

通過對(duì)比柱花草近緣物種大豆GmLysM-RLK家族基因序列，在鑒定到的11個(gè)成員中，發(fā)現(xiàn)7條

SgLysM-RLKs具有完整的開放閱讀框（Openread-ingframe，ORF），而Unigene29792（SgLYK2a），Unigene43983（SgLYK2b），Unigenel5056（SgLYK3b）和Unigene20089（SgLYK7）這4條無完整ORF序列。因此使用快速擴(kuò)增cDNA末端聚合酶鏈?zhǔn)椒磻?yīng)（Rapid Amplification of cDNA Ends PolymeraseChainReaction，RACE-PCR）試劑盒（Accurate，貨號(hào)AG11618）進(jìn)行 試驗(yàn)，獲得其全長ORF序列。將獲取到Unigene29792（SgLYK2a），Unigene43983（SgLYK2b），Unigene15056（SgLYK3b）和Unigene20089（SgLYK7）的全長ORF序列上傳到NCBI數(shù)據(jù)庫中，GenBank登錄號(hào)為PV098358至PV098361。

由于SgLysM-RLKs尚未命名，因此利用柱花草、擬南芥、大豆和百脈根LysM-RLKs之間的進(jìn)化關(guān)系進(jìn)行命名。

采用鄰接法構(gòu)建上述4種植物共63個(gè)LysM-RLKs進(jìn)化樹，其系統(tǒng)發(fā)育分析表明，柱花草、擬南芥、大豆和百脈根的LysM-RLKs家族成員被分為7個(gè)亞組（圖1）。其中：亞組V只有大豆；亞組Ⅲ、Ⅳ有柱花草、大豆和百脈根，沒有擬南芥；其它4個(gè)亞組中4種植物均有LysM-RLKs家族成員分布。根據(jù)11個(gè)SgLysM-RLKs與GmLysM-RLKs的進(jìn)化關(guān)系，進(jìn)行一致的命名，如表2所示。

2.2SgLysM-RLKs家族成員的理化性質(zhì)分析及亞細(xì)胞定位預(yù)測

如表2所示，柱花草轉(zhuǎn)錄組SgLysM-RLK家族成員的理化性質(zhì)表明，其氨基酸數(shù)量為 個(gè)，相對(duì)分子質(zhì)量為 ，理論等電點(diǎn)為 ；約81. 82 % 的成員為酸性蛋白，不穩(wěn)定指數(shù)范圍為 ，僅SgLYKl，SgLYK5和SgLYK3b的不穩(wěn)定系數(shù)小于4O，為較穩(wěn)定的蛋白質(zhì)，其余蛋白質(zhì)（8個(gè)成員，占比 7 2 . 7 3 % 穩(wěn)定性較差；總平均親水性指數(shù)范圍為 僅 數(shù)值為正值，屬于疏水蛋白，其余的均為負(fù)值，屬于親水蛋白。亞細(xì)胞定位預(yù)測結(jié)果顯示，9個(gè)成員定位于膜上，1個(gè)成員定位于葉綠體，1個(gè)成員定位于細(xì)胞質(zhì)。

2.3SgLysM-RLKs家族成員的保守基序和結(jié)構(gòu)域分析

利用MEME在SgLysM-RLKs中共搜索10個(gè)保守基序（表3），結(jié)果（圖2）表明：SgLysM-RLKs保守基序的分布相似，motif1，motif3，motif4，motif5，motif7，motif8和motif10在所有成員中均存在，保守程度較高，其中motif1的保守性得分最高；SgNFR5中不存在motif2，SgLYK2a中不存在motif6和motif9，SgLYK3a和SgLYK3b中均不存在motif9，其中SgLYK2a的motif數(shù)量最少，僅8個(gè)。

結(jié)構(gòu)域分析共鑒定出7種保守結(jié)構(gòu)域，如圖3所示，所有SgLysM-RLKs都含有典型的胞外LysM結(jié)構(gòu)域、跨膜結(jié)構(gòu)域和胞內(nèi)激酶結(jié)構(gòu)域；定位在質(zhì)膜的成員均含有信號(hào)肽，而定位在葉綠體和細(xì)胞質(zhì)的成員則沒有信號(hào)肽。

2.4SgLysM-RLKs家族成員的蛋白質(zhì)二、三級(jí)結(jié)構(gòu)預(yù)測

蛋白二級(jí)結(jié)構(gòu)分析結(jié)果如表4所示：SgLysM-RLKs蛋白的二級(jí)結(jié)構(gòu)均含有 α 螺旋、β折疊、無規(guī)則卷曲和延伸鏈形式，其中 α 螺旋的占比為 3 3 . 1 4 % ～ 40 . 6 % ， β 折疊的占比為 2 . 9 9 % ～ 9 . 0 2 % ，無規(guī)卷曲的占比為 3 5 . 8 0 % ～ 4 7 . 7 2 % ，延伸鏈的占比為

1 4 . 1 7 % ～ 1 9 . 2 5 % 。因此，SgLysM-RLKs的二級(jí)結(jié)構(gòu)以 α 螺旋和無規(guī)卷曲為主，延伸鏈為輔。

為了觀察蛋白空間結(jié)構(gòu)的變化，使用AlphaFold2對(duì)SgLysM-RLKs蛋白的三級(jí)結(jié)構(gòu)進(jìn)行預(yù)測（圖4），結(jié)果顯示：所有成員的LysM結(jié)構(gòu)域均存在 α 螺旋和β折疊結(jié)構(gòu)，跨膜結(jié)構(gòu)均為 α 螺旋結(jié)構(gòu)，激酶結(jié)構(gòu)域均存在 α 螺旋和β折疊結(jié)構(gòu)。

2.5SgLysM-RLKs家族成員響應(yīng)炭疽菌及低磷脅迫的轉(zhuǎn)錄組學(xué)分析

SgLysM-RLKs家族成員在響應(yīng)炭疽菌及低磷脅迫處理下的轉(zhuǎn)錄組學(xué)分析結(jié)果如圖5所示。

由圖5A可知，在接種炭疽菌的處理下，S g L Y K I ， S g L Y K 3 a 和 S g L Y K 8 基因在 時(shí)表達(dá)量顯著下調(diào)， S g L Y K 4 a 基因的表達(dá)量在各個(gè)時(shí)間點(diǎn)無顯著差異， S g L Y K 5 基因在 時(shí)表達(dá)量顯著上調(diào)， 則顯著下調(diào)。

由圖5B和圖5C可見，在磷處理18d后，SgLysM-RLKs家族成員在根中均無顯著差異，在葉片中， S g L Y K 4 a 顯著下調(diào)， S g L Y K 4 b 顯著上調(diào)，剩余的9個(gè)成員在葉片中均無顯著差異。

2.6SgLysM-RLKs家族成員的組織特異性表達(dá)分析

SgLysM-RLKs家族成員在根、莖、葉組織的表達(dá)模式如圖6所示，SgLysM-RLKs家族成員在3種組織中均有表達(dá)， S g L Y K 4 b ， S g L Y K 5 ， S g L Y K 7 主要在根中表達(dá)，其中 S g L Y K 2 b ， S g L Y K 4 a ， S g L Y K 7 在莖中的表達(dá)量也相對(duì)較高，表明它們可能主要在根或莖中發(fā)揮功能； S g L Y K 2 a ， S g L Y K 3 a ， S g L Y K 3 b ， S g L Y K 8 S g N F R 5 主要在莖中表達(dá)，其中 S g L Y K 3 a 在葉中的表達(dá)量也相對(duì)較高，表明它們可能主要在莖或葉中發(fā)揮功能。

2.7SgLysM-RLKs家族成員響應(yīng)不同生物脅迫的表達(dá)模式分析

柱花草葉片SgLysM-RLKs響應(yīng)幾丁質(zhì)處理的表達(dá)模式分析結(jié)果如圖7所示。幾丁質(zhì)處理后，S g L Y K I ， S g L Y K 2 a ， S g L Y K 3 a ， S g L Y K 3 b ， S g L Y K 4 a ， S g L Y K 5 ， S g L Y K 7 ， S g L Y K 8 和 S g N F R 5 的表達(dá)量隨處理后時(shí)間的延長顯著上調(diào)（與 相比）；S g L Y K 2 b 的表達(dá)量先顯著下調(diào)后上調(diào)；而 S g L Y K 4 b 在任何時(shí)間點(diǎn)的表達(dá)量與 相比均無顯著差異。

柱花草葉片SgLysM-RLKs響應(yīng)炭疽菌DZ-19接種的表達(dá)模式分析結(jié)果如圖8所示。與 相比，S g L Y K 4 a ， S g L Y K 5 ， S g L Y K 7 ， S g L Y K 8 的表達(dá)量隨處理后時(shí)間的延長顯著上調(diào)； S g L Y K I 的表達(dá)量先顯著上調(diào)后下調(diào)； S g L Y K 3 b 的表達(dá)量先顯著下調(diào)后上調(diào)； S g L Y K 2 b ， S g L Y K 4 b ， S g L Y K 3 a ， S g N F R 5 的表達(dá)量顯著下調(diào)；而 S g L Y K 2 a （圖中未列出）在處理后不表達(dá)。

柱花草葉片SgLysM-RLKs響應(yīng)炭疽菌WC-02接種的表達(dá)模式分析結(jié)果如圖9所示。與 相比，S g L Y K 1 ， S g L Y K 4 a ， S g L Y K 5 ， S g L Y K 7 ， S g L Y K 8 的表達(dá)隨處理后時(shí)間的延長顯著上調(diào)； S g L Y K 3 a SgLYK3b的表達(dá)量先顯著上調(diào)后下調(diào)； S g L Y K 2 b SgLYK4b的表達(dá)量顯著下調(diào)； S g N F R 5 在任何時(shí)間點(diǎn)的表達(dá)量與0h相比無顯著差異；而 圖中未列出）在處理后不表達(dá)。另外，接種炭疽菌后柱花草的生長變化為：第3d柱花草葉尖稍顯發(fā)黃，第4d僅有少量葉片出現(xiàn)病斑，第7d時(shí)，病斑擴(kuò)散，逐漸出現(xiàn)掉葉癥狀。

綜上所述，在幾丁質(zhì)和兩種炭疽菌的處理中，表達(dá)量均上調(diào)的基因有5個(gè)： S g L Y K 1 ， S g L Y K 4 a S g L Y K 5 ， S g L Y K 7 ， S g L Y K 8 。表明這5個(gè)基因在柱花草抗病中可能起到重要作用。

2.8SgLysM-RLKs家族成員響應(yīng)低磷脅迫的表達(dá)模式分析

柱花草根部和葉片SgLysM-RLKs響應(yīng)低磷脅迫表達(dá)模式分析結(jié)果如圖10所示，進(jìn)行低磷處理10d后，根部組織只有 S g L Y K 1 ， S g L Y K 3 a 和 S g L Y K 5 的表達(dá)量顯著上調(diào)，其他家族成員的表達(dá)量無顯著變化甚至不表達(dá)；葉片組織中 S g L Y K 2 a 和S g L Y K 3 a 的表達(dá)量顯著上調(diào)， S g L Y K I 和 S g L Y K 4 a 的表達(dá)量顯著下調(diào)，其他家族成員的表達(dá)量無明顯變化甚至不表達(dá)。因此，只有 S g L Y K 3 a 基因在低磷處理1Od后的根和葉組織中顯著上調(diào)表達(dá)，表明其可能在低磷脅迫中起到重要作用。另外，低磷處理后柱花草的生長變化為：與正常供磷處理相比，低磷處理的柱花草地上部矮小、葉片數(shù)減少、分枝發(fā)育受阻，根系則表現(xiàn)為側(cè)根更長、密度更大。

Fig.5Transcriptome analysisof SgLysM-RLKs family in response to C gloeosporioidesand lowphosphorusstress 注：柱花草接種炭疽菌的轉(zhuǎn)錄組分析（A）；柱花草根部在低磷脅迫下的轉(zhuǎn)錄組分析（B）；柱花草葉片在低磷脅迫下的轉(zhuǎn)錄組分析（C）。High-Pi為 高磷 處理下的樣品；Low-Pi為低磷 處理下的樣品。橫坐標(biāo)的數(shù)字1，2，3代表每個(gè)生物學(xué)重復(fù)標(biāo)記 Note：Transcriptome analysis of S. guianensis in response to ：gloeosporioides（A）；Transcriptome analysis of stylo’sroots under low phosphorus stress（B）；Transcriptome analysis ofstylo’sleavesunderlowphosphorus stress（C）.High-Pirefers tosampleunderHP treatment；Low-Pi refers to sample under LP treatment. The numbers 1，2 and 3 in the absciss represent each biological replicate

3 討論與結(jié)論

3.1SgLysM-RLKs家族鑒定及生物信息學(xué)分析

由于尚無柱花草的基因組發(fā)表，因此本研究利用柱花草轉(zhuǎn)錄組數(shù)據(jù)鑒定其LysM-RLKs家族成員，首次得到11個(gè)成員，并進(jìn)行了命名。雖然可能未得到完整的LysM-RLKs家族成員，但鑒定到的成員數(shù)量多于擬南芥（5個(gè)）]，說明LysM-RLKs家族在柱花草中發(fā)生了明顯擴(kuò)增，且這些成員可能增加了新的功能。已有研究報(bào)道，亞組I中的AtLYK1/AtCERK1，LjCERK6和亞組V中的

AtLYK5具有識(shí)別幾丁質(zhì)功能[4.25-26]，參與植物免疫反應(yīng)；亞組I中的LjNFR1，GmNFRlα， 和亞組ⅢI中的LjNFR5具有識(shí)別Nod因子（NFs）功能[27-28]，參與根瘤或叢枝菌根共生。此外，在系統(tǒng)發(fā)育樹上，柱花草、百脈根和大豆的成員均能聚類到同一小分支下，說明同為豆科植物的它們親緣關(guān)系相近，在進(jìn)化過程中相對(duì)保守，功能類似的可能性也更大。據(jù)此推測，亞組I中的SgLYK1可能具有識(shí)別幾丁質(zhì)和Nod因子的功能，亞組V中的SgLYK5可能具有識(shí)別幾丁質(zhì)功能，亞組Ⅲ中的SgNFR5可能具有識(shí)別Nod因子功能。

注：結(jié)果代表3次實(shí)驗(yàn)的平均值，誤差線表示3次生物學(xué)重復(fù)表達(dá)值的標(biāo)準(zhǔn)誤差。使用 T -test方差分析比較不同組織之間的差異顯著性， 0.05， 和

Note：Theulsteexptssatedrfpefhlepl cates.Significantdifference between different tissues were compared using T -test analysis of variance， ， and

Fig.7Analysis of the expresson pattern ofSgLysM-RLKs in response to chitin treatment in stylo’s leaves注：結(jié)果代表3次實(shí)驗(yàn)的平均值，誤差線表示3次生物學(xué)重復(fù)表達(dá)值的標(biāo)準(zhǔn)誤差。使用 T -test方差分析比較0h和其他處理時(shí)間之間的差異顯著性， 5 和 ，下同

Note：TheresultsepreentteragefromtheexpemetsErrorbarsidicatestadardoofthexpressoaluesfromthological replicates. Significant difference between Oh and other times were compared using T -test analysis of variance， ， and 0.001，similarlyhereinafter

Fig.8AnalysisofteexpressonpatteofSgLysM-RLKsinesponsetoinoculationwithColetotrichumgloeosporiodesstrair DZ-19 in stylo'sleaves

發(fā)揮模式識(shí)別受體功能的LysM-RLKs蛋白絕大部分定位在質(zhì)膜上，例如蘋果（Malusdomestica）MdCERK1[22]、葡萄（Vitis vinifera）VvLYK5-1[29]和煙草（NicotianatabacumL.）NbLYK4.2[30]等植物的LysM-RLKs成員均已證實(shí)定位在質(zhì)膜上。本研究亞細(xì)胞定位預(yù)測分析表明，柱花草大多數(shù)LysM-

RLKs蛋白定位在質(zhì)膜上，說明LysM-RLKs蛋白主要在質(zhì)膜上發(fā)揮功能，這與葡萄、蘋果、梨（Pyrusbretschneideri）和獼猴桃（Actinidiachinensis）的亞細(xì)胞定位預(yù)測結(jié)果類似[31-33]。本研究表明，SgLysM-RLKs家族成員的多個(gè)保守基序和結(jié)構(gòu)域的排序相似，且都含有典型的胞外LysM結(jié)構(gòu)域、跨膜結(jié)構(gòu)域和胞內(nèi)激酶結(jié)構(gòu)域（圖2，圖3），數(shù)量和分布也高度相似，因此，推測LysM-RLKs家族成員在柱花草響應(yīng)炭疽菌侵染和磷脅迫的過程中發(fā)揮重要的激發(fā)、轉(zhuǎn)導(dǎo)和調(diào)控作用。

SgLYK1 SgLYK2a SgLYK2b SgLYK3a 皖 主 * 十 3 華 1.0- 5 工 工 2↓ 王 0.5 士 不表達(dá)不表達(dá)不表達(dá) 05 王 0 不表達(dá) 0 expesedpsee HP-R LP-R HP-R LP-L HP-R LP-R HP-R LP-L HP-R LP-R HP-R LP-L HP-RLP-RHP-RLP-L SgLYK3b SgLYK4a SgLYK4b SgLYK5 1.5- T Barsgeseeeneree 1.5 1.0 王 華 1.0- 工 士 “ 1.0 士 擊 No不表達(dá) Note 不表達(dá) 5 不表達(dá) □ HP-RLP-RHP-RLP-L HP-RLP-R HP-R LP-L HP-RLP-RHP-RLP-L HP-RLP-RHP-RLP-L SgLYK7 SgLYK8 SgNFR5 出 1.0 0.5- 0.5 0.5 No不表達(dá) ot不表達(dá) 0 0 0 HP-RLP-RHP-RLP-L HP-RLP-RHP-R LP-L HP-RLP-RHP-R LP-L

stressinstylo’srootsand leaves 注： 為高磷 處理下的根部/葉片樣品； 為低磷 處理下的根部/葉片樣品。結(jié)果 代表3次實(shí)驗(yàn)的平均值，誤差線表示3次生物學(xué)重復(fù)表達(dá)值的標(biāo)準(zhǔn)誤差。使用 T -test方差分析比較高磷和低磷處理之間的差異顯著性， 0.05， 和 （204 Note： refers to roots/leaves under HP treatment；LP-R/L refers to roots/leaves under LP （204 ）treatment.Theresultsrepresent theaveragefromthreexperimentsErrorbarsindicatestandarderoroftheexpressionvaluesfrom three biological replicates.Significant diference between highand low P treatment were compared using T -test analysis of variance， ， and （20

圖4顯示，SgLysM-RLKs家族成員的蛋白三級(jí)結(jié)構(gòu)存在差異，特別是SgLYK7與其他成員表現(xiàn)出了明顯差異。據(jù)報(bào)道，植物L(fēng)ysM-RLKs家族成員除參與免疫和共生外，還參與細(xì)胞死亡、耐鹽性、脫落酸信號(hào)傳導(dǎo)等[31.34]，而SgLYK7只在低磷處理的轉(zhuǎn)錄組中被鑒定到，因此推測SgLysM-RLKs家族成員在長期進(jìn)化的過程中產(chǎn)生了明顯的功能分化。另外，SgLYK1和AtLYK1/AtCERK1的三級(jí)結(jié)構(gòu)類似，因而可能具有AtLYK1/AtCERK1的類似功能，這進(jìn)一步加大了SgLYK1具有識(shí)別幾丁質(zhì)功能的可能性。

3.2SgLysM-RLKs家族參與調(diào)控抗病功能

LysM-RLKs是一類在植物抗病中發(fā)揮重要作用的基因家族，在很多植物中已得到證實(shí)，例如擬南芥[35]、水稻[36]、番茄（Solanum lycopersicum）[37]和蘋果[22]等，同時(shí)，基因的表達(dá)情況與其功能密切相關(guān)。本研究表明，在SgLysM-RLKs家族響應(yīng)幾丁質(zhì)和炭疽菌的過程中，RT-qPCR結(jié)果與轉(zhuǎn)錄組數(shù)據(jù)存在一定差異：在接種炭疽菌的轉(zhuǎn)錄組數(shù)據(jù)中，S g L Y K I 基因在 時(shí)表達(dá)量顯著下調(diào)，而RT-qPCR結(jié)果則在 或 時(shí)表達(dá)量已經(jīng)顯著上調(diào)， 后表達(dá)量已無顯著差異。這可能是由于LysM-RLKs家族成員響應(yīng)病原體侵染主要在前期，后期停止過度表達(dá)以避免過度的免疫反應(yīng)而損傷植物本身。由于轉(zhuǎn)錄組數(shù)據(jù)取樣時(shí)間在 后，此時(shí)已處于響應(yīng)后期，因此可能在轉(zhuǎn)錄組的時(shí)段內(nèi)基因表達(dá)無顯著差異。同理，這也可能是 S g L Y K 4 a 和 S g L Y K 8 基因轉(zhuǎn)錄組數(shù)據(jù)與RT-qPCR結(jié)果存在差異的主要原因。

綜上所述，在幾丁質(zhì)和炭疽菌處理中均上調(diào)的基因SgLYK1， S g L Y K 4 a ，SgLYK5， S g L Y K 7 和S g L Y K 8 在柱花草抗病中可能起到重要作用。

3.3SgLysM-RLKs家族參與調(diào)控低磷脅迫功能

本研究表明，在低磷脅迫處理的轉(zhuǎn)錄組數(shù)據(jù)中，SgLysM-RLKs家族成員在根中的表達(dá)量均無顯著差異， S g L Y K 4 a 和 S g L Y K 4 b 在葉片中的表達(dá)量分別顯著下調(diào)和上調(diào)。RT-qPCR結(jié)果顯示S g L Y K I ， S g L Y K 3 a 和 S g L Y K 5 在根中的表達(dá)量顯著上調(diào)；葉片組織中 S g L Y K 2 a 和 S g L Y K 3 a 的表達(dá)量顯著上調(diào)， S g L Y K I 和 S g L Y K 4 a 的表達(dá)量顯著下調(diào)，即只有 S g L Y K 4 a 符合轉(zhuǎn)錄組的表達(dá)趨勢，其他基因存在差異，這可能是處理時(shí)間點(diǎn)的不同所致。總體上，SgLysM-RLKs家族成員在磷脅迫中也存在差異表達(dá)基因，因此部分SgLysM-RLKs家族成員可能在低磷脅迫中起到作用，且可能參與調(diào)控柱花草響應(yīng)低磷脅迫，但其具體分子機(jī)制有待進(jìn)一步研究。

綜上所述，本研究利用轉(zhuǎn)錄組數(shù)據(jù)和生物信息學(xué)方法共鑒定到11個(gè)LysM-RLKs家族成員 （ S g L Y K I S g L Y K 2 a ， S g L Y K 2 b ， S g L Y K 3 a ， S g L Y K 3 b ， S g L Y K 4 a ， S g L Y K 4 b ， S g L Y K 5 ， S g L Y K 7 ， S g L Y K 8 和 S g N F R 5 ），其蛋白大多數(shù)定位于質(zhì)膜上，含有典型的胞外LysM結(jié)構(gòu)域、跨膜結(jié)構(gòu)域和胞內(nèi)激酶結(jié)構(gòu)域，蛋白二級(jí)結(jié)構(gòu)以 α 螺旋和無規(guī)卷曲為主，蛋白三維結(jié)構(gòu)方面存在差異。表達(dá)模式分析表明，該家族成員均存在組織表達(dá)特異性，其中 S g L Y K I ， S g L Y K 4 a S g L Y K 5 ， S g L Y K 7 和 S g L Y K 8 可能在響應(yīng)炭疽菌侵染和幾丁質(zhì)誘導(dǎo)的免疫反應(yīng)中發(fā)揮重要功能，S g L Y K 3 a 則可能參與調(diào)控柱花草響應(yīng)低磷脅迫。

參考文獻(xiàn)

[1]YANG C，WANG E T，LIU J. CERK1，more than a co-receptorin plant-microbe interactions[J].New Phytologist，2022， 234（5）：1606-1613

[2] PETUTSCHNIGEK，JONESAME，SERAZETDINOVA L，etal.The lysin motif receptor-like kinase（LysM-RLK） CERK1 isamajorchitin-bindingprotein inArabidopsis thalianaand subjectto chitin-induced phosphorylation[J].Journal ofBiological Chemistry，2010，285（37）：28902-28911

[3]BROWN H E，ESHER S K，ANDREW ALSPAUGH J. Chitin：a“hidden figure”in the fungal cell wall[J].Current TopicsinMicrobiologyand Immunology，2020，425：83-111

[4]MIYAA，ALBERTP，SHINYAT，etal.CERK1，aLysM receptorkinase，isessential forchitinelicitorsignaling inArabidopsis[J].Proceedingsof theNational Academy of Sciencesof theUnitedStates of America，2007，104（49）：19613-19618

[5]YAO KJ，WANGYL，LIX，etal.Genome-wide identificationof thesoybeanLysM-RLK family genesand itsnitrogen response[J].InternationalJournalofMolecularSciences，2023， 24（17）：13621

[6]RUMAN H，KAWAHARADA Y.A new classification of lysinmotif receptor-likekinasesinLotusjaponicus[J].Plantamp; CellPhysiology，2023，64（2）：176-190

[7]ZHANG X C，WU X L，F(xiàn)INDLEY S，et al. Molecular evolution of lysin motif-type receptor-like kinases in plants[J].Plant Physiology，2007，144（2）：623-636

[8]ZHANG C，HE JM，DAI HL，et al. Discriminating symbiosis and immunity signals by receptor competition in rice[J].Proceedings of the National Academy of Sciences of the United States of America，2021，118（16）：e2023738118

[9]GIBELIN-VIALA C，AMBLARD E，PUECH-PAGES V，et al. The Medicago truncatula LysM receptor-like kinase LYK9 plays a dual role in immunity and thearbuscular mycorrhizal symbiosis[J].The New Phytologist，2019，223（3）：1516-1529

[10]易克賢，鄭金龍，習(xí)金根，等．柱花草炭疽病流行因子初探[J]. 熱帶農(nóng)業(yè)科學(xué)，2014，34（2）：74-78

[11]CHANDRA A，PATHAK P，BHATT R. Stylosanthes research in India： Prospects and challenges ahead[J]. Current Science，2006，90：915-921

[12]DE SOUZAJM，DA SILVA FRAGOSO KN，ORILIO AF， et al. Two viruses from Stylosanthes guianensis may represent a new genus within Potyviridae[J].Virus Research，2021，293： 198257

[13］劉國道，白昌軍，何華玄，等．熱研5號(hào)柱花草選育研究[J]．草 地學(xué)報(bào)，2001，9（1）：1-7

[14]易克賢，鄭金龍，習(xí)金根，等．中國柱花草炭疽病病原菌遺傳 多態(tài)性的AFLP分析[J]．熱帶作物學(xué)報(bào)，2014，35（5）： 909-917

[15］鄭麗，易小平，文衍堂，等．省熱帶牧草病害調(diào)查及病原 菌初步鑒定[J].熱帶作物學(xué)報(bào)，2014，35（5）：967-973

[16] CHEN H L，HU C P，YI K X，et al. Cloning of insertion site flanking sequence and construction of transfer DNA insert mutant library in Stylosanthes Colletotrichum[J].PLoS One， 2014，9（10）：e111172

[17]DU Y M，TIAN J，LIAO H，et al. Aluminium tolerance and high phosphorus efficiency helps Stylosanthes better adapt to low-Pacid soils[J].Annals of Botany，2009，103（8）：1239- 1247

[18] ZHANG S Z，XU YF，WANG F，et al. Transcriptomic and physiological analysis of the effects of exogenous phloretin and pterostilbene on resistance responses of Stylosanthes against anthracnose[J]. International Journal of Molecular Sciences， 2024，25（5）：2701

[19]FAMOSO A N，CLARK R T，SHAFF JE，et al. Development of a novel aluminum tolerance phenotyping platform used forcomparisons of cereal aluminum tolerance and investigations into rice aluminum tolerance mechanisms[J]. Plant Physiology， 2010，153（4）：1678-1691

[20] JIANG L Y，WU PP，YANG L Y，et al. Transcriptomics and metabolomics reveal the induction of flavonoid biosynthesis pathway in the interaction of Stylosanthes-Colltotrichum gloeosporioides[J]. Genomics，2021，113（4）：2702-2716

[21] LUO J J， CHEN Z J，HUANG R，et al. Multi-omics analysis reveals the roles of purple acid phosphatases in organic phosphorus utilization by the tropical legume Stylosanthes guianensis[J]. ThePlantJournal，2024，117（3）：729-746

[22]高靜，楊麗云，張世子，等．柱花草葉肉原生質(zhì)體提取與瞬時(shí) 表達(dá)體系構(gòu)建[J].草地學(xué)報(bào)，2023，31（3）：893-902

[23］高夢(mèng)澤．柱花草膠孢炭疽菌與擬南芥互作體系的建立及 SgPALs的初步功能探究[D].：大學(xué)，2021：22

[24]LIUPD，XUEYB，CHENZJ，etal.Characterizationof purpleacid phosphatases involved in extracellular dNTP utilizationinStylosanthes[J].Journalof Experimental Botany，2016， 67（14）：4141-4154

[25]BOZSOKI Z，CHENG J，F(xiàn)ENG F，et al. Receptor-mediated chitin perception in legume roots is functionally separable from Nod factor perception[J]. Proceedings of the National Academyof Sciencesof theUnited StatesofAmerica，2Ol7，114 （38）：E8118-E8127

[26］閆一嘉.擬南芥LYKs受體激酶介導(dǎo)幾丁質(zhì)激發(fā)免疫反應(yīng)的 分子機(jī)制研究[D].武漢：華中農(nóng)業(yè)大學(xué)，2020：49-64

[27]RADUTOIUS，MADSENLH，MADSENEB，etal.Plant recognition of symbiotic bacteriarequirestwo LysM receptorlikekinases[J].Nature，2003，425（6958）：585-592

[28]INDRASUMUNARA，SEARLEI，LINMH，etal.Nodulation factor receptor kinase lα controls nodule organ number in soybean（GlycinemaxL.Merr）[J].ThePlant Journal，2011， 65（1）：39-50

[29]ROUDAIRE T，MARZARI T，LANDRY D，et al. The grapevineLysMreceptor-like kinase VvLYK5-1 recognizeschitin oligomers through itsassociation with VvLYK1-1[J]. Frontiers inPlant Science，2023，14：1130782

[30]WANG T M，GASCIOLLI V，GASTON M，et al. LysM receptor-like kinases involved inimmunity perceive lipo-chitooligosaccharides in mycotrophic plants[J].Plant Physiology， 2023，192（2）：1435-1448

[31]褚明宇，李婉瑩，左存武，等.葡萄LysM類受體激酶基因家族 的鑒定及表達(dá)分析[J].農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，2023，31（2）： 282-297

[32］王東東，溫有福，趙丹，等．蘋果和梨LysM-RLK響應(yīng)腐爛病 信號(hào)成員鑒定[J].植物生理學(xué)報(bào)，2022，58（2）：381-392

[33]何蕾，王文霞，胡容平，等.獼猴桃LysM型類受體激酶基因家 族生物信息學(xué)分析[J].分子植物育種，2023，21（10）：3191- 3198

[34]HUSP，LIJJ，DHARN，et al.Lysin motif（LysM）proteins： interlinking manipulation of plant immunity and fungi[J].InternationalJournal ofMolecular Sciences，2O21，22（6）：3114

[35]WANJR，ZHANGXC，NEECED，etal.ALysMreceptorlikekinase playsa critical role in chitin signaling and fungal resistanceinArabidopsis[J].ThePlant Cell，2Oo8，20（2）： 471-481

[36]SHIMIZUT，NAKANOT，TAKAMIZAWAD，etal.Two LysM receptor molecules，CEBiPand OsCERK1，cooperatively regulatechitin elicitorsignaling in rice[J].ThePlant Journal， 2010，64（2）：204-214

[37]ZENGLR，VELASQUEZAC，MUNKVOLDKR，etal.A tomato LysM receptor-like kinase promotesimmunity and its kinase activityis inhibited by AvrPtoB[J].ThePlantJournal， 2012，69（1）：92-103

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