-
景宁木兰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通过促进AtFT和AtTSF的表达从而促进拟南芥开花[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蛋白质序列(http://www.plantgdb.org/ZmGDB/)[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 MsbHLH68 F-CCTCAGGCTCTGTCATTGGG 59.5 90.86 0.998 R-AGCGGTTGGGCTTCTTCACG 59.5 MsPIF4 F-GCCGAAGAATCTGAGTGCCA 57.4 96.98 0.999 R-CTATCCCTCCGTCTCCTTTC 57.4 MsbHLH1 F-CTTTGTGCTATCTTCGGGAACG 57.7 104.48 0.997 R-GGGACCCACCAATCAACGAC 59.5 MsPIF8 F-AAGCGAGGATGAGAACAAGG 55.4 93.01 0.983 R-TGCGTTCGGATTGGTTATGG 55.4 MsbHLH66 F-AGCAGTAACGGCACGCAGAC 59.5 101.57 0.995 R-GAAATGGGCATGAGGCAGAG 57.4 MsPIF3 F-CAGACTCAGCCGTCAACTCA 57.4 92.93 0.991 R-GCAGGCCCACTTCCACCAAT 59.5 MsPIF1 F-CCACATATTCCTCCATTTGAT 51.7 100.90 0.994 R-GAAGGCTGAGGTAGTGTTGATA 55.8 MsbHLH48 F-CCGTCGAGTCTCCAGTGGTC 61.6 96.92 0.993 R-CGTTTGAAGGGAAGGTAGGG 57.4 MsbHLH23 F-GCAGACAAAGGACCCAAGGA 57.4 93.01 0.993 R-CACGAGGAGTGAGGAACAAGAA 57.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 氨基酸数/个 等电点 不稳定系数 脂肪系数 平均亲水性 亚细胞定位 磷酸化位点 Ser Thr Tyr MsbHLH68 53 604.79 506 5.18 58.19 77.51 −0.447 细胞核(N) 47 10 1 MsPIF4 57 883.06 527 5.96 56.05 64.57 −0.611 细胞核(N) 40 12 3 MsbHLH1 74 557.30 671 5.52 53.08 79.66 −0.452 细胞核(N) 34 14 6 MsPIF8 50 048.29 456 8.22 42.09 58.29 −0.686 细胞核(N) 36 20 4 MsbHLH66 45 976.46 430 6.79 60.45 69.72 −0.558 细胞核(N) 45 14 1 MsPIF3 78 957.02 735 6.31 67.37 63.80 −0.513 细胞核(N) 74 13 3 MsPIF1 20 314.56 188 5.93 60.78 60.32 −0.341 细胞核(N) 4 2 1 MsbHLH48 40 255.00 366 6.43 59.85 70.90 −0.639 细胞核(N) 40 6 2 MsbHLH23 39 480.07 362 5.82 72.07 70.14 −0.577 细胞核(N) 43 9 2 说明：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
基序 宽度/个 最佳匹配 Motif1 40 MKALQELIPNSNKTDKASMLDEAIEYLKSLQLQVQMLSM Motif2 21 RGRAAEVHNLSERRRRERINE Motif3 21 MMFPGVQQYMPPMGMGMGMGM Motif4 36 FEQZIVKLMEEDMGSAMQYLQGKGLCLMPIALAAAI Motif5 8 ELLWENGQ Motif6 8 CVPBWBTD 
图 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家族转录因子在遮阴调控网络中发挥的作用还需要进一步研究。
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 Abstract:Objective The shading caused by the community is one of the important factors that lead to the endangerment of Magnolia sinostellata. Therefore, it is important to conduct a systematic analysis and research of PIF family transcription factors which play an important role in light signal transduction and plant growth. Also, such analysis will help lay a foundation for the exploration of its role in the light signal transduction mechanism of M. sinostellata. Method With the transcriptome data of M. sinostellata collected, transcription factors of PIF family were identified and analyzed by bioinformatic while the expression patterns were analyzed employing the qRT-PCR technology. Result The nine MsPIFs transcription factor genes screened from the M. sinostellata enjoyed a length of 188−735 aa, a protein size of 20314.56−78957.02 Da, and a theoretical isoelectric point range of 5.18−8.22. The proteins encoded by MsPIFs gene were unstable proteins, and all proteins were hydrophilic proteins localized in the nucleus as was demonstrated in the subcellular localization prediction. All nine proteins have Ser, Thr and Try phosphorylation sites. As was shown in the qRT-PCR results, under extreme shading conditions, there was the occurrence of changes of different degrees in the gene expression of 9 MsPIFs families of which MSBHLH23 has an expression level that was 52.77 and 20.03 times higher than that of the control at 5 d and 10 d after shading treatment. Conclusion The PIF transcription factor family of M. sinostellata can respond to shading and this study has laid a foundation for the identification of MsPIFs biological function. [Ch, 7 fig. 3 tab. 32 ref. ] -
图 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
图 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 MsbHLH68 F-CCTCAGGCTCTGTCATTGGG 59.5 90.86 0.998 R-AGCGGTTGGGCTTCTTCACG 59.5 MsPIF4 F-GCCGAAGAATCTGAGTGCCA 57.4 96.98 0.999 R-CTATCCCTCCGTCTCCTTTC 57.4 MsbHLH1 F-CTTTGTGCTATCTTCGGGAACG 57.7 104.48 0.997 R-GGGACCCACCAATCAACGAC 59.5 MsPIF8 F-AAGCGAGGATGAGAACAAGG 55.4 93.01 0.983 R-TGCGTTCGGATTGGTTATGG 55.4 MsbHLH66 F-AGCAGTAACGGCACGCAGAC 59.5 101.57 0.995 R-GAAATGGGCATGAGGCAGAG 57.4 MsPIF3 F-CAGACTCAGCCGTCAACTCA 57.4 92.93 0.991 R-GCAGGCCCACTTCCACCAAT 59.5 MsPIF1 F-CCACATATTCCTCCATTTGAT 51.7 100.90 0.994 R-GAAGGCTGAGGTAGTGTTGATA 55.8 MsbHLH48 F-CCGTCGAGTCTCCAGTGGTC 61.6 96.92 0.993 R-CGTTTGAAGGGAAGGTAGGG 57.4 MsbHLH23 F-GCAGACAAAGGACCCAAGGA 57.4 93.01 0.993 R-CACGAGGAGTGAGGAACAAGAA 57.7 
下载: 导出CSV
表 2 MsPIF基因家族
Table 2. MsPIF gene family of M. sinostellata
基因名称 蛋白质分子量/Da 氨基酸数/个 等电点 不稳定系数 脂肪系数 平均亲水性 亚细胞定位 磷酸化位点 Ser Thr Tyr MsbHLH68 53 604.79 506 5.18 58.19 77.51 −0.447 细胞核(N) 47 10 1 MsPIF4 57 883.06 527 5.96 56.05 64.57 −0.611 细胞核(N) 40 12 3 MsbHLH1 74 557.30 671 5.52 53.08 79.66 −0.452 细胞核(N) 34 14 6 MsPIF8 50 048.29 456 8.22 42.09 58.29 −0.686 细胞核(N) 36 20 4 MsbHLH66 45 976.46 430 6.79 60.45 69.72 −0.558 细胞核(N) 45 14 1 MsPIF3 78 957.02 735 6.31 67.37 63.80 −0.513 细胞核(N) 74 13 3 MsPIF1 20 314.56 188 5.93 60.78 60.32 −0.341 细胞核(N) 4 2 1 MsbHLH48 40 255.00 366 6.43 59.85 70.90 −0.639 细胞核(N) 40 6 2 MsbHLH23 39 480.07 362 5.82 72.07 70.14 −0.577 细胞核(N) 43 9 2 说明：Ser为丝氨酸，Thr为苏氨酸，Tyr为酪氨酸
下载: 导出CSV
表 3 景宁木兰PIF家族蛋白质的主要MEME基序
Table 3. Major MEME motif sequences in M. sinostellata PIF protein
基序 宽度/个 最佳匹配 Motif1 40 MKALQELIPNSNKTDKASMLDEAIEYLKSLQLQVQMLSM Motif2 21 RGRAAEVHNLSERRRRERINE Motif3 21 MMFPGVQQYMPPMGMGMGMGM Motif4 36 FEQZIVKLMEEDMGSAMQYLQGKGLCLMPIALAAAI Motif5 8 ELLWENGQ Motif6 8 CVPBWBTD
下载: 导出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. -
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20200488
点击查看大图
图(7) / 表(3)
计量
- 文章访问数: 1219
- HTML全文浏览量: 370
- PDF下载量: 75
- 被引次数: 0
