小花吊兰叶绿体基因组密码子偏好性分析

    Analysis of Codon Bias in Chlorophytum laxum Chloroplast Genome

    • 摘要:
      目的 为明确药用植物小花吊兰(Chlorophytum laxum)叶绿体基因组中密码子的使用偏好性,分析自然选择与突变压力对密码子偏好性形成的影响,并筛选最优密码子,为后续叶绿体基因工程及系统进化研究提供理论依据。
      方法 从NCBI数据库获取小花吊兰叶绿体全基因组序列(OR078538),经严格过滤剔除短序列(< 300 bp)及含有中间终止密码子的序列后,获得49条有效蛋白编码序列(CDS)。利用CUSP、CodonW 1.4.2等生物信息学程序计算各位置碱基的GC含量(GC1、GC2、GC3和GCall)、有效密码子数(ENC)及相对同义密码子使用度(RSCU)。随后,结合中性绘图、ENC-plot以及PR2-plot等相关性分析模型,综合评估突变压力和自然选择在密码子偏好性形成中的相对贡献。最后,通过构建高、低表达基因库筛选确定该物种的最优密码子。
      结果 小花吊兰49条叶绿体CDS的平均总GC含量(GCall)为38.10%,且密码子不同位置的GC含量呈现明显的不均衡性:GC1(46.92%) > GC2(39.21%) > GC3(28.17%)。ENC值介于37.26~60.01之间,平均值为46.70(大于45的基因占65.31%),显示其整体密码子偏好性较弱。相关性与中性绘图分析显示,GC12与GC3无显著相关性(Radj2=-0.011);ENC-plot及PR2-plot图谱表明,绝大多数基因偏离标准预期曲线,且第3位碱基中G和T的使用频率分别高于C和A。综合筛选鉴定出17个最优密码子(如UUU、UUA、AUU等),均以A或U结尾。
      结论 小花吊兰叶绿体基因组密码子偏好性整体较弱,且强烈偏好使用以A或U结尾的同义密码子。这种偏好模式主要受自然选择压力的主导驱动,而基因突变压力的影响相对微弱。鉴定的17个最优密码子为小花吊兰及相关百合科植物的叶绿体基因工程密码子改造提供了参考。

       

      Abstract:
      Objective To clarify the codon usage bias in the chloroplast genome of the medicinal plant Chlorophytum laxum, analyze the impacts of natural selection and mutation pressure on the formation of codon bias, and screen the optimal codons, so as to provide a theoretical basis for subsequent chloroplast genetic engineering and phylogenetic studies.
      Method The complete chloroplast genome sequence of C. laxum (OR078538) was obtained from the NCBI database. After strict filtering to exclude short sequences (< 300 bp) and sequences containing internal stop codons, 49 valid protein-coding sequences (CDS) were obtained. Bioinformatics programs including CUSP and CodonW 1.4.2 were used to calculate the GC content at different codon positions (GC1, GC2, GC3, and GCall), the effective number of codons (ENC), and the relative synonymous codon usage (RSCU). Subsequently, combined with correlation analysis models such as neutrality plot, ENC-plot, and PR2-plot, the relative contributions of mutation pressure and natural selection to the formation of codon usage bias were comprehensively evaluated. Finally, the optimal codons of this species were identified by constructing high- and low-expression gene libraries.
      Result The average overall GC content (GCall) of the 49 chloroplast CDSs of C. laxum was 38.10%, and the GC content at different codon positions exhibited significant imbalance: GC1 (46.92%) > GC2 (39.21%) > GC3 (28.17%). The ENC values ranged from 37.26 to 60.01, with an average of 46.70 (genes with ENC > 45 accounted for 65.31%), indicating a weak overall codon usage bias. Correlation and neutrality plot analyses revealed no significant correlation between GC12 and GC3 (R2 = -0.011). The ENC-plot and PR2-plot indicated that the vast majority of genes deviated from the standard expected curve, and the usage frequencies of G and T at the third base position were higher than those of C and A, respectively. Through comprehensive screening, 17 optimal codons (e.g., UUU, UUA, AUU) were identified, all of which end in A or U.
      Conclusion The overall codon usage bias of the C. laxum chloroplast genome is relatively weak, with a strong preference for synonymous codons ending in A or U. This bias pattern is primarily driven by natural selection pressure, while the influence of gene mutation pressure is relatively weak. The 17 identified optimal codons can provide a reference for codon modification in the chloroplast genetic engineering of C. laxum and related Liliaceae plants.