肖曉燕, 朱成磊, 楊克彬, 劉燕, 李紫陽, 郭棟, 高志民

低溫促進(jìn)毛竹葉片類黃酮合成及相關(guān)基因的表達(dá)模式分析

肖曉燕, 朱成磊, 楊克彬, 劉燕, 李紫陽, 郭棟, 高志民*

(國際竹藤中心竹藤資源基因科學(xué)與產(chǎn)業(yè)化研究所，國家林業(yè)和草原局/北京市共建竹藤科學(xué)與技術(shù)重點實驗室，北京 100102)

類黃酮在植物耐低溫脅迫方面發(fā)揮著重要作用，為揭示低溫對毛竹()葉片中類黃酮合成的影響，采用分光光度法測定了不同生長時期和低溫脅迫下毛竹幼苗葉片中的類黃酮含量，通過生物信息學(xué)方法對毛竹類黃酮早期生物合成關(guān)鍵酶基因進(jìn)行了鑒定，并用qPCR方法分析了其表達(dá)模式。結(jié)果表明，隨著葉片的生長，類黃酮含量呈現(xiàn)先升高后降低的趨勢，而低溫下，功能葉片中類黃酮含量則呈現(xiàn)上調(diào)趨勢，且在8 h時達(dá)極顯著水平。在毛竹基因組中鑒定了類黃酮早期生物合成3個酶基因家族共29個成員，包括20個查爾酮合酶基因(s)、8個查爾酮異構(gòu)酶基因(s)和1個黃酮-3-烴化酶基因()，這些基因的啟動子中均含有響應(yīng)低溫及其他非生物脅迫的調(diào)控元件。s傾向在根和葉中表達(dá), 而s為組成型表達(dá)。在不同生長時期的葉片中，僅表達(dá)與類黃酮的含量變化趨勢一致；而低溫脅迫下，3個s、2個s和在功能葉片中呈上調(diào)表達(dá)，與類黃酮含量變化趨勢一致。因此，毛竹可能通過提高類黃酮早期生物合成酶基因的表達(dá)量促進(jìn)類黃酮的合成來響應(yīng)低溫脅迫。

毛竹；類黃酮；分子特征；低溫脅迫；基因表達(dá)

溫度影響植物的生長發(fā)育，低溫脅迫是植物生長發(fā)育中最關(guān)鍵的非生物脅迫之一。低溫會直接影響植物的光合系統(tǒng)，使植物的生長發(fā)育受到抑制, 產(chǎn)生凍害和冷害[1]，嚴(yán)重時甚至死亡。葉片是植物最重要的光合器官，同時也是感知外界環(huán)境變化較為敏感的器官，能夠及時對低溫環(huán)境做出響應(yīng)。隨溫度降低，冬油菜()的葉面積、氣孔導(dǎo)度、光合速率、蒸騰速率等均逐漸下降[2]。在低溫下，菠菜()葉片的類黃酮總量增加，增強了植株的抗逆性并使糖分積累，葉片中的類黃酮含量與抗寒性的強弱有顯著相關(guān)性[3]。枳()實生苗的根系和葉片中的黃烷酮的含量隨著低溫脅迫時間的延長而增加，從而增強抗氧化脅迫的能力[4]。因此，類黃酮在植物應(yīng)對低溫脅迫方面發(fā)揮著重要作用。

類黃酮是植物中一類重要的次生代謝產(chǎn)物，其生物合成途徑前期需要查爾酮合酶(CHS)、查爾酮異構(gòu)酶(CHI)和黃烷酮-3-羥化酶(F3H)等3個酶家族，后期需要二氫黃酮醇-4-還原酶(DFR)、花青素合成酶(ANS)，無色花青素還原酶(ANR)和花青素還原酶(LAR)等4個家族參與，經(jīng)過一系列反應(yīng)生成不同種類的類黃酮[5]。研究表明，隨茉莉花()花瓣的顏色變淺，的表達(dá)量降低[6]；通過沉默非洲菊()中的，花瓣顏色由粉色變?yōu)榧儼咨玔7]；通過過表達(dá)牡丹()可以改變煙草()花瓣的著色[8]。除了參與花色的形成，類黃酮合成相關(guān)酶基因也參與響應(yīng)逆境脅迫。低溫脅迫下泡核桃()中的顯著上調(diào)表達(dá)[9]；芹菜()葉柄中的CHI和CHS酶活性在低溫下顯著升高[10]；紅花草莓(×)的表達(dá)會隨著溫度的升高而受到抑制[11]；低溫下擬南芥()中、和顯著上調(diào)，誘導(dǎo)了花青素的積累[12]。過表達(dá)番茄()的能增強煙草的耐冷性[13]。

竹子是重要的森林資源，具有生長快速的特點，其光合作用、初生增粗和筍的木質(zhì)化等已開展了廣泛研究[14–16]。毛竹()是我國栽培面積最大的竹種，已從葉片組織結(jié)構(gòu)、代謝組和轉(zhuǎn)錄組方面對其響應(yīng)低溫進(jìn)行了分析，找到了類黃酮代謝途徑的差異表達(dá)基因[17–18]。另外，研究表明毛竹中可能受不同激素的調(diào)控[19]。對七彩紅竹(‘Rainbow’)的研究表明,和在幼嫩紅稈中顯著表達(dá)[20–21]。因此，對竹子類黃酮代謝及其酶基因的研究尚不夠深入。本研究以毛竹葉片為材料，研究不同生長時期及低溫脅迫條件下類黃酮含量的變化趨勢，通過生物信息學(xué)的方法鑒定毛竹類黃酮早期生物合成的酶基因，明確其分子特征與系統(tǒng)進(jìn)化關(guān)系，并通過qPCR驗證關(guān)鍵酶基因的表達(dá)模式，以期為深入研究基因的功能提供參考依據(jù)。

1 材料和方法

1.1 材料

將毛竹()種子播種于草炭?蛭石=3:2的基質(zhì)中，置于約25 ℃的人工氣候室(光16 h/暗8 h，相對濕度約80%)中培養(yǎng)。選擇株高和長勢一致的3個月大小的毛竹實生苗，一部分用于采集不同生長時期的葉片，包括卷葉(S1)、半卷葉(S2)、展葉(S3)和老葉(S4)；另一部分隨機分為4組，每組4盆，置于4 ℃培養(yǎng)箱中模擬低溫脅迫處理。根據(jù)葉片類黃酮的含量變化，選擇含量最高時期的葉片作為低溫處理材料，分別在低溫處理0 (對照)、4和8 h后采集葉片，液氮速凍后置于-80 ℃冰箱中用于后續(xù)實驗。

1.2 類黃酮含量測定

類黃酮的基本結(jié)構(gòu)是A和C環(huán)、B和C環(huán)2個共軛體系，用蘆丁標(biāo)準(zhǔn)品繪制標(biāo)準(zhǔn)曲線。首先制備蘆丁標(biāo)準(zhǔn)品母液(1.0 mg/mL)，用60%乙醇稀釋為0、0.2、0.4、0.6、0.8、1.0 mg/mL共5個濃度梯度,在510 nm處測定吸光值，以蘆丁濃度為橫坐標(biāo), 吸光值為縱坐標(biāo)繪制標(biāo)準(zhǔn)曲線。采用亞硝酸鈉-硝酸鋁比色法提取葉片類黃酮并測定吸光度后，根據(jù)標(biāo)準(zhǔn)曲線計算類黃酮含量[22]。類黃酮含量(mg/g FW)= (A+0.0049)×÷(1.6277×)×D, 式中,為提取液體積(1.5 mL)；為樣品質(zhì)量(g)；為稀釋倍數(shù), 未稀釋即為1。

1.3 基因查找與生物信息學(xué)分析

利用毛竹最新基因組數(shù)據(jù)[23]，以水稻()和擬南芥的3個類黃酮生物合成酶家族(CHS、CHI和F3H)的基因序列為誘餌，利用TBtools v1.092軟件[24]內(nèi)置的BLASN和BLASTX進(jìn)行同源序列比對，閾值E<10–20?；贜CBI的CDD在線軟件進(jìn)行保守結(jié)構(gòu)域分析，保留具有完整類黃酮早期生物合成酶基因保守結(jié)構(gòu)域的序列。并對3個類黃酮早期生物合成酶基因家族成員命名。

利用ProtParam (http://web.Expasy.org/protparam/)對類黃酮早期生物合成酶的理化性質(zhì)進(jìn)行分析。運用TBtools軟件以及毛竹基因組注釋文件對毛竹類黃酮早期生物合成酶基因的結(jié)構(gòu)進(jìn)行分析；使用Plant-mPLoc server在線軟件(http://www.csbio.sjtu.edu. cn/bioinf/plant-multi/)進(jìn)行亞細(xì)胞定位的預(yù)測[25]。使用PlantCARE (http://bioinformatics.Psb.Ugent.be/web tools/plantcare/html/)對其啟動子序列(2 000 bp)中所含調(diào)控元件進(jìn)行分析預(yù)測[26]。

為探究毛竹類黃酮合成酶基因的潛在功能及不同物種間的進(jìn)化關(guān)系，從Phytozome V.12數(shù)據(jù)庫(https://phytozome.jgi.doe.gov/)獲得水稻、擬南芥、玉米()和高粱()等物種類黃酮早期生物合成酶的氨基酸序列，利用MEGA 7.0，基于p-distance模型采用鄰接法對毛竹和其他物種的序列進(jìn)行比對分析，生成系統(tǒng)進(jìn)化樹[27]，校驗參數(shù)為1 000次重復(fù)，氨基酸替換模型選用No. of differences[28]。

1.4 組織表達(dá)模式分析

從NCBI Short Read Archive下載毛竹組織(鞭、筍、根、葉片、芽和鞘)的26個轉(zhuǎn)錄組數(shù)據(jù)(登錄號：SRX2408703～SRX2408728)[23]，篩選其中類黃酮早期合成酶基因的FPKM (fragments per kilobase of exon model per million mapped fragments)，并用Log2(FPKM+1)的值作為基因的表達(dá)量，用TBtools軟件進(jìn)行表達(dá)譜熱圖繪制，可視化分析各基因的表達(dá)情況。

1.5 葉片中的定量表達(dá)分析

選擇關(guān)鍵的毛竹類黃酮早期生物合成酶基因,根據(jù)其基因序列，使用Primer Premier 5.0軟件設(shè)計特異定量引物(表1)。用植物RNA提取試劑盒(天漠, TR205-50, 中國)提取RNA，用反轉(zhuǎn)錄試劑盒反轉(zhuǎn)錄為cDNA (TaKaRa, RR037A, 日本)。qPCR試驗在qTOWER2.2系統(tǒng)上進(jìn)行，體系參照Roche Light Cycler?480 SYBR Green I Master kit試劑盒(Roche, 4707516001, 德國)，擴增程序：95 ℃預(yù)變性10 min，然后95 ℃變性10 s，60 ℃解鏈10 s，40個循環(huán)。選擇作為內(nèi)參基因[29]，相對表達(dá)量采用2–ΔΔCT法[30]計算。

2 結(jié)果和分析

2.1 毛竹葉片中類黃酮含量的變化

用蘆丁標(biāo)準(zhǔn)品建立標(biāo)準(zhǔn)曲線，回歸方程為= 1.6277-0.0049，相關(guān)系數(shù)為2=0.999 9，表明蘆丁在一定的范圍內(nèi)具有良好的線性關(guān)系(圖1: A)。根據(jù)蘆丁標(biāo)準(zhǔn)曲線，分別得到不同生長時期和低溫脅迫下毛竹葉片的類黃酮含量。在葉片的4個生長時期，類黃酮的含量存在顯著差異，隨葉片的生長呈先增加后降低的趨勢，且在S3時達(dá)到最高(圖1: B)。低溫處理下毛竹S3葉片的類黃酮含量隨脅迫時間的延長而增加，在8 h時顯著高于對照，差異達(dá)極顯著水平(圖1: C)，這表明毛竹葉片的生長時期和低溫均會影響類黃酮含量。

表1 實時定量PCR引物

圖1 毛竹葉片的類黃酮含量。A: 蘆丁標(biāo)準(zhǔn)曲線; B: 不同生長時期; C: 低溫處理; S1: 卷葉期; S2: 半卷葉期; S3: 展葉期; S4: 老葉期; **: P<0.01。

2.2 類黃酮早期生物合成酶基因鑒定

毛竹葉片類黃酮的形成受到其生物合成酶基因表達(dá)的影響，因此對毛竹中類黃酮早期生物合成酶基因進(jìn)行了鑒定。經(jīng)過在毛竹基因組數(shù)據(jù)庫中比對分析，共獲得29個編碼具有完整保守結(jié)構(gòu)域的類黃酮早期生物合成酶的基因，分別屬于CHS、CHI和F3H家族，分別命名為～、～和?；蚪Y(jié)構(gòu)分析表明，大部分s (14個)含有1個內(nèi)含子，的內(nèi)含子數(shù)量最多(3個)，、、、和均含有2個內(nèi)含子；s內(nèi)含子數(shù)量差異很大，有2～13個，其中最多(13個)，最少(2個)，有2個內(nèi)含子(圖2)。

啟動子序列分析結(jié)果表明，s、s和的啟動子序列中包含多種非生物脅迫響應(yīng)元件，如低溫的LTR、干旱的MBS和缺氧的ARE等(圖3: A)。s和s的啟動子中均有STRE (逆境響應(yīng)防御)，其中和的啟動子最多(9個)；多數(shù)s和s啟動子中含有低溫響應(yīng)防御元件(MYB-like sequence)、逆境傷害相關(guān)(WUN-motif)和真菌誘導(dǎo)(W box)等；4個s和的啟動子中都含有低溫響應(yīng)元件(LTR)；s和的啟動子中均含有逆境響應(yīng)防御元件(TC-rich repeats)。此外，在這些序列中還鑒定到其他調(diào)控元件，如MYB轉(zhuǎn)錄因子的結(jié)合元件，與激素相關(guān)的響應(yīng)元件等(圖3: B)。這說明類黃酮早期生物合成基因的轉(zhuǎn)錄可能會受到低溫及其他各種非生物脅迫的影響。

s編碼蛋白長度為368～461 aa，分子量為39.94～50.37 kDa，理論等電點為5.54～9.46，亞細(xì)胞定位預(yù)測在葉綠體和細(xì)胞質(zhì)中。s編碼蛋白長度為164～400 aa，分子量為18.75～44.92 kDa, 理論等電點為4.83～9.26，預(yù)測亞細(xì)胞定位均在葉綠體。編碼蛋白長度為363 aa，分子量為40.23 kDa，理論等電點為5.38，預(yù)測亞細(xì)胞定位在細(xì)胞質(zhì)(表2)?；蚪Y(jié)構(gòu)和蛋白理化性質(zhì)的差異可能是導(dǎo)致基因功能不同的主要原因。

2.3 PeCHSs和PeCHIs的系統(tǒng)進(jìn)化

根據(jù)水稻中類黃酮早期生物合成CHS的分類構(gòu)建，毛竹20個PeCHSs聚類在2個組，分別含有2個和18個成員，其中分支I中PeCHS14和PeCHS16分別與已知功能的OsCHS24和OsCHS8高度同源，推測它們可能具有合成查爾酮的功能(圖4: A)；分支II中PeCHSs與OsCHSs各成員存在不同程度的同源性，但功能尚不確定。通過構(gòu)建毛竹與其他物種CHI的進(jìn)化分類，PeCHIs在Group I、Group III和Group IV 3個組中(圖4: B)，PeCHI1和PeCHI5與ZmCHI聚類在Group I分支，推測它們具有合成柚皮素的功能。在F3H的進(jìn)化樹中，毛竹PeF3H1與水稻OsF3H同源性最高，推測它們可能具有相似的功能。

2.4 PeCHSs、PeCHIs和PeF3H1的組織表達(dá)特性

組織表達(dá)分析在一定程度上能夠驗證基因的功能，根據(jù)毛竹26個組織的轉(zhuǎn)錄組數(shù)據(jù)分析了類黃酮早期合成酶基因的轉(zhuǎn)錄水平。3個基因家族的各個成員在同一組織的不同生長時期以及不同組織中的表達(dá)模式存在差異(圖5)。s在毛竹的根、筍和葉的部分組織中能檢測到表達(dá)，更傾向在根和葉中表達(dá)，如、、和僅在葉中高表達(dá)，、、和僅在根中高表達(dá)；大部分s基因在筍中(除和外)的表達(dá)水平較低；而在26個組織中均沒有檢測到。s (除和外)在各個組織中均有表達(dá)，它們可能與毛竹整個生長發(fā)育過程密切相關(guān)。傾向在根和籜片中表達(dá)。

圖2 PeCHSs、PeCHIs和PeF3H1的基因結(jié)構(gòu)

圖3 PeCHSs、PeCHIs和PeF3H1啟動子中的順式作用元件(A)和應(yīng)答元件(B)

表2 PeCHSs、PeCHIs和PeF3H1的基本特征

圖4 基于不同種植物的CHS (A)和CHI (B)家族成員構(gòu)建的系統(tǒng)發(fā)育樹。Pe: 毛竹; Os: 水稻; Zm: 玉米; At: 擬南芥; Gm: 大豆; La: 羽扇豆; Lj: 百脈根; Ma: 紫苜蓿; Ah: 落花生; Ph: 矮牽牛; Dc: 康乃馨; Cs: 茶樹。

2.5 PeCHSs、PeCHIs和PeF3H1的定量分析

植物在不同生長時期和低溫脅迫下的類黃酮含量會發(fā)生顯著變化，以適應(yīng)外界環(huán)境。根據(jù)系統(tǒng)進(jìn)化、啟動子元件和轉(zhuǎn)錄組數(shù)據(jù)分析，篩選獲得6個代表基因(、、、、和)進(jìn)行qPCR分析。結(jié)果表明，6個基因隨著葉片的生長表達(dá)量變化差異較大，其中呈現(xiàn)先升高后降低的趨勢，與類黃酮含量的趨勢相同，在S3達(dá)到最高峰，表達(dá)量是S1的687.7倍；、和的表達(dá)量顯著下調(diào)；呈現(xiàn)先下調(diào)后上調(diào)的表達(dá)趨勢，在S4顯著升高，是S1的41.58倍；而的表達(dá)則呈現(xiàn)波動變化的趨勢，在S2和S3變化不顯著，在S4時，其表達(dá)量顯著上調(diào)，是S1的4.57倍(圖6: A)。經(jīng)過低溫脅迫后，這6個基因均呈現(xiàn)逐漸上升的表達(dá)趨勢，且均在8 h達(dá)到峰值, 與類黃酮含量表達(dá)趨勢一致(圖6: B)。因此，推測他們在毛竹葉片中的類黃酮合成和抵御低溫脅迫中發(fā)揮著關(guān)鍵作用。

圖5 PeCHSs、PeCHIs和PeF3H1在毛竹26個組織中的表達(dá)模式

3 結(jié)論和討論

類黃酮化合物廣泛存在于植物體內(nèi)，是一類重要的次生代謝產(chǎn)物，在植物遭受低溫、干旱等非生物脅迫時，類黃酮合成的相關(guān)基因被誘導(dǎo)表達(dá)，通過清除活性氧在植物的生長發(fā)育和逆境脅迫中發(fā)揮重要功能。毛竹廣泛分布于我國長江流域及南方各省區(qū)，具有重要的經(jīng)濟價值，但是溫度制約了毛竹生長和地理分布，尤其是冬季低溫和冰凍雨雪災(zāi)害的影響尤為嚴(yán)重。因此，提高毛竹抵御低溫脅迫的能力是其育種的重要目標(biāo)之一。本研究通過檢測毛竹葉片中類黃酮含量，包括不同生長時期和4 ℃低溫處理下的功能葉片，發(fā)現(xiàn)類黃酮含量發(fā)生了明顯變化，表明在毛竹生長發(fā)育和低溫條件下類黃酮可能發(fā)揮著重要作用。通過分析類黃酮早期生物合成途徑中涉及的酶基因家族成員發(fā)現(xiàn)，毛竹CHS家族成員與其它物種相似，不同分支的家族成員可能具有不同的功能。雖然PeCHS13和OsCHS27聚類在一起，但水稻中OsCHS27缺乏2個保守的氨基酸而不具有催化合成類黃酮的功能[31]，因此PeCHSs的具體功能需要進(jìn)一步研究。PeCHIs家族的成員多于水稻和擬南芥，這與該家族基因發(fā)生過復(fù)制和擴張有關(guān)[23]。Group I中PeCHI1和PeCHI5與玉米和水稻中具有合成類黃酮功能的CHI在同一組；Group III中的CHI屬于脂肪酸結(jié)合蛋白(FAPs),PeCHI2、PeCHI3、PeCHI4、PeCHI7和PeCHI8可能影響毛竹細(xì)胞中脂肪酸的生物合成及儲存[32]; Group IV作為類黃酮產(chǎn)生的增強劑，PeCHI6可能間接促進(jìn)類黃酮的生物合成[33]。F3H家族在毛竹中只有1個基因，這與F3H以單拷貝形式存在，在進(jìn)化上保守[34]一致。

圖6 PeCHSs、PeCHIs和PeF3H1在不同時期(A)和低溫處理(B)葉片中的相對表達(dá)量。*: P<0.05; **: P<0.01。

研究表明，類黃酮生物合成酶基因的表達(dá)具有發(fā)育時期特異性。在大豆中，傾向于在綠葉中表達(dá)，而傾向于在子葉中表達(dá)[35]。葡萄的在幼葉、幼根、果實和種子中均有表達(dá)[36]。銀杏()葉片幼葉時期，雖然分生組織活躍，類黃酮合成快，但由于葉片自身的合成能力有限，因此類黃酮的含量較低；在展葉時期，葉片自身的合成能力強，類黃酮含量較高[37]。在毛竹卷葉中的類黃酮含量最低，功能葉片中的最高，不同時期葉片中僅的表達(dá)與類黃酮含量一致，其它類黃酮早期生物合成基因的表達(dá)變化差異明顯但沒有規(guī)律，和在卷葉時期表達(dá)量最高，和在老葉時期表達(dá)量最高，這可能與葉片處于不同生長時期以及類黃酮生物合成途徑存在分支有關(guān)。

當(dāng)受到低溫脅迫時，植物體內(nèi)通過上調(diào)表達(dá)類黃酮合成酶基因促進(jìn)其積累來增強植物對低溫的耐受性。對低溫處理的桃樹()葉片進(jìn)行轉(zhuǎn)錄組分析表明，類黃酮合成相關(guān)基因顯著上調(diào)[38]; 低溫促進(jìn)煙草葉片和苦蕎麥()中類黃酮相關(guān)基因上調(diào)表達(dá)，從而提高類黃酮的含量[39–40]。本研究對低溫處理毛竹葉中類黃酮含量和基因表達(dá)的分析也獲得了類似的結(jié)果，類黃酮含量及其合成酶基因(3個s、2個s和)的表達(dá)隨著低溫處理時間的延長而逐漸增加，且在8 h時達(dá)到峰值。然而，除外, 其他類黃酮早期生物合成的關(guān)鍵酶基因在不同生長階段葉片中的表達(dá)與類黃酮含量變化并不一致，其原因可能與葉片生長發(fā)育的復(fù)雜調(diào)控有關(guān)，有待于進(jìn)一步研究。

本研究在毛竹中共鑒定到3個酶基因家族的29個成員，分別為20個查爾酮合酶基因(s)、8個查爾酮異構(gòu)酶基因(s)和1個黃酮-3-烴化酶基因()。在基因結(jié)構(gòu)、進(jìn)化分析、理化性質(zhì)、順式作用元件和轉(zhuǎn)錄組數(shù)據(jù)等生物學(xué)分析的基礎(chǔ)上，通過低溫下葉片中類黃酮含量的變化和不同基因的表達(dá)差異，初步確定了、、、、和可能是毛竹類黃酮早期生物合成的關(guān)鍵酶基因，為進(jìn)一步解析類黃酮早期合成關(guān)鍵酶基因的功能和抗逆基因資源挖掘提供了參考依據(jù)。

[1] DING Y H, WANG M, XIE J, et al. Response to low temperature stress in plants and advances in research methods [J]. Jiangsu Agric Sci, 2019, 47(14): 31–36. [丁紅映, 王明, 謝潔, 等. 植物低溫脅迫響應(yīng)及研究方法進(jìn)展 [J]. 江蘇農(nóng)業(yè)科學(xué), 2019, 47(14): 31–36. doi: 10.15889/j. issn.1002-1302.2019.14.007.]

[2] LI P, MA L, XU F, et al. Comparison of morphological and physio- logical responses of two different species of northern winter rapeseed before overwintering [J]. Agric Res Arid Areas, 2022, 40(5): 42–51. [李鵬, 馬驪, 徐芳, 等. 越冬前北方不同類型強冬性冬油菜形態(tài)及生理響應(yīng)研究 [J]. 干旱地區(qū)農(nóng)業(yè)研究, 2022, 40(5): 42–51. doi: 10. 7606/j.issn.1000-7601.2022.05.05.]

[3] WATANABE M, AYUGASE J. Effect of low temperature on flavor- noids, oxygen radical absorbance capacity values and major compo- nents of winter sweet spinach (L.) [J]. J Sci Food Agric, 2015, 95(10): 2095–2104. doi: 10.1002/jsfa.6925.

[4] YAN X. The effects of environmental stresses on the compont and gene expression of flovoinds in[D]. Wuhan: Huazhong Agricultural University, 2011. [晏校. 逆境脅迫對枳實生苗類黃酮組分含量及關(guān)鍵酶基因表達(dá)量的影響 [D].武漢: 華中農(nóng)業(yè)大學(xué), 2011.]

[5] XU W J, GRAIN D, BOBET S, et al. Complexity and robustness of the flavonoid transcriptional regulatory network revealed by comprehensiveanalyses of MYB-bHLH-WDR complexes and their targets inseed [J]. New Phytol, 2014, 202(1): 132–144. doi: 10.1111/nph.12620.

[6] LI M, CAO Y T, YE S R, et al. Isolation ofgene fromflowers and its regulation in anthocyanin biosynthesis [J]. Molecules, 2016, 22(1): 44. doi: 10.3390/molecules22010044.

[7] HELARIUTTA Y, ELOMAA P, KOTILAINEN M, et al. Chalcone synthase-like genes active during corolla development are differentially expressed and encode enzymes with different catalytic properties in(Asteraceae) [J]. Plant Mol Biol, 1995, 28(1): 47–60. doi: 10.1007/BF00042037.

[8] ZHOU L, WANG Y, REN L, et al. Overexpression of, a homologue of the chalcone isomerase gene from tree peony (), reduces the intensity of flower pigmentation in transgenic tobacco [J]. Plant Cell Tiss Organ Cult, 2014, 116(3): 285–295. doi: 10.1007/s11240-013-0403-2.

[9] WANG Y, XIAO L Y, MA T, et al. Cloning and functional analysis of chalcone synthase gene in cold-induced[J]. Mol Plant Breed, 2018, 16(2): 386–391. [王毅, 肖良俊, 馬婷, 等. 低溫誘導(dǎo)泡核桃中查爾酮合成酶基因克隆及功能分析 [J]. 分子植物育種, 2018, 16(2): 386–391. doi: 10.13271/j.mpb.016.000386.]

[10] ZHANG M. Effects of temperature, drought and UV-B irradiation time on celery growth and flavonoids accumulation [D]. Nanjing: Nanjing Agricultural University, 2017. [張萌. 溫度、干旱和UV-B照射時間對芹菜生長及黃酮類物質(zhì)積累的影響 [D]. 南京: 南京農(nóng)業(yè)大學(xué), 2017.]

[11] LIANG J C. Response of flower color to temperature and analysis of flower color related gene expression in red-flowered strawberry [D]. Shenyang: Shenyang Agricultural University, 2018. [梁吉昌. 紅花草莓花色對溫度的響應(yīng)及花色相關(guān)基因表達(dá)分析 [D]. 沈陽: 沈陽農(nóng)業(yè)大學(xué), 2018.

[12] ZHANG Y Q, ZHENG S, LIU Z J, et al. Both HY5 and HYH are necessary regulators for low temperature-induced anthocyanin accumu- lation inseedlings [J]. J Plant Physiol, 2011, 168(4): 367– 374. doi: 10.1016/j.jplph.2010.07.025.

[13] JIA H F, JIU S, ZHANG C, et al. Abscisic acid and sucrose regulate tomato and strawberry fruit ripening through the abscisic acid-stress- ripening transcription factor [J]. Plant Biotechnol J, 2016, 14(10): 2045–2065. doi: 10.1111/pbi.12563.

[14] ZHAO H S, LOU Y F, SUN H Y, et al. Transcriptome and comparative gene expression analysis ofin response to high light [J]. BMC Plant Biol, 2016, 16: 34. doi: 10.1186/s12870-016-07 20-9.

[15] WEI Q, JIAO C, GUO L, et al. Exploring key cellular processes and candidate genes regulating the primary thickening growth of moso underground shoots [J]. New Phytol, 2017, 214(1): 81–96. doi: 10. 1111/nph.14284.

[16] YANG K B, LI L C, LOU Y F, et al. A regulatory network driving shoot lignification in rapidly growing bamboo [J]. Plant Physiol, 2021, 187(2): 900–916. doi: 10.1093/plphys/kiab289.

[17] WANG H Y, GUO L, ZHA R F, et al. Histological, metabolomic, and transcriptomic analyses reveal mechanisms of cold acclimation of moso bamboo () leaf [J]. Tree Physiol, 2022, 20: tpac064. doi: 10.1093/treephys/tpac064.

[18] HUANG Z, ZHU P L, ZHONG X J, et al. Transcriptome analysis of Moso bamboo () reveals candidate genes involved in response to dehydration and cold stresses [J]. Front Plant Sci, 2022, 13: 960302. doi: 10.3389/fpls.2022.960302.

[19] ZHANG N, BAO H, CUI B L, et al. Genomic analysis of chalcone synthase gene family of[J]. Mol Plant Breed, 2022, 20(3): 817–825. [章妮, 暴涵, 崔博亮, 等. 毛竹查爾酮合成酶基因家族全基因組分析 [J]. 分子植物育種, 2022, 20(3): 817–825. doi: 10.13271/j.mpb.020.000817.]

[20] ZHOU X, WANG C C, BI W, et al. Cloning and sequence analysis of chalcone synthase genefromcv. Rainbow [J]. J Bamboo Res, 2014, 33(3): 11–16. [周旭, 王晨晨, 畢瑋, 等. 七彩紅竹查爾酮合成酶基因的克隆與分析 [J]. 竹子研究匯刊, 2014, 33(3): 11–16. doi: 10.3969/j.issn.1000-6567.2014.03.002.]

[21] WANG C C, BI W, WANG J, et al. Cloning and expression analysis offromcv. Rainbow [J]. J W China For Sci, 2014, 43(2): 85–90. [王晨晨, 畢瑋, 王娟, 等. 七彩紅竹查爾酮異構(gòu)酶基因的克隆與表達(dá)分析 [J]. 西部林業(yè)科學(xué), 2014, 43(2): 85–90. doi: 10.16473/j.cnki.xblykx1972.2014.02.018.]

[22] XIA G H, LI X H, WANG X Y. Effect of different processing on the antioxidant activity of flavonoids from[J]. Food Sci Technol, 2018, 43(3): 181–185. [夏光輝, 李新華, 王曉雅. 不同加工處理對玉竹黃酮抗氧化活性的影響研究 [J]. 食品科技, 2018, 43(3): 181–185. doi: 10.13684/j.cnki.spkj.2018.03.034.]

[23] PENG Z H, LU Y, LI L B, et al. The draft genome of the fast-growing non-timber forest species moso bamboo () [J]. Nat Genet, 2013, 45(4): 456–461. doi: 10.1038/ng.2569.

[24] CHEN C J, CHEN H, ZHANG Y, et al. TBtools: An integrative toolkit developed for interactive analyses of big biological data [J]. Mol Plant, 2020, 13(8): 1194–1202. doi: 10.1016/j.molp.2020.06.009.

[25] CHOU K C, SHEN H B. Cell-PLoc 2.0: An improved package of web- servers for predicting subcellular localization of proteins in various organisms [J]. Nat Sci, 2010, 2(10): 1090–1103. doi: 10.4236/ns.2010. 210136.

[26] LESCOT M, DéHAIS P, THIJS G, et al. PlantCARE, a database of plant-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences [J]. Nucl Acids Res, 2002, 30(1): 325– 327. doi: 10.1093/nar/30.1.325.

[27] KUMAR S, STECHER G, TAMURA K. MEGA7: Molecular evolu- tionary genetics analysis version 7.0 for bigger datasets [J]. Mol Biol Evol, 2016, 33(7): 1870–1874. doi: 10.1093/molbev/msw054.

[28] JONES D T, TAYLOR W R, THORNTON J M. The rapid generation of mutation data matrices from protein sequences [J]. Comput Appl Biosci, 1992, 8(3): 275–282. doi: 10.1093/bioinformatics/8.3.275.

[29] FAN C J, MA J M, GUO Q R, et al. Selection of reference genes for quantitative real-time PCR in bamboo () [J]. PLoS One, 2013, 8(2): e56573. doi: 10.1371/journal.pone.0056573.

[30] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2?ΔΔCTMethod [J]. Methods, 2001, 25(4): 402–408. doi: 10.1006/meth.2001.1262.

[31] PARK H L, YOO Y, BHOO S H, et al. Two chalcone synthase isozymes participate redundantly in UV-induced sakuranetin synthesis in rice [J]. Int J Mol Sci, 2020, 21(11): 3777. doi: 10.3390/ijms2111 3777.

[32] NGAKI M N, LOUIE G V, PHILIPPE R N, et al. Evolution of the chalcone-isomerase fold from fatty-acid binding to stereospecific catalysis [J]. Nature, 2012, 485(7399): 530–533. doi: 10.1038/nature 11009.

[33] MORITA Y, TAKAGI K, FUKUCHI-MIZUTANI M, et al. A chalcone isomerase-like protein enhances flavonoid production and flower pigmentation [J]. Plant J, 2014, 78(2): 294–304. doi: 10.1111/tpj.12469.

[34] PELLETIER M K, SHIRLEY B W. Analysis of flavanone 3-hydro- xylase inseedlings: Coordinate regulation with chalcone synthase and chalcone isomerase [J]. Plant Physiol, 1996, 111(1): 339– 345. doi: 10.1104/pp.111.1.339.

[35] ANGURAJ VADIVEL A K, KRYSIAK K, TIAN G, et al. Genome- wide identification and localization of chalcone synthase family in soybean ((L.) Merr) [J]. BMC Plant Biol, 2018, 18(1): 325. doi: 10.1186/s12870-018-1569-x.

[36] ZHOU J, YAO Q H, PENG R H, et al. Cloning and expression analysis ofofgrape by semi-quantity RT-PCR [J]. Acta Bot Boreali- Occid Sin, 2009, 29(9): 1723–1729. [周軍, 姚泉洪, 彭日荷, 等. 巨峰葡萄查爾酮異構(gòu)酶基因克隆及表達(dá)分析 [J]. 西北植物學(xué)報, 2009, 29(9): 1723–1729. doi: 10.3321/j.issn:1000-4025.2009.09.001.]

[37] VOGT T, GUL P G. Accumulation of flavonoids during leaf develop- ment in[J]. Phytochemistry, 1994, 36(3): 591–597. doi: 10.1016/S0031-9422(00)89780-0.

[38] ZHOU P, LIN Z K, GUO R, et al. Effects of low temperature treatment on gene expression and flavonoids biosynthesis metabolism in peach () leaves [J]. J Agric Biotechnol, 2021, 29(7): 1283– 1294. [周平, 林志楷, 郭瑞, 等. 低溫處理對桃樹葉片基因表達(dá)及類黃酮合成代謝的影響 [J]. 農(nóng)業(yè)生物技術(shù)學(xué)報, 2021, 29(7): 1283– 1294. doi: 10.3969/j.issn.1674-7968.2021.07.005.]

[39] ZHANG G, PANG T, WANG S S, et al. Regulation of different growth temperatures on flavonoids content and their metabolic pathways in tobacco leaves [J]. Genom Appl Biol, 2016, 35(12): 3527–3537. [張剛, 逄濤, 王莎莎, 等. 不同生長溫度對煙草葉片黃酮類化合物含量及其代謝途徑的調(diào)控 [J]. 基因組學(xué)與應(yīng)用生物學(xué), 2016, 35(12): 3527–3537. doi: 10.13417/j.gab.035.003527.]

[40] SONG Y J, FENG J C, LIU D M, et al. Different phenylalanine pathway responses to cold stress based on metabolomics and transcript- tomics in tartary buckwheat landraces [J]. J Agric Food Chem, 2022, 70(2): 687–698. doi: 10.1021/acs.jafc.1c06915.

Low Temperature Promotes Flavonoid Synthesis and Expression Pattern Analysis of the Related Genes in Leaves of

XIAO Xiaoyan, ZHU Chenglei, YANG Kebin, LIU Yan, LI Ziyang, GUO Dong, GAO Zhimin*

(Institute of Gene Science and Industrialization for Bamboo and Rattan Resources, International Center for Bamboo and Rattan, Key Laboratory of National Forestry and Grassland Administration/Beijing for Bamboo & Rattan Science and Technology, Beijing 100102, China)

Flavonoids play an important role in plant resistance to low temperature stress. To reveal the effects of low temperature on flavonoid synthesis in leaves of, the content of flavonoid in leaves of bamboo seedlings at different growth stages and under low temperature stress was determined by spectrophoto- metric method, the key enzyme genes of early biosynthesis of flavonoid in bamboo were identified by bioinfor- matics methods and their expression patterns were analyzed using qPCR. The results showed that the flavonoid content increased at first and then decreased, while that in the functional leaves increased under low temperature and reached significant level at 8 hours. There were a total of 29 members, belonging to three gene families, involved in early biosynthesis of flavonoid in., including 20 chalcone synthase genes (s), eight chalcone isomerase genes (s) and one flavanone-3-hydroxylase gene (), the promoters of these genes all containedregulatory elements in response to low temperature and other abiotic stress responses. Thes tended to be expressed in roots and leaves, whiles expressed constitutively. Only the expression trend ofwas consistent with the content of flavonoid in.leaves at different growth stages, while six key enzyme genes (3s, 2s and) were upregulated continuously in functional leaves under low temperature, which was consistent with the change trend of flavonoid content. Therefore,might respond to low temperature stress by increasing the expression of genes involved in early biosynthesis of flavonoid.

; Flavonoid; Molecular characteristics; Low temperature stress; Gene expression

10.11926/jtsb.4717

2022-08-16

2022-10-02

國家重點研發(fā)計劃項目(2021YFD2200502)資助

This work was supported by the National Key Research and Development Program of China (Grant No. 2021YFD2200502).

肖曉燕(1996年生)，女，碩士研究生，研究方向為毛竹生長發(fā)育的分子基礎(chǔ)研究。E-mail: xiaoxiaoyan@icbr.ac.cn

* 通訊作者 Corresponding author. E-mail: gaozhimin@icbr.ac.cn