參與連翹苯丙烷合成途徑的WRKY轉錄因子的篩選與鑒定

張佳琪，劉曉冉，陳佳茜，郭廣洋，譚新杰，趙杏利，胥華偉，劉 浩，侯典云*

參與連翹苯丙烷合成途徑的WRKY轉錄因子的篩選與鑒定

張佳琪1, 2，劉曉冉1, 2，陳佳茜1, 2，郭廣洋1, 2，譚新杰1, 2，趙杏利1, 2，胥華偉1, 2，劉 浩1, 2，侯典云1, 2*

1. 河南科技大學農學院，河南 洛陽 471023 2. 洛陽市道地藥材繁育與創新利用工程技術研究中心，河南 洛陽 471023

篩選出參與連翹苯丙烷合成途徑的WRKY轉錄因子并進行生物信息學分析。通過對從連翹轉錄組數據中篩選出的62條連翹WRKY基因序列進行鑒定；通過采用蛋白理化性質和基序分析、蛋白結構分析、系統進化分析、功能注釋、蛋白互作分析、外源激素處理后實時熒光定量分析等方式進行相關分析。最終篩選出52條具有WRKY結構域的連翹WRKY轉錄因子。WRKY轉錄因子編碼蛋白的氨基酸數目在105～728，相對分子質量在12 277.95～80 321.99。蛋白基序分析顯示其均有WRKY結構域，與擬南芥WRKY轉錄因子構建系統進化樹將連翹52個WRKY轉錄因子進一步分為3大類；篩選出5個可能參與連翹苯丙烷合成途徑的連翹WRKY轉錄因子，并推測其可能以被MAPK級聯反應所調控的方式介導連翹中苯丙烷合成途徑。通過分析轉錄組測序結果，從52個連翹WRKY轉錄因子中篩選出5個可能參與連翹苯丙烷合成途徑的轉錄因子，其表現出組織特異性表達，且對不同濃度茉莉酸甲脂表現出不同的響應，為進一步研究其作用機制提供研究方向。

連翹；次生代謝產物；WRKY；轉錄因子；苯丙烷合成途徑

連翹(Thunb.) Vahl為木犀科連翹屬落葉灌木，多見于中國河南、河北、山東、陜西、山西等地區[1]，其干燥果實入藥，具有清熱解毒、消腫散結、疏散風熱等功效，用于風熱感冒、溫病初起、高熱煩渴、神昏發斑、熱淋澀痛等癥治療[2]。據報道[3]，連翹作為傳統中藥，在新型冠狀病毒肺炎的預防和治療中，也具有較好的療效。研究表明，連翹中有多種次生代謝產物發揮藥理作用，如類黃酮及萜類、酚酸類等[4]，其中，通過苯丙烷途徑產生的木脂素類（如連翹苷）和苯乙醇苷類（如連翹脂苷）是連翹的主要活性成分[3]。苯丙烷途徑含有豐富的苯丙烷類化合物和天然化合物。超過70%的醫藥產品是以植物代謝物為基礎的，其中許多來自苯丙烷途徑，包括木質素、黃酮類、香豆素、醌和木脂素，它們是許多結構聚合物的基本成分，也是調節生物相互作用的花色素、香味化合物或信號分子等[5]。在許多中藥植物中，如唇形科植物丹參、豆科植物黃芪等，其主要生物活性成分為苯丙烷類化合物，已被證明具有抗氧化劑、自由基清除劑、抗炎和抗癌化合物的作用[6]。

茉莉酸甲酯（methyl jasmonate，MeJA）被證實可作為一種外源誘導因子對藥用植物體內藥理成分的積累產生影響，如雷公藤經MeJA誘導之后，顯著提高了愈傷中雷公藤紅素的合成[7]，另外在連翹中也證明了茉莉酸甲脂噴施后連翹葉中連翹苷、連翹脂苷等都有不同程度影響[8]。

WRKY轉錄因子作為植物中最大的轉錄因子家族之一，僅存在于高等植物中，自1994年Ishiguro等[9]在甘薯中第一次發現WRKY轉錄因子以來，越來越多的WRKY轉錄因子在多種植物中被鑒定和驗證。據報道，WRKY轉錄因子家族在多種植物中參與苯丙烷代謝途徑，在次生代謝物質合成中起到重要的調控作用，如生姜WRKY轉錄因子Zo WRKY1[10]、葡萄WRKY轉錄因子Vv WRKY2[11]和Vv WRKY26[12]、苜蓿WRKY轉錄因子Mt WRKY 100630和Mt WRKY 108715[13]均參與調控苯丙烷代謝途徑，但在連翹中未見報道。

本研究在課題組對連翹噴施200 μmol/L茉莉酸甲脂12 h和未做處理的連翹果實、葉片進行比較轉錄組測序的基礎上，通過采用蛋白理化性質和基序分析、蛋白結構分析、系統進化分析等方式最終篩選出52條連翹WRKY轉錄因子，此外，通過進化樹分析、功能注釋、蛋白互作分析、外源激素處理后實時熒光定量分析等方式，篩選出5個可能參與連翹苯丙烷合成途徑的連翹WRKY轉錄因子，并推測其可能以被絲裂原活化蛋白激酶（mitogen- activated protein kinase，MAPK）級聯反應所調控的方式介導連翹中次生代謝物質的合成，而其表現出組織特異性表達，且對不同濃度茉莉酸甲脂表現出不同的響應，為進一步研究其作用機制提供研究方向。

1 材料與儀器

1.1 材料

研究所用的樣品采自河南洛陽市宜陽縣連翹種植基地，由侯典云教授鑒定為連翹(Thunb.) Vahl。

1.2 儀器

5415R型低溫高速離心機（Eppendorf公司）、JY04S-3D型凝膠成像分析系統（君意電泳公司）、實時熒光定量PCR儀（Lightcycler96, Roche公司，瑞士）等。

2 方法

2.1 實驗設計

7月下旬，選擇生長環境一致，生長狀況良好的連翹植株掛簽標明噴施濃度和日期，于7: 00時對葉片進行噴施不同濃度茉莉酸甲脂，噴施濃度為50、200、400 μmol/L，每濃度噴施10株連翹。用于轉錄組測序的材料為200 μmol/L噴施12 h后的連翹葉片和果實及其對照。擬南芥WRKY轉錄因子家族序列來源于TAIR（https://www.arabidopsis.org/ index.jsp）。

2.2 連翹WRKY轉錄本的篩選與功能注釋

從連翹轉錄組序列中篩選獲得62條WRKY轉錄因子序列，通過結構域在線分析軟件SMART（http://smart.embl.de/）對保守結構域進行分析，使用ORF Finder在線分析網站（https://www.ncbi. nlm.nih.gov/orffinder/）分析預測轉錄因子的開放閱讀框（open reading frame，ORF），獲得連翹WRKY轉錄因子家族全長氨基酸序列。將從轉錄組數據中篩選到的連翹WRKY轉錄因子比對到基因本體（gene ontology，GO）、京都基因與基因組百科全書（Kyoto encyclopedia of genes and genomes，KEGG）數據庫進行注釋與分析。

2.3 連翹WRKY家族蛋白理化性質和蛋白基序分析鑒定

利用在線工具Expasy ProtParam tool（https:// web.expasy.org/protparam/）分析連翹WRKY蛋白的基本理化性質。利用SOPMA和signaIP 5.0對蛋白質二級結構特征和蛋白信號肽進行預測。通過CELLO在線預測蛋白亞細胞定位。利用The MEME Suite 5.4.1（https://meme-suite.org/meme/tools/ meme）在線分析連翹WRKY家族蛋白的基序，基序數量設置為10，其余參數為默認值。

2.4 系統進化樹的構建

從TAIR擬南芥數據庫中下載完整的擬南芥WRKY轉錄因子氨基酸序列。利用MEGA7.0軟件構建連翹與擬南芥WRKY轉錄因子蛋白序列進化樹，利用Clustal W進行多序列比對分析，采用鄰接法（neighbor-joining，NJ）構建系統進化樹，設置重復次數為1000，其余參數選擇默認值。

2.5 連翹苯丙烷合成途徑相關WRKY轉錄因子功能分析

通過查閱文獻獲得與苯丙烷合成途徑相關的WRKY轉錄因子基因序列，與連翹WRKY轉錄因子構建系統進化樹，推測連翹WRKY轉錄因子的功能。利用String在線蛋白互作數據庫對篩選出的連翹WRKY轉錄因子進行蛋白互作網絡分析，物種參數選擇擬南芥。

2.6 連翹苯丙烷合成途徑相關WRKY轉錄因子的驗證

使用ChamQ Universal SYBR qPCR Master Mix（南京諾維贊生物科技股份有限公司）對以上預測出的苯丙烷合成途徑相關WRKY轉錄因子進行熒光定量PCR驗證，分析連翹葉片、果實中WRKY轉錄因子的表達量隨著茉莉酸甲脂噴施不同濃度12 h后的變化趨勢。熒光定量反應程序：預變性95 ℃、30 s，變性95 ℃、10 s，退火60 ℃、30 s，40個循環。反應體系20 μL：2×ChamQ Universal SYBR qPCR Master Mix 10 μL，上、下游引物各0.5 μL，cDNA 1 μL，用ddH2O補至20 μL。所用引物見表1。

3 結果與分析

3.1 連翹WRKY轉錄因子篩選

從本實驗室連翹轉錄組序列中篩選出62條WRKY轉錄因子序列，通過ORF預測全長氨基酸及SMART分析序列保守結構域結果見圖1，發現ORF長度為321～2304 bp，10條不具有完整的WRKY結構域的序列，52條具有完整ORF和完整WRKY結構域的序列。

3.2 連翹WRKY家族蛋白理化性質和結構分析

連翹WRKY家族蛋白理化性質見表2，連翹WRKY蛋白的氨基酸數量為105～728，相對分子質量介于12 277.95～80 321.99。51.9%的蛋白等電點小于7.00，且親水性系數均小于零，說明該家族蛋白大多數為酸性、親水蛋白。同時，82.7%的蛋白不穩定系數大于40，該家族蛋白大多數為不穩定蛋白。蛋白質二級結構預測結果顯示，連翹WRKY家族蛋白均有α-螺旋、β-轉角、無規則卷曲和延伸鏈，其中，無規則卷曲占64.1%，是最主要的二級結構形式。信號肽和亞細胞定位預測結果表明52個連翹WRKY家族蛋白均不具有信號肽，且除TRINITY_DN1857_c0_g1（定位在葉綠體）、TRINITY_DN9153_c0_g2（定位在葉綠體）、TRINITY_DN37721_c0_g1（定位在細胞外間隙）外，其它的均定位在細胞核中。

表1 qRT-PCR引物

Table 1 Primers for qRT-PCR

基因ID引物序列（5’-3’） TRINITY_DN14046_c0_g1CCTGCCCGTGACACTAACAT CACCTGCTTCTGCCCATACT TRINITY_DN25488_c0_g1ATCCCGACAGAGTTCCACCG GGTTCATCGCCTTCTGCCTC TRINITY_DN28667_c0_g1TTTTCAGGGCTCCAGATGTC GCTAAGCGTGTCGGTAAAGG TRINITY_DN3231_c0_g1TGCTAGTCTACTGCCCAAGGAT CATAGATGAAGAAGTGGGTGGC TRINITY_DN6456_c0_g2CAAAGGTTCTCCTCACCCAAG GGATCATCCAAGCATTTCTCA UKNCAGACCAGCTTTGAGGAGTATC GGCCAGAAACCAGTAGTCAATA

3.3 連翹WRKY轉錄因子保守結構域鑒定

利用在線軟件MEME對連翹WRKY家族蛋白序列進行分析（圖2），發現5個蛋白序列具有5個motif，52個連翹WRKY蛋白均具有WRKY結構域（motif 1、motif 3），進一步證實其均為WRKY家族蛋白。

3.4 連翹WRKY轉錄因子的功能注釋與分析

3.4.1 GO數據庫注釋結果統計與分析 GO數據庫注釋與分析結果表明，連翹WRKY轉錄因子中有52個注釋到細胞組分（cellular component，CC），2個注釋到分子功能（molecular function，MF），1個注釋到生物過程（biological process，BP），見圖3。

3.4.2 KEGG數據庫注釋結果統計與KEGG代謝通路分析 由圖4可以看出有6條序列被注釋到環境信息處理（environmental information processing， EIP）中的植物MAPK信號通路（MAPK signaling pathway-plant）中，11條序列被注釋到OS中的植物-病原相互作用（plant-pathogen interaction）中。

圖1 連翹WRKY蛋白家族結構域預測

表2 連翹WRKY蛋白理化性質分析

Table 2 Analysis of physical and chemical properties of WRKY protein from F. suspensa

UnigenORF長度/bp氨基酸數目等電點不穩定系數脂肪系數親水平均系數二級結構 α-螺旋/%β-轉角/%無規則卷曲/%延伸鏈/% TRINITY_DN32638_c0_g110473486.6547.1159.43?0.96928.451.4460.34 9.77 TRINITY_DN15730_c0_g1 7682574.9953.3158.02?0.60512.063.1173.9310.89 TRINITY_DN92_c0_g110473487.0749.1863.05?0.480 5.752.8775.8615.52 TRINITY_DN130_c0_g121877285.5352.2252.90?0.79410.992.8875.6910.44 TRINITY_DN31960_c0_g1 8672885.2070.5744.69?1.09121.182.7865.9710.07 TRINITY_DN17865_c0_g1 9033006.0656.6452.00?0.96420.001.3371.00 7.67 TRINITY_DN5195_c0_g216625536.1163.0865.37?0.97412.842.7170.7113.74 TRINITY_DN26114_c0_g1 6692228.9956.0252.30?0.75631.982.2550.4515.32 TRINITY_DN7278_c0_g2 5491828.6838.7747.64?0.86226.374.9552.2016.48 TRINITY_DN3781_c0_g113444476.8460.0261.25?0.61510.512.0176.7310.74 TRINITY_DN19770_c0_g1 5191729.6735.0465.06?0.81925.006.4051.7416.86 TRINITY_DN36133_c1_g210623476.1454.4969.71?0.77438.901.4448.9910.66 TRINITY_DN32493_c0_g1 6992328.2046.5451.21?0.99311.217.3361.6419.83 TRINITY_DN15948_c0_g1 5701898.2350.8055.71?1.02820.115.8258.2015.87 TRINITY_DN2185_c0_g117103045.9075.8844.93?1.043 7.572.6381.58 8.22 TRINITY_DN9835_c0_g1 4291844.7567.8352.34?1.09229.351.6359.78 9.24 TRINITY_DN34934_c0_g1 3631217.7831.5378.18?0.88944.635.7940.50 9.09 TRINITY_DN11182_c0_g2 8072688.6350.4063.69?0.74731.723.3654.4810.45 TRINITY_DN17876_c0_g211494536.5956.5061.06?0.93515.452.8771.0810.60 TRINITY_DN1015_c1_g2 7772655.4460.7052.72?0.70910.193.0276.6010.19 TRINITY_DN10738_c0_g1 8942977.6869.7044.31?1.00524.921.6863.64 9.76 TRINITY_DN130_c0_g221217065.8254.5658.16?0.71811.612.8374.65 2.83 TRINITY_DN28667_c0_g114974986.0950.8565.30?0.58930.322.6157.23 9.84 TRINITY_DN13376_c0_g111731469.3047.0272.19?0.22426.034.1147.2622.60 TRINITY_DN3231_c0_g123042098.9248.6043.35?1.362 8.134.7875.6011.48 TRINITY_DN17876_c0_g111493976.00.60.4265.04?0.88217.382.0269.5211.08 TRINITY_DN25488_c0_g112394128.9134.9162.48?0.84513.112.6772.3311.89 TRINITY_DN6456_c0_g210383459.1956.7864.70?0.77522.614.0664.35 8.99 TRINITY_DN3962_c0_g1 5761719.2554.6768.95?0.61733.336.4350.29 9.94 TRINITY_DN3381_c0_g1 4171859.8969.1162.70?0.90815.683.2474.05 7.03 TRINITY_DN17815_c0_g110503494.9747.7563.41?0.69720.064.0167.91 8.02 TRINITY_DN33631_c0_g111281649.1466.9363.54?0.55215.242.4475.00 7.32 TRINITY_DN1051_c0_g1 5821934.6952.9267.15?0.52518.652.0769.95 9.33 TRINITY_DN53314_c0_g1 5371799.5444.3655.98?0.93030.177.8253.07 8.94 TRINITY_DN18337_c0_g1 9063015.7358.1256.38?0.81027.912.3362.46 7.31 TRINITY_DN7399_c0_g1 6902337.6761.5542.66?0.95920.603.4362.6613.30 TRINITY_DN11067_c0_g3 3691229.3833.5141.39?1.13322.957.3850.8218.85

續表2

UnigenORF長度/bp氨基酸數目等電點不穩定系數脂肪系數親水平均系數二級結構 α-螺旋/%β-轉角/%無規則卷曲/%延伸鏈/% TRINITY_DN56780_c0_g11725330 9.2542.0065.18?0.52221.821.2166.0610.91 TRINITY_DN5090_c0_g11443325 5.2643.5156.74?0.71616.922.7768.9211.38 TRINITY_DN18709_c0_g3 981326 6.3750.2951.75?0.85822.091.8469.63 6.44 TRINITY_DN37721_c0_g1 52516110.0369.8746.65?1.027 7.452.4875.7814.29 TRINITY_DN37541_c0_g12187728 5.8150.1859.05?0.75011.542.7575.6910.03 TRINITY_DN48409_c0_g2 351105 9.4077.4449.24?1.109 8.571.9066.6722.86 TRINITY_DN9153_c0_g2 579156 9.6861.1286.92?0.42432.693.8550.6412.82 TRINITY_DN13886_c0_g1 573190 4.8444.9650.84?0.73124.212.6367.89 5.26 TRINITY_DN1857_c0_g11299432 5.7959.2248.10?0.923 7.643.2479.86 9.26 TRINITY_DN41966_c0_g11002333 6.6875.6652.43?0.84919.222.1069.07 9.61 TRINITY_DN9289_c0_g1 489162 9.1536.3078.77?0.28333.956.1740.1219.75 TRINITY_DN14046_c0_g11587528 6.5861.2755.80?0.89413.263.2270.8312.69 TRINITY_DN20171_c0_g1 573190 5.8435.5586.00?0.16324.212.6350.5322.63 TRINITY_DN49536_c0_g1 573190 9.2140.0068.26?0.58930.006.3252.1111.58 TRINITY_DN10845_c0_g1 591196 5.0534.5465.10?0.61425.004.0861.22 9.69

圖2 連翹WRKY蛋白家族保守基序和基因結構分析

Fig.2 Analysis of conserved motifs and gene structure of WRKY protein family in F. suspensa

3.5 連翹WRKY蛋白進化樹分析

利用MEGA 7.0構建連翹和擬南芥WRKY蛋白系統進化樹，結果如圖5所示，參考擬南芥WRKY分類依據，將連翹WRKY家族分為3大類，第I類包含12個連翹WRKY家族蛋白，其中有9個蛋白具有2個WRKY結構域；第II類包含33個連翹WRKY家族蛋白，根據其鋅指結構的不同進一步分為IIb（5個）、IIc（17個）、IId（1個）、IIe（10個），連翹WRKY蛋白不具有IIa成員；第Ⅲ類包含7個連翹WRKY家族蛋白，且皆與擬南芥對應類蛋白在同一分支。

圖3 連翹WRKY轉錄因子GO二級分類圖

圖4 連翹WRKY的KEGG代謝通路富集圖

3.6 連翹苯丙烷合成途徑相關WRKY轉錄因子功能分析

基于NJ法構建系統進化樹結果表明，連翹TRINITY_DN14046_c0_g1、TRINITY_DN3231_c0_g1、TRINITY_DN25488_c0_g1、TRINITY_DN28667_c0_ g1、TRINITY_DN6456_c0_g2分別與Vv WRKY2、Mt WRKY 108715、Vv WRKY26、Mt WRKY 100630和Zo WRKY1聚為一支（圖6），預測這些連翹WRKY蛋白也可能參與苯丙烷合成途徑，繼而調控連翹次生代謝物質的合成。

圖5 連翹和擬南芥WRKY蛋白系統進化分析

圖6 連翹與其他參與苯丙烷途徑WRKY轉錄因子家族系統進化樹

3.7 連翹WRKY蛋白互作分析

利用STRING在線數據庫檢索構建了篩選出的5個連翹WRKY轉錄因子互作蛋白網絡（圖7）。與其互作的蛋白有MPK3、MPK4、SIB1、SIB2、MKS1等，其中MPK3在丹參中被證實Sm MPK3（可能與At MPK3同源）通過MAPK級聯調控Sm AREB1、Sm MYB36、Sm MYB39、Sm MYB111、Sm PAP1參與酚酸合成途徑[14]。因此推測篩選出的與MAPK3互作的3個連翹WRKY轉錄因子可能參與苯丙烷途徑的連翹WRKY轉錄因子與MPK3相互作用來調控連翹中次生代謝物質的形成。

3.8 連翹苯丙烷合成途徑相關WRKY轉錄因子的驗證分析

為驗證5個可能參與連翹苯丙烷合成途徑相關的連翹WRKY轉錄因子，分析了不同濃度茉莉酸甲脂處理12 h后，5個連翹WRKY轉錄因子對不同濃度茉莉酸甲脂處理12 h后的相對表達量變化。結果見圖8，連翹的5個WRKY轉錄因子在果實和葉片中均有表達，除TRINITY_DN3231_c0_g1外，其余4個轉錄因子在葉片中的表達量均高于果實。

在50、200、400 μmol/L茉莉酸甲脂處理12 h后，5個轉錄因子均表現出對不同濃度茉莉酸甲脂的響應，結果見圖9，除TRINITY_DN14046_c0_g1在葉片中表現出隨著激素濃度的升高，表達量也升高外，其余均為低濃度促進，高濃度抑制，圖10在果實中也呈現出不同程度的低濃度促進、高濃度抑制現象，且有研究表明，在外源激素處理后，參與植物次生代謝的基因均對其有不同程度的響應。因此，進一步推測5個轉錄因子可能參與連翹苯丙烷合成途徑。

圖7 連翹5個WRKY轉錄因子蛋白互作分析

不同字母表示差異顯著，P＜0.05，下同

圖9 5個連翹WRKY轉錄因子在不同濃度MeJA處理12 h后葉片中相對表達量

圖10 5個連翹WRKY轉錄因子在不同濃度MeJA處理12 h后果實中相對表達量

4 討論

植物可以采用多種方法應對生物和非生物脅迫，其中許多防御機制是基于次生代謝物的產生[15]。其中一些次生代謝物質還具有營養或藥理活性（抗生素、抗癌、抗糖尿病等）[16]。據估計，目前使用的治療化合物中至少有30%來自藥用植物的次生代謝物[17]。越來越多的證據表明，WRKY轉鐵蛋白調節植物的抗逆性和次生代謝的同時也可調節植物對非生物脅迫的相應，從而影響植物的生長發育[18]。且有研究表明，通過調控植物轉錄因子，可以改變次生代謝物的含量[19]。

MAPK級聯普遍存在于真核細胞中。這些級聯反應涉及4種不同類型的激酶：MAPKs、MAPKKs、MAPKKKs和MAPKKKs[20]。MAPK級聯對于植物的生長、發育以及防御生物脅迫和非生物脅迫都很重要[21]。此外，MAPK級聯很可能參與次生代謝，包括駱駝蓬蛋白和吲哚硫代葡萄糖苷、尼古丁、花青素、植保素[22]。同時，MAPK級聯調控SA、JA、生長素（in，AUX）和乙烯[23]等植物激素的生物合成和信號轉導途徑。Yang等[6]研究表明，連翹噴施不同濃度茉莉酸甲脂后，可對連翹中次生代謝物質產生不同程度的影響，即總黃酮含量隨著MeJA濃度的升高而增長；木脂素類化合物連翹苷含量表現出低濃度抑制，高濃度增加；苯乙醇苷類化合物連翹脂苷A積累整體上呈現抑制作用。因此推測，MAPK級聯可能會對外源MeJA處理后連翹的次生代謝物合成途徑產生影響。

連翹作為我國傳統中藥材，以果實入藥，產生的化合物能夠有效的治療一些疾病。本研究基于該條件通過對連翹轉錄組篩選出52個WRKY轉錄因子，進行蛋白結構域分析，蛋白一、二級結構預測等一系列生物信息學分析，篩選出5個可能參與連翹苯丙烷合成途徑的連翹WRKY轉錄因子，并通過STRING分析得到了5個轉錄因子可能被MAPK級聯反應所調控，而具體蛋白互作情況仍不清楚，需后續實驗驗證其作用方式。

利益沖突 所有作者均聲明不存在利益沖突

[1] 吳孟華, 石上梅, 曹暉. 青翹與老翹優劣的本草探析 [J]. 中國中藥雜志, 2019, 44(24): 5508-5512.

[2] 中國藥典[S]. 一部. 2015: 26.

[3] 齊麗娜, 陳炫好, 金華, 等. 中藥連翹化學成分及藥理活性研究進展 [J]. 天津中醫藥大學學報, 2021, 40(2): 168-175.

[4] 段文娟, 耿巖玲, 祝賀, 等. 中藥連翹化學成分和分析方法的研究進展 [J]. 山東科學, 2010, 23(2): 33-37.

[5] VOGT, T. Phenylpropanoid biosynthesis. (Special Issue: Plant metabolism.) [J]., 2010, 3(1):2-20.

[6] Yang L, Yang C Q, Li C Y,. Recent advances in biosynthesis of bioactive compounds in traditional Chinese medicinal plants [J]., 2016, 61(1): 3-17.

[7] Mosongo I M. 茉莉酸甲脂誘導雷公藤愈傷組織中雷公藤紅素合成的影響研究[D]. 北京: 中國科學院大學,2019.

[8] 陶娟, 宋小鋒, 周圓圓. MeJA對連翹葉片光合作用及次生代謝的影響 [J]. 河南農業科學, 2021, 50(6): 54-60.

[9] Ishiguro S, Nakamura K. Characterization of a cDNA encoding a novel DNA-binding protein, SPF1, that recognizes SP8 sequences in the 5' upstream regions of genes coding for sporamin and beta-amylase from sweet potato [J]., 1994, 244(6): 563-571.

[10] 唐建民, 齊力旺, 賀小婷, 等. 生姜ZoWRKY1基因的克隆及其與6-姜酚生物合成的相關性分析 [J]. 湖南農業科學, 2020(10): 1-5.

[11] Guillaumie S, Mzid R, Méchin V,. The grapevine transcription factor WRKY2 influences the lignin pathway and xylem development in tobacco [J]., 2010, 72(1/2): 215-234.

[12] Amato A, Cavallini E, Zenoni S,. A grapevine TTG2-like WRKY transcription factor is involved in regulating vacuolar transport and flavonoid biosynthesis [J]., 2016, 7: 1979.

[13] Naoumkina M A, He X Z, Dixon R A. Elicitor-induced transcription factors for metabolic reprogramming of secondary metabolism in[J]., 2008, 8: 132.

[14] Xie Y F, Ding M L, Zhang B,. Genome-wide characterization and expression profiling of MAPK cascade genes inreveals the function of SmMAPK3 and SmMAPK1 in secondary metabolism [J]., 2020, 21(1): 630.

[15] Hidalgo W, Chandran J N, Menezes R C,. Phenylphenalenones protect banana plants from infection byand are deactivated by metabolic conversion [J]., 2016, 39(3): 492-513.

[16] Salam A M, Quave C L. Opportunities for plant natural products in infection control [J]., 2018, 45: 189-194.

[17] Thiericke R. Drug discovery from Nature: Automated high-quality sample preparation [J]., 2000, 22(5): 149-157.

[18] Zhang Y, Butelli E, Alseekh S,. Multi-level engineering facilitates the production of phenylpropanoid compounds in tomato [J]., 2015, 6: 8635.

[19] Yousfi F E, Makhloufi E, Marande W,. Comparative analysis ofgenes potentially involved in salt stress responses in. ssp.[J]., 2016, 7: 2034.

[20] Zhang M X, Chiang Y H, Toru?o T Y,. The MAP4 kinase SIK1ensures robust extracellular ROS burst and antibacterial immunity in plants [J]., 2018, 24(3): 379-391.

[21] Ding H D, He J, Wu Y,. The tomato mitogen-activated protein kinase SlMPK1 is as a negative regulator of the high-temperature stress response [J]., 2018, 177(2): 633-651.

[22] Genot B, Lang J, Berriri S,. Constitutively activeMAP kinase 3 triggers defense responses involving salicylic acid and SUMM2 resistance protein [J]., 2017, 174(2): 1238-1249.

[23] Mine A , Berens M L , Nobori T ,. Pathogen exploitation of an abscisic acid- and jasmonate-inducible MAPK phosphatase and its interception by Arabidopsisimmunity [J]., 2017, 22: 236.

Screening and identification of transcription factors WRKY involved in phenylpropane synthesis pathway in

ZHANG Jia-qi1, 2, LIU Xiao-ran1, 2, CHEN Jia-xi1, 2, GUO Guang-yang1, 2, TAN Xin-jie1, 2, ZHAO Xing-li1, 2, XU Hua-wei1, 2, LIU Hao1, 2, HOU Dian-yun1, 2

1. College of Agricultural, Henan University of Science and Technology, Luoyang 471023, China 2. Luoyang Engineering Research Center of Breeding and Utilization of Dao-di Herbs, Luoyang 471023, China

To screen out the WRKY transcription factors involved in the phenylpropane synthesis pathway ofand conduct bioinformatics analysis.A total of 62 WRKY gene sequences screened fromtranscriptome data were identified. Then the correlation analysis was carried out by means of protein physical and chemical properties and motif analysis, protein structure analysis, phylogenetic analysis, functional annotation, protein interaction analysis, real-time fluorescence quantitative analysis after exogenous hormone treatment.A total of 52 WRKY transcription factors with WRKY domain were screened. The number of amino acids of the protein encoded by WRKY transcription factor ranged from 105 to 728, and its molecular weight ranged from 12 277.95 to 80 321.99. The analysis of protein motif showed that all of them had WRKY domains, and the phylogenetic tree was constructed withWRKY transcription factors to further divide 52 WRKY transcription factors ofinto three categories; FiveWRKY transcription factors that may participate in the phenylpropane synthesis pathway ofwere screened, and it was speculated that they may mediate the phenylpropane synthesis pathway ofin a way regulated by MAPK cascade reaction.By analyzing the transcriptome sequencing results, this study screened five transcription factors that may participate in the phenylpropane synthesis pathway ofWRKY from 52 transcription factors, which showed tissue-specific expression and different responses to different concentrations of methyl jasmonate, providing a research direction for the further study of its mechanism of action.

(Thunb.) Vahl; secondary metabolites; WRKY; transcription factor; phenylpropane synthesis pathway

R286.12

A

0253 - 2670(2023)18 - 6055 - 10

10.7501/j.issn.0253-2670.2023.18.023

2023-02-03

國家自然科學基金項目（U1404829）；中央本級重大增減支項目“名貴中藥資源可持續利用能力建設項目”（2060302）；河南省自然科學基金項目（202102110156）；河南省中藥材產業科技特派員服務團、河南省中藥產業技術體系建設專項資金資助

張佳琪（1998—），碩士研究生，主要從事藥用植物分子生物學研究。E-mail: 799535911@qq.com

侯典云（1975—），教授，博士生導師，主要從事藥用植物資源評價與利用。Tel: (0379)6428230 E-mail: dianyun518@163.com

[責任編輯 時圣明]