青花菜P450CYP79F1全长克隆、表达及其与不同器官中莱菔硫烷含量的相关性分析

李占省，刘玉梅，方智远，杨丽梅，庄木，张扬勇，吕红豪



青花菜P450全长克隆、表达及其与不同器官中莱菔硫烷含量的相关性分析

李占省，刘玉梅，方智远，杨丽梅，庄木，张扬勇，吕红豪

（中国农业科学院蔬菜花卉研究所/农业部园艺作物生物学与种质创制重点实验室，北京 100081）

【目的】克隆青花菜细胞色素家族P450，分析其序列特征，阐明其在不同发育时期器官中的表达情况及其与莱菔硫烷含量的关系，为进一步揭示莱菔硫烷合成机理提供科学依据，为青花菜品种改良和新品种选育奠定基础。【方法】通过5′和3′端RACE克隆技术获得青花菜全长序列，利用DNAMAN 6.0进行基因拼接、氨基酸序列分析和蛋白二级结构预测，结合在线分析程序和软件进行基因和编码蛋白生物信息学分析；利用荧光定量PCR技术研究在不同发育器官中的表达情况；结合HPLC方法测定各相应时期器官中莱菔硫烷含量，包括青花菜发育时期的根、茎、叶，成熟花球，抽薹期花蕾（顶端蕾、成熟蕾、开花前1 d的蕾和花）和种子；对表达量与和莱菔硫烷含量进行相关性分析。【结果】获得了全长序列，GenBank登录号为MG012890，该基因全长2 014 bp，包含一个1 620 bp的开放读码框（ORF），编码540个氨基酸；编码的蛋白与甘蓝、白菜、芥蓝和油菜等同源蛋白的相似性均在95%以上；生物信息学分析表明该基因编码的酶蛋白有2个信号肽和2个跨膜结构域，为亲水性稳定蛋白，Wolf Psort预测其亚细胞定位于细胞质中；在根和茎中的表达量较高，叶中表达量最低，在花器官发育时期，自顶端蕾到花发育过程中呈现逐渐降低的表达趋势，种子中该基因表达量与花处于同一水平。相关性分析显示，青花菜发育时期的根、茎、叶、成熟花球、抽薹期花蕾和种子中莱菔硫烷含量与表达量无显著相关关系，但抽薹期花蕾自顶端花蕾到花中莱菔硫烷生成量与表达量呈显著正相关关系（=0.96，<0.05），在调控莱菔硫烷的生成方面起重要作用，尤其是在生殖器官花蕾的发育过程中，能够显著上调莱菔硫烷的生成量。【结论】在青花菜不同发育器官中对莱菔硫烷的生成发挥着重要调控作用，可能与青花菜不同器官中莱菔硫烷含量多样性密切相关。

青花菜；莱菔硫烷；；表达分析；HPLC

0 引言

【研究意义】青花菜（var.）又名西兰花、绿菜花等，十字花科芸薹属甘蓝类蔬菜变种，一、二年生植物。青花菜营养全面，富含多种矿物质和维生素，同时富含4-甲酰基亚磺酰基丁基硫甙（glucoraphanin，GRA），4-甲酰基亚磺酰基丁基硫甙能被位于液泡中的黑芥子酶水解生成异硫氰酸盐类物质——莱菔硫烷（sulforaphane，SF或SFN）。流行病和临床医学研究表明，莱菔硫烷能够显著降低肝癌、胃癌、肺癌、膀胱癌、乳腺癌和结肠直肠癌等多种癌症的患病率[1-7]，在降低糖尿病、高血压和心脑血管疾病等发病率方面具有很好的医学前景[8-11]。据报道，青花菜不同器官中莱菔硫烷含量差异较大，不同基因型材料也存在显著差异，种子中莱菔硫烷含量最高，芽苗在前期（5—7 d）含量较高，后期降低，花球中莱菔硫烷含量高于幼茎，叶片最低。另有研究表明，青花菜叶片中莱菔硫烷含量处于较高水平，这可能由不同发育时期的样品和基因型造成[12-18]。拟南芥研究结果表明，为莱菔硫烷代谢的上游调控基因，在莱菔硫烷的合成方面可能具有重要作用[19-22]。因此克隆青花菜该基因全长，分析其在青花菜不同发育时期各器官中的表达量与莱菔硫烷含量的相关性，将为进一步阐明青花菜莱菔硫烷含量分布的器官多样性提供科学依据，对青花菜不同器官的科学利用和新品种选育都具有重要意义。【前人研究进展】拟南芥硫甙代谢研究表明，能够将长链或短链甲硫氨酸催化形成相应的醛肟，尤其是短链脂肪酸（1C-4C）的合成尤为重要[19,21]。莱菔硫烷的前体为4-甲基亚磺酰基丁基硫甙，为短链脂肪酸合成途径，因此，在调控4-甲基亚磺酰基丁基硫甙的生成方面尤为重要。另外，过表达能够调控拟南芥株高（变矮）和植株大小，且与根部发育也具有直接的调控作用[23-25]。王红波等[26]通过转基因技术研究指出，的过表达能够显著增加4-甲基亚磺酰基丁基硫甙含量，吲哚组硫甙成分则相对降低，暗示该基因与脂肪族硫甙和莱菔硫烷代谢密切相关。过表达能够增加脂肪族硫甙成分在拟南芥上同样得到了证实[25,27-28]。【本研究切入点】青花菜富含抗癌活性成分莱菔硫烷已逐渐被大众认知，青花菜不同器官中莱菔硫烷含量多样性分布的原因目前尚未见报道。【拟解决的关键问题】从青花菜中克隆获得全长，分析其家族序列特征，研究其在青花菜不同发育时期各器官中的表达特性与莱菔硫烷含量相关性，为深入研究青花菜不同器官中莱菔硫烷代谢和转运机理奠定基础。

1 材料与方法

试验于2015年8月至2017年7月在中国农业科学院蔬菜花卉研究所进行。

1.1 试验材料

青花菜材料B691为纯合自交系材料（自交16代），由中国农业科学院蔬菜花卉研究所甘蓝青花菜课题组选育。

试验材料于2015年8月中旬定植于露地，每份材料设置3次重复，每重复种植11株，行间距为45 cm×55 cm。2015年10月初开始对试验材料的根、茎、叶、发育时期的花球进行取样，待植株移栽到温室后，对抽薹期的花蕾（顶端蕾、成熟蕾、开花前1 d蕾和花）进行取样。相同材料保存2份，一份于超低温保存后真空冷冻干燥，机械粉碎后密封保存；另一份提取总RNA进行定量表达分析。

1.2 青花菜CYP79F1克隆

青花菜总RNA的提取采用Trizol法，RNA提取试剂盒和反转录试剂盒（MMLV反转录试剂盒）均由日本TaKaRa公司购置。按照TaKaRa试剂盒要求除去痕量DNA的总RNA，进行cDNA第一条链的合成。5′-RACE cDNA第一链的合成：根据Invitrogen公司Gene- RacerTM试剂盒要求进行总RNA的脱磷、脱帽处理，然后在5′端加上RNA Oligo，利用SuperScript反转录酶Ⅲ，以基因特异性引物进行反转录（表1）。3′-RACE cDNA第一链合成仍参照Gene-RacerTM试剂盒要求进行，所需引物信息见表1。

cDNA核心克隆则参照报道的已知序列（Genebank：KP693683）设计核心引物（表1）。反应体系为：Taq 2X Master Mix（百泰克）25 μL，正向引物（10 μmol∙L-1）1 μL，反向引物（10 μmol∙L-1）1 μL，cDNA模板（80—100 ng）1 μL，双蒸水22 μL，共50 μL。PCR扩增条件为：95℃变性5 min；95℃变性30 s，58℃退火30 s，72℃延伸45 s，35个循环；72℃延伸10 min。从反应产物中取10 μL在1%的琼脂糖凝胶上进行电泳分析。

将所有扩增产物进行1.5%TAE琼脂糖凝胶电泳检测后，使用1%琼脂糖凝胶电泳，回收PCR产物并与Blunt Zero Cloning载体连接，构建重组质粒转化DH5感受态细胞，筛选阳性克隆，送至华大生物技术有限公司测序，将测序结果与保守区序列进行比对拼接。

1.3 生物信息学分析

利用DNAMAN 6.0进行基因拼接、氨基酸序列分析和蛋白二级结构预测，在线NCBI进行序列处理（https://www.ncbi.nlm.nih.gov/tools/vecscreen/，Vecscreen）与分析，Blastx和Blastp分析，蛋白质一级结构分析（http://web.expasy.org/protparam/，ProtParam）和修饰分析（http://myhits.isb-sib.ch/ cgibin/motif_scan，MotifScan），跨膜预测分析（http:// www.cbs.dtu.dk/services/TMHMM-2.0/，TMHMM），系统进化树构建在线分析（NCBI），亚细胞定位在线分析（https:// omictools.com/wolf-psort-tool，WoLF PSORT）。

1.4 青花菜CYP79F1表达分析

利用软件Primer Premier 5.0软件设计扩增引物，引物序列为双向引物为F（5′-3′）：GTCACGCCAGACG AAATCAAA和R（5′-3′）：GCACAAGCCTGTCTTTT CCAACT，目标长度为169 bp，内参基因为-12，引物序列F（5′-3′）：GGCTCTATCTTGGCTTCTCTCA GT和（5′-3′）：CCAGATTCATCATACTCGGCTTT。参照TaKaRa SYBR® Green PCR Master Mix说明进行操作，在ABI 7500实时定量PCR仪上进行PCR扩增。反应体系为50 μL：25 μL SYBR® Premix，5 μL cDNA模板，双向引物各1 μL，1 μL Rox，加蒸馏水至50 μL。PCR反应程序为：95℃预变性30 s，95℃变性5 s，58℃退火30 s，40个循环，每样品各设置3次重复。反应完成后进行融解曲线和荧光值变化曲线分析，采用2-ΔΔCt法计算目的基因的相对表达情况。

表1 CYP79F1扩增引物信息

采用实时定量PCR时各样品加样量均为2 µL，然而由于受RNA浓度定量误差和RNA逆转录效率误差等影响，每个样品2 µL的cDNA含量并不完全相同，为校正此差异，使用内参基因（不同样品间表达量基本恒定）进行校正。根据Real-time PCR原始检测结果，按照2-△△ct相对定量法计算出各样品的目的基因相对定量结果[16]。

1.5 莱菔硫烷含量的提取与HPLC分析

青花菜不同器官中莱菔硫烷含量的提取与HPLC分析参照李占省等[17,29]的方法进行。

1.6 相关性分析与数据处理

利用统计学方法，将青花菜不同发育时期各器官中基因的表达量与莱菔硫烷含量进行相关性分析，具体采用相关性分析。试验数据采用Excel 2003进行基本处理，方差分析和卡方检验采用SPSS 17.0软件，试验结果用平均值±标准误（mean±SD，n=3）表示，one-way ANOVA用于显著性检验，多重比较采用Duncan’s新复极差法分析（<0.05）。

2 结果

2.1 青花菜CYP79F1全长克隆

根据PCR产物5′-RACE和3′-RACE扩增产物测序并拼接，长度为2 014 bp的序列（图1），前边包含62个非编码区，从63处开始翻译，到1 685 bp处终止，将拼接后的序列在线Blastx分析（NCBI），结果如图2所示，其中的小写字母通常被设定为没有三维结构的氨基酸。核心cDNA序列包含一个完整的1 620 bp的开放阅读框（open reading flame，ORF），5′端有62 bp的非翻译区（untranslate region，UTR），富含A和T，A（597）和T（541）的总含量达56.5%，G（459）和C（417）含量占总含量的43.5%。经DNAMAN 6.0分析和NCBI在线Blastx分析（图2），该基因编码区编码一条540个氨基酸残基的多肽，分子量61.434 kD，等电点8.29，pH 7.0时的带电荷数（ch）为6.63，3′端具有完整的Poly A尾巴。该基因检测到2个甲基化位点，DAM甲基化酶对应位点为GATC 2-A，DCM甲基化酶对应位点为CCWGG 2-C。该基因序列已在GenBank注册，登记号为MG012890。

2.2 CYP79F1编码蛋白质的特性分析

该基因编码区发现22个抗原肽段，其分值和长度变化分别为1.03—1.19和6—31个氨基酸。

1—4代表3端′端RANCE克隆，6—7为核心区克隆，8—9为5′端RANCE克隆，M为5000 bp

图2 CYP79F1与甘蓝型油菜blastx比对分析结果

经氨基酸特性分析，构成蛋白质中共存在20种氨基酸，赖氨酸（Leu）含量最高（10.3%），其次是精氨酸（Arg）和谷氨酸（Glu），分别占6.9%，然后依次是甘氨酸（Gly）（6.5%）、异亮氨酸（Ile）和丙氨酸（Ala）（6.3%），色氨酸（Trp）含量最低，为1.5%。就构成氨基酸的特性而言，该蛋白包含82个强碱性氨基酸（H、K和R），66个强酸性氨基酸（D和E），240个疏水氨基酸（A、F、I、L、M、P、V和W）和152个极性氨基酸（C、G、N、Q、S、T和Y），预测该蛋白为亲水性蛋白（亲水平均值-0.11，疏水平均值-0.26）。

Blastp搜索显示，CYP79F1属于细胞色素家族P450超级家族，编码540个氨基酸。Motif Scan 分析显示，该蛋白存在细胞色素P450C半胱氨酸血红素-iron配体信号（Cytochrome P450 cysteine heme-iron ligand signature），位于第469—478氨基酸处。PROSITE数据库分析表明，该蛋白存在2处酰胺化位点（AMIDATION Amidation site），207—210和471—474；6处酪蛋白激酶II磷酸化位点（CK2_ PHOSPHO_SITE Casein kinase II phosphorylation site），246—249、307—310、318—321、344—347、458—461和509—512；4处N-豆蔻酰化位点（MYRISTYL N-myristoylation site），129—134、293—298、428—433和478—483。此外发现，6处蛋白激酶C磷酸化位点（PKC_PHOSPHO_SITE Protein kinase C），156—158、202—204、307—309、390—392、432—434和471—473；1处酪氨酸激酶磷酸化位点（TYR_PHOSPHO_SITE Tyrosine kinase phosphorylation site），434—442。

TMHMM显示，该蛋白为跨膜蛋白，存在2处跨膜区，第一个位于编码氨基酸7—35处，包含29个氨基酸；第二个位于编码氨基酸472—500处，包含29个氨基酸，其结构详见图3。该蛋白二级结构组成为：45%的H（alpha-helix）、6%的E（beta-sheet）和47%的C（loop或coil）。利用SWISSMODEL对青花菜CYP79F1蛋白的三级结构进行预测，详见图4。亚细胞定位预测位于细胞质中。

2.3 CYP79F1 蛋白的同源比对分析

同源分析表明，该CYP79F1蛋白与十字花科作物中的甘蓝型油菜（）、野生甘蓝（）、芥蓝（var）、白菜（）和萝卜（）等氨基酸序列高度同源，相似性分别为99%、99%、99%、97%和95%，其次是十字花科的菘蓝（），又名北板蓝根，相似性为94%，这为分析两种基于该基因代谢中的相似成分与调控机理提供了依据[30]。

2.4 CYP79F1的表达分析

定量表达分析显示，在青花菜的根和茎中表达量最高，其次是顶端蕾、开花前1 d蕾、成蕾、花、叶和花球，种子中该基因的表达量与抽薹期花中的表达量处于同一水平（图5）。从花蕾发育期可以看出，该材料自顶端蕾到花的发育过程中基本呈现逐渐降低的表达趋势，该结果已经多次验证并报道[16,31]，暗示与花器官发育和激素代谢密切相关，其调控机理需要深入研究。

2.5 青花菜中莱菔硫烷含量HPLC分析

HPLC分析表明，青花菜不同发育时期各器官中莱菔硫烷含量差异显著（<0.05）（图6）。结果表明，种子中莱菔硫烷含量最高，平均含量为4 855.64 mg∙kg-1DW，与当前报道一致[14-18]。花器官发育时期总体含量处于较高水平，平均值为1 027.35 mg∙kg-1DW，含量最高的为顶端花蕾（1 327 mg∙kg-1DW），然后依次是成蕾（1 001.81 mg∙kg-1DW）、开花前1 d蕾（994.39 mg∙kg-1DW）和花（786.20 mg∙kg-1DW），其余器官中，花球（456.45 mg∙kg-1DW）高于幼茎（282.86 mg∙kg-1DW），叶片（160.46 mg∙kg-1DW）中莱菔硫烷含量处于较低水平，根中莱菔硫烷含量最低（17.73 mg∙kg-1DW）。同时，莱菔硫烷含量检出值在花器官发育时期呈现出一定的变化规律，表现为顶端花蕾发育到花的过程中呈逐渐降低的趋势，顶端蕾中含量最高，是根中莱菔硫烷检出值的75倍，暗示莱菔硫烷含量合成与植物分化可能存在一定的关系。

图3 CYP79F1蛋白跨膜结构预测

A和B分别代表不同角度的3D图 Different points of 3D view shown by A and B

不同小写字母表示差异显著（P＜0.05）。下同 Different lowercases indicated significant differences at P＜0.05 level. The same as below

图6 青花菜不同器官中莱菔硫烷含量HPLC分析

2.6 青花菜不同发育器官中莱菔硫烷含量与CYP79F1表达量相关性分析

相关性分析表明，表达量与根、茎、叶、花球、顶端蕾、成蕾、开花前1 d蕾、花和种子中莱菔硫烷生成量间的相关系数为-0.273，负相关性不显著（=0.478）（表2）。为了明确顶端花蕾发育到花的过程中莱菔硫烷含量逐渐降低是否与表达量具有显著相关性，将该发育时期的莱菔硫烷含量与对应器官的基因表达量进行相关性分析。结果发现，两者呈显著正相关关系，相关系数为0.959（=0.041）（表3）。表明在青花菜花器官发育中与莱菔硫烷的积累具有明显相关性，可能发挥了直接调控作用。

表2 不同器官中莱菔硫烷含量与基因表达量相关性分析

表3 抽薹期花蕾莱菔硫烷含量与基因表达量相关性分析

*表示显著相关（＜0.05）

* indicate significant correlation (＜0.05)

3 讨论

在调控脂肪族硫甙合成方面发挥着重要作用，与青花菜中莱菔硫烷生成量密切相关，同时在植物发育过程中起着重要的调控作用[14,20-21,23]。目前，拟南芥中已经克隆该基因（AT1G16410），该基因属于细胞色素家族P450成员，在基因组中的长度为3 383 bp，其cDNA为1 200 bp。油菜中该基因位于基因组scaffold_526区域内，该scaffold DNA全长为436 263 bp。青花菜基因组不同于拟南芥和油菜，存在多倍化和重组事件，同时不同变种之间往往存在基因变异的可能，而青花菜中有关该基因与莱菔硫烷方面的研究鲜有报道。因此，本研究通过RACE技术，从青花菜中获得mRNA全长为2 014 bp，同源比对结果表明，该序列与拟南芥、油菜基因相似性分别为86%和99%，与油菜亲缘关系更近，表明该基因在芸薹属异源四倍体基因组中能够稳定遗传，暗示了其在进化和重组中的保守性。经在线分析（NCBI），与该基因同源的还有芥蓝和野生甘蓝，其相似性均为99%，再次验证了该克隆的准确性和保守性[32]。

目前，功能已在拟南芥和部分芸薹属作物中获得验证[19,21]。本研究预测青花菜中亚细胞定位于细胞质中，与当前报道一致，同时，编码蛋白为亲水性跨膜蛋白，存在2处跨膜区，各含有3个P450信号序列，含有2个信号肽序列，这与拟南芥上有关该基因在内质网上发挥功能的定位与结论一致[21,33]。此外，本研究揭示了该编码蛋白存在2处酰胺化位点，6处酪蛋白激酶II磷酸化位点，4处N-豆蔻酰化位点，6处蛋白激酶C磷酸化位点和1处酪氨酸激酶磷酸化位点，这方面研究尚未见报道，该发现为深入研究CYP79F1高级结构和功能修饰位点提供了依据和支持[34]。

本研究发现，在青花菜不同发育器官中存在明显的表达特异性，莱菔硫烷的生成量在不同器官存在多样性。在根中和茎中表达量较高，但相应器官中莱菔硫烷含量并不是处于较高水平，种子中莱菔硫烷含量处于最高水平，但其表达量并不高。相关性分析结果表明，青花菜发育不同发育时期的器官中，表达量与莱菔硫烷的生成量没有显著相关性，但在抽薹期发育的花蕾中则呈现出显著的正相关关系，暗示该基因在花器官发育过程中可能直接参与了莱菔硫烷的代谢调控，该发现也为少数报道该基因间接参与植物激素代谢和细胞分化途径提供了参考[35-37]。

由此可见，青花菜发育时期不同器官中莱菔硫烷含量多样性是一个复杂的调控网络和运输机制。本研究发现，在生殖器官中总体表达较为活跃，而在营养器官种呈现出不同的表达水平，营养器官中莱菔硫烷HPLC检测值明显低于生殖器官（花器官和种子），暗示在调控莱菔硫烷的生成量方面发挥了重要调控作用，但营养器官和生殖器官间是否存在莱菔硫烷的主动运输关系，还需要进一步验证。

此外，为脂肪族硫甙上游调控基因，尤其是能够直接调控3C和4C脂肪族硫甙生成，而对应的下游硫甙包括GRA、PRO、NAP等，这些硫甙成分与青花菜等十字花科蔬菜（作物）的风味、植保和营养密切相关，因此该基因可能在十字花科芸薹属作物中发挥了重要的调控作用[38-41]。本研究结论为深入研究与芸薹属作物硫甙成分与风味奠定了基础。与细胞分化和细胞增殖密切相关，这与拟南芥中报道的该基因能够调控拟南芥株高和大小的结论相一致[36-37]，也暗示了可能参与了多个代谢途径的调控。

4 结论

利用RACE技术从青花菜中获得了全长（2 014 bp），ORF为1 620 bp，编码540个氨基酸。该基因属于P450基因超家族，各含有3个P450信号序列，含有2个信号肽序列，属于亲水性跨膜蛋白，亚细胞定位于细胞质中，其编码的蛋白与油菜、芥蓝和野生甘蓝等芸薹属蔬菜具有高度同源性。在青花菜不同发育器官中存在明显的表达特异性，在营养器官根和茎中表达量最高，其次是花器官发育时期，总体处于较高水平，与花器官发育呈负相关关系。结合各时期发育器官中莱菔硫烷含量变化，推测该基因通过组织表达特异性，可能直接参与调控青花菜不同发育器官中莱菔硫烷的生成和积累。

[1] Abdull Razis A F, Iori R, Ioannides C. The natural chemopreventive phytochemical R-sulforaphane is a far more potent inducer of the carcinogen-detoxifying enzyme systems in rat liver and lung than the S-isomer., 2011, 128(12): 2775-2782.

[2] Fahey J W, Haristoy X, Dolan P M, Kensler T W, Scholtus I, Stephenson K K, Talalay P, Lozniewski A. Sulforaphane inhibits extracellular, intracellular, and antibiotic- resistant strains of Helicobacter pylori and prevents benzo[a]pyrene- induced stomach tumors., 2002, 99(11): 7610-7615.

[3] Berger M S, Cooley M E, Abrahm J L. A pain syndrome-associated with large adrenal metastases in patients with lung-cancer., 1995, 10(2): 161-166.

[4] Abbaoui B, Riedl K M, Ralston R A, Thomas-Ahner J M, Schwartz S J, Clinton S K, Mortazavi A. Inhibition of bladder cancer by broccoli isothiocyanates sulforaphane and erucin: Characterization, metabolism, and interconversion., 2012, 56(11): 1675-1687.

[5] Dang Y M, Huang G, Chen Y R, Dang Z F, Chen C, Liu F L, Guo Y F, Xie X D. Sulforaphane inhibits the proliferation of the BIU87 bladder cancer cell line via IGFBP-3 elevation., 2014, 15(4): 1517-1520.

[6] Agyeman A S, Chaerkady R, Shaw P G, Davidson N E, Visvanathan K, Pandey A, Kensler T W. Transcriptomic and proteomic profiling of KEAP1 disrupted and sulforaphane-treated human breast epithelial cells reveals common expression profiles., 2012, 132(1): 175-187.

[7] Chung F L, Conaway C C, Rao C V, Reddy B S. Chemoprevention of colonic aberrant crypt foci in Fischer rats by sulforaphane and phenethyl isothiocyanate., 2000, 21(12): 2287-2291.

[8] Jiang X, Bai Y, Zhang Z, Xin Y, Cai L. Protection by sulforaphane from type 1 diabetes-induced testicular apoptosis is associated with the up-regulation of Nrf2 expression and function., 2014, 279(2): 198-210.

[9] Miao X, Bai Y, Sun W, Cui W, Xin Y, Wang Y, Tan Y, Miao L, Fu Y, Su G. Sulforaphane prevention of diabetes- induced aortic damage was associated with the up-regulation of Nrf2 and its down-stream antioxidants., 2012, 9(1): 84.

[10] Li Z, Galli U, Becker L E, Bruns H, Nickkolgh A, Hoffmann K, Karck M, Schemmer P. Sulforaphane protects hearts from early injury after experimental transplantation., 2013, 18: 558-566.

[11] Piao C S, Gao S, Lee G H, Kim do S, Park B H, Chae S W, Chae H J, Kim S H. Sulforaphane protects ischemic injury of hearts through antioxidant pathway and mitochondrial K(ATP) channels., 2010, 61(4): 342-348.

[12] Brown A F, Yousef G G, Jeffrey E H, Klein B P, Wallig M A, Kushad M M, Juvik J A. Glucosinolate profiles in broccoli: Variation in levels and implications in breeding for cancer chemoprotection., 2002, 127(5): 807-813.

[13] Farnham M W, Wilson P E, Stephenson K K, Fahey J W. Genetic and environmental effects on glucosinolate content and chemoprotective potency of broccoli., 2004, 123(1): 60-65.

[14] Fahey J W, Holtzclaw W D, Wehage S L, Wade K L, Stephenson K K, Talalay P. Sulforaphane bioavailability from glucoraphanin-rich broccoli: Control by active endogenous myrosinase., 2015, 10(11): e0140963.

[15] Liang H, Yuan Q P, Dong H R, Liu Y M. Determination of sulforaphane in broccoli and cabbage by high-performance liquid chromatography., 2006, 19(5): 473-476.

[16] Li Z S, Liu Y M, Fang Z Y, Yang L M, Zhuang M, Zhang Y Y, Zhao W, Sun P T. Variation of sulforaphane levels in Broccoli (var.) during flower development and the role of gene Aop2., 2014, 37(9): 1199-1211.

[17] 李占省, 刘玉梅, 方智远, 杨丽梅, 庄木, 张扬勇, 吕红豪, 孙培田. 甘蓝和青花菜不同器官中莱菔硫烷含量差异研究. 核农学报, 2017, 31(3): 447-454.

LI Z S, LIU Y M, FANG Z Y, YANG L M, ZHUANG M, ZHANG Y Y, LÜ H H, SUN P T. Study on difference of sulforaphane in cabbage head and different organs of broccoli., 2017, 31(3): 447-454. (in Chinese)

[18] 姚雪琴, 谢祝捷, 李光庆, 邱海荣, 李媛. 青花菜不同器官中4-甲基亚磺酰丁基硫苷及萝卜硫素含量分析. 中国农业科学, 2011, 44(4): 851-858.

YAO X Q,XIE Z J,LI G Q,QIU H R,LI Y. Analysis of glucoraphanin and sulforaphane contents in different organs of broccoli., 2011, 44(4): 851-858 . (in Chinese)

[19] Sonderby I E, Geu-Flores F, Halkier B A. Biosynthesis of glucosinolates-gene discovery and beyond., 2010, 15(5): 283-290.

[20] Falk K L, Vogel C, Textor S, Bartram S, Hick A, Pickett JA, Gershenzon J. Glucosinolate biosynthesis: demonstration and characterization of the condensing enzyme of the chain elongation cycle in., 2004, 65(8): 1073-1084.

[21] Grubb C D, Abel S. Glucosinolate metabolism and its control., 2006, 11(2): 89-100.

[22] Wu S H, Lei J J, Chen G J, Chen H C, Cao B H, Chen C M. De novo Transcriptome assembly of Chinese kale and global expression analysis of genes involved in glucosinolate metabolism in multiple tissue., 2017, 8(27): 1-16.

[23] Chen S X, Glawischnig E, Jorgensen K, Naur P, Jorgensen B, Olsen C E, Hansen C H, Rasmussen H, Pickett J A, Halkier B A. CYP79F1 and CYP79F2 have distinct functions in the biosynthesis of aliphatic glucosinolates in., 2003, 33(5): 923-937.

[24] Reintanz B, Lehnen M, Reichelt M, Gershenzon J, Kowalczyk M, Sandberg G, Godde M, Uhl R, Palme K. Bus, a bushy arabidopsis CYP79F1 knockout mutant with abolished synthesis of short-chain aliphatic glucosinolates., 2001, 13(2): 351-367.

[25] Hansen C H, Wittstock U, Olsen C E, Hick A J, Pickett J A, Halkier B A. Cytochrome P450CYP79F1 from Arabidopsis catalyzes the conversion of dihomomethionine and trihomomethionine to the corresponding aldoximes in the biosynthesis of aliphatic glucosinolates., 2001, 276(14): 11078-11085.

[26] 王红波, 李成, 曹阳理惠, 李晶. CYP79F1基因对西兰花的遗传转化的研究. 植物研究, 2014, 34(4): 516-523.

WANG H B, LI C, CAO Y L H, LI J. Genetic Transformation of CYP79F1 Gene invar., 2014, 34(4): 516-523. (in Chinese)

[27] Brown P D, Tokuhisa J G, Reichelt M, Gershenzon J. Variation of glucosinolate accumulation among different organs and developmental stages of., 2003, 62(3): 471-481.

[28] Sharma M, Mukhopadhyay A, Gupta V, Pental D, Pradhan A K. Bju B. CYP79F1 Regulates synthesis of propyl fraction of aliphatic glucosinolates in oilseed mustard: Functional validation through genetic and transgenic approaches., 2016, 11(2): e0150060.

[29] 李占省, 刘玉梅, 方智远, 杨丽梅, 庄木, 张扬勇, 袁素霞, 赵文, 刘二艳, 孙培田. 青花菜DH群体花球中莱菔硫烷含量的遗传效应分析. 园艺学报, 2012, 39(1): 101-108.

Li Z S, Liu Y M, Fang Z Y, Yang L M, Zhuang M, Zhang Y Y, Yuan S X, Zhao W, Liu E Y, Sun P T. Determination of sulforaphane by high performance liquid chromatography and genetic analysis of DH population in broccoli florets., 2012, 39(1): 101-108. (in Chinese)

[30] Tantikanjana T, Mikkelsen M D, Hussain M, Halkier B A, Sundaresan V. Functional analysis of the tandem-duplicated P450 genes SPS/BUS/CYP79F1 and CYP79F2 in glucosinolate biosynthesis and plant development by Ds transposition-generated double mutants., 2004, 135(2): 840-848.

[31] Li Z S, Liu Y M, Fang Z Y, Yang L M, Zhuang M, Zhang Y Y, LV H H. Development and identification of anti-cancer component of sulforaphane in developmental stages of broccoli (var.L.)., 2016, 4(8): 490-497.

[32] Guo L P, Yang R Q, Gu Z X. Cloning of genes related to aliphatic glucosinolate metabolism and the mechanism of sulforaphane accumulation in broccoli sprouts under jasmonic acid treatment., 2016, 96(13): 4329-4336.

[33] Yan X F, Chen S X. Regulation of plant glucosinolate metabolism., 2007, 226(6): 1343-1352.

[34] Kopriva S, Gigolashvili T. Glucosinolate synthesis in the context of plant metabolism., 2016, 80: 99-124.

[35] Luck K, Jia Q D, Huber M, Handrick V, Wong G K S, Nelson D R, Chen F, Gershenzon J, Kollner T G. CYP79 P450 monooxygenases in gymnosperms: CYP79A118 is associated with the formation of taxiphyllin in Taxus baccata., 2017, 95(1-2): 169-180.

[36] Chaban C, Waller F, Furuya M, Nick P. Auxin responsiveness of a novel cytochrome p450 in rice coleoptiles., 2003, 133(4): 2000-2009.

[37] Bak S, Feyereisen R. The involvement of two P450 enzymes, CYP83B1 and CYP83A1, in auxin homeostasis and glucosinolate biosynthesis., 2001, 127(1): 108-118.

[38] Baik H Y, Juvik J, Jeffery E H, Wallig M A, Kushad M, Klein B P. Relating glucosinolate content and flavor of broccoli cultivars., 2003, 68(3): 1043-1050.

[39] Rohr F, Ulrichs C, Mewis I. Variability of aliphatic glucosinolates in(L.) - Impact on glucosinolate profile and insect resistance., 2009, 82(2): 131-135.

[40] Rohr F, Ulrichs C, Schreiner M, Nguyen C N, Mewis I. Effects of two folivorous insect species on the aliphatic glucosinolate profile inecotype crosses., 2012, 18: 219.

[41] Wittstock U, Kliebenstein D J, Lambrix V, Reichelt M, Gershenzon J. Glucosinolate hydrolysis and its impact on generalist and specialist insect herbivores., 2003, 37: 101-125.

（责任编辑 赵伶俐）

Full Length Cloning, Expression and Correlation Analysis of P450Gene with Sulforaphane Content in Different Broccoli Organs

LI ZhanSheng, LIU YuMei, FANG ZhiYuan, YANG LiMei, ZHUANG Mu, ZHANG YangYong, LÜ HongHao

(Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences/Key Laboratory of Biological and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, Beijing 100081)

【Objective】The aim of this study was to clone the full length of P450 familygene from broccoli, investigate its sequence characteristic and analyze its expression and correlation with sulforaphane content in different developmental broccoli organs. These findings will provide a scientific basis for further revealing the mechanism of sulforaphane metabolism in broccoli and lay a practical foundation for the genetic improvement and breeding of broccoli varieties. 【Method】 The full length sequence ofgene was obtained by RACE cloning (5′ and 3′). According to the software of DNAMAN 6.0, gene splicing, amino acid sequence analysis and protein secondary structure prediction ofgene, at the same time, bioinformatics information was also predicted by online program and software. The quantitative real-time PCR (qRT-PCR) method was used to study the expression ofgene in different developmental organs of broccoli. The sulforaphane contents of different developmental organs were detected by high performance liquid chromatography (HPLC) method, and they were roots, stems, leaves, developmental buds at bolting stage (top buds, mature buds, buds one day before flowering and flowers), and seeds.analysis was carried out by the software SPSS 10.0.【Result】Full length sequence ofgenewere cloned, and the gene ofwas 2 014 bp in length, containing a 1 620 bp opening reading frame (ORF), encoded a polypeptide of 540 amino acids. They shared over 95% identity with the homologous proteins fromplants, such as cabbage, Chinese cabbage, Chinese kale and rapeThe bioinformatics analysis showed that the CYP79F1 protein was hydrophilic protein with two signal-peptides and two transmembrane regions, and the Wolf Psort protection indicated that they were located in the cytoplasm. The expression ofgene was higher in the roots and the stems, and the lowest in the leaves. During developmental stage of buds, the expression ofgenewas up-regulated in the early (young buds), and then decreased (mature buds to flowers). There was no difference significance in gene expression between roots, stems, leaves, developmental buds at bolting stage and seeds. Correlation analysis showed that the expression level ofgeneshowed an extremely significant positive correlation with sulforaphane content in developmental buds at bolting stage (=0.96，<0.05), and the expression ofgenemight influence the generation of sulforahane in different organs of broccoli, especially in developmental buds.【Conclusion】Thegenewere comprehensively obtained and characterized from broccoli,mayplay an important role in sulforaphane metabolism and revealing the diversity of sulforaphane content in different organs of broccoli.

broccoli; sulforaphane;;expression analysis; HPLC

2017-10-24；

2018-01-15

国家自然科学基金（31372067，31501761）、国家大宗蔬菜产业技术体系（CARS-23-A）、中国农业科学院科技创新工程（CAAS-ASTIP-IVFCAAS）

李占省，E-mail：lizhansheng@caas.cn。

刘玉梅，E-mail：liuyumei@caas.cn

doi: 10.3864/j.issn.0578-1752.2018.12.011