亚氨基酸编码方法及其应用

赋予氨基酸编码方法下除终止子之外的密码子突变为终止子时每一位发生的变化权值,利用矩阵来表示所有的突变方式和难易程度,综合亲水性与疏水性理化性质,提出亚氨基酸编码方法,并给出该编码方法下同义密码子的相对使用度fSubtypesRelativeSynonymousCodonUsage,SRSCU）．然后选取15条H5N1序列,使用MEGA4．0分析它们的同源性,并分别在氨基酸编码、拟氨基酸编码、亚氨基酸编码这三种环境下研究所选序列使用密码子的偏好性,对比结果验证,亚氨基酸编码方法具有相应的优越性．     

沙枣叶绿体全基因组序列及其使用密码子偏性分析

该研究以2株野生沙枣(Elaeagnus angustifolia Linn.)嫩枝经温室水培后的嫩叶为材料,采用CTAB法分别提取总DNA,并利用第二代测序技术进行总DNA从头测序,组装后得到2株沙枣叶绿体基因组全序列,并详细分析了其蛋白质编码基因密码子使用的偏好性及其原因,为沙枣叶绿体基因工程和分子系统进化等研究奠定基础。结果显示:(1)组装得到沙枣叶绿体基因组序列全长150 546 bp,由长度为81 113 bp的长单拷贝(LSC)区域和25 494 bp的短单拷贝(SSC)区域,以及1对分隔开它们的长18 445 bp的反向重复序列(IRS)组成;注释共得到132个基因,包括86个蛋白编码基因、38个tRNA基因和8个rRNA基因。(2)沙枣叶绿体基因组蛋白编码基因密码子的第三位碱基GC含量(GC_3)为28.47%,明显低于整个叶绿体基因组GC含量(37%),也低于第一位(GC_1)和第二位(GC_2)碱基的GC含量,说明密码子对AT碱基结尾有偏好性;其中, UCU、CCU、UGU、GCU、CUU、GAU、UCA和UAA为最优密码子。(3)同义密码子相对使用频率(RSCU)分析发现,影响密码子使用模式的因素并不单一,密码子的偏好性受到突变、选择及其他因素的共同影响,并且自然选择表达引起的序列差异比突变对密码子偏好性的影响要显著;中性绘图分析、有效密码子数(ENC-plot)分析和奇偶偏好性(PR2-plot)分析表明,沙枣叶绿体基因组使用密码子的偏性受选择的影响更大。(4)通过最大似然法、最大简约法和贝叶斯方法对胡颓子科6个物种和1个枣的叶绿体基因序列构建系统发育树,与它们使用密码子偏性聚类的结果一致,表明叶绿体基因组使用密码子偏性与物种的亲缘关系相关。     

杜仲基因密码子使用模式分析

为了解杜仲基因密码子使用模式,该文以杜仲基因组密码子为研究对象,运用CodonW软件对杜仲的320个蛋白编码基因进行同义密码子相对使用频率(RSCU)分析、ENC-GC3s关联分析编码基因的密码子ENC值、PR2-plot偏倚分析编码基因的密码子碱基使用频率,并运用CUSP软件与Codon Usage Database...     

Synonymous Codon Usage in TTSuV2: Analysis and Comparison with TTSuV1

Two species of the DNA virus Torque teno sus virus (TTSuV), TTSuV1 and TTSuV2, have become widely distributed in pig-farming countries in recent years. In this study, we performed a comprehensive analysis of synonymous codon usage bias in 41 available TTSuV2 coding sequences (CDS), and compared the codon usage patterns of TTSuV2 and TTSuV1. TTSuV codon usage patterns were found to be phylogenetically conserved. Values for the effective number of codons (ENC) indicated that the overall extent of codon usage bias in both TTSuV2 and TTSuV1 was not significant, the most frequently occurring codons had an A or C at the third codon position. Correspondence analysis (COA) was performed and TTSuV2 and TTSuV1 sequences were located in different quadrants of the first two major axes. A plot of the ENC revealed that compositional constraint was the major factor determining the codon usage bias for TTSuV2. In addition, hierarchical cluster analysis of 41 TTSuV2 isolates based on relative synonymous codon usage (RSCU) values suggested that there was no association between geographic distribution and codon bias of TTSuV2 sequences. Finally, the comparison of RSCU for TTSuV2, TTSuV1 and the corresponding host sequence indicated that the codon usage pattern of TTSuV2 was similar to that of TTSuV1. However the similarity was low for each virus and its host. These conclusions provide important insight into the synonymous codon usage pattern of TTSuV2, as well as better understangding of the molecular evolution of TTSuV2 genomes.     

Gene optimization mechanisms: A multi-gene study reveals a high success rate of full-length human proteins expressed in Escherichia coli

The genetic code is universal, but recombinant protein expression in heterologous systems is often hampered by divergent codon usage. Here, we demonstrate that reprogramming by standardized multi‐parameter gene optimization software and de novo gene synthesis is a suitable general strategy to improve heterologous protein expression. This study compares expression levels of 94 full‐length human wt and sequence‐optimized genes coding for pharmaceutically important proteins such as kinases and membrane proteins in E. coli. Fluorescence‐based quantification revealed increased protein yields for 70% of in vivo expressed optimized genes compared to the wt DNA sequences and also resulted in increased amounts of protein that can be purified. The improvement in transgene expression correlated with higher mRNA levels in our analyzed examples. In all cases tested, expression levels using wt genes in tRNA‐supplemented bacterial strains were outperformed by optimized genes expressed in non‐supplemented host cells.     

The compositional transition of vertebrate genomes: an analysis of the secondary structure of the proteins encoded by human genes

Fluctuations and increments of both C(3) and G(3) levels along the human coding sequences were investigated comparing two sets of Xenopus/human orthologous genes. The first set of genes shows minor differences of the GC(3) levels, the second shows considerable increments of the GC(3) levels in the human genes. In both data sets, the fluctuations of C(3) and G(3) levels along the coding sequences correlated with the secondary structures of the encoded proteins. The human genes that underwent the compositional transition showed a different increment of the C(3) and G(3) levels within and among the structural units of the proteins. The relative synonymous codon usage (RSCU) of several amino acids were also affected during the compositional transition, showing that there exists a correlation between RSCU and protein secondary structures in human genes. The importance of natural selection for the formation of isochore organization of the human genome has been discussed on the basis of these results.     

Detecting overlapping coding sequences in virus genomes

Background  

Detecting new coding sequences (CDSs) in viral genomes can be difficult for several reasons. The typically compact genomes often contain a number of overlapping coding and non-coding functional elements, which can result in unusual patterns of codon usage; conservation between related sequences can be difficult to interpret – especially within overlapping genes; and viruses often employ non-canonical translational mechanisms – e.g. frameshifting, stop codon read-through, leaky-scanning and internal ribosome entry sites – which can conceal potentially coding open reading frames (ORFs).     

The Synthetic Gene Designer: a flexible web platform to explore sequence manipulation for heterologous expression

"Codon optimization" is a general approach to improving heterologous expression where genes are moved from their native genomes into alternatives that exhibit different patterns of codon usage. However, despite reports of successful manipulations and the existence of stand-alone codon optimization software packages or commercial services that offer to redesign genes, the scientific community lacks any systematic understanding of what exactly it means to optimize codon usage. Thus we present a bona fide web application, the "Synthetic Gene Designer," which contrasts with existing software by providing a centralized, free, and transparent platform for the broader scientific community to develop knowledge about synthetic gene design. Consistent with this goal, our software is associated with a moderated e-forum that promotes discussion of synthetic gene design and offers technical support. In addition, the Synthetic Gene Designer presents enhanced functionality over existing software options: for example, it enables users to work with non-standard genetic codes, with user-defined patterns of codon usage and an expanded range of methods for codon optimization. The Synthetic Gene Designer, together with on-line tutorials and the forum, is available at .     

Codon usage bias in human cytomegalovirus and its biological implication

Human cytomegalovirus (HCMV) infection, a worldwide contagion, causes a serious disorder in infected individuals. Analysis of codon usage can reveal much molecular information about this virus. The effective number of codon (ENC) values, relative synonymous codon usage (RSCU) values, codon adaptation index (CAI), and nucleotide contents was investigated in approximately 160 coding sequences (CDS) among 17 human cytomegalovirus genomes using the software CodonW. Linear regression analysis and logistic regression were performed to explore the preliminary data. The results showed that, overall, HCMV genomes had low codon usage bias (mean ENC = 47.619). However, the ENC of individual CDS varied widely and was distributed unevenly between host-related genes and viral-self-function genes (P = 0.002, odds ratio (OR) = 3.194), as did the GC content (P = 0.016, OR = 2.178). The ENC values correlated with CAI, GC content, and the nucleotide composing at the 3rd codon position (GC3s) (P < 0.001). There was a significant variation in the codon preference that depended on the RSCU data. The predicted ENC curve suggested that mutational pressure, rather than natural selection, was one of the main factors that determined the codon usage bias in HCMV. Among 123 genes with known function, the genes related to viral self-replication and viral–host interaction showed different ENC and CAI values, and GC and GC3s contents. In conclusion, the detailed codon usage bias theoretically revealed information concerning HCMV evolution and could be a valuable additional parameter for HCMV gene function research.     

Analysis of synonymous codon usage in 11 human bocavirus isolates

Human Bocavirus (HBoV) is a novel virus which can cause respiratory tract disease in infants or children. In this study, the codon usage bias and the base composition variations in the available 11 complete HBoV genome sequences have been investigated. Although, there is a significant variation in codon usage bias among different HBoV genes, codon usage bias in HBoV is a little slight, which is mainly determined by the base compositions on the third codon position and the effective number of codons (ENC) value. The results of correspondence analysis (COA) and Spearman's rank correlation analysis reveals that the G + C compositional constraint is the main factor that determines the codon usage bias in HBoV and the gene's function also contributes to the codon usage in this virus. Moreover, it was found that the hydrophobicity of each protein and the gene length are also critical in affecting these viruses’ codon usage, although they were less important than that of the mutational bias and the genes’ function. At last, the relative synonymous codon usage (RSCU) of 44 genes from these 11 HBoV isolates is analyzed using a hierarchical cluster method. The result suggests that genes with same function yet from different isolates are classified into the same lineage and it does not depend on geographical location. These conclusions not only can offer an insight into the codon usage patterns and gene classification of HBoV, but also may help in increasing the efficiency of gene delivery/expression systems.     

Comparative analysis of human coronaviruses focusing on nucleotide variability and synonymous codon usage patterns

The prevailing COVID-19 pandemic has drawn the attention of the scientific community to study the evolutionary origin of Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2). This study is a comprehensive quantitative analysis of the protein-coding sequences of seven human coronaviruses (HCoVs) to decipher the nucleotide sequence variability and codon usage patterns. It is essential to understand the survival ability of the viruses, their adaptation to hosts, and their evolution.The current analysis revealed a high abundance of the relative dinucleotide (odds ratio), GC and CT pairs in the first and last two codon positions, respectively, as well as a low abundance of the CG pair in the last two positions of the codon, which might be related to the evolution of the viruses. A remarkable level of variability of GC content in the third position of the codon among the seven coronaviruses was observed. Codons with high RSCU values are primarily from the aliphatic and hydroxyl amino acid groups, and codons with low RSCU values belong to the aliphatic, cyclic, positively charged, and sulfur-containing amino acid groups. In order to elucidate the evolutionary processes of the seven coronaviruses, a phylogenetic tree (dendrogram) was constructed based on the RSCU scores of the codons. The severe and mild categories CoVs were positioned in different clades. A comparative phylogenetic study with other coronaviruses depicted that SARS-CoV-2 is close to the CoV isolated from pangolins (Manis javanica, Pangolin-CoV) and cats (Felis catus, SARS(r)-CoV). Further analysis of the effective number of codon (ENC) usage bias showed a relatively higher bias for SARS-CoV and MERS-CoV compared to SARS-CoV-2. The ENC plot against GC3 suggested that the mutational bias might have a role in determining the codon usage variation among candidate viruses.A codon adaptability study on a few human host parasites (from different kingdoms), including CoVs, showed a diverse adaptability pattern. SARS-CoV-2 and SARS-CoV exhibit relatively lower but similar codon adaptability compared to MERS-CoV.     

甜瓜蔗糖转化酶基因的密码子偏好性分析

以甜瓜蔗糖转化酶基因序列为材料,研究甜瓜蔗糖转化酶基因密码子偏好性,为改良甜瓜风味与品质提供理论依据。运用在线分析软件Codon W对甜瓜蔗糖转化酶的编码序列(Coding sequence,CDS)进行密码子分析,利用Mobyle在线工具分析同义密码子相对使用度(RSCU)、有效密码子数(ENC)、GC及GC1s、GC2s、GC3s含量。甜瓜蔗糖转化酶基因偏好于以A或T结尾的密码子。密码子ATT、GTT和AGA的RSCU值都大于1,属于共同偏好使用的密码子,而密码子GCG、CGG的RSCU值小于1,属于使用频率较低的密码子。发现密码子偏好性与亲缘关系的远近有一定的关系。要实现目的基因在外源表达系统中的成功表达和提高其表达量,可通过增加目的基因剂量,目的基因密码子优化,改善培养条件等方法实现,其中目的基因密码子优化起到了关键作用。     

Optimization of codon usage of the envelope protein E2 gene from various genotypes of hepatitis C virus to predict the expression level in <Emphasis Type="Italic">Pichia pastoris</Emphasis>

Hepatitis C virus infection (HCV) alarmingly increases worldwide; it causes chronic hepatitis, liver cirrhosis and hepatocellular carcinoma, so there is urgent need of developing effective and sufficient quantity of vaccine. HCV envelope protein E2 is the main target for developing as a vaccine candidate. Presently recombinant proteins can successfully be used as a vaccine for many diseases. This concern, it is challenging to produce sufficient quantities of many recombinant proteins from their expression hosts. One of the main factors affecting the success of expression of foreign genes in heterologous hosts is the divergence of codon usage of the target gene from that used in the expression system. In this study, we optimized the various genotypes of HCV envelope protein E2 gene according to the codon usage of Pichia pastoris and predicted the expression level. Synonymous codon usage of E2 adapted to that used by P. pastoris was estimated using the relative synonymous codon usage value (RSCU), codon adaptation index (CAI) and effective number of codon (ENC). The CAI of optimized HCV E2 sequences was enhanced from 0.638 to 0.833 and %GC was decreased from 56.05 to 44.05; this was significantly (p < 0.01) different from the native sequences. Codon with RSCU value less than one was replaced with most preferred synonymous codons. The ENC values of optimized HCV E2 sequences varied from 47.00 to 47.50, with a mean value of 47.15 and an SD of 0.14. Our study suggested that, from the measured values of predicted expression level, the codon optimized HCV E2 protein could be produced in sufficient quantity in the expression host; knowledge of the codon usage patterns of E2 of various genotypes facilitate the production of a promising unique vaccine candidate for HCV.     

Insight into pattern of codon biasness and nucleotide base usage in serotonin receptor gene family from different mammalian species

5-HT (5-Hydroxy-tryptamine) or serotonin receptors are found both in central and peripheral nervous system as well as in non-neuronal tissues. In the animal and human nervous system, serotonin produces various functional effects through a variety of membrane bound receptors. In this study, we focus on 5-HT receptor family from different mammals and examined the factors that account for codon and nucleotide usage variation. A total of 110 homologous coding sequences from 11 different mammalian species were analyzed using relative synonymous codon usage (RSCU), correspondence analysis (COA) and hierarchical cluster analysis together with nucleotide base usage frequency of chemically similar amino acid codons. The mean effective number of codon (ENc) value of 37.06 for 5-HT(6) shows very high codon bias within the family and may be due to high selective translational efficiency. The COA and Spearman's rank correlation reveals that the nucleotide compositional mutation bias as the major factors influencing the codon usage in serotonin receptor genes. The hierarchical cluster analysis suggests that gene function is another dominant factor that affects the codon usage bias, while species is a minor factor. Nucleotide base usage was reported using Goldman, Engelman, Stietz (GES) scale reveals the presence of high uracil (>45%) content at functionally important hydrophobic regions. Our in silico approach will certainly help for further investigations on critical inference on evolution, structure, function and gene expression aspects of 5-HT receptors family which are potential antipsychotic drug targets.     

Comparison of Correspondence Analysis Methods for Synonymous Codon Usage in Bacteria

Synonymous codon usage varies both between organisms and among genes within a genome, and arises due to differences in G + C content, replication strand skew, or gene expression levels. Correspondence analysis (CA) is widely used to identify major sources of variation in synonymous codon usage among genes and provides a way to identify horizontally transferred or highly expressed genes. Four methods of CA have been developed based on three kinds of input data: absolute codon frequency, relative codon frequency, and relative synonymous codon usage (RSCU) as well as within-group CA (WCA). Although different CA methods have been used in the past, no comprehensive comparative study has been performed to evaluate their effectiveness. Here, the four CA methods were evaluated by applying them to 241 bacterial genome sequences. The results indicate that WCA is more effective than the other three methods in generating axes that reflect variations in synonymous codon usage. Furthermore, WCA reveals sources that were previously unnoticed in some genomes; e.g. synonymous codon usage related to replication strand skew was detected in Rickettsia prowazekii. Though CA based on RSCU is widely used, our evaluation indicates that this method does not perform as well as WCA.Key words: correspondence analysis, synonymous codon usage, horizontal gene transfer, strand-specific mutational bias, translational selection     

Nucleotide sequence from the genetic left end of bacteriophage T7 DNA to the beginning of gene 4

The nucleotide sequence running from the genetic left end of bacteriophage T7 DNA to within the coding sequence of gene 4 is given, except for the internal coding sequence for the gene 1 protein, which has been determined elsewhere. The sequence presented contains nucleotides 1 to 3342 and 5654 to 12,100 of the approximately 40,000 base-pairs of T7 DNA. This sequence includes: the three strong early promoters and the termination site for Escherichia coli RNA polymerase: eight promoter sites for T7 RNA polymerase; six RNAase III cleavage sites; the primary origin of replication of T7 DNA; the complete coding sequences for 13 previously known T7 proteins, including the anti-restriction protein, protein kinase, DNA ligase, the gene 2 inhibitor of E. coli RNA polymerase, single-strand DNA binding protein, the gene 3 endonuclease, and lysozyme (which is actually an N-acetylmuramyl-l-alanine amidase); the complete coding sequences for eight potential new T7-coded proteins; and two apparently independent initiation sites that produce overlapping polypeptide chains of gene 4 primase. More than 86% of the first 12,100 base-pairs of T7 DNA appear to be devoted to specifying amino acid sequences for T7 proteins, and the arrangement of coding sequences and other genetic elements is very efficient. There is little overlap between coding sequences for different proteins, but junctions between adjacent coding sequences are typically close, the termination codon for one protein often overlapping the initiation codon for the next. For almost half of the potential T7 proteins, the sequence in the messenger RNA that can interact with 16 S ribosomal RNA in initiation of protein synthesis is part of the coding sequence for the preceding protein. The longest non-coding region, about 900 base-pairs, is at the left end of the DNA. The right half of this region contains the strong early promoters for E. coli RNA polymerase and the first RNAase III cleavage site. The left end contains the terminal repetition (nucleotides 1 to 160), followed by a striking array of repeated sequences (nucleotides 175 to 340) that might have some role in packaging the DNA into phage particles, and an A · T-rich region (nucleotides 356 to 492) that contains a promoter for T7 RNA polymerase, and which might function as a replication origin.     

Genome-wide comparative analysis of the codon usage patterns in plants

Codon usage analysis has been a classical area of study for decades and is important for evolution, mRNA translation, and new gene discovery. Recently, genome sequencing has made it possible to perform studies of the entire genome in plant kingdoms. The base composition of the coding sequence, codon usage pattern, codon pairs, and related indicators of relative synonymous codon usage (RSCU), including the Fop, Nc, RSCU, CAI and GC contents, were analyzed. We found that the GC content of single-celled algae is the highest, whereas dicotyledons are the lowest. Moreover, the base composition of plants is similar within the same family. In addition, the GC content of the second base of the codon is lower than the first and third base. In conclusion, the codon usage characteristics are opposite in Gramineae, single-celled algae, fern and dicotyledon, moss, and Pinaceae. Furthermore, the degree of codon usage bias is decreasing with evolution. Therefore, we hypothesize that the lower the plants, the more that they must optimize codons and that higher plants no longer need to optimize codons.     

Analysis of synonymous codon usage patterns of HPRT1 gene across twelve mammalian species

Genetic changes in Hypoxanthine guanine phosphoribosyltransferace (HPRT1) gene can alter the expression of the dopamine neurotransmitter leads to abnormal neuron function, a disease called Lesch-Nyhan syndrome (LNS). Although different studies were conducted on LNS, information on codon usage bias (CUB) of HPRT1 gene is limited. The present study examines the genetic determinants of CUB in HPRT1 gene using twelve mammalian species. In the coding sequence of HPRT1 genes, A/T ending codons was most frequently used. A higher ENC value was observed indicating lower HPRT1 gene expression in the selected mammalian species. Correlation analysis indicates that compositional constraints under mutation pressure can involve in CUB of HPRT1 genes among the selected mammalian species. Relative synonymous codon usage (RSCU) value revealed that the codons such as ACT, AGG, ATT and AGC were over-represented in each of the mammalian species. Result from the analysis of the RSCU indicates that compositional constraint is a key driver for the variation in codon usage. Ratio of nonsynonymous (dN) and synonymous (dS) substitution further suggested that purifying selection occurs among the HPRT1 gene of studied mammals to maintain its protein function under the process of evolution. Our findings report an insight into the codon usage patterns of HPRT1 gene and will be useful for LNS management.