流感病毒基因的密码子偏好性及聚类分析

流行性感冒病毒是一种造成人类及动物患流行性感冒的RNA病毒,它造成急性上呼吸道感染,并由空气迅速传播,在世界各地常有周期性的大流行。根据该病毒的基因组CDS序列,探讨了基因组序列密码子的使用模式和特性,并进行了病毒间的聚类分析。结果表明:流感病毒的G+C含量均低于A+U含量,偏向使用以A、U结尾的密码子的程度比使用以G、C结尾的较高,CUG、UCA、AGU、AGC、AGA、AGG、GUG、CCA、ACA、GGA、GCA、AUU、UGA、CAU、CAA、AAU、AAA、GAA等18个密码子为流感病毒共有的偏好性密码子,且以A结尾的居多,尤其偏爱AGA、GGA。聚类结果表明首先亚洲流感病毒H2N2和香港流感病毒H2N2聚为一类,亚洲流感病毒H1N1和俄罗斯流感病毒H1N1聚为一类,1997年和2003年~2004年发生的人禽流感聚为一类,说明它们的密码子使用的偏好性相似;而2009年爆发的甲型H1N1流感和任何一个流感的距离都比较远,说明甲型H1N1流感病毒是一种新型的病毒,不同于以往任何一种流感病毒。     

Analysis of the codon usage of the cephamycin C producer Nocardia lactamdurans

Abstract The G + C content in a sequenced region of 27 kb of the Nocardia lactamdurans genome is 70.4 and 70.6% in the 14 characterized ORFs, showing an extreme average G + C content (94.9%) in the third codon position. The codon usage parameters of the N. lactamdurans genes studied are closely related and depart weakly from the values of other species of the genus Nocardia . The homologies and differences in the codon usage between N. lactamdurans and Streptomyces sp. or other high-G + C Gram-positive genera are analysed.     

Cluster analysis of genes for nitrogen fixation from several diazotrophs

Hierarchical clustering and similarity coefficients of pairwise alignments of the published nucleotide sequences of 27nifH genes suggest thatnif genes are as ancient as the archaebacteria and clostridia. The positions ofnifHl ofMethanococcus thermolithotrophicus, nifH3 ofClostridium pasteurianum, nifH3 ofAzotobacter vinelandii andnifH ofFrankia suggest that a variety of lateral transfers may have occurred during evolution ofnifH gene. The genes for type 3 nitrogenase ofA. vinelandii may have diverged early from methanogens and clostridia. A high similarity coefficient with the derived amino acid sequence of type 3 nitrogenase suggests the presence of a functionally similar enzyme inC. pasteurianum. The type 2 nitrogenase genenifH2 of azotobacters seems to have originated recently from the genenifHl for conventional type I nitrogenase. RhizobialnifH genes comprise two closely related but discrete clusters that are in consonance with the plasmid or chromosomal location ofnif genes. The chromosomal and plasmid locatednifH of rhizobia seem to have evolved independently but contemporaneously.     

Nature of selective constraints on synonymous codon usage of rice differs in GC-poor and GC-rich genes

Synonymous codon usage and cellular tRNA abundance are thought to be co-evolved in optimizing translational efficiencies in highly expressed genes. Here in this communication by taking the advantage of publicly available gene expression data of rice and Arabidopsis we demonstrated that tRNA gene copy number is not the only driving force favoring translational selection in all highly expressed genes of rice. We found that forces favoring translational selection differ between GC-rich and GC-poor classes of genes. Supporting our results we also showed that, in highly expressed genes of GC-poor class there is a perfect correspondence between majority of preferred codons and tRNA gene copy number that confers translational efficiencies to this group of genes. However, tRNA gene copy number is not fully consistent with models of translational selection in GC-rich group of genes, where constraints on mRNA secondary structure play a role to optimize codon usage in highly expressed genes.     

An analysis of the codon usage of Pasteurella haemolytica A1

Analysis of approximately 17 kbp of nucleotide sequences from three different regions of the genome of Pasteurella haemolytica A1 showed that the mol% G+C of P. haemolytica A1 DNA is 38.5%. When only the coding sequences (approx. 10 kbp) were analysed, a similar value of 38.8% was obtained. A comparison of the relative synonymous codon usage values of the cloned genes showed that P. haemolytica A1 has a very different codon usage pattern from that of Escherichia coli.     

Nucleic acid composition,codon usage,and the rate of synonymous substitution in protein-coding genes

Summary Based on the rates of synonymous substitution in 42 protein-codin gene pairs from rat and human, a correlation is shown to exist between the frequency of the nucleotides in all positions of the codon and the synonymous substitution rate. The correlation coefficients were positive for A and T and negative for C and G. This means that AT-rich genes accumulate more synonymous substitutions than GC-rich genes. Biased patterns of mutation could not account for this phenomenon. Thus, the variation in synonymous substitution rates and the resulting unequal codon usage must be the consequence of selection against A and T in synonymous positions. Most of the varition in rates of synonymous substitution can be explained by the nucleotide composition in synonymous positions. Codon-anticodon interactions, dinucleotide frequencies, and contextual factors influence neither the rates of synonymous substitution nor codon usage. Interestingly, the nucleotide in the second position of codons (always a nonsynonymous position) was found to affect the rate of synonymous substitution. This finding links the rate of nonsynonymous substitution with the synonymous rate. Consequently, highly conservative proteins are expected to be encoded by genes that evolve slowly in terms of synonymous substitutions, and are consequently highly biased in their codon usage.     

Trends in the codon usage patterns of Chromohalobacter salexigens genes

Chromohalobacter salexigens, a Gammaproteobacterium belonging to the family Halomonadaceae, shows a broad salinity range for growth. In order to reveal the factors influencing architecture of protein coding genes in C. salexigens, pattern of synonymous codon usage bias has been investigated. Overall codon usage analysis of the microorganism revealed that C and G ending codons are predominantly used in all the genes which are indicative of mutational bias. Multivariate statistical analysis showed that the genes are separated along the first major explanatory axis according to their expression levels and their genomic GC content at the synonymous third positions of the codons. Both NC plot and correspondence analysis on Relative Synonymous Codon Usage (RSCU) indicates that the variation in codon usage among the genes may be due to mutational bias at the DNA level and natural selection acting at the level of mRNA translation. Gene length and the hydrophobicity of the encoded protein also influence the codon usage variation of genes to some extent. A comparison of the relative synonymous codon usage between 10% each of highly and lowly expressed genes determines 23 optimal codons, which are statistically over represented in the former group of genes and may provide useful information for salt-stressed gene prediction and gene-transformation. Furthermore, genes for regulatory functions; mobile and extrachromosomal element functions; and cell envelope are observed to be highly expressed. The study could provide insight into the gene expression response of halophilic bacteria and facilitate establishment of effective strategies to develop salt-tolerant crops of agronomic value.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻(Oryza rufipogon Griff.)线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻(O. sativa L.)的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示:普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数(Nc)分布于45.32～61.00之间,其密码子偏性较弱; Nc值仅与GC_3呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数(CBI)和最优密码子使用频率(Fop)的相关性均达到显著水平;而GC_3和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度(RSCU)、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻(O. sativa L. subsp. japonica Kato)和籼稻(O. sativa L. subsp.indica Kato)内具有同质性。     

Nonrandom CpG mutations affect the synonymous codon usage of moderately GC-rich single copy actin genes

Summary In species where actin genes exist as single copies, analysis of their synonymous codon usage and of the substitutions occurring between the genes of closely related species shows that there is a positive selection for codons that do not have highly mutable CpG dinucleotides in codon positions 2 and 3 when the GC content of these genes is less than 57%.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻（Oryza rufipogon Griff.）线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻（O.sativa L.）的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示：普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数（Nc）分布于45.32~61.00之间,其密码子偏性较弱;Nc值仅与GC₃呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数（CBI）和最优密码子使用频率（Fop）的相关性均达到显著水平;而GC₃和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度（RSCU）、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻（O.sativa L.subsp.japonica Kato）和籼稻（O.sativa L.subsp.indica Kato）内具有同质性。     

CpG islands are the second main factor shaping codon usage in human genes

A correspondence analysis of codon usage in human genes revealed, as expected, that the first axis is strongly correlated with the base composition at synonymous third codon positions. At one extreme of the second axis were localized genes with a high frequency of NCG and CGN codons. The great majority of these sequences were embedded in CpG islands, while the opposite is true for the genes placed at the other extreme. The two main conclusions of this paper are: (1) the influence of CpG islands on codon usage, and (2) since the second axis is orthogonal (and therefore independent) of the first, GC3-rich genes are not necessarily associated with CpG islands.     

Cluster analysis of genes in codon space

Summary We construct a codon space in which a given DNA sequence can be plotted as a function of its base composition in each of the three codon positions. We demonstrate that the base composition is very highly nonrandom, with sequences from more primitive organisms having the least random compositions. By using cluster analysis on the points plotted in codon space we show that there is a strong correlation between base composition and type of organism, with the most primitive organisms having the highest A or T content in the second and third codon positions. A smooth transition toward lower A+T and higher G+C content is observed in the second and third codon positions as the evolutionary complexity of the organism increases. Besides this general trend, more detailed structure can be observed in the clustering that will become clearer as the data base is increased.     

Analysis of synonymous codon usage bias in 09H1N1

A novel subtype of influenza A virus 09H1N1 has rapidly spread across the world. Evolutionary analyses of this virus have revealed that 09H1N1 is a triple reassortant of segments from swine, avian and human influenza viruses. In this study, we investigated factors shaping the codon usage bias of 09H1N1 and carried out cluster analysis of 60 strains of influenza A virus from different subtypes based on their codon usage bias. We discovered that more preferentially used codons of 09H1N1 are A-ended or U-ended, and the intra-genomic codon usage bias of 09H1N1 is quite low. Base composition constraint, dinucleotide biases and translational selection are the main factors influencing the codon usage bias of 09H1N1. At the genome level, we find that the codon usage bias of 09H1N1 is similar to H1N1 (A/swine/Kansas/77778/2007H1N1), H9N2 from Asia, H1N2 from Asia and North America and H3N2 from North America. Our results provide insight for understanding the processes governing evolution, regulation of gene expression, and revealing the evolution of 09H1N1.     

Synonymous codon usage in forty staphylococcal phages identifies the factors controlling codon usage variation and the phages suitable for phage therapy

The immergence and dissemination of multidrug-resistant strains of Staphylococcus aureus in recent years have expedited the research on the discovery of novel anti-staphylococcal agents promptly. Bacteriophages have long been showing tremendous potentialities in curing the infections caused by various pathogenic bacteria including S. aureus. Thus far, only a few virulent bacteriophages, which do not carry any toxin-encoding gene but are capable of eradicating staphylococcal infections, were reported. Based on the codon usage analysis of sixteen S. aureus phages, previously three phages were suggested to be useful as the anti-staphylococcal agents. To search for additional S. aureus phages suitable for phage therapy, relative synonymous codon usage bias has been investigated in the protein-coding genes of forty new staphylococcal phages. All phages appeared to carry A and T ending codons. Several factors such as mutational pressure, translational selection and gene length seemed to be responsible for the codon usage variation in the phages. Codon usage indeed varied phage to phage. Of the phages, phages G1, Twort, 66 and Sap-2 may be extremely lytic in nature as majority of their genes possess high translational efficiency, indicating that these phages may be employed in curing staphylococcal infections.     

Analysis of codon usage pattern in the radioresistant bacterium Deinococcus radiodurans

The main factors shaping codon usage bias in the Deinococcus radiodurans genome were reported. Correspondence analysis (COA) was carried out to analyze synonymous codon usage bias. The results showed that the main trend was strongly correlated with gene expression level assessed by the "Codon Adaptation Index" (CAI) values, a result that was confirmed by the distribution of genes along the first axis. The results of correlation analysis, variance analysis and neutrality plot indicated that gene nucleotide composition was clearly contributed to codon bias. CDS length was also key factor in dictating codon usage variation. A general tendency of more biased codon usage of genes with longer CDS length to higher expression level was found. Further, the hydrophobicity of each protein also played a role in shaping codon usage in this organism, which could be confirmed by the significant correlation between the positions of genes placed on the first axis and the hydrophobicity values (r=-0.100, P<0.01). In summary, gene expression level played a crucial role, nucleotide mutational bias, CDS length and the hydrophobicity of each protein just in a minor way in shaping the codon usage pattern of D. radiodurans. Notably, 19 codons firstly defined as "optimal codons" may provide useful clues for molecular genetic engineering and evolutionary studying.     

Non-uniqueness of factors constraint on the codon usage in Bombyx mori

Background

The analysis of codon usage is a good way to understand the genetic and evolutionary characteristics of an organism. However, there are only a few reports related with the codon usage of the domesticated silkworm, Bombyx mori (B. mori). Hence, the codon usage of B. mori was analyzed here to reveal the constraint factors and it could be helpful to improve the bioreactor based on B. mori.

Results

A total of 1,097 annotated mRNA sequences from B. mori were analyzed, revealing there is only a weak codon bias. It also shows that the gene expression level is related to the GC content, and the amino acids with higher general average hydropathicity (GRAVY) and aromaticity (Aromo). And the genes on the primary axis are strongly positively correlated with the GC content, and GC3s. Meanwhile, the effective number of codons (ENc) is strongly correlated with codon adaptation index (CAI), gene length, and Aromo values. However, the ENc values are correlated with the second axis, which indicates that the codon usage in B. mori is affected by not only mutation pressure and natural selection, but also nucleotide composition and the gene expression level. It is also associated with Aromo values, and gene length. Additionally, B. mori has a greater relative discrepancy in codon preferences with Drosophila melanogaster (D. melanogaster) or Saccharomyces cerevisiae (S. cerevisiae) than with Arabidopsis thaliana (A. thaliana), Escherichia coli (E. coli), or Caenorhabditis elegans (C. elegans).

Conclusions

The codon usage bias in B. mori is relatively weak, and many influence factors are found here, such as nucleotide composition, mutation pressure, natural selection, and expression level. Additionally, it is also associated with Aromo values, and gene length. Among them, natural selection might play a major role. Moreover, the “optimal codons” of B. mori are all encoded by G and C, which provides useful information for enhancing the gene expression in B. mori through codon optimization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1596-z) contains supplementary material, which is available to authorized users.     

Nucleotide composition and synonymous codon usage of open reading frames in Norovirus GII.4 variants

Abstract

Norovirus GII.4 variants, a genotype in genogroup II belonging to the genus Norovirus, is a single-strand positive sense RNA containing three open reading frames (ORF1, ORF2 and ORF3) and is the most important pathogen causing nonbacterial gastroenteritis outbreaks. By using bioinformatic softwares such as Codon W, SPSS and so on, a total of 292 strains of the viruses isolated from 1974 to 2016 were analyzed for nucleotide composition and synonymous codon usage in each ORF. The result shows that it is enriched for A over the other bases in nucleotide composition, G behind the other bases in the 3^rd site of all synonymous codons in the three ORFs. The patterns of nucleotide composition and codon bias of ORF2 are similar to those of ORF3 and different from those of ORF1. There are generally UpA motif and CpG motif in the codons with the lowest proportion. Correspondence analysis indicates that the codon usage may be changing over a certain time period for ORF1 in 2006 and 2012, ORF2 in 2012, and ORF3 in 2013. ENC (effective number of codons) plot and other analyses indicate that both natural selection and mutational pressure play partly roles in the ORFs, but natural selection is more important for ORF2 and ORF3. Besides, we also found all optimal codons in the ORFs. The study provides a basic understanding of the mechanism for norovirus GII.4 codon usage bias. Abbreviations ORF Open Reading Frame

ENC Effective Number of Codons

COA correspondence analysis

RSCU Relative Synonymous Codon Usage

CAI Codon Adaptation Index

CBI Codon Bias Index

Fop frequency of optimal codons

L_sym number of synonymous codons

L_aa length amino acids

GRAVY grand average of hydropathicity

Aroma aromaticity

Communicated by Ramaswamy H. Sarma