-
光是影响植物生长发育的最重要的环境因子之一。为了感知生长环境中的光信号并对各种光信号的变化及时做出响应,植物进化出了不同的光受体,通过光受体的调节使环境因子达到最适宜植物生长发育的状态[1]。先前的研究表明,植物共进化出四种光受体,分别包括感受红光和远红光的光敏色素(Phy)[2],感受蓝光的隐花色素(Cry)和向光蛋白(Phot)[3-4]以及感受紫外光-蓝光的UV-B受体(UV-B receptor)[4-5], 其中,感受红光和远红光的光受体-光敏色素是目前研究较清楚的光受体,参与了植物从种子萌发到成熟的整个生长发育过程[6-7]。光敏色素是一种可溶性色素蛋白,在植物体内存在两种状态,分别是具有生物活性的远红光吸收型(FR-absorbing form, Pfr)和不具有生物活性的红光吸收型(R-absorbing form, Pr),且这两种光敏色素可以在红光和远红光的照射下进行相互转换[8]。
在拟南芥基因组中共发现了5个光敏色素基因(PhyA-E)[9-10], 根据其在光下的稳定性不同,拟南芥光敏色素可以分为Ⅰ型光敏色素和Ⅱ型光敏色素。Ⅰ型光敏色素主要存在于黄化幼苗中,遇光极易降解,PhyA属于Ⅰ型光敏色素;Ⅱ型光敏色素主要存在于绿色组织中,遇光稳定,PhyB-E属于Ⅱ型光敏色素[8, 11]。在暗生长的黄化幼苗中,光敏色素以不具有生物活性的状态定位于细胞质中;感受到红光后,光敏色素转变为具有生物活性的远红光吸收型定位于细胞核中,并与细胞核内的其它蛋白质发生相互作用[12-14]。研究表明,光敏色素相互作用因子(PIF)能够直接与具有生物活性的光敏色素结合而启始光信号的传导[15-16]。
PIFs属于bHLH(basic helix-loop-helix)转录因子家族的一个亚家族,也称作PILs(phytochrome interacting factor-like),能够与具有生物活性型的光敏色素相互作用[17]。研究表明,大部分PIFs蛋白质在N端具有保守的APB结构域(active phytochrome-binding domain),APB结构域是PIFs与具有生物活性型光敏色素相结合的不可或缺的元件[18]。拟南芥中至少有7个PIFs成员,PIF1、PIF3、PIF4、PIF5、PIF6、PIF7和PIF8[18-24], 其中,PIF3是利用PHYB作为诱饵通过酵母双杂交筛选出来的最早的PIFs家族成员[20]。免疫共沉淀结果表明,PIF3能与具有生物活性型的PhyA和PhyB的C端相结合。拟南芥中只有PIF1和PIF3能够同时与PhyA和PhyB相互作用,而PIF4-PIF8只能与PhyB相互作用[25]。拟南芥PIF基因家族成员在下胚轴伸长和幼苗避荫反应等过程发挥了重要作用[26-29]。PIF3参与黄花幼苗的形成,在幼苗由黑暗至光照和由光照至黑暗过程中起着重要作用[30]。PIF1参与光诱导的种子萌发,抑制下胚轴伸长,负调控叶绿素合成[21]。PIF4正调控气孔的发育,PIF4、PIF5和PIF7受光周期的调控而影响下胚轴延伸[31]。此外,PIF类的转录因子还具有与靶基因启动子区顺式调控元件G-box (CACGTG)相结合的DNA结合结构域[21]。这些研究结果表明,PIF类转录因子在调控种子萌发、幼苗脱黄化以及开花时间等发育过程具有重要作用。
杨树是多年生木本植物,光是调控其生长发育的重要环境信号之一,光信号的变化使杨树需要具备适应不同发育状态的能力。PIF作为光敏色素相互作用因子,对光环境的变化做出相应的改变调整其达到最适宜树木生长的状态;但目前对于杨树的PtPIF基因家族成员在基因组中分布及其表达模式却鲜见报道。为此,本研究利用生物信息学对杨树PtPIF基因家族成员进行了系统进化分析,鉴定其基因结构、保守结构域,分析了不同家族成员的组织表达特性及其在不同非生物胁迫下的表达规律,研究结果为以后深入解析杨树PtPIF基因的功能奠定了基础。
HTML
-
实验所用材料毛果杨(Populus trichocarpa Torr. & Gray)由中国林业科学研究院林木遗传育种国家重点实验室保存。
-
RNA提取试剂盒RNeasy Plant;Mini与RNase-free DNase Ⅰ均购买于Qiagen公司;反转录试剂盒SuperscriptⅢ First-strand Kit以及PCR所用试剂均购买于Takara公司。采用Primer 5.0软件设计引物,实验所需引物在Invitrogen公司合成。
-
首先获取TAIR(http://www.arabidopsis.org/)数据库中已公布的拟南芥PIF基因及其编码的蛋白质序列,并用得到的拟南芥PIF蛋白序列作为诱饵,在杨树基因组数据库(http://www.phytozome.net/poplar.php)中通过blastp比对检索,获得PIF直系同源序列。使用MEGA5.0软件中的Neighbor-Joining(NJ)法构建系统发育进化树。杨树PtPIF蛋白的基本理化性质采用ProtParam(http://web.expasy.org/portparam/)分析完成。使用Gene Structure Display Server(http://gsds.cbi.pku.edu.cn/)软件分析杨树PtPIF基因的外显子及内含子结构。使用在线软件MEME(http://meme-suite.org/)分析杨树PtPIF蛋白的结构域。
-
使用RNA提取试剂盒(RNeasy Plant Kit)提取毛果杨的幼叶、成熟叶、初生茎、次生茎和根的总RNA,反转录合成cDNA第一链,并以cDNA为模板进行qRT-PCR。反应体系为:LightCycler®80 2×SYBR®Premix Ex TaqTM 10 μL,5′-引物(10 μmol·L-1)0.8 μL,3′-引物(10 μmol·L-1)0.8 μL,cDNA2 μL,超纯水6.4 μL;反应程序为:94℃30 s;94℃12 s,58℃30 s,72℃40 s,45个循环。引物序列见表 1。
引物
Primer引物序列(5′-3′)
Primer sequence(5′-3′)PtPIF1-F TACCTGGTACTAAAGAGATCCGATC PtPIF1-R GAGGATATAGCAAATCGGCG PtPIF3a-F CAAGTAGTAGCAATCCGGCTG PtPIF3a-R GGAATTGTAGCCTCAGGTAGCT PtPIF3b-F AATAATCTATGTGCAGGATGCTGA PtPIF3b-R CTCATTGATCCTATCCCTTCG PtPIF4/5a-F CATAGAGCCAGACAATGAGCTAGTA PtPIF4/5a-R CTTGAATCAAATGACTTGAATTCAG PtPIF4/5b-F TGGAAGGAAAGTCGGCTACA PtPIF4/5b-R TGAGGTTTATTGTTGTTGAGTTGTC PtPIF8a-F AATGGGAATGGGAATGGG PtPIF8a-R GGCATTACCGTGGCTGAA PtPIF8b-F CAGTACTGGTTCGGACCATGA PtPIF8b-R GCTGGGCCGGAGTCTATC PtPIF9a-F ATCAACGTCTTCTTCCTCATTTATA PtPIF9a-R GTCTCGTTCTCACATTCCCAG PtPIF9b-F AGCTCAGAGCAAAATTCTCTTAATC PtPIF9b-R CGTGGACAGATTGAAAGGTG PtPIF10-F TTCAATGACCAACAATATCTAATGG PtPIF10-R ACCTCTTACTGGTAGTGGCATAGTT PtActin-F GTCCTTCTAACTTCCCAACAGTGC PtActin-R GACTACCAAAGTGTCTGACCACCA Table 1. The primer pair sequences used for quantitative real-time PCR
-
毛果杨组培苗培养1个月后,转移至人工气候箱中用Hogaland营养液水培。水培15 d后,分别进行4℃低温处理、干旱处理(10% PEG6000)及盐处理(150 mmol·L-1 NaCl),于处理6、12、24、36 h后取由顶端数的第4片叶片,每个处理设置3个重复,每次重复3株,以未处理的毛果杨作为对照。
-
使用RNeasy Plant Mini试剂盒与RNase-free DNase Ⅰ试剂盒提取处理后杨树叶片的RNA,使用反转录试剂盒SuperscriptⅢ First-strand Kit合成cDNA第一链。使用Primer 5.0在杨树PIF基因核苷酸序列的非保守区设计定量PCR扩增引物(表 1),引物序列的GC含量低于30%,Tm值为58~60℃,并由Invitrogen公司合成。实时荧光定量PCR根据SYBR®Premix Ex TaqTMⅡkit(TaKaRa Dalian, China)和Roche480实时荧光定量PCR仪的使用说明进行操作。反应体系为:LightCycler®480 2×SYBR®Premix Ex TaqTM 10 μL,上下游引物各0.8 μL,cDNA 2 μL,ddH2O 6.4 μL。选取杨树Actin(基因序列号:EF418792.1)作为内参基因。每个cDNA模板设置四个重复,所有实验重复三次,采用2-ΔΔCt计算方法进行数据分析。
1.1. 实验材料
1.2. 实验试剂
1.3. 杨树PIF基因序列的获取与分析
1.4. 杨树PIF基因组织特异性表达分析
1.5. 材料处理
1.6. 实时荧光定量PCR分析
-
通过序列比对分析,在杨树基因组数据库中共筛选得到10个PtPIF基因成员(表 2)。这10个PIF基因长度为2 810~6 064 bp,所编码的蛋白序列长度为310~742 aa,蛋白质分子质量介于33~79 kDa,等电点为5.42~9.17,既有酸性蛋白,又有碱性蛋白。使用NJ法构建基于杨树和拟南芥PIF成员序列的系统发育进化树(图 1),根据聚类结果,将杨树10个PIF家族成员基因分别命名为PtPIF1、PtPIF3a、PtPIF3b、PtPIF4/5a、PtPIF4/5b、PtPIF8a、PtPIF8b、PtPIF9a、PtPIF9b和PtPIF10。
基因名
Gene基因号
Gene number基因在染色体上的位置
Location基因长度
Length of gene/bp蛋白质长度
Length of protein/aa等电点/分子量
pI/(MW/kDa)PtPIF1 Potri.002G252800 Chr02(-) 24187022-24193085 6 064 615 7.29/(67.366 00) PtPIF3a Potri.005G001800 Chr05 (-) 107969-111992 4 024 742 5.89/(79.784 41) PtPIF3b Potri.013G001300 Chr13 (-) 86122-91428 5 307 719 5.61/(77.620 79) PtPIF4/5a Potri.002G055400 Chr02 (-) 3705303-3709545 4 243 583 6.13/(64.164 79) PtPIF4/5b Potri.005G207200 Chr05 (+) 22153104-22156744 3 641 478 9.17/(52.311 36) PtPIF8a Potri.002G143300 Chr02 (+) 10651063-10655266 4 204 484 8.80/(52.715 58) PtPIF8b Potri.014G066500 Chr14 (-) 5402159-5406848 4 690 471 8.43/(51.139 76) PtPIF9a Potri.005G139700 Chr05 (+) 11648940-11653997 5 058 327 6.13/(35.941 34) PtPIF9b Potri.014G025800 Chr14 (+) 2209016-2211913 2 898 310 5.42/(33.741 64) PtPIF10 Potri.014G111400 Chr14 (-) 8714454-8717263 2 810 561 6.97/(61.476 18) Table 2. Basic information of PtPIF family members in Populus
通过对杨树PIF家族各成员的基因结构(外显子/内含子)分析,发现其中PtPIF10由5个外显子和4个内含子组成,PtPIF4/5b、PtPIF8a、PtPIF8b和PtPIF9b均由6个外显子和5个内含子组成,PtPIF3a、PtPIF3b和PtPIF9a均由7个外显子和6个内含子组成,PtPIF4/5a由8个外显子和7个内含子组成,而PtPIF1由10个外显子和9个内含子组成(图 2)。
-
使用软件对杨树PIF蛋白进行保守结构域分析。结果表明:杨树PIF蛋白与拟南芥PIF蛋白在功能结构域上具有较高的保守性(图 3),都具有bHLH DNA结合结构域,推测它们可能具有相似的功能。
-
为了探究杨树中PIF基因成员的组织表达特性,利用qRT-PCR方法分析了PtPIF在幼叶、成熟叶、初生茎、次生茎及根等不同组织中的表达模式, 结果(图 4)表明:PtPIF1、PtPIF8a与PtPIF9b在各组织中普遍表达,且PtPIF1与PtPIF8a均在成熟叶的表达相对较高,而PtPIF9b在幼叶的表达相对较高。PtPIF3a与PtPIF8b在根中检测不到表达,但在其他部位均有表达,且在幼叶和成熟叶中高表达。PtPIF4/5b只在成熟叶中表达,在其它部位检测不到表达。PtPIF3b、PtPIF4/5a、PtPIF9a和PtPIF10在各组织中的表达水平较低。杨树PIF基因家族成员在各组织中的表达差异,说明不同成员在杨树发育过程中可能具有不同的功能。
-
为揭示杨树PIF基因在非生物胁迫下的响应机制,利用qRT-PCR方法检测其在低温、干旱和NaCl胁迫下的表达模式(图 5)。在盐(150 mmol ·L-1 NaCl)、干旱(10% PEG)及4℃低温处理下,杨树PtPIF基因家族各成员间在响应时间、程度上各有不同。低温处理后,PtPIF1、PtPIF3a、PtPIF4/5a、PtPIF4/5b、PtPIF9b和PtPIF10响应时间较早,表达量均在6 h较对照显著改变,其中, PtPIF9b和PtPIF10分别是对照的0.3和4.3倍。之后PtPIF9b的表达量基本维持这个水平,而PtPIF10的表达量则呈先上升后下降的趋势,在24 h达到对照的8.2倍,36 h仅为对照的2.3倍。PtPIF1、PtPIF4/5a和PtPIF8b的表达量均随着处理时间的延伸而降低,且在36 h达到最低值,分别是对照的0.12、0.04、0.23倍。PtPIF3b和PtPIF9a的表达量呈现先下降后上升的趋势,在36 h均到达对照的2.5倍。干旱处理下,PtPIF1、PtPIF3a、PtPIF3b、PtPIF4/5a、PtPIF4/5b、PtPIF9b和PtPIF10的表达量较对照组均下调,PtPIF8a和PtPIF9a的表达量上调。盐处理下,PtPIF3a、PtPIF3b、PtPIF8b、PtPIF9a和PtPIF10的表达量均上调,PtPIF3b、PtPIF9a和PtPIF10的表达量在36 h达到对照的4.3倍。PtPIF4/5a的表达量随着处理时间的延伸而降低,在36 h为对照的0.25倍。