刺柏屬4種藥用植物葉綠體基因組密碼子偏好性分析

舒軍霞，楊 林，周 濤，沈蓮文，王大瑋*

刺柏屬4種藥用植物葉綠體基因組密碼子偏好性分析

舒軍霞1，楊 林1，周 濤1，沈蓮文2，王大瑋1*

1. 西南林業大學云南省高校林木遺傳改良與繁育重點實驗室，云南 昆明 650224 2. 西南林業大學西南山地森林資源保育與利用教育部重點實驗室，云南 昆明 650224

明確刺柏屬圓柏、垂枝香柏、昆明柏和鋪地柏4種藥用植物葉綠體基因組密碼子的使用偏好性及影響其密碼子偏好性的因素。利用CodonW、CUSP、SPSS等軟件對4種刺柏屬藥用植物葉綠體基因組密碼子偏好性進行分析。4種刺柏屬藥用植物葉綠體基因組密碼子的GC含量在35.8%～37.3%，有效密碼子數（effective number of codon，ENC）值在47.10～47.79，表明其密碼子偏好性較弱。中性繪圖、ENC-plot、PR2-plot分析表明影響4種刺柏屬藥用植物葉綠體基因組密碼子使用偏性的因素有選擇和突變。同時利用相對同義密碼子使用度（relative synonymous codon usage，）值和ENC值篩選出共65個最優密碼子，其中有9個密碼子為4種刺柏屬植物共有的最優密碼子。4種刺柏屬藥用植物密碼子偏好使用A/T結尾，其密碼子使用偏性主要受到自然選擇的影響。通過對4種刺柏屬植物葉綠體基因組密碼子的使用偏性進行分析及外類群聚類分析驗證，揭示影響其密碼子使用偏性的主要因素，為后續刺柏屬葉綠體基因組外源蛋白表達載體的構建及優化提供理論依據。

刺柏屬；圓柏；垂枝香柏；昆明柏；鋪地柏；密碼子偏好性；最優密碼子；自然選擇

葉綠體基因組是高等植物細胞質基因組的組成成分之一，攜帶大量自我調節的編碼基因。與核基因組相比具有序列長度適中、基因直系同源、進化速率適中等優勢[1-2]。這些優勢可以用于揭示物種起源、進化、演變及比較基因組學等研究方向[3-5]。隨著高通量測序技術的發展，越來越多的植物葉綠體基因組得到了測序和組裝。因此，葉綠體基因組已經廣泛應用于植物系統發育研究、物種鑒定、葉綠體基因工程和遺傳結構分析等方面的研究，為植株的多抗性轉化、農作物代謝途徑改造等方面奠定理論依據[6-7]。

密碼子是生物體內mRNA（或DNA）上的3聯體核苷酸殘基序列，是遺傳信息的序列單位，參與mRNA翻譯成蛋白質的過程[8]。密碼子的使用偏好性是指某一物種或某一基因在蛋白質翻譯過程中傾向于使用一種或幾種特定的同義密碼子的現象[9]。目前對密碼子偏好性的研究主要集中在分子進化、翻譯調控等方面[10-11]。研究表明，密碼子偏好性可以通過改變蛋白質結構、核苷酸序列等方式影響外源基因的表達[12-13]。因此，對密碼子偏好性的分析可提供密碼子偏好性優化策略，進而為利用基因工程技術改良植物重要性等研究奠定理論基礎。

圓柏L.、垂枝香柏W. C. Cheng ex Ferré、昆明柏W. C. Cheng和鋪地柏(Endlicher) Siebold ex Miquel均屬于柏科刺柏屬植物，主要分布于熱帶及亞熱帶地區。刺柏屬植物是具有藥用、景觀、生態價值的經濟樹種。刺柏屬植物的枝葉中含有具抗炎、抑菌效果的總酚、總黃酮等化合物，是中醫常用來止血、散腫毒的傳統藥物[14-15]；同時研究發現圓柏葉提取物制備的生物功能銀納米顆粒對肺癌細胞有抗轉移和抗增殖的作用[16]。此外，刺柏屬植物還具有耐寒耐旱及改善土壤環境等特點，是我國西北干旱以及被重金屬污染地區的主要栽培樹種[17-18]。刺柏屬植物在多個領域都有開發利用價值，因此加強刺柏屬葉綠體基因組的基礎研究對其高效利用具有重要意義。本研究分析了4種刺柏屬葉綠體基因組密碼子使用模式并通過外類群聚類分析驗證，確定了影響其密碼子使用偏好性的因素以及最優密碼子，為4種刺柏屬植物葉綠體基因組后續進行異源基因表達載體設計和基因工程應用和研究提供理論依據。

1 材料

所用的植物樣本采于云南省昆明市（102°10'E，24°23'N，平均海拔1891 m）西南林業大學（圓柏、昆明柏、鋪地柏，2020年7月3日）和中國科學院昆明植物研究所（垂枝香柏，2020年7月11日），所有樣本均由西南林業大學劉江華副教授分別鑒定柏科刺柏屬植物圓柏L.、垂枝香柏W. C. Cheng ex Ferré、昆明柏W. C. Cheng和鋪地柏(Endlicher) Siebold ex Miquel。樣本送至安諾優達基因科技（北京）有限公司進行葉綠體基因組測序并將測序結果上傳至NCBI數據庫（https://www.ncbi. nlm.nih.gov/），圓柏、昆明柏、垂枝香柏、鋪地柏登錄號分別為MZ151419、MZ151420、MZ151421、MZ151422。

2 方法

2.1 數據處理方法

剔除CDS中基因長度小于300 bp、重復和不完整的序列，選擇起始密碼子為ATG、終止密碼子為TAA、TAG、TGA的序列，對符合條件的CDS進行下一步的分析。以4種刺柏屬植物篩選得到的CDS序列進行分析，使用CodonW1.4.2、CUSP、SPSS26.0、Excel、CUSP（http://imed.med.ucm.es/ EMBOSS/）等軟件進行數據分析以及圖表繪制。

2.1.1 ENC-plot繪圖分析 以GC3為橫坐標（），有效密碼子數（effective number of codon，ENC）值為縱坐標（），繪制散點圖和ENC值的標準曲線圖（ENC＝2＋GC3＋29/[GC32＋(1－GC3)2]）。若基因落點位于標準曲線上或附近表明密碼子偏好主要受突變影響，反之則主要受到選擇影響。同時為了更準確地評估觀察值（ENCobs）與期待值（ENCexp）之間的差異，根據公式（ENCexp－ENCobs）/ENCexp計算ENC的比值頻數，并結合比值頻數對差異進行量化分析。

2.1.2 PR2-plot分析 以G3/(G3＋C3)為橫坐標（），以A3/(A3＋T3)為縱坐標（），繪制散點圖，圖中心點位置表明堿基含量A＝T、G＝C，即密碼子不具有偏好性，其余點與中心點的矢量距離則代表其偏倚程度和方向。

2.1.3 中性繪圖分析 以GC3為橫坐標（），GC12（GC1、GC2的平均值）為縱坐標（），繪制散點圖并做直線擬合分析。若直線斜率接近1，表明密碼子在堿基組成上無差異，則密碼子偏好性主要決定因素為突變；若直線斜率接近0，則更多地受到選擇的影響。

2.2 高表達基因庫和低表達基因庫

以密碼子的ENC值作為偏好性參考的重要指標，選取升序排列后的ENC值的前后兩端10%的基因，并依此建立高、低表達基因庫。根據CodonW軟件計算2個表達庫中密碼子的相對同義密碼子使用度（relative synonymous codon usage，）值和RSCU值（RSCU＝RSCU高表達－RSCU低表達），確定4種刺柏屬植物葉綠體基因組的高表達密碼子（RSCU＞1）和高頻密碼子（?RSCU≥0.08），同時滿足2個條件的密碼子為最優密碼子。

2.3 密碼子偏性驗證分析

從NCBI數據庫中下載擬南芥L.、香子含笑L.、杉木L.的葉綠體基因組序列，在4種刺柏屬葉綠體基因組基礎上篩選相同基因的CDS并計算密碼子的RSCU值，進行外類群聚類分析，以此來驗證4種刺柏屬葉綠體基因組密碼子偏好性。

3 結果與分析

3.1 密碼子組成特征分析

通過CodonW1.4.2和在線CUSP程序對圓柏（MZ151419）、昆明柏（MZ151420）、垂枝香柏（MZ151421）、鋪地柏（MZ151422）葉綠體基因組進行分析（表1）。對密碼子各位置的GC含量和ENC值進行相關性分析（表2），并對4種刺柏屬葉綠體基因組密碼子參數分析（表3）。結果表明GCall和GC1、GC2、GC3的相關性均達到極顯著水平（0.543～0.774）；GC3與GC1、GC2的相關性不顯著，但與GCall（0.543～0.743）和ENC（0.391～0.552）呈顯著相關性，說明密碼子的第3位與第1、2位的堿基組成不同。

表1 4種刺柏屬植物葉綠體基因組

表2 4種刺柏屬植物葉綠體基因的GC含量及ENC值

表3 4種刺柏屬葉綠體基因組密碼子參數的相關性分析

*< 0.05**< 0.01

3.2 密碼子偏好性影響因素分析

3.2.1 ENC-plot分析 通過對葉綠體基因組中各基因的ENC值和GC3的關聯分析，然后根據各基因與標準曲線的距離來衡量影響密碼子偏好性的主要因素。落點位于標準曲線上或附近表明密碼子偏好性主要受突變影響，反之則主要受到選擇影響。

根據ENC-plot分析圖（圖1），大部分落點位于標準曲線下方，極少部分位于標準曲線下方較遠位置。ENC值在31.65～59.39，且絕大部分ENC值大于45，說明這4種刺柏屬植物葉綠體基因組密碼子偏好性均較弱。結合ENC比值頻數分布表（表4）發現，4種刺柏屬葉綠體基因組分別有14、16、17、14個基因位于標準曲線附近，分別有15、22、25、23個基因距標準曲線較遠，表明突變影響了密碼子的使用模式。

3.2.2 PR2-plot分析 以G3/(G3＋C3)為橫坐標，A3/(A3＋T3)為縱坐標進行PR2-plot分析來研究葉綠體基因組密碼子的第3位上（A/T）與（C/G）之間的關系。結果表明在4個區域內的點分布并不均勻（圖2）。4種刺柏屬植物的點都集中在右半部分，這說明堿基使用情況T＞A，G＞C。結果表明4個刺柏屬葉綠體基因組密碼子在第3位上A/T的使用頻率較高，且密碼子使用模式受到多種因素的影響。

圖1 ENC-plot分析圖

表4 ENC比值頻率分布

圖2 PR2-plot分析

3.2.3 中性繪圖分析 本研究對4種刺柏屬植物葉綠體基因組做中性繪圖分析以進一步確定影響其密碼子偏好性的因素。以GC3為橫坐標（），GC12為縱坐標（）繪制散點圖。若基因落點沿對角線分布，斜率接近1，則3個位置上的密碼子堿基組成無明顯差異，密碼子偏好性主要受突變影響，反之則為自然選擇影響。

如圖3所示，4個葉綠體基因組密碼子的GC12值在0.33～0.52，GC3值在0.15～0.39，斜率在0.010～0.049，相關性不顯著，說明密碼子1、2位堿基和第3位堿基的組成存在明顯差異，密碼子的偏好性受自然選擇影響較大，這與ENC-plot分析結果一致。

圖3 中性繪圖分析

3.3 最優密碼子確定

通過構建高表達、低表達基因庫，分別計算4種刺柏屬植物在這2個基因庫中的RSCU和?RSCU，見表5。選取RSCU＞1的密碼子為高頻密碼子，?RSCU≥0.08的密碼子為高表達密碼子。同時滿足高頻、高表達的密碼子選為最優密碼子。結果表明（表6），這4個刺柏屬葉綠體基因組中以G結尾的最優密碼子只占極小部分，以A和U結尾的密碼子占絕大部分，其中又以U結尾的密碼子數最多。

表5 4種刺柏屬植物葉綠體基因組密碼子RSCU值

3.4 外類群聚類分析

為驗證本研究中所得出的4種刺柏屬葉綠體基因組密碼子的使用偏性能否用于后續對密碼子的研究，使用擬南芥（NC-000932.1）、香子含笑（MW470947）、杉木（KC427270.1）公開發表的葉綠體基因組數據篩選相同基因，計算密碼子的RSCU值（表7），根據RSCU值做聚類分析。結果表明（圖4），7個物種可以分為3大類：擬南芥、香子含笑和針葉樹種。其中鋪地柏單獨聚為1支，表明鋪地柏與其他3種刺柏屬植物密碼子使用模式具有明顯的差異。

4 討論

通常密碼子第3位堿基受到的自然選擇壓力較小，不會影響與氨基酸的對應關系，因此對研究密碼子偏好性具有重要意義[19]。對4種刺柏屬植物葉綠體基因組進行相關性分析、中性繪圖分析，結果顯示GC3和GC1、GC2相關性不顯著，且GC3＜GC1、GC2，表明密碼子GC3與GC1、GC2有差異，且密碼子偏向以A/G堿基結尾，這與巴山松[20]、臭柏[21]等大部分裸子植物葉綠體基因組密碼子偏好性的研究結論相似，表明親緣關系近的高等植物密碼子偏好使用模式具有一定的相似性。

表6 4種刺柏屬植物葉綠體基因組最優密碼子

*＜0.05 **＜0.01 ***＜0.001

表7 7個物種葉綠體基因組密碼子RSCU值

續表7

圖4 外類群聚類分析圖

影響密碼子偏好性形成的因素有很多，這其中最主要的因素是自然選擇和突變[22]。根據ENC-plot分析、PR2-plot分析的結果顯示4種刺柏屬植物葉綠體基因組密碼子的使用頻率T＞A、G＞C，說明影響密碼子偏好性形成的主要因素是自然選擇，這與馬尾松[23]、思茅松[24]等針葉樹種的結論一致；但與杜梨[25]、銀白楊[26]等主要受到突變影響以及巨桉[27]、芒果[28]等受突變和自然選擇的影響相對均衡的結論不一致，推測可能存在多種因素影響植物密碼子偏好性，不同物種間影響密碼子偏好性的主要因素也有所不同[29]。

利用與外類群葉綠體基因組密碼子的RSCU值進行聚類分析來對4種刺柏屬植物葉綠體基因組密碼子的使用偏性結果進行群體驗證。結果顯示，4種刺柏屬植物中鋪地柏單獨聚為1支；昆明柏、圓柏、垂枝香柏聚為另外1支。鋪地柏原產于日本，為匍匐灌木；垂枝香柏和昆明柏是中國特有種，為喬木或小喬木；圓柏雖分布廣泛，但屬于喬木[30]。研究結果與4種刺柏屬植物的地理分布及形態學特征一致，證明本研究中4種刺柏屬植物的密碼子偏好性分析的結果可用于后續對異源基因表達載體的設計和基因工程的應用和研究中。

本研究在4種刺柏屬植物葉綠體基因組密碼子中篩選出共65個最優密碼子，確定了GUA、UCU、CCU、CCA、GCU、AAA、GAU、CGU、GGA等9個密碼子為4種刺柏屬葉綠體基因組共有最優密碼子，大部分最優密碼子以A或U結尾，尤其以U結尾，這與杉木[31]、云南藍果樹[32]、梧桐[33]等高等植物研究結果類似，說明高等植物葉綠體基因組的進化相對保守。

本研究對4個刺柏屬葉綠體基因組密碼子的特征及偏性進行了分析和驗證，找出最優密碼子，對制定密碼子優化策略來提高異源蛋白在刺柏屬植物中的表達量從而改良刺柏屬植物重要性狀奠定堅實基礎。

利益沖突 所有作者均聲明不存在利益沖突

[1] 張悅. 紅豆杉科葉綠體基因組結構多樣化研究 [D]. 北京: 北京林業大學, 2019.

[2] 任婷. 十種報春花屬植物的葉綠體基因組研究 [D]. 西安: 西北大學, 2018.

[3] Li S F, Su T, Cheng G Q,. Chromosome evolution in connection with repetitive sequences and epigenetics in plants [J].(), 2017, 8(10): E290.

[4] 姜汶君, 郭夢月, 龐曉慧. 葉綠體基因組在藥用植物鑒定及系統進化研究中的應用 [J]. 世界中醫藥, 2020, 15(5): 702-708.

[5] 丁彥強, 方揚, 靳艷玲, 等. 基于葉綠體基因組的浮萍亞科系統進化 [J]. 應用與環境生物學報, 2017, 23(2): 215-219.

[6] 李泳潭, 張軍, 黃亞麗, 等. 杜梨葉綠體基因組分析 [J]. 園藝學報, 2020, 47(6): 1021-1032.

[7] 母連勝, 何勇, 田志宏. 植物葉綠體遺傳轉化技術及應用研究進展 [J]. 長江大學學報: 自科版, 2017, 14(14): 52-57.

[8] 周承哲, 朱晨, 李小楨, 等. 茶樹密碼子使用偏性分析方法及其研究進展 [J]. 分子植物育種, 2020, 18(5): 1480-1488.

[9] 馬孟莉, 張薇, 孟衡玲, 等. 豆蔻屬藥用植物葉綠體基因組密碼子偏性分析 [J]. 中草藥, 2021, 52(12): 3661-3670.

[10] Yu C H, Dang Y K, Zhou Z P,.usage influences the local rate of translation elongation to regulate co-translational protein folding [J]., 2015, 59(5): 744-754.

[11] 吳憲明, 吳松鋒, 任大明, 等. 密碼子偏性的分析方法及相關研究進展 [J]. 遺傳, 2007, 29(4): 420-426.

[12] 唐曉芬, 陳莉, 馬玉韜. 密碼子使用偏性量化方法研究綜述 [J]. 基因組學與應用生物學, 2013, 32(5): 660-666.

[13] 楊云彭, 馬曉焉, 霍毅欣. 密碼子優化策略在異源蛋白表達中的應用 [J]. 生物工程學報, 2019, 35(12): 2227-2237.

[14] 符文豪. 刺柏木材中抗炎活性成分的研究 [D]. 廈門: 廈門大學, 2019.

[15] 趙俊淇. 柏木、圓柏的精油提取及其抑菌性研究 [D]. 雅安: 四川農業大學, 2018.

[16] Noorbazargan H, Amintehrani S, Dolatabadi A,. Anti-cancer & anti-metastasis properties of bioorganic-capped silver nanoparticles fabricated fromextract against lung cancer cells [J]., 2021, 11(1): 61.

[17] 張翼, 來靜. 內蒙古中西部地區圓柏移植季節分析與選擇 [J]. 農業與技術, 2021, 41(5): 109-111.

[18] 施辰陽, 沙旭明, 馬丹丹, 等. 基于外源硅處理下刺柏對重金屬鎘污染土壤的修復治理研究 [J]. 杭州師范大學學報: 自然科學版, 2019, 18(1): 46-52.

[19] 王曉雙，楊芳，羅茜，等. 石斛屬葉綠體基因組密碼子使用偏性及系統發育分析[J/OL].分子植物育種, 2022, http://kns.cnki.net/kcms/detail/46.1068.S.20210303.1340.012.html.

[20] 徐貝. 巴山松及其近緣種的葉綠體基因組比較分析 [D]. 西安: 西北大學, 2021.

[21] 路東曄, 張磊, 郝蕾, 等. 臭柏葉綠體基因組結構與系統進化分析 [J]. 西北植物學報, 2018, 38(8): 1464-1475.

[22] Rao Y S, Wu G Z, Wang Z F,. Mutation bias is the driving force ofusage in thegenome [J]., 2011, 18(6): 499-512.

[23] 葉友菊, 倪州獻, 白天道, 等. 馬尾松葉綠體基因組密碼子偏好性分析 [J]. 基因組學與應用生物學, 2018, 37(10): 4464-4471.

[24] 原曉龍, 陳偉, 王毅, 等. 思茅松葉綠體基因組密碼子偏好性分析 [J]. 基因組學與應用生物學, 2021, 40(3): 1278-1283.

[25] 辛雅萱, 董章宏, 瞿紹宏, 等. 杜梨葉綠體基因組密碼子偏好性分析 [J]. 河北農業大學學報, 2020, 43(6): 51-59.

[26] Zhou M, Long W, Li X. Analysis of synonymoususage in chloroplast genome of[J]., 2008, 19(4): 293-297.

[27] 王鵬良, 吳雙成, 楊利平, 等. 巨桉葉綠體基因組密碼子偏好性分析 [J]. 廣西植物, 2019, 39(12): 1583-1592.

[28] 唐玉娟, 趙英, 黃國弟, 等. 芒果葉綠體基因組密碼子使用偏好性分析 [J]. 熱帶作物學報, 2021, 42(8): 2143-2150.

[29] 秦政, 鄭永杰, 桂麗靜, 等. 樟樹葉綠體基因組密碼子偏好性分析 [J]. 廣西植物, 2018, 38(10): 1346-1355.

[30] 中國科學院中國植物志編輯委員會. 中國植物志(第7卷) [M]. 北京: 科學出版社, 1978: 362-363.

[31] 鄭薇瑋, 陳金慧, 郝兆東, 等. 杉木等5種針葉樹葉綠體密碼子偏好性分析 [J]. 分子植物育種, 2016, 14(5): 1091-1097.

[32] 原曉龍, 康洪梅, 王毅. 云南藍果樹葉綠體基因組密碼子偏好性分析 [J]. 西北林學院學報, 2020, 35(4): 26-31.

[33] 陸奇豐, 駱文華, 黃至歡. 兩種梧桐葉綠體基因組密碼子使用偏性分析 [J]. 廣西植物, 2020, 40(2): 173-183.

Analysis of codon bias in the chloroplast genome of four medicinal plants of

SHU Jun-xia1, YANG Lin1, ZHOU Tao1, SHEN Lian-wen2, WANG Da-wei1

1. Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming 650224, China 2. Key Lab of Ministry of Education for Conservation and Utilization of Mountain Forest Resources in Southwest China, Southwest Forestry University, Kunming 650224, China

To determine the codon usage preference of the chloroplast genome of four medicinal plants ofgenus (,,, and) and the factors that affecting its codon preference.CodonW, CUSP, SPSS and other software was used to analyze the codon preference of the chloroplast genome of four medicinal plants of.The GC content of the chloroplast genome codons of the four medicinal plants ofwas from 35.8% to 37.3%, and the effective number of codon (ENC) value was from 47.10 to 47.79, indicating that their codon preference was weak. Neutral plotting, ENC-plot, PR2-plot analysis showed that the factors affecting the codon usage bias of the chloroplast genome of the four medicinal plants ofspecies were selection and mutation. At the same time, the relative synonymous codon usage (RSCU) value and ENC value were used to filter out a total of 65 optimal codons, of which nine codons are the optimal codons shared by four species of.The codon preference of four species ofspecies uses A/T ending, and its codon usage bias is mainly affected by natural selection. This study analyzes the codon usage bias of the chloroplast genome of four species ofplants and verified it by cluster analysis of the outer group cluster analysis, revealing the main factors affecting the codon usage bias ofgenus, and it is the foreign protein of the subsequentchloroplast genome. The construction and optimization of expression vectors provide a theoretical basis.

L.;L.;W. C. Cheng ex Ferré;W. C. Cheng;(Endlicher) Siebold ex Miquel; codon preference; optimal codon; natural selection

R286.12

A

0253 - 2670(2022)23 - 7507 - 09

10.7501/j.issn.0253-2670.2022.23.022

2022-03-01

云南森林資源培育與利用協同創新中心開放基金項目（201903611）

舒軍霞（1998—），女，碩士研究生，研究方向為林木遺傳育種。Tel: 18008713477 E-mail: shujunxia@swfu.edu.cn

通信作者：王大瑋，碩士生導師，副教授，研究方向為林木遺傳育種。Tel: 13888915161 E-mail: wangdawei@swfu.edu.cn

[責任編輯 時圣明]