Ribosome traffic in E. coli and regulation of gene expression

Analysis of codon usage bias of WRKY transcription factors in Helianthus annuus

Background

The phenomenon of codon usage bias is known to exist in many genomes and is mainly determined by mutation and selection. Codon usage bias analysis is a suitable strategy for identifying the principal evolutionary driving forces in different organisms. Sunflower (Helianthus annuus L.) is an annual crop that is cultivated worldwide as ornamentals, food plants and for their valuable oil. The WRKY family genes in plants play a central role in diverse regulation and multiple stress responses. Evolutionary analysis of WRKY family genes of H. annuus can provide rich genetic information for developing hybridization resources of the genus Helianthus.

Results

Bases composition analysis showed the average GC content of WRKY genes of H. annuus was 43.42%, and the average GC3 content was 39.60%, suggesting that WRKY gene family prefers A/T(U) ending codons. There were 29 codons with relative synonymous codon usage (RSCU) greater than 1 and 22 codons ending with A and U base. The effective number of codons (ENC) and codon adaptation index (CAI) in WRKY genes ranged from 43.47–61.00 and 0.14–0.26, suggesting that the codon bias was weak and WRKY genes expression level was low. Neutrality analysis found a significant correlation between GC12 and GC3. ENC-plot showed most genes on or close to the expected curve, suggesting that mutational bias played a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that the usage of AT and GC was disproportionate. A total of three codons were identified as the optimal codons.

Conclusion

Apart from natural selection effects, most of the genetic evolution in the H. annuus WRKY genome might be driven by mutation pressure. Our results provide a theoretical foundation for elaborating the genetic architecture and mechanisms of H. annuus and contributing to enrich H. annuus genetic resources.

Analysis of codon usage patterns and influencing factors in rice tungro bacilliform virus

Rice tungro bacilliform virus (RTBV) belongs to genus Tungrovirus within the family Caulimoviridae harbors circular double-stranded DNA (dsDNA). Rice tungro disease (RTD) caused by RTBV, responsible for severe rice yield losses in South and Southeast Asia. Here, we performed a systematic evolutionary and codon usage bias (CUB) analysis of RTBV genome sequences. We analysed different bioinformatics techniques to calculate the nucleotide compositions, the relative synonymous codon usage (RSCU), and other indices. The results indicated slightly or low codon usage bias in RTBV isolates. Mutation and natural selection pressures have equally contributed to this low codon usage bias. Additionally, multiple factors such as host, geographical distribution also affect codon usage patterns in RTBV genomes. RSCU analysis revealed that RTBV shows mutation bias and prefers A and U ended codons to code amino acids. Codon usage patterns of RTBV were also found to be influenced by its host. This indicates that RTBV have evolved codon usage patterns that are specific to its host. The findings from this study are expected to increase our understanding of factors leading to viral evolution and fitness with respect to hosts and the environment.

Composition, codon usage pattern, protein properties, and influencing factors in the genomes of members of the family Anelloviridae

The present study was carried out on 62 genome sequences of members of the family Anelloviridae, as there have been no reports of genome analysis of these DNA viruses using a bioinformatics approach. The genes were found to be rich in AC content with low codon usage bias (CUB). Relative synonymous codon usage (RSCU) values identified the preferred codons for each amino acid in the family. The codon AGA was overrepresented, while the codons TCG, TTG, CGG, CGT, ACG, GCG and GAT were underrepresented in all of the genomes. A significant correlation was found between the effective number of codons (ENC) and base constraints, indicating that compositional properties might have influenced the CUB. A highly significant correlation was observed between the overall base content and the base content at the third codon position, indicating that mutations might have affected the CUB. A highly significant positive correlation was observed between GC12 and GC3 (r = 0.904, p < 0.01), which indicated that directional mutation pressure influenced all three codon positions. A neutrality plot revealed that the contribution of mutation and natural selection in determining the CUB was 58.6% and 41.4%, respectively.

Analysis of Synonymous Codon Usage Bias in Flaviviridae Virus

Background

Flaviviridae viruses are single-stranded, positive-sense RNA viruses, which threat human constantly mediated by mosquitoes, ticks, and sandflies. Considering the recent increase in the prevalence of the family virus and its risk potential, we investigated the codon usage pattern to understand its evolutionary processes and provide some useful data to develop the medications for most of Flaviviridae viruses.

Results

The overall extent of codon usage bias in 65 Flaviviridae viruses is low with the average value of GC contents being 50.5% and the highest value being 55.9%; the lowest value is 40.2%. ENC values of Flaviviridae virus genes vary from 48.75 to 57.83 with a mean value of 55.56. U- and A-ended codons are preferred in the Flaviviridae virus. Correlation analysis shows that the positive correlation between ENC value and GC content at the third nucleotide positions was significant in this family virus. The result of analysis of ENC, neutrality plot analysis, and correlation analysis revealed that codon usage bias of all the viruses was affected mainly by natural selection. Meanwhile, according to correspondence analysis (CoA) based on RSCU and phylogenetic analysis, the Flaviviridae viruses mainly are made up of two groups, Group I (Yellow fever virus, Apoi virus, Tembusu virus, Dengue virus 1, and others) and Group II (West Nile virus lineage 2, Japanese encephalitis virus, Usutu virus, Kedougou virus, and others).

Conclusions

All in, the bias of codon usage pattern is affected not only by compositional constraints but also by natural selection. Phylogenetic analysis also illustrates that codon usage bias of virus can serve as an effective means of evolutionary classification in Flaviviridae virus.

Analysis of the codon usage pattern in Middle East Respiratory Syndrome Coronavirus

Middle East Respiratory Syndrome Coronavirus (MERS-CoV), which first broken out in Jeddah in 2012, causes a severe acute respiratory illness with a high mortality rate. To better understand the molecular characteristics of isolated MERS-CoV genomes, we first analysed the codon usage pattern of the zoonotic MERS-CoV strains comprehensively to gain an insight into the mechanism of cross-species transmission. We found that MERS human/camel isolates showed a low codon usage bias. Both mutation and nature selection pressure have contributed to this low codon usage bias, with the former being the main determining factor. We also observed that gene function, evolution time and the different host species of the virus all contributed to the bias of MERS-CoV, to some extent. Additionally, the codon usage pattern of MERS-CoV isolates is different from other related Nidovirales viruses isolated from bats and hedgehogs. In the future, more epidemiological surveys are required to examine the factors that resulted in the emergence and outbreak of this virus.

Analysis of codon usage pattern evolution in avian rotaviruses and their preferred host

Rotavirus infection is a worldwide problem, with occurrence of highly divergent viruses classified in 8 species (A-H). We report here the evolution assessment of codon usage patterns in virus-host system in avian rotavirus (AvRV) of species RVA, RVD, RVF and RVG (preferentially affecting birds). The nucleotide contents, codon usage bias (CUB), relative synonymous codon usage (RSCU), and effective number of codons (ENCs) values were investigated targeting overexpressing major inner capsid viral protein (VP6) of these AvRV species. The results confirm that the evolutionary characteristics influences the rotavirus (RV) genetic diversity and impact of host's natural selection on the AvRVs codons. Synonymous codon usage patterns were evaluated following multivariate statistical procedures on all available AvRV coding gene sequences. RSCU trees accommodated all AvRV species and preferred host sequences in one topology confirming greater imminence of AvRVs with the host chicken cell genes. Similarly, the codon adaptation index (CAI) results also displayed a higher adaptation of AvRVs to its chicken host. The codon preference analysis of RVs revealed that VP6 gene express more proficiently in the yeast system, whereas, codon optimization might be required for the effectual expression in Escherichia coli and Homo sapiens. The findings provide basic evidence on the dynamics of AvRV evolution and its host adaptation, which could be exploited for additional research on avian species in future.

Analysis of codon usage preference in hemagglutinin genes of the swine-origin influenza A (H1N1) virus

BACKGROUND

The swine-origin influenza A (H1N1) virus (S-OIV) has come to the forefront since 2009 and was identified as a new reassortant strain. The hemagglutinin (HA) glycoprotein mediates virus binding, contains antigenic regions recognized by neutralizing antibodies, and is associated with viral cross-species infection and adaption. The comparison study of codon usage preferences in influenza viral genomes was less extensive. In this study, we used codon usage pattern analyses to validate the adaption and origins of S-OIV.

METHODS

Codon usage pattern was used to estimate the host adaption of S-OIVs. Phylogenetic analysis of the HA gene was conducted to understand the phylogeny of H1N1 viruses isolated from different hosts. Amino acid signature pattern on antigenic sites of HA was analyzed to understand the antigenic characteristics.

RESULTS

Results of phylogenetic analyses of HA gene indicate that S-OIVs group in identical clusters. The synonymous codon usage pattern analyses indicate that the effective number of codons versus GC content at the third codon position in the HA1 gene slightly differ from those in swine H1N1 and gradually adapted to human. Our data indicate that S-OIV evolution occurred according to positive selection within these antigenic regions. A comparison of antigenic site amino acids reveals similar signature patterns between S-OIV and 1918 human influenza strains.

CONCLUSION

This study proposes a new and effective way to gain a better understanding of the features of the S-OIV genome and evolutionary processes based on the codon usage pattern. It is useful to trace influenza viral origins and cross-species virus transmission.

Analysis of synonymous codon usage in Newcastle disease virus hemagglutinin-neuraminidase (HN) gene and fusion protein (F) gene

Newcastle disease virus (NDV) hemagglutinin-neuraminidase (HN) is a multifunctional protein, which possesses both the receptor recognition and neuraminidase activities. The fusion (F) protein is a type I membrane glycoprotein that mediates the merger of the viral envelope to the host cell membrane. Although the functions of the HN and F proteins have been well studied, however, the factors shaping synonymous codon usage bias and nucleotide composition in HN and F genes have been few reported. In our study, we analyzed synonymous codon usage using the 69 NDV HN and F genes, respectively. The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in HN and F genes. In addition, other factors, such as the aromaticity and hydrophobicity, also influence the codon usage variation among HN and F genes. This study represents the most comprehensive analysis to date of NDV HN and F genes codon usage patterns and provides a basic understanding of the mechanisms for codon usage bias.

Large-scale analysis of conserved rare codon clusters suggests an involvement in co-translational molecular recognition events

MOTIVATION

An increasing amount of evidence from experimental and computational analysis suggests that rare codon clusters are functionally important for protein activity. Most of the studies on rare codon clusters were performed on a limited number of proteins or protein families. In the present study, we present the Sherlocc program and how it can be used for large scale protein family analysis of evolutionarily conserved rare codon clusters and their relation to protein function and structure. This large-scale analysis was performed using the whole Pfam database covering over 70% of the known protein sequence universe. Our program Sherlocc, detects statistically relevant conserved rare codon clusters and produces a user-friendly HTML output.

RESULTS

Statistically significant rare codon clusters were detected in a multitude of Pfam protein families. The most statistically significant rare codon clusters were predominantly identified in N-terminal Pfam families. Many of the longest rare codon clusters are found in membrane-related proteins which are required to interact with other proteins as part of their function, for example in targeting or insertion. We identified some cases where rare codon clusters can play a regulating role in the folding of catalytically important domains. Our results support the existence of a widespread functional role for rare codon clusters across species. Finally, we developed an online filter-based search interface that provides access to Sherlocc results for all Pfam families.

AVAILABILITY

The Sherlocc program and search interface are open access and are available at http://bcb.med.usherbrooke.ca

Characterization of Codon Usage Bias in UL21gene from Duck Enteritis Virus

In this paper, the comprehensive analysis of codon usage bias of Duck enteritis virus (DEV) UL21 gene was performed by using CAI, CHIPS and CUSP program of EMBOSS. Our study showed that codon usage bias of DEV UL21 had strong bias towards the A-ended or T-ended codons, and GC3s contents of the codon usage bias in DEV UL21 gene were significantly varied compared with those of other 27 reference herpesviruses. The CAI, ENC value of DEV CHv strain UL21 gene is 0.615 and 55.167, respectively, indicating that the codon usage bias of this gene is weak and lowly expressed. The plot of ENC versus GC3S content revealed that DEV UL21 gene is subject to GC compositional constraints. The phylogentic analysis about amino acids codon usage bias of DEV UL21 and the27 reference herpesviruses showed that DEV was evolutionarily closer to herpesviruses Mardivirus. In addition, the codon usage bias of DEV UL21 gene was compared with those of E. coli, yeast and humans. There are 42, 45, 39 same codons usage bias between the DEV UL21 to E.coli, Yeast, H.sapiens, respectively, indicaiting that UL21 gene of DEV may be more efficiently expressed in the yeast system.

The characteristics of the synonymous codon usage in hepatitis B virus and the effects of host on the virus in codon usage pattern

BACKGROUND

Hepatitis B virus (HBV) infection is one of the main human health problem and causes a large-scale of patients chronic infection worldwide.. As the replication of HBV depends on its host cell system, codon usage pattern for the viral gene might be susceptible to two main selections, namely mutation pressure and translation selection. In this case, a deeper investigation between HBV evolution and host adaptive response might assist control this disease.

RESULT

Relative synonymous codon usage (RSCU) values for the whole HBV coding sequence were studied by Principal component analysis (PCA). The characteristics of the synonymous codon usage patterns, nucleotide contents and the comparison between ENC values of the whole HBV coding sequence indicated that the interaction between virus mutation pressure and host translation selection exists in the processes of HBV evolution. The synonymous codon usage pattern of HBV is a mixture of coincidence and antagonism to that of host cell. But the difference of genetic characteristic of HBV failed to be observed to its different epidemic areas or subtypes, suggesting that geographic factor is limited to influence the evolution of this virus, while genetic characteristic based on HBV genotypes could be divided into three groups, namely (i) genotyps A and E, (ii) genotype B, (iii) genotypes C, D and G.

CONCLUSION

Codon usage patterns from PCA for identification of evolutionary trends in HBV provide an alternative approach to understand the evolution of HBV. Further more, a combined selection of mutation pressure with translation selection on codon usage might shed a light on understanding the evolutionary trends of HBV genotypes.

Analysis of synonymous codon usage inShigella flexneri2a strain 301 and otherShigellaandEscherichia colistrains

In this study, we analysed synonymous codon usage in Shigella flexneri 2a strain 301 (Sf301) and performed a comparative analysis of synonymous codon usage patterns in Sf301 and other strains of Shigella and Escherichia coli . Although there was a significant variety in codon usage bias among different Sf301 genes, there was a slight but observable codon usage bias that could primarily be attributable to mutational pressure and translational selection. In addition, the relative abundance of dinucleotides in Sf301 was observed to be independent of the overall base composition but was still caused by differential mutational pressure; this also shaped codon usage. By comparing the relative synonymous codon usage values across different Shigella and E. coli strains, we suggested that the synonymous codon usage pattern in the Shigella genomes was strain specific. This study represents a comprehensive analysis of Shigella codon usage patterns and provides a basic understanding of the mechanisms underlying codon usage bias.

Bioinformatics analysis of rabbit haemorrhagic disease virus genome

Background

Rabbit haemorrhagic disease virus (RHDV), as the pathogeny of Rabbit haemorrhagic disease, can cause a highly infectious and often fatal disease only affecting wild and domestic rabbits. Recent researches revealed that it, as one number of the Caliciviridae, has some specialties in its genome, its reproduction and so on.

Results

In this report, we firstly analyzed its genome and two open reading frameworks (ORFs) from this aspect of codon usage bias. Our researches indicated that mutation pressure rather than natural is the most important determinant in RHDV with high codon bias, and the codon usage bias is nearly contrary between ORF1 and ORF2, which is maybe one of factors regulating the expression of VP60 (encoding by ORF1) and VP10 (encoding by ORF2). Furthermore, negative selective constraints on the RHDV whole genome implied that VP10 played an important role in RHDV lifecycle.

Conclusions

We conjectured that VP10 might be beneficial for the replication, release or both of virus by inducing infected cell apoptosis initiate by RHDV. According to the results of the principal component analysis for ORF2 of RSCU, we firstly separated 30 RHDV into two genotypes, and the ENC values indicated ORF1 and ORF2 were independent among the evolution of RHDV.

Analysis of synonymous codon usage in hepatitis A virus

BACKGROUND

Hepatitis A virus is the causative agent of type A viral hepatitis, which causes occasional acute hepatitis. Nevertheless, little information about synonymous codon usage pattern of HAV genome in the process of its evolution is available. In this study, the key genetic determinants of codon usage in HAV were examined.

RESULTS

The overall extent of codon usage bias in HAV is high in Picornaviridae. And the patterns of synonymous codon usage are quite different in HAV genomes from different location. The base composition is closely correlated with codon usage bias. Furthermore, the most important determinant that results in such a high codon bias in HAV is mutation pressure rather than natural selection.

CONCLUSIONS

HAV presents a higher codon usage bias than other members of Picornaviridae. Compositional constraint is a significant element that influences the variation of synonymous codon usage in HAV genome. Besides, mutation pressure is supposed to be the major factor shaping the hyperendemic codon usage pattern of HAV.

A detailed comparative analysis on the overall codon usage patterns in hepatitis A virus

Hepatitis A virus (HAV) is a hepatotropic member of the family Picornaviridae. HAV has several unique biological characteristics that distinguish it from other members of this family. Recent and previous studies revealed that codon usage plays a key role in HAV replication and evolution. In this study, the patterns of synonymous codon usage in HAV have been studied through multivariate statistical methods on 30 complete open reading frames (ORFs) from the available 30 full-length HAV sequences. Effective number of codons (ENC) indicates that the overall extent of codon usage bias in HAV genomes is significant. The relative dinucleotide abundances suggest that codon usage in HAV can also be strongly influenced by underlying biases in dinucleotide frequencies. These factors strongly correlated with the first major axis of correspondence analysis (COA) on relative synonymous codon usage (RSCU). The distribution of the HAV ORFs along the plane defined by the first two major axes in COA showed that different genotypes are located at different places in the plane, suggesting that HAV codon usage is also reflecting an evolutionary process. It has been very recently described that fine-tuning translation kinetics selection also contributes to codon usage bias of HAV. The results of these studies suggest that HAV genomic biases are the result of the co-evolution of genome composition, controlled translation kinetics and probably the ability to escape the antiviral cell responses.

The characteristics of the synonymous codon usage in enterovirus 71 virus and the effects of host on the virus in codon usage pattern

To give a new perspective on the evolutionary characteristics shaping the genetic diversity of enterovirus 71 (EV71) and the effects of natural selection from its host on the codon usage pattern of the virus, the relative synonymous codon usage (RSCU) values, codon usage bias (CUB) values, effective number of codons (ENCs) values and nucleotide contents were calculated to implement a comparative analysis to evaluate the dynamics of the virus evolution. The characteristics of the synonymous codon usage patterns and nucleotide contents of EV71 and the comparison between ENC values for the whole coding sequence of EV71 and that of coding sequences for viral proteins of EV71 all indicate that the interaction between mutation pressure from virus and natural selection from host exists in the processes of evolution of EV71. The synonymous codon usage pattern of EV71 is a mixture of coincidence and antagonism to that of host cell. In addition, the genetic diversity of EV71 strains and the preferential selection of some synonymous codons in EV71 strains based on the different epidemic areas were observed, suggesting that geographic and social factors may play roles in influencing the evolution of this virus.

Adaptive evolution of the Chlamydia trachomatis dominant antigen reveals distinct evolutionary scenarios for B- and T-cell epitopes: worldwide survey

Background

Chlamydia trachomatis is one of the most disseminated human pathogens, for which no vaccine is available yet. Understanding the impact of the host pressure on pathogen antigens is crucial, but so far it was only assessed for highly-restricted geographic areas. We aimed to evaluate the evolutionary picture of the chlamydial key antigen (MOMP), which is one of the leading multi-subunit vaccine candidates, in a worldwide basis.

Methodology/Principal Findings

Using genetics, molecular evolution methods and mathematical modelling, we analyzed all MOMP sequences reported worldwide, composed by 5026 strains from 33 geographic regions of five continents. Overall, 35.9% of variants were detected. The evolutionary pattern of MOMP amino acid gains/losses was found to differ from the remaining chromosome, reflecting the demanding constraints of this porin, adhesin and dominant antigen. Amino acid changes were 4.3-fold more frequent in host-interacting domains (P<10⁻¹²), specifically within B-cell epitopes (P<10⁻⁵), where 25% of them are at fixation (P<10⁻⁵). According to the typical pathogen-host arms race, this rampant B-cell antigenic variation likely represents neutralization escape mutants, as some mutations were previously shown to abrogate neutralization of chlamydial infectivity in vitro. In contrast, T-cell clusters of diverse HLA specificities are under purifying selection, suggesting a strategy that may lead to immune subversion. Moreover, several silent mutations are at fixation, generating preferential codons that may influence expression, and may also reflect recombination-derived ‘hitchhiking-effect’ from favourable nonsilent changes. Interestingly, the most prevalent C. trachomatis genotypes, E and F, showed a mutation rate 22.3-fold lower than that of the remainder (P<10⁻²⁰), suggesting more fitted antigenic profiles.

Conclusions/Significance

Globally, the adaptive evolution of the C. trachomatis dominant antigen is likely driven by its complex pathogenesis-related function and reflects distinct evolutionary antigenic scenarios that may benefit the pathogen, and thus should be taking into account in the development of a MOMP-based vaccine.

Analysis of synonymous codon usage in porcine reproductive and respiratory syndrome virus

In this study, we calculated the relative synonymous codon usage (RSCU) values and codon usage bias (CUB) values to implement a comparative analysis of codon usage pattern of open reading frames (ORFs) which belong to the two main genotypes of porcine reproductive and respiratory syndrome virus (PRRSV). By analysis of synonymous codon usage values in each ORF of PRRSV, the optimal codons for most amino acids were all C or G-ended codons except GAU for Asp, CAU for His, UUU for Phe and CCU for Pro. The synonymous codon usage patterns in different ORFs of PRRSV were different and genetically conserved. Among them, ORF1a, ORF4, ORF5 and ORF7 could cluster these strains into the two main serotypes (EU and US). Due to mutational pressure, compositional constraint played an important role in shaping the synonymous codon usage pattern in different ORFs, and the synonymous codon usage diversity in ORFs was correlated with gene function. The degree of CUB for some particular amino acids under strong selection pressure probably served as a potential genetic marker for each ORF in PRRSV. However, gene length and translational selection in nature had no effect on the synonymous codon usage pattern in PRRSV. These conclusions could not only offer an insight into the synonymous codon usage pattern and differentiation of gene function, but also assist in understanding the discrepancy of evolution among ORFs in PRRSV.

Analysis of synonymous codon usage in foot-and-mouth disease virus

In this study, we calculate the relative synonymous codon usage (RSCU) values and codon usage bias (CUB) values to carry out a comparative analysis of codon usage pattern for open reading frames (ORFs) among 85 samples which belong to all seven serotypes of foot-and-mouth disease virus (FMDV). Although the degree of CUB for ORFs is a relatively slight, there is a significant variation for CUB among different serotypes, which is mainly determined by codon usage pattern depending on RSCU. By comparison with RSCU values for all samples, although RSCU values fail to show the relationship of specific-lineage serotype, there are two main genetic populations existing in FMDV, namely (i) serotypes Asia 1, A, C & O; (ii) serotypes SAT 1, 2 & 3. This interesting characteristic may be formed by the mechanism of RNA virus recombination. The analysis of quantitative & qualitative evaluation based on CUB indicates interesting characteristic of codon usage, which suggests that more FMDV genome diversity may exist in specific-lineage serotypes rather than exist randomly. Furthermore, the relationship between amino acids and codon usage pattern indicates that mutation pressure rather than translational selection in nature is the important determinant of the codon usage bias observed. Our work might give some sight into some characteristics of FMDV ORF and some evolutionary information of this virus.

Analysis of synonymous codon usage in Aeropyrum pernix K1 and other Crenarchaeota microorganisms

In this study, a comparative analysis of the codon usage bias was performed in Aeropyrum pernix K1 and two other phylogenetically related Crenarchaeota microorganisms (i.e., Pyrobaculum aerophilum str. IM2 and Sulfolobus acidocaldarius DSM 639). The results indicated that the synonymous codon usage in A. pernix K1 was less biased, which was highly correlated with the GC(3S) value. The codon usage patterns were phylogenetically conserved among these Crenarchaeota microorganisms. Comparatively, it is the species function rather than the gene function that determines their gene codon usage patterns. A. pernix K1, P. aerophilum str. IM2, and S. acidocaldarius DSM 639 live in differently extreme conditions. It is presumed that the living environment played an important role in determining the codon usage pattern of these microorganisms. Besides, there was no strain-specific codon usage among these microorganisms. The extent of codon bias in A. pernix K1 and S. acidocaldarius DSM 639 were highly correlated with the gene expression level, but no such association was detected in P. aerophilum str. IM2 genomes.

Analysis of synonymous codon usage in classical swine fever virus

Using the complete genome sequences of 35 classical swine fever viruses (CSFV) representing all three genotypes and all three kinds of virulence, we analyzed synonymous codon usage and the relative dinucleotide abundance in CSFV. The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in CSFV. Furthermore, we observed that the relative abundance of dinucleotides in CSFV is independent of the overall base composition but is still the result of differential mutational pressure, which also shapes codon usage. In addition, other factors, such as the subgenotypes and aromaticity, also influence the codon usage variation among the genomes of CSFV. This study represents the most comprehensive analysis to date of CSFV codon usage patterns and provides a basic understanding of the mechanisms for codon usage bias.

Analysis of synonymous codon usage in the UL24 gene of duck enteritis virus

The analysis on codon usage bias of UL24 gene of duck enteritis virus (DEV) may improve our understanding of the evolution and pathogenesis of DEV and provide a basis for understanding the relevant mechanism for biased usage of synonymous codons and for selecting appropriate expression systems to improve the expression of target genes. The codon usage bias of UL24 genes of DEV and 27 reference herpesviruses were analyzed. The results showed that codon of UL24 gene of DEV was strong bias toward the synonymous codons with A and T at the third codon position. A high level of diversity in codon usage bias existed, and the effective number of codons used in a gene plot revealed that the genetic heterogeneity in UL24 gene of herpesviruses was constrained by the G + C content. The phylogentic analysis suggested that DEV was evolutionarily closer to Alphaherpesvirinae and that there was no significant deviation in codon usage in different virus strains. There were 20 codons showing distinct usage differences between DEV and Escherichia coli, 23 between DEV and Homo sapiens, but only 16 codons between DEV and yeast. Therefore the yeast expression system may be more suitable for the expression of DEV genes.

Analysis of codon usage bias and base compositional constraints in iridovirus genomes

The codon usage bias and the base composition variations in the available 12 complete iridovirus genome sequences have been investigated. We re-evaluated the number of open reading frames (ORFs) in each published iridovirus genome and analyzed its correlation against the genome size. The result shows that there is a direct relationship between the number of ORFs and the genome size. The codon usage patterns of these iridoviruses are found to be phylogenetically conserved. A significant variation in the base content among the 12 iridovirus genomes has been observed, with G+C content ranges widely from 27 to 55%. Moreover, the preferential use of bases in codons is different among higher and lower G+C content genomes. A preferential codon usage among viral genomes is also noticed. Effective number of codon (Enc) plot reveals that the G+C compositional constraint is the main factor that determines the codon usage bias. Relative synonymous codon usage analysis of methyltransferase containing as well as lacking viruses suggests that the codon usage is not influenced by the methylation-mediated mutation. In addition, the comparison of the codon usage of iridovirus hosts and the iridovirus genomes reveals that the host tRNA pool may be responsible for the base compositional constraint. This study represents the most comprehensive analysis to date for iridovirus codon usage patterns.

Synonymous codon usage bias in 16 Staphylococcus aureus phages: implication in phage therapy

To reveal the factors influencing architecture of protein-coding genes in staphylococcal phages, relative synonymous codon usage variation has been investigated in 920 protein-coding genes of 16 staphylococcal phages. As expected for AT rich genomes, there are predominantly A and T ending codons in all 16 phages. Both Nc plot and correspondence analysis on relative synonymous codon usage indicates that mutation bias influences codon usage variation in the 16 phages. Correspondence analysis also suggests that translational selection and gene length also influence the codon usage variation in the phages to some extent and codon usage in staphylococcal phages is phage-specific but not S. aureus-specific. Further analysis indicates that among 16 staphylococcal phages, 44AHJD, P68 and K may be extremely virulent in nature as most of their genes have high translation efficiency. If this is true, then above three phages may be useful for curing staphylococcal infections.

Modelling of translation of human protein disulfide isomerase in Escherichia coli—A case study of gene optimisation

Recombinant human protein disulfide isomerase (PDI) was expressed in vivo in Escherichia coli using a non-optimised gene sequence and an optimised sequence with four 5' codons substituted by synonymous codons that take less time to translate. The optimisation resulted in a 2-fold increase of total PDI concentration and by successive optimisation with expression at low temperature in a 10-fold increase of the amount of soluble PDI in comparison with the original wild-type construct. The improvement can be due to a faster clearing of the ribosome binding site on the mRNA, elevating the translation initiation rate and resulting in higher ribosome loading and better ribosome protection of the PDI mRNA against endonucleolytic cleavage by RNase. This hypothesis was supported by a novel computer simulation model of E. coli translational ribosome traffic based upon the stochastic Gillespie algorithm. The study indicates the applicability of such models in optimisation of recombinant protein sequences.

Factors influencing the synonymous codon and amino acid usage bias in AT-rich Pseudomonas aeruginosa phage PhiKZ

To reveal how the AT-rich genome of bacteriophage PhiKZ has been shaped in order to carry out its growth in the GC-rich host Pseudomonas aeruginosa, synonymous codon and amino acid usage bias of PhiKZ was investigated and the data were compared with that of P. aeruginosa. It was found that synonymous codon and amino acid usage of PhiKZ was distinct from that of P. aeruginosa. In contrast to P. aeruginosa, the third codon position of the synonymous codons of PhiKZ carries mostly A or T base; codon usage bias in PhiKZ is dictated mainly by mutational bias and, to a lesser extent, by translational selection. A cluster analysis of the relative synonymous codon usage values of 16 myoviruses including PhiKZ shows that PhiKZ is evolutionary much closer to Escherichia coli phage T4. Further analysis reveals that the three factors of mean molecular weight, aromaticity and cysteine content are mostly responsible for the variation of amino acid usage in PhiKZ proteins, whereas amino acid usage of P. aeruginosa proteins is mainly governed by grand average of hydropathicity, aromaticity and cysteine content. Based on these observations, we suggest that codons of the phage-like PhiKZ have evolved to preferentially incorporate the smaller amino acid residues into their proteins during translation, thereby economizing the cost of its development in GC-rich P. aeruginosa.

Comparative Analysis of the Base Composition and Codon Usages in Fourteen Mycobacteriophage Genomes

To study the possible codon usage and base composition variation in the bacteriophages, fourteen mycobacteriophages were used as a model system here and both the parameters in all these phages and their plating bacteria, M. smegmatis had been determined and compared. As all the organisms are GC-rich, the GC contents at third codon positions were found in fact higher than the second codon positions as well as the first + second codon positions in all the organisms indicating that directional mutational pressure is strongly operative at the synonymous third codon positions. Nc plot indicates that codon usage variation in all these organisms are governed by the forces other than compositional constraints. Correspondence analysis suggests that: (i) there are codon usage variation among the genes and genomes of the fourteen mycobacteriophages and M. smegmatis, i.e., codon usage patterns in the mycobacteriophages is phage-specific but not the M. smegmatis-specific; (ii) synonymous codon usage patterns of Barnyard, Che8, Che9d, and Omega are more similar than the rest mycobacteriophages and M. smegmatis; (iii) codon usage bias in the mycobacteriophages are mainly determined by mutational pressure; and (iv) the genes of comparatively GC rich genomes are more biased than the GC poor genomes. Translational selection in determining the codon usage variation in highly expressed genes can be invoked from the predominant occurrences of C ending codons in the highly expressed genes. Cluster analysis based on codon usage data also shows that there are two distinct branches for the fourteen mycobacteriophages and there is codon usage variation even among the phages of each branch.

Analysis of synonymous codon usage in H5N1 virus and other influenza A viruses

In this study, we calculated the codon usage bias in H5N1 virus and performed a comparative analysis of synonymous codon usage patterns in H5N1 virus, five other evolutionary related influenza A viruses and a influenza B virus. Codon usage bias in H5N1 genome is a little slight, which is mainly determined by the base compositions on the third codon position. By comparing synonymous codon usage patterns in different viruses, we observed that the codon usage pattern of H5N1 virus is similar with other influenza A viruses, but not influenza B virus, and the synonymous codon usage in influenza A virus genes is phylogenetically conservative, but not strain-specific. Synonymous codon usage in genes encoded by different influenza A viruses is genus conservative. Compositional constraints could explain most of the variation of synonymous codon usage among these virus genes, while gene function is also correlated to synonymous codon usages to a certain extent. However, translational selection and gene length have no effect on the variations of synonymous codon usage in these virus genes.

The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens

The role of silent position in the codon on the protein structure is an interesting and yet unclear problem. In this paper, 563 Homo sapiens genes and 417 Escherichia coli genes coding for proteins with four different folding types have been analyzed using variance analysis, a multivariate analysis method newly used in codon usage analysis, to find the correlation between amino acid composition, synonymous codon, and protein structure in different organisms. It has been found that in E. coli, both amino acid compositions in differently folded proteins and synonymous codon usage in different gene classes coding for differently folded proteins are significantly different. It was also found that only amino acid composition is different in different protein classes in H. sapiens. There is no universal correlation between synonymous codon usage and protein structure in these two different organisms. Further analysis has shown that GC content on the second codon position can distinguish coding genes for different folded proteins in both organisms.

Analysis of synonymous codon usage in SARS Coronavirus and other viruses in the Nidovirales

In this study, we calculated the codon usage bias in severe acute respiratory syndrome Coronavirus (SARSCoV) and performed a comparative analysis of synonymous codon usage patterns in SARSCoV and 10 other evolutionary related viruses in the Nidovirales. Although there is a significant variation in codon usage bias among different SARSCoV genes, codon usage bias in SARSCoV is a little slight, which is mainly determined by the base compositions on the third codon position. By comparing synonymous codon usage patterns in different viruses, we observed that synonymous codon usage pattern in these virus genes was virus specific and phylogenetically conserved, but it was not host specific. Phylogenetic analysis based on codon usage pattern suggested that SARSCoV was diverged far from all three known groups of Coronavirus. Compositional constraints could explain most of the variation of synonymous codon usage among these virus genes, while gene function is also correlated to synonymous codon usages to a certain extent. However, translational selection and gene length have no effect on the variations of synonymous codon usage in these virus genes.

Mean-field approaches to the totally asymmetric exclusion process with quenched disorder and large particles

The process of protein synthesis in biological systems resembles a one-dimensional driven lattice gas in which the particles (ribosomes) have spatial extent, covering more than one lattice site. Realistic, nonuniform gene sequences lead to quenched disorder in the particle hopping rates. We study the totally asymmetric exclusion process with large particles and quenched disorder via several mean-field approaches and compare the mean-field results with Monte Carlo simulations. Mean-field equations obtained from the literature are found to be reasonably effective in describing this system. A numerical technique is developed for computing the particle current rapidly. The mean-field approach is extended to include two-point correlations between adjacent sites. The two-point results are found to match Monte Carlo simulations more closely.

Totally asymmetric exclusion process with extended objects: a model for protein synthesis

The process of protein synthesis in biological systems resembles a one dimensional driven lattice gas in which the particles have spatial extent, covering more than one lattice site. We expand the well studied totally asymmetric exclusion process, in which particles typically cover a single lattice site, to include cases with extended objects. Exact solutions can be determined for a uniform closed system. We analyze the uniform open system through two approaches. First, a continuum limit produces a modified diffusion equation for particle density profiles. Second, an extremal principle based on domain wall theory accurately predicts the phase diagram and currents in each phase. Finally, we briefly consider approximate approaches to a nonuniform open system with quenched disorder in the particle hopping rates and compare these approaches with Monte Carlo simulations.

Cis control of gene expression in E.coli by ribosome queuing at an inefficient translational stop signal

An UGA stop codon context which is inefficient because of the 3'-flanking context and the last two amino acids in the gene protein product has a negative effect on gene expression, as shown using a model protein A' gene. This is particularly true at low mRNA levels, corresponding to a high intracellular ribosome/mRNA ratio. The negative effect is smaller if this ratio is decreased, or if the distance between the initiation and termination signals is increased. The results suggest that an inefficient termination codon can cause ribosomal pausing and queuing along the upstream mRNA region, thus blocking translation initiation of short genes. This cis control effect is dependent on the stop codon context, including the C-terminal amino acids in the gene product, the translation initiation signal strength, the ribosome/mRNA ratio and the size of the mRNA coding region. A large proportion of poorly expressed natural Escherichia coli genes are small, and the weak termination codon UGA is under-represented in small, highly expressed E.coli genes as compared with the efficient stop codon UAA.

An insight into the possible mechanism of working of two-cistronic gene expression systems and rational designing of newer systems

The initial attempts at hyper-expressing buffalo/goat growth hormone (GH)-ORFs in Escherichia coli directly under various strong promoters were not successful despite the presence of a functional gene. High level expression of GH was achieved as a fusion protein with glutathione-S-transferase (GST). To produce native GH in an unfused state, we adapted an established strategy of two-cistronic approach in our system. In this strategy, utilizing one of the highly efficient reported sequences as the first cistron led to a nearly 1000-fold enhancement in the level of expression under an E. coli promoter (trc). In search of a newer first-cistron sequence as well as to see the generality of the two-cistronic approach, we explored the ability of different lengths of a highly expressing natural gene to act as an efficient first cistron. Surprisingly, GST, which is naturally highly expressible in E. coli, could not be fitted into a successful two-cistronic construct. In addition, placement of the entire two-cistronic expression cassette (which had earlier given high-level GH expression under trc promoter) under the T7 promoter in E. coli failed to hyper-express GH. These results suggest that the successful exploitation of the two-cistron arrangement for hyper-expression of eukaryotic ORFs in bacteria is not as straightforward as was previously thought. It appears probable that factors such as the sequence context, together with the length and codons used in the first cistron are important as well.

Conformational toxicity and sporadic conformational diseases

Spontaneous, so-called 'conformational' diseases, specially of the neurodegenerative type like Alzheimer's, are linked to certain protein types which have the normal amino-acid sequence but are misfolded and accumulate due to resistance to proteolysis. In the case of prion diseases, the 'protein only' hypothesis assumes that the misconformation of a native protein could be initiated upon interaction with a sister-protein already in the misfolded state. There is an alternative to this sister protein contamination scheme, which assumes that the misconformation is acquired upon protein synthesis, that is de novo. Misfoldling and resistance to proteolysis could result from defects responsible for shortage or inactivity of the cellular factors in charge of protein folding and degradation. The defects could have a genetic origin (the gene of the faulty factor involved could have been mutated, or control and regulation of its expression could have been altered, etc.). Alternatively, the cell's actual biosynthetic and/or proteolytic resources could have become overloaded and unavailable, due to unscheduled mass-production of proteins resulting from unscheduled cell growth or proliferation, cell stress, etc. Xenobiotics, active for instance as endocrine proliferators, stressors, or inducing copious, unscheduled gene expression, etc. could give rise to shortage of cellular factors necessary for the production of native proteins and for proteolysis. Alternatively, xenobiotics could alter expression or activity of some of these factors. In both cases, the xenobiotic could be a 'conformational toxicant' by inducing misfolding of selected proteins. The xenobiotic could trigger some conformational disease if it targets a specific protein and tissue.