小黑麥葉綠體基因組密碼子偏好性分析與基因功能相關(guān)性分析

摘要：為揭示小黑麥（Triticale）的密碼子使用偏好性及其影響因素，用CodonW 1.4.2和在線軟件CUSP（Codon Usage and Synonymous Condon Usage Profile）對小黑麥葉綠體基因組中52個編碼序列（Coding sequence，CDS）特征進行分析。結(jié)果表明：小黑麥葉綠體基因組密碼子第三位主要偏向以A/U為主；有效密碼子數(shù)（Effective number of codon，ENC）值在37.64～61.00之間，平均值為47.31，大于45的有42個；密碼子適應(yīng)指數(shù)（Codon adaptation index，CAI）在0.12～0.33之間，平均值為0.17，說明小黑麥葉綠體基因組密碼子偏好性較弱；16個均以A/U結(jié)尾的密碼子被篩選為最優(yōu)密碼子。中性繪圖、ENC-GC_3s和偏倚性繪圖分析（Parity Rule 2，PR2-plot）結(jié)果表明多數(shù)功能基因的密碼子使用偏好性受自然選擇壓力影響為主。小黑麥葉綠體基因組中與光合作用密切相關(guān)的基因可能受到了突變和選擇的強烈作用，使這些基因的密碼子使用頻率具有共同特征，表現(xiàn)為密碼子使用模式與基因功能有較大的相關(guān)性。

關(guān)鍵詞：小黑麥；葉綠體基因組；密碼子偏好性；最優(yōu)密碼子

中圖分類號：S544.9 """""""文獻標識碼：A """""""文章編號：1007-0435（2025）03-0707-12

Correlation Between Codon Usage Patterns and Gene Functions in the Chloroplast Genome of Triticale

TIAN Tian，"ZHANG Yin-xia，"DU Wen-hua^*

（College of Grassland Science，"Gansu Agricultural University/Key Laboratory of Grassland Ecosystem，"Gansu Agricultural University，"Ministry of Education/Pratacultural Engineering Laboratory of Gansu Province/Sino-U.S. Centers for Grazingland Ecosystem Sustainability，"Lanzhou，"Gansu Province 730070，"China）

Abstract：To investigate the codon usage bias and its influencing factors in triticale （Triticale），"CodonW 1.4.2 and the online tool CUSP （Codon Usage and Synonymous Condon Usage Profile）"were used to analyze the characteristics of 52 coding sequences （CDS）"in the chloroplast genome of triticale. The results showed that the third codon position in the triticale chloroplast genome predominantly favored A/U bases. The effective number of codons （ENC）"ranged from 37.64 to 61.00，"with an average of 47.31，"and 42 sequences had ENC values greater than 45. The codon adaptation index （CAI）"ranged 0.12-0.33，"with an average of 0.17，"indicating a weak codon usage bias in the triticale chloroplast genome. Sixteen codons ending in A/U were identified as optimal codons. Neutral plotting，"ENC-GC3s analysis，"and PR2-plot analysis （Parity Rule 2）"revealed that the codon usage bias of most functional genes was mainly influenced by natural selection. Genes in the triticale chloroplast genome closely related to photosynthesis were likely subjected to strong mutation and selection pressures，"resulting in shared codon usage patterns among these genes. This suggests a significant correlation between codon usage patterns and gene functions.

Key words：Triticale；Chloroplast genome；"Codon usage bias；"Optimal codon

葉綠體是一種半自主細胞器，通過光合作用將太陽能轉(zhuǎn)化為生物能，參與了許多其他化合物的合成^［1-4］，葉綠體基因組是其實現(xiàn)活躍的代謝功能的必要條件^［3-4］。單個葉綠體基因組的詳細表征與比較葉綠體基因組學(xué)可以作為功能基因組研究的基礎(chǔ)^［5］，也可以指導(dǎo)葉綠體基因工程的試驗^［6-7］。葉綠體基因轉(zhuǎn)化具有外源基因表達效率高、穩(wěn)定遺傳、定點整合以及沒有位置效應(yīng)現(xiàn)象等優(yōu)點^［8-9］，可以用來開發(fā)分子標記、系統(tǒng)發(fā)育分析、種群水平研究、基因分型和定位等研究^［10］。密碼子是生物體遺傳信息和蛋白質(zhì)翻譯之間轉(zhuǎn)化的橋梁，在實現(xiàn)遺傳信息的傳遞和表達上起到了至關(guān)重要的作用。蛋白質(zhì)翻譯過程中，翻譯為同種氨基酸的不同密碼子稱為同義密碼子，同義密碼子在組合頻率上也有所差異，這種現(xiàn)象稱為密碼子偏好性（Codon usage bias，"CUB）^［11］。同義密碼子偏好性可能是由進化過程中基因上的不同因素造成的^［12］，是基因表達的一種微調(diào)方式^［13］。分析目標基因或受體植物的基因密碼子使用模式，可以為構(gòu)建植物成熟穩(wěn)定的葉綠體轉(zhuǎn)基因系統(tǒng)奠定基礎(chǔ)^［14-15］。研究表明，密碼子偏好性由自然選擇和突變共同決定，前人提出了“選擇-突變-漂移”假說和“中性理論”假說^［16］，前者認為密碼子的使用偏倚性是由最優(yōu)密碼子的選擇壓力和非同義密碼子的突變-漂移所決定，后者認為同義密碼子的突變是中性選擇，與選擇壓力無關(guān)^［17-18］。除此之外，GC含量的高低、上下游依賴性核苷酸的組成、種群大小以及基因表達水平和長度等因素也會對密碼子的偏倚性產(chǎn)生影響^［11］。同時，密碼子的使用偏好會影響mRNA折疊、翻譯延伸和蛋白折疊，從而對基因的表達水平產(chǎn)生不同程度的影響^{［13，19］}。

小黑麥（Triticale）是由小麥屬（Triticum）和黑麥屬（Secale）經(jīng)過屬間雜交和染色體數(shù)加倍而人工育成的新品種^［20］。小黑麥雜種優(yōu)勢強，具有產(chǎn)量高、品質(zhì)好以及抗逆性強等優(yōu)勢，是繼苜蓿、飼用玉米后的在西北地區(qū)的優(yōu)質(zhì)飼草作物^［21］。同時，小黑麥越冬性能好，在高寒農(nóng)牧交錯區(qū)能夠錯季種植，保證了應(yīng)季作物的正常播種，同時減少了冬春季節(jié)農(nóng)田裸露帶來的環(huán)境問題^［22-25］。小黑麥六倍體葉綠體基因組已完成測序（GenBank登錄號："ON422331.1），但其密碼子偏好性研究尚未有見報道。明確小黑麥葉綠體基因的密碼子偏好性是開展小黑麥葉綠體基因工程的基礎(chǔ)之一，此工程在改良植物性狀和提高抵抗病蟲害等方面具有重要作用^［26］。培育高產(chǎn)、抗病性強、抗逆性強的小黑麥一直是遺傳改良的目標，但受到小黑麥遺傳背景復(fù)雜、倍性多樣和傳統(tǒng)育種周期長等問題的限制。分子育種一定程度上可以緩解這些問題，但外源基因的表達水平不穩(wěn)定，會受到密碼子偏好性的影響^［27］。另外，規(guī)律間隔成簇短回文重復(fù)序列（Clustered regularly interspaced short palindromic repeats，CRISPR）基因編輯技術(shù)中，Cas9蛋白的表達量在基因編輯中起著關(guān)鍵作用。Cas9主要通過sgRNA的引導(dǎo)定位到目標DNA序列進行特異性切割DNA雙鏈從而在DNA斷裂時引入突變實現(xiàn)基因的敲除、插入或替換，因此，Cas9的定位效率直接影響到編輯效率，而Cas9的定位效率與其表達量有關(guān)^［28-30］。Ma等人^［31］研究發(fā)現(xiàn)，影響Cas9定位效率的因素主要有以下幾點：1）不同啟動子驅(qū)動的sg RNA靶向的位點突變率相似；2）高GC含量（高于50%）具有更高的編輯效率；3）所有堿基都具有相似的編輯效率。這一發(fā)現(xiàn)擴大了有效基因組編輯目標的選擇范圍，不限于U3啟動子的A和U6啟動子的G的常規(guī)靶標^［28］。研究表明，密碼子的使用模式會影響基因的表達^［13］，Shan等^［32］對密碼子進行優(yōu)化了Cas9，靶向水稻OsPDS基因的突變率達14.5%～20%，Zhang等^［33］使用了兩種密碼子類型的Cas9（h Sp Cas9，p Sp Cas9）研究靶向水稻除草劑抗性基因BEL基因，結(jié)果表明，p Sp Cas9比h Sp Cas9的編輯效率高了7.42倍，差異可能與p Sp Cas9轉(zhuǎn)基因植物中較高水平的Cas9蛋白積累相關(guān)。因此，Cas9的表達量影響其定位效率，從而影響基因編輯效率。而當前尚未有任何報道關(guān)于Cas9在小黑麥葉綠體基因組密碼子偏好性中的關(guān)系。因此，通過研究小黑麥葉綠體基因組密碼子的使用特征，確定最優(yōu)密碼子，并探討植物葉綠體基因組的密碼子使用模式與基因功能的相關(guān)性，該研究結(jié)果能夠為提高基因在小黑麥中表達、預(yù)測及外源受體的選擇提供一定參考，同時，在小黑麥基因編輯的研究中，能夠為提高小黑麥Cas9表達從而提高編輯效率提供一定的研究價值。

1 材料與方法

1.1 序列獲取

從NCBI數(shù)據(jù)庫下載小黑麥的葉綠體基因組序列（GenBank登錄號："ON422331.1）87條，剔除長度小于300 bp的序列、去除重復(fù)序列、去除非ATG作為起始密碼子的序列、刪掉異常終止密碼子和CDS內(nèi)部存在終止密碼子的序列后，得到52條序列用于后續(xù)分析。

1.2 密碼子組成分析

用CodonW 1.4.2對小黑麥葉綠體基因組的52條CDS序列進行同義密碼子相對使用度（Relative usage of synonymous codons，RSCU）、密碼子適應(yīng)指數(shù)（Codon adaptation index，CAI）、密碼子偏好指數(shù)（Codon bias index，CBI）和最優(yōu)密碼子使用頻率值（Frequency of optimal codons，F(xiàn)OP）分析；使用在線軟件CUSP（http："http://emboss.toulouse.inra.fr/cgi-bin/emboss/cusp）統(tǒng)計得到的52條序列的密碼子數(shù)量（Codon Number，CN）、密碼子不同位置堿基的GC含量，包括第1位GC含量（GC₁）、第2位GC含量（GC₂）、第3位GC含量（GC₃）以及GC總含量（GC_all），同時計算出GC₁和GC₂的平均值（GC₁₂）利用R語言進行相關(guān)性分析。

1.3 中性繪圖

以GC₃為橫坐標，GC₁和GC₂的平均值GC₁₂為縱坐標繪制中性圖。若密碼子偏好性僅受到突變的影響，則基因分布于對角線上；若受到選擇壓力的影響更大，則基因分布于對角線周圍^［34］。同時，可通過橫縱坐標GC₃和GC₁₂的相關(guān)性來判斷密碼子偏好性受到的影響因素，若兩者顯著相關(guān)，回歸系數(shù)（對角線斜率）較大且接近于1，則密碼子受到的影響因素主要為突變，反之（回歸系數(shù)接近于0）則為選擇壓力^［35］。

1.4 ENC-plot分析

有效密碼子數(shù)（Effective Number of Codon，ENC）用于度量同義密碼子的偏倚性，當ENC值=20時，表示每個氨基酸只使用一個密碼子，為極端偏好；ENC值=61，表示密碼子偏向隨機使用，不存在使用偏好。以公式（1）為標準曲線^［36］，即期望值，每點代表一個基因。

公式：

（1）

每個基因與標準曲線（即期望值）之間的距離為判斷影響密碼子偏倚性因素的標準。當ENC實際值與期望值相差較大，即距離標準曲線較遠時，偏倚性受選擇壓力的影響，反之（即ENC實際值距離標準曲線較近時）則主要由堿基組成決定，即受突變的影響較大。為了更好地反映ENC實際值與期望值之間的差異，可以通過計算ENC比值頻率分布，以-0.05～0.05的范圍為界對差異進行量化分析。在此范圍內(nèi)則接近期望值，即密碼子偏好性受到突變的影響較大；在此范圍外則影響密碼子偏好性的主要因素為選擇壓力。

1.5 PR2-plot繪圖分析

以A₃/（A₃+T₃）為縱坐標，以G₃/（G₃+C₃）為橫坐標，繪制PR2偏倚平面圖，每個基因的堿基組成在一個平面上表示，中心點表示密碼子無偏好性使用時的狀態(tài)（A=T，C=G），而其他各點與中心點的矢量距離表示偏倚的方向和程度^［37］，以此來判斷突變壓力和自然選擇是否影響了DNA雙鏈上核苷酸的組成。如果4種堿基的使用頻率相似，則基因坐落于中線或靠近中心點，表示密碼子偏好性只受到突變的影響；反之坐落于中心點四周或遠離中心點，表示同時受到了突變和自然選擇的影響。

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

為了得到小黑麥葉綠體基因的最優(yōu)密碼子，以ENC值為參考指標，篩選前后各10%的基因構(gòu)建高低基因表達庫，通過Codon W1.4.2計算其RSCU值，并計算高低表達庫之間的差值（△RSCU），將△RSCU≥0.08且RSCUgt;1的密碼子篩選出，確定為最優(yōu)密碼子^［38］。

2 結(jié)果與分析

2.1 密碼子組成分析

小黑麥葉綠體基因組總長為135 899 bp，從中得到87個CDS（不含假基因），經(jīng)篩選而得到52個CDS序列用于后續(xù)分析。結(jié)果如表1所示，52個CDS的總密碼子數(shù)為17 152。總GC含量分布范圍為32.55%～49.20%，平均值為38.91%。GC₁，GC₂和GC₃含量分布范圍分別在32.75%～56.28%，29.87%～55.56%，22.22%～48.28%之間，平均值為GC₁（46.83%）gt;GC₂（39.62%）gt;GC₃（30.28%），表明密碼子分布的不均勻性，GC含量較低，主要以A/U結(jié)尾。其中，GC₂和GC_all的差異不大（39.62%和38.91%），而GC₃最?。?0.28%），所有GC含量均在50%以下，說明小黑麥葉綠體基因組密碼子向A/U偏倚。

根據(jù)偏倚性分析結(jié)果（表1），ENC在37.64～61.00之間（平均值為47.31），其中有42個基因的ENC值大于45，表明小黑麥葉綠體基因組密碼子偏好性相對較弱。CAI范圍為0.12～0.33，平均值為0.17，趨近于0，表明小黑麥葉綠體基因組密碼子偏好性弱。CBI取值為-0.21～0.24，均值為-0.09，則最優(yōu)密碼子使用次數(shù)小于密碼子平均使用次數(shù)；FOP取值為0.27～0.56，均值為0.35，則最優(yōu)密碼子占同義密碼子數(shù)量的比值較小。以上結(jié)果表明小黑麥葉綠體基因組密碼子使用偏好性較弱。

各參數(shù)相關(guān)性分析表明（圖1），GC_all與GC₁，GC₂以及GC₃呈極顯著相關(guān)關(guān)系（Plt;0.01），GC₁和GC₂呈顯著相關(guān)關(guān)系（Plt;0.05），但GC₃和GC₁，GC₂相關(guān)性差，表明小黑麥葉綠體基因組密碼子第1，2位堿基組成上較相似，而與第3位存在較大差異。ENC與GC₃呈極顯著相關(guān)關(guān)系，而與GC₂呈極顯著負相關(guān)關(guān)系，與GC₁沒有相關(guān)性，表明密碼子第2位和第3位對ENC具有偏性影響。密碼子數(shù)CN與ENC，GC₁，GC₂，GC₃以及GC_all均無相關(guān)性，說明小黑麥密碼子堿基受偏性影響較小。

密碼子總體ESCU分析表明（圖2），Leu，Ser和"Arg分別由6個密碼子編碼，Val，Pro，Thr，Ala和Gly分別由四個密碼子編碼，其余除了Met和Trp兩種氨基酸分別由一個密碼子編碼以及Lie由3個密碼子編碼之外，剩余所有的氨基酸均由兩個密碼子編碼。

由圖3可知，編碼Ile的密碼子AUU數(shù)量最多，有711個；編碼Ter的密碼子UGA最少，僅有12個。小黑麥葉綠體基因組中密碼子RSCUgt;1的有33個，其中29個以A/U結(jié)尾，RSCU＜1的有31個，具有明顯的A/U偏向性。

2.2 中性繪圖分析

通過中性繪圖來判斷小黑麥葉綠體基因組密碼子偏好性的主要影響因素，結(jié)果如圖4。GC₁₂的取值范圍在0.3577～0.5348之間，GC₃取值范圍在0.2222～0.4828之間，僅有1個基因在對角線位置下方，其余所有基因均位于對角線位置上方，說明主要的影響來源是自然選擇；GC₁₂與GC₃回歸系數(shù)為0.119，相關(guān)系數(shù)R為0.0671，相關(guān)性弱，另外，GC含量僅為47.43，進一步表明自然選擇是主要的影響因素。同時，涉及保守假設(shè)性葉綠體ORF、核糖體大亞基蛋白和翻譯起始因子3種功能中的4個基因位于對角線附近甚至下方，說明在一定程度上，個別基因受到突變的影響較大。

2.3 ENC-plot分析

ENC-plot繪圖（圖5）分析發(fā)現(xiàn)，實際ENC值與期望值存在較大差異，除了涉及NADH脫氫酶亞基、核糖體大亞基蛋白和核糖體小亞基蛋白3種功能中的4個基因ENC值接近期望值（標準曲線）以外，其他基因的ENC值均高于或低于期望值，即大多數(shù)與標準曲線的距離較遠，說明自然選擇是密碼子偏好性的主要影響因素。后續(xù)通過統(tǒng)計小黑麥葉綠體基因ENC比值頻率以更為準確地反映差異，結(jié)果表明（圖6），僅有11個基因分布于-0.05～0.05之間，占比21.15%，其余41個基因（占比78.85%）全部在此范圍外，距離ENC標準期望值較遠，進一步說明小黑麥葉綠體基因組密碼子偏好性的主要影響因素是選擇壓力。

2.4 PR2-plot繪圖分析

通過PR2-plot分析密碼子偏倚性，當突變?yōu)槲ㄒ挥绊懸蛩貢r，四種堿基在圖中的表現(xiàn)為均衡，即集中于中心點，而當有選擇壓力同時存在時，根據(jù)影響大小而不同程度偏離中心點，向四周不均勻排布。小黑麥葉綠體基因組PR2-plot分析結(jié)果表明（圖7），涉及保守假設(shè)性葉綠體ORF和RNA聚合酶亞基兩種功能中的2個基因接近于中心點，偏向程度小，表明其密碼子偏好性受到突變的影響較多。其他大多數(shù)基因均坐落于平面圖的下方，尤其以右下方基因數(shù)量最多，說明在堿基出現(xiàn)的頻率上G₃gt;C₃，T₃gt;A₃；RuBisCO酶大亞基全部基因落于左下方，即密碼子偏倚性為C₃gt;G₃，T₃gt;A₃。圖中可看出各基因分布不均勻，以平面圖左下和右下居多，說明4種堿基的使用頻率存在一定的差異，且同時受到突變和選擇的多重因素的影響。

整體而言，多數(shù)基因在A₃/U₃偏倚表現(xiàn)出了功能相關(guān)性：光合系統(tǒng)I亞基、光合系統(tǒng)II亞基、細胞色素c基因、NADH脫氫酶亞基、成熟酶等基因的密碼子都更偏好使用U₃；RNA聚合酶亞基的眾多基因中，僅有一個未表現(xiàn)出明顯的A₃/U₃偏倚，其余基因都更偏好使用U₃；核糖體小亞基蛋白的眾多基因中，多數(shù)都更偏好使用A₃。然而，在C₃/G₃偏倚方面，整體上未表現(xiàn)出明顯的基因功能相關(guān)性。小黑麥葉綠體基因組密碼子偏好性受到了突變和選擇的共同影響，而這種影響可能與相應(yīng)基因的功能是相關(guān)的。

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

通過構(gòu)建高、低表達基因庫，將△RSCU≥0.08且RSCU≥1的密碼子選出作為最優(yōu)密碼子，共得到GCU，AGA和CGU等16個最優(yōu)密碼子（GCU，AGA，CGU，CAA，GGU，CAU，AUU，UUA，AAA，UUU，CCU，AGU，UCU，ACA，ACU，GUU），全部以A/U結(jié)尾（表2），這與GC₃、RSCU分析結(jié)果一致。

3 討論

研究在獲得了小黑麥的完整葉綠體基因組序列之后，探究了52個CDS的密碼子使用特征及其與基因功能的關(guān)系。首先，將小黑麥葉綠體基因組中的蛋白編碼序列作為一個整體進行分析而不考慮其功能基因時，密碼子組成特征分析表明，GC含量較低，且GC₁（46.83%）gt;GC₂（39.62%）gt;GC₃（30.28%），說明密碼子分布不均，偏向于A/U結(jié)尾，這與前人在高粱（Sorghum bicolor"L.）、葦狀羊茅（Festuca arundinacea"Schreb）、大麥（Hordeum vulgare"subsp. vulgare）以及多年生黑麥草（Lolium perenne"L.）等23種禾本科（Gramineae）植物中的研究結(jié)果一致^［39-40］。根據(jù)RSCU分析結(jié)果，小黑麥葉綠體基因組的編碼氨基酸的密碼子中，偏好使用和非偏好使用的密碼子種類幾乎持平，但總體偏好使用的密碼子多以A/U結(jié)尾，結(jié)果與此前研究所報道的雙子葉植物葉綠體密碼子使用偏好特征一致^{［14，41-42］}，符合“高等植物密碼子傾向于使用A/U結(jié)尾”的假設(shè)^［43］，說明植物葉綠體密碼子的總體使用模式在進化過程中比較保守。密碼子使用特征的數(shù)據(jù)主要應(yīng)用于優(yōu)化或重新設(shè)計基因，以便在異源系統(tǒng)中獲得高表達或最佳表達，本研究確定了小黑麥葉綠體基因組的16個最優(yōu)密碼子，在進行基因工程研究時，針對最優(yōu)密碼子的優(yōu)化將有利于調(diào)控相關(guān)基因的表達^［38］。

將密碼子與基因功能相聯(lián)合分析時，不同功能的基因使用的密碼子具有較大的區(qū)別。研究表明，CAI值越高，基因的密碼子使用偏好性越強，在生物體內(nèi)的表達水平越高^［44］。小黑麥葉綠體基因組的52個CDS中，各基因的CAI值在0～1.0之間，平均值為0.17，密碼子適應(yīng)指數(shù)CAI總體偏小，說明該基因組的密碼子使用偏好性總體偏弱。相較而言，小黑麥的rpl16，rps8，ycf68，ndhG和rpl20等基因的密碼子使用偏好性較弱，而psbB，rps3，psbD，rbcL和psbA等基因的密碼子使用偏好性較強，說明psbB，rps3，psbD等基因可能在小黑麥的葉綠體中表達水平較高。事實上，光合作用功能是葉綠體的最重要的功能^［3］，這意味著與光合作用相關(guān)的psbB，rps3，psbD等基因必然有較高的表達水平。

對于各基因，ENC的值介于20（當每個氨基酸只使用一個密碼子時存在極端偏差）和61（當所有密碼子都被隨機使用時）之間，小黑麥葉綠體CDS的ENC保持在較高水平（均值為47.31gt;45.00），說明密碼子使用的偏好性總體較弱^［45］，與CAI的分析結(jié)果一致。值得注意的是，ENC值高達61的inFA屬于翻譯起始因子，通過與核糖體、mRNA和tRNA的相互作用，幫助識別起始密碼子并啟動蛋白質(zhì)合成，但在小黑麥的葉綠體基因組上所有密碼子都被隨機使用，說明該基因表達量低，起到特定識別的能力較弱。中性繪圖、ENC-GC_3s和PR2-plot均顯示絕大多數(shù)功能基因的密碼子使用的偏好性受突變和自然選擇壓力的共同影響，且以自然選擇的影響為主。生物現(xiàn)有同義密碼子使用模式的形成是突變、自然選擇和基因翻譯效率的遺傳漂變之間平衡的結(jié)果^［46-47］，小黑麥葉綠體基因組中與光合作用密切相關(guān)的重要功能基因必然會受到突變和選擇的強烈作用，也許正是這些原因?qū)е铝诵『邴溔~綠體基因組上與光合作用相關(guān)基因的密碼子使用常表現(xiàn)出更鮮明的共同特征。

本研究結(jié)果可以為小黑麥葉綠體基因組在基因工程中的應(yīng)用提供依據(jù)。密碼子使用模式及其影響因素的分析結(jié)果表明，小黑麥葉綠體基因，尤其是與光合作用密切相關(guān)的基因，可能是因受到選擇突變和選擇的強烈影響，其密碼子使用模式與基因功能的有較大的相關(guān)性。將來的研究可以通過結(jié)合生物信息學(xué)和系統(tǒng)生物學(xué)等技術(shù)手段，進行大規(guī)模的葉綠體基因組密碼子偏好性研究，以探究葉綠體基因組密碼子偏好性與基因表達調(diào)控之間的關(guān)系，從而揭示更多的密碼子使用模式規(guī)律和生物學(xué)意義，以便更好地了解葉綠體基因表達的分子機制和密碼子使用模式在生物演化中的作用和意義。

4 結(jié)論

通過對小黑麥葉綠體基因組密碼子偏好性特征及其影響因素進行多角度分析，篩選出16個最優(yōu)密碼子，明確了密碼子偏好性較弱，但對A/U具有明顯的偏倚性，造成偏倚性的主要因素為選擇壓力。研究表明小黑麥葉綠體基因組密碼子使用模式與功能基因之間具有較大的相關(guān)性，尤其是光合作用相關(guān)基因的密碼子使用模式上，更偏向于使用A₃/U₃。研究結(jié)果可為小黑麥基因工程提供理論基礎(chǔ)，特別是在提高基因表達量的設(shè)計上或能提供研究基礎(chǔ)。

參考文獻

［1］"VIDI P A，KANWISCHER M，BAGINSKY S，et al. Tocopherol cyclase （VTE1）"localization and vitamin E accumulation in chloroplast plastoglobule lipoprotein particles［J］. Journal of Biological Chemistry，2006，281（16）：11225-11234

［2］"LIPPOLD F，VOM DORP K，ABRAHAM M，et al. Fatty acid phytyl ester synthesis in chloroplasts of Arabidopsis［J］. The Plant Cell，2012，24（5）：2001-2014

［3］"DANIELL H，JIN S X，ZHU X G，et al. Green giant—a tiny chloroplast genome with mighty power to produce high-value proteins：history and phylogeny［J］. Plant Biotechnology Journal，2021，19（3）：430-447

［4］"DOBROGOJSKI J，ADAMIEC M，LUCI?SKI R. The chloroplast genome：a review［J］. Acta Physiologiae Plantarum，2020，42（6）：98

［5］"SUGIURA M. History of chloroplast genomics［J］. Photosynthesis Research，2003，76（1/2/3）：371-377

［6］"GREVICH J J，DANIELL H. Chloroplast genetic engineering：recent advances and future perspectives［J］. Critical Reviews in Plant Sciences，2005，24（2）：83-107

［7］"DANIELL H，LIN C S，YU M，et al. Chloroplast genomes：diversity，evolution，and applications in genetic engineering［J］. Genome Biology，2016，17（1）：134

［8］"KWAK S Y，LEW T T S，SWEENEY C J，et al. Chloroplast-selective gene delivery and expression in planta using chitosan-complexed single-walled carbon nanotube carriers［J］. Nature Nanotechnology，2019，14（5）：447-455

［9］"RUF S，HERMANN M，BERGER I J，et al. Stable genetic transformation of tomato plastids and expression of a foreign protein in fruit［J］. Nature Biotechnology，2001，19（9）：870-875

［10］"RAUBESON L A，PEERY R，CHUMLEY T W，et al. Comparative chloroplast genomics：analyses including new sequences from the angiosperms Nuphar advena and Ranunculus macranthus［J］. BMC Genomics，2007，8：174

［11］"高守輿，李鈺瑩，楊志青，等. 白羊草葉綠體基因組密碼子使用偏好性分析［J］. 草業(yè)學(xué)報，2023，32（7）：85-95

［12］"IKEMURA T. Codon"usage and tRNA content in unicellular and multicellular organisms［J］. Molecular Biology and Evolution，1985，2（1）：13-34

［13］"QUAX T E F，CLAASSENS N J，S?LL D，et al. Codon"bias as a means to fine-tune gene expression［J］. Molecular Cell，2015，59（2）：149-161

［14］"WANG Z J，XU B B，LI B，et al. Comparative analysis of Codon"usage patterns in chloroplast genomes of six Euphorbiaceae"species［J］. PeerJ，2020，8：e8251

［15］"SHENG J J，SHE X，LIU X Y，et al. Comparative analysis of Codon"usage patterns in chloroplast genomes of five Miscanthus species and related species［J］. PeerJ，2021，9：e12173

［16］"MAZUMDAR P，BINTI OTHMAN R，MEBUS K，et al. Codon"usage and Codon"pair patterns in non-grass monocot genomes［J］. Annals of Botany，2017，120（6）：893-909

［17］"CHEN S L，LEE W，HOTTES A K，et al. Codon"usage between genomes is constrained by genome-wide mutational processes［J］. Proceedings of the National Academy of Sciences of the United States of America，2004，101（10）：3480-3485

［18］"劉慧，王夢醒，岳文杰，等. 糜子葉綠體基因組密碼子使用偏性的分析［J］. 植物科學(xué)學(xué)報，2017，35（3）：362-371

［19］"原曉龍，李云琴，張勁峰，等. 降香黃檀葉綠體基因組密碼子偏好性分析［J］. 廣西植物，2021，41（4）：622-630

［20］"劉晶，趙明月，郭宏玉，等. 高寒牧區(qū)小黑麥研究進展［J］. 青海科技，2021，28（4）：23-26

［21］"王偉強，劉晶，田新會，等. ‘甘農(nóng)4號’小黑麥品種在青海省不同區(qū)域的適應(yīng)性評價［J］. 草地學(xué)報，2020，28（6）：1626-1634

［22］"裴亞斌，杜文華，劉漢成，等. 甘南高寒牧區(qū)3種飼草不同種植模式下的生產(chǎn)性能及經(jīng)濟效益［J］. 草業(yè)科學(xué)，2020，37（4）：791-799

［23］"楊曉，宋謙，余小亮，等. 隴東旱塬區(qū)秋播小黑麥與青貯玉米的復(fù)種效果［J］.草業(yè)科學(xué)，2019，36（8）：2127-2134

［24］"楊永強，尹玉林，杜文華. 青藏高原川西北牧區(qū)一年兩季飼草種植模式對土壤理化性質(zhì)及酶活性的影響［J］. 草原與草坪，2024，44（2）：57-66

［25］"宋謙，馬淑梅，張?zhí)m，等. 隴東地區(qū)秋播小黑麥茬后復(fù)種青貯玉米品種的篩選［J］. 中國飼料，2024（7）：140-145

［26］"倪典墨，李冉郡，郭忠靜，等. 百合科植物psbA基因密碼子偏好性分析［J］. 生物學(xué)雜志，2020，37（6）：68-72

［27］"BIN PEK H，KLEMENT M，ANG K S，et al. Exploring Codon"context bias for synthetic gene design of a thermostable invertase in Escherichia coli［J］. Enzyme and Microbial Technology，2015，75：57-63

［28］"時歡，林玉玲，賴鐘雄，等. CRISPR/Cas9介導(dǎo)的植物基因編輯技術(shù)研究進展［J］. 應(yīng)用與環(huán)境生物學(xué)報，2018，24（3）：640-650

［29］"DOUDNA J A，CHARPENTIER E. The new frontier of genome engineering with CRISPR-Cas9［J］. Science，"2014，346（6213）：1258096

［30］"JIANG F G，DOUDNA J A. CRISPR–Cas9 structures and mechanisms［J］. Annual Review of Biophysics，"2017，46（1）：505-529

［31］"MA X L，ZHANG Q Y，ZHU Q L，et al. A robust CRISPR/Cas9 system for convenient，"high-efficiency multiplex genome editing in monocot and dicot plants［J］. Molecular Plant，2015，8（8）：1274-1284

［32］"SHAN Q W，WANG Y P，LI J，et al. Targeted genome modification of crop plants using a CRISPR-Cas system［J］. Nature Biotechnology，2013，31（8）："686-688

［33］"ZHANG H，ZHANG J S，WEI P L，et al. The CRISPR/Cas9 system produces specific and homozygous targeted gene editing in rice in one generation［J］. Plant Biotechnology Journal，2014，12（6）：797-807

［34］"SUEOKA N. Directional mutation pressure and neutral molecular evolution［J］. Proceedings of the National Academy of Sciences of the United States of America，1988，85（8）：2653-2657

［35］"KAWABE A，MIYASHITA N T. Patterns of Codon"usage bias in three dicot and four monocot plant species［J］. Genes amp; Genetic Systems，2003，78（5）：343-352

［36］"WRIGHT F. The ‘effective number of codons’"used in a gene［J］. Gene，1990，87（1）：23-29

［37］"SUEOKA N. Translation-coupled violation of Parity Rule 2 in human genes is not the cause of heterogeneity of the DNA G+C content of third Codon"position［J］. Gene，1999，238（1）：53-58

［38］"WANG L J，ROOSSINCK M J. Comparative analysis of expressed sequences reveals a conserved pattern of optimal Codon"usage in plants［J］. Plant Molecular Biology，2006，61（4/5）：699-710

［39］"洪森榮，林順來，李盈萍，等. 甜高粱葉綠體基因組特征及密碼子偏好性分析［J］.草地學(xué)報，2023，31（12）：3636-3650

［40］"ZHANG Y R，NIE X J，JIA X O，et al. Analysis of Codon"usage patterns of the chloroplast genomes in the Poaceae family［J］. Australian Journal of Botany，2012，60（5）：461

［41］"季凱凱，宋希強，陳春國，等. 木蘭科葉綠體基因組的密碼子使用特征分析［J］. 中國農(nóng)業(yè)科技導(dǎo)報，2020，22（11）：52-62

［42］"YANG S W，LI G，LI H. Molecular characterizations of genes in chloroplast genomes of the genus Arachis"L.（Fabaceae）"based on the Codon"usage divergence［J］. PLoS One，2023，18（3）：e0281843

［43］"CAMPBELL W H，GOWRI G. Codon"usage in higher plants，green algae，and cyanobacteria［J］. Plant Physiology，1990，92（1）：1-11

［44］"SHARP P M，LI W H. An evolutionary perspective on synonymous Codon"usage in unicellular organisms［J］. Journal of Molecular Evolution，1986，24（1/2）：28-38

［45］"SALIM H M W，CAVALCANTI A R O. Factors influencing Codon"usage bias in genomes［J］. Journal of the Brazilian Chemical Society，2008，19（2）：257-262

［46］"LIU Q P. Mutational bias and translational selection shaping the Codon"usage pattern of tissue-specific genes in rice［J］. PLoS One，2012，7（10）：e48295

［47］"PLOTKIN J B，KUDLA G. Synonymous but not the same：the causes and consequences of Codon"bias［J］. Nature Reviews Genetics，2011，12（1）：32-42

（責任編輯""閔芝智）

引用格式：田甜， 張銀霞， 杜文華.小黑麥葉綠體基因組密碼子偏好性分析與基因功能相關(guān)性分析[J].草地學(xué)報，2025，33（3）：707-718

Citation：TIAN Tian， ZHANG Yin-xia， DU Wen-hua.Correlation Between Codon Usage Patterns and Gene Functions in the Chloroplast Genome of"Triticale[J].Acta Agrestia Sinica，2025，33（3）：707-718

基金項目：農(nóng)業(yè)生物育種重大專項（2023ZD0402605-02）；國家自然科學(xué)基金（32260339）；甘肅省高等學(xué)校產(chǎn)業(yè)支撐計劃（2022CYZC-49）；甘肅省重大專項（21ZD4NA012）；西藏重大專項（XZ202101ZD003N）資助

作者簡介：田甜（1997-），女，苗族，貴州銅仁人，博士研究生，主要從事草地生物多樣性研究，E-mail：1138115730@qq.com；

* 通信作者"Author for correspondence，E-mail：duwh@gasu.edu.cn