杉木木质素合成酶基因CCR的克隆和序列分析_陈喜
亚麻木质素合成关键酶CAD基因的克隆与分析的开题报告

亚麻木质素合成关键酶CAD基因的克隆与分析的开题报告1. 研究背景亚麻是一种经济作物,其纤维可以用于生产纺织品,种子可以提取出亚麻油。
亚麻木中含有大量木质素,木质素是一种复杂的天然高分子化合物,具有很强的稳定性和抗生物降解性,在工业、医药、环保等领域有重要的应用价值。
木质素的合成主要是由亚麻木细胞壁中存在的木质素合成关键酶催化完成的。
CAD是木质素合成的一个关键酶,为木质素单体的预芳香化的限速酶,其合成量和合成活性直接影响亚麻木质素的含量与质量。
因此研究CAD基因的克隆与分析对于了解亚麻木质素生物合成机制,提高亚麻木质素的含量和质量具有重要意义。
2. 研究目的本研究的主要目的是克隆亚麻木质素合成关键酶CAD基因,分析其结构和功能特点,进一步了解亚麻木质素生物合成的机制。
3. 研究内容和方法3.1 克隆CAD基因从亚麻木材中提取总RNA,通过RT-PCR方法克隆CAD基因,利用PCR扩增方法连接载体并进行重组,将重组载体转化到大肠杆菌中进行扩增。
3.2 分析CAD基因结构和功能特点通过生物信息学方法对CAD基因进行序列分析,包括ORF长度、氨基酸序列、保守区域等。
同时进一步分析CAD的保守结构域、亲水性、疏水性等结构特征。
对CAD基因进行功能现象研究,如酶活检测等实验。
4. 研究意义通过本研究,能够更深入地了解亚麻木质素合成的生物合成机制,探究CAD基因的结构和功能特点,为亚麻木质素的生产和发展提供理论根据,同时在木质素生物制造方面具有积极的应用价值。
5. 研究进展和计划目前已完成了从亚麻木材中提取总RNA和克隆CAD基因的实验工作。
接下来的工作将会进行CAD基因序列分析及其功能实验,进一步深入了解CAD基因的结构和功能特点。
最后,将进一步研究亚麻木质素合成的生物合成机制。
预计本研究将在2年内完成。
一个新杉木纤维素合酶ClCesA基因的克隆与植物表达载体构建

一个新杉木纤维素合酶ClCesA基因的克隆与植物表达载体构建魏晓骁;王士亚;陈潇潇;汪凤林;曹光球;曹世江【摘要】In order to explore the molecular mechanism of Chinese fir wood formation, total RNA was extracted from Chinese fir new leaves and cDNA was synthesized. A new ClCesA was cloned using RT-PCR. The length of open reading frame of ClCesA is 2955 bp, with 984 amino acids. Based on TMHMM Server v. 2.0 and MEGA5.1 software, ClCesA proteins consist of 8 across-membrane structures, with average 18 amino acid residues and zinc finger in N end. Also, a conservative DDDQXXLRW structure domain was found in ClCesA protein. Phylogenetic tree analysis showed that ClCesA may be associated with secondary cell wall formation. The flu-orescent quantitative PCR analysis showed that ClCesA can be expressed in organsof root, stem and leaf in Chinese fir, with the highest expression in the root and the least in leaf. Lastly, a plant overexpression vector pGWB506-35s-ClCesA-GFP was successful-ly constructed from pGWB506 vector which contains GFP tag and hygromycin resistance gene by Gateway technology.%为探索杉木木材形成的分子机制,以杉木嫩叶总RNA反转录的cDNA为模板,运用RT-PCR技术克隆出一个新杉木纤维素合酶基因ClCesA.该基因的开放阅读框(ORF)为2955 bp,共编码984个氨基酸.通过TMHMM Server v.2.0和MEGA 5.1软件预测ClCesA蛋白的跨膜结构和功能,发现该蛋白N端含有锌指结构,共有8个跨膜结构,平均跨膜区域18个氨基酸残基,并具有保守的DDDQXXLRW结构域.系统进化树分析结果表明,ClCesA可能与细胞次生壁形成有关.荧光定量PCR分析表明,ClCseA基因在杉木的不同器官(根、茎、叶)中均有表达,但在杉木根中的表达含量最高,茎中次之,叶中最低.随后,为后续该基因功能验证提供基础,将ClCesA 克隆到含有GFP标签和潮霉素抗性的pGWB506载体,构建植物过表达载体pGWB506-35S-ClCesA-GFP.【期刊名称】《福建农林大学学报(自然科学版)》【年(卷),期】2017(046)001【总页数】7页(P66-72)【关键词】杉木;纤维素合酶;基因克隆;表达分析;载体构建【作者】魏晓骁;王士亚;陈潇潇;汪凤林;曹光球;曹世江【作者单位】福建农林大学林学院;国家林业局杉木工程技术研究中心,福建福州350002;福建农林大学林学院;国家林业局杉木工程技术研究中心,福建福州350002;福建农林大学林学院;国家林业局杉木工程技术研究中心,福建福州350002;福建农林大学林学院;国家林业局杉木工程技术研究中心,福建福州350002;福建农林大学林学院;国家林业局杉木工程技术研究中心,福建福州350002;福建农林大学林学院;国家林业局杉木工程技术研究中心,福建福州350002【正文语种】中文【中图分类】Q785;S791.27纤维素作为地球上最丰富的可再生资源,是木材次生木质部细胞壁的主要组成成分之一,约占木材干重的40%,是决定木材品质的重要因子之一.为此揭示植物纤维素生物合成机制对木材材质改良、木材定向培育都具有十分重要的意义.纤维素合酶(CesA)是纤维素生物合成的关键酶之一[1,2].研究表明,CesA基因通过在转录水平上对纤维素合成进行调控,对植物细胞壁发育起到重要作用[3-5].美国密西根大学在2000年初次从欧洲颤杨(Populus tremuloides)中克隆得到林木的纤维素合酶基因PtrCesA1[6],并且研究发现该基因主要在植物茎中表达,表达高峰期主要出现在次生壁形成期间.此后,在毛果杨(Populus trichocarpa)[7]、大青杨(Populus ussuriensis)[8]、松树(Pinus radiata)[9]、苎麻(Boehmeria nivea)[10]、桉树(Eucalyptus)[11]、马尾松(Pinus massoniana)[12]、杉木(Cunninghamia lanceolata)[13]、毛竹(Phyllostachys edulis)[14]等木本植物中都陆续克隆分离出CesA基因,其中毛果杨中鉴定出18个CesA基因[7],杂交杨(Populus tremula (L)×P.tremuloides)中鉴定出10个CesA基因,绿竹(Bambusa bambusaoldhami)中鉴定出8个CesA基因[14].这些基因共同构成了林木的纤维素合酶基因家族.对ClCesA基因功能的研究结果表明,PtrCesA1、PtrCesA2、PtrCesA3这3个基因可能组成同一个纤维素合成酶复合体,共同完成杨树次生壁的合成[15];而PtrCesA4、PtrCesA5、PtrCesA7这3个基因可能与初生壁的形成有关[16].在毛果杨中,PtCesA4、PtCesA5、PtCesA7、PtCesA8、PtCesA17和PtCesA18在木质部特异高表达[17].Song et al[18]通过研究杨树的纤维素合酶发现,杨树的木质部细胞膜上至少具有2种纤维素合酶复合体,复合体Ⅰ由CesA4、CesA7、CesA8、CesA17和CesA18组成,参与细胞次生壁的合成;复合体Ⅱ由CesA3、CesA10、CesA11、CesA13、CesA15和CesA16组成,参与细胞初生壁和次生壁的合成.杉木作为中国南方重要的速生用材树种之一,用途非常广泛.近几年,与杉木材性相关的分子机制研究取得一定进展[19,20],而有关杉木的纤维素生物合成分子机制的研究报道甚少.彭沙沙等[21]克隆了杉木纤维素合成酶类似蛋白ClCslD1;而庞景等[13]克隆了杉木纤维素合成酶基因CesA,通过将本研究克隆出的基因ClCseA与庞景等克隆的CesA基因的DNA和蛋白序列进行比较,发现2个基因在核苷酸的编码区有11处位点存在差异,分别为773、821、1 249、1 537、1 545、1 554、2 166、2 436、2 580、2 823、2 856位点;在氨基酸的编码区序列也发现3处位点存在差异,分别为258、274、513位点.由于杉木基因组测序未完成,关于杉木纤维素合成酶基因家族的数量、各个基因的确切功能,以及各基因间的相互作用关系均不清楚.虽然克隆出的2个基因存在相似性,但也存在差异,有可能为2个不同的基因.关于基因的功能尚需进一步验证.探索杉木中纤维素合酶基因家族的功能,特别是相关次生细胞壁形成的纤维素合酶基因,对于利用现代分子生物学手段改良杉木木材品质具有重要意义.本研究以杉木三代种子园种子发育植株为材料,利用RT-PCR技术分离克隆出1个新的杉木ClCesA基因,长度为2 955 bp,编码984个氨基酸;同时进行了生物信息学分析以及不同器官表达情况的研究;采用Invitrogen公司的gateway特异性位点重组技术将ClCesA基因构建到植物表达载体中,以期为进一步的功能解析提供依据,同时也为杉木材性的改良提供重要的基因资源.1.1 材料收集与RNA提取试验材料由国家林业局杉木工程技术研究中心提供,为杉木种子发育而成的杉木幼苗.取长势好的杉木幼嫩叶片组织0.2 g,按照TIANGEN公司的总RNA提取试剂盒操作说明书提取杉木总RNA,通过1%(质量分数)的琼脂糖凝胶电泳和超微量分光光度计Nanodrop分别检测样品的RNA质量及浓度.1.2 ClCesA基因克隆与序列分析利用NCBI网上的GenBank数据库查询到已知杉木的DNA序列,采用gateway 技术设计特异引物(表1).取检测合格的杉木总RNA,利用TIANScript cDNA第一链合成试剂盒进行第1条链cDNA的合成.合成的cDNA稀释10倍后进行PCR反应,具体反应体系见表2.反应程序为:95 ℃预变性3 min;94 ℃变性1 min,53 ℃退火3 min,72 ℃延伸3 min,共30循环;最后72 ℃延伸10 min.扩增后的产物经1%(质量分数)的琼脂糖凝胶电泳检测后,通过BP连接到pDONR207载体上,热激转化到大肠杆菌DH-5α感受态细胞,进行Gen抗性筛选.选取单菌落做菌落PCR,将检测为阳性的菌落进行液体LB过夜培养,使用天根质粒试剂盒提取质粒,并送铂尚生物公司进行全长测序.利用实时荧光定量PCR仪研究ClCesA基因在不同组织的表达情况.qPCR扩增试剂为SYBR Premix Ex Taq,采用的荧光染料为SYBR GreenⅠ,以杉木EF1α基因作为内参基因.在进行qCR分析前,先使用普通PCR检测体系的反应条件,如退火温度、模板量等,并将PCR 产物进行测序,以保证qPCR结果的准确性.1.3 植物表达载体的构建使用Gateway BP反应试剂盒,将回收带有attB位点的ClCesA基因PCR产物与入门载体pDONR207在室温下进行BP反应1 h.反应体系为:PCR产物0.6 μL,入门载体0.4 μL,BP酶0.25 μL.随后将反应产物全部转化到E.coli Trans 5α感受态细胞,涂布于含Gen(50 mg·L-1)的LB平板;次日随机挑选阳性克隆摇菌并提取质粒,经PCR检测验证正确后进行测序,即可获得入门克隆载体pDONR207-ClCesA.参考Gateway LR反应试剂盒的操作说明书,将入门载体与目的载体pGWB506进行LR反应,室温25 ℃,反应后转化大肠杆菌DH5α感受态细胞,并涂布于LB平板(含50 mg·L-1潮霉素),37 ℃培养过夜.次日挑单菌落,摇菌后用试剂盒提取质粒,经PCR检测载体构建的正确性,并进行测序验证.1. 4 纤维素合酶基因的序列分析在DNAMAN V6软件上设计引物;通过Vecter NTI软件进行序列比对;在TMHMM Server v.2.0软件上进行蛋白质结构分析;运用MEGA 5.1软件通过邻接法构建不同植物的基因进化树;采用Microsoft Excel 2007软件对基因表达进行作图.2.1 杉木总RNA提取杉木叶片总RNA经1%(质量分数)琼脂糖凝胶电泳,结果见图1.从图1可以看出,RNA条带完整,28 S rRNA和18 S rRNA 2条主带明显,28 S rRNA的亮度约为18 S rRNA的2倍,且无DNA污染.采用超微量分光光度计进行浓度测定和质量分析、检测,结果显示,样品总RNA浓度为231 μg·μL-1,D260 nm/D280nm=1.86,其浓度与质量均满足反转录的要求,可用于后续试验.最后,将总RNA 样品置于-80 ℃冰箱保存.2.2 扩增产物的碱基序列测定与分析取上述总RNA,采用天根公司的逆转录酶试剂盒获得cDNA,通过RT-PCR反应得到1个约3 000 bp的产物(图2),将该产物连接到pDONR207载体且转化DH-5α,挑取阳性克隆测序.测序结果显示该序列与CesA基因的同源性很高,并含有保守结构域.其氨基酸序列及结构如图3所示.2.3 蛋白质跨膜结构的预测从图4可以看出,新的ClCesA蛋白具有8个跨膜区,分别为175~197、204~223、760~782、794~816、831~853、878~900、910~932、945~964,最大跨膜区域为22个氨基酸残基,最小跨膜区域为12氨基酸残基,平均跨膜区域为18个氨基酸残基.2.4 基因序列系统发生树本研究运用MEGA 5.1软件构建系统进化树,将克隆到的杉木ClCesA基因编码区与拟南芥(Arabidopsis thaliana, At)、颤杨(P.tremuloides, Ptr)、欧洲山杨×颤杨(Populus tremula×P.tremuloides, Ptt)、火炬松(Pinus taeda L., Pt)、银灰杨(Populus canescens (Ait.) Smith, Pca)、毛白杨(Populus tomentosa Carrière, Pto)的初生壁和次生壁形成的CesA基因进行系统进化树分析.初生壁有关基因有AtCesA1(O48946)、AtCesA2(O48947)、AtCesA5(Q8L778)、AtCesA9(Q9SJ22)、PtrCesA7(AAO25581)、PttCesA2(AAT09895)、PtrCesA4(AAO25536)、PtrCesA6(AAP40636).次生壁有关基因有PtCesA1(AAX18647)、PtCesA2(AAX18648、PtCesA3(AAX18649)、PtrCesA2(AAM26299)、PtrCesA7(AAO25581)、PtrCesA3(AAQ08987)、PcaCesA1(AAC78476)、PttCesA1(AAT09894)、PttCesA3-1(AAT09896)、PttCesA3-2(AAT09897)、AtCesA4(Q84JA6)、AtCesA8(Q8LPK5)、PtoCesA1(AAY21910).从图5可以看出ClCesA分布在次生壁较密集的分支,且与杉木ClCesA1在一个小分支,说明两者的亲缘关系较近,可能执行类似的功能,由此推测ClCesA基因可能与次生壁合成有关.2.5 ClCesA基因在不同组织的表达杉木纤维素合酶基因在不同器官表达情况的分析结果(图6)显示,杉木的根、茎、叶3个器官均有纤维素合酶基因表达,而且相对表达量在杉木根中最高,茎中次之,叶中最低.因此,可以推测ClCesA基因与杉木木材生长发育的调控有着密切关系.2.6 植物过表达载体的构建将检测正确的入门载体pDONR207-ClCesA和目的载体pGWB506进行LR反应,获得植物表达载体pGWB506-35S-ClCesA-GFP,载体构建策略如图7所示.使用特异性的载体引物和基因特异性引物做菌落PCR,扩增出1 000 bp的片段(图8),表明LR反应成功.重组质粒的测序结果显示,ClCesA并未产生突变,并且和原序列相同,说明本试验的植物表达载体连接正确.通过序列比对与分析,Pear et al[22]首先从棉花中克隆得到与细菌CelA基因同源的纤维素合酶基因GhCesA1、2.通过基因组学分析研究,Richmond et al[23]得到了拟南芥基因组中的10个CesA基因的蛋白质序列结构特点.Suzuki et al[17]分析了杨树基因组的18个CesA基因的结构特点.目前为止,研究发现CesA基因编码的蛋白质有着相类似的结构,共有8个跨膜结构域, N端有2个, C端有6个.1个锌指结构域在N端,可能和蛋白质之间的相互作用有关.CesA蛋白都含有较典型的结构如DDDQxxRW(天冬氨酸、天冬氨酸、天冬氨酸、谷氨酞胺-x-x精氨酸—色氨酸)[24].通过对CesA基因功能的研究,已陆续鉴定出与初生壁和次生壁形成相关的基因.在欧洲颤杨中,PtrCesA1、PtrCesA2、PtrCesA3可能形成一个纤维素合成酶复合体,共同完成杨树次生壁的合成[15,24],而PtrcesA4、PtrCesA5、PtrCesA7可能与初生壁的形成有关[16].本研究克隆到一个新杉木ClCesA基因,通过系统进化树,该基因与聚类在一个进化分支并在根中表达量最高,推测该基因可能参与植物次生细胞壁的形成.而庞景等克隆到的2个CesA基因,ClcesAl与PtrCesAl的同源性最高;ClCesA2与ZmCesA2聚类在一起,两者在茎中高度表达.本研究中所用材料为培养30 d左右的杉木种子萌发而成的幼苗,而庞景所用的杉木材料为1年生杉木幼苗,且纤维素合酶CesA在植物不同发育阶段、不同器官中的表达量不同,所以这可能是造成ClCesA在杉木器官中表达量不同的原因.植物中除了CesA基因外,还有类纤维素合酶(cellulose synthase-like, CSL)基因家族参与纤维素的生物合成.拟南芥CSL蛋白有9个家族,分别为CSL(A/B/C/D/E/F/G/H/J)[25];在基因序列上,CSL基因与CesA基因具有部分同源,主要参与各种β-聚糖链的合成[23,25].彭莎莎等[21]也从杉木中克隆到了CSL基因,并分析其序列结构特点.【相关文献】[1] KUMAR M, TURNER S. Plant cellulose synthesis: CESA proteins crossing kingdoms[J]. Phytochemistry, 2015,112:91-99.[2] KAUR S, DHUGGA K, GILl K. Novel structural and functional motifs in cellulose synthase (CesA) genes of bread wheat (Triticum aestivum L.)[J]. PLoS ONE, 2016,11(1):e0147046. [3] DESPREZ T, JURANIEC M, CROWELLE E F. Organization of cellulose synthase complexesinvolved in primary cell wall synthesis in Arabidopsis thaliana[J]. Proceedings of the National Academy of Sciences, 2007,104(39):15 572-15 577.[4] TAYLOR N G. Cellulose biosynthesis and deposition in higher plants[J]. New Phytologist, 2008,178(2):239-252.[5] TANAKA K, MURATA K, YAMAZAKI M. Three distinct rice cellulose synthase catalytic subunit genes required for cellulose synthesis in the secondary wall[J]. Plant Physiology, 2003,133(1):73-83.[6] WU L, JOSHI C P, CHIANG V L. A xylem-specific cellulose synthase gene from aspen (Populus tremuloides) is responsive to mechanical stress[J]. The Plant Journal,2000,22(6):495-502.[7] DJERBIJ S, LINDSKOG M, ARVESTAD L. The genome sequence of black cottonwood (Populus trichocarpa) reveals 18 conserved cellulose synthase (CesA) genes[J]. Planta, 2005,221(5):739-746.[8] 许雷,刘一星,方连玉.大青杨纤维素合成酶PuCesA6基因cDNA的克隆及序列分析[J].西南林业大学学报,2012,32(5):26-32.[9] KRAUSKOPF E, HARRIS P J, PUTTERILL J. The cellulose synthase gene PrCESA10 is involved in cellulose biosynthesis in developing tracheids of the gymnosperm Pinus radiata[J]. Gene, 2005,350(2):107-116.[10] 刘昱翔,陈建荣,彭彦,等.苎麻纤维素合成酶基因BnCesA4 cDNA序列的克隆与表达分析[J].作物研究,2014,28(5):472-478.[11] RANIK M, MYBURG A A. Six new cellulose synthase genes from Eucalyptus are associated with primary and secondary cell wall biosynthesis[J]. Tree Physiology, 2006,26(5):545-556.[12] 阮维程,潘婷,季孔庶.马尾松纤维素合成酶基因PmCesA1的克隆及其分析[J].分子植物育种,2015,13(4):861-870.[13] 庞景,童再康,黄华宏,等.杉木纤维素合成酶基因CesA的克隆及表达分析[J].浙江农林大学学报,2015,32(1):40-46.[14] 邓小波,胡尚连,曹颖,等.绿竹与毛竹纤维素合成酶(CesA)的生物信息学分析[J].福建林学院学报,2011,31(1):84-90.[15] SAMUGA A, JOSHI C P. Cloning and characterization of cellulose synthase-like gene, PtrCSLD2 from developing xylem of aspen trees[J]. Physiologia Plantarum,2004,120(4):631-641.[16] KALLURI U C, JOSHI C P. Differential expression patterns of two cellulose synthase genes are associated with primary and secondary cell wall development in aspen trees[J]. Planta, 2004,220(1):47-55.[17] SUZUKI S, LI L, SUN Y H. The cellulose synthase gene superfamily and biochemical functions of xylem-specific cellulose synthase-like genes in populus trichocarpa[J]. PlantPhysiology, 2006,142(3):1 233-1 245.[18] SONG D, SHEN J, LI L. Characterization of cellulose synthase complexes in Populus xylem differentiation[J]. New Phytologist, 2010,187(3):777-790.[19] 蒋向辉,佘朝文,许栋,等.杉木CCoAOMT基因部分cDNA克隆与序列分析[J].中南林业科技大学学报,2009,29(6):24-28.[20] 吕运舟,郑佳,陈金慧,等.参与杉木次生壁合成调控的转录因子ClMYB4的克隆及在大肠杆菌中表达[J].分子植物育种,2012,10(5):512-519.[21] 彭沙沙,童再康,黄华宏,等.杉木纤维素合成酶类似蛋白基因ClCslD1的克隆及其生物信息学分析[J].浙江农林大学学报,2012,29(1):1-6.[22] PEAR JR, KAWAGOE Y, SCHRECKENGOET W E. Higher plants contain homologs of the bacterial celA genes encoding the catalytic subunit of cellulose synthase[J]. Proceedings of the National Academy of Sciences, 1996,93(22):12 637-12 642.[23] RICHMOND T A, SOMERVILLE C R. The cellulose synthase superfamily[J]. Plant Physiology, 2000,124(2):495-498.[24] SAMUGA A, JOSHI CP. A new cellulose synthase gene (PtrCesA2) from aspen xylem is orthologous to Arabidopsis AtCesA7 (irx3) gene associated with secondary cell wall synthesis[J]. Gene, 2002,296(1):37-44.[25] LIEPMAN A H, CAVALIER D M. The cellulose synthase-like A and cellulose synthase-like C families: recent advances and future perspectives[J]. Front Plant Sci, 2012,3(109):1-7.。
杉木涩籽形成过程的花青素还原酶基因表达

收稿日期:2018-06-12修回日期:2018-09-05基金项目:国家重点研发计划(2016YFD0600300)ꎻ国家林业局杉木工程技术研究中心平台建设项目(ptjh130002)ꎮ第一作者简介:王培(1993-)ꎬ女ꎬ硕士研究生ꎬ从事林木遗传育种研究ꎮEmail:790646055@qq.comꎮ通信作者:陈宇(1983-)ꎬ男ꎬ研究实习员ꎬ从事林木遗传育种研究ꎮEmail:28811852@qq.comꎮDOI:10.13324/j.cnki.jfcf.2019.01.001杉木涩籽形成过程的花青素还原酶基因表达王㊀培1ꎬ2ꎬ理㊀挪1ꎬ2ꎬ张家君1ꎬ2ꎬ马志慧2ꎬ3ꎬ陈㊀宇1ꎬ2ꎬ林思祖1ꎬ2(1.福建农林大学林学院ꎬ福建福州350002ꎻ2.国家林业局杉木工程技术研究中心ꎬ福建福州350002ꎻ3.福建农林大学生命科学学院ꎬ福建福州350002)摘要:为探究杉木涩籽形成过程中ꎬ花青素还原酶(ANR)基因在类黄酮生物合成途径中的分子调控机制ꎬ以种子为试验材料ꎬ通过RT ̄PCR结合RACE的方法对杉木ANR基因进行全长克隆ꎬ运用qPCR技术分析不同涩籽率的杉木家系中ANR基因在种子不同发育时期的表达模式ꎮ结果显示ꎬ克隆到的ANR基因全长cDNA序列为1357bpꎬ具有一个1056bp的开放阅读框ꎬ编码351个氨基酸ꎮ保守域预测显示ANR蛋白属于SDR超家族ꎬ具有一个FR_SDR_e特异性结构域ꎮqPCR结果发现ꎬ在高涩籽率杉木中ꎬ90㊁100d的表达量一直显著高于中㊁低涩籽率杉木ꎻ而在中㊁低涩籽率杉木中ꎬANR基因的表达模式仅仅是种子发育80d时较高ꎬ之后明显下调ꎮ表明杉木涩籽率可能与其ANR基因表达相关ꎮ关键词:杉木ꎻANR基因ꎻ基因克隆ꎻ序列分析ꎻ表达模式中图分类号:Q945.78ꎻS791.27㊀㊀㊀文献标识码:A文章编号:2096-0018(2019)01-0001-08ANRgeneexpressionduringtheformationofsterileseedsofCunninghamialanceolataWANGPei1ꎬ2ꎬLINuo1ꎬ2ꎬZHANGJiajun1ꎬ2ꎬMAZhihui1ꎬ2ꎬCHENYu1ꎬ2ꎬLINSizu1ꎬ2(1.CollegeofForestryꎬFujianAgricultureandForestryUniversityꎬFuzhouꎬFujian350002ꎬChinaꎻ2.StateForestryAdministrationEngineeringResearchCenterofChineseFirꎬFuzhouꎬFujian350002ꎬChinaꎻ3.CollegeofLifeSciencesꎬFujianAgricultureandForestryUniversityꎬFuzhouꎬFujian350002ꎬChina)Abstract:Toexploretheroleofanthocyaninreductase(ANR)geneinflavonoidbiosynthesispathwayduringtheformationofsterileseedsofChinesefirꎬthefull ̄lengthANRgenewasclonedbyRT ̄PCRcombinedwithRACEꎬandtheANRgeneexpressionindifferentdevelopmentstagesofseedswasanalyzedbyqPCR.Theresultsshowedthatthefull ̄lengthcDNAsequenceoftheclonedANRgenewas1357bpꎬitcontainsa1056bpopenreadingframeandencodes351aminoacids.ConservativedomainpredictionindicatedthatANRproteinbelongstoSDRsuperfamilyandhasaFR_SDR_especificdomain.qPCRresultsshowedthatANRgeneexpressioninfamilywithhighsterileseedratewassignificantlyhigherat90and100dofseeddevelopmentthanthefamilieswithmiddleandlowsterileseedrate.HoweverꎬtheANRgeneexpressioninfamilieswithmiddleandlowsterileseedratewasonlyhigherat80dofseeddevelopmentꎬandthendecreased.TheresultsshowedthatthesterileseedrateofChinesefirmightberelatedtoANRgeneexpression.Keywords:Chinesefir(Cunninghamialanceolata)ꎻANRgeneꎻgenecloningꎻsequenceanalysisꎻexpressionpattern杉木[Cunninghamialanceolata(Lamb.)Hook.]是我国南方地区常用的经济用材树种ꎬ在全国林业生产中占有举足轻重的地位ꎮ目前ꎬ已建立了杉木第3代种子园作为杉木的良种选育基地ꎬ但同时面临着败育率居高不下的现实问题ꎮ杉木种子败育受多种因素的影响ꎬ比如自身因素包括开花习性㊁球果的形态和自交的不亲和等ꎬ外界因素包括雨量㊁光照㊁温度和病虫害等[1]ꎮ败育种子可分为瘪籽和涩籽ꎬ瘪籽是仅有外种皮和一层珠心薄膜的种子ꎻ涩籽则与健籽的外观并无太大的差别ꎬ但剥去种皮后会发现种胚内充满了褐色物质ꎬ后证实这种褐色物质为单宁类的多酚物质[2]ꎮ在涩籽与健籽的研究中发现涩籽多形成于后胚期ꎬ成熟的干涩籽中单宁的含量比健籽高4~6倍ꎬ且过氧化酶活性显著降低ꎬMDA含量升高进一步加重了涩籽的形成[3-4]ꎮ基于前人的结论ꎬ林小琴等[5]测定了杉木成熟的涩籽与健籽的多酚含量ꎬ发现涩籽中多酚高达383.22mg g-1ꎬ而健籽中多酚仅18.36mg g-1ꎬ表明涩籽的形成与多酚物质具有一定的相关性ꎮ㊀㊀以涩籽和健籽为材料ꎬ本实验室对多酚物质形成途径ꎬ即类黄酮生物合成途径的相关基因进行了转录组测序ꎬ发现参与类黄酮物质合成的多个编码基因的表达量具有显著性的差异ꎮ其中ꎬ类黄酮生物合成途森林与环境学报㊀2019ꎬ39(1):1-8第39卷第1期JournalofForestandEnvironment2019年1月2 森㊀林㊀与㊀环㊀境㊀学㊀报第39卷㊀径中与缩合单宁合成有关的关键酶ꎬ花青素还原酶(anthocyanidinreducetaseꎬANR)其编码基因BANYULS(BANꎬ也称ANR)的表达量在涩粒中相对较高ꎮ研究表明拟南芥ANR基因参与种皮颜色的形成ꎬANR基因缺失型的拟南芥由于花青素无法转化为缩合单宁ꎬ种皮颜色加深变红ꎬ在野生型拟南芥中随着种子发育ANR基因调控种子中缩合单宁的积累[6-7]ꎮ在桑葚发育过程中ꎬMaANR基因在生长初期表达量较低ꎬ中期迅速生长且MaANR优势表达ꎬ但成熟后期表达急速下调ꎬ以达到花青素的积累[8]ꎬ说明ANR基因的表达具有一定的时序性且与花青素的积累呈负相关ꎮ然而ꎬ杉木ANR基因的相关研究尚未见报道ꎬ综合前人的研究以及本实验室前期转录组数据ꎬ推测杉木涩籽中单宁物质的积累可能是因为ANR基因抑制机制出现异常ꎬ导致花青素不断地向缩合单宁方向衍化ꎮ为了探究杉木涩籽形成过程中ꎬ花青素还原酶基因的表达相关性ꎬ以杉木种子为试验材料ꎬ运用分子克隆的手段对ANR基因进行克隆及生物信息学分析ꎬ采用qPCR技术进行表达特性分析ꎬ探究杉木涩籽形成过程中ANR基因在类黄酮生物合成途径中的分子调控机制及其作用机理ꎬ以期为以后杉木的分子辅助育种提供相关依据ꎮ1㊀材料与方法1.1㊀试验材料与试剂㊀㊀以福建省尤溪林场杉木3代种子园中的M家系(涩籽率50%)种子为试验材料ꎬ采集授粉后140d的未成熟种子ꎬ-80ħ冰箱中保存备用ꎮ由于种子发育早期种胚发育较小ꎬ很难分离出涩籽与健籽ꎬ所以采用不同涩籽率的杉木家系为试验对象ꎬ利用H家系(涩籽率90%)㊁M家系和L家系(涩籽率24%)授粉后不同时期(授粉后70㊁80㊁90㊁100d)的完整种子为试验材料ꎬ取大小一致无病害的球果进行剥离ꎬ每两个球果的种子混装于同一离心管ꎬ随后称取0.1g种子ꎬ4个生物学重复ꎬ进行基因实时荧光定量PCR分析ꎮ㊀㊀试验过程中所用到的试剂盒如下:天根多糖多酚RNA提取试剂盒ꎬ全式金公司的2ˑTransTaq ̄TPCRSuperMix和TransScriptone ̄stepgDNARemovalandcDNASynthesissuperMix试剂盒ꎬTakara公司的SMARTer®RACE5ᶄ/3ᶄKit㊁PrimeScriptTM1stStrandcDNASynthesisKit㊁胶回收试剂盒和pMD18 ̄T载体等ꎮ引物为擎科生物公司合成ꎬ大肠杆菌DH5α菌种为本实验室保存ꎬ抗生素及其他生化试剂均购自上海生工生物工程股份有限公司ꎮ1.2㊀试验方法1.2.1㊀总RNA提取和cDNA的合成㊀根据天根多糖多酚RNA提取试剂盒方法ꎬ从0.1g的杉木种子中提取总RNAꎮ分别采用PrimeScriptTM1stStrandcDNASynthesisKit合成用于基因克隆的cDNAꎬTransScriptone ̄stepgDNARemovalandcDNASynthesissuperMix合成用于qPCR扩增的cDNAꎮ1.2.2㊀杉木ANR基因的扩增㊀根据本实验室杉木转录组测序得到的ANR基因的部分序列ꎬ在其高度保守的区域进行特异性引物设计ꎬ正反向引物分别是ANR ̄F和ANR ̄Rꎮ以cDNA原液为模板进行RT ̄PCR扩增ꎬPCR反应条件为94ħ预变性2minꎻ94ħ变性30sꎬ56ħ退火30sꎬ72ħ延伸30sꎬ反应35个循环ꎻ72ħ延伸10minꎬ4ħ保存ꎮ1%琼脂糖凝胶电泳检测扩增产物ꎬ紫外切胶仪回收扩增条带ꎬ胶回收试剂盒进行纯化后连接PCR产物到pMD ̄18T载体ꎬ随后转入DH5α大肠杆菌ꎬ菌液PCR筛选阳性克隆后送至华大基因公司测序ꎮ基于测序的保守片段ꎬ设计特异引物扩增ANR基因的5ᶄ和3ᶄ端序列ꎬ以上试验参照了SMARTer®RACE5ᶄ/3ᶄKit试剂盒说明书ꎮ利用巢式PCR进行ANR5ᶄRACE和3ᶄRACE实验ꎬPCR产物回收测序ꎬ将克隆到的5ᶄ㊁3ᶄ端和保守区序列用DNAMAN软件进行拼接ꎬ获得ANR基因的全长cDNA序列ꎮ利用ORFfinder工具确定ANR基因的ORF区域ꎬ再根据全长cDNA序列设计扩增ORF序列的特异性引物ꎬPCR反应后回收产物ꎬ再次测序验证ORF序列是否准确ꎮ以上实验所用到的引物如表1所示ꎮ表1㊀引物设计Table1㊀Primerdesign引物名称Primername引物序列PrimersequencePrimername引物序列Primersequence引物名称ANR ̄F5ᶄGGATTCATCGGCTCCTCCCTTGT3ᶄANR ̄GSP55ᶄGGATCACCCGTCTCACTGTTTTTG3ᶄANR ̄R5ᶄGCCCTTCCCTTCATATTCCTTCA3ᶄANR ̄NGSP55ᶄCTAGAGTTCCATCAATGGCAGGCT3ᶄANR ̄GSP35ᶄGGTAACAGTGTTGGCAGCTCTTGTGGG3ᶄANR ̄ORF ̄F5ᶄCCAAAGAAGGTGAGTCGAAAGG3ᶄANR ̄NGSP35ᶄCCTTCATTCACGCCCACATTTCCCCAC3ᶄANR ̄ORF ̄R5ᶄGAAGCATAAATTGATAGGCGAA3ᶄ1.2.3㊀生物信息学分析㊀利用ORFfinder分析ANR基因全长序列的ORF区域ꎬ并将其氨基酸序列在Blastp上进行同源性比对ꎬNCBIconserveddomains在线网站进行蛋白结构域预测ꎬ利用MEGA5软件构建系统进化树ꎻExPASyProtParam软件推测蛋白质的基本理化性质ꎬProtScale软件预测蛋白亲水性ꎬTMHMM软件进行蛋白质跨膜区预测ꎮ并用SignalP4.1Server软件预测蛋白质信号肽ꎬPredictProtein软件预测蛋白质的二级结构ꎬPhyre2软件预测蛋白质三级结构ꎬProtComp9.0软件进行蛋白质亚细胞定位预测ꎬMotifScan进行蛋白质的功能位点预测ꎮ1.2.4㊀ANR基因实时荧光定量PCR分析㊀由克隆的cDNA序列设计qPCR特异性引物ANR ̄1:5ᶄCGTG ̄GCAAAGACATTAGCGGA3ᶄꎬANR ̄2:5ᶄACTCTCTGGGTTGCCTGTGATTG3ᶄꎬ选择杉木EF1α为内参基因ꎬ设计引物EF1:5ᶄGCACTGTTATTGATGCTCCT3ᶄꎬEF2:5ᶄCCAGTGGTGGAGTCAATGAT3ᶄꎮ以稀释10倍的cDNA为模板进行qPCR扩增ꎬPCR反应条件为94ħ预变性30sꎻ94ħ变性5sꎬ54ħ退火15sꎬ72ħ延伸10sꎬ反应40个循环ꎻ72ħ延伸10sꎮ运用2-әәCt计算方法进行ANR基因的相对表达量的分析ꎮ2㊀结果与分析2.1㊀总RNA提取和cDNA的合成㊀㊀采用1.0%琼脂糖凝胶电泳法和酶标仪法ꎬ检测杉木种子总RNA的完整性和纯度ꎮ凝胶电泳结果表明28S和18S条带清晰ꎬ说明RNA无降解ꎬ完整性良好果均[图1(a)]ꎻ测定的OD260/280的值在1.8~2.0范围内ꎬ说明无DNA污染ꎮ2.2㊀杉木ANR基因的扩增㊀㊀合成cDNA第一条链ꎬ以cDNA为模板进行保守区中间序列扩增ꎬPCR产物电泳检测结果如图1(b)所示ꎬ得到一条长度约850bp的条带ꎬ回收产物ꎬTA克隆测序ꎮ按照SMARTer®RACE5ᶄ/3ᶄKitUserManual的方法ꎬ分别扩增了5ᶄ端和3ᶄ端基因序列ꎬ1%琼脂糖凝胶电泳检测RACE产物ꎮ㊀㊀3ᶄRACE和5ᶄRACE扩增分别得到两条长度约820和450bp的目的基因条带[图1(c)]ꎬ回收产物连接到克隆载体进行测序ꎮ将克隆到的基因片段进行拼接ꎬ获得了杉木ANR基因的全长cDNA序列ꎬ其长度为1357bpꎬ具有一个1056bp的开放阅读框(ORF)ꎬ编码351个氨基酸ꎮ㊀㊀根据全长序列ꎬ设计引物ꎬ进行ORF序列的PCR扩增[图1(d)]ꎬ得到一条长度约1100bp的条带ꎬ回收测序验证ORF序列的准确性ꎮ将ORF编码的氨基酸序列(图2)在NCBI上进行Blastp分析ꎬ结果显示:杉木ANR蛋白与银杏的ANR蛋白比对ꎬ同源性达到48%ꎻ与胡杨㊁毛果杨和黑杨的同源性均达到42%ꎻ与无油樟和蓖麻的同源性达到43%ꎮ同时ꎬ比对结果显示ꎬANR蛋白与银杏㊁大豆和拟南芥的DFR蛋白具有同源性ꎬ分别是44%㊁40%和41%ꎮ在DevicM对花青素还原酶基因的研究中ꎬ发现BAN基因与DFR基因具有高度的相似性ꎬ而且是种皮早期分化和发育的一种标记基因[9]ꎮThirug ̄nanasambanthamK选用VvDFR的F链为模板ꎬ运用SP2server构建CsANR模型ꎬ发现CsANR与VvDFR模型叠加[10]ꎮ表明杉木ANR可能与DFR蛋白的功能相似ꎮ图1㊀ANR基因的克隆结果Figure1㊀CloningresultsoftheANRgene3 ㊀第1期王培ꎬ等:杉木涩籽形成过程的花青素还原酶基因表达图2㊀核苷酸和推导的氨基酸序列Figure2㊀Nucleotideanddeducedaminoacidsequenc2.3㊀生物信息学分析2.3.1㊀蛋白质基本性质及结构预测㊀根据推导的氨基酸序列ꎬ采用ExPASyProtParam软件推测ANR蛋白的理化性质ꎬ结果显示:351个氨基酸编码的ANR蛋白分子量为38824.37g mol-1ꎬ等电点为5.70ꎬ半衰期为30hꎬ不稳定指数44.24ꎬ属不稳定蛋白质ꎮ其中亮氨酸含量最高ꎬ达10.0%ꎬ其次ꎬ丙氨酸为7.7%㊁丝氨酸7.4%ꎬ带负电残基的总数为38个ꎬ带正电残基的总数为28个ꎬ预测亲水性平均值为-0.068ꎬ表示ANR蛋白为亲水性蛋白ꎮProtScale预测结果显示[图3(a)]负值多于正值ꎬ也证明了该蛋白的亲水性质ꎮ㊀㊀运用SignalP4.1Server信号肽预测如图3(b)所示ꎬ发现ANR蛋白N末端并没有信号肽ꎬTMHMM跨膜区预测如图3(c)所示ꎬ也未发现跨膜结构ꎬ说明ANR蛋白不属于分泌或跨膜蛋白ꎬ有可能是以游离态形式存在ꎮ运用ProtComp9.0亚细胞定位预测结果如表2所示ꎬ发现ANR蛋白最有可能存在于细胞外ꎬ其次是在细胞器过氧化物酶体和叶绿体中ꎮ㊀㊀蛋白质结构预测结果显示:ANR蛋白具有12个蛋白质-蛋白质相互作用位点ꎬ二级结构中无规则卷曲(Cc)和α螺旋(Hh)占比较大ꎬ分别是40.17%和38.18%ꎬ其次是延伸链(Ee)占14.81%ꎬβ转角(Tt)占6.84%[图3(d)]ꎮ三级结构预测结果[图3(e)]与二级结构预测结果基本一致ꎮ2.3.2㊀ANR蛋白结构域及功能位点分析㊀经过NCBIconserveddomains网站搜索ꎬANR蛋白属于SDR超家族成员(图4)ꎬ具有一个FR_SDR_e特异性结构域ꎬ以及还原酶典型的NADP结合位点㊁活动位点以及底物结合位点ꎮSDR是属于功能多样的氧化还原酶家族ꎬ主要有经典型㊁扩展型㊁中间型㊁复合型和扩散型[11]ꎮFR_SDR_e是包含了扩展型SDR特性的FR蛋白ꎬ在类黄酮物质发生NADP依赖性还原反应中发挥着重要作用ꎮ用MotifScan软件预测蛋白的功能位点ꎬ显示ANR蛋白具有4个N ̄肉豆蔻酰化位点(N ̄myris ̄toylationsite)ꎬ5个蛋白激酶C磷酸化位点(ProteinkinaseCphosphorylationsite)和6个酪蛋白激酶磷酸化位点(Caseinkinasephosphorylationsite)ꎬ具体位置如表3所示ꎮ 4 森㊀林㊀与㊀环㊀境㊀学㊀报第39卷㊀图3㊀ANR蛋白质基本性质及结构预测Figure3㊀BasicpropertiesandstructurepredictionofANRprotein表2㊀ANR蛋白亚细胞定位结果Table2㊀SubcellularlocalizationofANRprotein亚细胞定位Subcellularlocalization细胞外Extracellular过氧化物酶体Peroxisomal叶绿体Chloroplast线粒体Mitochondrial细胞质Cytoplasmic高尔基Golgi内质网Endoplasmretic存在概率Probability2.482.271.921.060.970.960.34表3㊀ANR蛋白功能位点预测Table3㊀PredictionoffunctionalsitesofANRprotein功能位点Functionalsites氨基酸位置Aminoacidlocation氨基酸序列Aminoacidsequence功能位点Functionalsites氨基酸位置Aminoacidlocation氨基酸序列AminoacidsequenceN ̄myristoylation13 18GGFIGSproteinkinaseCphosphorylation311 313SYKN ̄myristoylation68 73GSFYSAcaseinkinasephosphorylation34 37TCRDN ̄myristoylation216 221GNPESLcaseinkinasephosphorylation72 75SAIEN ̄myristoylation323 328GIQYAKcaseinkinasephosphorylation85 88TSMDproteinkinaseCphosphorylation4 6TKKcaseinkinasephosphorylation169 172TLAEproteinkinaseCphosphorylation34 36TCRcaseinkinasephosphorylation267 270SIADproteinkinaseCphosphorylation117 119TVRcaseinkinasephosphorylation345 348SSNDproteinkinaseCphosphorylation302 304SNK 5㊀第1期王培ꎬ等:杉木涩籽形成过程的花青素还原酶基因表达图4㊀ANR蛋白保守域分析Figure4㊀AnalysisofconserveddomainofANRprotein2.3.3㊀ANR蛋白系统进化树的构建㊀根据Blastp的蛋白序列比对结果ꎬ选取9种不同物种的蛋白序列ꎬ其中包括可可㊁木本棉㊁苦荞麦㊁毛果杨㊁黑杨㊁茶树㊁拟南芥㊁银杏㊁北美云杉等物种ꎮ运用MEGA5软件中的NJ法进行ANR蛋白的系统进化树构建(图5)ꎮ从进化关系上可以看出ꎬ杉木ANR蛋白与北美云杉位于同一分支ꎬ表明两者的亲缘关系最为接近ꎬ但北美云杉的蛋白质特征还未知ꎬ仅具有花青素还原酶的CDD区ꎬ区域特征属于FR_SDR_e结构域ꎮ银杏和拟南芥则各位于一支ꎮ目前裸子植物的ANR研究报道少之又少ꎬ2006年首次对银杏的ANR基因进行分离与鉴定ꎬ发现其在叶子中表达最强[12]ꎮ而拟南芥中ꎬ最早在种皮中发现BAN基因的存在ꎬ证明了其在类黄酮合成途径中作为一种负调节剂来防止种皮色素的积累[13]ꎮ㊀㊀其他物种聚为一个大分支ꎬ分别包括3组ꎬ一组茶树ꎬ一组毛果杨和黑杨ꎬ另一组可可㊁木本棉和苦荞麦ꎬ根据距离远近可以看出这些物种与杉木的亲缘关系相对较远ꎮ蒋春丽等[14]的研究表明ꎬ在苦荞麦不同发育时期的种子中ꎬANR基因的表达受不同光质的影响ꎬ蓝光和远红光可导致ANR表达量的上调ꎮ在可可中ꎬ种子的ANR表达量仅次于外果皮ꎬ但PA含量却最低ꎬ结果表明PA的积累与ANR不相关ꎬ而与LAR相关ꎻ对拟南芥BAN突变体进行TcANR基因过表达时ꎬ发现突变体种子恢复到野生型的表型ꎬ种皮内有缩合单宁的积累ꎬ同时下胚轴花青素的含量有所降低[15]ꎮ缩合单宁在种皮细胞内积累的主要作用则是保护胚和胚乳的正常发育[16]ꎮ在过表达CsANR的烟草中ꎬ发现CsANR提高了烟草根对铝毒害的耐受性[17]ꎮ结合以上研究发现ꎬANR基因可能具有一定的物种特异性ꎬ其表达特性受多种因素干扰ꎬ但进化树亲缘关系的远近一定程度上表明了杉木ANR可能具备的蛋白功能ꎬ亦可为后期ANR蛋白功能的验证提供一定的理论参考ꎮ图5㊀杉木ANR蛋白系统进化树Figure5㊀PhylogenetictreeofANRproteininCunninghamialanceolata2.4㊀ANR基因相对表达量的分析㊀㊀对不同涩籽率的杉木家系不同发育时期的种子进行实时荧光定量PCR分析ꎬ结果显示(图6)ꎬ在70~80d之间ꎬ3个家系的ANR基因的相对表达量都显著上升ꎬ且在70d时ꎬ各家系中ANR基因的相 6 森㊀林㊀与㊀环㊀境㊀学㊀报第39卷㊀对表达量大小依次是H>M>Lꎬ差异显著ꎻ80d时ꎬM家系的表达量最大ꎬL和H家系大致相同ꎮ80~90d之间ꎬL和M家系的ANR基因相对表达量迅速降低ꎬH家系的相对表达量变化则没有明显的差异ꎮ90~100d之间ꎬL家系的ANR基因相对表达量有所上升ꎬM家系的继续下降ꎬ到达100d时两个家系的ANR基因相对表达量基本一致ꎬ但是H家系的相对表达量仍然处于高水平状态ꎬ且无降低的趋势ꎮ表明在H家系种子中ꎬANR基因在90㊁100d的表达量一直显著高于M和L家系ꎻ而在M和L家系中ꎬANR基因的表达模式仅仅是种子发育80d时较高ꎬ之后明显下调ꎮ说明杉木涩籽率可能与其ANR基因表达相关ꎮ3㊀讨论㊀㊀本试验通过RT ̄PCR结合RACE的方法ꎬ克隆了杉木花青素还原酶ANR基因的全长cDNA序列ꎬ㊀㊀注:大写字母表示同一家系不同发育时期差异显著ꎬ小写字母表示同一发育时期不同家系差异显著(P<0.05)Note:capitallettersindicatesignificantdifferenceamongdifferentdevelopmentstagesofthesamefamilyꎬlowerlettersindicatesignificantdifferenceamongdifferentfamiliesatthesamedevelopmentstage(P<0.05).图6㊀种子不同发育时期ANR基因的相对表达量Figure6㊀RelativeexpressionofANRgeneindifferentdevelopmentstagesofseeds其长度为1357bpꎬ具有一个1056bp的开放阅读框(ORF)ꎬ编码351个氨基酸ꎮ蛋白质基本性质预测:ANR蛋白分子量38824.37Daꎬ理论等电点5.70ꎬ预测的半衰期为30hꎬ不稳定指数44.24ꎬ亲水性平均值-0.068ꎬ属于不稳定的亲水性蛋白质ꎬ且不具备信号肽和跨膜区结构ꎮ保守域预测ANR蛋白为SDR超家族ꎬ具有一个FR_SDR_e特异性结构域ꎮ在进行同源性分析时ꎬ发现杉木ANR蛋白和其他物种的ANR蛋白和DFR蛋白都具有一定的同源性ꎮ㊀㊀本研究对不同涩籽率的杉木家系不同发育时期的种子进行ANR基因实时荧光定量PCR分析结果显示ꎬL家系㊁M家系和H家系的种子在70~80d的过程中ꎬANR基因均是显著上调表达ꎬL家系和M家系在80d之后ANR基因表达水平显著降低ꎬ但H家系在80d以后其表达量一直保持较高水平ꎮ㊀㊀有研究表明ꎬ黄醇酮对矮牵牛和玉米的花粉萌发和花粉管发育具有至关重要的作用[18-19]ꎬ在烟草中沉默黄酮醇合酶(FLS)基因ꎬ则使得烟草的黄酮醇含量降低ꎬ黄烷 ̄3 ̄醇的含量升高ꎬ且花粉萌发受阻ꎬ种子数量减少[20]ꎮ在类黄酮生物合成途径中ꎬFLS和DFR蛋白酶是以二氢黄酮醇作为共同底物ꎬ分别进行黄酮醇和无色花青素的合成ꎬ无色花青素又经ANS和ANR的催化转化为黄烷 ̄3 ̄醇ꎬ即单宁单体[21]ꎮ由于ANR是催化花青素向缩合单宁衍化的关键酶之一ꎬANR基因的表达对黄烷 ̄3 ̄醇和PA的生物合成起着重要作用[22]ꎮ在杉木涩籽形成过程中ꎬ可能因为杉木ANR与DFR具有一定的相似性ꎬ促使二氢黄酮醇向PA途径的方向衍化ꎬANR基因的高度表达使得黄酮醇的合成受到抑制ꎬ进而导致花粉育性降低引发了种子败育现象ꎮ因此推测杉木H家系的高败育率可能跟ANR基因的异常表达具有一定的关联性ꎮ但在BAUDRYetal[23]研究中发现ꎬTT2(MYB转录因子)㊁TT8(bHLH转录因子)和TTG1(WDR蛋白)基因是BAN基因正确表达的必要条件ꎮ戴鹏辉等[24]对珙桐的DiMYB1进行表达分析时ꎬ发现DiMYB1在败育种子中的表达量显著高于正常种子ꎬ而且DiMYB1在发育中期的败育种子中表达量最高ꎮ所以在杉木涩籽中ꎬANR基因的异常表达可能是这类基因单一或共同调控的结果ꎮ而杉木ANR基因的调控机制到底是如何进行的ꎬ还需后期进行实验加以验证ꎮ参考文献[1]符梅忠ꎬ傅远志ꎬ蒋国洪.杉木籽粒败育的研究[J].浙江林业科技ꎬ1989ꎬ9(2):15-19.[2]管康林.杉木种子空瘪涩粒成因的综合分析[J].浙江林学院学报ꎬ1997ꎬ14(1):88-93.[3]黄烺增ꎬ谢永强ꎬ齐清琳ꎬ等.杉木种子涩籽胚败育过程及其生理生化指标[J].福建林学院学报ꎬ1993ꎬ13(4):401-406.[4]黄烺增.杉木种子涩籽胚败育与过氧化物酶同工酶活性及脂质过氧化作用的关系[J].西北植物学报ꎬ1996ꎬ16(5):83-88. 7 ㊀第1期王培ꎬ等:杉木涩籽形成过程的花青素还原酶基因表达8 森㊀林㊀与㊀环㊀境㊀学㊀报第39卷㊀[5]林小琴ꎬ马志慧ꎬ何宗明ꎬ等.杉木种子多酚提取工艺优化[J].森林与环境学报ꎬ2017ꎬ37(1):47-53. [6]XIEDYꎬSHARMASBꎬPAIVANLꎬetal.RoleofanthocyanidinreductaseꎬencodedbyBANYULSinplantflavonoidbiosynthesis[J].Scienceꎬ2003ꎬ299(5605):396-399.[7]DEBEAUJONIꎬNESINꎬPEREZPꎬetal.Proanthocyanidin ̄accumulatingcellsinArabidopsistesta:regulationofdifferenti ̄ationandroleinseeddevelopment[J].ThePlantCellꎬ2003ꎬ15(11):2514-2531.[8]李军ꎬ梁燕梅ꎬ赵爱春ꎬ等.桑树花青素还原酶基因MaANR的克隆和表达分析[J].蚕业科学ꎬ2016ꎬ42(4):570-575. [9]DEVICMꎬGUILLEMINOTJꎬDEBEAUJONIꎬetal.TheBANYULSgeneencodesaDFR ̄likeproteinandisamarkerofearlyseedcoatdevelopment[J].ThePlantJournalꎬ1999ꎬ19(4):387-398.[10]THIRUGNANASAMBANTHAMKꎬMURALIDARANSꎬMANDALAKA.MolecularcloningꎬcomputationalandexpressionanalysisofanthocyanidinreductaseinTea(Camelliasinensis)[J].AppliedBiochemistryandBiotechnologyꎬ2014ꎬ174(1):130-145.[11]OPPERMANNUꎬFILLINGCꎬHULTMꎬetal.Short ̄chaindehydrogenases/reductases(SDR):the2002update[J].Chemico ̄BiologicalInteractionsꎬ2003ꎬ143/144:247-253.[12]SHENGAꎬPANGYZꎬWUWSꎬetal.IsolationandcharacterizationofaputativeanthocyanidinreductasegenefromGinkgobiloba[J].JournalofPlantPhysiologyꎬ2006ꎬ163(2):224-227.[13]ALBERTSꎬDELSENYMꎬDEVICM.BANYULSꎬanovelnegativeregulatorofflavonoidbiosynthesisintheArabidopsisseedcoat[J].ThePlantJournalꎬ1997ꎬ11(2):289-299.[14]姜春丽ꎬ苏新堯ꎬ许亚春ꎬ等.光照对药食同源植物苦荞麦种子中原花青素积累关键酶FtANRꎬFtLAR基因表达的影响[J].中国中药杂志ꎬ2018ꎬ43(3):469-477.[15]LIUYꎬSHIZꎬMAXIMOVASꎬetal.ProanthocyanidinsynthesisinTheobromacacao:genesencodinganthocyanidinsynthaseꎬanthocyanidinreductaseꎬandleucoanthocyanidinreductase[J].BMCPlantBiologyꎬ2013ꎬ13:202. [16]LEPINIECLꎬDEBEAUJONIꎬROUTABOULJMꎬetal.Geneticsandbiochemistryofseedflavonoids[J].AnnualReviewofPlantBiologyꎬ2006ꎬ57:405-430.[17]KUMARVꎬYADAVSK.OverexpressionofCsDFRandCsANRenhancedflavonoidsaccumulationandantioxidantpotentialofrootsintobacco[J].PlantRootꎬ2013ꎬ7:65-76.[18]MOYꎬNAGELCꎬTAYLORLP.Biochemicalcomplementationofchalconesynthasemutantsdefinesaroleforflavonolsinfunctionalpollen[J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmericaꎬ1992ꎬ89(15):7213-7217.[19]YLSTRABꎬTOURAEVAꎬMORENORMBꎬetal.Flavonolsstimulatedevelopmentꎬgerminationꎬandtubegrowthoftobaccopollen[J].PlantPhysiologyꎬ1992ꎬ100(2):902-907.[20]MAHAJANMꎬAHUJAPSꎬYADAVSK.Post ̄transcriptionalsilencingofflavonolsynthasemRNAintobaccoleadstofruitswitharrestedseedset[J].PLoSOneꎬ2011ꎬ6(12):e28315.[21]乔小燕ꎬ马春雷ꎬ陈亮.植物类黄酮生物合成途径及重要基因的调控[J].天然产物研究与开发ꎬ2009ꎬ21(2):354-360ꎬ207.[22]ZHUYꎬWANGHYꎬPENGQZꎬetal.Functionalcharacterizationofananthocyanidinreductasegenefromthefibersofuplandcotton(Gossypiumhirsutum)[J].Plantaꎬ2015ꎬ241(5):1075-1089.[23]BAUDRYAꎬHEIMMAꎬDUBREUCQBꎬetal.TT2ꎬTT8ꎬandTTG1synergisticallyspecifytheexpressionofBANYULSandproanthocyanidinbiosynthesisinArabidopsisthaliana[J].ThePlantJournalꎬ2004ꎬ39(3):366-380.[24]戴鹏辉ꎬ任锐ꎬ曹福祥ꎬ等.珙桐MYB转录因子DiMYB1基因的克隆及表达分析[J].植物生理学报ꎬ2016ꎬ52(8):1255-1262.(责任编辑:叶宝鉴)㊀㊀。
木质素生物合成酶CCR基因的生物信息学分析

木 质 素 生 物合 成 酶 C C R基 因 的 生 物 信 息 学 分 析
陈 刚 , 徐秉 良 , 白江 平
( 1 .甘 肃 农业 大学 草 业 学 院 、 草 业 生态 系统 教 育 部 重 点 实 验 室 、 中 一美 草 地 畜 牧 业 可持 续 发 展 研 究 中 心 、 甘 肃 省草 业 工 程 实 验 室
l i g n i n b i o s y n t h e s i s i n p l a n t s
CHE N G a n g , ’ XU B i n g — l i a n g ’ 。 B A I J i a n g - p i n g 。 ’
( 1 .C o l l e g e o f G r a s s l a n d S c i e n c e f o G a n s u A g r i c u l t u r a l U n i v e r s i t y / K e y L a b o r a t o r y f o G r a s s l a n d E c o s y s t e m o f MO E /S i n o— S .C e n t e r s f o r G r a z i n g L a n d E c o s y s t e m S u s t a i n a b i l i t y / P r a t a c u h u r a l E n g i n e e r i n g ab L o r a t o r y o f G a n s u P r  ̄i n c e , L a n z h o u 7 3 0 0 7 0, C h i n a ;
中 图分 类号 : Q 5 1 8 . 2 文 献标 识 码 : C 文章编号 : 1 6 7 2— 5 5 6 5 ( 2 0 1 3 ) 一 0 1 — 0 5 0— 0 8
杉木天冬酰胺合成酶基因的克隆与生物信息学分析

中丝氨酸 15 个,苏氨酸 7 个,酪氨酸 3 个,其二级及三级结构主要由无规则卷曲和α-螺旋组成。关键词:杉Βιβλιοθήκη ;天冬酰胺合成酶基因;克隆;生物信息学
中图分类号 S791.27
文献标识码 A
文章编号 1007-7731(2017)22-0028-03
Cloning and Bioinformatics Analysis of Asparagine Synthetase Gene from Cunninghamia lanceolata
天冬酰胺合成酶是广泛存在于植物体内的一类由小 基因家族编码的氨基转移酶,主要以氨或谷氨酰胺及天 冬氨酸为底物催化合成天冬酰胺。根据其底物的不同, 天冬酰胺合成酶可分为 AS-A 和 AS-B,前者为氨依赖性, 仅存于原核生物;后者为谷氨酰胺依赖性,广泛存在于原 核和真核生物[1]。因此,普遍认为植物天冬酰胺主要由 AS-B 催化合成。目前,已从大豆[2-3]、小黑杨[4]、玉米[5]、海 岸松[6-8]、桑树[9]、小黑麦[10]、向日葵[11]等多种植物中成功 分离 ASN 基因,且进一步的功能研究发现 ASN 基因具有 一定的抗逆功能。杉木作为我国南方重要的造林树种, 已成功分离出一个 ASN1 基因[12],本研究将在此基础上进 一步分离杉木 ASN 基因其他新成员,以期为后续系统研 究 ASN 基因在杉木生长发育及逆境过程中的调控作用奠 定基础。
1 材料与方法
1.1 材料 供试材料为国家林业局杉木工程技术研究中 心保存的杉木组培苗 105。 1.2 方法 1.2.1 RNA 提取和 cDNA 合成 参考官清娜[12]方法进行 杉木总 RNA 提取,经酶标仪和琼脂糖凝胶(1%)检测合格
的 RNA 用于 cDNA 合成。根据 TIANScript cDNA 第一链 合成试剂盒(TIANGEN)说明书反转录 cDNA。 1.2.2 引物设计与 PCR 扩增 将已扩杉木 ASN1 基因序 列与杉木转录组中预测的 ASN 基因序列比对,设计上游 引物 5’-ATGAAAGCTTTGCACGATGACTG-3’和下游引 物 5’-TTAACCCTGAATGACAACTCCCCT-3’用于 PCR 扩 增。扩增反应程序:94℃预变性 3min;94℃变性 30s,60℃ 退火 30s,72℃延伸 2min,35 个循环;72℃延伸 7min。将 PCR 扩增产物送北京六合华大基因生物技术公司测序。 1.2.3 生物信息学分析 利用在线软件对目的基因所编 码蛋白的理化性质(ExPASy)、信号肽(SignaIP 4.1 Serv⁃ er)、跨膜结构(ExPASy)、亚细胞定位(PSORT)、磷酸化位 点(NetPhos 2.0 Sever)、卷曲螺旋结构(EMBnet Coils)、二 级结构(GOR IV)和三级结构(SWISSMODEL)进行预测 和分析。
杉木纤维素合成酶类似蛋白基因ClCslD1的克隆及其生物信息学分析

的 c N 序 列 , 长 度 为 4 10b ,开 放 阅 读 框 架 为 336b ,编 码 l12个 氨 基 酸 并 含 有 保 守 的 D,D DA 5 p 9 p 3 .D.
Q R 天冬氨 酸 残基 序列 应 用 美国 国家 生物技 术信 息 中心 ( B ) XX W NC I数据 库 对获得 的基 因进 行序 列分析 比对 .发现 它属 于 f D基 因家族 ,命 名 为 CC I 基 因。多重序 列 比对 结果 分析表 明 ,该蛋 白与 来 自颤杨 P p lst muo e . Is D1 o uu r li s e d 水稻 O yast a rz ai .拟 南芥 Arbdp i tain v io s l a等 的 同源基 因相似 性 高达 7 %以上 利 用生物 软件 对 其进行 生 物 学分 a sh a l
mo o o s s q e c s o h e A g n s A l p e c mp rs n s q e c n l sswa o d ce o c re t l g u e u n e ft e C s e e . mu t l o a io e u n e a a y i s c n u t d c n u r n l i y
w t iifr t t o st n lz h b an d s q e c . s l ftes q e c n lsss o dt a hs i bonomai meh d o a ay ete o tie e u n e Re ut o e u n ea ay i h we h tti h c s h
析 .为进 一 步研 究其在 植 物纤 化 中的功 能奠 定基 础 图 5参 l O
关键 词 :林木 育种 学 ;杉木 ;纤 维素合 成 酶类 似蛋 白;克 隆 ;表达 中图分 类 号 :¥ 2 .;Q 8 7 21 7 5 文 献标 志码 :A 文章 编号 :2 9 .7 6 2 1 ) 10 0 —6 0 50 5 (0 2 0 .0 10
木质素合成酶基因F5H的克隆及其鉴定

木质素合成酶基因F5H的克隆及其鉴定
付月;李学龙;薛永常
【期刊名称】《安徽农业科学》
【年(卷),期】2006(034)008
【摘要】通过RT-PCR从正在分化的2年生欧美杨107号茎的次生木质部提取mRNA,扩增出一基因片段,与pGM-T载体连接,重组质粒经EcoRI酶切、特异性引物扩增、测序鉴定.结果表明,该基因片段长度为867 bp,与EMBL核酸序列库中Van Doorsselaere等发表的白杨杂种F5H的cDNA序列(PAJ10324)相比,同源性高达99.19%,可以断定所扩增的DNA片段即为F5H基因片段.
【总页数】3页(P1553-1555)
【作者】付月;李学龙;薛永常
【作者单位】大连轻工业学院生物与食品工程学院,辽宁,大连,116034;大连轻工业学院生物与食品工程学院,辽宁,大连,116034;大连轻工业学院生物与食品工程学院,辽宁,大连,116034
【正文语种】中文
【中图分类】Q785
【相关文献】
1.苎麻木质素合成酶CAD基因的cDNA克隆及表达分析 [J], 朱伟溢;陈建荣;彭彦;张学文;郭清泉;赵燕
2.毛白杨木质素合成酶基因F5H克隆与生物信息学分析 [J], 王轶男;陈雪;盖颖
3.柳枝稷木质素合成酶基因F5H核心片段的克隆与表达特性的分析 [J], 李园;王珊;方程
4.欧美杨107木质素生物合成酶CCR基因的克隆及其鉴定 [J], 薛永常;赵一玲;赵文超;王雪霞
5.木质素合成酶基因F5H的RT-PCR引物筛选 [J], 付月;薛永常;姚奇
因版权原因,仅展示原文概要,查看原文内容请购买。
欧美杨107木质素生物合成酶CCR基因的克隆及其鉴定

T栽 体连接 , 重组 质粒 经 限制性 内切酶 酶切 、 异 引物 P R扩增 和测序 鉴 定。 结果表 明该扩 增 片 特 C
段 长 为 9 1p含 一 个编码 31 基酸 的 完整开放 阅读框 , 6b , 0氨 其核 苷酸序 列与 Lpe 发表 的 白杨 杂种 el等 Ppl ioa aC R的 e N ouu tc cr C srh p D A序 列 (J296 A248 )同 源性 达 到 9 . 82% , 并且反 向插 入 到 p 2 MD 0一T 载体 , 因此构 建 了 C R的反 义大肠杆 茵表 达 载体 , C 为后 续的 基 因操 作 打 下基础 。 该基 因序 列 已提 交 国际核 苷酸序 列库 , 录号 为 A 9 19 。 登 M 268 关 键 词 : C R P R; 隆 ; C R;T— C 克 木质 素
‘47” dvl i cn y m R A, db l e t p D 0一Tvc rte et e yr tco ny eP R a — ' /6 )ee pn s odxl m N a ec ndio M 2 7 o ge e n o n et , ni nid b sii ezm ,C m o h d f i e r tn
p t a o d sq e cn .T es qln eo tea pi e DN f g e t a 9 1b s an a d c nan a 0R i ei n l l t n a e u n ig h e l c f h m l d e i f A i m n w s 6 aep i n o tie n ' a d Fw iht n lt hc r s e a a a p dn w t 0 m n cd .A di w sc n it MD 0一T v co nr v r mt i 3 1 a i o a is n t a l e no p 2 h od e tri e es d e 0 .A i me t i D A , q// e o n l n n w t e N  ̄ l l f g h ec
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第30卷第6期浙江林业科技Vol. 30 No.6 2 0 1 0年11月JOUR. OF ZHEJIANG FOR. SCI. & TECH. Nov., 2 0 1 0 文章编号:1001-3776(2010)06-0006-06杉木木质素合成酶基因CCR的克隆和序列分析陈喜,童再康*,黄华宏,林二培,程龙军(浙江农林大学林业与生物技术学院,浙江临安 311300)摘要:以杉木(Cunninghamia lanceolata)茎叶为材料,根据已报道的肉桂酰辅酶A还原酶(Cinnamoyl-CoA reductase,CCR)基因保守区设计简并引物,采用RT-PCR结合RACE的策略克隆获得杉木CCR基因的全长cDNA (命名为CLCCR),该基因cDNA序列全长1 379 bp,具有一个975 bp的开放阅读框(Open Reading Frame,简称ORF),编码324个氨基酸,序列分析表明:杉木CCR基因核苷酸序列与已报道的欧洲云杉(Picea abies,CAK18610.1)CCR基因具有87%的相似性;与其他植物氨基酸序列聚类分析显示杉木与松科植物欧洲云杉、马尾松(Pinus massoniana)、火炬松(P. taeda)的亲缘关系较近。
关键词:木质素;cDNA克隆;杉木;序列分析中图分类号:S791.27 文献标识码:AGene Cloning and Sequence Analysis of CCR inCunninghamia lanceolata Lignin SynthesisCHEN Xi,TONG Zai-kang,HUANG Hua-hong,LIN Er-pei,CHENG Long-jun(School of Forestry and Bio-technology, Zhejiang A & F University, Lin’an 311300, China)Abstract: Cinnamoyl-CoA reductase(CCR) is one of the key enzymes in the process of lignin synthesis. A full length cDNA from the Cunninghamia lanceolata was amplified by RT-PCR and RACE. The full cDNA named CLCCR which was 1379 bp in length, and with an open reading frame of 975 bp encoding 324 amino acids. Sequence analysis demonstrated that nucleotide sequences of C. lanceolata CCR had 87% of genetic relationship with that of Picea abies one. Cluster analysis on CCR amino acids sequence of other plants showed that C. lanceolata had close relationship with P. abies, Pinus massoniana and P. taeda.Key words: lignin; cDNA cloning; Cunninghamia lanceolata; sequence analysis木质素(lignin)是一种由肉桂醇等单体聚合而成的酚类多聚体,是维管植物细胞壁的重要组成成分,在木材植物中含量较高,一般约占干重的15% ~ 35%,为地球上仅次于纤维素的有机组分[1]。
木质素的生物合成是在一系列酶催化下使苯丙氨酸或酪氨酸逐步转化为木质素单体,最终聚合成木质素的过程。
肉桂酰CoA还原酶(Cinnamyol Co-A Reductase,CCR)是影响木质素合成的一个关键酶,负责将各种羟基肉桂酰—辅酶A(CoA)酯最终还原成各种木质素单体,对木质素的合成与代谢也起关键作用。
木质素总含量和单体组成对植物体的最终用途起着决定性作用。
CCR基因除了影响木质素合成,还影响木材的微纤夹角(Microfibral Angle,mfa或MFA),进而导致材质密度和硬度的变异[2],因此对CCR基因进行克隆和研究具有重要的实际意义[3~4]。
杉木(Cunninghamia lanceolata)属杉科(Taxodiaceae)杉木属(Cunninghamia),是我国特有的常绿针叶6期陈喜,等:杉木木质素合成酶基因CCR的克隆和序列分析7树种之一。
杉木生长快,木材纹理直,结构细,耐腐力强,用途广泛,又因其萌芽力强,成林迅速,故成为重要的商品用材和造林树种。
建国后,我国加强了杉木科研工作,并取得了重要和丰硕成果[5~8],涉及杉木森林培育、经营、遗传育种、计测、木材加工利用等领域[9~11],其中杉木遗传改良研究最为活跃,研究成果及报道文献最多。
对杉木的CCR基因进行克隆,并与近缘属树种进行对比分析发现CCR基因变异情况,以期为杉木育种、遗传进化及转基因研究提供基础遗传学数据。
本文应用RACE(Rapid Amplification of cDNA Ends)等分子生物学方法,获得杉木木质素生物合成酶CCR基因的全长cDNA序列,并进行生物信息学分析,为CCR的进一步研究奠定基础。
1 材料与方法1.1 材料1.1.1 实验材料杉木材料采集于浙江林学院教学实习基地杉木种子园内。
样品采集后,直接于液氮速冻,带回实验室,用于mRNA的提取。
1.1.2 试剂 Taq DNA聚合酶、T4 DNA连接酶、DNA凝胶回收试剂盒、RNA逆转录试剂盒、克隆载体均购自TaKaRa公司;RACE试剂盒购自Clontech公司;其他试剂为进口或国产分析纯;宿主菌为E.coli DH5α,为本实验室保存。
测序由南京金思特公司完成。
1.2 方法1.2.1 杉木RNA的提取参考相关的文献,采用改良的CTAB法进行提取[12~16],具体方法如下:取1 g杉木茎和叶,液氮研磨后用药匙将粉末加入到65℃水浴的700μL提取液中,用手摇动混匀。
然后65℃水浴20 min后冷却至室温。
加等体积(700μL)的氯仿/异戊醇(24:1),混匀,10 000 r/min,4℃离心10 min。
小心吸取上清液,重复上述步骤,直至中间层无白色絮状沉淀。
吸取上清液,加1/4体积的10M LiCl混合,4℃沉淀过夜,然后10 000 r/min离心20 min。
弃上清液,用500μL SSTE溶解沉淀,加等体积氯仿/异戊醇混匀(24:1),10 000 r/min,4℃离心10 min。
离心吸取上清液,加2倍体积的无水乙醇,-70℃沉淀至少30 min,或-20℃下沉淀2 h。
12 000 r/min,4℃,离心20 min,去上清液,用75%的乙醇洗2次12 000 r/min,4℃离心5 min,小心弃上清液,吹干沉淀。
用40μL DEPC处理的水溶解沉淀,得到RNA样品。
1.2.2 CCR基因cDNA部分序列的克隆 cDNA合成:按照反转录TaKaRa PrimscriptTM RT-PCR Kit试剂盒(TaKaRa)的说明进行总RNA的反转录反应,获取cDNA。
1.2.2.1 RT-PCR扩增①引物设计。
根据GenBank中各种植物的CCR基因的核酸序列进行同源性比较,找出保守的区段,依据同源性高和简并性低的原则,设计一对简并引物:CCR-D-U 5’-GWASAAATGGTGGAGCCGC-3’CCR-D-L 5’-CCACGGTGGAGVACRCTCTC-3’引物由上海生工合成;②进行RT-PCR扩增,扩增体系为20μL体系,其中cDNA 2μL,10×PCR buffer 2μL,2.5 mM dNTP Mixture 0.8μL,5 U/μL TaKaRa Ex Taq HS 0.1μL,10μM Forward Primer 0.2μL,10μM Reverse Primer 0.2μL,ddH2O14.7μL。
扩增条件采用Touchdown PCR进行,反应程序:94℃预变性5 min;94℃变性30 s,54℃ 30 s,72℃延伸1 min,共14个循环,每个循环下降1℃;94℃变性30 s,54℃ 30s,72℃延伸1 min,30个循环;最后72℃延伸10 min。
RT-PCR产物经1%琼脂糖凝胶电泳检测,将回收产物进行目的片断的克隆和重组子筛选,并将目的基因片段测序鉴定。
1.2.3 杉木CCR基因cDNA片段3’端和5’端的扩增RACE采用BD SMARTTM RACE cDNA Amplification Kit(BD Biosciences Clontech),方法按照说明书再稍加改良。
①RACE引物设计及扩增。
反应利用Primer 5.0软件设计基因特异引物GSP1:5’>TGCTGGCATTGATGGCAGATGCAG<3’,用于3’-RACE;GSP2:5’> GGTACGCCGTGTGGTGTTTACATCT < 3’,用于5’-RACE。
按Clontech Smart RACE试剂盒说明操作,首先制备RACE-cDNA,然后进行RACE-PCR。
用1%琼脂糖凝胶电泳检测PCR结果。
②RACE产物的克隆测序和cDNA8 浙 江 林 业 科 技 30卷 全长获得。
3’RACE 产物、5’RACE 产物分别用pGEM-T 载体连接,经Amp 抗性筛选、X-gal/IPTG 蓝白斑显色筛选、重组质粒酶切鉴定后,送生物技术公司测序。
将3’RACE 产物测序结果和5’RACE 产物测序结果拼接为杉木CCR 基因全长cDNA 序列,将此cDNA 序列GenBank 数据库进行核苷酸同源性比较分析。
利用DNA-STAR 软件分析测序结果的序列与拼接结果的序列,最终将所得到的序列在NCBI 的Blast 进行比对,核苷酸序列的相似性比对在NCBI 的Blastn (/blastn )中分析,并确定开放阅读框分析(ORF ),将其推导相应的氨基酸序列。
蛋白质氨基酸相似性比对在NCBI 的Blastp 中分析;同样,运用Mega3.0软件对推导的氨基酸序列与已知报道植物推导的氨基酸序列进行系统进化分析。
2 结果与分析2.1 杉木次生木质部总RNA 的提取采用CTAB-LiCl 法提取杉木次生木质部总RNA ,经琼脂糖凝胶电泳检测,测得结果见图1。