留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式

陆丹迎 程少禹 章颖佳 刘志高 金梦婷 董彬 张寿洲 彭豪 戴梦怡 王卓为 赵宏波 申亚梅

陆丹迎, 程少禹, 章颖佳, 等. 景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式[J]. 浙江农林大学学报, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
引用本文: 陆丹迎, 程少禹, 章颖佳, 等. 景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式[J]. 浙江农林大学学报, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
LU Danying, CHENG Shaoyu, ZHANG Yingjia, et al. Bioinformatics analysis of PIF transcription factors in Magnolia sinostellata and expression pattern analysis under extreme shading condition[J]. Journal of Zhejiang A&F University, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
Citation: LU Danying, CHENG Shaoyu, ZHANG Yingjia, et al. Bioinformatics analysis of PIF transcription factors in Magnolia sinostellata and expression pattern analysis under extreme shading condition[J]. Journal of Zhejiang A&F University, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488

景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式

DOI: 10.11833/j.issn.2095-0756.20200488
基金项目: “十三五”浙江省林木新品种选育专项(2016C02056-1);浙江省重点研发项目(2019C02023)
详细信息
    作者简介: 陆丹迎(ORCID: 0000-0002-4182-0131),从事园林植物与遗传育种研究。E-mail: ludanying0713@163.com
    通信作者: 申亚梅(ORCID: 0000-0002-3079-2817),教授,从事园林植物与遗传育种研究。E-mail: yameishen@zafu.edu.cn
  • 中图分类号: S718.3

Bioinformatics analysis of PIF transcription factors in Magnolia sinostellata and expression pattern analysis under extreme shading condition

  • 摘要:   目的  群落所造成的遮阴是导致景宁木兰Magnolia sinostellata濒危的重要因素之一。PIF家族转录因子在光信号传导和植物生长发育中起到重要作用。对PIF家族转录因子进行系统分析和研究,为探究其在景宁木兰光信号转导机制中的作用奠定基础。  方法  从景宁木兰转录组数据中鉴定获得PIF家族转录因子并进行生物信息学分析,利用实时荧光定量聚合酶链式反应(qRT-PCR)技术对其在极端遮阴条件下的表达模式进行分析。  结果  从景宁木兰转录组中共筛选出9个MsPIFs转录因子基因,其编码的蛋白质长度为188~735 个氨基酸,蛋白质大小为20 314.56~78 957.02 Da,理论等电点范围为5.18~8.22。MsPIFs基因编码的蛋白质均为不稳定蛋白质,所有蛋白质均为亲水性蛋白质。亚细胞定位预测结果显示所有蛋白质均定位于细胞核。9个蛋白质均具有丝氨酸(Ser)、苏氨酸(Thr)和酪氨酸(Try)磷酸化位点。qRT-PCR结果表明:极端遮阴条件下,9个MsPIFs家族基因表达均发生不同程度的变化。其中,MsbHLH23的表达变化较其他基因更为明显,遮阴处理5和10 d时的表达量分别上调为对照的52.77与20.03倍。  结论  景宁木兰PIF转录因子家族均能响应遮阴,为后续对MsPIFs进行生物学功能鉴定奠定了基础。图7表3参32
  • 景宁木兰Magnolia sinostellata是木兰科Magnoliaceae木兰属Magnolia植物,属于典型的极小种群濒危树种。景宁木兰幼苗生长需要一定的荫庇环境[1],但过度遮阴会抑制景宁木兰幼苗的光合速率[2]。自然状态下,群落所形成的遮阴是影响群落内濒危植物生长发育的重要因素[3]。景宁木兰多生长于光照充足的针叶群落,在针阔混交群落和阔叶群落分布较少[4]。据此推断,群落所形成的遮阴会削弱景宁木兰的生长势,是导致景宁木兰濒危的重要因素。光是影响植物生长发育的重要环境因子。植物的整个生命周期都受到不断变化的光环境影响。光敏色素是植物重要的光感受器,能够感受环境中光信号的变化。光敏色素作用因子(phytochrome interacting factors, PIFs)能够感知红光和远红光信号,在光信号传导和植物生长发育中起到重要作用。光信号在通过信号转导后,能够诱导一系列生理生化反应。光敏色素作用因子是光信号转导中的重要调控因子,属于bHLH超家族的第15亚族,其家族成员都具有高度保守的bHLH结构域。bHLH结构域由约15个氨基酸的碱性结构域(basic region)和60个氨基酸的HLH(helix-loop-helix)结构域组成[5]。目前对PIF家族转录因子的研究主要集中在模式植物拟南芥Arabidopsis thaliana中。PIF转录因子广泛参与了拟南芥的光调控反应,并在光敏色素介导的光信号网络调控中处于主导地位[6]。黑暗条件下,AtPIF1能够抑制拟南芥种子的萌发,AtPIF3、AtPIF4、AtPIF5能够促进叶绿素分解并抑制叶绿体活性,从而最终导致叶片衰老[7]。遮阴条件下,AtPIF1能够通过调控AtHB1的表达来促进下胚轴的伸长[8]AtPIF4、AtPIF5和AtPIF7通过促进AtFTAtTSF的表达从而促进拟南芥开花[9]AtPIF4、AtBZR1和AtARF6通过整合生长素信号与BR和光敏色素通路来协同促进拟南芥下胚轴的伸长[10]。由此可见,PIF家族转录因子是遮阴条件下光信号转导途径中的关键调控因子。据此推测,PIF家族转录因子在光照强度对木兰属濒危树种生长发育的影响及光受体信号转导途径中起到重要作用。因此,本研究对景宁木兰PIF家族转录因子进行分析鉴定,探究其在极端遮阴条件下的表达模式,从而为揭示光介导的景宁木兰生长发育的分子机制、群落引起的遮阴而导致景宁木兰濒危的内在机制奠定基础,最终为珍稀植物景宁木兰的保育提供理论依据。

    景宁木兰来自浙江省杭州市青山湖花园中心苗圃培育的3年生扦插苗。2019年6月下旬,选择长势一致且无病虫害的扦插苗作为材料,置于温度为(25±2) ℃、光照强度为2 000~2 500 lx、相对湿度为40%~60%的人工气候室中进行培养,光照与黑暗时间各12 h,定期浇水保持土壤湿润。

    3层遮阳网形成的光照强度是景宁木兰自然生境光照条件之一,会显著削弱景宁木兰的生长势[1]。2019年7月上旬开始遮阴处理,本研究于浙江省杭州市临安区浙江农林大学人工气候室内,用3层黑色遮阳网及竹竿搭建极端遮阴处理(25%光照,ST)装置,以模拟遮阴条件下的景宁木兰野外生存环境。以100%全光照为对照(ck)。所有盆栽苗叶片之间无重叠。每个处理4个生物学重复。分别自实验开始处理第0、1、3、5、10、20、30天时,拍照记录植株表型变化、采集景宁木兰植株的中层叶片装入锡纸袋中,并迅速置于液氮中速冻,随后转入−80 ℃冰箱中储存,用于植物RNA的提取及后续实时荧光定量聚合酶链式反应(qRT-PCR)实验。

    1.3.1   景宁木兰PIF家族转录因子的筛选

    基于实验室已有的转录组测序数据,在NR、NT、Swiss-Prot这3个数据库中进行注释,初步注释到46个PIF转录因子相关基因,将这些基因编码蛋白质逐个进行NCBI Blast和NCBI CDD预测,去除重复序列及冗余转录本后,最终得到9个PIF家族转录因子的蛋白质序列。利用MEGA 5.0对景宁木兰PIF转录因子序列与拟南芥PIF转录因子序列进行比对分析,根据TOLEDO等[5]命名拟南芥PIF家族的方法对9个景宁木兰PIF转录因子进行命名。拟南芥PIF转录因子序列从拟南芥转录因子数据库(https://www.arabidopsis.org/)中获得。

    1.3.2   景宁木兰PIF家族转录因子的理化性质分析

    采用Compute pI/Mw在线工具预测景宁木兰PIF转录因子蛋白质的等电点及分子量,采用在线网站Cell-PLoc 2.0预测亚细胞定位,利用NetPhos分析蛋白质磷酸化位点;利用MEME分析蛋白质保守结构域;利用DNAMAN比对分析保守结构域;利用在线软件heatmapper(http://www2.heatmapper.ca/expression/)绘制热图;为研究PIF转录因子家族的进化关系,利用MEGA5.0将9个景宁木兰MsPIF蛋白质序列、15个拟南芥AtPIF蛋白质序列[11]、6个玉米Zea mays ZmPIF蛋白质序列[12]、10个杨树Populus trichocar PtPIF蛋白质序列(http://www.Phytozome.net/poplar.php)[13]进行序列比对,并用邻接法构建系统进化树,进行Bootstrap测试,重复设置为1 000 次。

    1.3.3   极端遮阴条件下景宁木兰PIF家族转录因子的表达模式

    采用植物RNA提取试剂盒提取样品总RNA。取1 μg RNA进行第1链 cDNA的合成。所用定量PCR试剂为BCG qPCR Master Mix。根据景宁木兰转录组数据设计特异性引物(表1)。以景宁木兰的EF1-α基因作为内参基因[1],利用Light Cycler 480II(Roche)仪器进行目的基因qRT-PCR表达分析。反应程序结束后,65~95 ℃每隔0.2 s作溶解曲线,采用${2^{ - \Delta \Delta {C_{\rm{t}}}}} $法计算目的基因的相对表达量。每个样品重复3次。

    表 1  实时荧光定量PCR引物
    Table 1  Primers used for qRT-PCR analysis
    基因名称正向(反向)引物序列(5'→3')退火温度/℃扩增效率/%相关系数R2
    MsbHLH68F-CCTCAGGCTCTGTCATTGGG59.590.860.998
    R-AGCGGTTGGGCTTCTTCACG59.5
    MsPIF4F-GCCGAAGAATCTGAGTGCCA57.496.980.999
    R-CTATCCCTCCGTCTCCTTTC57.4
    MsbHLH1F-CTTTGTGCTATCTTCGGGAACG57.7104.480.997
    R-GGGACCCACCAATCAACGAC59.5
    MsPIF8F-AAGCGAGGATGAGAACAAGG55.493.010.983
    R-TGCGTTCGGATTGGTTATGG55.4
    MsbHLH66F-AGCAGTAACGGCACGCAGAC59.5101.570.995
    R-GAAATGGGCATGAGGCAGAG57.4
    MsPIF3F-CAGACTCAGCCGTCAACTCA57.492.930.991
    R-GCAGGCCCACTTCCACCAAT59.5
    MsPIF1F-CCACATATTCCTCCATTTGAT51.7100.900.994
    R-GAAGGCTGAGGTAGTGTTGATA55.8
    MsbHLH48F-CCGTCGAGTCTCCAGTGGTC61.696.920.993
    R-CGTTTGAAGGGAAGGTAGGG57.4
    MsbHLH23F-GCAGACAAAGGACCCAAGGA57.493.010.993
    R-CACGAGGAGTGAGGAACAAGAA57.7
    下载: 导出CSV 
    | 显示表格
    1.3.4   数据分析

    采用SPSS 19.0进行方差分析并用SigmaPlot 14绘图,对平均值采用独立样本T检验进行比较,显著水平为0.05。

    利用ExPASy对9个基因编码的氨基酸进行理化性质进行分析(表2),发现9个基因编码的蛋白质大小不等,最小的为MsPIF1,编码188个氨基酸,最大的为MsPIF3,编码735个氨基酸。9个蛋白质的理论等电点为5.18(MsbHLH68)~8.22(MsPIF8),表明不同的PIF家族蛋白质在不同的微环境中发挥不同的生物学功能。9个PIF家族蛋白质均属于不稳定、亲水性蛋白质。

    表 2  MsPIF基因家族
    Table 2  MsPIF gene family of M. sinostellata
    基因名称蛋白质分子量/Da氨基酸数/个等电点不稳定系数脂肪系数平均亲水性亚细胞定位磷酸化位点
    SerThrTyr
    MsbHLH6853 604.795065.1858.1977.51−0.447细胞核(N)47101
    MsPIF457 883.065275.9656.0564.57−0.611细胞核(N)40123
    MsbHLH174 557.306715.5253.0879.66−0.452细胞核(N)34146
    MsPIF850 048.294568.2242.0958.29−0.686细胞核(N)36204
    MsbHLH6645 976.464306.7960.4569.72−0.558细胞核(N)45141
    MsPIF378 957.027356.3167.3763.80−0.513细胞核(N)74133
    MsPIF120 314.561885.9360.7860.32−0.341细胞核(N) 4 21
    MsbHLH4840 255.003666.4359.8570.90−0.639细胞核(N)40 62
    MsbHLH2339 480.073625.8272.0770.14−0.577细胞核(N)43 92
      说明:Ser为丝氨酸,Thr为苏氨酸,Tyr为酪氨酸
    下载: 导出CSV 
    | 显示表格

    通过在线软件Cell-PLoc 2.0(http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/)对9个PIF基因家族成员的氨基酸序列进行亚细胞定位预测(表2),9个基因编码的蛋白质均定位于细胞核内,符合其转录因子的功能,可能通过调控细胞核基因的转录来发挥功能。

    蛋白质磷酸化在很多植物信号转导过程中起着重要作用[14]。蛋白质磷酸化的可逆性使得植物细胞能够快速响应外部环境变化[15]。在正常光条件下,PIF家族转录因子与光敏色素相互作用,从而导致PIF蛋白质的磷酸化和降解[16]。而在低光照条件下,PIF蛋白质的磷酸化和降解受到抑制[17]。因此,分析PIF蛋白质的磷酸化位点对研究其在信号转导中的功能有重要意义。利用NetPhos对9个PIF家族蛋白质的磷酸化位点进行分析(表2),发现9个蛋白质均具有丝氨酸(Ser)、苏氨酸(Thr)和酪氨酸(Tyr)磷酸化位点,且Ser的磷酸化位点最多,Tyr的磷酸化位点最少。其中磷酸化位点最多的是MsPIF3,为94个,MsPIF1的磷酸化位点最少,仅7个。

    进一步利用MEME在线软件分析MsPIF蛋白质[16],共鉴定到6个基序(表3)。结果显示(图1):9个蛋白质含有不同数目的基序,MsPIF4和MsPIF8含有的基序最多,为5个,而MsbHLH48中含有的基序最少,仅2个。bHLH家族蛋白质均含有bHLH基序,该基序由位于N端的basic基序(图2A)和位于C端的HLH基序(图2B)组成[5]。9个PIF家族蛋白质均含有basic基序及HLH基序。不同蛋白质的同源异型基序之间存在一定差异。MsPIF4、MsPIF3和MsPIF8中含有APB基序(active phyB-binding motif)(图2C),表明其可以与PhyB(phytochrome B)相结合,从而对下游基因的表达进行调控[16]。此外,MsPIF4和MsPIF3中还含有APA基序(active phyA-binding motif)(图2D),表明其既可以与PhyA(phytochrome A)相结合也可以与PhyB相结合。

    表 3  景宁木兰PIF家族蛋白质的主要MEME基序
    Table 3  Major MEME motif sequences in M. sinostellata PIF protein
    基序宽度/个最佳匹配
    Motif140MKALQELIPNSNKTDKASMLDEAIEYLKSLQLQVQMLSM
    Motif221RGRAAEVHNLSERRRRERINE
    Motif321MMFPGVQQYMPPMGMGMGMGM
    Motif436FEQZIVKLMEEDMGSAMQYLQGKGLCLMPIALAAAI
    Motif5 8ELLWENGQ
    Motif6 8CVPBWBTD
    下载: 导出CSV 
    | 显示表格
    图 1  MsPIF蛋白质的基序
    Figure 1  Motif of MsPIF protein
    图 2  MsPIF家族蛋白质的basic、HLH、APB和APA基序蛋白质序列比对
    Figure 2  Sequence alignment of basic, HLH, APB and APA motif protein sequences of MsPIF

    运用MEGA5对所得景宁木兰的蛋白质序列与拟南芥、玉米、杨树中的序列对比分析(图3),采用邻接法及Bootstrap分析(1 000次重复)构建系统进化树[18]。结果表明:植物PIF家族蛋白质表现出多样性,可以分为7个进化分枝。其中,9个MsPIF蛋白质分布在α、β、γ、δ、η、ε等6个进化枝。15个AtPIF蛋白质分布在α、β、γ、δ、ζ、ε等6个进化枝,其中ζ为拟南芥所特有。ZmPIF蛋白质仅分布于β和δ进化枝,10个PtPIF蛋白质分布在α、β、γ、δ、ε等5个进化枝。从进化树可以看出,MsPIF8、MsPIF4、MsPIF1、MsPIF3及MsbHLH23均与相应的拟南芥蛋白质与杨树蛋白质聚类在一起,表明其亲缘关系较为接近。MsbHLH68、MsbHLH66、MsbHLH48汇聚于η枝,为景宁木兰所特有。

    图 3  景宁木兰、拟南芥、玉米和杨树的PIF家族蛋白质的系统进化分析
    Figure 3  Phylogenetic tree analysis of MsPIF protein and AtPIF protein sequences in M. sinostellata, A. thaliana, Z. mays, P. trichocarpa

    光照是植物生长发育过程中不可或缺的因素,极端遮阴处理过程中,景宁木兰叶片表型变化十分明显(图4)。处理10 d,景宁木兰叶片出现褐色斑点;处理20 d,褐色斑点面积增大,植株开始萎蔫;处理30 d,大量叶片干枯脱落,植株萎蔫严重。

    图 4  遮阴处理过程中景宁木兰的表型变化
    Figure 4  Phenotypic changes of M. sinostellata during shading treatment

    PCR扩增琼脂糖凝胶电泳检测发现9个基因均扩出预期条带(图5)。qRT-PCR检测结果(图6)表明:遮阴处理1 d时,MsbHLH23和MsbHLH66的表达水平显著上调,MsbHLH1、MsbHLH48、MsPIF8、MsPIF3的表达水平均显著下调,其余基因表达水平无显著变化。遮阴处理3 d时,MsPIF8的表达显著低于对照,其余基因的表达均显著上调。遮阴处理5 d时,MsPIF8的表达仍显著降低,MsbHLH68、MsbHLH23、MsbHLH66、MsbHLH1、MsbHLH48的表达显著上调,其余基因的表达无显著差异。遮阴处理10 d时,MsPIF8、MsPIF7、MsPIF3的表达均显著下调,MsbHLH68、MsbHLH23、MsbHLH66、MsPIF4的表达显著上调。遮阴处理20 d时,MsPIF8、MsbHLH23、MsPIF1的表达相对于对照显著下调,MsbHLH68和MsbHLH66的表达显著上调,其余基因无显著差异。遮阴处理30 d时,MsPIF8、MsbHLH23的表达相对于对照显著下调,MsbHLH68、MsbHLH1、MsbHLH66、MsbHLH48、MsPIF3的表达显著上调。其中,MsbHLH23的表达变化比其他基因更为明显,遮阴处理5和10 d时的表达量分别上调为对照的52.77与20.03倍。

    图 5  qRT-PCR扩增产物特异性
    Figure 5  Specificity of qRT-PCR amplification products
    图 6  遮阴条件下景宁木兰PIF转录因子的qRT-PCR分析
    Figure 6  qRT-PCR analysis of PIF transcription factor under shading treatment of M. sinostellata

    基因的表达模式与基因功能关系密切,基于获得的景宁木兰处理组的qRT-PCR数据,采用heatmap在线软件对这些基因的差异表达进行分层聚类分析(图7)。结果发现:MsPIF1、MsPIF3、MsPIF4、MsPIF8的表达模式十分相似,其表达均在处理3 d时明显升高,随后逐渐下降。MsbHLH48的表达在处理0~30 d逐渐上调,其中处理30 d时的表达上调最为明显。MsbHLH1、MsbHLH68、MsbHLH66、MsbHLH23的表达模式均为先逐渐上调后逐渐下调,其中MsbHLH1和MsbHLH23的表达在处理5 d显著上调,MsbHLH68和MsbHLH66的表达分别在处理10和20 d明显上调。

    图 7  长期极端遮阴条件下景宁木兰PIF转录因子的表达模式
    Figure 7  Pattern of expression of PIF family during shading treatment period of M. sinostellata

    PIF家族转录因子广泛存在于植物中。目前PIF家族转录因子在拟南芥中研究的最为透彻,共有15个[5]。在其他物种中,杨树至少存在10个[13],玉米有6个[12],水稻Oryza sativa有6个[19]。本研究从景宁木兰中仅筛选到9个,这一方面是由于景宁木兰本身具有特异性,另一方面可能是因为转录组数据库还不够丰富造成的。NCBI blast及NCBI CDD预测结果表明:9个PIF蛋白质均含bHLH结构域,属于PIF转录因子家族。同源性较高的蛋白质通常功能也相似[18]。因此,本研究用所得9个景宁木兰PIF家族转录因子编码蛋白质与6个ZmPIFs、10个PtPIFs 和15个功能明确的AtPIFs编码蛋白质作系统进化树,主要根据与15个拟南芥PIF家族转录因子编码蛋白质进化关系的远近推测景宁木兰PIF转录因子的功能。AtPIF3与PtPIF3a、PtPIF3b同源性较高,且与AtPIF3和ZmPIF3位于同一进化枝,据此推测MsPIF3与AtPIF3功能相似。MsPIF4与PtPIF4、ZmPIF4、AtPIF4位于同一进化枝。AtPIF3蛋白质在植株刚刚暴露到光下时会急剧减少,在黑暗条件下会重新积累,在拟南芥幼苗阶段,AtPIF3在光诱导的叶绿体发育中起到积极的作用[20]。AtPIF3、AtPIF4和AtPIF5在年龄引起的和黑暗引起的植物衰老中起到重要作用,其中AtPIF4通过调控叶绿素降解和维持叶绿体活性的相关基因从而促进植株衰老[7]。AtPIF4能够调控下游基因从而促进拟南芥开花[21]

    此外,AtPIF4能与AtBZR1共同调节下游光响应的相关基因。AtPIF7与MsbHLH1同源性较高,推测其功能可能类似。遮阴条件下,AtPIF7是调控拟南芥茎伸长的重要因子[17]。AtPIF7是调节植物去黄化的弱负反馈因子,与PhyB结合的过程中通常不伴随着可以监测到的磷酸化和分解[22]。MsPIF8与AtPIF8同源,推测其功能相似。与AtPIF3不同,AtPIF8蛋白质在远红光下的积累比在黑暗条件中或红光条件下更多,AtPIF8能够抑制由PhyA介导的种子萌发、下胚轴伸长等生命活动[23]。MsbHLH23与AtbHLH24功能相似,AtbHLH24在抑制种子萌发和维持种子休眠状态中具有重要作用[24]。AtbHLH24能够激活细胞分裂素信号,同时激活生长素生物合成和雄蕊内部结构域的转运基因,并且在花瓣和果实发育中起到重要作用[25]。MsbHLH68,MsbHLH66和MsbHLH48共同汇集于同一进化枝,其功能还需要进一步挖掘。但是,即使是同一分支内蛋白质功能也可能存在差异,并且不同物种之间具有较大的差异[18]。因而关于MsPIF家族蛋白质的具体功能还需要通过后续研究。

    群落层片造成的遮阴会降低下层植被的生长率与成活率[26]。本研究极端遮阴处理组的景宁木兰在处理过程中逐步萎蔫死亡。光敏色素作用因子(PIFs)在植物的生长发育中起到重要调节作用。PIFs作为光信号转导途径中的负调控因子,对植物生长发育的多个过程都有调控作用[27-32]。PIF家族转录因子是光敏色素介导的光信号转导途径中的负反馈因子[28]。遮阴条件下,PIF蛋白质大量积累,从而调控植物发生一系列生理生化反应。考虑到植物的个体差异和基因表达差异,本研究的每个取样时间点都作了相应的空白对照。本研究中,长期极端遮阴条件下,景宁木兰的PIF类家族转录因子的表达水平多数显著上调,只有MsPIF8的表达水平显著下调。在拟南芥中,与其他PIF家族成员相比,AtPIF8能够与PhyB结合但效果微弱[23]MsPIF8的表达模式与前人的研究结论存在差异,可能与物种特异性相关,有关MsPIF8的功能还需要进一步研究。MsPIF1、MsPIF3和MsPIF4的表达模式与前人的研究结论[29]吻合。黑暗条件下,AtPIF1的表达上调,从而抑制叶绿素的合成与种子萌发;AtPIF3、AtPIF4的表达显著上调,从而促进叶绿素的降解与叶片衰老[29]。植物群落造成的遮阴使得AtPIF4的表达上调,从而使成花诱导相关基因的表达上调,最终诱导拟南芥提早开花[21]。PIF功能十分复杂,景宁木兰PIF家族转录因子在遮阴调控网络中发挥的作用还需要进一步研究。

  • 图  1  MsPIF蛋白质的基序

    Figure  1  Motif of MsPIF protein

    图  2  MsPIF家族蛋白质的basic、HLH、APB和APA基序蛋白质序列比对

    Figure  2  Sequence alignment of basic, HLH, APB and APA motif protein sequences of MsPIF

    图  3  景宁木兰、拟南芥、玉米和杨树的PIF家族蛋白质的系统进化分析

    Figure  3  Phylogenetic tree analysis of MsPIF protein and AtPIF protein sequences in M. sinostellata, A. thaliana, Z. mays, P. trichocarpa

    图  4  遮阴处理过程中景宁木兰的表型变化

    Figure  4  Phenotypic changes of M. sinostellata during shading treatment

    图  5  qRT-PCR扩增产物特异性

    Figure  5  Specificity of qRT-PCR amplification products

    图  6  遮阴条件下景宁木兰PIF转录因子的qRT-PCR分析

    Figure  6  qRT-PCR analysis of PIF transcription factor under shading treatment of M. sinostellata

    图  7  长期极端遮阴条件下景宁木兰PIF转录因子的表达模式

    Figure  7  Pattern of expression of PIF family during shading treatment period of M. sinostellata

    表  1  实时荧光定量PCR引物

    Table  1.   Primers used for qRT-PCR analysis

    基因名称正向(反向)引物序列(5'→3')退火温度/℃扩增效率/%相关系数R2
    MsbHLH68F-CCTCAGGCTCTGTCATTGGG59.590.860.998
    R-AGCGGTTGGGCTTCTTCACG59.5
    MsPIF4F-GCCGAAGAATCTGAGTGCCA57.496.980.999
    R-CTATCCCTCCGTCTCCTTTC57.4
    MsbHLH1F-CTTTGTGCTATCTTCGGGAACG57.7104.480.997
    R-GGGACCCACCAATCAACGAC59.5
    MsPIF8F-AAGCGAGGATGAGAACAAGG55.493.010.983
    R-TGCGTTCGGATTGGTTATGG55.4
    MsbHLH66F-AGCAGTAACGGCACGCAGAC59.5101.570.995
    R-GAAATGGGCATGAGGCAGAG57.4
    MsPIF3F-CAGACTCAGCCGTCAACTCA57.492.930.991
    R-GCAGGCCCACTTCCACCAAT59.5
    MsPIF1F-CCACATATTCCTCCATTTGAT51.7100.900.994
    R-GAAGGCTGAGGTAGTGTTGATA55.8
    MsbHLH48F-CCGTCGAGTCTCCAGTGGTC61.696.920.993
    R-CGTTTGAAGGGAAGGTAGGG57.4
    MsbHLH23F-GCAGACAAAGGACCCAAGGA57.493.010.993
    R-CACGAGGAGTGAGGAACAAGAA57.7
    下载: 导出CSV

    表  2  MsPIF基因家族

    Table  2.   MsPIF gene family of M. sinostellata

    基因名称蛋白质分子量/Da氨基酸数/个等电点不稳定系数脂肪系数平均亲水性亚细胞定位磷酸化位点
    SerThrTyr
    MsbHLH6853 604.795065.1858.1977.51−0.447细胞核(N)47101
    MsPIF457 883.065275.9656.0564.57−0.611细胞核(N)40123
    MsbHLH174 557.306715.5253.0879.66−0.452细胞核(N)34146
    MsPIF850 048.294568.2242.0958.29−0.686细胞核(N)36204
    MsbHLH6645 976.464306.7960.4569.72−0.558细胞核(N)45141
    MsPIF378 957.027356.3167.3763.80−0.513细胞核(N)74133
    MsPIF120 314.561885.9360.7860.32−0.341细胞核(N) 4 21
    MsbHLH4840 255.003666.4359.8570.90−0.639细胞核(N)40 62
    MsbHLH2339 480.073625.8272.0770.14−0.577细胞核(N)43 92
      说明:Ser为丝氨酸,Thr为苏氨酸,Tyr为酪氨酸
    下载: 导出CSV

    表  3  景宁木兰PIF家族蛋白质的主要MEME基序

    Table  3.   Major MEME motif sequences in M. sinostellata PIF protein

    基序宽度/个最佳匹配
    Motif140MKALQELIPNSNKTDKASMLDEAIEYLKSLQLQVQMLSM
    Motif221RGRAAEVHNLSERRRRERINE
    Motif321MMFPGVQQYMPPMGMGMGMGM
    Motif436FEQZIVKLMEEDMGSAMQYLQGKGLCLMPIALAAAI
    Motif5 8ELLWENGQ
    Motif6 8CVPBWBTD
    下载: 导出CSV
  • [1] YU Qin, SHEN Yamei, WANG Qianying, et al. Light deficiency and waterlogging affect chlorophyll metabolism and photosynthesis in Magnolia sinostellata [J]. Trees, 2019, 33(1): 11 − 22.
    [2] 俞芹, 王倩颖, 刘志高, 等. 光强与水分处理下景宁木兰光合光响应模型拟合比较[J]. 生态学杂志, 2018, 37(3): 898 − 905.

    YU Qin, WANG Qianying, LIU Zhigao, et al. Comparison of the light response models of photosynthesis in leaves of Magnolia sinostellata under different light intensity and moisture conditions [J]. Chin J Ecol, 2018, 37(3): 898 − 905.
    [3] TAMAKI I, NOMURA K, NOMURA R, et al. Evaluation of a field experiment for the conservation of a Magnolia stellata stand using clear-cutting [J]. Landscape Ecol Eng, 2018, 14(2): 269 − 276.
    [4] 余泽智, 陈翔翔, 卢璐, 等. 玉景宁玉兰种群分布与群落结构研究[J]. 浙江林业科技, 2015, 35(3): 47 − 52.

    YU Zezhi, CHEN Xiangxiang, LU Lu, et al. Distribution and community structure ofMagnolia sinostellata [J]. J Zhejiang For Sci Technol, 2015, 35(3): 47 − 52.
    [5] TOLEDO-ORTIZ G, HUQ E, QUAIL P H. The Arabidopsis basic/helix-loop-helix transcription factor family [J]. Plant Cell, 2003, 15(8): 1749 − 1770.
    [6] CASTILLON A, SHEN Hui, HUQ E. Phytochrome interacting factors: central players in phytochrome-mediated light signaling networks [J]. Trends Plant Sci, 2007, 12(11): 514 − 521.
    [7] SONG Yi, YANG Chuangwei, GAO Shan, et al. Age-triggered and dark-induced leaf senescence require the bHLH transcription factors PIF3, 4, and 5 [J]. Mol Plant, 2014, 7(12): 1776 − 1787.
    [8] CAPELLA M, RIBONE P A, ARCE A L, et al. Arabidopsis thaliana HomeoBox 1 (AtHB1), a homedomain-leucine zipper I (HD-Zip I) transcription factor, is regulated by PHYTOCHROME-INTERACTING FACTOR 1 to promote hypocotyl elongation [J]. New Phytol, 2015, 207(3): 669 − 682.
    [9] GALVĀO V C, FIORUCCI A S, TREVISAN M, et al. PIF transcription factors link a neighbor threat cue to accelerated reproduction in Arabidopsis[J]. Nat Commun, 2019, 10(1): 4005. doi: 10.1038/s41467-019-11882-7.
    [10] EUNKYOO O, ZHU Jiaying, BAI Mingyi, et al. Cell elongation is regulated through a central circuit of interacting transcription factors in the Arabidopsis hypocotyl[J]. eLife, 2014, 3: e03031. doi: 10.7554/eLife.03031.
    [11] HASAN M, RASHID M, KHATUN S, et al. Computational identification of microbial phosphorylation sites by the enhanced characteristics of sequence information [J]. Sci Rep, 2019, 9: 8258. doi: 10.1038/s41598-019-44548-x.
    [12] 陈亮, 崔芬芬, 王勇飞, 等. 基于转录组金银花WRKY转录因子的挖掘与分析[J]. 分子植物育种, 2019, 17(6): 1780 − 1787.

    CHEN Liang, CUI Fenfen, WANG Yongfei, et al. Mining and analysis of WRKY transcription factors in transcriptome-based Lonicera japonica [J]. Mol Plant Breed, 2019, 17(6): 1780 − 1787.
    [13] ESSER D, HOFFMANN L, PHAM T K, et al. Protein phosphorylation and its role in archaeal signal transduction [J]. FEMS Microbiol Rev, 2016, 40(5): 625 − 647.
    [14] PHAM V N, KATHARE P K, HUQ E. Phytochromes and phytochrome interacting factors [J]. Plant Physiol, 2018, 176(2): 1025 − 1038.
    [15] HUANG Xu, ZHANG Qian, JIANG Yupei, et al. Shade-induced nuclear localization of PIF7 is regulated by phosphorylation and 14-3-3 proteins in Arabidopsis[J]. eLife, 2018, 7: e31636. doi: 10.7554/eLife.31636.
    [16] 庄黎丽, 王剑, 杨志民. 基于转录组数据库的高羊茅HD-Zip I转录因子的鉴定及表达模式解析[J]. 草业学报, 2017, 27(3): 67 − 77.

    ZHUANG Lili, WANG Jian, YANG Zhimin. Transcriptome-wide identification and expression analysis of HD-Zip I transcription factors in Festuca arundinacea [J]. Acta Prat Sin, 2017, 27(3): 67 − 77.
    [17] 徐向东, 任逸秋, 张利, 等. 杨树PIF基因家族成员表达模式研究[J]. 林业科学研究, 2018, 31(2): 19 − 25.

    XU Xiangdong, REN Yiqiu, ZHANG Li, et al. Analysis of expression pattern of PIF family members in Populus [J]. For Res, 2018, 31(2): 19 − 25.
    [18] GAO Yong, REN Xiaoyun, QIAN Jingjie, et al. The phytochrome-interacting family of transcription factors in maize (Zea mays L.): identification, evolution, and expression analysis[J]. Acta Physiol Plant, 2019, 41(1): 8. doi: 10.1007/s11738-018-2802-9.
    [19] NAKAMURA Y, KATO T, YAMASHINO T, et al. Characterization of a set of phytochrome-interacting factor-like bHLH proteins in Oryza sativa [J]. Biosci Biotechnol Biochem, 2007, 71(5): 1183 − 1191.
    [20] MONTE E, TEPPERMAN J M, AL-SADY B, et al. The phytochrome-interacting transcription factor, PIF3, acts early, selectively, and positively in light-induced chloroplast development [J]. Proc Natl Acad Sci USA, 2004, 101(46): 16091 − 16098.
    [21] KUMAR S V, LUCYSHYN D, JAEGER K E, et al. Transcription factor PIF4 controls the thermosensory activation of flowering [J]. Nature, 2012, 484(7393): 242 − 245.
    [22] LEIVAR P, MONTE E, AL-SADY B, et al. The Arabidopsis phytochrome-interacting factor PIF7, together with PIF3 and PIF4, regulates responses to prolonged red light by modulating phyB levels [J]. Plant Cell, 2008, 20(2): 337 − 352.
    [23] OH J, PARK E, SONG K, et al. PHYTOCHROME INTERACTING FACTOR8 Inhibits phytochrome a-mediated far-red light responses in Arabidopsis [J]. Plant Cell, 2020, 32(1): 186 − 205.
    [24] PENFIELD S, JOSSE E M, KANNANGARA R, et al. Cold and light control seed germination through the bHLH transcription factor SPATULA [J]. Curr Biol, 2005, 15(22): 1998 − 2006.
    [25] REYES-OLALDE J I, ZÚÑIGA-MAYO V M, SERWATOWSKA J, et al. The bHLH transcription factor SPATULA enables cytokinin signaling, and both activate auxin biosynthesis and transport genes at the medial domain of the gynoecium [J]. PLoS Genet, 2017, 13(4): e1006726. doi: 10.1371/journal.pgen.1006726.
    [26] GROSZMANN M, PAICU T, SMYTH D R. Functional domains of SPATULA, a bHLH transcription factor involved in carpel and fruit development in Arabidopsis [J]. Plant J, 2008, 55(1): 40 − 52.
    [27] LU Deliang, WANG G G, ZHANG Jinxin, et al. Converting larch plantations to mixed stands: effects of canopy treatment on the survival and growth of planted seedlings with contrasting shade tolerance [J]. For Ecol Manage, 2018, 409: 19 − 28.
    [28] DONG Jie, SUN Ning, YANG Jing, et al. The transcription factors TCP4 and PIF3 antagonistically regulate organ-specific light induction of SAUR genes to modulate cotyledon opening during de-etiolation inArabidopsis [J]. Plant Cell, 2019, 31(5): 1155 − 1170.
    [29] OH E, YAMAGUCHI S, HU Jianhong, et al. PIL5, a phytochrome-interacting bHLH protein, regulates gibberellin responsiveness by binding directly to the GAI and RGA promoters in Arabidopsis seeds [J]. Plant Cell, 2007, 19(4): 1192 − 1208.
    [30] OH E, ZHU Jiaying, WANG Zhiyong. Interaction between BZR1 and PIF4 integrates brassinosteroid and environmental responses [J]. Nat Cell Biol, 2012, 14(8): 802 − 809.
    [31] PARK E, PARK J, KIM J, et al. Phytochrome B inhibits binding of phytochrome-interacting factors to their target promoters [J]. Plant J, 2012, 72(4): 537 − 546.
    [32] 任小芸, 吴美琴, 陈建民, 等. 光敏色素作用因子PIFs参与植物激素信号转导的分子机制[J]. 植物生理学报, 2016, 52(10): 1466 − 1473.

    REN Xiaoyun, WU Meiqin, CHEN Jianmin, et al. The molecular mechanisms of phytochrome interacting factors (PIFs) in phy-tohormone signaling transduction [J]. Plant Physiol Commun, 2016, 52(10): 1466 − 1473.
  • [1] 王书伟, 周明兵.  毛竹ICE基因家族的全基因组鉴定及低温胁迫下的表达模式分析 . 浙江农林大学学报, 2024, 41(3): 568-576. doi: 10.11833/j.issn.2095-0756.20230445
    [2] 郭玉婷, 杜长霞.  黄瓜R2R3-MYB亚家族鉴定及生物信息学分析 . 浙江农林大学学报, 2024, 41(2): 286-296. doi: 10.11833/j.issn.2095-0756.20230278
    [3] 孟超敏, 耿翡翡, 卿桂霞, 张富厚, 李雪林, 刘逢举.  陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析 . 浙江农林大学学报, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
    [4] 刘俊, 李龙, 陈玉龙, 陈随清.  杜仲WOX家族基因鉴定及在叶片发育中的表达 . 浙江农林大学学报, 2023, 40(1): 1-11. doi: 10.11833/j.issn.2095-0756.20210725
    [5] 郑晨璐, 王倩颖, 陆丹迎, 申亚梅, 马晶晶, 刘璐, 王云.  高温胁迫下景宁木兰淀粉与蔗糖代谢途径转录分析 . 浙江农林大学学报, 2023, 40(1): 55-63. doi: 10.11833/j.issn.2095-0756.20220170
    [6] 刘俊, 李龙, 陈玉龙, 刘燕, 吴耀松, 任闪闪.  杜仲CONSTANS-like全基因组鉴定、系统进化及表达模式分析 . 浙江农林大学学报, 2022, 39(3): 475-485. doi: 10.11833/j.issn.2095-0756.20210385
    [7] 王绍良, 张雯宇, 高志民, 周明兵, 杨克彬, 宋新章.  毛竹磷转运蛋白Ⅰ家族基因鉴定及表达模式 . 浙江农林大学学报, 2022, 39(3): 486-494. doi: 10.11833/j.issn.2095-0756.20210471
    [8] 赖梦霞, 杜长霞, 樊怀福.  黄瓜AQP基因家族的鉴定与生物信息学分析 . 浙江农林大学学报, 2022, 39(2): 318-328. doi: 10.11833/j.issn.2095-0756.20210361
    [9] 陈雪冰, 刘聪, 程赫, 姜廷波, 夏德安, 魏志刚.  毛果杨ZHD家族全基因组水平鉴定及在干旱胁迫下的表达分析 . 浙江农林大学学报, 2022, 39(3): 465-474. doi: 10.11833/j.issn.2095-0756.20210373
    [10] 黄奕孜, 钱旺, 邱姗, 王文新, 黄华宏, 林二培.  光皮桦AP2/ERF基因家族鉴定与表达分析 . 浙江农林大学学报, 2022, 39(6): 1183-1193. doi: 10.11833/j.issn.2095-0756.20220331
    [11] 缪云锋, 周丹, 董彬, 赵宏波.  桂花OfNAC转录因子鉴定及在花开放阶段的表达分析 . 浙江农林大学学报, 2021, 38(3): 433-444. doi: 10.11833/j.issn.2095-0756.20200474
    [12] 黄元城, 郭文磊, 王正加.  薄壳山核桃全基因组LBD基因家族的生物信息学分析 . 浙江农林大学学报, 2021, 38(3): 464-475. doi: 10.11833/j.issn.2095-0756.20200454
    [13] 原晓龙, 李娟, 李云琴, 王毅.  1个含有SDR结构域PKS/NRPS基因的克隆 . 浙江农林大学学报, 2019, 36(6): 1247-1253. doi: 10.11833/j.issn.2095-0756.2019.06.024
    [14] 王倩颖, 常鹏杰, 申亚梅, 张超, 董彬, 时宝柱.  景宁木兰热胁迫下实时荧光定量PCR内参基因的筛选 . 浙江农林大学学报, 2019, 36(5): 935-942. doi: 10.11833/j.issn.2095-0756.2019.05.012
    [15] 程占超, 侯丹, 马艳军, 高健.  毛竹生长素反应因子基因的生物信息学分析及差异表达 . 浙江农林大学学报, 2017, 34(4): 574-580. doi: 10.11833/j.issn.2095-0756.2017.04.002
    [16] 刘健健, 刘俊丽, 季敏杰, 陈家栋, 杨晓峰, 陈爱群.  番茄质膜H+-ATPase家族基因的鉴定和表达分析 . 浙江农林大学学报, 2016, 33(5): 734-741. doi: 10.11833/j.issn.2095-0756.2016.05.002
    [17] 安苗苗, 刘静, 郦元, 周明兵.  花叶矢竹转录组中的转座子表达分析 . 浙江农林大学学报, 2016, 33(6): 935-943. doi: 10.11833/j.issn.2095-0756.2016.06.003
    [18] 王庆灵, 刘文鑫, 赵嘉平.  山海关杨PdERF-18转录因子的表达特征分析 . 浙江农林大学学报, 2014, 31(5): 716-723. doi: 10.11833/j.issn.2095-0756.2014.05.009
    [19] 彭沙沙, 童再康, 黄华宏, 周厚君, 时剑, 林二培.  杉木纤维素合成酶类似蛋白基因ClCslD1的克隆及其生物信息学分析 . 浙江农林大学学报, 2012, 29(1): 1-6. doi: 10.11833/j.issn.2095-0756.2012.01.001
    [20] 吴家胜, 汪旭升.  数量性状位点(QTLs)内候选基因的生物信息学分析方法 . 浙江农林大学学报, 2008, 25(1): 104-108.
  • 期刊类型引用(2)

    1. 徐建恩,石从广,申亚梅,仇建习,李因刚,诸葛菲,周伟琴,杨伟丽. 景宁木兰保护生物学研究进展. 浙江林业科技. 2022(03): 97-100 . 百度学术
    2. 郝新迪,许梦男,董波波,陈丽,刘超杰,范双喜,韩莹琰. 基于转录组信息的叶用莴苣MADS-box转录因子家族分析. 北京农学院学报. 2021(04): 14-19 . 百度学术

    其他类型引用(1)

  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20200488

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2021/3/445

图(7) / 表(3)
计量
  • 文章访问数:  1608
  • HTML全文浏览量:  462
  • PDF下载量:  79
  • 被引次数: 3
出版历程
  • 收稿日期:  2020-07-28
  • 修回日期:  2021-01-23
  • 网络出版日期:  2021-06-09
  • 刊出日期:  2021-06-09

景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式

doi: 10.11833/j.issn.2095-0756.20200488
    基金项目:  “十三五”浙江省林木新品种选育专项(2016C02056-1);浙江省重点研发项目(2019C02023)
    作者简介:

    陆丹迎(ORCID: 0000-0002-4182-0131),从事园林植物与遗传育种研究。E-mail: ludanying0713@163.com

    通信作者: 申亚梅(ORCID: 0000-0002-3079-2817),教授,从事园林植物与遗传育种研究。E-mail: yameishen@zafu.edu.cn
  • 中图分类号: S718.3

摘要:   目的  群落所造成的遮阴是导致景宁木兰Magnolia sinostellata濒危的重要因素之一。PIF家族转录因子在光信号传导和植物生长发育中起到重要作用。对PIF家族转录因子进行系统分析和研究,为探究其在景宁木兰光信号转导机制中的作用奠定基础。  方法  从景宁木兰转录组数据中鉴定获得PIF家族转录因子并进行生物信息学分析,利用实时荧光定量聚合酶链式反应(qRT-PCR)技术对其在极端遮阴条件下的表达模式进行分析。  结果  从景宁木兰转录组中共筛选出9个MsPIFs转录因子基因,其编码的蛋白质长度为188~735 个氨基酸,蛋白质大小为20 314.56~78 957.02 Da,理论等电点范围为5.18~8.22。MsPIFs基因编码的蛋白质均为不稳定蛋白质,所有蛋白质均为亲水性蛋白质。亚细胞定位预测结果显示所有蛋白质均定位于细胞核。9个蛋白质均具有丝氨酸(Ser)、苏氨酸(Thr)和酪氨酸(Try)磷酸化位点。qRT-PCR结果表明:极端遮阴条件下,9个MsPIFs家族基因表达均发生不同程度的变化。其中,MsbHLH23的表达变化较其他基因更为明显,遮阴处理5和10 d时的表达量分别上调为对照的52.77与20.03倍。  结论  景宁木兰PIF转录因子家族均能响应遮阴,为后续对MsPIFs进行生物学功能鉴定奠定了基础。图7表3参32

English Abstract

陆丹迎, 程少禹, 章颖佳, 等. 景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式[J]. 浙江农林大学学报, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
引用本文: 陆丹迎, 程少禹, 章颖佳, 等. 景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式[J]. 浙江农林大学学报, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
LU Danying, CHENG Shaoyu, ZHANG Yingjia, et al. Bioinformatics analysis of PIF transcription factors in Magnolia sinostellata and expression pattern analysis under extreme shading condition[J]. Journal of Zhejiang A&F University, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
Citation: LU Danying, CHENG Shaoyu, ZHANG Yingjia, et al. Bioinformatics analysis of PIF transcription factors in Magnolia sinostellata and expression pattern analysis under extreme shading condition[J]. Journal of Zhejiang A&F University, 2021, 38(3): 445-454. DOI: 10.11833/j.issn.2095-0756.20200488
  • 景宁木兰Magnolia sinostellata是木兰科Magnoliaceae木兰属Magnolia植物,属于典型的极小种群濒危树种。景宁木兰幼苗生长需要一定的荫庇环境[1],但过度遮阴会抑制景宁木兰幼苗的光合速率[2]。自然状态下,群落所形成的遮阴是影响群落内濒危植物生长发育的重要因素[3]。景宁木兰多生长于光照充足的针叶群落,在针阔混交群落和阔叶群落分布较少[4]。据此推断,群落所形成的遮阴会削弱景宁木兰的生长势,是导致景宁木兰濒危的重要因素。光是影响植物生长发育的重要环境因子。植物的整个生命周期都受到不断变化的光环境影响。光敏色素是植物重要的光感受器,能够感受环境中光信号的变化。光敏色素作用因子(phytochrome interacting factors, PIFs)能够感知红光和远红光信号,在光信号传导和植物生长发育中起到重要作用。光信号在通过信号转导后,能够诱导一系列生理生化反应。光敏色素作用因子是光信号转导中的重要调控因子,属于bHLH超家族的第15亚族,其家族成员都具有高度保守的bHLH结构域。bHLH结构域由约15个氨基酸的碱性结构域(basic region)和60个氨基酸的HLH(helix-loop-helix)结构域组成[5]。目前对PIF家族转录因子的研究主要集中在模式植物拟南芥Arabidopsis thaliana中。PIF转录因子广泛参与了拟南芥的光调控反应,并在光敏色素介导的光信号网络调控中处于主导地位[6]。黑暗条件下,AtPIF1能够抑制拟南芥种子的萌发,AtPIF3、AtPIF4、AtPIF5能够促进叶绿素分解并抑制叶绿体活性,从而最终导致叶片衰老[7]。遮阴条件下,AtPIF1能够通过调控AtHB1的表达来促进下胚轴的伸长[8]AtPIF4、AtPIF5和AtPIF7通过促进AtFTAtTSF的表达从而促进拟南芥开花[9]AtPIF4、AtBZR1和AtARF6通过整合生长素信号与BR和光敏色素通路来协同促进拟南芥下胚轴的伸长[10]。由此可见,PIF家族转录因子是遮阴条件下光信号转导途径中的关键调控因子。据此推测,PIF家族转录因子在光照强度对木兰属濒危树种生长发育的影响及光受体信号转导途径中起到重要作用。因此,本研究对景宁木兰PIF家族转录因子进行分析鉴定,探究其在极端遮阴条件下的表达模式,从而为揭示光介导的景宁木兰生长发育的分子机制、群落引起的遮阴而导致景宁木兰濒危的内在机制奠定基础,最终为珍稀植物景宁木兰的保育提供理论依据。

    • 景宁木兰来自浙江省杭州市青山湖花园中心苗圃培育的3年生扦插苗。2019年6月下旬,选择长势一致且无病虫害的扦插苗作为材料,置于温度为(25±2) ℃、光照强度为2 000~2 500 lx、相对湿度为40%~60%的人工气候室中进行培养,光照与黑暗时间各12 h,定期浇水保持土壤湿润。

    • 3层遮阳网形成的光照强度是景宁木兰自然生境光照条件之一,会显著削弱景宁木兰的生长势[1]。2019年7月上旬开始遮阴处理,本研究于浙江省杭州市临安区浙江农林大学人工气候室内,用3层黑色遮阳网及竹竿搭建极端遮阴处理(25%光照,ST)装置,以模拟遮阴条件下的景宁木兰野外生存环境。以100%全光照为对照(ck)。所有盆栽苗叶片之间无重叠。每个处理4个生物学重复。分别自实验开始处理第0、1、3、5、10、20、30天时,拍照记录植株表型变化、采集景宁木兰植株的中层叶片装入锡纸袋中,并迅速置于液氮中速冻,随后转入−80 ℃冰箱中储存,用于植物RNA的提取及后续实时荧光定量聚合酶链式反应(qRT-PCR)实验。

    • 基于实验室已有的转录组测序数据,在NR、NT、Swiss-Prot这3个数据库中进行注释,初步注释到46个PIF转录因子相关基因,将这些基因编码蛋白质逐个进行NCBI Blast和NCBI CDD预测,去除重复序列及冗余转录本后,最终得到9个PIF家族转录因子的蛋白质序列。利用MEGA 5.0对景宁木兰PIF转录因子序列与拟南芥PIF转录因子序列进行比对分析,根据TOLEDO等[5]命名拟南芥PIF家族的方法对9个景宁木兰PIF转录因子进行命名。拟南芥PIF转录因子序列从拟南芥转录因子数据库(https://www.arabidopsis.org/)中获得。

    • 采用Compute pI/Mw在线工具预测景宁木兰PIF转录因子蛋白质的等电点及分子量,采用在线网站Cell-PLoc 2.0预测亚细胞定位,利用NetPhos分析蛋白质磷酸化位点;利用MEME分析蛋白质保守结构域;利用DNAMAN比对分析保守结构域;利用在线软件heatmapper(http://www2.heatmapper.ca/expression/)绘制热图;为研究PIF转录因子家族的进化关系,利用MEGA5.0将9个景宁木兰MsPIF蛋白质序列、15个拟南芥AtPIF蛋白质序列[11]、6个玉米Zea mays ZmPIF蛋白质序列[12]、10个杨树Populus trichocar PtPIF蛋白质序列(http://www.Phytozome.net/poplar.php)[13]进行序列比对,并用邻接法构建系统进化树,进行Bootstrap测试,重复设置为1 000 次。

    • 采用植物RNA提取试剂盒提取样品总RNA。取1 μg RNA进行第1链 cDNA的合成。所用定量PCR试剂为BCG qPCR Master Mix。根据景宁木兰转录组数据设计特异性引物(表1)。以景宁木兰的EF1-α基因作为内参基因[1],利用Light Cycler 480II(Roche)仪器进行目的基因qRT-PCR表达分析。反应程序结束后,65~95 ℃每隔0.2 s作溶解曲线,采用${2^{ - \Delta \Delta {C_{\rm{t}}}}} $法计算目的基因的相对表达量。每个样品重复3次。

      表 1  实时荧光定量PCR引物

      Table 1.  Primers used for qRT-PCR analysis

      基因名称正向(反向)引物序列(5'→3')退火温度/℃扩增效率/%相关系数R2
      MsbHLH68F-CCTCAGGCTCTGTCATTGGG59.590.860.998
      R-AGCGGTTGGGCTTCTTCACG59.5
      MsPIF4F-GCCGAAGAATCTGAGTGCCA57.496.980.999
      R-CTATCCCTCCGTCTCCTTTC57.4
      MsbHLH1F-CTTTGTGCTATCTTCGGGAACG57.7104.480.997
      R-GGGACCCACCAATCAACGAC59.5
      MsPIF8F-AAGCGAGGATGAGAACAAGG55.493.010.983
      R-TGCGTTCGGATTGGTTATGG55.4
      MsbHLH66F-AGCAGTAACGGCACGCAGAC59.5101.570.995
      R-GAAATGGGCATGAGGCAGAG57.4
      MsPIF3F-CAGACTCAGCCGTCAACTCA57.492.930.991
      R-GCAGGCCCACTTCCACCAAT59.5
      MsPIF1F-CCACATATTCCTCCATTTGAT51.7100.900.994
      R-GAAGGCTGAGGTAGTGTTGATA55.8
      MsbHLH48F-CCGTCGAGTCTCCAGTGGTC61.696.920.993
      R-CGTTTGAAGGGAAGGTAGGG57.4
      MsbHLH23F-GCAGACAAAGGACCCAAGGA57.493.010.993
      R-CACGAGGAGTGAGGAACAAGAA57.7
    • 采用SPSS 19.0进行方差分析并用SigmaPlot 14绘图,对平均值采用独立样本T检验进行比较,显著水平为0.05。

    • 利用ExPASy对9个基因编码的氨基酸进行理化性质进行分析(表2),发现9个基因编码的蛋白质大小不等,最小的为MsPIF1,编码188个氨基酸,最大的为MsPIF3,编码735个氨基酸。9个蛋白质的理论等电点为5.18(MsbHLH68)~8.22(MsPIF8),表明不同的PIF家族蛋白质在不同的微环境中发挥不同的生物学功能。9个PIF家族蛋白质均属于不稳定、亲水性蛋白质。

      表 2  MsPIF基因家族

      Table 2.  MsPIF gene family of M. sinostellata

      基因名称蛋白质分子量/Da氨基酸数/个等电点不稳定系数脂肪系数平均亲水性亚细胞定位磷酸化位点
      SerThrTyr
      MsbHLH6853 604.795065.1858.1977.51−0.447细胞核(N)47101
      MsPIF457 883.065275.9656.0564.57−0.611细胞核(N)40123
      MsbHLH174 557.306715.5253.0879.66−0.452细胞核(N)34146
      MsPIF850 048.294568.2242.0958.29−0.686细胞核(N)36204
      MsbHLH6645 976.464306.7960.4569.72−0.558细胞核(N)45141
      MsPIF378 957.027356.3167.3763.80−0.513细胞核(N)74133
      MsPIF120 314.561885.9360.7860.32−0.341细胞核(N) 4 21
      MsbHLH4840 255.003666.4359.8570.90−0.639细胞核(N)40 62
      MsbHLH2339 480.073625.8272.0770.14−0.577细胞核(N)43 92
        说明:Ser为丝氨酸,Thr为苏氨酸,Tyr为酪氨酸
    • 通过在线软件Cell-PLoc 2.0(http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/)对9个PIF基因家族成员的氨基酸序列进行亚细胞定位预测(表2),9个基因编码的蛋白质均定位于细胞核内,符合其转录因子的功能,可能通过调控细胞核基因的转录来发挥功能。

    • 蛋白质磷酸化在很多植物信号转导过程中起着重要作用[14]。蛋白质磷酸化的可逆性使得植物细胞能够快速响应外部环境变化[15]。在正常光条件下,PIF家族转录因子与光敏色素相互作用,从而导致PIF蛋白质的磷酸化和降解[16]。而在低光照条件下,PIF蛋白质的磷酸化和降解受到抑制[17]。因此,分析PIF蛋白质的磷酸化位点对研究其在信号转导中的功能有重要意义。利用NetPhos对9个PIF家族蛋白质的磷酸化位点进行分析(表2),发现9个蛋白质均具有丝氨酸(Ser)、苏氨酸(Thr)和酪氨酸(Tyr)磷酸化位点,且Ser的磷酸化位点最多,Tyr的磷酸化位点最少。其中磷酸化位点最多的是MsPIF3,为94个,MsPIF1的磷酸化位点最少,仅7个。

    • 进一步利用MEME在线软件分析MsPIF蛋白质[16],共鉴定到6个基序(表3)。结果显示(图1):9个蛋白质含有不同数目的基序,MsPIF4和MsPIF8含有的基序最多,为5个,而MsbHLH48中含有的基序最少,仅2个。bHLH家族蛋白质均含有bHLH基序,该基序由位于N端的basic基序(图2A)和位于C端的HLH基序(图2B)组成[5]。9个PIF家族蛋白质均含有basic基序及HLH基序。不同蛋白质的同源异型基序之间存在一定差异。MsPIF4、MsPIF3和MsPIF8中含有APB基序(active phyB-binding motif)(图2C),表明其可以与PhyB(phytochrome B)相结合,从而对下游基因的表达进行调控[16]。此外,MsPIF4和MsPIF3中还含有APA基序(active phyA-binding motif)(图2D),表明其既可以与PhyA(phytochrome A)相结合也可以与PhyB相结合。

      表 3  景宁木兰PIF家族蛋白质的主要MEME基序

      Table 3.  Major MEME motif sequences in M. sinostellata PIF protein

      基序宽度/个最佳匹配
      Motif140MKALQELIPNSNKTDKASMLDEAIEYLKSLQLQVQMLSM
      Motif221RGRAAEVHNLSERRRRERINE
      Motif321MMFPGVQQYMPPMGMGMGMGM
      Motif436FEQZIVKLMEEDMGSAMQYLQGKGLCLMPIALAAAI
      Motif5 8ELLWENGQ
      Motif6 8CVPBWBTD

      图  1  MsPIF蛋白质的基序

      Figure 1.  Motif of MsPIF protein

      图  2  MsPIF家族蛋白质的basic、HLH、APB和APA基序蛋白质序列比对

      Figure 2.  Sequence alignment of basic, HLH, APB and APA motif protein sequences of MsPIF

    • 运用MEGA5对所得景宁木兰的蛋白质序列与拟南芥、玉米、杨树中的序列对比分析(图3),采用邻接法及Bootstrap分析(1 000次重复)构建系统进化树[18]。结果表明:植物PIF家族蛋白质表现出多样性,可以分为7个进化分枝。其中,9个MsPIF蛋白质分布在α、β、γ、δ、η、ε等6个进化枝。15个AtPIF蛋白质分布在α、β、γ、δ、ζ、ε等6个进化枝,其中ζ为拟南芥所特有。ZmPIF蛋白质仅分布于β和δ进化枝,10个PtPIF蛋白质分布在α、β、γ、δ、ε等5个进化枝。从进化树可以看出,MsPIF8、MsPIF4、MsPIF1、MsPIF3及MsbHLH23均与相应的拟南芥蛋白质与杨树蛋白质聚类在一起,表明其亲缘关系较为接近。MsbHLH68、MsbHLH66、MsbHLH48汇聚于η枝,为景宁木兰所特有。

      图  3  景宁木兰、拟南芥、玉米和杨树的PIF家族蛋白质的系统进化分析

      Figure 3.  Phylogenetic tree analysis of MsPIF protein and AtPIF protein sequences in M. sinostellata, A. thaliana, Z. mays, P. trichocarpa

    • 光照是植物生长发育过程中不可或缺的因素,极端遮阴处理过程中,景宁木兰叶片表型变化十分明显(图4)。处理10 d,景宁木兰叶片出现褐色斑点;处理20 d,褐色斑点面积增大,植株开始萎蔫;处理30 d,大量叶片干枯脱落,植株萎蔫严重。

      图  4  遮阴处理过程中景宁木兰的表型变化

      Figure 4.  Phenotypic changes of M. sinostellata during shading treatment

      PCR扩增琼脂糖凝胶电泳检测发现9个基因均扩出预期条带(图5)。qRT-PCR检测结果(图6)表明:遮阴处理1 d时,MsbHLH23和MsbHLH66的表达水平显著上调,MsbHLH1、MsbHLH48、MsPIF8、MsPIF3的表达水平均显著下调,其余基因表达水平无显著变化。遮阴处理3 d时,MsPIF8的表达显著低于对照,其余基因的表达均显著上调。遮阴处理5 d时,MsPIF8的表达仍显著降低,MsbHLH68、MsbHLH23、MsbHLH66、MsbHLH1、MsbHLH48的表达显著上调,其余基因的表达无显著差异。遮阴处理10 d时,MsPIF8、MsPIF7、MsPIF3的表达均显著下调,MsbHLH68、MsbHLH23、MsbHLH66、MsPIF4的表达显著上调。遮阴处理20 d时,MsPIF8、MsbHLH23、MsPIF1的表达相对于对照显著下调,MsbHLH68和MsbHLH66的表达显著上调,其余基因无显著差异。遮阴处理30 d时,MsPIF8、MsbHLH23的表达相对于对照显著下调,MsbHLH68、MsbHLH1、MsbHLH66、MsbHLH48、MsPIF3的表达显著上调。其中,MsbHLH23的表达变化比其他基因更为明显,遮阴处理5和10 d时的表达量分别上调为对照的52.77与20.03倍。

      图  5  qRT-PCR扩增产物特异性

      Figure 5.  Specificity of qRT-PCR amplification products

      图  6  遮阴条件下景宁木兰PIF转录因子的qRT-PCR分析

      Figure 6.  qRT-PCR analysis of PIF transcription factor under shading treatment of M. sinostellata

      基因的表达模式与基因功能关系密切,基于获得的景宁木兰处理组的qRT-PCR数据,采用heatmap在线软件对这些基因的差异表达进行分层聚类分析(图7)。结果发现:MsPIF1、MsPIF3、MsPIF4、MsPIF8的表达模式十分相似,其表达均在处理3 d时明显升高,随后逐渐下降。MsbHLH48的表达在处理0~30 d逐渐上调,其中处理30 d时的表达上调最为明显。MsbHLH1、MsbHLH68、MsbHLH66、MsbHLH23的表达模式均为先逐渐上调后逐渐下调,其中MsbHLH1和MsbHLH23的表达在处理5 d显著上调,MsbHLH68和MsbHLH66的表达分别在处理10和20 d明显上调。

      图  7  长期极端遮阴条件下景宁木兰PIF转录因子的表达模式

      Figure 7.  Pattern of expression of PIF family during shading treatment period of M. sinostellata

    • PIF家族转录因子广泛存在于植物中。目前PIF家族转录因子在拟南芥中研究的最为透彻,共有15个[5]。在其他物种中,杨树至少存在10个[13],玉米有6个[12],水稻Oryza sativa有6个[19]。本研究从景宁木兰中仅筛选到9个,这一方面是由于景宁木兰本身具有特异性,另一方面可能是因为转录组数据库还不够丰富造成的。NCBI blast及NCBI CDD预测结果表明:9个PIF蛋白质均含bHLH结构域,属于PIF转录因子家族。同源性较高的蛋白质通常功能也相似[18]。因此,本研究用所得9个景宁木兰PIF家族转录因子编码蛋白质与6个ZmPIFs、10个PtPIFs 和15个功能明确的AtPIFs编码蛋白质作系统进化树,主要根据与15个拟南芥PIF家族转录因子编码蛋白质进化关系的远近推测景宁木兰PIF转录因子的功能。AtPIF3与PtPIF3a、PtPIF3b同源性较高,且与AtPIF3和ZmPIF3位于同一进化枝,据此推测MsPIF3与AtPIF3功能相似。MsPIF4与PtPIF4、ZmPIF4、AtPIF4位于同一进化枝。AtPIF3蛋白质在植株刚刚暴露到光下时会急剧减少,在黑暗条件下会重新积累,在拟南芥幼苗阶段,AtPIF3在光诱导的叶绿体发育中起到积极的作用[20]。AtPIF3、AtPIF4和AtPIF5在年龄引起的和黑暗引起的植物衰老中起到重要作用,其中AtPIF4通过调控叶绿素降解和维持叶绿体活性的相关基因从而促进植株衰老[7]。AtPIF4能够调控下游基因从而促进拟南芥开花[21]

      此外,AtPIF4能与AtBZR1共同调节下游光响应的相关基因。AtPIF7与MsbHLH1同源性较高,推测其功能可能类似。遮阴条件下,AtPIF7是调控拟南芥茎伸长的重要因子[17]。AtPIF7是调节植物去黄化的弱负反馈因子,与PhyB结合的过程中通常不伴随着可以监测到的磷酸化和分解[22]。MsPIF8与AtPIF8同源,推测其功能相似。与AtPIF3不同,AtPIF8蛋白质在远红光下的积累比在黑暗条件中或红光条件下更多,AtPIF8能够抑制由PhyA介导的种子萌发、下胚轴伸长等生命活动[23]。MsbHLH23与AtbHLH24功能相似,AtbHLH24在抑制种子萌发和维持种子休眠状态中具有重要作用[24]。AtbHLH24能够激活细胞分裂素信号,同时激活生长素生物合成和雄蕊内部结构域的转运基因,并且在花瓣和果实发育中起到重要作用[25]。MsbHLH68,MsbHLH66和MsbHLH48共同汇集于同一进化枝,其功能还需要进一步挖掘。但是,即使是同一分支内蛋白质功能也可能存在差异,并且不同物种之间具有较大的差异[18]。因而关于MsPIF家族蛋白质的具体功能还需要通过后续研究。

    • 群落层片造成的遮阴会降低下层植被的生长率与成活率[26]。本研究极端遮阴处理组的景宁木兰在处理过程中逐步萎蔫死亡。光敏色素作用因子(PIFs)在植物的生长发育中起到重要调节作用。PIFs作为光信号转导途径中的负调控因子,对植物生长发育的多个过程都有调控作用[27-32]。PIF家族转录因子是光敏色素介导的光信号转导途径中的负反馈因子[28]。遮阴条件下,PIF蛋白质大量积累,从而调控植物发生一系列生理生化反应。考虑到植物的个体差异和基因表达差异,本研究的每个取样时间点都作了相应的空白对照。本研究中,长期极端遮阴条件下,景宁木兰的PIF类家族转录因子的表达水平多数显著上调,只有MsPIF8的表达水平显著下调。在拟南芥中,与其他PIF家族成员相比,AtPIF8能够与PhyB结合但效果微弱[23]MsPIF8的表达模式与前人的研究结论存在差异,可能与物种特异性相关,有关MsPIF8的功能还需要进一步研究。MsPIF1、MsPIF3和MsPIF4的表达模式与前人的研究结论[29]吻合。黑暗条件下,AtPIF1的表达上调,从而抑制叶绿素的合成与种子萌发;AtPIF3、AtPIF4的表达显著上调,从而促进叶绿素的降解与叶片衰老[29]。植物群落造成的遮阴使得AtPIF4的表达上调,从而使成花诱导相关基因的表达上调,最终诱导拟南芥提早开花[21]。PIF功能十分复杂,景宁木兰PIF家族转录因子在遮阴调控网络中发挥的作用还需要进一步研究。

参考文献 (32)

目录

/

返回文章
返回