普通油茶葉綠體基因組密碼子偏好性分析

王鵬良 楊利平 吳紅英 農有良 吳雙成 肖玉菲 覃子海 王華宇 劉海龍

摘 要： 為了利用葉綠體基因工程技術改良普通油茶的重要經濟性狀，該研究以普通油茶葉綠體全基因組序列為材料，從中篩選出51條長度大于300 bp且以ATG起始的非重復CDS （Coding DNA Sequence）為對象，利用CodonW軟件分析其密碼子偏好性。結果表明：密碼子第三位GC含量為27.55%，ENC范圍在35.23～56.67之間，平均值為46.09；RSCU值大于1.00的密碼子數目為30個，其中29個第三位堿基以U或A結尾；中性繪圖表明GC12與GC3的相關系數為0.143，相關性不顯著，回歸系數為0.0573；頻數分布顯示，55%基因的ENC比值集中分布在0～0.1，25%基因的ENC比值分布在0.1～0.2之間；對應分析結果表明，第一向量軸占10.12%的差異，第二向量軸占9.36%的差異，其余兩軸分別占7.97%和7.46%，前4軸累計差異為34.91%。中性繪圖、ENC-plot和對應性分析均表明普通油茶葉綠體基因密碼子偏好受突變作用，更多受選擇的影響。最終取高表達優越密碼子和高頻密碼子共有的CUU、AUU、GUU、GUA、UAA、CAA、AAA、GAC、GAA、CCU、ACU、GCU、GCA、UGU、CGU、AGU、UUG、GGU等18個密碼子作為最優密碼子。該研究結果為利用葉綠體基因工程技術改良普通油茶重要經濟性狀奠定了基礎。

關鍵詞： 普通油茶， 葉綠體， 密碼子偏好， ENC， RSCU

中圖分類號： Q943.2

文獻標識碼： A

文章編號： 1000-3142（2018）02-0135-10

Condon preference of chloroplast genome in Camellia oleifera

WANG Pengliang1，3*， YANG Liping2， WU Hongying2， NONG Youliang1， WU Shuangcheng1， XIAO Yufei3， QIN Zihai3， WANG Huayu1，2， LIU Hailong3

（ 1. Guangxi Key Laboratory of Beibu Gulf Marine Biodiversity Conservation， Qinzhou University， Qinzhou 535011， Guangxi， China； 2. Qinzhou Key Laboratory of Plant Biotechnology， Qinzhou Forestry Science Insititute， Qinzhou 535099， Guangxi， China； 3. Guangxi Key Laboratory of Special Non-Wood Forest Cultivation and Utilization， Guangxi Zhuang Autonomous Region Forestry Science Institute， Nanning 530002， China ）

Abstract： With the rapid development of the industry and increase of planting area of Camellia oleifera， the harvest was still affected by the factors such as the tolerance of disease， insect and drought as well as mating compatibility. Genetic engineering of chloroplast was an effective way to improve the important traits of plant. In order to improve the important traits of C. oleifera， it was necessary to carry out the analysis on codon preference of its chloroplast genome. The whole genome sequence of chloroplast in C. oleifera taken as raw material， 51 CDS who were longer than 300 bps， started with ATG and not repeated were screened for future analysis. The codon preference was conducted by the CodonW software. The results indicated that the GC content of the third base of genetic code was 27.55% and effective number of codons ranged from 35.23 to 56.67 with an average of 46.09 and then 29 codons of 30 whose RSCU were more than 1 ended with A or U. Neutral plot analysis showed the correlation between GC12 and GC3 was not significant and the coefficient of correlation and regression was 0.143 and 0.057 3， respectively. Frequency analysis showed 55% and 25% of the ENC ratio of genes focused on the region of 0-0.1 and 0.1-0.2. Correspondence analysis indicated the first axis accounted for 10.12 variation and the rested three axises accounted for 9.36%， 7.97% and 7.46% and the first four axises accounted for 34.91% in total. All the methods used in this study showed that codon preference was affected by not only mutation and but also selection. Finally the 18 optimal codons including CUU， AUU， GUU， GUA， UAA， CAA， AAA， GAC， GAA， CCU， ACU， GCU， GCA， UGU， CGU， AGU， UUG， GGU were determined. This work provides plenty of information for improving important traits of C. oleifera utilizing the genetic engineering of chloroplast.

Key words： Camellia oleifera， chloroplast， codon preference， effective number of codon， relative synonymous codon usage

普通油茶（Camellia oleifera）是山茶科（Theacae）山茶屬（Camellia）一類植物，是我國南方最重要的木本食用植油料物種之一（張恩慧等， 2016）。目前，油茶最具價值的部分為油茶果。茶籽油富含多種對人體健康有益的成分，對于提高人體的健康水平有很大影響，尤其是在人體免疫力和預防高血壓等方面，因此是一種非常優質的食用油，被譽為“東方的橄欖油”；同時茶籽油還含有多種有益于皮膚的成分，易于吸收，無刺激性，所以常被用于調制藥膏和高級化妝品的用油（莊瑞林， 2007）；其果殼由木質素、半纖維素和纖維素組成，可以用來制作糠醛、木糖醇、培養基、活性炭等；茶籽粕可以直接用于制作肥料及提取茶皂素、鞣質（單寧）、黃酮、生物堿等生物活性物質。

當前由于我國人口規模增長和居民收入水平提高，我國食用植物油消費總量快速增長，盡管近年來植物油生產穩步增長，然而消費增速高于生產增速，缺口仍然較大（張雯麗， 2016）。國家主動調整油料物種產業結構，由于油茶油質好、生長快、壽命長、適應性強、不占用耕地等優點，油茶產業得到重視，油茶產業獲得快速發展，新品種數量增加，推廣面積持續擴大。到目前為止，全國油茶種植面積約為426.67萬hm2（羅晨，2016）。

油茶產量持續增長，但油茶生產中還會遇到以下問題：首先，近年來大面積油茶病害的發生，如炭疽病，軟腐病和葉腫病（徐麗萍和檀根甲， 2015）；其次，大規模蟲害的出現，如蚧殼蟲，象鼻蟲等；第三，由于油茶的抗旱性弱導致低齡油茶植株死亡率高，引起油茶落花、落果，降低產量（左繼林等， 2012； 陳永忠等， 2014）；最后，油茶的交配親和性也嚴重影響了油茶的產量（常維霞和姚小華， 2016； 常維霞等， 2016）。

基因工程技術是改良植物重要性狀的有效手段 （王關林和方宏筠， 2014）。葉綠體基因組是理想的轉化載體，高表達優勢明顯，同時能將轉化基因的擴散控制在最低程度（Daniell & Chase， 2004）。葉綠體基因組中存在密碼子偏好，而且不同物種間密碼子偏好差異顯著（Morton， 1999； Liu & Xue， 2005）。為了更好地設計葉綠體高效表達載體，進一步利用葉綠體基因工程技術改良油茶重要性狀，有效降低油茶產量損失，分析普通油茶葉綠體基因組密碼子偏好十分必要。

1 材料與方法

1.1 序列材料

根據已報道的葉綠體基因組的GenBank登錄號NC_023084.1（Shi et al， 2013），登陸NCBI數據庫（https：//www.ncbi.nlm.nih.gov/），搜索并下載完整普通油茶葉綠體基因組序列及其CDS（Coding DNA Sequence），序列總長度為156 971 bp，共包含87條編碼蛋白的CDS。由于短序列無法正確計算有效密碼子數 （Wright， 1990），為了減少樣本誤差，本研究選擇以ATG為起始密碼子的非重復且長度大于300 bp的序列進行分析。

1.2 方法

1.2.1 密碼子偏好參數計算 以選取的CDS為研究對象，使用CodonW 1.4.2軟件和在線程序CUSP（http：//imed.med.ucm.es/EMBOSS/）分析密碼子使用偏好參數，即有效密碼子數（Effective Number of Codon， ENC）、相對同義密碼子使用度（Relative Synonymous Codon Usage， RSCU）、GC、GC1、GC2、GC3、GC3s。其中，GC、GC1、GC2、GC3分別代表了基因CDS序列的GC總含量和密碼子的對應位置上的GC含量，GC3s為同義密碼子第三位的GC含量。

1.2.2 中性繪圖分析 中性繪圖是一種可初步判斷影響密碼子使用偏好因素的方法。首先計算密碼子GC1和GC2的平均值GC12，再以計算出的平均值為縱坐標，以GC3為橫坐標作散點圖，圖中每一個點就代表一個基因的位置。由中性繪圖可分析密碼子的前兩位堿基和第三位堿基的關系，從而判斷密碼子使用偏好由突變或選擇造成。

1.2.3 ENC-plot繪圖 以GC3為橫坐標，ENC為縱坐標，構建二維散點圖，并在圖中構建出ENC值的標準曲線。每個基因點分布在標準曲線附近，這代表密碼子偏好性只是受到了突變作用的影響，如果各基因點集中在預期曲線下方，表明密碼子偏好更多受選擇影響。

標準曲線的公式如下（Wright， 1990）：

ENC=2+GC3+29GC23+（1-GC3）2

1.2.4 對應性分析 對應性分析是一種利用多元統計探究不同基因密碼子變異趨勢的方法。利用CodonW軟件的RSCU對應分析功能，將所有基因分布到59維向量空間，從而反映出密碼子使用偏好特點。通過分析基因位置探究基因向量及基因間的值的變異程度，推測導致密碼子偏好發生的原因。

1.2.5 最優密碼子確定 以ENC偏好性為標準排序，兩端各選出10%基因，分別建立高低偏性庫，取兩庫中ΔRSCU>0.08的密碼子作為高表達優越密碼子（楊惠娟等， 2012； 楊國鋒等， 2015）。將高表達優越密碼子與高頻密碼子相比較，確定出最優密碼子（劉慶坡和薛慶中， 2004； 續晨等， 2010； 羅洪等， 2015； 胡莎莎等， 2016）。

2 結果與分析

2.1 密碼子組成分析

根據密碼子使用偏好分析要求，本文篩選到51條CDS，再利用CodonW軟件和CUSP程序對篩選出的CDS進行分析。表1結果表明，密碼子第一位的GC平均含量為46.86%、第二位為39.28%、第三位為27.55%，第三位的GC含量明顯小于前兩位的GC含量，而且相差較大，說明普通油茶葉綠體基因密碼子偏好于A和U結尾。 ENC范圍在35.23～56.67之間，平均值為46.09。ENC值小于35說明密碼子偏性強，高于35則說明偏性比較弱的標準（Jiang et al， 2008），說明普通油茶葉綠體基因密碼子偏性較弱。

相關分析結果（表2）顯示，GC和GC1、GC2、GC3的相關性均達到極顯著水平，其相關系數分別為0.829、0.766、0.379。GC1和GC2的相關性也達到了極顯著水平，其相關系數為0.445，但GC3與GC1，GC2的相關性均未達到顯著相關水平，說明密碼子前兩位的GC含量比較相似，第三位GC含量和前兩位差異較大。ENC與GC3呈顯著相關，其相關系數為0.350，說明第三位堿基的組成會對密碼子的使用偏性有影響。密碼子數目N與GC3、ENC顯著相關，表明基因序列長度對GC3和密碼子的偏好性有一定的影響。

RSCU分析結果（表3）表明，RSCU值大于1.00的密碼子數目為30個。其中，16個第三位堿基為U；13個以A結尾；1個以G結尾。這說明普通油茶葉綠體基因組密碼子偏愛以A或U（T）結尾，不偏好以G或C結尾的密碼子。

2.2 中性繪圖分析

普通油茶葉綠體基因中性繪圖分析（圖1）顯示， GC12的取值范圍為0.319 6～0.550 4， GC3的范圍為0.191 2～0.371 0。GC12與GC3的相關系數為0.143，相關性未達到顯著水平，回歸系數為0.057 3，說明GC12和GC3 的相關性微弱。突變對密碼子第一位、第二位和對第三位堿基組成的影響不同。各基因都落在對角線的上方，并且沒有對角線上或沿著對角線分布。突變不是密碼子偏性形成的主要原因。因此，普通油茶葉綠體基因密碼子的使用更多地受選擇的影響。

2.3 ENC-plot繪圖

以GC3為X軸，ENC為Y軸作圖（圖2），圖2中各點為基因分布情況。代表基因的各點基本都落在了標準曲線的下方。為了能更直接地顯示出

差異，根據ENCexp-ENCobsENCexp公式，計算ENC比值，再統計ENC比值的分布（表4）。頻數分布顯示，55%的基因的ENC比值集中分布在0～0.1，25%的基因ENC比值分布在0.1～0.2的區間，結果表明ENCexp與ENCobs有差異。假如ENCexp與ENCobs接近時，表明基因的密碼子受突變影響；本文中ENCexp與ENCobs有差異較大，所以普通油茶葉綠體基因組密碼子偏好更多受選擇的影響（胡莎莎等， 2016）。

2.4 對應性分析

RSCU的對應分析結果表明，第一向量軸占10.12%的差異，第二向量軸占9.36%的差異，其余兩軸分別為7.97%和7.46%，前4軸累計差異為34.91%，說明第一軸對普通油茶葉綠體基因組密碼子偏好性造成的影響最大。第一軸與GC3的相關性為0.215和ENC的相關性為0.104，均相關不顯著。這說明GC3是對普通油茶葉綠體基因組密碼子偏好性造成影響的不大，密碼子的偏好性受到其他多方面因素的影響。以Axis 1為橫坐標，Axis 2為縱坐標，從各基因在第一、第二軸的分布（圖3）看出，與遺傳結構相關的基因分布相對比較集中，說明兩個類基因的密碼子使用是比較接近的；而其他類別基因可能是因為密碼子使用模式不同，分布比較分散。

2.5 最優密碼子分析

以密碼子的有效數量ENC為標準對所有基因的CDS序列進行排序，先從排序結果的兩端各選取5個基因分別建立高表達基因庫和低表達基因庫；然后使用CodonW計算兩個表達庫中的RSCU值，并計算高低庫ΔRSCU值（表5）；最后確定22個密碼子（表中帶星號的）為普通油茶葉綠體基因組高表達優越密碼子，10個以U結尾、6個以A結尾、5個以C結尾和1個以G結尾。

先將22個高表達優越密碼子和表3得出的30個高頻密碼子相比較，然后選取共有的密碼子作為最優密碼子，最后確定了18個密碼子為普通油茶葉綠體基因組最優密碼子，分別為CUU、AUU、GUU、GUA、UAA、CAA、AAA、GAC、GAA、CCU、ACU、GCU、GCA、UGU、CGU、AGU、UUG、GGU，其中有2個密碼子以G或C結尾，其余的16個密碼子均以A或U結尾。

3 討論

在20種氨基酸中，除了甲硫氨酸和色氨酸只有單個密碼子外，其他18種氨基酸均有多種密碼子，即遺傳密碼子存在簡并性。密碼子簡并性常常表現在第三位密碼子上（朱圣庚和徐長發， 2016），因此在普通油茶葉綠體基因組堿基組成中發現，GC1、GC2顯著相關，GC3與GC1、GC2相關不顯著，且GC3小于GC1和GC2，說明普通油茶葉綠體基因組的堿基組成偏向于A和T，第3位密碼子通常以A、T結尾。這與多數植物，如文心蘭（李冬梅等， 2012）、甜蕎（羅洪等， 2015）、苦蕎（胡莎莎等， 2016）、蝴蝶蘭（續晨等， 2010）、蒺藜苜蓿（楊國鋒等， 2015）等物種的特點相一致。

突變和選擇是密碼子偏好形成的主要原因。ENC-plot分析結果表明，部分基因落點位于標準曲線附近，基因的實際ENC值與理論ENC值基本一致，說明這些密碼子偏好受到GC含量（突變）的影響較大；而大多數基因的落點離標準曲線較遠，實際ENC值與理論ENC值的差異較大，受選擇的影響較大。因此，普通油茶葉綠體密碼子偏好受到突變的影響，更多的是受到了選擇的作用。

中性繪圖和對應性分析均支持普通油茶密碼偏好受突變和選擇共同作用。這與已報道的文心蘭（李冬梅等， 2012）、蝴蝶蘭（續晨等， 2010）、銀白楊（Zhou et al， 2008）、甜蕎（羅洪等， 2015）等有相似的結論；本文結果也印證了選擇對葉綠體密碼子使用起顯著影響（Morton， 1998）。

本研究采用由高表達優越密碼子和高頻密碼子中選出兩者共有的密碼子，可以避免由于單方法存在問題而導致結果出錯，最終確定了CUU、AUU、GUU、GUA、UAA、CAA、AAA、GAC、GAA、CCU、ACU、GCU、GCA、UGU、CGU、AGU、UUG、GGU等18個密碼子為普通油茶葉綠體基因的最優密碼子。在此基礎上可以通過對目的基因的密碼子進行優化，提高該基因的表達效率，從而為改良普通油茶重要性狀奠定堅實基礎。

參考文獻：

CHANG WX， YAO XH， 2016. Study onself-compatibility of Camellia oleifera clones [J]. For Res， 29（4）： 508-514. [常維霞，姚小華， 2016. 油茶無性系自交親和性分析 [J]. 林業科學研究，29（4）： 508-514.]

CHANG WX， YAO XH， LONG W， et al， 2016. Cross-compatibility of four kinds of Camellia species [J]. Bull Bot Res， 36（4）： 527-534. [常維霞，姚小華，龍偉，等，2016. 4個油茶物種雜交親和性分析 [J]. 植物研究，36（4）： 527-534.]

CHEN YZ， LUO J， CHEN LS， et al， 2014. Drought injury and countermeasure of Camellia oleifera in Hunan Province [J]. Nonwood For Res， 32（3）： 22-29. [陳永忠，羅健，陳隆升，等，2014. 湖南省油茶旱害及防治對策 [J]. 經濟林研究，32（3）： 22-29.]

DANIELL H， CHASE C， 2004. Molecular biology and biotechnology of plant organelles [M]. Dordrecht， Springer.

HU SS， LUO H， WU Q， et al， 2016. Analysis of codon bias of chloroplast genome of Fagopyrum tataricum [J]. Mol Plant Breed， 14（2）： 309-317. [胡莎莎，羅洪，吳琦，等，2016. 苦蕎葉綠體基因組密碼子偏愛性分析 [J]. 分子植物育種，14（2）： 309-317.]

JIANG Y， DENG F， WANG H， et al， 2008. An extensive analysis on the global codon usage pattern of baculoviruses [J]. Arch Virol， 153（12）： 2273-2282.

LI DM， L FB， ZHU GF， et al， 2012. Analysis on codon usage of chloroplast genome of Oncidium Gower Ramsey [J]. Guangdong Agric Sci， 10： 61-65. [李冬梅，呂復兵，朱根發，等，2012. 文心蘭葉綠體基因組密碼子使用的相關分析 [J]. 廣東農業科學，10： 61-65.]

LIU Q， XUE Q， 2005. Comparative studies on codon usage pattern of chloroplasts and their host nuclear genes in four plant species [J]. J Gene， 84（1）： 55-62.

LIU QP， XUE QZ， 2004. Codon usage in the chloroplast genome of rice （Oryza sativa L. ssp. japonica） [J]. Acta Agronom Sin， 30（12）： 1220-1224. [劉慶坡，薛慶中，2004. 粳稻葉綠體基因組的密碼子用法 [J]. 作物學報，30（12）： 1220-1224.]

LUO C， 2016. The area of Camellia achieved 64 million acres [EB/OL]. WWW. cnfood. cn/n/2016/1219/98628. html. [羅晨， 2016. 全國油茶面積已達6400多萬畝 [EB/OL]. WWW. cnfood. cn/n/2016/1219/98628. html.]

LUO H， HU SS， WU Q， et al， 2015. Analysis of buckwheat chloroplast gene codon bias [J]. Genomics Appl Biol， 34（11）： 2457-2464. [羅洪，胡莎莎，吳琦，等，2015. 甜蕎葉綠體基因密碼子偏愛性分析 [J]. 基因組學與應用生物學，34（11）： 2457-2464.]

MORTON BR， 1998. Selection on the codon bias of chloroplast and cyanelle genes in different plant and algal lineages [J]. J Mol Evol， 46（4）： 449-459.

MORTON BR， 1999. Strand asymmetry and codon usage bias in the chloroplast genome of Euglena gracilis [J]. Proc Nat Acad Sci USA， 96（9）： 5123-5128.

SHI C， LIU Y， HUANG H， et al， 2013. Contradiction between plastid gene transcription and function due to complex posttranscriptional splicing： An exemplary study of ycf15 function and evolution in angiosperms [J]. PLoS ONE， 8（3）： e59620.

WANG GL， FANG HY， 2014. Plant genetic engineering [M]. 2nd ed. Beijing： Science Press. [王關林，方宏筠， 2014. 植物基因工程 [M]. 第2版. 北京： 科學出版社.]

WRIGHT F， 1990. The ‘effective number of codons used in a gene [J]. Gene， 87（1）： 23-29.

XU LP， TAN GJ， 2015. Epidemic of Camellia oleifera diseases and their relations to the ecological conditions [J]. J Anhui Agric Univ， 42（2）： 272-275. [徐麗萍，檀根甲， 2015. 油茶主要病害流行與生態條件的關系和生態調控技術 [J]. 安徽農業大學學報，42（2）： 272-275.]

XU C， BEN AL， CAI XN， 2010. Analysis of synonymous codon usage in chloroplast genome of Phalaenopsis aphrodite subsp. formosana [J]. Mol Plant Breed， 8（5）： 945-950. [續晨，賁愛玲，蔡曉寧，2010. 蝴蝶蘭葉綠體基因組密碼子使用的相關分析 [J]. 分子植物育種，8（5）： 945-950.

YANG GF， SU KL， ZHAO YR， et al， 2015. Analysis of codon usage in the chloroplast genome of Medicago truncatula [J]. Acta Pratacul Sin， 171-179. [楊國鋒，蘇昆龍，趙怡然，等，2015. 蒺藜苜蓿葉綠體密碼子偏好性分析 [J]. 草業學報， 24（12）： 171-179.]

YANG HJ， LIU GS， ZHANG ST， et al， 2012. Codon usage bias studies and cluster analysis on chloroplastic genes in tobacco [J]. Acta Tabac Sin， 18（2）： 37-43. [楊惠娟，劉國順，張松濤，等， 2012. 煙草葉綠體密碼子的偏好性及聚類分析 [J]. 中國煙草學報，18（2）： 37-43.]

ZHANG EH， WANG XY， QIN ZH， et al， 2016. Genetic diversity analysis of Camellia oleifera in Guangxi using SSR markers [J]. Guihaia， 36（7）： 806-811. [張恩慧，王曉云，覃子海，等， 2016. 廣西普通油茶種質資源遺傳多樣性的SSR分析 [J]. 廣西植物，36（7）： 806-811.]

ZHANG WL， 2016. Supply and demand of Chinas oilseeds and edible vegetable oil and their prospects during the 13th five-year plan period [J]. Agric Outlook， （2）： 4-8. [張雯麗，2016. 中國油料及食用植物油供需形勢及“十三五”展望 [J]. 農業展望，（2）： 4-8.]

ZHU SG， XU CF， 2016. Biochemistry [M]. 4th ed. Beijing： Higher Education Press. [朱圣庚，徐長發，2016. 生物化學 [M]. 第4版. 北京： 高等教育出版社.]

ZHUANG RL， 2007. Chinese oil tea [M]. 2nd ed. Beijing： Chinese Forestry Publishing House. [莊瑞林，2007. 中國油茶 [M]. 第2版. 北京： 中國林業出版社.]

ZHOU M， LONG W， LI X， 2008. Analysis of synonymous codon usage in chloroplast genome of Populus alba [J]. J For Res， 19（4）： 293-297.

ZUO JL， WANG YJ， GONG C， et al， 2012. Effects of different management measures on the fruit shape growth of high yield Camellia oleifera under during the summer drought period [J]. J Centr Southern Univ For & Technol， 32（4）： 15-20. [左繼林，王玉娟，龔春，等，2012. 高產油茶夏旱期不同經營措施對其果形生長的影響 [J]. 中南林業科技大學學報，32（4）： 15-20.]