摘要：为了解猕猴桃查尔酮合成酶基因（Chalcone Synthase，CHS）的密码子使用特征，采用CodonW、SPSS、MEGA软件和EMBOSS在线程序等分析猕猴桃CHS密码子偏好性，并构建单双子叶植物CHS系统发育树。研究结果表明，猕猴桃CHS有效密码子数（ENc）、总G和C含量（GC）和密码子第3位上G或C含量（GC3s）分别为5721、0533和0607；同义密码子相对使用度（RSCU）大于10的密码子数为30，但不存在偏好性极强的密码子；聚类分析结果发现，基于密码子使用偏好性分类可将单双子叶进行准确归类，而CHS序列相似性聚类结果并不理想；密码子使用频率分析表明，猕猴桃CHS与模式生物拟南芥和大肠杆菌基因组间密码子使用偏好性差异相对较小。研究认为，猕猴桃CHS密码子偏好性较弱，但倾向使用G和C，且以G或C结尾的密码子；单双子叶植物CHS间存在严格的密码子使用法则；此外，拟南芥和大肠杆菌可能是猕猴桃CHS遗传转化和异源表达较为理想的受体系统。
关键词：猕猴桃；查尔酮合成酶；密码子使用偏好性；聚类分析
中图分类号：S 6634文献标识码：A文章编号：1008-0384（2017）11-1213-05
Abstract：To investage codon usage characteristic of chalcone synthase gene of kiwifruit， CodonW， SPSS， MEGA software and EMBOSS online program were used to analyse codon usgae bias （CUB） of CHS， and phylogenetic tree of CHS between monocotyledon and dicotyledon was constructed in this paper The results showed that， the effective number of codons （ENc）， total GC content （GC） and GC content on the 3rd site （GC3s） was 5721， 0533 and 0607，respectively The number of condons，the number of relative synonymous codon usage （RSCU） of which was more than 10，was 30， but no codons was with strongly CUB The cluster tree of CHS based on CUB was closer to the real category system than it did by sequence similarity According to codon usage frequency， the CUB between CHS of kiwifruit and Arabidopsis thaliana， Escherichia coli were more similar These results indicated that， the CUB level of CHS was low in kiwifruit， but it was bias toward GC rich codons and synonymous codons with G or C at 3rd codon position And CHS in monocotyledons and dicotyledons had the distinctive CUB rules In addition， Arabidopsis thaliana and Escherichia coli migth be the best receptors for genetic transformation and heterologous expression of CHS in kiwifruit.
Key words：Kiwifruit； chalcone synthase； codon usage bias； cluster analysis
獼猴桃是猕猴桃科猕猴桃属多年生藤本果树，其果实中富含酚类化合物、微量元素、氨基酸和维生素C等多种营养和保健成分。猕猴桃栽培主要包括红肉、黄肉和绿肉等品种，尤其红肉猕猴桃因其果肉颜色鲜艳，风味优良而广受欢迎。花青素的积累是果实着色的重要环节。研究发现，红肉猕猴桃内果皮中的红色物质也由花青素构成[1]。目前，植物花青素的分子研究较为广泛，花青素的合成和转运参与的功能基因较多，主要包括查尔酮合成酶CHS[2]、查尔酮异构酶CHI[3]、黄烷酮羟基化酶F3H[4]和二氢黄酮醇还原酶DFR[5]等。其中，查尔酮合成酶CHS是植物聚酮合成酶家族的一员，能够催化3个分子的丙二酰辅酶A和1个分子的对香豆酰辅酶A形成查尔酮，在植物黄酮类化合物合成途径中发挥重要作用[6]。研究表明，查尔酮合成酶基因CHS在猕猴桃果实转色开始期表达急剧升高，而注射CHSRNAi干扰载体后果肉颜色变浅，花青素含量降低[2，7]。可见，CHS在猕猴桃果实着色过程中可能发挥重要的功能。
自然界长期进化过程中，生物体或基因会逐渐形成适应环境所需的密码子使用法则，在这一过程中编码同一氨基酸的多个同义密码子并非均等随机使用，即存在密码子使用偏好性（Codon usage bias，CUB）。密码子使用偏好性能够反映基因的起源和进化规律，且与基因的转录水平[8]、内含子含量[9]、蛋白质翻译效率[10]和DNA甲基化[11]等密切相关。开展密码子使用偏好性研究对基因功能和分子进化机理的阐明具有重要参考价值。鉴于此，本试验以GenBank中登录的猕猴桃查尔酮合成酶基因CHS为研究对象，采用CodonW、SPSS190、MEGA52软件和EMBOSS在线程序分析其密码子使用偏好性，同时基于密码子使用偏好性和核酸序列相似性构建CHS系统发育树，以期为猕猴桃CHS后续研究提供科学依据。
1材料与方法
11基因序列获得
试验所研究猕猴桃CHS从GenBank中搜索并下载，全长1 475 bp，5′UTR、3′UTR和ORF分别124 bp、184 bp和1 170 bp，编码389个氨基酸，GenBank登录号为KF1573941。以CHS为检索名从NCBI中下载其他单子叶植物和双子叶植物查尔酮合成酶基因序列。
12方法
121密码子使用偏好性参数分析
采用CodonW分析CHS有效密码子数（Effective number of codons，ENc）、密码子总G和C含量（GC）、密码子第1、2和3位上G和C含量（GC1s、GC2s和GC3s）、密码子适应指数（Codon adaptation index，CAI）和同义密码子相对使用度（Relative synonymous codon usage，RSCU）等密码子偏好性参数。采用EMBOSS在线程序中的CUSP模块分析CHS密码子使用频率。模式生物拟南芥Arabidopsis thaliana、烟草Nicotiana sylvestris、大肠杆菌Escherichia coli以及酵母菌Saccharomyces cerevisiae的基因组密码子使用偏好性来源于Codon Usage Database 在线数据库。
122聚类分析
以121分析获得的物种间CHS同义密码子相对使用度RSCU为基础，基于SPSS190采用离差和平方法（Word法）构建单双子叶植物间CHS系统发育树，同时采用MEGA52组间邻接法进行CHS序列相似性聚类分析。
2结果与分析
21密码子使用偏好性参数分析
有效密码子数ENc是密码子偏好性的重要衡量标准。试验中，猕猴桃CHS ENc值为5721，远远高于35[12]而接近61，表明该基因密码子使用偏好性较低，氨基酸编码过程中对碱基的选择较为随机。而在密码子组成成分上，猕猴桃CHS的GC值为0533，说明其偏好使用含有较多G和C的密码子。此外，GC1s、GC2s和GC3s分别为0572、0403和0604，可见猕猴桃CHS密码子第1位和第3位倾向使用G或C，尤其对G或C结尾密码子有很强的选择偏好性。而猕猴桃CHS密码子适应指数CAI仅为0224，远小于10，也说明其密码子使用偏好性较弱。
22同义密码子相对使用度分析
同义密码子相对使用度RSCU表示特定密码子在对应同义密码子中的相对概率。猕猴桃CHS RSCU分析结果如表1所示，偏好性较强（RSCU>10）的密码子个数达到30个。在这些偏好性较强的密码子中，RSCU值均介于10～20，并不存在偏好性极强的密码子，由此也进一步反映猕猴桃CHS密码子使用偏好性较弱。此外，在编码半胱氨酸（Cys）时，猕猴桃CHS对TGC和TGT均等选择，而AGA则完全不参与精氨酸（Arg）的编码过程。23聚类分析
基于密码子偏好性进行单双子叶植物CHS聚类分析，结果如图1A所示。13个物种按照单双子叶的不同分为2个大类，表明单双子叶植物间CHS具有严格的密码子偏好模式差异。此外，双子叶植物类群又可以分为2个亚类群，类群Ⅰ中包含5种木本植物，类群Ⅱ包含2种草本植物，可见，双子叶植物中草本植物和木本植物之间CHS密码子使用规则也存在明显界限，而猕猴桃与甜橙Citrus sinensis的密码子使用特性较为相似。基于序列相似性的CHS聚类结果中，除拟南芥Arabidopsis thaliana外，其余双子叶植物均能准确归类，而单子叶植物CHS间的相似性较弱（图1B）。可见，双子叶植物间CHS的进化方向和动力可能较为一致，而在长期进化过程中部分单子叶植物开始朝不同的方向进行特征或功能分化。综合分析，单双子叶植物间CHS具有严格的密码子使用法则，可为物种间真实系统分类提供重要参考标准。
24异源转化受体系统筛选
密码子偏好性是物种间基因异源转化效率的重要影响因子，通常认为，二者密码子使用频率比值介于05～20则使用偏好性差异较小，否则差异明显。猕猴桃CHS与模式生物拟南芥（CHS/At）、烟草（CHS/Ns）、大肠杆菌（CHS/Ec）和酵母菌（CHS/Sc）基因组之间的密码子使用频率比值如图2所示，与4种模式基因组密码子使用偏好性差异较大的密码子个数分别为17、20、17和20个，表明在模式植物遗传转化过程中拟南芥可作为猕猴桃CHS较为理想的载体，而大肠杆菌可能更有利于该基因异源表达。然而，猕猴桃CHS与这些模式基因组仍存在部分偏好性差异极显著的密码子，在后续研究中仍需对这些密码子进行优化和定向改造。
3讨论
密码子使用偏好性是物种或基因长期进化和自然环境适应性的结果，对揭示基因的起源和进化具有重要参考价值。然而，猕猴桃中相关方面的研究尚未见报道。本试验中，猕猴桃CHS有效密码子数达到5721，密码子适应指数为0224，结果均表明猕猴桃CHS密码子使用偏好性较弱，在编码氨基酸过程中对密码子的选择具有较高的随机性；而在同义密码子相对使用度分析中，所有密码子的RSCU值均较小，预示该基因可能并不存在极强选择偏好的密码子或最优密码子群。密码子成分分析发现，猕猴桃CHS仍倾向使用G和C，且以G或C结尾的密码子，并不符合Murray[13]等公布的大部分双子叶植物密码子特性，可见猕猴桃CHS存在特殊的密码子使用和进化方式。此外，根据Novembre[14]和赖瑞联[15]等的观点推测，猕猴桃CHS进化过程中可能以由AT向GC突变为主。
亲缘关系较近的密码子往往表现相似的密码子使用法则，密码子使用偏好性分析可为物种间遗传关系和系统分类提供理论依据[15-16]。试验中，基于核酸序列相似性的CHS聚类結果发现，无论单子叶植物类群或双子叶植物类群CHS均难以准确归类，而密码子偏好性聚类结果更符合CHS真实系统分类学。结果预示着物种间CHS在长期进化过程中可能因为不同的进化方向和动力而逐渐产生序列特征差异，而单双子叶植物间仍存在的严格的密码子使用法则，与龙眼ARF3[17]、野生蕉FAD7[18]和橄榄CHI[19]研究结果较一致。由此也说明，采用序列相似性聚类和密码子使用偏好性聚类相结合更有利于提高分类结果的准确性。
采用遗传背景更为明确、技术手段更为成熟的模式生物系统进行基因异源转化是目前基因功能鉴定的常用方法之一。而目标基因和模式基因组之间的密码子使用偏好性差异会影响转化和表达效率，试验中往往有必要对其密码子进行优化和定向改造[20-21]。本研究表明，拟南芥和大肠杆菌可能作为猕猴桃CHS遗传转化和异源表达较为理想的受体，然而在具体试验中仍需对部分密码子进行优化。研究结果为猕猴桃CHS后续功能验证提供了重要理论依据。
参考文献：[1]
MONTEFIORI M， MCGHIE T K， COSTA G， et al. Pigments in the fruit of redfleshed kiwifruit （Actinidia chinensis and Actinidia deliciosa）[J]. Journal of Agricultural & Food Chemistry， 2005， 53（24）： 9526-9530.
[2]張亮， 满玉萍， 姜正旺， 等. 猕猴桃花青素合成途径基因AcCHS和AcLDOX的克隆与表达分析[J]. 园艺学报， 2012， 39（11）： 2124-2132.
[3]黄春辉， 葛翠莲， 张晓慧， 等. ‘红阳’猕猴桃突变体果实花青苷合成相关结构基因的表达分析[J]. 果树学报， 2014， 31（2）： 169-174.
[4]杨红丽， 王彦昌， 姜正旺， 等. ‘红阳’猕猴桃cDNA文库构建及F3H基因的表达初探[J]. 遗传， 2009， 31（12）： 1265-1272.
[5]杨俊， 姜正旺， 王彦昌. 红肉猕猴桃DFR基因的克隆及表达分析[J]. 植物科学学报， 2010， 28（6）： 673-681.
[6]王燕， 许锋， 程水源. 植物查尔酮合成酶分子生物学研究进展[J]. 河南农业科学， 2007， 36（8）： 5-9.
[7]SUN J， TANG W， LIU Y. Transient expression of CHSRNAi effectively influences the accumulation of anthocyanin in fruit of kiwifruit （Actinidia chinensis）[J]. Chinese Journal of Applied & Environmental Biology， 2014， 20（5）： 929-933.
[8]PEK H B， KLEMENT M， ANG K S， et al. Exploring codon context bias for synthetic gene design of a thermostable invertase in Escherichia coli[J]. Enzyme & Microbial Technology， 2015， 75-76： 57-63.
[9]RENSING S A， FRITZOWSKY D， LANG D， et al. Protein encoding genes in an ancient plant： analysis of codon usage， retained genes and splice sites in a moss， Physcomitrella patens[J]. BMC Genomics， 2005， 6（1）： 1-13.
[10]TULLER T， WALDMAN Y Y， KUPIEC M， et al. Translation efficiency is determined by both codon bias and folding energy[J]. Proceedings of the National Academy of Sciences of the United States of America， 2010， 107（8）： 3645-3650.
[11]ZHEN Q， CAI Z， XIA G， et al. Synonymous codon usage bias is correlative to intron number and shows disequilibrium among exons in plants[J]. BMC Genomics， 2013， 14（1）： 1-11.
[12]SHARP P M， LI W H. An evolutionary perspective on synonymous codon usage in unicellular organisms[J]. Journal of Molecular Evolution， 1986， 24（1）： 28-38.
[13]MURRAY E E， LOTZER J， EBERLE M. Codon usage in plant genes[J]. Nucleic Acids Research， 1989， 17（2）： 477-498.
[14]NOVEMBRE J A. Accounting for background nucleotide composition when measuring codon usage bias[J]. Molecular Biology and Evolution， 2002， 19（8）： 1390-1394.
[15]赖瑞联， 林玉玲， 钟春水， 等. 龙眼生长素受体基因TIR1密码子偏好性分析[J]. 园艺学报， 2016， 43（4）： 771-780.
[16]SHARP P M， COWE E， HIGGINS D G， et al. Codon usage patterns in Escherichia coli， Bacillus subtilis， Saccharomyces cerevisiae， Schizosaccharomyces pombe， Drosophila melanogaster and Homo sapiens； a review of the considerable withinspecies diversity[J]. Nucleic Acids Research， 1988， 16（17）： 8207-8211.
[17]赖瑞联， 林玉玲， 林梦桦， 等. 龙眼生长素响应因子ARF3密码子偏好性分析[EB/OL]. 北京： 中国科技论文在线[2016-04-30].
[18]赖志宸， 林玉玲. 野生蕉ω3脂肪酸去饱和酶基因FAD7密码子偏性分析[J]. 福建农业学报， 2017， 32（5）： 503-507.
[19]赖瑞联， 冯新， 陈瑾， 等. 橄榄查尔酮异构酶基因CHI的密码子偏好模式[J]. 应用与环境生物学报， 2017，23 （ 5 ） ： 0945-0951
[20]DUAN H， WANG H， MA B， et al. Codon optimization and expression of irisin in Pichia pastoris GS115[J]. International Journal of Biological Macromolecules， 2015， 79（1）： 21-26.
[21]Li C， Pan L L， Wang Y， et al. Codon bias of the gene for chloroplast glycerol3phosphate acyltransferase in Camellia sinensis （L.） O. Kuntze[J]. Biochemical Systematics & Ecology， 2014， 55（2）： 212-218.
（责任编辑：张梅）