兩種梧桐葉綠體基因組密碼子使用偏性分析

陸奇豐　駱文華　黃至歡

摘要：? 為了分析美麗梧桐、云南梧桐葉綠體基因組密碼子的使用偏性，該研究通過篩選美麗梧桐、云南梧桐葉綠體基因組中各52條蛋白編碼序列，并利用CodonW、CUSP和SPSS軟件對(duì)其密碼子使用模式及偏性進(jìn)行了分析。結(jié)果表明：（1）美麗梧桐、云南梧桐的GC含量分別為38.12%、38.05%，表明葉綠體基因組內(nèi)富含A/T堿基。（2）有效密碼子數(shù)（ENC）范圍為36.91～56.46、36.55～58.04，表明多數(shù)密碼子偏性較弱。（3）相對(duì)同義密碼子（RSCU）分析顯示，RSCU>1的密碼子各有29個(gè)，其中28個(gè)以A、U結(jié)尾。（4）中性繪圖顯示，GC3與GC12的相關(guān)性不顯著，回歸曲線斜率分別為0.195和0.304，說(shuō)明密碼子偏好性主要受到自然選擇的影響。（5）ENC-plot分析中大部分基因分布于曲線的周圍和下方，ENC比值多分布于-0.04～0.10之間，表明突變會(huì)影響密碼子偏性的形成。此外，17、18個(gè)密碼子分別被鑒定為美麗梧桐、云南梧桐的最優(yōu)密碼子。以上結(jié)果說(shuō)明美麗梧桐、云南梧桐葉綠體基因組的密碼子使用偏性可能受選擇和突變共同作用，且使用模式較為相似，但具有一定的差異，可能與適應(yīng)環(huán)境的進(jìn)化機(jī)制有關(guān)。

關(guān)鍵詞： 美麗梧桐， 云南梧桐， 葉綠體基因組， 密碼子偏性， 最優(yōu)密碼子

中圖分類號(hào)：? Q943.2文獻(xiàn)標(biāo)識(shí)碼：? A

文章編號(hào)：? 1000-3142（2020）02-0173-11

Abstract：? In order to analyze the codon usage patterns of the chloroplast genome of Firmiana pulcherrima and F. major. Fifty-two genes selected from chloroplast genome （Cp genome）? of these two species? were analyzed using CodonW， CUSP and SPSS software. The results showed that： （1） GC content in Firmiana pulcherrima was 38.12%， and that in F. major was 38.05%， suggesting that A/T base was rich in chloroplast genome. （2） The effective number of codons in F. pulcherrima ranged from 36.91 to 56.46， and that in F. major ranged from 36.55 to 58.04， meaning the codon bias was weak in F. pulcherrima and F. major. （3） Both F. pulcherrima and F. major? had 29 codons with relative synonymous codon usage （RSCU） greater than 1 and 28 codons ending with A and U base. （4） Neutrality plot showed that there was not significant correlation between GC3 and GC12， and the regression slope was 0.195 in F. pulcherrima， and 0.304 in F. major， suggesting that codon bias was mainly affected by natural selection. （5） ENC-plot analysis showed that most of genes fell below and followed the standard curve， and ENC ratio of most genes ranged from -0.04 to 0.10， suggesting that the mutation also affected the formation of codon bias in F. pulcherrima and F. major. In addition， 17 codons were identified as optimal codons in F. pulcherrima， and 18 codons in F. major. We suggest that the codon usage bias of F. pulcherrima and F. major are formed under effect of mutation and selection. Although the codon usage pattern of F. pulcherrima and F. major are very similar， there are some differences， which may be related to the evolutionary me-chanism of adaptive environment.

Key words： Firmiana pulcherrima， F. major， chloroplast genome， codon usage bias， optimal codons

自然界中，基因序列根據(jù)三聯(lián)密碼子法則翻譯成氨基酸序列，并進(jìn)一步組裝修飾形成具有特定空間結(jié)構(gòu)和功能活性的蛋白質(zhì)。密碼子由3個(gè)相連的核苷酸堿基所組成，一共存在64個(gè)密碼子，除去3個(gè)終止密碼子 （UAA、UAG、UGA），61個(gè)密碼子對(duì)應(yīng)特定的氨基酸。天然氨基酸有20種，除了甲硫氨酸（AUG）、色氨酸（UGG）對(duì)應(yīng)單一密碼子編碼以外，其余18種氨基酸還均對(duì)應(yīng)2～6個(gè)編碼密碼子。編碼相同氨基酸的不同密碼子又稱為同義密碼子，同義密碼子使用頻率在基因內(nèi)不是隨機(jī)的，而是優(yōu)先使用其中一些密碼子，這一現(xiàn)象被稱為密碼子使用偏性（Gustafsson et al.， 2004;Li et al.， 2017）。影響密碼子偏性的因素較多，如自然選擇下DNA的變異（Ashesh， 2002;Suzuki， 2011）、tRNA豐度（Buchan et al.， 2006）、基因序列長(zhǎng)度（Gabriel et al.， 2001;Stoletzki et al.， 2011;Ribeiro et al.， 2012）、蛋白質(zhì)翻譯效率（Haupt et al.， 2009;Li et al.， 2017）等。研究密碼子偏性不僅有助于解釋物種間進(jìn)化規(guī)律，而且還對(duì)基因表達(dá)、載體構(gòu)建、未知功能基因的分析具有重要意義。

葉綠體屬于半自主性細(xì)胞器，是進(jìn)行光合作用和固定CO2的主要場(chǎng)所，參與植物生長(zhǎng)環(huán)境脅迫（干旱、鹽、光等）相關(guān)的重要代謝過程（Liu et al.， 2018）。此外，葉綠體具有母系遺傳、結(jié)構(gòu)穩(wěn)定和基因組較小等特性，從而廣泛運(yùn)用于植物多樣性、系統(tǒng)發(fā)育、DNA 條碼開發(fā)、植物適應(yīng)性和遺傳工程研究中（Sajjad et al.， 2017;Liu et al.， 2018）。不同物種葉綠體基因組密碼子使用偏性存在一定差異（Liu & Xue， 2005;鄭薇瑋等，2016）。密碼子使用模式對(duì)基因功能和表達(dá)有一定的影響作用，基因的密碼子偏性越強(qiáng)，其表達(dá)水平越高（Zhou et al.，2016）。因此，為提高葉綠體遺傳轉(zhuǎn)化效率及目的基因表達(dá)水平，很有必要開展葉綠體全基因組密碼子使用模式的研究。

美麗梧桐（Firmiana pulcherrima）、云南梧桐（F. major）是梧桐屬內(nèi)的重要物種，具有較高的研究?jī)r(jià)值和經(jīng)濟(jì)價(jià)值。但是，由于生態(tài)環(huán)境的破壞，美麗梧桐、云南梧桐野外資源日益減少，已被列為瀕危樹種（中國(guó)科學(xué)院，1984;余志祥等，2013;周亞東等，2013）。目前，對(duì)美麗梧桐、云南梧桐的研究主要集中于資源調(diào)查（黎國(guó)運(yùn)等，2011;Yang et al.， 2018）、種苗繁育（黎國(guó)運(yùn)等，2011;余志祥等，2013;周亞東等，2013）等方面，兩種梧桐葉綠體全基因組雖已完成測(cè)序 （Wang et al.， 2018;Ya et al.， 2018），但關(guān)于兩者的葉綠體基因組編碼基因的密碼子使用模式比較的研究卻尚未見有報(bào)道。因此，本研究在前人研究的基礎(chǔ)上，對(duì)美麗梧桐和云南梧桐的葉綠體基因組中52個(gè)基因的密碼子堿基組成特點(diǎn)及使用偏性進(jìn)行比較分析，揭示其密碼子偏性的影響因素，并分析最優(yōu)密碼子，旨在為美麗梧桐、云南梧桐制定全面的保育策略以及為梧桐科葉綠體基因組學(xué)研究提供參考。

1材料與方法

1.1 葉綠體基因序列數(shù)據(jù)

從NCBI 分別檢索下載美麗梧桐、云南梧桐葉綠體全基因組序列信息，登錄號(hào)分別為NC_036395.1（F. pulcherrima 美麗梧桐）、NC_037242.1（F. major 云南梧桐），分別獲得82條和85條編碼基因。對(duì)所獲得的編碼基因按以下條件進(jìn)行篩選：基因長(zhǎng)度大于300 bp，以ATG為起始密碼子，以TAA、TAG、TGA為終止密碼子，并剔除重復(fù)的編碼基因，最終獲得各52條編碼基因進(jìn)行后續(xù)分析。

1.2 方法

1.2.1 密碼子相關(guān)參數(shù)的計(jì)算分別利用軟件CodonW（http：//codonw.sourceforge.net/）和CUSP（http：//imed.med.ucm.es/EMBOSS/）分析CDS的有效密碼子數(shù)（ENC）、密碼子第3位堿基組成特性，以CT3、CC3、CA3、CG3表示，統(tǒng)計(jì)密碼子第1位、第2位、第3位堿基的GC含量及52個(gè)基因總GC含量，分別用CGC1、CGC2、CGC3、GC表示。

1.2.2 中性繪圖分析計(jì)算各CDS的CGC1、CGC2的平均值，記為CGC12，以此為縱坐標(biāo)，以CGC3為橫坐標(biāo)繪制散點(diǎn)圖，分析CGC3與CGC12的相關(guān)性，并推測(cè)影響密碼子偏性的因素。若CGC3與CGC12呈顯著相關(guān)，表明密碼子對(duì)堿基使用無(wú)差異，突變是密碼子偏性的主要影響因素;若CGC3與CGC12相關(guān)性不顯著，且斜率接近0，表明密碼子第1位、第2位堿基與第3位堿基使用存在差異，密碼子偏性與第3位堿基顯著相關(guān)，選擇壓力對(duì)密碼子偏性的影響較強(qiáng)（尚明照等，2011; 秦政等，2018）。

1.2.3 PR2-plot繪圖分析PR2-plot分析是考察突變和自然選擇對(duì)密碼子偏性影響的方法之一（尚明照等，2011;劉慧等，2017）。為避免密碼子第3位堿基上的A/T和G/C突變的不平衡，通過統(tǒng)計(jì)兩種梧桐葉綠體52條基因中編碼丙氨酸、甘氨酸、亮氨酸（CTT、CTC、CTA、CTG）、脯氨酸、精氨酸（CGT、CGC、CGA、CGG）、絲氨酸（TCT、TCC、TCA、TCG）、蘇氨酸、纈氨酸的密碼子第3位堿基的T3、C3、A3、G3的組成情況，并以A3/（A3+T3）為縱坐標(biāo)，以G3/（G3+C3）為橫坐標(biāo)繪制散點(diǎn)圖。

1.2.4 ENC-plot繪圖分析為進(jìn)一步了解兩種梧桐堿基組成對(duì)密碼子偏性的影響，以兩種梧桐葉綠體中各基因的ENC、GC3值為縱、橫坐標(biāo)做散點(diǎn)圖，并與ENC期望值（ENC =2+GC3+29/（[GC32+ （1-GC3） 2]） 進(jìn)行比較，期望值代表密碼子偏性僅受堿基組成的影響。若基因分布于期望值曲線的上方，其密碼子偏性僅受突變影響;若基因分布于曲線下方較遠(yuǎn)位置，表明選擇是影響密碼子偏性的主要因素（牛元等，2018）。

1.2.5 相對(duì)同義密碼子使用度分析為考察同義密碼子偏好性情況，以相對(duì)同義密碼子使用度（RSCU）表示特定密碼子在編碼對(duì)應(yīng)氨基酸的同義密碼子間的相對(duì)使用概率（楊國(guó)鋒等，2015）。當(dāng)RSCU>1時(shí)，表示該密碼子使用頻率較高，偏好性較強(qiáng);當(dāng)RSCU<1時(shí)，表示該密碼子使用頻率較低，偏好性較弱;當(dāng)RSCU=1時(shí)，該密碼子無(wú)偏好性。

1.2.6 最優(yōu)密碼子分析對(duì)CDS的ENC值進(jìn)行排序，從兩極各選出5個(gè)基因作為高、低表達(dá)組，分別計(jì)算對(duì)應(yīng)的密碼子RSCU值，篩選高低表達(dá)組內(nèi)對(duì)應(yīng)密碼子△RSCU之差大于0.08作為優(yōu)越高表達(dá)密碼子，并與整體密碼子RSCU進(jìn)行結(jié)合比較，將△RSCU≥0.08，且整體RSCU>1的密碼子定義為最優(yōu)密碼子（秦政等，2018;石元豹等，2018 ）。

2結(jié)果與分析

2.1 密碼子組成分析

通過CodonW和在線CUSP對(duì)美麗梧桐、云南梧桐葉綠體編碼基因組的52條CDS進(jìn)行分析，結(jié)果如表1所示，美麗梧桐、云南梧桐葉綠體編碼基因密碼子CT3和CA3的含量分別為46.56%、46.67%與43.04%、43.17%，表明兩種梧桐葉綠體編碼基因第3位堿基以A/T結(jié)尾為主。密碼子上不同位置堿基GC含量有所差異，第1位堿基的GC含量最高，分別為46.10%、46.17%;其次為第2位，分別為38.14%、38.06%;最低為第3位堿基，分別為30.11%、29.92%，且兩種梧桐葉綠體基因組GC含量分別為38.12%、38.05%，表明葉綠體編碼基因富含A/T堿基。ENC值可反映密碼子偏性程度，其取值范圍為20～61，ENC值越小表明偏性越強(qiáng)，ENC值越大則密碼子偏性越弱，一般以ENC的值為35時(shí)區(qū)分密碼子偏性強(qiáng)弱（Jiang et al.，2008）。本研究中，兩種梧桐葉綠體基因組的ENC值為50左右，大于35表明其密碼子偏性較弱。

利用在線CUSP對(duì)各基因密碼子不同位置堿基的GC含量進(jìn)行分析，結(jié)果如表2所示，不同基因的GC含量存在一定的差異，CGC1、CGC2與CCG3差異較大，兩種梧桐葉綠體編碼基因ENC范圍分別為36.91～56.46、36.55～58.04，且大部分基因的ENC大于44，表明葉綠體編碼基因密碼子偏性較弱。美麗梧桐、云南梧桐中atpE基因的ENC值最大，分別為56.46和58.04，說(shuō)明美麗梧桐中編碼atpE基因密碼子偏性強(qiáng)于云南梧桐，同時(shí)兩種梧桐葉綠體編碼基因GC含量、ENC值之間也存在差異，但差異不大。此外，本研究還發(fā)現(xiàn)，美麗梧桐中不存在clpP基因及云南梧桐中的ndhK基因?yàn)榧倩颍砻鲀煞N梧桐葉綠體基因組具有較高的保守性，但種間卻存在一定的差異。

CGC1、CGC2、CCG3、GC、ENC及序列長(zhǎng)度（sequence length， SL）相關(guān)性分析結(jié)果如表3所示，CGC1與CGC2極顯著相關(guān)，CGC1、CGC2均與CCG3相關(guān)性不顯著，表明兩種梧桐葉綠體CDS密碼子的第1位、第2位堿基組成相似，且與第3位堿基組成存在較大差異，GC與CGC1、CGC2、CCG3均呈極顯著相關(guān)，SL與CGC2呈負(fù)相關(guān)。此外，美麗梧桐的ENC與CCG3呈顯著相關(guān)，與CGC1、CGC2、GC、SL均未達(dá)到顯著水平，表明美麗梧桐葉綠體編碼基因密碼子第3位堿基組成對(duì)密碼子偏性影響較大。云南梧桐中ENC與CCG1、CGC3顯著相關(guān)，相關(guān)值分別為0.29、0.32，與CGC2、GC、SL相關(guān)性較低，表明CCG3對(duì)云南梧桐葉綠體密碼子偏性的影響強(qiáng)于CGC1。兩種梧桐密碼子參數(shù)的相關(guān)性分析結(jié)果存在差異，其中 CGC1與云南梧桐葉綠體基因組密碼子偏性顯著相關(guān)，而美麗梧桐葉綠體基因組密碼子偏性與CGC1相關(guān)不顯著，進(jìn)一步說(shuō)明兩種植物葉綠體基因組密碼子偏性的影響因素有所差異。

2.2 中性繪圖分析

如圖1所示，美麗梧桐的葉綠體CDS的CGC3與CGC12含量分布范圍分別為0.223 0～0.370 5、0.329 2～0.568 4，云南梧桐為0.223 0～0.369 7、0.328 8～0.564 8，兩者的分布范圍幾乎重合，說(shuō)明兩種梧桐葉綠體基因組較為保守。兩種梧桐的CGC3與CGC12的Pearson相關(guān)系數(shù)分別為0.121和0.206，進(jìn)一步說(shuō)明CGC1和CGC2均與CCG3相關(guān)性不顯著，密碼子的第1位、第2位堿基與第3位堿基組成存在顯著差異，密碼子偏性與第3位堿基顯著相關(guān)。同時(shí)，回歸曲線斜率分別為0.195和0.304，表明選擇壓力對(duì)兩種梧桐屬葉綠體基因組密碼子偏性的影響較強(qiáng)。此外，突變對(duì)密碼子偏性也有一定的影響，且云南梧桐葉綠體密碼子偏性受到突變的影響高于美麗梧桐。

2.3 PR2-plot繪圖分析

PR2-plot繪圖分析是指通過對(duì)密碼子第3位堿基組成進(jìn)行分析，探討突變和自然選擇對(duì)密碼子使用偏性的影響，當(dāng)密碼子第3位堿基表現(xiàn)為A=T和C=G時(shí)，表明DNA兩條互補(bǔ)鏈上不存在選擇效應(yīng)。突變是影響密碼子偏性的唯一因素（Sueoka， 2001;Zhang et al.， 2018）。本研究中，兩種梧桐大部分基因分布在左下方區(qū)域（圖2），T、C堿基頻率高于A、G堿基，說(shuō)明密碼子第3位堿基偏愛使用嘧啶堿基，表明兩種梧桐葉綠體編碼基因密碼子使用偏性不僅受到突變的影響，而且還受到其他因素的影響，如自然選擇壓力等。

2.4 ENC-plot繪圖分析

從圖3可以看出，兩種梧桐的大部分基因較集中分布于標(biāo)準(zhǔn)曲線的兩側(cè)，說(shuō)明其ENC值與ENC期望值較為接近，而rps8、rps14、rps18、petD、rpl16、psbA基因偏離標(biāo)準(zhǔn)曲線較遠(yuǎn)，且ENC值較小，密碼子偏性較強(qiáng)。根據(jù)公式（ENC期望-ENC實(shí)際）/ENC期望統(tǒng)計(jì)基因ENC的分布頻數(shù)（表4），結(jié)果顯示大部分基因的ENC頻數(shù)集中于-0.04～0.10之間，進(jìn)一步表明大部分基因的ENC值與ENC期望值差異較小，表明這部分基因密碼子偏性主要受到突變的影響。rps8、rps14、rps18、petD、rpl16、psbA分布于0.18～0.33之間，與0點(diǎn)相差較遠(yuǎn)，表明這部分基因密碼子偏性主要受到選擇壓力的影響。

2.5 RSCU分析

利用CondonW程序?qū)煞N梧桐葉綠體基因組52個(gè)編碼基因進(jìn)行RSCU分析，結(jié)果如表5所示，兩種梧桐的高頻密碼子（RSCU>1）各有29個(gè)，除了UUG以外，其余密碼子均以A/U結(jié)尾，表明兩種梧桐密碼子偏好使用以A/U結(jié)尾的密碼子，且兩種梧桐植物對(duì)應(yīng)密碼子RSCU值極為接近，進(jìn)一步說(shuō)明其葉綠體基因組密碼子使用偏性較為接近，葉綠體基因組較為保守，種間差異較小。

2.6 最優(yōu)密碼子分析

對(duì)兩種梧桐植物葉綠體基因組基因的ENC值進(jìn)行排序，選出兩極的10%基因作為高、低表達(dá)基因庫(kù)，其中兩種梧桐屬葉綠體基因組高表達(dá)基因?yàn)閞ps8、rps14、petD、rps18、rpl16，平均ENC值分別為38.41（F. pulcherrima）、38.02（F. major）;F. pulcherrima低表達(dá)基因?yàn)閥cf2、petA、 ycf3、ndhJ、atpE，圖 3ENC-plot繪圖分析

Fig. 3Analysis of ENC-plot （ENC value against CGC3）而F. major為ycf2、clpP、petA、ycf3、ndhJ、atpE，平均ENC值分別為55.21（F. pulcherrima）、55.24（F. major）。分別計(jì)算高、低表達(dá)基因庫(kù)的RSCU值，并計(jì)算對(duì)應(yīng)密碼子的△RSCU值（高表達(dá)庫(kù)RSCU值-低表達(dá)庫(kù)RSCU值），同時(shí)與整體密碼子的RSCU值結(jié)合比較，篩選出△RSCU值≥0.08，且整體RSCU>1的密碼子，定義為優(yōu)越密碼子。最終分別確定了17、18個(gè)表達(dá)優(yōu)越密碼子（表4），分別為GCA、UGU、GAA、UUU、GGA、GGU、CAU、AUU、AAA、UUA、CAA、CGA、CGU、AGU、ACA、ACU（為F. major優(yōu)越密碼子）、GUA、UAU，這些優(yōu)越密碼子均以A/U結(jié)尾。此外， 我們觀察到兩種梧桐植物葉綠體基因組使用的優(yōu)越密碼子基本相同，但其△RSCU值存在一定的差異，F(xiàn). pulcherrima 中的GAA、 CGA密碼子△RSCU>0.08且<0.30， 而F.

3討論與結(jié)論

密碼子是遺傳信息的基本單位，是核酸與蛋白質(zhì)的聯(lián)系紐帶，研究其使用情況對(duì)了解物種基因組進(jìn)化具有重要意義（李平等，2011）。密碼子使用偏性是物種不斷進(jìn)化適應(yīng)的結(jié)果，其形成機(jī)制受多種因素的影響，其中突變作用和自然選擇壓力是其重要影響因素（王鵬良等， 2018）。不同物種、基因家族中密碼子使用偏性存在一定的差異，因此密碼子使用偏性是研究物種基因組進(jìn)化關(guān)系及影響的重要指標(biāo)之一。

堿基組成差異對(duì)密碼子的使用有重要的影響。密碼子在第1位、第2位堿基上的突變通常會(huì)導(dǎo)致編碼氨基酸的改變，第3位堿基存在簡(jiǎn)并性，其突變對(duì)編碼氨基酸的改變通常無(wú)影響（尚明照等，2011）。因此，密碼子第3位堿基的使用模式與密碼子偏性密切相關(guān)。有研究表明，植物葉綠體基因組中密碼子第3位堿基偏好使用A/T（Zhou et al.， 2008;秦政等，2018 ）。本研究中，兩種梧桐植物CT3和CA3高于CG3、CC3，且CGC1、CGC2均與CCG3相關(guān)性不顯著，CCG3比例最低，說(shuō)明梧桐屬植物葉綠體基因組密碼子偏好使用A/T堿基結(jié)尾，與大花香水月季（Rosa odorata var. gigantea）、表 5兩種梧桐密碼子RSCU及最優(yōu)密碼子樟樹（Cinnamomum camphora）、甜蕎（Fagopyrum esculentum） （羅洪等，2015;牛元等，2018;秦政等，2018）等植物的研究結(jié)果一致。

自然選擇和突變是影響密碼子使用偏性的重要因素（Sharp et al.，2010;Rao et al.，2011）。本研究通過對(duì)兩種梧桐葉綠體基因組52個(gè)編碼基因CGC3與CGC12相關(guān)分析表明，密碼子第3位堿基GC含量與第1位、第2位堿基相關(guān)性弱，且線性回歸斜率較小，表明兩種梧桐密碼子的使用受到選擇壓力因素較強(qiáng)。此外，PR2-plot分析發(fā)現(xiàn)，密碼子第3位堿基上A/T堿基和C/G堿基使用頻率存在差異，更偏好使用T、C堿基，進(jìn)一步表明兩種梧桐葉綠體基因組密碼子使用不僅與堿基突變的有關(guān)，而且還更多受到選擇的影響。ENC-plot分析顯示：大部分基因的ENC值與期望ENC值較為接近，其分布頻數(shù)集中于-0.04～0.10之間，表明這部分基因受密碼子偏性與第3位堿基的GC含量相關(guān)，突變是主要影響因素;而rps8、rps14、rps18、petD、rpl16、psbA分布于標(biāo)準(zhǔn)曲線的下方較遠(yuǎn)位置，自然選擇影響其密碼子偏性的主要因素。綜合中性繪圖分析、ENC-plot、PR2-plot等分析發(fā)現(xiàn)，突變和自然選擇等多重因素共同影響著兩種梧桐葉綠體密碼子使用偏性。這與柿屬植物、普通油茶、苦蕎、甜蕎等物種的葉綠體基因組密碼子偏性研究結(jié)果相似（羅洪等，2015;胡莎莎等，2016;傅建敏等，2017;王鵬良等，2018）。

本研究發(fā)現(xiàn)，兩種梧桐中29個(gè)高頻密碼子（RSCU>1）第3位堿基主要以A/U結(jié)尾，與前人在其他植物中的研究結(jié)果相似（尚明照等，2011;秦政等，2018）。與高表達(dá)密碼子進(jìn)行綜合分析，最終確定美麗梧桐的最優(yōu)密碼子為17個(gè)，云南梧桐18個(gè)，分別為GCA、UGU、GAA、UUU、GGA、GGU、CAU、AUU、AAA、UUA、CAA、CGA、CGU、AGU、ACA、ACU（為F. major優(yōu)越密碼子）、GUA、UAU，這些最優(yōu)密碼子以A/U結(jié)尾，與大多植物葉綠體基因組最優(yōu)密碼子相似，表明梧桐與其他雙子葉植物葉綠體基因組密碼子第3位堿基偏愛使用A/U結(jié)尾（Kawabe & Miyashita， 2003）。有研究表明，物種間的親緣關(guān)系越近，其密碼子使用模式越相似（張?zhí)龋?017;張孟偉等，2018）。本研究中，兩種梧桐葉綠體基因組密碼子使用偏性極為相似，但存在一定的差異。此外，本研究還發(fā)現(xiàn)，美麗梧桐中葉綠體clpP編碼基因丟失，云南梧桐的ndhK為假基因，推測(cè)兩種梧桐適應(yīng)環(huán)境改變的進(jìn)化趨勢(shì)有所不同，有待深入研究。

本研究分析了美麗梧桐、云南梧桐葉綠體基因組編碼基因的密碼子偏性形成的影響因素和特點(diǎn)，以及兩種梧桐葉綠體基因組的最優(yōu)密碼子，為外源基因密碼子改造及梧桐屬葉綠體基因組進(jìn)化機(jī)制研究提供一定的理論基礎(chǔ)。

參考文獻(xiàn)：

ASHESH N， 2002. Investigations on evolutionary changes in base distributions in gene sequences? [J]. Internet Electr J Mol Design， 1： 545-558.

BUCHAN JR， AUCOTT LS， STANSFIELD I， 2006. tRNA properties help shape codon pair preferences in open reading frames? [J]. Nucl Acids Res， 34（3）： 1015-1027.

Chinese Academy of Sciences， 1984. Flora Reipublicae Popularis Sinicae [M]. http： //frps.iplant.cn/， 49（2）： 136-137. [中國(guó)科學(xué)院， 1984. 中國(guó)植物志 [M]. http： //frps.iplant.cn/， 49（2）： 136-137.]

FU JM， SUO YJ， LIU HM， et al.， 2017. Analysis on codon usage in the chloroplast protein coding genes of Diospyros spp.? [J]. Nonw For Res， 35（2）：38-44.? [傅建敏， 索玉靜， 劉慧敏， 等， 2017. 柿屬植物葉綠體蛋白質(zhì)編碼基因密碼子用法 [J]. 經(jīng)濟(jì)林研究， 35（2）： 38-44.]

GABRIEL M， LAURENT D，2001. Synonymous codon usage， accuracy of translation， and gene length in Caenorhabditis elegans? [J]. J Mol Evol， 52： 275-280.

GUSTAFSSON C， GOVINDARAJAN S， MINSHULL J，2004. Codon bias and heterologous protein expression? [J]. Trends Biotechnol， 22（7）： 346-253.

HAUPT S， ZIEGLER A， COWAN G，et al.， 2009. Studies of the role and function of barley stripe mosaic virus encoded proteins in replication and movement using GFP fusions? [M]. Methods Mol Biol， 515： 287-297.

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

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

KAWABE A， MIYASHITA NT， 2003. Patterns of codon usage bias in three dicot and four moncot plant species? [J]. Genes Genet Syst， 78（5）： 343-352.

LI GY， SHENG XH， XU PL， 2011. A preliminary report on wild resource survey of Firmiana pulcherrima? [J]. Torp For， 39（1）： 50-52.? [黎國(guó)運(yùn)， 盛小彬， 徐佩玲， 2011. 美麗梧桐野生資源調(diào)查初報(bào) [J]. 熱帶林業(yè)， 39（1）： 50-52.]

LI GY， GUAN XJ， SHENG XB， et al.， 2011. Research on shoot cutting propagation of Firmiana pulcherrima Hsue? [J]. Torp For， 39（4）： 36-38.? [黎國(guó)運(yùn)， 關(guān)秀娟， 盛小彬， 等， 2011. 美麗梧桐扦插育苗技術(shù)研究 [J]. 熱帶林業(yè)， 39（4）： 36-38.]

LI J， LI HY， ZHI JK， et al.， 2017. Codon usage of expansin genes in Populus trichocarpa? [J]. Current Bioinform， 12： 452-461.

LI P， BAI YF， FENG RY， et al.， 2011. Analysis of codon bias of NAD-ME gene in Amaranthus hypochondriacus? [J]. Chin J Appl Environ Biol， 1： 12-17.? [李平， 白云鳳， 馮瑞云， 等， 2011. 籽粒莧蘋果酸酶（NAD-ME）基因密碼子偏好性分析 [J]. 應(yīng)用與環(huán)境生物學(xué)報(bào)， 1： 12-17.]

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

LIU X， ZHOU BY， YANG HY， et al.， 2018. Sequencing and analysis of Chrysanthemum carinatum Schousb and Kalimeris indica. The complete chloroplast genomes reveal two inversions and rbcL as barcoding of the vegetable? [J]. Molecules， 23（6）： 1358.

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

NIU Y， XU Q， WANG YD， et al.， 2018. An analysis on codon usage bias of chloroplast genmoe of Rosa odorata var. gigantea? [J]. J NW For Univ， 33（3）： 123-130.? [牛元， 徐瓊， 王崳德， 等， 2018. 大花香水月季葉綠體基因組密碼子使用偏性分析 [J]. 西北林學(xué)院學(xué)報(bào)， 33（3）： 123-130.]

QIN Z， ZHENG YJ， GUI LJ， et al.， 2018. Codon usage bias analysis of the chloroplast genome of camphor tree （Cinnamomum camphora）? [J]. Guihaia， 38（10）： 1346-1355.? [秦政， 鄭永杰， 桂麗靜， 等， 2018. 樟樹葉綠體基因組密碼子偏好性分析 [J]. 廣西植物， 38（10）： 1346-1355.]

RAO Y， WU G， WANG Z， et al.， 2011. Mutation bias is the driving force of codon usage in the gallus gallus genome? [J]. DNA Res， 18（6）： 499-512

RIBEIRO AS， HKKINEN A， LLOYD-PRICE J， 2012. Effects of gene length on the dynamics of gene expression? [J]. Comput Biol Chem， 41（4）： 1-9.

SAJJAD A， MUHAMMAD W， ABDUL LK， et al.， 2017. The complete chloroplast genome of wild rice （Oryza minuta） and its comparison to related species? [J]. Front Plant Sci， 8： 304.

SHANG MZ， LIU F， HUA JP， et al.， 2011. Analysis on codon usage of chloroplast genome of Gossypium hirsutum? [J]. Sci Agric Sin， 44（2）：245-253.? [尚明照， 劉方， 華金平， 等， 2011. 陸地棉葉綠體基因組密碼子使用偏性的分析 [J]. 中國(guó)農(nóng)業(yè)科學(xué)， 44（2）：245-253.]

SHARP PM EMERY LR， ZENG K， 2010. Forces that influence the evolution of codon bias? [J]. Phil Trans R Soc B， 365（1544）： 1203-1212.

SHI YB， WANG GB， YANG XM， et al.， 2018. Analysis of codon usage of WRKY transcription factors in Ginkgo biloba? [J]. Mol Plant Breed， http：kns.cnki.net/kcms/detail/46.1068.S.20180809.1423.008.html.? [石元豹， 汪貴斌， 楊曉明， 等， 2018. 銀杏WRKY家族基因密碼子使用偏向性分析 [J]. 分子植物育種， http：//kns.cnki.net/kcms/detail/46.1068.S.20180809.1423.008.html.]

STOLETZKI N， 2011. The surprising negative correlation of gene length and optimal codon use-disentangling translational selection from GC-biased gene conversion in yeast? [J]. BMC Evol Biol， 11（1）： 1179-1187.

SUEOKA N， 2001. Near homogeneity of PR2-bias fingerprints in the human genome and their implications in phyloge-netic analyses? [J]. J Mol Evol， 53： 469 –476.

SUZUKI Y， 2011. Statistical methods for detecting natural selection from genomic data? [J]. Genes Genet Syst， 85（6）： 359-376.

WANG JH， CAI YC， ZHAO KK， et al.， 2018. Characterization of the complete chloroplast genome sequence of Firmiana pulcherrima （Malvaceae）? [J]. Conserv Genet Resour， 10（3）： 445-448.

WANG PL， YANG LP， WU HY， et al.， 2018. Codon preference of chloroplast genome in Camellia oleifera? [J]. Guihaia， 38（2）： 135-144.? [王鵬良， 楊利平， 吳紅英， 等， 2018. 普通油茶葉綠體基因組密碼子偏好性分析 [J]. 廣西植物， 38（2）： 135-144.]

YA JD， YU ZX， YANG YQ， et al.， 2018. Complete chloroplast genome of Firmiana major （Malvaceae）， a critically endangered species endemic to Southwest China? [J]. Conserv Genet Resour， 10（4）： 713-715.

YANG GF， SU KL， ZHAO YR， et al.， 2015. Analysis of condon usage in the chloroplast genome of Medicago truncatula? [J]. Atca Pratac Sin， 24（12）： 171-179.? [楊國(guó)鋒， 蘇昆龍， 趙怡然， 等， 2015. 蒺藜苜蓿葉綠體密碼子偏好性分析 [J]. 草業(yè)學(xué)報(bào)， 24（12）： 171-179.]

YANG J， CHEN G， SUN WB， 2018. Conserving Firmiana major， a tree species endemic to China? [J]. Oryx， 52（2）：? 211-213.

YU ZX， YANG YQ， MO X， et al.， 2013. A shudy of artificial breeding of an endangered panlt （Firmiana major）? [J]. J Sichuan For Sci Technol， 34（2）： 62-65.? [余志祥， 楊永瓊， 莫旭， 等， 2013. 瀕危植物云南梧桐（Firmiana major）人工繁育研究 [J]. 四川林業(yè)科技， 34（2）： 62-65.]

ZHANG H， LI J， WANG RX， et al.， 2018. Comparative analysis of expansin gene codonusage patterns among eight plant species [J]. J Biomol Struct Dyn， 7： 1-8.

ZHANG MW， ZHANG TK， LIU CY， et al.， 2018. Biosysthsis pathway evolution and codon usage bias analysis of lignin in Punica granatunm? [J]. Mol Plant Breed， 16（7）： 2131-2138.? [張孟偉， 張?zhí)?劉翠玉， 等， 2018. 石榴木質(zhì)素合成途徑進(jìn)化與密碼子使用偏好性分析 [J]. 分子植物育種， 16（7）： 2131-2138.]

ZHANG TK， QI GH， YE HL， et al.， 2017. Codon usage bias in pomegranate transcriptome? [J]. Acta Hortic Sin， 44 （4）： 675-690.? [張?zhí)?起國(guó)海， 葉紅蓮， 等， 2017. 石榴轉(zhuǎn)錄組密碼子使用偏向性 [J]. 園藝學(xué)報(bào)， 44 （4）： 675-690.]

ZHENG WW， CHEN JH， HAO ZD， et al.， 2016. Analysis on codon usage bias of chloroplast genes from five conifers including Cunninghamia lanceolata? [J]. Mol Plant Breed， 14（5）： 1091-1097.? [鄭薇瑋， 陳金慧， 郝兆東， 等， 2016. 杉木等5種針葉樹葉綠體密碼子偏好性分析 [J]. 分子植物育種， 14（5）： 1091-1097.]

ZHOU YD， ZHOU YB， YAN SL，et al.， 2013. Effects of different substrates on seed germination and seeding growth of Firmiana pulcherrima Hsue? [J]. J Trop Biol， 4（4）： 322-326.? [周亞東， 周有標(biāo)， 顏速亮， 等， 2013. 不同基質(zhì)對(duì)美麗梧桐種子育苗的影響 [J]. 熱帶生物學(xué)報(bào)， 4（4）： 322-326.]

ZHOU ZP， DANG YK， ZHOU M， et al.， 2016. Codon usage is an important determinant of gene expression levels largely through its effects on transcription [J]. Proc Natl Acad Sci USA， 26： e6117-e6125.

（責(zé)任編輯 蔣巧媛）