-
药用植物中的次生代谢物是中药药效物质的主要来源,已知的植物次生代谢物生物合成途径有乙酸-丙二酸途径、异戊二烯途径、莽草酸途径等[1],探究植物生物合成的调控因素不仅能提升药材的品质,也为中药有效成分体外合成的工业化提供可能。
植物昼夜节律钟是植物体内应对光照、温度等外界因素随昼夜节律性改变而进化出的一套适应机制[2],对植物生长发育具有不可或缺的作用。大量研究表明,如黄酮类化合物合成的相关结构基因表达,也具有明显的昼夜节律性特点,受昼夜节律钟调控[3-4]。昼夜节律钟的核心部分中央振荡器是MYB蛋白LHY(late elongated hypocotyl)、CCA1(circadian clock associated 1)和伪应答调控蛋白家族(PRRs,pseudo-response regulators)组成,对维持植物昼夜节律的稳定至关重要[5-6]。
PRRs基因都带有2个保守的结构域,氨基端的响应接受结构域(receiver-like domain,RLD),其结构上与磷酸接受域相似,羧基端带有的CCT(Constans/Constans-like/TOC1)结构域,这2个结构域被一个保守程度不高的“可变域”所分隔[7-9]。目前研究发现在CO、CO-like以及TOC1基因中也有带此类结构域,对于植物开花进程有着重要作用[10-11]。
研究表明PRRs家族基因具有增加植物抗逆性,影响植物生物量的积累[12-14],以及调控花发育及衰老等作用[15]。目前对模式植物中PRRs基因的研究较多,如在拟南芥以及水稻中的PRRs基因证明具有调控开花周期的作用[17],但药用植物中PRRs基因的研究则罕见报道。
中药红花(Carthami flos)是菊科植物红花(Carthamus tinctorius L.)的干燥花,有活血化瘀的功效。研究表明,红花主要药效物质为黄酮类化合物,如羟基红花黄色素A(HSYA)、红花素、槲皮素、山奈酚、野黄芩苷[18]等,目前已有对红花黄酮类化合物生物合成的关键基因查尔酮合酶、查尔酮异构酶、糖基转移酶等多种研究[19-20],但调控红花中黄酮类化合物生物合成途径的基因未完全明确。验证昼夜节律钟调控红花中黄酮类化合物的生物合成对提升红花品质意义重大。
本研究依据前期红花转录组数据库基因注释筛选昼夜节律相关基因,并与红花黄酮类化合物的积累量进行相关分析,得到具有调控红花黄酮类化合物生物合成功能的昼夜节律基因,通过qPCR、液质联用(UHPLC-MS)等方法以期阐释PRRs基因的特征与功能,为进一步研究昼夜节律钟调控红花黄酮类化合物的生物合成积累资料。
-
植物材料:云南巍山红花品系(编号ZHH0119),种植于海军军医大学药学院温室,温室条件:恒温25 ℃,昼夜节律为16 h光照/8 h黑暗。仪器与试剂:荧光qPCR仪:ABI7500;Phanta Max Master Mix高保真酶,Trans Top green qPCR super mix,Trans one-step cDNA synthesis super mix逆转录试剂盒(北京全式金公司);Meitler Toledo电子天平(十万分之一量程);高效液相色谱仪:Agilent1290 Infinty LC system; 质谱仪:Agilent 6538 Accurate Mass。
-
取用新鲜红花的花冠,根、茎、叶约100 mg研磨成粉。依据Plant Zol说明书提取总RNA,并用紫外分光光度计测定总RNA浓度与质量,样品的A260/A280在1.9~2.1之内可认为符合后续实验要求。将其作为模板逆转录合成cDNA第一链,−20 ℃保存备用。
-
基于红花转录组数据库,结合基因注释筛选出昼夜节律相关基因,通过红花花冠表达谱获取表达量,与红花代谢组数据库中不同花期黄酮类化合物的积累量进行Pearson分析,获取与黄酮类化合物积累量相关系数r≥0.7的昼夜节律基因进行生物信息学分析,设计引物(表1)进行克隆。PCR产物进行凝胶电泳,回收含有目的条带的凝胶,连接载体后转化大肠杆菌感受态细胞,LBA平板培养,挑取阳性克隆送至上海生工生物工程有限公司测序。
表 1 全长克隆引物
引物名称 引物序列 PRR1-F ACCTCAAGGGCCACTGGTTC PRR1-R GTAACAAAAGACTTTCTGAA -
在NCBI网站用BLAST在线分析红花PRR1全长序列以及编码蛋白的氨基酸序列进行比对以及同源性分析;在ExPASy在线工具SOMPA得出所编码蛋白的二级结构特征;通过Swiss-Model同源建模预测蛋白三级结构;在ClustalX 2.1软件对其编码蛋白的氨基酸序列与同源蛋白的氨基酸序列进行多重比对分析;使用相邻节点算法(Neighbor-Joining)构建系统进化树,自展分析法(Bootstrap analysis)进行1 000次重复验证进化树可靠性。
-
以红花根、茎、叶、花4个部位;花冠开花前3 d(Ⅰ期)、开花1 d(Ⅱ期)、开花3 d(Ⅲ期)、开花5 d(Ⅳ期)4个时期[12]以及1 d中8个时间点(6:00、9:00、12:00、15:00、18:00、21:00、0:00、3:00)的花冠第一链cDNA为模板,设计RT-qPCR引物,以Ct60s(KJ634810)作为内参基因,进行qPCR实验,每个样品设3个复孔。
-
取红花Ⅲ期单日内4个时间点(9:00、12:00、15:00、21:00)的花冠烘干至恒重,精密称取5.00 mg置1.5 ml离心管,再加入精密量取的1 ml 60%甲醇,室温放置12 h后超声处理40 min,12 000 r/min,10 min离心取上清液,UHPLC-MS在正、负离子模式下进行检测(表2),参比离子:正离子模式为121.050 9,922.009 8;负离子模式为119.036 3,966.000 7,数据采集与分析使用Agilent MassHunter Analysis4.0软件。以柚皮素、芹菜素、槲皮素、HSYA、山奈酚、红花素、山奈酚-3-O-葡萄糖苷以及野黄芩素为对照品。
表 2 UHPLC-MS分析方法
色谱条件 质谱条件 色谱柱:Agilent ZORBAX C18(3.5 μm,2.1 mm ×100 mm) 毛细管电压:3.5 kV 流动相:A:0.1%甲酸-水,B:乙腈 锥孔电压:60 V 流速:0.4 ml/min 离子源温度:350 ℃ 柱温:40 ℃ 碎裂电压:120 V 进样量:0.4 ml 采集范围:50~1 100 m/z 梯度洗脱:0 min,5%B; 4 min,20%B;
6 min, 21% B; 9 min, 26% B; 11 min, 40% B; 15 min, 80% B; 17 min, 95% B -
使用STRING数据库对CtPRR1基因进行蛋白互作网络分析,可信区间设置为0.4。利用Python下的networkX将互作基因进行可视化。
-
将红花cDNA均一化处理得到红花酵母双杂cDNA文库,与载体PGADT7进行同源克隆,产物转化大肠杆菌感受态细胞,涂布LBA平板,37 ℃培养过夜,挑取单克隆菌落进行菌液PCR验证
-
根据PRR1 ORF设计同源引物引入BamHI、XhoI酶切位点。PCR产物经电泳后回收,连接酶切载体pGBKT7转入大肠杆菌,测序验证后完成pGBKT7-PRR1 bait载体的构建。
将pGBKT7-PRR1bait质粒通过PEG/LiAc法转化酵母感受态细胞,涂布SD平板,28 ℃培养,挑取单克隆进行菌液PCR鉴定。将阳性单克隆酵母菌液划线涂布含有X-α-gal和Aba的平板。
-
将酵母感受态细胞、载体DNA、红花酵母双杂cDNA文库质粒按共转法混合处理后均匀涂布于SD/-His/-Trp/-Ura三缺平板培养,挑取单克隆转移至含有X-α-gal的SD/-His/-Trp/-Ura培养基,扩增阳性质粒用于测序以及酵母双杂验证。
-
实验结果经SPSS19.0软件处理,计量资料统一表示为(
${{\bar x}} \pm {{s}}$ ),组间比较采用ANOVA分析,以P<0.05为差异具有统计学意义。 -
基于红花花冠转录组数据库的基因注释,检索出PRR1、PRR2、ELF3、FT、PHYB、GI、ZTL7个昼夜节律基因,将其在花冠不同时期的表达量与红花代谢组数据库中柚皮素、HYSA、苯丙氨酸、山奈酚-3-O-葡萄糖苷、芹菜素、槲皮苷、野黄芩素、槲皮素-3-葡萄糖苷7个化合物含量进行PEARSON相关性分析,如图1所示,PRR1与红花中主要黄酮类成分的积累量具有良好的相关性(r≥0.7)。
-
在红花全长转录组数据库中得到的PRR1基因全长序列3 201 bp,ORF FINDER结果显示开放阅读框1 549~2 814 bp,编码421个氨基酸,命名为CtPRR1(GenBank登录号:MW492035),将全长序列进行BLAST,系统发育进化树分析如图2A,CtPRR1氨基酸序列与水稻(Oryza sativa L.)OsPRR73氨基酸序列(A2XFB7.2)同源性最高。Prot-param分析PRR1基因所编码的蛋白质分子式C1900H3039N611O653S15,分子量为45 300,理论pI=8.52,对PRR蛋白质三维结构预测如图2B,Prot Scale分析表明预测PRR1蛋白为亲水性蛋白,无信号肽属非分泌蛋白;蛋白跨膜性分析预测PRR不含有跨膜区域,为非跨膜蛋白。
-
qPCR结果表明(图3A)在红花不同部位中,CtPRR1基因在花中的表达量最高,且与根、茎、叶都有显著性差异(P<0.05),对红花花冠Ⅰ期、Ⅱ期、Ⅲ期、Ⅳ期的CtPRR1基因的表达情况进行分析,CtPRR1在Ⅰ期花冠中表达含量最低;Ⅱ期表达水平略有上升,与Ⅰ期无显著性差异;Ⅲ期花冠中表达量明显上升,相比Ⅰ,Ⅱ期具有显著性差异,Ⅳ期花冠中的CtPRR1基因表达水平略微下降(图3B),表明CtPRR1基因在Ⅲ期花冠中转录水平最高,CtPRR1在不同时间表达量的变化与红花黄酮类化合物积累量的变化相符[22]。对红花8个时间点的花冠CtPRR1表达量分析发现CtPRR1表达量在日间6:00至18:00逐渐升高,18:00达峰值,晚间21:00至2:00逐渐下降,3:00最低,表达水平在白天与夜晚两个连续的时间周期内呈现昼夜节律性(图3C)。
-
UHPLC-MS检测黄酮类化合物在单日内不同时间的含量变化(图4),除柚皮素外的7个化合物在白天含量逐渐降低,而在晚间含量又升至较高的水平,柚皮素则相反在白天升高,晚间下降,但8个化合物含量变化也呈现昼夜节律性,与同时间点CtPRR1的表达量进行PEARSON分析得到红花素(r=−0.9485)、山奈酚(r=−0.9423)、野黄芩苷(r=−0.9504)、HSYA(r=−0.8372)、山奈酚-3-O-葡萄糖苷(r=−0.8792)、柚皮素(r=0.7415)、芹菜素(r=−0.6652)、槲皮素(r=−0.4876),目前研究认为柚皮素在红花黄酮类化合物生物合成途径的上游[20],而CtPRR1与柚皮素呈正相关,与其他化合物呈负相关,说明CtPRR1对整体黄酮类化合物生物合成具有负调控作用而导致柚皮素的积累量增多,进一步说明CtPRR1参与调控红花黄酮类化合物的合成。
-
STRING数据库预测红花转录组数据库中CtPRR1的互作蛋白,获得互作蛋白TOC1、PIF3、COP1、ZTL、LHY、ELF4、ELF3、CCA1、GI、PRR5、GRP7、TIC、FKF1、RVE114个,互作关系网络如图5。与PRRs基因共同参与调控植物生长发育,调控植物次生代谢途径。由此获得了较为完善的红花昼夜节律核心元件系统,为进一步探索昼夜节律调控红花次生代谢的分子机制提供了依据。
通过酵母双杂实验发现,转入pGBKT7-PRR bait质粒的酵母细胞在28 ℃生长4 d后情况正常。带有pGBKT7-PRR bait载体的酵母细胞与构建的红花cDNA文库质粒共转后,28 ℃培养2 d后可以观察到蓝色菌斑出现,4 d后出现直径约2 mm的蓝色单克隆菌斑(图6),菌液PCR鉴定为阳性克隆。结果如表3所示,有2个热休克蛋白,3个AP2转录因子,热休克蛋白具有增强植物抗逆性的功能,AP2转录因子广泛参与植物生长发育,调控体内次生代谢。
表 3 酵母双杂筛选互作蛋白
序号 蛋白家族名称 预测功能 1 热休克蛋白 广泛参与植物响应外界环境胁迫,生物合成以及生长发育相关的转录调控 2 热休克蛋白 3 AP2转录因子 4 AP2转录因子 5 AP2转录因子
Characterization and function study of circadian rhythm gene CtPRR1 in Carthamus tinctorius L.
-
摘要:
目的 探究昼夜节律基因对红花黄酮类物质生物合成的影响及机制。 方法 基于红花花冠转录组及代谢组数据库筛选潜在调控红花黄酮类化合物生物合成的昼夜节律基因;用qPCR测定红花各部位以及花冠单日不同时间点昼夜节律基因的表达量,液质联用测定黄酮类化合物的积累量,并分析二者的相关性;酵母双杂交实验验证昼夜节律基因的互作蛋白。 结果 获得7个昼夜节律基因PRR1、PRR2、ELF3、FT、PHYB、GI、ZTL,其中PRR1基因与黄酮类化合物积累量呈正相关(r≥0.7)。PRR1全长3 201 bp,编码421个氨基酸,与水稻OsPRR73基因高度同源,将其命名为CtPRR1(GenBank登录号:MW492035)。CtPRR1主要在花中表达,表达量在日间逐渐升高,晚间逐渐下降;黄酮类化合物芹菜素、槲皮素、HSYA、山奈酚、Carthamin、山奈酚-3-O-葡萄糖苷以及野黄芩素的含量为白天逐渐降低,晚间逐渐升高,二者都有昼夜节律性且呈负相关(r≥−0.7)。酵母双杂实验得到2个热休克蛋白、3个AP2转录因子。 结论 CtPRR1对红花黄酮类成分的昼夜节律性积累起负调节作用;CtPRR1可能受这些互作蛋白的影响调控红花黄酮类成分的昼夜节律性积累。 Abstract:Objective To explore the effect of circadian rhythm genes on flavonoids biosynthesis in safflower and its molecular mechanism. Methods Based on the transcriptome and metabolomic database of safflower corolla, we screened the circadian rhythm genes that correlate with biosynthesis of flavonoids in safflower. qPCR was used to quantify the expressions of circadian rhythm genes in different flowering stages at different time points in a single day. LC-MS was performed to determine the accumulation of flavonoids. The correlation between them was analyzed as well. Yeast Two-Hybrid experiment was used to verify the interactive proteins of these genes. Results Seven circadian rhythm genes PRR1, PRR2, ELF3, FT, PHYB, GI and ZTL were obtained. PRR1 gene was positively correlated with flavonoids accumulation (r≥0.7). The full length of PRR1 is 3 201 bp, encoding 421 amino acids, which is highly homologous with rice OsPRR73 gene and named as CtPRR1 (GenBank accession number: MW492035). CtPRR1 was mainly expressed in flowers, and the expression level increased in the daytime and declined in the evening gradually. Correspondingly, the content of flavonoids showed an opposite variation. Both of them displayed a circadian rhythm with a negative correlation (r≥−0.7). In addition, 2 heat shock proteins along with 3 AP2 transcription factors interacting with CtPRR1 protein were obtained via Yeast Two-Hybrid experiment. Conclusion CtPRR1 negatively regulated the safflower flavonoids accumulation in a circadian rhythm way, which may be affected by these interacting proteins. -
Key words:
- Carthami flos /
- circadian gene /
- flavonoids /
- interaction protein
-
表 1 全长克隆引物
引物名称 引物序列 PRR1-F ACCTCAAGGGCCACTGGTTC PRR1-R GTAACAAAAGACTTTCTGAA 表 2 UHPLC-MS分析方法
色谱条件 质谱条件 色谱柱:Agilent ZORBAX C18(3.5 μm,2.1 mm ×100 mm) 毛细管电压:3.5 kV 流动相:A:0.1%甲酸-水,B:乙腈 锥孔电压:60 V 流速:0.4 ml/min 离子源温度:350 ℃ 柱温:40 ℃ 碎裂电压:120 V 进样量:0.4 ml 采集范围:50~1 100 m/z 梯度洗脱:0 min,5%B; 4 min,20%B;
6 min, 21% B; 9 min, 26% B; 11 min, 40% B; 15 min, 80% B; 17 min, 95% B表 3 酵母双杂筛选互作蛋白
序号 蛋白家族名称 预测功能 1 热休克蛋白 广泛参与植物响应外界环境胁迫,生物合成以及生长发育相关的转录调控 2 热休克蛋白 3 AP2转录因子 4 AP2转录因子 5 AP2转录因子 -
[1] 王升, 蒋待泉, 康传志, 等. 药用植物次生代谢在中药材生态种植中的作用及利用[J]. 中国中药杂志, 2020, 45(9):2002-2008. [2] DUNLAP J C. Molecular bases for circadian clocks[J]. Cell,1999,96(2):271-290. doi: 10.1016/S0092-8674(00)80566-8 [3] BENDIX C, MARSHALL C M, HARMON F G. Circadian clock genes universally control key agricultural traits[J]. Mol Plant,2015,8(8):1135-1152. doi: 10.1016/j.molp.2015.03.003 [4] 徐小冬, 谢启光. 植物生物钟研究的历史回顾与最新进展[J]. 自然杂志, 2013, 35(2):118-126. [5] BECKWITH E J, YANOVSKY M J. Circadian regulation of gene expression: at the crossroads of transcriptional and post-transcriptional regulatory networks[J]. Curr Opin Genet Dev,2014,27:35-42. doi: 10.1016/j.gde.2014.03.007 [6] ITO S, MATSUSHIKA A, YAMADA H, et al. Characterization of the APRR9 pseudo-response regulator belonging to the APRR1/TOC1 quintet in Arabidopsis thaliana[J]. Plant Cell Physiol,2003,44(11):1237-1245. doi: 10.1093/pcp/pcg136 [7] COCKRAM J, THIEL T, STEUERNAGEL B, et al. Genome dynamics explain the evolution of flowering time CCT domain gene families in the Poaceae[J]. PLoS One,2012,7(9):e45307. doi: 10.1371/journal.pone.0045307 [8] 陈华夏, 申国境, 王磊, 等. 4个物种CCT结构域基因家族的序列进化分析[J]. 华中农业大学学报, 2010, 29(6):669-676. [9] MATSUSHIKA A, MAKINO S, KOJIMA M, et al. Circadian waves of expression of the APRR1/TOC1 family of pseudo-response regulators in Arabidopsis thaliana: insight into the plant circadian clock[J]. Plant Cell Physiol,2000,41(9):1002-1012. doi: 10.1093/pcp/pcd043 [10] KURUP S, JONES H D, HOLDSWORTH M J. Interactions of the developmental regulator ABI3 with proteins identified from developing Arabidopsis seeds[J]. Plant J,2000,21(2):143-155. doi: 10.1046/j.1365-313x.2000.00663.x [11] LEGNAIOLI T, CUEVAS J, MAS P. TOC1 functions as a molecular switch connecting the circadian clock with plant responses to drought[J]. EMBO J,2009,28(23):3745-3757. doi: 10.1038/emboj.2009.297 [12] GRUNDY J, STOKER C, CARRÉ I A. Circadian regulation of abiotic stress tolerance in plants[J]. Front Plant Sci,2015,6:648. [13] SEO P J, MAS P. STRESSing the role of the plant circadian clock[J]. Trends Plant Sci,2015,20(4):230-237. doi: 10.1016/j.tplants.2015.01.001 [14] SALOMÉ P A, WEIGEL D, MCCLUNG C R. The role of the Arabidopsis morning loop components CCA1, LHY, PRR7, and PRR9 in temperature compensation[J]. Plant Cell,2010,22(11):3650-3661. doi: 10.1105/tpc.110.079087 [15] NAKAMICHI N, KIBA T, HENRIQUES R, et al. PSEUDO-RESPONSE REGULATORS 9, 7, and 5 are transcriptional repressors in the Arabidopsis circadian clock[J]. Plant Cell,2010,22(3):594-605. doi: 10.1105/tpc.109.072892 [16] MIZUNO T. Two-component phosphorelay signal transduction systems in plants: from hormone responses to circadian rhythms[J]. Biosci Biotechnol Biochem,2005,69(12):2263-2276. doi: 10.1271/bbb.69.2263 [17] 王晗. 生物钟生物学及其研究进展[J]. 生命科学, 2015, 27(11):1313-1319. [18] 杨志福, 梅其炳, 蒋永培. 红花有效成分及药理作用[J]. 西北药学杂志, 2001, 16(3):131-133. doi: 10.3969/j.issn.1004-2407.2001.03.025 [19] 刘飞. 红花黄酮类化合物生物合成途径关键酶基因的克隆与功能验证[D]. 上海: 第二军医大学, 2014. [20] GUO D D, GAO Y, LIU F, et al. Integrating molecular characterization and metabolites profile revealed CtCHI1& apos; s significant role in Carthamus tinctorius L[J]. BMC Plant Biol,2019,19(1):1-13. doi: 10.1186/s12870-018-1600-2 [21] BROWN B A, JENKINS G I. UV-B signaling pathways with different fluence-rate response profiles are distinguished in mature Arabidopsis leaf tissue by requirement for UVR8, HY5, and HYH[J]. Plant Physiol,2008,146(2):576-588. [22] FARRÉ E M, WEISE S E. The interactions between the circadian clock and primary metabolism[J]. Curr Opin Plant Biol,2012,15(3):293-300. doi: 10.1016/j.pbi.2012.01.013 [23] 何贝轩, 郭丹丹, 贾鑫磊, 等. 昼夜节律钟调控花青素类成分的生物合成[J]. 药学学报, 2018, 53(3):345-355.