趙 英，付海天，金 剛，徐 釧，宋恩亮，羅燕春，莫永龍

（廣西壯族自治區(qū)亞熱帶作物研究所，廣西 南寧 530001）

木薯（Manihot esculenta Crantz）起源于南美洲亞馬遜河流域，是世界三大薯類作物之一，有“地下糧倉”和“淀粉之王”的美稱［1］，肩負著世界近8億人口的口糧任務［2］，具有高光效、高生物量、耐干旱貧瘠、易種植等特點。木薯除可以作為糧食外，也被用作飼料和重要的工業(yè)原料［3］，是我國南方熱帶亞熱帶旱地坡地主要的經(jīng)濟作物之一［4］。當前，對于木薯的研究相對其他作物而言還比較薄弱，選育種和高產(chǎn)高效栽培是其研究的主要方向［5］；木薯淀粉積累的機制尚不明確［6］；其功能基因組學的研究也相對滯后，限制了木薯遺傳育種的發(fā)展［7］。

將DNA或RNA序列以3個核苷酸為一組的密碼子翻譯為蛋白質的氨基酸序列，稱為遺傳密碼，又稱密碼子或者三聯(lián)體，其具有簡并性；編碼同一氨基酸的密碼子稱為同義密碼。密碼子使用頻率的不同形成了密碼子的偏好性，不同基因組有其特有的密碼子使用偏好性，表現(xiàn)在物種進化過程中受到基因突變壓力和自然選擇的影響不同［8］。在特定物種或基因中使用頻率最高的密碼子稱為最優(yōu)密碼子。研究目標物種的密碼子使用偏好性，根據(jù)該物種的最優(yōu)密碼子對外源基因的表達載體進行優(yōu)化設計，對提高外源基因組在葉綠體中的表達量具有重要的作用［9］。

葉綠體基因組屬單親遺傳，突變率低，進化速度慢［10］。對葉綠體基因組進化的研究可加深對植物間親緣關系、葉綠體基因組及核基因組間信息交流等的了解，并為設計載體、轉化后同質化效率的提高等提供有益的信息［11］，以促進目的基因高效穩(wěn)定表達，進而更好地了解基因［12］。本研究對木薯葉綠體基因組的密碼子組成及其偏好性進行了分析，篩選出了木薯的最優(yōu)密碼子，并根據(jù)其密碼子的偏好性，對目的基因的密碼子進行了改良，以期使盡可能多的最優(yōu)密碼子參與到外源基因的表達中，為木薯遺傳性狀的改良、系統(tǒng)發(fā)育及功能基因組學的研究提供科學依據(jù)。

1 材料和方法

1.1 試驗材料

木薯葉綠體基因組系列信息從NCBI數(shù)據(jù)庫獲取，GenBank登錄號為NC_035239，共83個基因序列。剔除重復長度小于300 bp，且起始密碼子為ATG，終止密碼子為TAA、TAG、TGA的序列，最終選出51條序列用于分析。

1.2 分析方法

木薯基因組密碼子的組成及相對同義密碼子使用度（RSCU）利用Codon W軟件參照吉德娟等［13］的方法進行分析。最優(yōu)密碼子的篩選與確定參照原曉龍等［14］的方法。中性繪圖分析參照唐玉娟等［15］的方法。ENC-plot繪圖分析參照吉德娟、李路麗等［13，16］的方法。PR2-plot繪圖分析參照李路麗等［16］的方法。

2 結果與分析

2.1 木薯基因組密碼子的組成分析

從表1可以看出：木薯葉綠體基因組密碼子基因編碼區(qū)GC的含量在28.31%～45.56%之間，平均含量為37.81%；GC在密碼子中的分布不均勻，表現(xiàn)為GC1＞GC2＞GC3，其含量平均值分別為46.64%、39.19%和27.58%；GC3s的平均含量為24.46%，表明木薯葉綠體密碼子中A和T的使用頻率高于G和C。

有效密碼子數(shù)（ENC）是衡量密碼子偏好性的重要指標之一。ENC的理論取值范圍為21～61，其值越接近20則偏好性越強；當ENC=61時，說明每個密碼子被均衡使用［17］。由表1可見，木薯葉綠體基因組的ENC均值為47.57，分布于36.55～60.29之間，故可以認為木薯葉綠體基因組密碼子的偏好性較弱。

表1 木薯葉綠體基因組密碼子的主要參數(shù)

木薯葉綠體基因組密碼子的相關分析結果（表2）顯示：GCall與GC1、GC2和GC3之間的相關性均極顯著；GC1與GC2極顯著相關；GC1、GC2均與GC12呈極顯著相關；GC1、GC2、GC12與GC3的相關性均不顯著，說明第1位與第2位堿基的組成相似，與第3位堿基的組成不同，木薯葉綠體基因組密碼子主要受選擇作用的影響；ENC值與GC1間相關性不顯著，但與GC2和GC3呈顯著或極顯著相關，表明木薯葉綠體基因密碼子的使用偏好性受第2位和第3位堿基組成的影響。

2.2 相對同義密碼子的使用度分析

基因表達水平與密碼子使用偏好性受同義密碼子第3位堿基含量的直接影響［18］。RSCU分析結果（表3）顯示，在RSCU＜1的密碼子中，以G或C結尾的占91%。RSCU＞1的密碼子共有30個，其中偏好性最強的是TTA，其RSCU值為2.05；偏好性最弱的是CTC和TCA，其RSCU值均為0.34；在這30個密碼子中，有97%的密碼子以A（13個）或T（16個）結尾，有3%的密碼子以G（1個）結尾，沒有以C結尾的密碼子。木薯葉綠體基因組密碼子第3位A或T的使用偏好性比C或G強，說明木薯偏好使用以A或T結尾的密碼子。

2.3 最優(yōu)密碼子的篩選

從木薯葉綠體基因組中篩選出atpF、atpE、cemA、rpl20、ycf2等5個高表達基因，以及rps14、petB、petD、rps8、ndhA等5個低表達基因。計算它們的△RSCU值，發(fā)現(xiàn)△RSCU＞0.08的高表達密碼子有31個。結合木薯葉綠體基因的相對同義密碼子使用度和同義密碼子相對使用度（表3和表4），篩選出最優(yōu)密碼子9個，分別是TTG、TCA、CCA、CAT、AAT、GAT、CGA、AGA、GGA。其中5個以A結尾，3個以T結尾，1個以G結尾，再次說明木薯葉綠體基因組偏好使用以A和T結尾的密碼子。

表3 木薯葉綠體基因組各氨基酸的RSCU分析結果

表4 木薯葉綠體基因組最優(yōu)密碼子的確定

2.4 密碼子使用偏好性的影響因素

2.4.1 中性繪圖分析 由圖1可見：除ycf2和cemA基因靠近對角線外，其他基因均距離對角線較遠；GC12的值在0.32～0.54之間，GC3的值在0.22～0.36之間，兩者的分布范圍均較窄；線性回歸方程的回歸系數(shù)為0.3483，接近于0，說明GC12與GC3間相關不顯著。另外，在表2中，GC12與GC3間的相關系數(shù)為0.277，顯示第1、2位與第3位堿基組成間的相關性不顯著。上述結果表明影響木薯葉綠體基因組密碼子偏性使用的主要因素是自然選擇。

圖1 中性繪圖分析結果

表2 木薯葉綠體基因組密碼子參數(shù)間的相關系數(shù)

2.4.2 ENC-plot繪圖分析 從圖2可以看出，GC3s值的分布范圍在0.22～0.36之間，分布較廣，說明不同基因密碼子的偏好性相差較大。ENC頻數(shù)分布結果（圖4）顯示：45.1%共23個基因分布在-0.05～0.05區(qū)間內，這23個基因距離標準曲線較近，其ENC的實際值與預期值接近；54.9%共28個基因分布在-0.05～0.05區(qū)間之外，這28個基因距離標準曲線較遠，其ENC的實際值與預期值偏差較大。上述結果同樣說明影響木薯葉綠體基因組密碼子偏性使用的主要因素是選擇。

圖2 ENC-plot繪圖分析結果

圖4 ENC比值頻數(shù)分布圖

2.4.3 PR2-plot繪圖分析 由圖3可以看出，木薯葉綠體基因組密碼子第3位堿基在4個區(qū)域內分布不均勻，幾乎全部點落在圖3的左上方，且堿基A的使用頻率高于T，堿基C的使用頻率高于G，因此可以推斷木薯葉綠體基因組密碼子的偏好性受多種因素的影響，但主要影響因素是選擇。

圖3 木薯PR2-plot繪圖分析結果

3 討論與結論

植物葉綠體基因組密碼子在基因突變、自然選擇和蛋白質結構等諸多因素的協(xié)同作用下，形成了適應自身偏好的模式以適應環(huán)境的變化，其中最重要的影響因素是基因突變和自然選擇。DNA序列中GC含量的變化可反映密碼子使用偏性突變的整體趨勢［19］；其中第3位堿基的突變不影響編碼氨基酸的變化，可作為密碼子偏好性分析的主要依據(jù)［20］。本研究首次對木薯葉綠體基因組51條密碼序列進行了分析，發(fā)現(xiàn)木薯葉綠體密碼子的GC含量在28.31%～45.56%之間，平均含量為37.81%；GC3s的平均含量為24.46%，表明木薯葉綠體密碼子偏好使用末位堿基為A/T的密碼子，與芒果［15］、紫花苜蓿［21］、劍麻［22］等雙子葉植物一樣偏好使用以A/T結尾的密碼子。ENC是衡量密碼子偏倚程度的重要指標之一，其值越小、越接近20，說明其偏好性越強，對應基因的表達量也越高。在本研究中，木薯葉綠體基因組的ENC均值為47.57，說明木薯葉綠體基因組密碼子的偏好性較弱。木薯葉綠體基因第1位與第2位堿基的組成相似，但與第3位堿基的組成不同，據(jù)此可以初步判斷木薯葉綠體基因組密碼子的偏好性主要受選擇作用的影響；中性繪圖、ENC-plot繪圖與PR2-plot繪圖分析結果再次證實木薯葉綠體基因組密碼子的使用偏好性主要受選擇的影響，該研究結果與其他作物如霸王［13］、美國紅梣［23］、杜梨［24］等相似，說明堿基組成成分的不同導致了密碼子使用頻率的差異。

基因的表達受同義密碼子偏好使用的影響，最優(yōu)密碼子可使翻譯的效率和準確率得到提高［25］。最優(yōu)密碼子在正選擇和突變壓力的作用下得以大量形成，而在純化選擇和突變壓力作用下則形成較少［26］。本研究通過高頻密碼子和高表達密碼子作為篩選最優(yōu)密碼子的標準，最終篩選出TTG、TCA、CCA、CAT、AAT、GAT、CGA、AGA、GGA等9個最優(yōu)密碼子；最優(yōu)密碼子數(shù)量較少，可能與木薯葉綠體基因組在進化過程中受純化選擇和突變壓力作用有關。最優(yōu)密碼子的確定，可以為木薯葉綠體改造外源基因密碼子、提高轉入表達量等提供參考。

筆者對木薯葉綠體基因組密碼子的組成與使用偏好性進行了系統(tǒng)研究，結果認為木薯葉綠體基因組密碼子的偏好性較弱，偏好使用以A/T結尾的密碼子，偏好性使用主要受選擇作用的影響；共篩選出9個最優(yōu)密碼子。本研究結果可以為木薯基因工程組學、異源表達基因改造、遺傳性狀的改良、種源鑒定等方面的研究提供科學依據(jù)。