毛仁俊 白朕卿 何志貴 李秋月 吳佳文 陳浩 閻巖

摘 要 為了解析含羞草中黃酮類物質(zhì)的生物合成途徑，利用Illumina platform 2000TM測(cè)序平臺(tái)對(duì)含羞草葉片進(jìn)行轉(zhuǎn)錄組測(cè)序，共獲得94 182個(gè)Unigene，平均長(zhǎng)度為695 bp，N50值為1 159 bp。共有30 243個(gè)Unigene注釋于50個(gè)GO功能組中，其中“代謝過程”“催化活性”以及“細(xì)胞”注釋的Unigene數(shù)量較多。KEGG通路分析鑒定出49個(gè)Unigene注釋在黃酮生物合成途徑，分別編碼查爾酮合成酶（CHS，12個(gè)Unigene）、查爾酮異構(gòu)酶（CHI，4）、黃烷酮-3-羥化酶（F3H，3）、黃烷酮-3-羥化酶（F3H，9）、黃烷酮-3，5-羥化酶（F3，5H，1）、二氫黃酮醇4-還原酶（DFR，5）、黃酮醇合成酶（FLS，3）以及無色花色素還原酶（LAR，12）基因。從轉(zhuǎn)錄組數(shù)據(jù)庫中鑒定出7 382個(gè)以單核苷酸和三核苷酸為主要類型的SSR標(biāo)記。隨機(jī)選出10對(duì)SSR引物進(jìn)行擴(kuò)增，其中有7對(duì)能夠擴(kuò)增出清晰的條帶。

關(guān)鍵詞 含羞草; 轉(zhuǎn)錄組測(cè)序; 代謝途徑; 黃酮生物合成

含羞草（Mimosa pudica Linn）是豆科含羞草屬的一種多年生草本或小灌木，主要分布于熱帶和亞熱帶地區(qū)，具有較高的觀賞價(jià)值。含羞草以全草入藥，在印度、孟加拉國(guó)、菲律賓以及中國(guó)廣西、云南等少數(shù)民族地區(qū)是一種常用藥材[1-2]。中國(guó)古代醫(yī)學(xué)典籍《生草藥性備要》、《本草求原》以及《嶺南采藥錄》中均記載了含羞草味甘、性寒，具有清熱解毒、消腫止痛的功效?，F(xiàn)代藥理學(xué)研究表明含羞草具有愈合創(chuàng)傷[3]、抗菌消炎[4]、安神止痛[5]、降低血壓[6]以及緩解糖尿病病癥[7]等多種功效。此外，含羞草作為觸感性研究的模式植物，具有顯著的觸感特性，即在受到強(qiáng)光、震動(dòng)、風(fēng)吹、碰觸等外界刺激時(shí)會(huì)迅速閉合葉片、收緊枝條（圖1）。

含羞草中含有兒茶素、槲皮素、豆甾醇以及對(duì)香豆酸等黃酮類物質(zhì)，也含有生物堿、皂苷以及酚酸類化合物[8-9]，這些化合物是其發(fā)揮藥效功能的物質(zhì)基礎(chǔ)。Zhang等[1]通過核磁共振氫譜（1H-NMR）、核磁共振碳譜（13C-NMR）以及質(zhì)譜（MS）的方法從含羞草中分離出異紅草素、異牡荊素及牡荊素等黃酮類物質(zhì)，并對(duì)其活性進(jìn)行分析。Yuan等[10]利用紅外光譜（IR）、質(zhì)譜（MS）以及核磁共振（NMR）的方法首次從含羞草中分離出黃酮類化合物6，7，3′，4′-四羥基-8-C-［α-L-鼠李糖-（1→2）］-β-D-葡糖黃酮碳苷和5，7，3′，4′-四羥基-8-C-［β-D-芹菜糖-（1→4）］-β-D-葡糖黃酮碳苷。喬文濤等[11]通過正交分析法對(duì)含羞草中黃酮類物質(zhì)的研究表明，以10倍量的80%乙醇提取，利用大孔吸附樹脂進(jìn)行分離，再以不同濃度的乙醇進(jìn)行梯度洗脫，最后經(jīng)過濃縮得到高濃度的黃酮。

目前，關(guān)于含羞草的研究主要集中于藥效功能[1，12]、觸感性效應(yīng)[13-14]以及生理活性等方面[15-16]，而分子水平的研究匱乏，組學(xué)研究鮮見報(bào)道。早期對(duì)于含羞草觸感性的研究已經(jīng)證明含羞草枝條基部的葉枕是調(diào)節(jié)其葉片開合的主要器官[17]，也鑒定出鉀-L-蘋果酸，鎂鉀反式烏頭酸等與其觸感性相關(guān)的化合物[13-14]。然而，近年來觸感性相關(guān)研究進(jìn)展緩慢，制約其發(fā)展的重要因素是分子水平研究不足，次生代謝物合成途徑尚未解析。因此，本研究利用轉(zhuǎn)錄組測(cè)序技術(shù)對(duì)含羞草葉片進(jìn)行分析，解析黃酮生物合成途徑，注釋途徑中的關(guān)鍵酶基因，以期為深入研究黃酮途徑中關(guān)鍵酶基因的功能以及進(jìn)一步開發(fā)含羞草藥用價(jià)值提供參考。

1 材料與方法

1.1 植物材料

含羞草樣品于2019年6月采集于廣西藥用植物園。試驗(yàn)樣品經(jīng)過延安大學(xué)生命科學(xué)學(xué)院閻巖教授鑒定為豆科（Leguminosae）含羞草屬（Mimosa Linn）植物含羞草（Mimosa pudica L.）。采集6株長(zhǎng)勢(shì)良好且無病蟲害的含羞草植株地上部分，選取葉片用于后續(xù)試驗(yàn)分析。將含羞草葉片分為兩份，一份用錫箔紙包裹后立即置于液氮中速凍，隨后置于超低溫冰箱保存，另一份樣品陰干后用于黃酮類物質(zhì)的測(cè)定。

1.2 總黃酮和兒茶素含量的測(cè)定

以亞硝酸鈉-氯化鋁-氫氧化鈉顯色法測(cè)定含羞草葉片中總黃酮的含量[18]。精確稱取過40目篩的含羞草葉片粉末1.0 g，置于具塞燒瓶中，加入25 mL體積分?jǐn)?shù)為 70%的乙醇，稱量。加熱回流1 h，冷卻后再次稱量，并用70%乙醇補(bǔ)足質(zhì)量，過濾后待用。取1 mL樣品溶液置于25 mL容量瓶中，加入6 mL超純水和1 mL 5%的亞硝酸鈉溶液，混勻后靜置6 min，再加入10% 硝酸鋁溶液1 mL，加入10 mL 4.3%氫氧化鈉，用70%乙醇定容，靜置15 min。利用分光光度計(jì)（UV-2450，Shimadzu）測(cè)定波長(zhǎng)510 nm處的吸光度值，以蘆?。兌?98%，Sigma，USA）為對(duì)照品制作標(biāo)準(zhǔn)曲線計(jì)算總黃酮含量。以香莢蘭素比色法[19]測(cè)定兒茶素的含量。取含羞草葉片粉末1.0 g置于具塞燒瓶中，加入25 mL 70%甲醇，稱量。加熱回流0.5 h，放冷后再次稱量，用70%甲醇補(bǔ)足質(zhì)量，過濾，取續(xù)濾液待測(cè)。以兒茶素（純度>98%，Sigma，USA）為對(duì)照品繪制標(biāo)準(zhǔn)曲線計(jì)算兒茶素含量。

1.3 RNA的提取和測(cè)序

利用RNeasy Plant Mini試劑盒（Qiagen，Hilden，Germany）提取含羞草葉片的總RNA，提取方法參照試劑盒說明。通過Nanodrop核酸蛋白檢測(cè)儀（Thermo Scientific，MA，USA）和Qubit 2.0核酸定量?jī)x（Life Technologies，CA，USA）分析樣品RNA的質(zhì)量。cDNA文庫的構(gòu)建參照RNA-Seq文庫構(gòu)建試劑盒（Illumina，USA）說明書進(jìn)行。測(cè)序工作由北京百邁客生物科技有限公司利用Illumina HiSeq2000TM測(cè)序平臺(tái)完成。

1.4 轉(zhuǎn)錄組數(shù)據(jù)的拼接、組裝及功能注釋

剔除測(cè)序結(jié)果中的接頭讀序、低質(zhì)量讀序以及包含ploy-N的讀序，獲得高質(zhì)量讀序。利用Trinity軟件[20]對(duì)序列進(jìn)行組裝，獲得Unigene和轉(zhuǎn)錄本（Transcript）。利用Blast2GO進(jìn)行GO功能注釋[21]，利用KOBAS 2.0進(jìn)行KEGG代謝通路注釋[22]。通過BLASTX將Unigene分別與Nr（NCBI非冗余蛋白數(shù)據(jù)庫）、Swiss-Prot（Swiss-蛋白序列數(shù)據(jù)庫）、eggNOG（基因進(jìn)化譜系：無監(jiān)督的同源群體）、KEGG（京都基因和基因組百科全書）和KOG數(shù)據(jù)庫（蛋白直系同源數(shù)據(jù)庫）進(jìn)行相似性比對(duì)，獲取Unigene的功能注釋和分類信息。

1.5 SSR標(biāo)記的開發(fā)與鑒定

使用MicroSatellite（http：//pgrc.ipk-gatersleben.de/misa/）軟件對(duì)轉(zhuǎn)錄組中的序列長(zhǎng)度大于1 kb的Unigene進(jìn)行串聯(lián)重復(fù)位點(diǎn)檢測(cè)。以單核苷酸重復(fù)不少于10次、二核苷酸重復(fù)不少于6次、三核苷酸、四核苷酸、五核苷酸、六核苷酸不少于4次為標(biāo)準(zhǔn)，統(tǒng)計(jì)分析SSR分子標(biāo)記的種類與數(shù)量[23]。從轉(zhuǎn)錄組開發(fā)出的SSR數(shù)據(jù)庫中隨機(jī)選出10對(duì)引物進(jìn)行PCR擴(kuò)增驗(yàn)證，引物序列信息見表1。

2 結(jié)果與分析

2.1 總黃酮與兒茶素含量測(cè)定

利用分光光度法測(cè)定含羞草葉片中總黃酮和兒茶素的含量。以不同濃度蘆丁為橫坐標(biāo)，對(duì)應(yīng)的吸光度值為縱坐標(biāo)繪制標(biāo)準(zhǔn)曲線，得到總黃酮標(biāo)準(zhǔn)曲線為Y=0.998 5X+0.008 3，R2=0.999 7，線性范圍0～128.9 mg/g（圖2-A）。測(cè)定結(jié)果表明3份含羞草樣品中總黃酮含量為75.71 mg/g±0.42 mg/g；以同樣方法獲得兒茶素標(biāo)準(zhǔn)曲線為Y=1.178 4X-0.002 4，R2=0.999 8，線性范圍0～1.1 mg/g，3份樣品中兒茶素含量為0.223 mg/g±0.012 mg/g（圖2-B）。

2.2 序列組裝

3個(gè)樣本共獲得67 400 933個(gè)高質(zhì)量讀序，對(duì)比成功的讀序數(shù)量為55 702 861，3個(gè)樣本中比對(duì)成功的讀序比例均超過82%（表2）。本次測(cè)序總數(shù)據(jù)量為20.08 GB。拼接后共獲得211 405個(gè)轉(zhuǎn)錄本，平均長(zhǎng)度為1 178 bp，N50值為1 962 bp。共獲得94 182個(gè)Unigene，平均長(zhǎng)度為695 bp，N50值為1 159 bp。樣本的平均GC含量為45.86%，Q30值均大于93.35%（表2）。測(cè)序錯(cuò)誤率小于0.1%的堿基數(shù)量百分比為92.48%。Unigene長(zhǎng)度分布結(jié)果表明，長(zhǎng)度為1 000～2 000 bp和300～500 bp的Unigene占比較大，分別為37.30%（35 134個(gè)）和11.46%（10 795個(gè)）。這些指標(biāo)表明本次測(cè)序結(jié)果良好，得到大量序列信息，能夠滿足后續(xù)生物信息學(xué)分析的要求。轉(zhuǎn)錄組數(shù)據(jù)已提交至美國(guó)國(guó)家生物技術(shù)信息中心數(shù)據(jù)庫，獲得序列號(hào)PRJNA634474。

2.3 Unigene的功能注釋

以E-value≤e-5且HMMER E-value≤e-10為標(biāo)準(zhǔn)進(jìn)行BLAST，共獲得47 090個(gè)有注釋信息的Unigene。進(jìn)一步分析表明分別有13 921、30 243、19 427、26 970、29 902、33 103、44 256和46 480個(gè)Unigene注釋在COG、GO、KEGG、KOG、Pfam、Swissprot、eggNOG和Nr數(shù)據(jù)庫中。注釋到長(zhǎng)度大于1 kb的Unigene共有15 775個(gè)，占所有注釋Unigene的33.5%（表3）。Nr注釋的物種分類表明，比對(duì)到木豆（Cajanus cajan）的Unigene最多，為3 207個(gè)，占比6.90%，隨后為大豆（Glycine max）3 147個(gè)、橡膠樹（Hevea brasiliensis）2 709個(gè)、羽扇豆（Lupinus angustifolius）2 608個(gè)以及麻風(fēng)樹（Jatropha curcas）2 144個(gè)（圖3）。

2.4 GO功能分類

按照GO功能分類，30 243個(gè)Unigene可分入細(xì)胞組成（CC）、分子功能（MF）和生物過程（BP）3大類，50個(gè)功能組（圖4）。細(xì)胞組成中的Unigene可進(jìn)一步分為15個(gè)功能組，其中以“細(xì)胞”（3 887個(gè)Unigene）和“細(xì)胞部分”（13 813）數(shù)量較多；分子功能中的Unigene同樣被分入15個(gè)功能組，主要注釋在“催化活性”（16 651）和“結(jié)合”（14 544）；在生物過程中，Unigene被分入20個(gè)功能組，以“代謝過程”（15 764）和“細(xì)胞過程”（1 416）中的Unigene數(shù)量較多（圖4）。

在KOG數(shù)據(jù)庫中比對(duì)Unigene信息，預(yù)測(cè) Unigene功能并進(jìn)行分類統(tǒng)計(jì)。結(jié)果表明共有94 182個(gè)Unigene獲得注釋信息，可分為25類（圖5）。其中“一般功能預(yù)測(cè)（R）”注釋的Unigene數(shù)量最多，為6 252個(gè)，占Unigene總數(shù)的23.18%。隨后依次為“翻譯后修飾、蛋白質(zhì)轉(zhuǎn)換、分子伴侶（O）”（2 979，11.05%）、“信號(hào)轉(zhuǎn)導(dǎo)機(jī)制（T）”? ?（2 689，9.97%）、“翻譯-核糖體結(jié)構(gòu)和生物發(fā)生（J）”（1 385，5.14%），而“細(xì)胞活動(dòng)性”（N）（9，0.03%）、“細(xì)胞外部結(jié)構(gòu)（W）”（89，0.33%）及“核結(jié)構(gòu)（Y）”（113，0.42%）注釋到的Unigene數(shù)量較少（圖5）。

2.5 KEGG途徑分類注釋

為進(jìn)一步研究含羞草的代謝途徑，將獲得功能注釋的Unigene進(jìn)行KEGG通路分析。結(jié)果表明共有9 155個(gè)Unigene注釋在“新陳代謝”、“遺傳信息處理”、“環(huán)境信息處理”、“細(xì)胞過程”以及“有機(jī)系統(tǒng)”5個(gè)大類、28個(gè)小類（圖6）。注釋在“新陳代謝”中的Unigene數(shù)量最多，共有3 885個(gè)，其中又以“碳代謝（700）、氨基酸生物合成（619）”、“淀粉及糖類代謝（508）”注釋到的Unigene數(shù)量較多。共有3 243個(gè)Unigene注釋在“遺傳信息處理”，其中有603個(gè)Unigene注釋在“核糖體”，數(shù)量最多，“內(nèi)質(zhì)網(wǎng)蛋白加工”中注釋的Unigene次之（482個(gè)）?！碍h(huán)境信息處理”中注釋到的Unigene較少，其中“ABC轉(zhuǎn)運(yùn)器”途徑僅有92個(gè)Unigene（圖6）。

2.6 含羞草黃酮生物合成途徑的解析

黃酮是含羞草中重要的藥效成分，本研究重點(diǎn)關(guān)注含羞草黃酮生物合成途徑（ko00941）。注釋分析的結(jié)果表明共有9 155個(gè)Unigene注釋在136條KEGG途徑，其中49個(gè)Unigene注釋在黃酮生物合成途徑，分別編碼查爾酮合成酶（CHS）、查爾酮異構(gòu)酶（CHI）、黃烷酮-3-羥化酶（F3H）、黃烷酮-3-羥化酶（F3H）、黃烷酮-3，5-羥化酶（F3，5H）、二氫黃酮醇4-還原酶（DFR）、黃酮醇合成酶（FLS）以及無色花色素還原酶（LAR），涵蓋了途徑中所有已知的關(guān)鍵酶。其中，編碼CHS和LAR基因的Unigene數(shù)量最多，均為12個(gè)，而編碼F3，5H的Unigene最少，僅有1個(gè)。如圖7所示，苯丙氨酸解氨酶（PAL）催化苯丙氨酸轉(zhuǎn)化為反式肉桂酸，經(jīng)過肉桂酸-4-羥化酶（C4H）、對(duì)香豆酸-輔酶A連接酶（4CL）的作用將苯丙烷代謝途徑（KEGG 數(shù)據(jù)庫通路 ID：ko00940）與黃酮生物合成途徑連接。香豆酰輔酶A在CHS的作用下形成查爾酮，隨后在CHI作用下形成柚皮素，再經(jīng)過F3H的催化形成二氫山奈酚。二氫山奈酚是黃酮合成途徑下游的一個(gè)節(jié)點(diǎn)物質(zhì)，能夠在不同酶的作用下生成多種黃酮類物質(zhì)。例如，二氫山奈酚在F3，5H的作用下生成二氫楊梅素，隨后在LAR的催化下生成兒茶酚。兒茶酚是含羞草中一種主要的黃酮類物質(zhì)。另外，二氫山奈酚在F3H的作用下生成二氫槲皮素，進(jìn)而經(jīng)過FLS的催化生成類黃酮（圖7）。本研究也鑒定出一些Unigene注釋在氮素代謝（ko00910）、萜類骨架的生物合成（ko00900）和苯丙素生物合成（ko00940）等次生代謝途徑中，其中注釋在苯丙素生物合成途徑的Unigene數(shù)量最多，為212條（表4）。

2.7 SSR標(biāo)記的開發(fā)

在含羞草轉(zhuǎn)錄組數(shù)據(jù)庫中共鑒定出7 382個(gè)SSR標(biāo)記。根據(jù)SSR的類型可分為單核苷酸（3 554）、二核苷酸（1 637）、三核苷酸（1 978）、四核苷酸（160）、五核苷酸（30）以及六核苷酸重復(fù)（23）。其中單核苷酸A/T（1 692/1 784）和三核苷酸GAA/TTC（172/125）為主要重復(fù)類型。從開發(fā)出的SSR標(biāo)記中隨機(jī)選取10對(duì)引物，并進(jìn)行PCR擴(kuò)增驗(yàn)證，結(jié)果如圖8所示，7對(duì)引物對(duì)能夠擴(kuò)增出清晰條帶，泳道4的引物擴(kuò)增出的條帶較暗，而泳道9、10的2對(duì)引物未能成功擴(kuò)增出條帶。

3 討論與結(jié)論

轉(zhuǎn)錄組測(cè)序具有分析效率高、覆蓋面廣，且不需要了解物種基因組信息便可獲得大量轉(zhuǎn)錄水平的數(shù)據(jù)，在植物代謝途徑解析[24]、SSR標(biāo)記位點(diǎn)挖掘[25]以及基因功能研究等[26]領(lǐng)域已有廣泛應(yīng)用。含羞草作為一種非模式植物，轉(zhuǎn)錄水平研究匱乏，參與黃酮生物合成的相關(guān)基因未見報(bào)道。本研究利用轉(zhuǎn)錄組測(cè)序技術(shù)對(duì)含羞草葉片進(jìn)行測(cè)序分析，共獲得20.08 GB數(shù)據(jù)，豐富了含羞草的轉(zhuǎn)錄組數(shù)據(jù)。原始數(shù)據(jù)經(jīng)過拼接后得到94 182個(gè)Unigene，GO富集分析將Unigene注釋于50個(gè)GO功能組。共有9 155個(gè)Unigene注釋在136條KEGG代謝通路，其中包括黃酮生物合成、苯丙素生物合成以及萜類骨架生物合成等重要的次生代謝途徑。

利用轉(zhuǎn)錄組數(shù)據(jù)開發(fā)SSR標(biāo)記具有廣闊的應(yīng)用前景。肖亮等[27]從葛根轉(zhuǎn)錄組中鑒定出25 452個(gè)SSR標(biāo)記，選取28對(duì)引物對(duì)不同葛根資源進(jìn)行了遺傳分析，表明不同葛根種質(zhì)具有很高的遺傳多樣性。梅利那等[28]基于轉(zhuǎn)錄組分析開發(fā)出大量馬尾松的SSR標(biāo)記。以多態(tài)性較高的24對(duì)引物對(duì)27份馬尾松進(jìn)行鑒定，發(fā)現(xiàn)供試樣品遺傳相似系數(shù)為0.32～0.92，UPGMA分析將樣品分為3大類。該研究結(jié)果為構(gòu)建馬尾松遺傳圖譜、繪制指紋圖譜提供了參考。本研究結(jié)果表明含羞草具有豐富的SSR位點(diǎn)，以單核苷酸（A/T）和三核苷酸（GAA/TTC）為主要類型。楊曉等[29]對(duì)秦艽，徐志軍等[30]對(duì)花生SSR的分析也得出相似的結(jié)果。開發(fā)出的SSR標(biāo)記能夠?yàn)榻窈笱芯亢卟葸z傳多樣性、種群結(jié)構(gòu)以及分子輔助育種提供參考，也為后續(xù)有針對(duì)性的開發(fā)藥用型、觀賞型含羞草資源提供了分子水平的依據(jù)。

黃酮生物合成途徑是植物體內(nèi)重要的次生代謝途徑之一。利用轉(zhuǎn)錄組測(cè)序能夠高效獲取黃酮合成相關(guān)的轉(zhuǎn)錄本，為解析黃酮合成途徑，分析基因功能提供參考。馬丹丹等[31]利用轉(zhuǎn)錄組測(cè)序技術(shù)對(duì)羅布麻進(jìn)行研究，鑒定出查爾酮合成酶、查爾酮異構(gòu)酶等參與黃酮合成的基因；姚運(yùn)法等[32]分析黃秋葵花和果莢中的類黃酮代謝途徑，結(jié)果表明黃秋葵花是通過F3H、DFR、ANS以及葡萄糖基轉(zhuǎn)移酶將柚皮素催化生成花青素苷，而黃秋葵果莢則是通過F3，5H將柚皮素催化成為二氫楊梅素，表明黃酮途徑在花和莢果中存在差異。本研究以轉(zhuǎn)錄組分析為基礎(chǔ)，鑒定出含羞草中黃酮合成途徑中的8個(gè)關(guān)鍵酶基因CHS、CHI、F3H、F3H、F3，5H、FLS、LAR和DFR的Unigene，初步構(gòu)建了黃酮生物合成途徑。此外，本研究也鑒定出參與異喹啉生物堿合成、萜類骨架生物合成以及苯丙素生物合成相關(guān)的Unigene。該結(jié)果與含羞草中含有L-含羞草堿、多酚類物質(zhì)的情況相符，為解析含羞草中其他藥效物質(zhì)的合成途徑提供了新的線索。

本研究對(duì)含羞草葉片進(jìn)行了轉(zhuǎn)錄組測(cè)序分析，獲得大量轉(zhuǎn)錄水平數(shù)據(jù)。GO富集分析表明Unigene主要注釋在“代謝過程”、“催化活性”以及“細(xì)胞”三個(gè)功能組。KEGG通路分析鑒定出49個(gè)Unigene注釋在黃酮生物合成途徑，編碼該途徑中的查爾酮合成酶、查爾酮異構(gòu)酶、黃烷酮-3-羥化酶、黃烷酮-3-羥化酶、黃烷酮-3，5-羥化酶、二氫黃酮醇4-還原酶、黃酮醇合成酶以及無色花色素還原酶基因。本研究結(jié)果豐富了含羞草的轉(zhuǎn)錄組數(shù)據(jù)，構(gòu)建了黃酮生物合成通路，為深入研究黃酮合成途徑關(guān)鍵酶基因的功能以及探索含羞草葉片開合的分子機(jī)制奠定了基礎(chǔ)。

參考文獻(xiàn) Reference：

PANG D D，LIU Y F，SUN Y N，et al.Full-length transcriptome and gene structure analysis of kucha（Camellia sinensis） ［J］.Acta Agriculturae Boreali-occidentalis Sinica，2021，30（4）：563-571.

［27］ 肖 亮，尚小紅，曹 升，等.基于轉(zhuǎn)錄組測(cè)序的葛根SSR標(biāo)記研究與利用［J］.西北植物學(xué)報(bào)，2019，39（1）：59-67.

XIAO L，SHANG X H，CAO SH，et al.Utilization of simple sequence repeat（SSR） markers development from a de novo transcriptome assembly in Pueraria thomsonii Benth［J］.Acta Botanica Boreali-Occidentalia Sinica，2019，39（1）：59-67.

［28］ 梅利那，范付華，崔博文，等.基于馬尾松轉(zhuǎn)錄組的SSR分子標(biāo)記開發(fā)及種質(zhì)鑒定［J］.農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，2017，25（6）：991-1002.

MEI L N，F(xiàn)AN F H，CUI B W，et al.Development of SSR molecular markers based on transcriptome sequences and germplasm identification in Masson pine（Pinus massoniana）［J］.Journal of Agricultural Biotechnology，2017，?25（6）：991-1002.

［29］ 楊 曉，馬子豪，馬 婕，等.粗莖秦艽種子萌發(fā)過程的轉(zhuǎn)錄組及關(guān)鍵因子分析［J］.中草藥，2021，52（1）：219-226.

YANG X，MA Z H，MA J，et al.Transcriptome and key factors analysis of Gentiana crasicaulis in seed germination process［J］.Chinese Traditional and Herbal Drugs，2021，52（1）：219-226.

［30］ 徐志軍，趙 勝，徐? 磊，等.基于 RNA-seq 數(shù)據(jù)的栽培種花生 SSR 位點(diǎn)鑒定和標(biāo)記開發(fā)［J］.中國(guó)農(nóng)業(yè)科學(xué)，2020，53（4）：695-706.

XU ZH J，ZHAO SH，XU L，et al.Discovery of microsatellite markers from RNA-seq data in cultivated peanut（Arachis hypogaea）［J］.Scientia Agricultura Sinica，2020，53（4）：695-706.

［31］ 馬丹丹，周佳熠，馬鵬舉，等.羅布麻轉(zhuǎn)錄組數(shù)據(jù)的注釋和黃酮合成相關(guān)基因的挖掘［J］.基因組學(xué)與應(yīng)用生物學(xué)，2019，38（10）：4580-4587.

MA D D，ZHOU J Y，MA P J，et al.Assembly on the data of Apocynum venetum transcriptome group and mining of genes involved in flavonoid synthesis［J］.Genomics and Applied Biology，2019，38（10）：4580-4587.

［32］ 姚運(yùn)法，張少平，練冬梅，等.黃秋葵花和果莢轉(zhuǎn)錄組測(cè)序及類黃酮代謝差異表達(dá)分析［J］.西北植物學(xué)報(bào)，2018， 38（11）：2000-2009.

YAO Y F，ZHANG SH P，LIAN D M，et al.Transcriptome sequencing and differential expression analysis of flavonoid metabolism in flowers and fruits of? okra［J］.Acta Botanica Boreali-Occidentalia Sinica，2018，38（11）：2000-2009.

Abstract To understand the flavonoid biosynthesis pathway of M.pudica，Illumina platform 2000TM platform was used to assemble the transcriptome of M.pudica leaf，a total of 94 182 unigene were obtained，with average length of 695 bp，the N50 value was 1 159 bp.Totally? of 30 243 unigene were annotated in 50 GO functional classification and dominated by “Metabolic process” “Catalytic activity” and “Cell”.A total of 7 382 SSR were explored from transcriptome database，among which mononucleotide and trinucleotide were the dominant types.Seven out of ten randomly selected primers could produce clear fragments.The results of KEGG analysis indicated that 49 unigene was identified in the flavonoid biosynthesis pathway，and chalcone synthase（CHS，12 unigene），chalcone isomerase（CHI，4），flavanone-3-hydroxylase（F3H，3），flavanone-3-hydroxylase（F3H，9），flavanone-3，5-hydroxylase（F3，5H，1），flavanone 4-reductase（DFR，5），flavonol synthase（FLS，3），and leucoanthocyanidin reductase（LAR，12） were encoded.

Key words Mimosa pudica Linn; Transcriptome sequencing; Metabolic pathway; Flavonoid biosynthesis