芍药花色调控基因的密码子使用模式及其影响因素分析

doi:10.3864/j.issn.0578-1752.2016.12.012

摘要/Abstract

摘要： 【目的】芍药花色的优劣影响其观赏价值和商业价值，研究芍药花色调控基因的密码子使用偏好性和密码子使用模式的影响因素，为芍药花色调控基因在mRNA翻译、转基因设计、新基因表达与功能预测以及分子生物进化研究提供参考。【方法】根据前期芍药花色嵌合体品种‘金辉’转录组测序筛选的6 345个芍药花色调控基因，并根据CDS序列特征和大于300 bp原则进行过滤后最终获得的2 234个基因序列作为研究对象，利用Mobyle软件计算GC含量、第1与2位密码子的平均GC含量（GC12）、第3位密码子的GC含量（GC3s）、有效密码子数ENC、密码子适应指数CAI、相对同义密码子使用度RSCU等密码子偏性指标，其次进行中性绘图（GC12 vs. GC3）、ENC-GC3s绘图以及PR2（Parity Rule 2）绘图分析，并运用多元统计分析方法探讨突变压力和选择作用对密码子使用模式的影响程度，最后以5%CAI值作为高、低表达样本组，计算这两个样本组的同义密码子相对使用度，利用卡方检验Chi-square test分析两组之间的显著性差异来确定最优密码子。【结果】芍药花色相关基因的密码子GC3s含量为46.37%，大部分基因GC含量主要分布在30%—55%；中性绘图分析表明GC3s与GC12呈极显著的正相关（R²=0.202，P＜0.01）；ENC-GC3s绘图表明大部分基因分布在标准曲线周围，也有一部分基因分布在标准曲线下方较远的位置，同时大部分基因（ENCexp-ENCobs）/ENCexp比值集中分布在0.0—0.4；PR2绘图分析显示密码子第三位T的使用频率高于A，C使用频率高于G，表明嘌呤（A和G）与嘧啶（T和C）的使用频率并不均衡；对应性分析COA（Correspondence Analysis）表明，第一轴上显示了38.09%的差异，其他3个轴分别为18.42%、15.09%、14.59%，表明芍药花色调控基因的密码子使用模式评价以第一轴（Axis 1）为主；突变压力和选择作用分析发现，第一主轴与GC3s、CAI的相关系数均达到极显著正相关（R²=0.736，P＜0.01；R²=0.286，P＜0.01）；利用△RSCU和卡方显著性检验的方法，确定了21个为芍药花色调控相关基因的最优密码子，其中18个以G或C结尾，仅CGU、GGU等2个密码子以U结尾。【结论】芍药花色调控基因的最优密码子多数以G/C结尾，并且密码子使用模式主要受到碱基差异（R²=0.736）和基因表达水平（R²=0.286）共同作用的影响，其中碱基差异占主导因素。本研究了解了芍药花色调控基因的密码子使用模式情况，为通过密码子改造开展芍药花色遗传改良以及分子进化研究提供了一定的理论依据。

关键词: 芍药, 花色调控基因, 密码子使用模式, 影响因素

Abstract: 【Objective】 The quality of Paeonia lactiflora flower color affects its ornamental value and the commercial value of ornamental plants. This study aims to understand the codon usage pattern of genes regulating flower color and probe into the main factors affecting the formation of codon bias, which has important biological significance for mRNA translation, design of transgenes, the prediction of expression level and functions of new genes, and studies of molecular biology and evolution, etc. 【Method】In a previous study, 6,345 differential genes were screened out by transcriptome sequencing of a flower color chimaera cultivar “Jinhui” with a consistent genetic background red outer-petal and yellow inner-petal, followed by a further filtering analysis according to the principle of CDS sequence characteristics and greater than 300 bp. We finally obtained 2,234 genes as our research object. Mobyle software was used to calculate different parameters for the codon usage, such as GC content, average GC content of the first and second positions (GC12), GC content of the third position (GC3s), effective number of codon (ENC), codon adaptation index (CAI), and relative synonymous codon usage (RSCU). Further analysis of a neutrality plot (GC12 vs. GC3), an ENC-GC3s plot, and a Parity Rule 2 (PR2) plot were performed. Additionally, we probed into the influence of mutational pressure and translational selection by a multivariate statistical analysis. Finally, we took 5% CAI value as high-expression and low-expression sample groups, then calculated the RSCU value, and analyzed the significant difference to determine the optimal codons by a chi-square test.【Result】The results showed that the GC content at the third position of codons was 46.37%. The GC content of most genes was mainly distributed between 30% and 55%. Neutrality analysis showed that there was a significant positive correlation (R²=0.202, P＜0.01) between GC3s and GC12 value. The ENC-plot showed most of the genes on or close to the expected curve, but also some points with low-ENC values were below it. The(ENCexp-ENCobs)/ENCexp ratio of most genes ranged from -0.05 to 0.05. The Parity Rule 2-plot showed that the frequency of T nucleotide at the third position was higher than A, and C was higher than G, suggesting that the use frequencies of four nucleotide were not balanced. Correspondence analysis showed that the first axis showed a 38.09% variation, while the other three axes showed 18.42%, 15.09%, and 14.59%, respectively, suggesting that the first axis was the main index evaluating the codon usage bias of Paeonia lactiflora genes regulating flower color. Mutation pressure and selection analysis showed there were significant negative correlations (R²=0.736, P＜0.01. R²=0.286, P＜0.01) between the first axis and GC3s, CAI value, respectively. Using the delta RSCU and significant chi-square test methods, we defined 21 codons as the major preference codons in the Paeonia lactiflora genes regulating flower color, of 18 codons ending with G or C, only CGU and GGU ending with U.【Conclusion】In conclusion, most of optimal codons ended with G or C. Meanwhile, the codon usage pattern of Paeonia lactiflora genes regulating flower color is formed under the effect of mutational pressure (R²=0.736) and translational selection (R²=0.286), but a mutational bias was the major influence on codon usage. This study not only preliminarily reveals the codon usage pattern of Paeonia lactiflora genes, but also provides a certain theoretical basis for further carrying out genetic improvement of Paeonia lactiflora flower color by codon reconstruction and analyzing the molecular evolution.

Key words: Paeony, flower color controlling gene, codon usage pattern, influence factor

吴彦庆，赵大球，陶俊. 芍药花色调控基因的密码子使用模式及其影响因素分析[J]. 中国农业科学, 2016, 49(12): 2368-2378.

WU Yan-qing, ZHAO Da-qiu, TAO Jun. Analysis of Codon Usage Pattern of Paeonia lactiflora Genes Regulating Flower Color and Its Influence Factors[J]. Scientia Agricultura Sinica, 2016, 49(12): 2368-2378.

0
/ / 推荐

导出引用管理器 EndNote|Reference Manager|ProCite|BibTeX|RefWorks

链接本文: https://www.chinaagrisci.com/CN/10.3864/j.issn.0578-1752.2016.12.012

https://www.chinaagrisci.com/CN/Y2016/V49/I12/2368

参考文献

[1] Ikemura T. Codon usage and tRNA content in unicellular and multicellular organisms. Molecular Biology and Evolution, 1985, 2: 13-34.

[2] Angellotti M C, Bhuiyan S B, Chen G, Wan X F. CodonO: Codon usage bias analysis within and across genomes. Nucleic Acids Research, 2007, 35: 132-136.

[3] LuH,ZhaoW M,ZhengY, Wang H, Qi M, Yu X P.AnalysisofsynonymouscodonusagebiasinChlamydia. Acta Biochimica Biophysica Sinica (Shanghai), 2005,37(1):1-10.

[4] 郭秀丽, 王玉, 杨路成, 丁兆堂. 茶树CBF1基因密码子使用特性分析. 遗传, 2012, 34(12): 1614-1623.

Guo X L, Wang Y, Yang L C, Ding Z T. Analysis of codon use features of CBF gene in Camellia sinensis. Hereditas, 2012, 34(12): 1614-1623. (in Chinese)

[5] 石秀凡, 黄京飞, 柳树群, 刘次全. 人类基因同义密码子偏好的特征以及与基因GC含量的关系. 生物化学与生物物理进展, 2002, 29(3): 411-414.

Shi X F, Huang J F, Liu S Q, Liu C Q. The features of synonymous codon bias and GC-content relationship in human genes. Progress in Biochemistry and Biophysics, 2002, 29(3): 411-414. (in Chinese)

[6] Fedorov A, Saxonov S, Gilbert W. Regular ties of context- dependent codon bias in eukaryotic genes. Nucleic Acids Research, 2002, 30: 1192-1197.

[7] Hiraoka Y, Kawamata K, Haraguchi T, Chikashige Y. Codon usage bias is correlated with gene expression levels in the fission yeast Schizosaccharomyces pombe. Genes to Cells, 2009, 14: 499-509.

[8] Moriyama E, Powell J. Codon usage bias and tRNA abundance in Drosophila. Journal of Molecular Evolution, 1997, 45(5): 514-523.

[9] Oresic M, Dehn M, Korenblum D, Shalloway D. Tracing specific synonymous codon-secondary structure correlations through evolution. Journal of Molecular Evolution, 2003, 56: 473-484.

[10] Sakai H, Washio T, Saito R, Shinagawa A, Itoh M, Shibata K, Carninci P, Konno H, Kawai J, Hayashizaki Y, Tomita M. Correlation between sequENCe conservation of the 5′untranslated region and codon usage bias in Mus musculus genes. Gene, 2001, 276(1/2): 101-105.

[11] Gustafsson C, Govindarajan S, Minshull J. Codon bias and heterologous protein expression. Trends in Biotechnology, 2004, 22(7): 346-353.

[12] 吴正常, 王靖, 赵乔辉, 朱世平, 訾臣, 吴圣龙, 包文斌. 猪脂多糖结合蛋白基因(LBP)的密码子偏好性分析. 农业生物技术学报, 2013, 21(10): 1135-1144.

Wu Z C, Wang J, Zhao Q H, Zhu S P, Zi C, Wu S L, Bao W B. Analysis of codon bias of lipopolysaccharide-binding protein gene (LBP) in pigs (Sus scrofa). Journal of Agricultural Biotechnology, 2013, 21(10): 1135-1144. (in Chinese)

[13] Sharp P M, Bailes E, Grocock R J, Peden J F, Sockett R E. Variation in the strength of selected codon usage bias among bacteria. Nucleic Acids Research, 2005, 33: 1141-1153.

[14] Sharp P M, Tuohy T M, Mosurski K R. Codon usage in yeast: Cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Research, 1986, 14: 5125-5143.

[15] Zhou T, Sun X, Lu Z. Synonymous codon usage in environmental chlamydia UWE25 reflects an evolutional divergence from pathogenic chlamydiae. Gene, 2006, 368: 117-125.

[16] Duret L, Mouchiroud D. Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96: 4482-4487.

[17] Bulmer M. The selection-mutation-drift theory of synonymous codon usage. Genetics, 1991, 129(3): 897-907.

[18] 吴彦庆, 葛金涛, 陶俊. 芍药AP1(APETALA1)基因密码子使用的偏好性分析. 湖南农业大学学报(自然科学版), 2015, 41(6): 610-615.

Wu Y Q, Ge J T, Tao J. Analysis the codon bias of AP1 gene in Paeonia lactiflora. Journal of Hunan Agricultural University (Natural Sciences Edition), 2015, 41(6): 610-615. (in Chinese)

[19] Zhao D Q, Jiang Y, Ning C L, Meng J S, Lin S, Ding W, Tao J. Transcriptome sequencing of a chimaera reveals coordinated expression of anthocyanin biosynthetic genes mediating yellow formation in herbaceous peony (Paeonia lactiflora Pall.). BMC Genomics, 2014, 15(1): 689.

[20] Liu H M, He R, Zhang H Y, Huang Y B, Tian M L, Zhang J J. Analysis of synonymous codon usage in Zea mays. Molecular Biology Reports, 2010, 37(2): 677-684.

[21] Peden J F. CodonW. University of Nottingham, 1999.

[22] Sueoka N. Directional mutation pressure and neutral molecular evolution. Proceedings of the National Academy of Sciences, 1988, 85(8): 2653-2657.

[23] Wright F. The ‘effective number of codons’ used in a gene. Gene, 1990, 87: 23-29.

[24] Comeron J M, Aguade M. An evaluation of measures of synonymous codon usage bias. Journal of Molecular Evolution, 1998, 47: 268-274.

[25] Sharp P M, Li W H. The codon adaptation index-a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Research, 1987, 15: 1281-1295.

[26] SharpPM,LiWH.Anevolutionaryperspectiveonsynonymouscodonusageinunicellularorganisms. Journal of Molecular Evolution,1986, 24(1/2):28-38.

[27] Greenacre M J. Theory and applications of correspondence analysis. London: Academic Press 1984.

[28] Noboru S. Near homogeneity of PR2-bias fingerrints in the human genome and their implications in phylogenetic analyses. Journal of Molecular Evolution, 2001, 53: 469-476.

[29] Yang X, Luo X N, Cai X P. Analysis of codon usage pattern in Taenia saginata based on a transcriptome dataset. Parasites & Vectors, 2014, 7: 527.

[30] Zhou J H, Ding Y Z, He Y, Chu Y F, Zhao P, Zhao P, Ma L Y, Wang X J, Li X R, Liu Y S. The effect of multiple evolutionary selections on synonymous codon usage of genes in the Mycoplasma bovis Genome. Plos One, 2014, 9: e108949.

[31] Naya H, Romero H, Carels N, Zavala A, Musto H. Translational selection shapes codon usage in the GC-rich genome of Chlamydomonas reinhardtii. Febs Letters, 2001,501: 127-130.

[32] Gupta S, Bhattacharyya T, Ghosh T C. Synonymous codon usage in Lactococcus lactis: mutational bias versus translational selection. Journal of Biomolecular Structure & Dynamics, 2004, 21: 527-535.

[33] Bellgard M, Schibeci D, Trifonov E, Gojobori T. Early detection of G+C differences in bacterial species inferred from the comparative analysis of the two completely sequenced Helicobacter pylori strains. Journal of Molecular Evolution, 2001, 53(4/5): 465-468.

[34] Kawabe A, Miyashita N T. Patterns of codon usage bias in three dicot and four monocot plant species. Genes & Genetic Systems, 2003, 78: 343-352.

[35] Saul A, Battistutta D. Codon usage in Plasmodium falciparum. Molecular and Biochemical Parasitology, 1988, 27: 35-42.

[36] Muto A, Yamao F, Osawa S. The genome of Mycoplasma capricolum. Progress in Nucleic Acid Research and Molecular Biology, 1987, 34: 29-58.

[37] Dybvig K, Voelker L L. Molecular biology of mycoplasmas. Annual Reviews in Microbiology, 1996, 50(1): 25-57.

[38] 续晨, 蔡小宁, 钱保俐, 贲爱玲. 葡萄基因组密码子使用偏好模式研究. 西北植物学报, 2012, 32(2): 409-415.

Xu C, Cai X N, Qian B L, Ben A L. Codon usage bias in Vitis vinifera. Acta Botanica Boreali-Occidentalia Sinica, 2012, 32(2): 409-415. (in Chinese)

[39] 刘庆坡, 薛庆中. 粳稻叶绿体基因组的密码子用法. 作物学报, 2004, 30(12): 1220-1224.

Liu Q P, Xue Q Z. Codon Usage in the chloroplast genome of rice (Oryza sativa L.ssp. japonica). Acta Agronomica Sinica, 2004, 30(12): 1220-1224. (in Chinese)

[40] 续晨, 贲爱玲, 蔡晓宁. 蝴蝶兰叶绿体基因组密码子使用的相关分析. 分子植物育种, 2010, 8(5): 945-950.

Xu C, Ben A L, Cai X L. Analysis of synonymous codon usage in chloroplast genome of Phalaenopsis aphrodite subsp. Formosana. Molecular Plant Breeding, 2010, 8(5): 945-950. (in Chinese)

[41] 尚明照, 刘方, 华金平, 王坤波. 陆地棉叶绿体基因组密码子使用偏性的分析. 中国农业科学, 2011, 44(2): 245-253.

Shang M Z, Liu F, Hua J P, Wang K B. Analysis on codon usage of chloroplast genome of Gossypium hirsutum. Scientia Agricultura Sinica, 2011, 44(2): 245-253. (in Chinese)

[42] 李秀璋, 宋辉, 李春杰. 茄腐镰孢(Fusarium solani)线粒体基因组密码子偏好性分析. 基因组学与应用生物学, 2015, 34(11): 2465-2472.

Li X Z, Song H, Li C J. Analysis of codon usage bias in mitochondria genome of Fusarium solani. Genomics and Applied Biology, 2015, 34(11): 2465-2472. (in Chinese)

[43] Liu Q P, Feng Y, Zhao X A, Dong H, Xue Q Z. Synonymous codon usage bias in Oryza sativa. Plant Science, 2004, 167(1): 101-105.

[44] Zhang W J, Zhou J, Li Z F, Wang L, Gu X, Zhong Y. Comparative analysis of codon usage patterns among mitochondrion, chloroplast and nuclear genes in Triticum aestivum L. Journal of Integrative Plant Biology, 2007, 49(2): 246-254.

[45] Wang L J, Roossinck M J. Comparative analysis of expressed sequences reveals a conserved pattern of optimal codon usage in plants. Plant Molecular Biology, 2006, 61(4): 699-710.