‘懷玉山’高山馬鈴薯葉綠體基因組特征及密碼子使用偏好性分析

洪森榮，張牧彤，徐子林，張欽榮，羅雨欣，田文慧，王心雨

（1.上饒師范學院 生命科學學院，江西 上饒 334001；2.上饒農業技術創新研究院，江西 上饒 334001；3.上饒師范學院 上饒市藥食同源植物資源保護與利用重點實驗室，江西 上饒 334001；4.上饒師范學院 上饒市薯芋類作物種質保存與利用重點實驗室，江西 上饒 334001）

葉綠體是高等植物細胞內一種重要的與光合作用和物質代謝相關的細胞器，葉綠體基因組是一套具有母系遺傳特征的獨立基因組，是高等植物細胞質基因組的組成成分之一[5]。與核基因組相比，葉綠體基因組全長序列短、易測序獲得、基因直系同源、基因結構穩定、保守性較高、進化速率適中，目前已經廣泛應用于植物系統發育分析、物種分類鑒定及分子標記開發等研究中，在物種起源、進化、演變及比較基因組學等研究領域發揮著越來越大的作用[6]。密碼子是核酸和蛋白質之間遺傳信息傳遞的橋梁[7]，mRNA 上的遺傳信息以tRNA 三重密碼子傳遞。氨基酸一般對應≥1 的密碼子[8]，這些密碼子稱為同義密碼子[9]。在自然選擇或突變偏好的情況下，基因傾向于使用≥1 的同義密碼子，即同義密碼子使用偏好性[10-12]。目前，關于茄屬的葉綠體基因組研究已有報道[13-19]，而針對‘懷玉山’高山馬鈴薯的研究大多集中在基因克隆[20]、轉錄組分析[3]、遺傳多樣性[4]、脫毒快繁[2]、DNA 甲基化敏感擴增多態性(MSAP)分析[21]等方面，對‘懷玉山’高山馬鈴薯葉綠體全基因組及其密碼子使用偏好性方面的研究還未見系統報道。本研究通過對‘懷玉山’高山馬鈴薯葉綠體基因組進行測序和組裝，分析基于葉綠體基因組的‘懷玉山’高山馬鈴薯系統進化、結構解析和密碼子偏好性等，為‘懷玉山’高山馬鈴薯葉綠體基因組研究和應用提供科學依據，也為進一步研究‘懷玉山’高山馬鈴薯遺傳背景、種質資源保護與開發利用奠定基礎。

1 材料與方法

1.1 材料

由上饒市薯芋類作物種質保存與利用重點實驗室提供的‘懷玉山’高山馬鈴薯試管苗。

1.2 方法

1.2.1 DNA 提取和測序 選取‘懷玉山’高山馬鈴薯(MLS)試管苗葉片組織，利用植物基因組DNA 提取試劑盒(北京天根生化科技有限公司)提取‘懷玉山’高山馬鈴薯試管苗DNA，質量分數為1%瓊脂糖凝膠電泳檢測DNA 的完整性，NanoDrop 2000 分光光度計(Thermo Scientific 公司)檢測 DNA 濃度和純度，用超聲波將DNA 片段化，然后對片段化的DNA 進行片段純化、末端修復、3′端加A、連接測序接頭，再用瓊脂糖凝膠電泳進行片段大小選擇，進行聚合酶鏈式反應(PCR)擴增形成測序文庫。建好的文庫先進行文庫質檢，質檢合格的文庫用BGISEQ-500 平臺進行測序。

1.2.2 葉綠體全基因組的組裝與注釋 通過SOAPnuk 1.3.0 對raw data (測序下機的原始數據)進行數據過濾，去除其中的接頭序列及低質量reads (高通量測序中一個反應獲得的測序序列)，獲得高質量的clean data (對原始數據進行過濾后并剔除了低質量數據的剩余數據)。采用Noveplastys 軟件組裝葉綠體基因組核心模塊，以起始組裝序列為起點開始組裝葉綠體contigs (很多reads 根據序列拼接在一起拼出的片段)，如果contigs 未環化，則利用CAP 3 軟件連接多個contigs 為完整葉綠體基因組，并手動調整環狀葉綠體基因組起始位置。使用GeSeq、tRNAscan-SE 對葉綠體基因組進行注釋，再經過手工校正后得到最終的基因注釋結果。將注釋完成的‘懷玉山’高山馬鈴薯葉綠體基因組序列提交至美國國家生物信息中心(NCBI)，獲得登錄號：OP589401。使用OGDRAW 繪制葉綠體基因組圖譜。

1.2.3 葉綠體基因組特征分析 通過JSHYCloud 在線工具集分析并統計葉綠體基因組、大單拷貝區(LSC)、小單拷貝區(SSC)和反向重復區(IR)的鳥嘌呤和胞嘧啶所占的比例(GC 比例)；使用MISA 軟件進行簡單重復序列(SSR)分析，單核苷酸、二核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸的最小重復值分別設置為10、6、5、5、5、5；利用REPuter 軟件進行長重復序列(longrepeat)分析，查找正向重復(F)、反向重復(R)、互補重復(C)、回文重復(P)等4 種重復類型；通過Pasteur Galaxy 在線工具集中的CodonW 模塊分析密碼子使用情況，設置輸出結果為有效密碼子數(ENC)和相對同義密碼子使用頻率(RSCU)，其他參數設為默認值。將‘懷玉山’高山馬鈴薯葉綠體基因組序列上傳至美國國家生物技術信息中心(NCBI) 進行BLASTn 比對，選擇highly similar sequence (megablast)比較相似性在95 %以上的序列，檢索獲得‘懷玉山’高山馬鈴薯的近緣種。利用Gview、VISTA tools、IRscope 和DNADnaSP 6.0 軟件繪制‘懷玉山’高山馬鈴薯及其10 個近緣種(S.cochoaeNC_062512、多毛番茄S.habrochaitesNC_026879、潘那利番茄S.pennelliiNC_035742、S.bukasoviiMT120867、S.bolivienseNC_062870、S.trisectumNC_062469、S.salamancaeNC_062480、S.clivorumNC_062513、S.mortoniiNC_062426、S.insanumMW384851)的變異圈圖、mVIST 結構變異圖、IR 結構變異圖，計算‘懷玉山’高山馬鈴薯及其10 個近緣種的基因組核酸多樣性(Pi)，參數設置100 bp 滑窗，25 bp 的步長，并進行中性繪圖分析(GC3-GC12 分析)、ENC-plot 分析、PR2-bias-plot 分析和最優密碼子分析；對‘懷玉山’高山馬鈴薯葉綠體基因的ENC 進行排序，分別選取兩端基因各5 個，構建高表達基因庫(ENC 小)和低表達基因庫(ENC 大)，并計算兩者的RSCU 差值(ΔRSCU)。篩選ΔRSCU≥0.08 的高表達密碼子，且將RSCU＞1.00 的高頻率密碼子定義為‘懷玉山’高山馬鈴薯葉綠體基因組的最優密碼子；最后利用mafft 7.0 和fasttree 2.1.10 軟件分別對‘懷玉山’高山馬鈴薯和18 個近緣種以及煙草屬Nicotiana2 個外類群物種進行序列比對和構建進化樹。

2 結果與分析

2.1 葉綠體基因組序列特征分析

經過測序組裝的完整的葉綠體基因組長度為155 296 bp，圖1 顯示：‘懷玉山’高山馬鈴薯葉綠體基因組呈典型的四分體結構，包含1 個LSC、1 個SSC 和2 個將LSC 與SSC 分隔開的IR (IRa 和IRb)?；蚪M的總GC 比例為37.88%，A、T、C、G 比例分別為30.65%、31.47%、19.24%、18.65%。LSC、SSC 和IR 的長度分別為85 737、18 373、25 593 bp。LSC 的GC 比例為36.01%，A、T、C、G 比例分別為31.29%、32.70%、18.40%、17.61%；SSC 的GC 比例為32.09%，A、T、C、G 比例分別為33.78%、34.14%、16.69%、15.40%；IRb 的GC 比例為43.10%，A、T、C、G 比例分別為28.57%、28.33%、20.72%、22.39%；IRa 的GC 比例為43.10%，A、T、C、G 比例分別為28.33%、28.57%、22.39%、20.72%。表明‘懷玉山’高山馬鈴薯IR 的GC 比例最大，LSC 次之，SSC 最少；葉綠體基因組總GC 比例顯著低于AT 比例；葉綠體基因組各堿基比例從大到小依次為T、A、C、G。

圖1 ‘懷玉山’高山馬鈴薯葉綠體基因組圖譜Figure 1 Chloroplast genome map of S.tuberosum var.cormosus ‘Huaiyushan’

2.2 葉綠體基因類型分析

葉綠體基因組共注釋到光合作用基因、自我復制基因、其他基因和未知功能基因4 類，包括87 個編碼區(CDS)基因、37 個tRNA 基因、8 個rRNA 基因、1 個假基因，共133 個基因。對有多個外顯子的葉綠體基因進行結構分析，由2 個外顯子構成的基因有21 個，包括13 個CDS 基因和8 個tRNA 基因；由 3 個外顯子構成的基因有4 個，為clpP1、ycf3、rps12 (2 個)基因。LSC 的基因數量最多(81 個)，其中CDS 基因59 個、tRNA 基因22 個；SSC 的基因數量為11 個，其中CDS 基因10 個、tRNA 基因1 個；IR 的基因數量為17 個，其中CDS 基因6 個、rRNA 基因4 個、tRNA 基因7 個；SSC 與IRb 邊界(JSB)的基因數量為2 個(ndhF和ycf1)；LSC 與IRb 邊界(JLB)的基因數量為1 個(rps19)；SSC 與IRa 邊界(JSA)的基因數量為2 個(ycf1)；LSC 與IRa 邊界(JLA)的基因數量為0。rps12 有2 個拷貝，每個拷貝具有3 個外顯子，且2 個拷貝共享第1 個外顯子，第1 個外顯子位于LSC，另外2 個外顯子位于IR (表1)。

表1 ‘懷玉山’高山馬鈴薯葉綠體基因功能分類Table 1 Chloroplast gene functional classification of S.tuberosum var.cormosus ‘Huaiyushan’

2.3 葉綠體基因組重復序列分析

葉綠體基因組中共檢測到38 個SSR 位點，其中，單堿基重復有36 個，雙堿基重復有2 個。其中，重復單元為A/T，重復頻率為10 的SSR 位點數量最多(18 個)，重復頻率為11 的SSR 位點數量次之(11 個)；重復單元為AT/AT、重復頻率為6 的SSR 位點數量為2 個。

四川大學陳劍等[11]從醬香型酒糟中篩選得到1株耐酸產酯細菌ZP-28，此株細菌具有顯著耐高溫耐酸特性，且酯化力較高，將其應用在麩曲中具有明顯改善麩曲品質的作用。本課題在此基礎上以該菌株為研究對象，采用分子生物學方法對其進行鑒定，并利用活菌計數法和分光光度法兩種不同的方法測定其生長曲線，采用單因素實驗考察了溫度、培養時間、接種量等培養條件對酯化力的影響，并在此基礎上采用正交法對上述發酵條件進行了優化，對改進制曲工藝，提高白酒酯香提供了一定的理論依據。

2.4 葉綠體基因組長重復序列分析

葉綠體基因組共鑒定到32 個長重復序列，包括16 個正向重復(15 個30～39 bp，1 個40～49 bp)，16 個回文重復 (13 個30～39 bp，2 個40～49 bp，1 個50～59 bp)，無反向重復和互補重復。

2.5 葉綠體基因組反向重復區邊界結構差異分析

‘懷玉山’高山馬鈴薯及其10 個近緣種葉綠體基因組結構從LSC 中間呈線性展開，均由1 個LSC、1 個SSC 和2 個IR (IRa 和IRb) 4 部分組成?！畱延裆健呱今R鈴薯及其10 個近緣種rpl22、rps19、rpl2、ycf1、ndhF、trnH和psbA位置基本一致，但收縮和擴張的長度存在一些差異(圖2)。

圖2 ‘懷玉山’高山馬鈴薯及其10 個近緣種葉綠體基因組大單拷貝區、小單拷貝區和反向重復區邊界位置的比較Figure 2 Comparison on the boundary locations of large single copy region, small single copy region and inverted repeat region in chloroplast genomes of S.tuberosum var.cormosus ‘Huaiyushan’ and its 10 related species

2.6 葉綠體基因組比對分析

‘懷玉山’高山馬鈴薯及其10 個近緣種葉綠體基因組核苷酸多樣性的變化范圍為0～0.139 27，高變區主要分布在LSC 和SSC。LSC 的trnL-UAA-trnF-GAA、cemA、rps12-exon1-clpP1、clpP1 基因變異率最高；SSC 的rpl32-trnL-UAG、ycf1 基因變異率最高。

2.7 葉綠體基因組密碼子使用偏性分析

2.7.1 同義密碼子的偏性分析 葉綠體基因組87 個CDS 基因密碼子3 個位置GC 比例的平均值為38.38%，GC1、GC2、GC3 分別為45.98%、39.55%、29.60%，這說明GC 在密碼子3 個位點上的分布存在顯著差異，只有GC2 與平均GC 大致接近(圖3)。ENC 是密碼子偏性分析的重要指標，通常將35 作為區分值來評估密碼子偏倚的強度。葉綠體基因組87 個CDS 基因的平均ENC 為47.29，ENC＞45 的基因有60 個，ENC＞35 的基因有83 個，有4 個基因的ENC＜35，這表明葉綠體基因組的密碼子偏性較弱。通過SPSS 20.0 進行相關性分析，結果表明：密碼子總GC 比例(GCall)與GC1、GC2 在0.01 水平上均存在極顯著的正相關，GCall 與GC3 在0.05 水平上顯著相關；GC1 與GC2 在0.05 水平上存在顯著正相關，但兩者均與GC3 不相關。這表明葉綠體基因組密碼子前2 位的堿基組成相似，而與第3 位不相似。ENC 與GC1、GC2、GC3 均不相關，說明密碼子上第1 位、第2 位和第3 位的堿基組成對ENC 沒有顯著影響。葉綠體基因組 87 個CDS 基因序列共有31 個RSCU＞1 的密碼子。在這31 個密碼子中，除AUG、UUG 外，其余都以A、U 結尾，表明A、U 堿基在密碼子最后位點上出現的頻率最高?！畱延裆健呱今R鈴薯葉綠體基因組密碼子偏好以A、U 結尾(表2)。

表2 ‘懷玉山’高山馬鈴薯葉綠體基因組同義密碼子的使用頻率Table 2 Relative synonymous codon usage (RSCU) of chloroplast genome of S.tuberosum var.cormosus ‘Huaiyushan’

圖3 ‘懷玉山’高山馬鈴薯及其10 個近緣種葉綠體基因組密碼子組成成分分析Figure 3 Composition analysis of chloroplast genome codons of S.tuberosum var.cormosus ‘Huaiyushan’ and its 10 related species

2.7.2 GC3-GC12 分析 分析發現：‘懷玉山’高山馬鈴薯及其10 個近緣種植物葉綠體基因的GC3 比例分布為0.142 9～0.443 2，GC12 比例分布為0.285 7～0.658 5，兩者大多沿對角線上方分布。兩者的相關系數(r)為0.110 1 (R2=0.012 1)，相關不顯著(P＞0.05)，回歸斜率為0.117 5，說明GC12 與GC3 不相關(圖4A)。表明‘懷玉山’高山馬鈴薯葉綠體基因組密碼子使用偏性很大程度上受自然選擇的影響，而受突變壓力的影響小。

圖4 ‘懷玉山’高山馬鈴薯及其10 個近緣種葉綠體基因組密碼子GC3-GC12 分析(A)、ENC-plot 分析(B)和PR2-plot 分析(C)Figure 4 GC3-GC12 analysis (A), ENC-plot analysis (B) and PR2-plot analysis (C) of chloroplast genome codons of S.tuberosum var.cormosus‘Huaiyushan’ and its 10 related species

2.7.3 ENC-plot 分析 分析表明：分布在期望曲線上或曲線附近的基因較少，分布在期望曲線下方且遠離曲線的基因較多，說明大部分基因的實際ENC (ENCobs)與理論ENC (ENCexp)存在差異。為了解實際ENC 和理論ENC 的差異度，計算了‘懷玉山’高山馬鈴薯ENC 比值頻數，即(ENCexp-ENCobs)/ENCexp。結果表明：‘懷玉山’高山馬鈴薯葉綠體基因組基因中，有16.47%(14 個)的基因分布在0～0.1 區間，分布于期望曲線上或曲線附近，即ENCobs 接近于ENCexp 值，有83.53%的基因分布在0～0.1 區間外，遠離期望曲線分布，即ENCexp 和ENCobs 相差較大，表明自然選擇是影響‘懷玉山’高山馬鈴薯葉綠體基因組密碼子使用偏性的主要因素，而突變壓力的作用較小(圖4B)。

2.7.4 PR2-plot 分析 分析表明： A3/AU3 軸、G3/GC3 軸均以0.5 為界限，發現4 個平面內基因分布不均衡。從G3/GC3 軸看，多數基因位于上方(＞0.5)，少數基因位于下方(＜0.5)；從A3/AU3 軸看，多數基因位于左側(＜0.5)，少數基因位于右側(＞0.5)。這表明4 種堿基在同義密碼子第3 位上存在C＞G、T＞A 現象(圖4C)。當密碼子使用存在偏性完全受突變壓力影響時，C 和G 以及A 和T 同義密碼子在第3 位上的分布應相等。因此，‘懷玉山’高山馬鈴薯葉綠體基因組密碼子使用偏性主要受自然選擇等因素影響。

2.7.5 最優密碼子確定 RSCU 分析可知：同時滿足RSCU＞1 和ΔRSCU≥0.08 的密碼子共10 個，即CGU、AAA、CUU、GUU、GGA、GUA、GGU、UCA、GCU、CCU，這些密碼子都以A、U 結尾，被確定為‘懷玉山’高山馬鈴薯葉綠體基因組的最優密碼子(表3)。

表3 ‘懷玉山’高山馬鈴薯葉綠體基因最優密碼子篩選Table 3 Optimal codon screening of chloroplast genome of S.tuberosum var.cormosus ‘Huaiyushan’

2.8 系統發育分析

基于‘懷玉山’高山馬鈴薯和18 個近緣種以及煙草屬2 個外類群物種葉綠體基因組構建的系統發育樹分析可知：茄屬聚為一大類，煙草屬聚為另一大類。在茄屬中，MLS 與S.tuberosumNC_008096(‘Ddeiree’)聚為一小分支。說明‘懷玉山’高山馬鈴薯與S.tuberosum‘Ddeiree’親緣關系較近，兩者同源(圖5)。

圖5 基于葉綠體基因組的‘懷玉山’高山馬鈴薯及其18 個近緣種的系統發育樹Figure 5 Phylogenetic tree of S.tuberosum var.cormosus ‘Huaiyushan’ and its 18 related species based on chloroplast genome

3 討論與結論

葉綠體基因組結構保守、獨立母系遺傳，是被子植物基因組的重要組成部分，廣泛用于被子植物的生長發育、類群分析和進化分析[22]。被子植物葉綠體基因組大小一般為120～180 kb，IR 大小一般為20～30 kb[23-24]。在本研究中，‘懷玉山’高山馬鈴薯葉綠體基因組長度和IR 長度分別為155 296 和25 593 bp，與S.tuberosum‘Shepody’[16]葉綠體基因組長度和IR 長度一致，與其他馬鈴薯品種[13-15,17-18]相比，葉綠體基因組長度和IR 長度不超過500 bp，說明馬鈴薯各個品種的葉綠體基因組較為保守。

葉綠體的SSR 不僅與核基因組SSR 一樣，具有高多態性、多等位性、共顯性[25]，也具有單親遺傳模式，結構簡單、相對保守[26]，因此，葉綠體的SSR 有較好的種間、種內遺傳變異區分能力，已成為區分物種的重要分子標記而被廣泛應用[27]。關惜今等[13]研究表明：S.fernandezianum與其野生近緣種(S.phureja、S.palustre、S.etuberosum)葉綠體基因組中共檢測到36、36、42、40 個SSR，SSR 類型比較單一，只有單核苷酸和二核苷酸等2 種類型，單核苷酸為A和T等2 種類型，二核苷酸包括TA 和AT 等2 種類型，其數目比較少。本研究結果與此一致。在本研究中，在‘懷玉山’高山馬鈴薯葉綠體基因組中共檢測到38 個SSR 位點，其中，單堿基重復有36 個，雙堿基重復有2 個，較少的SSR 位點存在表明‘懷玉山’高山馬鈴薯葉綠體基因組可能不易發生重排。

IR 和單拷貝區(SC)邊界的膨脹和收縮被認為是被子植物葉綠體全基因組大小變化的主要機制[28]，同一屬不同品種葉綠體基因組IR/SC 邊界位置變化也不同[29]。關惜今等[13]研究表明：S.fernandezianum與其野生近緣種(S.phureja、S.palustre、S.etuberosum)葉綠體基因組rps19 基因均橫跨JLB，S.phureja的ndhF基因橫跨JSB，S.fernandezianum、S.palustre、S.etuberosum的ndhF基因均右移，分布在SSC，S.fernandezianum、S.phureja、S.palustre、S.etuberosum的ycf1 基因總長度為5 664 bp，均橫跨SSC 和IRa 區域。在本研究中，對‘懷玉山’高山馬鈴薯及其10 個近緣種葉綠體基因組 IR/SC 邊界區域的分析結果表明：這些葉綠體基因組的IR 都存在擴張或收縮的現象?！畱延裆健呱今R鈴薯的rps19 基因橫跨 JLB，橫跨 JLB 的左邊和右邊長度分別為209 和69 bp，在JSA，‘懷玉山’高山馬鈴薯ycf1 基因為5 663 bp，左邊和右邊長度分別為4 541 和1 122 bp。

許多植物存在密碼子偏好性(CUB)，即某一或幾種特定密碼子頻率超過其他同義密碼子。密碼子偏好性可用來評估基因組中蛋白質編碼區(CDS)的密碼子使用情況[30]。植物密碼子偏好性是物種不斷適應外界環境進化所導致的結果，生物獲得特定的密碼子使用模式以適應起源、進化、自然選擇和突變壓力等多種因素[31]。影響不同物種中密碼子偏好性差異的因素主要有堿基突變、基因表達水平、自然選擇等，自然選擇和突變壓力被認為是2 個最重要的因素[32-34]。密碼子第3 個堿基的同義突變不能改變氨基酸的類型，但被認為是決定氨基酸類型的重要特征，因此GC3 經常被用作密碼子偏向的重要指標[35-36]。本研究發現‘懷玉山’高山馬鈴薯葉綠體基因組的平均GC 比例為38.38%，GC3 為29.60%，更傾向于使用A/T 密碼子。RSCU 分析結果也證實了這一點?！畱延裆健呱今R鈴薯葉綠體基因組中存在A/T 密碼子使用偏向，這與大多數高等植物的模式一致[37]?！畱延裆健呱今R鈴薯葉綠體基因組平均ENC 為47.29，ENC＞35 的基因有83 個，有4 個基因的ENC＜35，表明‘懷玉山’高山馬鈴薯葉綠體基因組的密碼子偏性較弱。當密碼子的使用受到自然選擇的影響時，GC3 值往往分布在一個較小的范圍內，GC12 和GC3 之間沒有顯著的相關性[38]。密碼子偏好性可以通過調節基因翻譯的準確性和效率影響基因表達，基因表達水平越高，密碼子偏好性越強[39-40]。通過建立的高低基因表達庫，本研究挖掘到‘懷玉山’高山馬鈴薯葉綠體基因組10 個最優密碼子，即CGU、AAA、CUU、GUU、GGA、GUA、GGU、UCA、GCU、CCU，說明‘懷玉山’高山馬鈴薯葉綠體基因組密碼子更偏好于以A/U 結尾。篩選到的最優密碼子可以用于設計葉綠體基因表達載體，以提高葉綠體基因組中基因的表達水平，也可以利用已知密碼子的使用偏好來推測和預測未知基因的表達和功能，可為今后從遺傳水平上進行‘懷玉山’高山馬鈴薯育種改良提供參考。

含有足夠信息位點的葉綠體基因組已被證明可有效判斷系統發育關系，甚至是在較低的分類學水平下植物之間也有較強的分類學意義，為物種間系統發育的研究提供了新的思路[41]。在本研究中，在茄屬中‘懷玉山’高山馬鈴薯與S.tuberosum‘Desiree’單獨聚為一分支。說明‘懷玉山’高山馬鈴薯與S.tuberosum‘Desiree’親緣關系較近，表明兩者同源，推測‘懷玉山’高山馬鈴薯可能是S.tuberosum‘Desiree’從美國引種的。

綜上所述，本研究測序組裝了‘懷玉山’高山馬鈴薯葉綠體基因組全序列，分析了其編碼蛋白基因的密碼子使用特點，從高表達優越密碼子和高頻密碼子中選出兩者共有的密碼子，最終篩選得到了10 個葉綠體蛋白編碼基因的最優密碼子?！畱延裆健呱今R鈴薯密碼子的偏好性受到突變、選擇及其他多方面因素的共同影響，但自然選擇的影響更大，這為用基因工程手段改造外源基因密碼子，提高其在‘懷玉山’高山馬鈴薯葉綠體中的表達量提供了參考，也為在分子水平上研究茄科茄屬植物的系統進化提供參考。