Codon usage and gene function are related in sequences of Arabidopsis thaliana

Codon usage bias and genetic diversity in chloroplast genomes of Elaeagnus species (Myrtiflorae: Elaeagnaceae)

Codon usage bias (CUB) reveals the characteristics of species and can be utilized to understand their evolutionary relationship, increase the target genes' expression in the heterologous receptor plants, and further provide theoretic assistance for correlative study on molecular biology and genetic breeding. The chief aim of this work was to analyze the CUB in chloroplast (cp.) genes in nine Elaeagnus species to provide references for subsequent studies. The codons of Elaeagnus cp. genes preferred to end with A/T bases rather than with G/C bases. Most of the cp. genes were prone to mutation, while the rps7 genes were identical in sequences. Natural selection was inferred to have a powerful impact on the CUB in Elaeagnus cp. genomes, and their CUB was extremely strong. In addition, the optimal codons were identified in the nine cp. genomes based on the relative synonymous codon usage (RSCU) values, and the optimal codon numbers were between 15 and 19. The clustering analyses based on RSCU were contrasted with the maximum likelihood (ML)-based phylogenetic tree derived from coding sequences, suggesting that the t-distributed Stochastic Neighbor Embedding clustering method was more appropriate for evolutionary relationship analysis than the complete linkage method. Moreover, the ML-based phylogenetic tree based on the conservative matK genes and the whole cp. genomes had visible differences, indicating that the sequences of specific cp. genes were profoundly affected by their surroundings. Following the clustering analysis, Arabidopsis thaliana was considered the optimal heterologous expression receptor plant for the Elaeagnus cp. genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01289-6.

Genome-wide analysis of codon usage in sesame (Sesamum indicum L.)

Sesamum indicum is an ancient oil crop grown in tropical and subtropical areas of the world. We have analyzed 23,538 coding sequences (CDS) of S. indicum to understand the factors shaping codon usage in this important oil crop plant. We identified eleven highly preferred codons in S. indicum that have AT-endings. The slope of a neutrality plot was less than one while effective number of codons (ENC) plot showed distribution above and below the standard curve. There is a significant relationship between protein length and relative synonymous codon usage (RSCU) at the primary axis while there is a weak correlation between protein length and Nc values. Correspondence analysis conducted on RSCU values differentiated CDS based on their GC content and their characteristic feature and showed a discrete distribution. Moreover, by determining codon usage, we found out that majority of the lignan biosynthesis related genes showed a weaker codon usage bias. These results provide insights into understanding codon evolution in sesame.

Codon usage bias

Codon usage bias is the preferential or non-random use of synonymous codons, a ubiquitous phenomenon observed in bacteria, plants and animals. Different species have consistent and characteristic codon biases. Codon bias varies not only with species, family or group within kingdom, but also between the genes within an organism. Codon usage bias has evolved through mutation, natural selection, and genetic drift in various organisms. Genome composition, GC content, expression level and length of genes, position and context of codons in the genes, recombination rates, mRNA folding, and tRNA abundance and interactions are some factors influencing codon bias. The factors shaping codon bias may also be involved in evolution of the universal genetic code. Codon-usage bias is critical factor determining gene expression and cellular function by influencing diverse processes such as RNA processing, protein translation and protein folding. Codon usage bias reflects the origin, mutation patterns and evolution of the species or genes. Investigations of codon bias patterns in genomes can reveal phylogenetic relationships between organisms, horizontal gene transfers, molecular evolution of genes and identify selective forces that drive their evolution. Most important application of codon bias analysis is in the design of transgenes, to increase gene expression levels through codon optimization, for development of transgenic crops. The review gives an overview of deviations of genetic code, factors influencing codon usage or bias, codon usage bias of nuclear and organellar genes, computational methods to determine codon usage and the significance as well as applications of codon usage analysis in biological research, with emphasis on plants.

Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB

Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite-induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duodenalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G. duodenalis CUB formation. This study systematically analyzes G. duodenalis codon usage pattern and clarifies the mechanisms of G. duodenalis CUB. These results will be very useful to identify new genes, molecular genetic manipulation, and study of G. duodenalis evolution.

The Complete Chloroplast Genome of the Vulnerable Oreocharis esquirolii (Gesneriaceae): Structural Features, Comparative and Phylogenetic Analysis

Oreocharis esquirolii, a member of Gesneriaceae, is known as Thamnocharis esquirolii, which has been regarded a synonym of the former. The species is endemic to Guizhou, southwestern China, and is evaluated as vulnerable (VU) under the International Union for Conservation of Nature (IUCN) criteria. Until now, the sequence and genome information of O. esquirolii remains unknown. In this study, we assembled and characterized the complete chloroplast (cp) genome of O. esquirolii using Illumina sequencing data for the first time. The total length of the cp genome was 154,069 bp with a typical quadripartite structure consisting of a pair of inverted repeats (IRs) of 25,392 bp separated by a large single copy region (LSC) of 85,156 bp and a small single copy region (SSC) of18,129 bp. The genome comprised 114 unique genes with 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. Thirty-one repeat sequences and 74 simple sequence repeats (SSRs) were identified. Genome alignment across five plastid genomes of Gesneriaceae indicated a high sequence similarity. Four highly variable sites (rps16-trnQ, trnS-trnG, ndhF-rpl32, and ycf 1) were identified. Phylogenetic analysis indicated that O. esquirolii grouped together with O. mileensis, supporting resurrection of the name Oreocharis esquirolii from Thamnocharisesquirolii. The complete cp genome sequence will contribute to further studies in molecular identification, genetic diversity, and phylogeny.

Demographic history and adaptive synonymous and nonsynonymous variants of nuclear genes in Rhododendron oldhamii (Ericaceae)

Demographic events are important in shaping the population genetic structure and exon variation can play roles in adaptive divergence. Twelve nuclear genes were used to investigate the species-level phylogeography of Rhododendron oldhamii, test the difference in the average GC content of coding sites and of third codon positions with that of surrounding non-coding regions, and test exon variants associated with environmental variables. Spatial expansion was suggested by R₂ index of the aligned intron sequences of all genes of the regional samples and sum of squared deviations statistic of the aligned intron sequences of all genes individually and of all genes of the regional and pooled samples. The level of genetic differentiation was significantly different between regional samples. Significantly lower and higher average GC contents across 94 sequences of the 12 genes at third codon positions of coding sequences than that of surrounding non-coding regions were found. We found seven exon variants associated strongly with environmental variables. Our results demonstrated spatial expansion of R. oldhamii in the late Pleistocene and the optimal third codon position could end in A or T rather than G or C as frequent alleles and could have been important for adaptive divergence in R. oldhamii.

Transit From Autotrophism to Heterotrophism: Sequence Variation and Evolution of Chloroplast Genomes in Orobanchaceae Species

The family Orobanchaceae including autotrophic, hemiparasitic, and holoparasitic species, is becoming a key taxa to study the evolution of chloroplast genomes in different lifestyles. But the early evolutionary trajectory in the transit from autotrophism to hemiparasitism still maintains unclear for the inadequate sampling. In this study, we compared 50 complete chloroplast genomes in Orobanchaceae, containing four newly sequenced plastomes from hemiparasitic Pedicularis, to elucidate the sequence variation patterns in the evolution of plastomes. Contrasted to the sequence and structural hypervariabilities in holoparasites, hemiparasitic plastomes exhibited high similarity to those of autotrophs in gene and GC contents. They are generally characterized with functional or physical loss of ndh/tRNA genes and the inverted small-single-copy region. Gene losses in Orobanchaceae were lineage-specific and convergent, possibly related to structural reconfiguration and expansion/contraction of the inverted region. Pseudogenization of ndh genes was unique in hemiparasites. At least in Pedicularis, the ndhF gene might be most sensitive to the environmental factors and easily pseudogenized when autotrophs transit to hemiparasites. And the changes in gene contents and structural variation potentially deeply rely on the feeding type. Selective pressure, together with mutational bias, was the dominant factor of shaping the codon usage patterns. The relaxed selective constraint, potentially with genome-based GC conversion (gBGC) and preferential codon usage, drive the fluctuation of GC contents among taxa with different lifestyles. Phylogenetic analysis in Orobanchaceae supported that parasitic species were single-originated while holoparasites were multiple-originated. Overall, the comparison of plastomes provided a good opportunity to understand the evolution process in Orobanchaceae with different lifestyles.

Correlation between gene expression levels under drought stress and synonymous codon usage in rice plant by in-silico study

We studied the correlation of synonymous codon usage (SCU) on gene expression levels under drought stress in rice. Sixty genes related to drought stress (with high, intermediate and low expression) were selected from rice meta-analysis data and various codon usage indices such as the effective number of codon usage (ENC), codon adaptation index (CAI) and relative synonymous codon usage (RSCU) were calculated. We found that in genes highly expressing under drought 1) GC content was higher, 2) ENC value was lower, 3) the preferred codons of some amino acids changed and 4) the RSCU ratio of GC-end codons relative to AT-end codons for 18 amino acids increased significantly compared with those in other genes. We introduce ARSCU as the Average ratio of RSCUs of GC-end codons to AT-end codons in each gene that could significantly separate high-expression genes under drought from low-expression genes. ARSCU is calculated using the program ARSCU-Calculator developed by our group to help predicting expression level of rice genes under drought. An index above ARSCU threshold is expected to indicate that the gene under study may belong to the "high expression group under drought". This information may be applied for codon optimization of genes for rice genetic engineering. To validate these findings, we further used 60 other genes (randomly selected subset of 43233 genes studied for their response to drought stress). ARSCU value was able to predict the level of expression at 88.33% of the cases. Using third set of 60 genes selected amongst high expressing genes not related to drought, only 31.65% of the genes showed ARSCU value of higher than the set threshold. This indicates that the phenomenon we described in this report may be unique for drought related genes. To justify the observed correlation between CUB and high expressing genes under drought, possible role of tRNA post transcriptional modification and tRFs was hypothesized as possible underlying biological mechanism.

Comparative analysis of codon usage between Gossypium hirsutum and G. barbadense mitochondrial genomes

Gossypium hirsutum and G. barbadense mitochondrial genomes were analyzed to understand the factors shaping codon usage. While most analyses of codon usage suggest minimal to no bias, nucleotide composition, specifically GC content, was significantly correlated with codon usage. In general, both mitochondrial genomes favor codons that end in A or U, with a secondary preference for pyrimidine rich codons. These observations are similar to previous reports of codon usage in cotton nuclear genomes, possibly suggestive of a general bias spanning genomic compartment. Although evidence for codon usage bias is weak for most genes, we identified six genes (i.e. atp8, atp9, sdh3, sdh4, mttB and rpl2) with significant nonrandom codon usage. In general, we find multiple factors that influence cotton mitochondrial genome codon usage, which may include selection in a subset of genes.

Analysis of codon usage patterns in "Lonicerae Flos" (Lonicera macranthoides Hand. -Mazz.) based on transcriptome data

Lonicera macranthoides Hand. -Mazz. is an important medicinal and economical plant in China, however, the codon usage bias (CUB) in L. macranthoides genes is still unknown. In this study, L. macranthoides transcriptome sequencing has been completed, and codon usage patterns in 36,090 reconstructed genes from the L. macranthoides transcriptome were examined. The mean GC content and GC3 value is 44.9% and 43.1%, respectively, which indicates that nucleotide contents of L. macranthoides genome is somewhat AT rich, and its codon bias pattern tends to use A/T-ending codons. According to neutrality plot, ENC plot, PR2-Bias plot and correspondence analysis, we know that both compositional constraint under selection and mutation could affect the CUB in L. macranthoides, and the mutation is the most determinant factor. Meanwhile, gene expression levels can influence its codon usage pattern. Furthermore, we identified 29 optimal codons and most of them ended with A/U. The study will lay a foundation for future research on gene prediction, genetic engineering and molecular evolution in L. macranthoides.

A cross-talk on compositional dynamics and codon usage patterns of mitochondrial CYB gene in Echinodermata

Codon usage bias (CUB) refers to a phenomenon in which some synonymous codons are used in mature mRNA at a higher frequency than other members codifying the same amino acid. CUB is mainly determined by mutation pressure and natural selection. We used bioinformatic tools to analyze the protein coding sequences of mitochondrial CYB gene in different classes of Echinodermata to understand the patterns of codon usage. The ENC values of CYB gene in five different classes of Echinodermata were 41.64, 30.33, 43.63, 41.11, and 41.33, which suggested that the CUB of this gene was low. The relative synonymous codon usage (RSCU) values showed that the patterns of over-represented and under-represented codons were different among different classes. Correspondence analysis indicated that the plots of CYB gene were different across classes, suggesting that the pattern of codon usage was also different among five classes under study. Highly significant correlation (p < .01) between overall nucleotide composition and its 3rd codon position indicated that both mutational pressure and natural selection had an influence on the codon usage bias of CYB gene. Furthermore, PR-2 bias plot analysis showed that both mutation pressure and natural selection might have affected the pattern of codon usage in CYB gene of Echinodermata.

Differences in Codon Usage Bias between Photosynthesis-Related Genes and Genetic System-Related Genes of Chloroplast Genomes in Cultivated and Wild Solanum Species

Solanum is one of the largest genera, including two important crops-potato (Solanum tuberosum) and tomato (Solanum lycopersicum). In this study we compared the chloroplast codon usage bias (CUB) among 12 Solanum species, between photosynthesis-related genes (Photo-genes) and genetic system-related genes (Genet-genes), and between cultivated species and wild relatives. The Photo-genes encode proteins for photosystems, the photosynthetic electron transport chain, and RuBisCO, while the Genet-genes encode proteins for ribosomal subunits, RNA polymerases, and maturases. The following findings about the Solanum chloroplast genome CUB were obtained: (1) the nucleotide composition, gene expression, and selective pressure are identified as the main factors affecting chloroplast CUB; (2) all these 12 chloroplast genomes prefer A/U over G/C and pyrimidines over purines at the third-base of codons; (3) Photo-genes have higher codon adaptation indexes than Genet-genes, indicative of a higher gene expression level and a stronger adaptation of Photo-genes; (4) gene function is the primary factor affecting CUB of Photo-genes but not Genet-genes; (5) Photo-genes prefer pyrimidine over purine, whereas Genet-genes favor purine over pyrimidine, at the third position of codons; (6) Photo-genes are mainly affected by the selective pressure, whereas Genet-genes are under the underlying mutational bias; (7) S. tuberosum is more similar with Solanum commersonii than with Solanum bulbocastanum; (8) S. lycopersicum is greatly different from the analyzed seven wild relatives; (9) the CUB in codons for valine, aspartic acid, and threonine are the same between the two crop species, S. tuberosum and S. lycopersicum. These findings suggest that the chloroplast CUB contributed to the differential requirement of gene expression activity and function between Photo-genes and Genet-genes and to the performance of cultivated potato and tomato.

Association of EPHA3 Gene Polymorphisms with Nonsyndromic Cleft Lip With or Without Cleft Palate

AIMS

Nonsyndromic cleft lip with or without palate (NSCL/P) represents a complex condition caused by genetic and environmental factors. The aim of this study was to investigate the relationship between the EPHA3 polymorphisms and NSCL/P.

MATERIALS AND METHODS

To investigate the relationship between five EPHA3 single nucleotide polymorphisms (SNPs) and NSCL/P, we selected 180 affected patients and 167 normal controls from the Chinese Han Population. EPHA3 SNPs (rs7650466, rs1398197, rs17801309, rs1054750, and rs7632427) were genotyped using the SNaPshot technique; bioinformatic analyses were performed to determine if any of them were potentially functional SNPs.

RESULTS

The rs7650466 T allele was associated with the incidence of NSCL/P (OR, 0.211; 95% CI, 0.131-0.338; adjusted p = 4.881 × 10^-10) and cleft lip with or without palate (CL/P) (OR, 0.176; 95% CI, 0.104-0.297; adjusted p = 3.617 × 10^-10), as well as with protective and dominant effects in both conditions. The rs7650466 T allele could be associated with reduced risk of the malformation. In a bioinformatics analysis, we found potential matching sites (miR-1255a, miR-125a-3p, miR-143, and miR-552) for rs7650466 and preliminarily analyzed its potential function.

CONCLUSIONS

Collectively, our data suggest that the EPHA3 rs7650466 polymorphism confers genetic risk for NSCL/P in the Chinese Han Population. Furthermore, rs7650466 is associated with CL/P incidence in stratified analyses, but not with cleft palate only.

Compositional bias coupled with selection and mutation pressure drives codon usage in Brassica campestris genes

The plant Brassica campestris includes the vegetables turnip and Chinese cabbage, important plants of economic importance. Here, we have analysed the codon usage bias of B. campestris for 116 protein coding genes. Neutrality analysis showed that B. campestris had a wide range of GC3s, and a significant correlation was observed between GC12 and GC3. Nc versus GC3s plot showed a few genes on or proximate to the expected curve, but the majority of points were found to be scattered distantly from the expected curve. Correspondence analysis on codon usage revealed that the position preference of codons on multidimensional space totally depends on the presence of A and T at synonymous third codon position. These results altogether suggest that composition bias along with selection (major) and mutation pressure (minor) affects the codon usage pattern of the protein coding genes in Brassica campestris.

Analysis of synonymous codon usage bias and phylogeny of coat protein gene in banana bract mosaic virus isolates

The banana is one of the world's most important livelihood crops. Banana plants are principally infected by four virus species, Banana bunchy top virus (genus Babuvirus), Cucumber mosaic virus (genus Cucumovirus), Banana streak virus (genus Badnavirus) and Banana bract mosaic virus (genus Potyvirus). The objective of this study is to understand the codon usage pattern and phylogeny of coat protein gene in different banana bract mosaic virus (BBrMV) isolates. The BBrMV Coat Protein (CP) gene was amplified from BBrMV infected banana plant samples collected from different districts of Tamil Nadu and Karnataka, India. Six new BBrMV isolates were submitted to National Center for Biotechnology Information. Phylogenetic analysis and codon usage indices were studied along with other isolates of BBrMV. Phylogenetic analysis of CP genes shows that most of BBrMV isolates are closely related to each other except KF385484.1 and KF385478.1. Relative codon usage patterns among different BBrMV isolates were calculated by software CodonW version 1.4.2. In BBrMV, codons with A-ended or U ended are the most preferential except the Leu and Gln whose optimized codons are CAG and UUG ending by G. The codon usage patterns of BBrMV isolates are principally influenced by mutational bias; however, compositional constraints along with mutational bias also play a major role.

Analysis of transcriptome data reveals multifactor constraint on codon usage in Taenia multiceps

Background

Codon usage bias (CUB) is an important evolutionary feature in genomes that has been widely observed in many organisms. However, the synonymous codon usage pattern in the genome of T. multiceps remains to be clarified. In this study, we analyzed the codon usage of T. multiceps based on the transcriptome data to reveal the constraint factors and to gain an improved understanding of the mechanisms that shape synonymous CUB.

Results

Analysis of a total of 8,620 annotated mRNA sequences from T. multiceps indicated only a weak codon bias, with mean GC and GC3 content values of 49.29% and 51.43%, respectively. Our analysis indicated that nucleotide composition, mutational pressure, natural selection, gene expression level, amino acids with grand average of hydropathicity (GRAVY) and aromaticity (Aromo) and the effective selection of amino-acids all contributed to the codon usage in T. multiceps. Among these factors, natural selection was implicated as the major factor affecting the codon usage variation in T. multiceps. The codon usage of ribosome genes was affected mainly by mutations, while the essential genes were affected mainly by selection. In addition, 21codons were identified as “optimal codons”. Overall, the optimal codons were GC-rich (GC:AU, 41:22), and ended with G or C (except CGU). Furthermore, different degrees of variation in codon usage were found between T. multiceps and Escherichia coli, yeast, Homo sapiens. However, little difference was found between T. multiceps and Taenia pisiformis.

Conclusions

In this study, the codon usage pattern of T. multiceps was analyzed systematically and factors affected CUB were also identified. This is the first study of codon biology in T. multiceps. Understanding the codon usage pattern in T. multiceps can be helpful for the discovery of new genes, molecular genetic engineering and evolutionary studies.

Electronic supplementary material

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

Biased Nucleotide Composition and Differential Codon Usage Pattern in HIV-1 and HIV-2

HIV-1 and HIV-2 are closely related retroviruses with differences in pathogenicity and geographic distribution. HIV-2 infection is associated with slower disease progression and transmission, longer latency period, low or undetectable plasmatic viral loads, and reduced likelihood of progression to AIDS, compared to HIV-1. In this investigation, we analyzed HIV-2 genes and genomes and compared them with that of HIV-1 belonging to various subtypes. Comparative analysis of the effective number of codons (ENC) for each of the nine genes of the two viruses revealed that the tat gene of HIV-2 had a higher ENC value compared to HIV-1 tat, reflecting lower levels of expression of HIV-2 tat. Lower levels of tat protein particularly during the early stages of infection could result in a lower viral load, lower viral set point, and delayed progression of disease in HIV-2-infected individuals compared to HIV-1-infected subjects. Furthermore, the GC3 composition of the regulatory genes of HIV-2 was ≥50%, suggesting a firm effort by these viruses to adapt themselves to evolutionary survival. We hypothesize that differential codon usage could be one of the possible factors that could contribute to the diminished pathogenicity of HIV-2 in the host as compared to HIV-1.

Genome-Wide Analysis of Codon Usage Bias in Epichloë festucae

Analysis of codon usage data has both practical and theoretical applications in understanding the basics of molecular biology. Differences in codon usage patterns among genes reflect variations in local base compositional biases and the intensity of natural selection. Recently, there have been several reports related to codon usage in fungi, but little is known about codon usage bias in Epichloë endophytes. The present study aimed to assess codon usage patterns and biases in 4870 sequences from Epichloë festucae, which may be helpful in revealing the constraint factors such as mutation or selection pressure and improving the bioreactor on the cloning, expression, and characterization of some special genes. The GC content with 56.41% is higher than the AT content (43.59%) in E. festucae. The results of neutrality and effective number of codons plot analyses showed that both mutational bias and natural selection play roles in shaping codon usage in this species. We found that gene length is strongly correlated with codon usage and may contribute to the codon usage patterns observed in genes. Nucleotide composition and gene expression levels also shape codon usage bias in E. festucae. E. festucae exhibits codon usage bias based on the relative synonymous codon usage (RSCU) values of 61 sense codons, with 25 codons showing an RSCU larger than 1. In addition, we identified 27 optimal codons that end in a G or C.

Molecular Evolution and Functional Diversification of Replication Protein A1 in Plants

Replication protein A (RPA) is a heterotrimeric, single-stranded DNA binding complex required for eukaryotic DNA replication, repair, and recombination. RPA is composed of three subunits, RPA1, RPA2, and RPA3. In contrast to single RPA subunit genes generally found in animals and yeast, plants encode multiple paralogs of RPA subunits, suggesting subfunctionalization. Genetic analysis demonstrates that five Arabidopsis thaliana RPA1 paralogs (RPA1A to RPA1E) have unique and overlapping functions in DNA replication, repair, and meiosis. We hypothesize here that RPA1 subfunctionalities will be reflected in major structural and sequence differences among the paralogs. To address this, we analyzed amino acid and nucleotide sequences of RPA1 paralogs from 25 complete genomes representing a wide spectrum of plants and unicellular green algae. We find here that the plant RPA1 gene family is divided into three general groups termed RPA1A, RPA1B, and RPA1C, which likely arose from two progenitor groups in unicellular green algae. In the family Brassicaceae the RPA1B and RPA1C groups have further expanded to include two unique sub-functional paralogs RPA1D and RPA1E, respectively. In addition, RPA1 groups have unique domains, motifs, cis-elements, gene expression profiles, and pattern of conservation that are consistent with proposed functions in monocot and dicot species, including a novel C-terminal zinc-finger domain found only in plant RPA1C-like sequences. These results allow for improved prediction of RPA1 subunit functions in newly sequenced plant genomes, and potentially provide a unique molecular tool to improve classification of Brassicaceae species.

Gene Evolutionary Trajectories and GC Patterns Driven by Recombination in Zea mays

Recombination occurring during meiosis is critical for creating genetic variation and plays an essential role in plant evolution. In addition to creating novel gene combinations, recombination can affect genome structure through altering GC patterns. In maize (Zea mays) and other grasses, another intriguing GC pattern exists. Maize genes show a bimodal GC content distribution that has been attributed to nucleotide bias in the third, or wobble, position of the codon. Recombination may be an underlying driving force given that recombination sites are often associated with high GC content. Here we explore the relationship between recombination and genomic GC patterns by comparing GC gene content at each of the three codon positions (GC₁, GC₂, and GC₃, collectively termed GC_x) to instances of a variable GC-rich motif that underlies double strand break (DSB) hotspots and to meiocyte-specific gene expression. Surprisingly, GC_x bimodality in maize cannot be fully explained by the codon wobble hypothesis. High GC_x genes show a strong overlap with the DSB hotspot motif, possibly providing a mechanism for the high evolutionary rates seen in these genes. On the other hand, genes that are turned on in meiosis (early prophase I) are biased against both high GC_x genes and genes with the DSB hotspot motif, possibly allowing important meiotic genes to avoid DSBs. Our data suggests a strong link between the GC-rich motif underlying DSB hotspots and high GC_x genes.

GC-Content of Synonymous Codons Profoundly Influences Amino Acid Usage

Amino acids typically are encoded by multiple synonymous codons that are not used with the same frequency. Codon usage bias has drawn considerable attention, and several explanations have been offered, including variation in GC-content between species. Focusing on a simple parameter—combined GC proportion of all the synonymous codons for a particular amino acid, termed GC_syn—we try to deepen our understanding of the relationship between GC-content and amino acid/codon usage in more details. We analyzed 65 widely distributed representative species and found a close association between GC_syn, GC-content, and amino acids usage. The overall usages of the four amino acids with the greatest GC_syn and the five amino acids with the lowest GC_syn both vary with the regional GC-content, whereas the usage of the remaining 11 amino acids with intermediate GC_syn is less variable. More interesting, we discovered that codon usage frequencies are nearly constant in regions with similar GC-content. We further quantified the effects of regional GC-content variation (low to high) on amino acid usage and found that GC-content determines the usage variation of amino acids, especially those with extremely high GC_syn, which accounts for 76.7% of the changed GC-content for those regions. Our results suggest that GC_syn correlates with GC-content and has impact on codon/amino acid usage. These findings suggest a novel approach to understanding the role of codon and amino acid usage in shaping genomic architecture and evolutionary patterns of organisms.

Comparative Analysis of Codon Usage Bias Patterns in Microsporidian Genomes

The sub-3 Mbp genomes from microsporidian species of the Encephalitozoon genus are the smallest known among eukaryotes and paragons of genomic reduction and compaction in parasites. However, their diminutive stature is not characteristic of all Microsporidia, whose genome sizes vary by an order of magnitude. This large variability suggests that different evolutionary forces are applied on the group as a whole. In this study, we have compared the codon usage bias (CUB) between eight taxonomically distinct microsporidian genomes: Encephalitozoon intestinalis, Encephalitozoon cuniculi, Spraguea lophii, Trachipleistophora hominis, Enterocytozoon bieneusi, Nematocida parisii, Nosema bombycis and Nosema ceranae. While the CUB was found to be weak in all eight Microsporidia, nearly all (98%) of the optimal codons in S. lophii, T. hominis, E. bieneusi, N. parisii, N. bombycis and N. ceranae are fond of A/U in third position whereas most (64.6%) optimal codons in the Encephalitozoon species E. intestinalis and E. cuniculi are biased towards G/C. Although nucleotide composition biases are likely the main factor driving the CUB in Microsporidia according to correlation analyses, directed mutational pressure also likely affects the CUB as suggested by ENc-plots, correspondence and neutrality analyses. Overall, the Encephalitozoon genomes were found to be markedly different from the other microsporidians and, despite being the first sequenced representatives of this lineage, are uncharacteristic of the group as a whole. The disparities observed cannot be attributed solely to differences in host specificity and we hypothesize that other forces are at play in the lineage leading to Encephalitozoon species.

Mutational bias is the driving force for shaping the synonymous codon usage pattern of alternatively spliced genes in rice (Oryza sativa L.)

Alternative splicing plays important roles in diverse aspects of plant development, metabolism, and stress responses. However, the regulatory mechanisms of alternative splicing of genes still remain incompletely elucidated, especially in plants. In this study, the synonymous codon usage pattern of alternatively spliced (AS) genes in rice was firstly explored using the combination of correspondence analysis (CA), internal CA, correlation and ANOVA analyses. The results show that alternatively and non-alternatively spliced (non-AS) genes have similar tendency for overall codon usage, but exhibit significant difference in 58 out of 64 codons. AS and non-AS genes are both under strong purifying selection, but the former ones have significant lower mutation rate and are prone to be enriched towards the chromosomal ends. In the group of AS genes, the variability in synonymous codon usage between genes is mainly due to the variations in GC content, CDS length, as well as gene functions. Mutational bias that accounts for 25.85 % of the total codon usage variability plays a major role in shaping the codon usage pattern of AS genes. In contrast, no obvious evidence is found for the contributions of translational selection, AS types, the conservation of AS events, and numbers of AS variants to the codon usage divergence between AS genes. These findings may be useful for further understanding the mechanisms of origination, differentiation and regulation of alternatively spliced genes in plants.

Improved expression of recombinant plant-made hEGF

KEY MESSAGE

The yield of recombinant hEGF was increased approximately tenfold through a range of optimisations. Further, the recombinant protein was found to have biological activity comparable to commercial hEGF. Human epidermal growth factor (hEGF) is a powerful mitogen that can enhance the healing of a wide range of injuries, including burns, cuts, diabetic ulcers and gastric ulcers. However, despite its clinical value, hEGF is only consistently used for the treatment of chronic diabetic ulcers due to its high cost. In this study, hEGF was transiently expressed in Nicotiana benthamiana plants and targeted to the apoplast, ER and vacuole. Several other approaches were also included in a stepwise fashion to identify the optimal conditions for the expression of recombinant hEGF. Expression was found to be highest in the vacuole, while targeting hEGF to the ER caused a decrease in total soluble protein (TSP). Using a codon optimised sequence was found to increase vacuolar targeted hEGF yield by ~34 %, while it was unable to increase the yield of ER targeted hEGF. The use of the P19 silencing inhibitor was able to further increase expression by over threefold, and using 5-week-old plants significantly increased expression compared to 4- or 6-week-old-plants. The combined effect of these optimisations increased expression tenfold over the initial apoplast targeted construct to an average yield of 6.24 % of TSP. The plant-made hEGF was then shown to be equivalent to commercial E. coli derived hEGF in its ability to promote the proliferation of mouse keratinocytes. This study supports the potential for plants to be used for the commercial production of hEGF, and identifies a potential limitation for the further improvement of recombinant protein yields.

Function and localization of the Arabidopsis thaliana diacylglycerol acyltransferase DGAT2 expressed in yeast

Diacylglycerol acyltransferases (DGATs) catalyze the final and only committed step of triacylglycerol synthesis. DGAT activity is rate limiting for triacylglycerol accumulation in mammals, plants and microbes. DGATs belong to three different evolutionary classes. In Arabidopsis thaliana, DGAT1, encoded by At2g19450, is the major DGAT enzyme involved in triacylglycerol accumulation in seeds. Until recently, the function of DGAT2 (At3g51520) has remained elusive. Previous attempts to characterize its enzymatic function by heterologous expression in yeast were unsuccessful. In the present report we demonstrate that expression of a codon-optimized version of the DGAT2 gene is able to restore neutral lipid accumulation in the Saccharomyces cerevisiae mutant strain (H1246), which is defective in triacylglycerol biosynthesis. Heterologous expression of codon-optimized DGAT2 and DGAT1 induced the biogenesis of subcellular lipid droplets containing triacylglycerols and squalene. Both DGAT proteins were found to be associated with these lipid droplets. The fatty acid composition was affected by the nature of the acyltransferase expressed. DGAT2 preferentially incorporated C16:1 fatty acids whereas DGAT1 displayed preference for C16:0, strongly suggesting that these enzymes have contrasting substrate specificities.

Mutation and selection cause codon usage and bias in mitochondrial genomes of ribbon worms (Nemertea)

The phenomenon of codon usage bias is known to exist in many genomes and it is mainly determined by mutation and selection. To understand the patterns of codon usage in nemertean mitochondrial genomes, we use bioinformatic approaches to analyze the protein-coding sequences of eight nemertean species. Neutrality analysis did not find a significant correlation between GC12 and GC3. ENc-plot showed a few genes on or close to the expected curve, but the majority of points with low-ENc values are below it. ENc-plot suggested that mutational bias plays a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that GC and AT were not used proportionally and we propose that codons containing A or U at third position are used preferentially in nemertean species, regardless of whether corresponding tRNAs are encoded in the mitochondrial DNA. Context-dependent analysis indicated that the nucleotide at the second codon position slightly affects synonymous codon choices. These results suggested that mutational and selection forces are probably acting to codon usage bias in nemertean mitochondrial genomes.

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.

In Arabidopsis thaliana codon volatility scores reflect GC3 composition rather than selective pressure

BACKGROUND

Synonymous codon usage bias has typically been correlated with, and attributed to translational efficiency. However, there are other pressures on genomic sequence composition that can affect codon usage patterns such as mutational biases. This study provides an analysis of the codon usage patterns in Arabidopsis thaliana in relation to gene expression levels, codon volatility, mutational biases and selective pressures.

RESULTS

We have performed synonymous codon usage and codon volatility analyses for all genes in the A. thaliana genome. In contrast to reports for species from other kingdoms, we find that neither codon usage nor volatility are correlated with selection pressure (as measured by dN/dS), nor with gene expression levels on a genome wide level. Our results show that codon volatility and usage are not synonymous, rather that they are correlated with the abundance of G and C at the third codon position (GC3).

CONCLUSIONS

Our results indicate that while the A. thaliana genome shows evidence for synonymous codon usage bias, this is not related to the expression levels of its constituent genes. Neither codon volatility nor codon usage are correlated with expression levels or selective pressures but, because they are directly related to the composition of G and C at the third codon position, they are the result of mutational bias. Therefore, in A. thaliana codon volatility and usage do not result from selection for translation efficiency or protein functional shift as measured by positive selection.

Acyl-CoA:diacylglycerol acyltransferase: molecular biology, biochemistry and biotechnology

Triacylglycerol (TG) is a storage lipid which serves as an energy reservoir and a source of signalling molecules and substrates for membrane biogenesis. TG is essential for many physiological processes and its metabolism is widely conserved in nature. Acyl-CoA:diacylglycerol acyltransferase (DGAT, EC 2.3.1.20) catalyzes the final step in the sn-glycerol-3-phosphate pathway leading to TG. DGAT activity resides mainly in two distinct membrane bound polypeptides, known as DGAT1 and DGAT2 which have been identified in numerous organisms. In addition, a few other enzymes also hold DGAT activity, including the DGAT-related acyl-CoA:monoacylglycerol acyltransferases (MGAT). Progress on understanding structure/function in DGATs has been limited by the lack of detailed three-dimensional structural information due to the hydrophobic properties of theses enzymes and difficulties associated with purification. This review examines several aspects of DGAT and MGAT genes and enzymes, including current knowledge on their gene structure, expression pattern, biochemical properties, membrane topology, functional motifs and subcellular localization. Recent progress in probing structural and functional aspects of DGAT1 and DGAT2, using a combination of molecular and biochemical techniques, is emphasized. Biotechnological applications involving DGAT enzymes ranging from obesity therapeutics to oilseed engineering are also discussed.

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.

Patterns and influencing factor of synonymous codon usage in porcine circovirus

Background

Analysis of codon usage can reveal much about the molecular evolution of the viruses. Nevertheless, little information about synonymous codon usage pattern of porcine circovirus (PCV) genome in the process of its evolution is available. In this study, to give a new understanding on the evolutionary characteristics of PCV and the effects of natural selection from its host on the codon usage pattern of the virus, Patterns and the key determinants of codon usage in PCV were examined.

Methods

We carried out comprehensive analysis on codon usage pattern in the PCV genome, by calculating relative synonymous codon usage (RSCU), effective number of codons (ENC), dinucleotides and nucleic acid content of the PCV genome.

Results

PCV genomes have relatively much lower content of GC and codon preference, this result shows that nucleotide constraints have a major impact on its synonymous codon usage. The results of the correspondence analysis indicate codon usage patterns of PCV of various genotypes, various subgenotypes changed greatly, and significant differences in codon usage patterns of Each virus of Circoviridae.There is much comparability between PCV and its host in their synonymous codon usage, suggesting that the natural selection pressure from the host factor also affect the codon usage patterns of PCV. In particular, PCV genotype II is in synonymous codon usage more similar to pig than to PCV genotype I, which may be one of the most important molecular mechanisms of PCV genotype II to cause disease. The calculations results of the relative abundance of dinucleotides indicate that the composition of dinucleotides also plays a key role in the variation found in synonymous codon usage in PCV. Furthermore, geographic factors, the general average hydrophobicity and the aromaticity may be related to the formation of codon usage patterns of PCV.

Conclusion

The results of these studies suggest that synonymous codon usage pattern of PCV genome are the result of interaction between mutation pressure and natural selection from its host. The information from this study may not only have theoretical value in understanding the characteristics of synonymous codon usage in PCV genomes, but also have significant value for the molecular evolution of PCV.

Rates of evolution in stress-related genes are associated with habitat preference in two Cardamine lineages

Background

Elucidating the selective and neutral forces underlying molecular evolution is fundamental to understanding the genetic basis of adaptation. Plants have evolved a suite of adaptive responses to cope with variable environmental conditions, but relatively little is known about which genes are involved in such responses. Here we studied molecular evolution on a genome-wide scale in two species of Cardamine with distinct habitat preferences: C. resedifolia, found at high altitudes, and C. impatiens, found at low altitudes. Our analyses focussed on genes that are involved in stress responses to two factors that differentiate the high- and low-altitude habitats, namely temperature and irradiation.

Results

High-throughput sequencing was used to obtain gene sequences from C. resedifolia and C. impatiens. Using the available A. thaliana gene sequences and annotation, we identified nearly 3,000 triplets of putative orthologues, including genes involved in cold response, photosynthesis or in general stress responses. By comparing estimated rates of molecular substitution, codon usage, and gene expression in these species with those of Arabidopsis, we were able to evaluate the role of positive and relaxed selection in driving the evolution of Cardamine genes. Our analyses revealed a statistically significant higher rate of molecular substitution in C. resedifolia than in C. impatiens, compatible with more efficient positive selection in the former. Conversely, the genome-wide level of selective pressure is compatible with more relaxed selection in C. impatiens. Moreover, levels of selective pressure were heterogeneous between functional classes and between species, with cold responsive genes evolving particularly fast in C. resedifolia, but not in C. impatiens.

Conclusions

Overall, our comparative genomic analyses revealed that differences in effective population size might contribute to the differences in the rate of protein evolution and in the levels of selective pressure between the C. impatiens and C. resedifolia lineages. The within-species analyses also revealed evolutionary patterns associated with habitat preference of two Cardamine species. We conclude that the selective pressures associated with the habitats typical of C. resedifolia may have caused the rapid evolution of genes involved in cold response.

Interspecific and intragenic differences in codon usage bias among vertebrate myosin heavy-chain genes

Synonymous codon usage bias is a broadly observed phenomenon in bacteria, plants, and invertebrates and may result from selection. However, the role of selective pressures in shaping codon bias is still controversial in vertebrates, particularly for mammals. The myosin heavy-chain (MyHC) gene family comprises multiple isoforms of the major force-producing contractile protein in cardiac and skeletal muscles. Slow and fast genes are tandemly arrayed on separate chromosomes, and have distinct patterns of functionality and expression in muscle. We analyze both full-length MyHC genes (~5400 bp) and a larger collection of partial sequences at the 3' end (~500 bp). The MyHC isoforms are an interesting system in which to study codon usage bias because of their length, expression, and critical importance to organismal mobility. Codon bias and GC content differs among MyHC genes with regards to functional type, isoform, and position within the gene. Codon bias even varies by isoform within a species. We find evidence in favor of both chromosomal influences on nucleotide composition and selection against nonsense errors (SANE) acting on codon usage in MyHC genes. Intragenic variation in codon bias and elongation rate is significant, with a strong trend for increasing codon bias and elongation rate towards the 3' end of the gene, although the trend is dependent upon the degeneracy class of the codons. Therefore, patterns of codon usage in MyHC genes are consistent with models supporting SANE as a major force shaping codon usage.

Reduced efficacy of natural selection on codon usage bias in selfing Arabidopsis and Capsella species

Population genetic theory predicts that the efficacy of natural selection in a self-fertilizing species should be lower than its outcrossing relatives because of the reduction in the effective population size (N(e)) in the former brought about by inbreeding. However, previous analyses comparing Arabidopsis thaliana (selfer) with A. lyrata (outcrosser) have not found conclusive support for this prediction. In this study, we addressed this issue by examining silent site polymorphisms (synonymous and intronic), which are expected to be informative about changes in N(e). Two comparisons were made: A. thaliana versus A. lyrata and Capsella rubella (selfer) versus C. grandiflora (outcrosser). Extensive polymorphism data sets were obtained by compiling published data from the literature and by sequencing 354 exon loci in C. rubella and 89 additional loci in C. grandiflora. To extract information from the data effectively for studying these questions, we extended two recently developed models in order to investigate detailed selective differences between synonymous codons, mutational biases, and biased gene conversion (BGC), taking into account the effects of recent changes in population size. We found evidence that selection on synonymous codons is significantly weaker in the selfers compared with the outcrossers and that this difference cannot be fully accounted for by mutational biases or BGC.

Codon usage: nature's roadmap to expression and folding of proteins

Biomedical and biotechnological research relies on processes leading to the successful expression and production of key biological products. High-quality proteins are required for many purposes, including protein structural and functional studies. Protein expression is the culmination of multistep processes involving regulation at the level of transcription, mRNA turnover, protein translation, and post-translational modifications leading to the formation of a stable product. Although significant strides have been achieved over the past decade, advances toward integrating genomic and proteomic information are essential, and until such time, many target genes and their products may not be fully realized. Thus, the focus of this review is to provide some experimental support and a brief overview of how codon usage bias has evolved relative to regulating gene expression levels.

Complex mutation and weak selection together determined the codon usage bias in bryophyte mitochondrial genomes

Mutation and selection are two major forces causing codon usage biases. How these two forces influence the codon usages in green plant mitochondrial genomes has not been well investigated. In the present study, we surveyed five bryophyte mitochondrial genomes to reveal their codon usage patterns as well as the determining forces. Three interesting findings were made. First, comparing to Chara vulgaris, an algal species sister to all extant land plants, bryophytes have more G, C-ending codon usages in their mitochondrial genes. This is consistent with the generally higher genomic GC content in bryophyte mitochondria, suggesting an increased mutational pressure toward GC. Second, as indicated by Wright's Nc-GC3s plot, mutation, not selection, is the major force affecting codon usages of bryophyte mitochondrial genes. However, the real mutational dynamics seem very complex. Context-dependent analysis indicated that nucleotide at the 2nd codon position would slightly affect synonymous codon choices. Finally, in bryophyte mitochondria, tRNA genes would apply a weak selection force to fine-tune the synonymous codon frequencies, as revealed by data of Ser4-Pro-Thr-Val families. In summary, complex mutation and weak selection together determined the codon usages in bryophyte mitochondrial genomes.

GC3 biology in corn, rice, sorghum and other grasses

BACKGROUND

The third, or wobble, position in a codon provides a high degree of possible degeneracy and is an elegant fault-tolerance mechanism. Nucleotide biases between organisms at the wobble position have been documented and correlated with the abundances of the complementary tRNAs. We and others have noticed a bias for cytosine and guanine at the third position in a subset of transcripts within a single organism. The bias is present in some plant species and warm-blooded vertebrates but not in all plants, or in invertebrates or cold-blooded vertebrates.

RESULTS

Here we demonstrate that in certain organisms the amount of GC at the wobble position (GC3) can be used to distinguish two classes of genes. We highlight the following features of genes with high GC3 content: they (1) provide more targets for methylation, (2) exhibit more variable expression, (3) more frequently possess upstream TATA boxes, (4) are predominant in certain classes of genes (e.g., stress responsive genes) and (5) have a GC3 content that increases from 5'to 3'. These observations led us to formulate a hypothesis to explain GC3 bimodality in grasses.

CONCLUSIONS

Our findings suggest that high levels of GC3 typify a class of genes whose expression is regulated through DNA methylation or are a legacy of accelerated evolution through gene conversion. We discuss the three most probable explanations for GC3 bimodality: biased gene conversion, transcriptional and translational advantage and gene methylation.

Fine-tuning translation kinetics selection as the driving force of codon usage bias in the hepatitis A virus capsid

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness.

Each organism has a specific codon usage signature. Translational selection i.e., selection for the codon adaptation to the tRNA pools, is one of the driving forces of codon bias. In the virus world, this implies an adjustment of the virus codon usage to that of the host cell. Hepatitis A virus appears as an exception to the rule, with a highly deoptimized codon usage, suggesting that translational selection is not the underlying mechanism of its codon bias. However, since the virus lacks a mechanism of cellular protein synthesis inhibition, the deoptimized codon usage may be envisaged as a hawk (cell) and dove (hepatitis A virus) competition strategy for tRNAs and translational selection as well. To confirm this possibility, we artificially induced cell protein synthesis shut-off, thus increasing the tRNA pool availability for the virus, and we took advantage of the quasispecies dynamics to elucidate changes in its codon usage. Virus adaptation to the drug results in a re-deoptimization of codon usage in the capsid region, suggesting a requirement of a slow translation rate, i.e., a translation kinetic selection, instead of a translational selection associated with an optimization of the codon usage. Translation kinetics control is based on the right combination of codons (common and rare) that allows a regulated ribosome traffic rate ensuring the proper protein folding. Capsid folding is critical for a virus transmitted through the fecal-oral route with long extracorporeal periods.

Codon usage is associated with the evolutionary age of genes in metazoan genomes

Background

Codon usage may vary significantly between different organisms and between genes within the same organism. Several evolutionary processes have been postulated to be the predominant determinants of codon usage: selection, mutation, and genetic drift. However, the relative contribution of each of these factors in different species remains debatable. The availability of complete genomes for tens of multicellular organisms provides an opportunity to inspect the relationship between codon usage and the evolutionary age of genes.

Results

We assign an evolutionary age to a gene based on the relative positions of its identified homologues in a standard phylogenetic tree. This yields a classification of all genes in a genome to several evolutionary age classes. The present study starts from the observation that each age class of genes has a unique codon usage and proceeds to provide a quantitative analysis of the codon usage in these classes. This observation is made for the genomes of Homo sapiens, Mus musculus, and Drosophila melanogaster. It is even more remarkable that the differences between codon usages in different age groups exhibit similar and consistent behavior in various organisms. While we find that GC content and gene length are also associated with the evolutionary age of genes, they can provide only a partial explanation for the observed codon usage.

Conclusion

While factors such as GC content, mutational bias, and selection shape the codon usage in a genome, the evolutionary history of an organism over hundreds of millions of years is an overlooked property that is strongly linked to GC content, protein length, and, even more significantly, to the codon usage of metazoan genomes.

Analysis of synonymous codon usage patterns in different plant mitochondrial genomes

Codon usage in mitochondrial genome of the six different plants was analyzed to find general patterns of codon usage in plant mitochondrial genomes. The neutrality analysis indicated that the codon usage patterns of mitochondrial genes were more conserved in GC content and no correlation between GC12 and GC3. T and A ending codons were detected as the preferred codons in plant mitochondrial genomes. The Parity Rule 2 plot analysis showed that T was used more frequently than A. The EN(C)-plot showed that although a majority of the points with low EN(C) values were lying below the expected curve, a few genes lied on the expected curve. Correspondence analysis of relative synonymous codon usage yielded a first axis that explained only a partial amount of variation of codon usage. These findings suggest that natural selection is likely to be playing a large role in codon usage bias in plant mitochondrial genomes, but not only natural selection but also other several factors are likely to be involved in determining the selective constraints on codon bias in plant mitochondrial genomes. Meantime, 1 codon (P. patens), 6 codons (Z. mays), 9 codons (T. aestivum), 15 codons (A. thaliana), 15 codons (M. polymorpha) and 15 codons (N. tabacum) were defined as the preferred codons of the six plant mitochondrial genomes.

Ribosome and transcript copy numbers, polysome occupancy and enzyme dynamics in Arabidopsis

Plants are exposed to continual changes in the environment. The daily alternation between light and darkness results in massive recurring changes in the carbon budget, and leads to widespread changes in transcript levels. These diurnal changes are superimposed on slower changes in the environment. Quantitative molecular information about the numbers of ribosomes, of transcripts for 35 enzymes in central metabolism and their loading into polysomes is used to estimate translation rates in Arabidopsis rosettes, and explore the consequences for important sub-processes in plant growth. Translation rates for individual enzyme are compared with their abundance in the rosette to predict which enzymes are subject to rapid turnover every day, and which are synthesized at rates that would allow only slow adjustments to sustained changes of the environment, or resemble those needed to support the observed rate of growth. Global translation rates are used to estimate the energy costs of protein synthesis and relate them to the plant carbon budget, in particular the rates of starch degradation and respiration at night.