黃笑宇，許在恩，郭小勤

（浙江農(nóng)林大學(xué) 亞熱帶森林培育國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地，浙江 臨安 311300）

基于全基因組的毛竹同義密碼子使用偏好性分析

黃笑宇，許在恩，郭小勤

（浙江農(nóng)林大學(xué) 亞熱帶森林培育國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地，浙江 臨安 311300）

密碼子使用偏好性是物種在遺傳信息傳遞過程中的一個(gè)重要特點(diǎn)，分析物種的密碼子使用偏好性對(duì)于了解該物種遺傳信息的傳遞規(guī)律具有重要意義。應(yīng)用CodonW軟件對(duì)毛竹Phyllostachys edulis基因組中的26 103個(gè)蛋白質(zhì)編碼基因序列進(jìn)行了分析，計(jì)算了位于密碼子3個(gè)位置的G＋C含量、有效密碼子數(shù)、同義密碼子的使用頻率等，確定了毛竹的最優(yōu)密碼子。結(jié)果顯示：毛竹密碼子第1位和第3位的G＋C含量明顯高于第2位，表現(xiàn)出對(duì)以G或C堿基開頭和結(jié)尾的密碼子發(fā)生強(qiáng)烈偏向使用，且確定的26種最優(yōu)密碼子均以G/C堿基結(jié)尾。與模式動(dòng)植物酵母Saccharomys cerevisiae，大腸埃希菌Escherichia coli，果蠅Drosophila melanogaster，擬南芥Arabidopsis thaliana，煙草Nicotiana tabacum，水稻Oryza sativa，玉米Zea mays，小麥Triticum aestivum等8個(gè)代表性物種進(jìn)行比較，結(jié)果顯示：毛竹密碼子偏好性與模式動(dòng)植物存在不同程度的差異，與大腸埃希菌、酵母、果蠅、擬南芥、煙草和玉米的差異較大，差異分別為8個(gè)、11個(gè)、6個(gè)、7個(gè)、6個(gè)和6個(gè)，而與小麥的差異較?。?個(gè)），與水稻完全一致。密碼子偏好性差異大小在一定程度上反映物種間的進(jìn)化關(guān)系。圖1表4參29

植物學(xué)；毛竹；密碼子偏好性；最優(yōu)密碼子

三聯(lián)密碼子是整個(gè)生物王國(guó)的核心，作為最基本的編碼組分編碼特定的氨基酸。除了甲硫氨酸Met和色氨酸Trp外，同一個(gè)氨基酸會(huì)由2~6個(gè)同義密碼子編碼［1］。根據(jù)中心法則，盡管同義突變不會(huì)引起蛋白序列的變化，但同義密碼子使用偏好性在基因組內(nèi)和基因組間廣泛存在［2－4］。密碼子使用偏好性是物種在遺傳信息傳遞過程中的一個(gè)重要特點(diǎn)，分析物種的密碼子使用偏好性對(duì)于了解該物種遺傳信息的傳遞規(guī)律具有重要意義。密碼子使用偏好性的研究有助于更好地理解分子生物及進(jìn)化，信使核糖核酸（mRNA）翻譯，轉(zhuǎn)基因設(shè)計(jì)，新基因發(fā)現(xiàn)，以及其他生物應(yīng)用［3－6］。幾十年來，大量模式物種如擬南芥Arabidopsis thaliana，水稻Oryza sativa，果蠅Drosophila melanogaster，楊樹Populus trichocarpa等的測(cè)序產(chǎn)生了大量的開放閱讀框，這些全長(zhǎng)編碼序列作為密碼生物學(xué)的基礎(chǔ)，為研究密碼子使用模式提供了強(qiáng)有力的保障。隨著深度測(cè)序技術(shù)的快速發(fā)展，非模式植物也紛紛被測(cè)序，產(chǎn)生了大量的序列，有關(guān)密碼子使用模式的工作也逐漸拉開序幕［7］。毛竹Phyllostachys edulis為禾本科Gramineae多年生木本植物，地上部分可材用，地下部分發(fā)育的筍可食用。毛竹基因組序列測(cè)定［8］，產(chǎn)生了大量基因組及編碼區(qū)序列，為研究毛竹密碼子的使用情況提供了大量的信息。本研究通過分析毛竹全基因組編碼序列數(shù)據(jù)，了解毛竹基因密碼子用法特征，并與不同代表性物種進(jìn)行比較。這些分析有助于我們理解毛竹的密碼子模式，提升植物密碼子使用的研究，同時(shí)為毛竹基因選擇合適的表達(dá)系統(tǒng)，優(yōu)化密碼子提高基因表達(dá)量等提供重要理論基礎(chǔ)。

1 材料與方法

1.1 序列數(shù)據(jù)

從http://202.127.18.221/bamboo/index.php的毛竹基因組注釋數(shù)據(jù)中獲取了31 987條蛋白質(zhì)對(duì)應(yīng)的編碼基因序列（coding DNA sequence，CDS），從中挑選出以ATG為起始密碼子，以TAA，TAG或TGA為終止密碼子的，且CDS長(zhǎng)度大于300 bp的26 103個(gè)基因作為序列分析樣本［9］。這個(gè)數(shù)據(jù)庫中已經(jīng)去除所有假基因的信息。

1.2 序列處理

采用C語言編寫程序進(jìn)行序列篩選與處理。

1.3 同義密碼子使用偏好性分析

采用EMBOSS軟件包中的CHIPS和CUSP程序在線http://emboss.bioinformatics.nl及CodonW1.4.4（http://mobyle.pasteur.fr/cgi-bin/portal.py?#forms::codonw）對(duì)毛竹全基因編碼序列進(jìn)行分析，計(jì)算有效密碼子數(shù)（effective number of codons，Enc），CDS區(qū)的GC含量，密碼子中第3位堿基的GC含量（GC3s），同義密碼子相對(duì)使用頻率（relative synonymous codon usage，RSCU）及密碼子使用概率。

衡量同義密碼子使用偏好性參數(shù)的含義：①有效密碼子數(shù)（Enc）。該值被認(rèn)為是在評(píng)價(jià)基因整體密碼子偏好性用法中最具有參考價(jià)值的參數(shù)之一，目前被廣泛用于評(píng)價(jià)基因密碼子偏好性，其取值范圍為20（每個(gè)氨基酸只使用1個(gè)密碼子的極端情況）到61（各個(gè)密碼子均被平均使用）［10］。②同義密碼子相對(duì)使用頻率RSCU。該值的計(jì)算方法為某一密碼子所使用的頻率與其在無偏好使用時(shí)預(yù)期頻率之間的比值，若某一密碼子的RSCU值等于1，則表明該密碼子的使用沒有偏好性；RSCU值大于l，表明該密碼子的使用頻率相對(duì)較高，反之亦然。它去除了氨基酸組成對(duì)密碼子使用的影響，且直觀地反映了密碼子使用的偏好性［11］。③同義密碼子使用的絕對(duì)頻率（Fract）。該值表示各個(gè)密碼子在編碼該氨基酸的密碼子中所占的比例（各比例相加總和為1）［12］。

1.4 最優(yōu)密碼子的確定

采用STENICO等［13］的方法，把密碼子使用偏好性強(qiáng)和弱的2組基因之間相應(yīng)密碼子出現(xiàn)頻率之差達(dá)到統(tǒng)計(jì)學(xué)上顯著水平的密碼子定義為最優(yōu)密碼子。具體方法如下：通過計(jì)算樣本中每個(gè)基因的有效密碼子數(shù)，并按該值的大小對(duì)基因進(jìn)行排列，從這一排列的兩端各取基因樣本總數(shù)的5%，分別組成高、低表達(dá)樣本組。計(jì)算這2組基因的相對(duì)密碼子使用度，并進(jìn)行卡方檢驗(yàn)，確定最優(yōu)密碼子［14］。

1.5 毛竹與其他物種密碼子偏好性比較

運(yùn)用CUSP程序計(jì)算毛竹基因各密碼子的使用頻率，并與從Codon Usage Database（http://www.kazusa. or.jp/codon/）中獲得的果蠅，酵母Saccharomys cerevisiae，大腸埃希菌Escherichia coli，擬南芥，煙草Nicotiana tabacum，水稻，玉米Zea mays，小麥Triticum aestivum等的密碼子使用頻率進(jìn)行比較。密碼子使用頻率若為 0.5~2.0，表明這2個(gè)物種對(duì)該密碼子的偏好性較接近，若≥2.0或≤0.5，則表明偏好性差異較大［15］。

1.6 基于密碼子使用偏好性的聚類

利用SPSS 19.0對(duì)毛竹及其他9個(gè)物種進(jìn)行基于密碼子使用偏好性的聚類分析，方法參考文獻(xiàn)［15］。

2 結(jié)果與分析

2.1 有效密碼子數(shù)與密碼子的堿基組成

將毛竹基因組注釋數(shù)據(jù)中獲取到的31 987條蛋白質(zhì)對(duì)應(yīng)的編碼基因序列進(jìn)行篩選后，獲得26 103條有效序列，將這些序列作為一整體，在線計(jì)算了其有效密碼子數(shù)及密碼子第1位、第2位、第3位和3個(gè)位置平均的堿基GC百分率，結(jié)果見表1。毛竹基因整體的有效密碼子數(shù)為57.88，表明毛竹整體基因的密碼子存在一定程度偏好，但偏好性不強(qiáng)。從密碼子的GC含量來看，3個(gè)位置平均GC含量為0.52。其中，GC1含量為0.56，比GC2（0.44）高0.12，而GC3含量（0.57）比GC1略高0.01，表明選擇壓力使得毛竹密碼子的第1位傾向于選擇G/C，第2位傾向于選擇T/A，第3位可以有大幅擺動(dòng)。

表1 毛竹基因密碼子中3個(gè)位置的GC含量及有效密碼子數(shù)值Table 1 GC content of different positions and effective number of codons in Phyllostachys edulis

2.2 同義密碼子的使用頻率

經(jīng)軟件計(jì)算的同義密碼子使用次數(shù)及頻率結(jié)果見表2。在64個(gè)密碼子中，GAG是出現(xiàn)次數(shù)最高的密碼子，絕對(duì)頻率為39.21，是毛竹平均頻率的2.51倍；緊隨其后的是AAG，達(dá)到33.90；處于第3位的是GAU，為28.40；CGA的出現(xiàn)頻率最低，僅為5.34，是毛竹平均頻率的1/3；另有6個(gè)密碼子（UUA，CUA，GUA，ACG，UGU，CGU）的頻率小于10.00（表2）。

有34個(gè)密碼子的RSCU值大于1，這些密碼子為毛竹基因的偏好密碼子，其中約1/3的密碼子以A/U結(jié)尾，2/3的密碼子以G/C結(jié)尾。AGG（編碼Arg），CUC（編碼Leu）和GUG（編碼Val）的RSCU值處于前3位，分別為1.57，1.45和1.39。CUG（編碼Cys）以及AAG（編碼Lys）和GGC（編碼Gly）相對(duì)于其同義密碼子的使用頻率高，分別為1.36和1.31。這5個(gè)密碼子為本文的高頻率密碼子。

4個(gè)NUA密碼子的RSCU值最低，AUA為0.69，CUA為0.52，GUA為0.47，UUA為0.45，表明這幾個(gè)是毛竹基因避免使用的密碼子。4個(gè)NCG的RSCU值相對(duì)來說接近于平均水平甚至更低，CCG為0.99，GCG為0.91，ACG為0.76，UCG為0.75，表明毛竹體內(nèi)的甲基化水平可能較低或中等，這點(diǎn)從NCG∶NCC的比值（為0.82）也可看出。終止密碼子UGA在毛竹基因中的使用頻率較其余2個(gè)終止密碼子高，為1.02，其次是UAG，RSCU值為0.67，UAA的使用頻率最低，僅為0.55。

2.3 最優(yōu)密碼子的確定

不僅同義密碼子間存在偏好性，且密碼子本身的使用也存在偏好性。目前，關(guān)于毛竹基因表達(dá)的數(shù)據(jù)偏少，多數(shù)轉(zhuǎn)錄組測(cè)序的數(shù)據(jù)也基于幾個(gè)毛竹的特異組織。因此，本研究依據(jù)Enc值來衡量基因的表達(dá)量。表3中的結(jié)果是通過計(jì)算高表達(dá)/低表達(dá)基因之間同義密碼子相對(duì)使用頻率之差，經(jīng)卡方測(cè)驗(yàn)確定的毛竹中的最優(yōu)密碼子，用*號(hào)標(biāo)記，共26個(gè)。這些密碼子均以G/C結(jié)尾，表明在高表達(dá)基因中優(yōu)先使用這些密碼子。這些密碼子的使用頻率在高表達(dá)基因組與低表達(dá)基因組之間的差異達(dá)到極顯著水平。

2.4 毛竹與模式動(dòng)植物密碼子偏好性比較

將毛竹與3種模式生物大腸埃希菌、酵母和果蠅密碼子使用頻率比較，比值0.5~2.0表明2物種使用該密碼子的偏好性相似，比值小于0.5或大于2.0，表明該密碼子的使用偏好性差異較大。結(jié)果顯示：毛竹與大腸埃希菌、酵母和果蠅密碼子的比值中，分別有8，11，6個(gè)小于0.5或大于2.0，表明毛竹與這些模式生物之間的密碼子偏好性存在一定差異。

表2 毛竹基因同義密碼子的使用頻率Table 2 Frequency of synonymous codons in genes of Phyllostachys edulis

表4的結(jié)果顯示：毛竹與雙子葉植物的代表種擬南芥和煙草的密碼子偏好性差異性較大，比值大于2.0或小于0.5的分別有7個(gè)和6個(gè)，與同科植物相比，與C4植物玉米的密碼子偏好性差異也較大，有6個(gè)，而與C3植物水稻和小麥的偏好性一致。

2.5 不同物種間密碼子偏好性的聚類分析

根據(jù)各物種編碼序列密碼子的使用頻率，利用SPSS 19.0進(jìn)行聚類分析（圖1）。從圖1可以看出：雙子葉植物擬南芥和煙草密碼子使用偏好更相近，禾本科植物毛竹與水稻的最近，其次與小麥和玉米。利用密碼子使用頻率得出的聚類結(jié)果一定程度上反映了各物種間的進(jìn)化關(guān)系。

表3 毛竹中高/低表達(dá)樣本的密碼子用法Table 3 Codon usage of high/low expressed genes in Phyllostachys edulis

3 討論

在長(zhǎng)期的進(jìn)化過程中，不同物種對(duì)進(jìn)化環(huán)境和選擇壓力的適應(yīng)不同，因此，任何一個(gè)物種都會(huì)形成特定的密碼子用法以適應(yīng)其基因組環(huán)境，最終使其宿主適應(yīng)外界進(jìn)化環(huán)境。由此，不同物種就形成了各自特定的密碼子偏好性。若要通過基因工程技術(shù)改造某一物種或?qū)⒛骋换蛴糜隗w外表達(dá)，應(yīng)先按照宿主的密碼子使用偏好性對(duì)所導(dǎo)入的基因進(jìn)行優(yōu)化和改造。本研究在毛竹全基因組測(cè)序的基礎(chǔ)上，對(duì)編碼蛋白基因的密碼子偏好性進(jìn)行了分析，結(jié)果表明與很多物種包括人、細(xì)菌、酵母、果蠅、玉米、草菇Volvariella volvacea等相似［15－16］，毛竹基因密碼子偏好使用G/C結(jié)尾的密碼子，這種密碼子使用偏好性有利于保證翻譯的準(zhǔn)確性［17－18］，但與同為植物界的雙子葉植物相比，密碼子使用偏好性差異較大，雙子葉植物偏好使用以A或T結(jié)尾的密碼子［14,19－20］。

圖1 基于不同物種密碼子使用頻率的聚類分析Figure 1 Cluster analysis dendrogram of frenquency of codon usage of different species

表4 毛竹與模式植物的密碼子偏好性比較Table 4 Comparison of codon preference between bamboo and other model plants

表4 （續(xù)）Table 4 Continued

從GC含量上看，很多植物的密碼子的GC1含量均比GC2高，兩者含量的差異達(dá)到0.096（Medicago truncatula）~0.155（Micromonas pusilla RCC299）。裸子植物、單子葉植物、綠藻等物種GC3的含量一般來說略高于GC1［7］。本研究中的毛竹GC1含量比GC2高0.12，而GC3含量比GC1高0.01，表明選擇壓力對(duì)毛竹密碼子不同位置的堿基組成影響不同。不同物種中GC3會(huì)隨著進(jìn)化不同而發(fā)生變化［7,21－22］。一般來說，原始的單細(xì)胞或多細(xì)胞綠色植物GC3含量會(huì)比較高，為0.690~0.854，苔蘚植物為0.481~0.578，而被子植物GC3s的含量變異差異比較大，單子葉植物的變異范圍為0.581~0.609，優(yōu)等雙子葉植物的變異范圍為0.335~0.482。本研究中毛竹GC3s為0.52，超出了單子葉植物的變異范圍。這樣的特例在其他物種中也有發(fā)生，如萊茵衣藻Chlamydomonas reinhardtii，團(tuán)藻Volvox carteri和 細(xì)小微胞藻Micromonas pusilla［7］。

本研究用CodonW軟件分析了毛竹同義密碼子的RSCU值，發(fā)現(xiàn)AGG，CUC，GUG，AAG和UGC 5個(gè)密碼子為本文的高頻率密碼子。而4個(gè)NUA密碼子RSCU值較低，AUA為0.69，CUA為0.52，GUA為0.47，UUA為0.45，表明毛竹基因避免使用UA密碼子，同一現(xiàn)象在其他物種中也發(fā)現(xiàn)，可能因?yàn)榈秃康腢A抑制了mRNA的降解，提高蛋白產(chǎn)物或產(chǎn)量［23］。毛竹中終止密碼子的使用以UGA的使用頻率最高，與大多數(shù)植物相吻合［24］。

NCG∶NCC的比值已廣泛用于評(píng)估CpG抑制，反映了編碼區(qū)甲基化水平，尤其在真雙子葉植物。甲基化水平低的物種往往其 NCG∶NCC的比值相對(duì)較高，如擬南芥（0.921），深山南芥 Arabis lyrata（0.93）；而高甲基化水平的物種，該比值相對(duì)較低，如葡萄Vitis vinifera（0.414），楊樹（0.463）；甲基化程度中等的物種；該比值中等，如蘋果Malus×domestica（0.639），番茄Solanum lycopersicon（0.634）。毛竹中該比值為0.819 7，表明毛竹為低甲基化水平的物種。由此可以判斷：甲基化水平對(duì)毛竹的生長(zhǎng)發(fā)育過程影響有限［7］。

本研究使用同義密碼子相對(duì)使用頻率（RSCU）方法鑒定出26個(gè)最優(yōu)密碼子，全部以G/C結(jié)尾，毛竹編碼蛋白序列的GC含量平均為52.4%，因此，本研究結(jié)果符合一般規(guī)律，即富含GC堿基的基因組中最優(yōu)密碼子也富含GC［2,25］。通過比較某一特定基因與外源表達(dá)系統(tǒng)之間的密碼子使用偏好性差異，從而分析是否會(huì)引起甲基化，導(dǎo)致基因表達(dá)量下降或基因沉默［26］，從而改造密碼子以提高外源基因在宿主中的表達(dá)［27－29］。本研究將毛竹基因組密碼子的偏好性與模式動(dòng)植物大腸埃希菌、酵母、果蠅、擬南芥、煙草、玉米、水稻和小麥待密碼子偏好性相比，結(jié)果表明：毛竹與不同物種的差異程度不同，其中與大腸埃希菌和酵母的差異最大，而與同科C3植物水稻和小麥的偏好性一致。因此，要將毛竹基因進(jìn)行體外表達(dá)時(shí)，需要通過密碼子的改造，來提高表達(dá)效率。若要將毛竹基因用于水稻和小麥中表達(dá)時(shí)，可以不用經(jīng)密碼子優(yōu)化直接進(jìn)行外源基因表達(dá)。本研究的聚類結(jié)果表明：密碼子偏好性差異大小在一定程度上反映物種間的進(jìn)化關(guān)系，與傳統(tǒng)分類有一定的吻合性，但不完全吻合，這與其他物種的基于密碼子偏好性聚類的結(jié)果類似［4,15］，很可能是因?yàn)閰?shù)選擇單一造成的。該研究結(jié)果可為毛竹基因外源表達(dá)選擇合適的受體提供理論基礎(chǔ)，同時(shí)，還為將毛竹基因轉(zhuǎn)入模式生物中進(jìn)行功能驗(yàn)證提供基礎(chǔ)資料。

4 結(jié)論

本研究對(duì)毛竹基因組中的26 103個(gè)蛋白質(zhì)編碼基因序列進(jìn)行了分析，根據(jù)同義密碼子相對(duì)使用頻率（RSCU值）確定了毛竹中的最優(yōu)密碼子26個(gè)，且均以G/C結(jié)尾。同時(shí)與模式動(dòng)植物9個(gè)代表性物種進(jìn)行了比較，毛竹密碼子偏好性與水稻完全一致。

［1］ 宋輝，王鵬飛，馬登超，等.蒺藜苜蓿WRKY轉(zhuǎn)錄因子密碼子使用偏好性分析［J］.農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，2015, 23（2）：203-212.

SONG Hui,WANG Pengfei,MA Dengchao,et al.Analysis of codon usage bias of WRKY transcription factors in Medicago truncatula［J］.J Agric Biotechnol,2015,23（2）:203-212.

［2］ HERSHBERG R,PETROV D A.Selection on codon bias［J］.Annu Rev Genet,2008,42（42）:287-299.

［3］ 劉漢梅，何瑞，趙耀，等.玉米密碼子用法分析［J］.核農(nóng)學(xué)報(bào)，2008，22（2）：141-147.

LIU Hanmei,HE Rui,ZHAO Yao,et al.Analysis of codon usage in maize［J］.Acta Agric Nucl Sin,2008,22（2）:141 -147.

［4］ 晁岳恩，吳政卿，楊會(huì)民，等.11種植物psbA基因的密碼子偏好性及聚類分析［J］.核農(nóng)學(xué)報(bào)，2011，25（5）：927-932.

ZHAO Yueen,WU Zhengqing,YANG Huimin,et al.Cluster analysis and codon usage bias studies on psbA genes from 11 plant species［J］.Acta Agric Nucl Sin,2011,25（5）:927-932.

［5］ DURET L,MOUCHIROUD D.Expression pattern and,surprisingly,gene length shape codon usage in Caenorhabditis, Drosophila,and Arabidopsis［J］.Proc Nat Acad Sci U S A,1999,96（8）:4482-4487.

［6］ SAKAI H,WASHIO T,SAITO R,et al.Correlation between sequence conservation of the 5′untranslated region and codon usage bias in Mus musculus genes［J］.Gene,1998,431（3）:101-105.

［7］ FENG Chao,XU Changjie,WANG Yue,et al.Codon usage patterns in Chinese bayberry （Myrica rubra）based on RNA-Seq data［J］.BMC Genet,2013,14（6）:986-991.

［8］ PENG Zhenhua,LU Ying,LI Lubin,et al.The draft genome of the fast-growing non-timber forest species moso bamboo（Phyllostachys heterocycla）［J］.Nat Genet,2013,45（4）:456-461.

［9］ 劉慶坡，譚軍，薛慶中.秈稻品種93-11同義密碼子的使用偏好性［J］.遺傳學(xué)報(bào)，2003，30（4）：335-340.

LIU Qingpo,TAN Jun,XUE Qingzhong.Synonymous codon usage bias in the rice cultivar 93-11（Oryza sativa L.ssp. indica）［J］.Acta Genet Sin,2003,30（4）：335–340.

［10］ WRIGHT F.The “effective number of codons” used in a gene［J］.Gene,1990,87（1）:23-29.

［11］ SHARP P M，LI W H.An evolutionary perspective on synonymous codon usage in unicellular organisms［J］.J Mol Evol,1986,24（1/2）：28-38.

［12］ 時(shí)慧，王玉，楊路成,等.茶樹抗寒調(diào)控轉(zhuǎn)錄因子ICE1密碼子偏好性分析［J］.園藝學(xué)報(bào)，2012,39（7）：1341 -1352.

SHI Hui,WANG Yu,YANG Lucheng,et al.Analysis of codon bias of the cold regulated transcription factor ICE1 in tea plant［J］.Acta Hortic Sin,2012,39（7）:1341-1352.

［13］ STENICO M,LLOYD A T,SHARP P M.Codon usage in Caenorhabditis elegans:delineation of translational selection and mutational biases［J］.Science,2002,296（5576）:2174-2176.

［14］ KAWABE A,MIYASHITA N T.Patterns of codon usage bias in three dicot and four monocot plant species［J］. Genes Genet Syst,2003,78（5）:343-352.

［15］ 蔣瑋，呂貝貝，何建華，等.草菇密碼子偏好性分析［J］.生物工程學(xué)報(bào)，2014，30（9）：1424-1435.

JIANG Wei,Lü Beibei,HE Jianhua,et al.Codon usage bias in the straw mushroom Volvariella volvacea［J］.Chin J Biotech,2014,30（9）:1424-1435.

［16］ 劉漢梅，何瑞，張懷渝，等.玉米同義密碼子偏愛性分析［J］.農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，2010，18（3）：456-461.

LIU Hanmei,HE Rui,ZHANG Huaiyu,et al.Analysis of synonymous codon bias in maize［J］.J Agric Biotechnol, 2010,18（3）:456-461.

［17］ 石秀凡，黃京飛，柳樹群，等.人類基因同義密碼子偏好的特征以及與基因GC含量的關(guān)系［J］.生物化學(xué)與生物物理進(jìn)展，2002，29（3）：411-414.

SHI Xiufan,HUANG Jingfei,LIU Shuqun,et al.The features of synonymous codon bias and GC content relationship in human genes［J］.Prog Biochem Biophys,2002,29（3）:411-414.

［18］ 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 within-species diversity［J］.Nucl Acid Res,1988,16（17）:8207-8211.

［19］ WANG Liangjiang,ROOSSINCK M J.Comparative analysis of expressed sequences reveals a conserved pattern of optimal codon usage in plants［J］.Plant Mol Biol,2006,61（4/5）:699-710.

［20］ WANG Huaichun,HICKEY D A.Rapid divergence of codon usage patterns within the rice genome［J］.BMC Evolut Biol,2007,7（S1）:173-188.

［21］ LYNCH D B,LOGUE M E,BUTLER G,et al.Chromosomal G+C content evolution in yeasts:systematic interspecies differences,and GC-poor troughs at centromeres［J］.Gen Biol Evol,2010,2（1）:572-583.

［22］ EYRE-WALKER A,HURST LD.The evolution of isochors［J］.Nat Rev Genet,2001,2（7）:549-555.

［23］ MAHER A S,KHABAR K S A.UU/UA dinucleotide frequency reduction in coding regions results in increased mRNA stability and protein expression［J］.Mol Ther J Am Soc Gene Ther,2012,20（5）:954-959.

［24］ SUN Jingchun,CHEN Ming,XU Jinlin,et al.Relationships among stop codon usage bias,its context,isochores, and gene expression level in various eukaryotes［J］.J Mol Evol,2005,61（4）:437-444.

［25］ RAO Yousheng,WU Guozuo,WANG Zhangfeng,et al.Mutation bias is the driving force of codon usage in the Gallus gallus genome［J］.J Jpn Veter Med Ass,1986,39（6）:154-158.

［26］ 張樂，金龍國(guó)，羅玲，等.大豆基因組和轉(zhuǎn)錄組的核基因密碼子使用偏好性分析［J］.作物學(xué)報(bào)，2011，37（6）：965-974.

ZHANG Le,JIN Longguo,LUO Ling,et al.Analysis of nuclear gene codon bias on soybean genome and transcriptome［J］.Acta Agron Sin,2011,37（6）:965-974.

［27］ SHAO Zhuqing,ZHANG Yanmei,FENG Xueying,et al.Synonymous codon ordering:a subtle but prevalent strategy of bacteria to improve translational efficiency［J］.Plos One,2012,7（3）:e33547.doi:10.1371/journal.pone.0033547.

［28］ QIAN Wenfeng,YANG Jianrong,PEARSON N M,et al.Balanced codon usage optimizes eukaryotic translational efficiency［J］.Plos Genet,2012,8（3）:e1002603.doi:10.1371/journal.pgen.1002603.

［29］ KYOKO H T,MPANJA N,TADAYOSHI H,et al.High-level accumulation of recombinant miraculin protein in transgenic tomatoes expressing a synthetic miraculin gene with optimized codon usage terminated by the native miraculin terminator［J］.Plant Cell Rep,2011,30（1）:113-124.

Synonymous codon bias of Phyllostachys edulis

HUANG Xiaoyu,XU Zaien,GUO Xiaoqin
（The Nurturing Station for the State Key Laboratory of Subtropical Silviculture,Zhejiang A&F University,Lin’an 311300,Zhejiang,China）

Analysis of codon usage bias for different species,an important characteristic of genetic information transfer in organisms,is important for understanding the rules of genetic information transfer.To better understand the characteristics of Phyllostachys edulis,coding DNA sequences （CDS）of 26 103 proteins in this species were analyzed.The content of G+C at three positions of codons’the Effective number of codons（ENC）and frequency of synonym codon usage for genes were calculated the“optimal codons”were determined by Codon W,CHIPS and CUSP.Then,the frequency of codon usage for Ph.edulis with other organisms including five model value species （Drosophila melanogaster,Saccharomys cerevisiae,Escherichia coli,Arabidopsis thaliana,and Nicotiana tabacum）and three other Gramineae species （Oryza sativa,Zea mays,and Triticum aestivum）were compared.Results showed that the content of G+C at the first and third position of codons was much higher than that of the second positions,and genes preferred codons with C or G in the synonymous position.Meanwhile,26 codons,ending with G or C,were determined as the “optimal codons”.Frequency of codon usage also had fewer differences for the three Gramineae species than it did for the six model value species.To some extent,differences in the size of codon bias reflected the evolutionary relationships between species.These results provided the useful information for understanding the evolution of Ph.edulis ［Ch,1 fig.4 tab.29 ref.］

botany;Phyllostachys edulis;codon bias;optimal codon

S722；S795.7

A

2095-0756（2017）01-0120-09

2016-02-29；

2016-04-20

浙江省自然科學(xué)基金資助項(xiàng)目（Y307499，LY13C160011）；國(guó)家自然科學(xué)基金資助項(xiàng)目（30901155）；浙江省大學(xué)生新苗人才計(jì)劃項(xiàng)目（2015R412020）

黃笑宇，從事林學(xué)研究。 E-mail：593796107@qq.com。通信作者：郭小勤，副教授，博士，從事竹類植物遺傳育種等研究。E-mail：xqguo@zafu.edu.cn

浙 江 農(nóng) 林 大 學(xué) 學(xué) 報(bào)，2017，34（1）：129-136

Journal of Zhejiang A＆F University

10.11833/j.issn.2095-0756.2017.01.018