Variation in global codon usage bias among prokaryotic organisms is associated with their lifestyles

Background  

It is widely acknowledged that synonymous codons are used unevenly among genes in a genome. In organisms under translational selection, genes encoding highly expressed proteins are enriched with specific codons. This phenomenon, termed codon usage bias, is common to many organisms and has been recognized as influencing cellular fitness. This suggests that the global extent of codon usage bias of an organism might be associated with its phenotypic traits.     

Protein elemental sparing and codon usage bias are correlated among bacteria

Highly expressed proteins can exhibit relatively small material costs, in terms of the quantities of carbon (C), nitrogen (N) or sulphur (S) atoms they contain. This 'elemental sparing' probably reflects selection to reduce the quantities of potentially growth-limiting elements in abundant proteins, but the evolutionary mechanisms for adaptive elemental sparing are still poorly understood. Here, we predict that the extent of 'elemental sparing' in highly expressed proteins will vary among organisms, according to the effectiveness of selection in determining the fate of mutations. We test this hypothesis in bacteria by asking whether 'elemental sparing' is correlated with codon usage bias. Bacteria exhibit extraordinary variation in their life histories and demography and consequently in the effectiveness of selection in determining whether preferred codons are used in highly expressed genes. We find that C sparing and S sparing, but not N sparing, are significantly correlated with adaptive codon usage bias among 148 genera of bacteria, suggesting that selection for elemental sparing and codon bias are promoted by similar bacterial traits. Our study helps identify principles that determine how nutrient scarcity can shape the elemental composition of proteins.     

Gene expression,nucleotide composition and codon usage bias of genes associated with human Y chromosome

Analysis of codon usage pattern is important to understand the genetic and evolutionary characteristics of genomes. We have used bioinformatic approaches to analyze the codon usage bias (CUB) of the genes located in human Y chromosome. Codon bias index (CBI) indicated that the overall extent of codon usage bias was low. The relative synonymous codon usage (RSCU) analysis suggested that approximately half of the codons out of 59 synonymous codons were most frequently used, and possessed a T or G at the third codon position. The codon usage pattern was different in different genes as revealed from correspondence analysis (COA). A significant correlation between effective number of codons (ENC) and various GC contents suggests that both mutation pressure and natural selection affect the codon usage pattern of genes located in human Y chromosome. In addition, Y-linked genes have significant difference in GC contents at the second and third codon positions, expression level, and codon usage pattern of some codons like the SPANX genes in X chromosome.     

Analysis of synonymous codon usage bias in Chlamydia

Chlamydiae are obligate intracellular bacterial pathogens that cause ocular and sexuallytransmitted diseases,and are associated with cardiovascular diseases.The analysis of codon usage mayimprove our understanding of the evolution and pathogenesis of Chlamydia and allow reengineering of targetgenes to improve their expression for gene therapy.Here,we analyzed the codon usage of C.muridarum,C.trachomatis(here indicating biovar trachoma and LGV),C.pneumoniae,and C.psittaci using the codonusage database and the CUSP(Create a codon usage table)program of EMBOSS(The European MolecularBiology Open Software Suite).The results show that the four genomes have similar codon usage patterns,with a strong bias towards the codons with A and T at the third codon position.Compared with Homosapiens,the four chlamydial species show discordant seven or eight preferred codons.The ENC(effectivenumber of codons used in a gene)-plot reveals that the genetic heterogeneity in Chlamydia is constrained bythe G+C content,while translational selection and gene length exert relatively weaker influences.Moreover,mutational pressure appears to be the major determinant of the codon usage variation among the chlamydialgenes.In addition,we compared the codon preferences of C.trachomatis with those of E.coli,yeast,adenovirus and Homo sapiens.There are 23 codons showing distinct usage differences between C.trachomatisand E.coli,24 between C.trachomatis and adenovirus,21 between C.trachomatis and Homo sapiens,butonly six codons between C.trachomatis and yeast.Therefore,the yeast system may be more suitable for theexpression of chlamydial genes.Finally,we compared the codon preferences of C.trachomatis with those ofsix eukaryotes,eight prokaryotes and 23 viruses.There is a strong positive correlation between the differ-ences in coding GC content and the variations in codon bias(r=0.905,P<0,001).We conclude that thevariation of codon bias between C.trachomatis and other organisms is much less influenced by phylogeneticlineage and primarily determined by the extent of disparities in GC content.     

IVT-seq reveals extreme bias in RNA sequencing

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.     

Variation in evolutionary processes at different codon positions

Evolutionary studies commonly model single nucleotide substitutions and assume that they occur as independent draws from a unique probability distribution across the sequence studied. This assumption is violated for protein-coding sequences, and we consider modeling approaches where codon positions (CPs) are treated as separate categories of sites because within each category the assumption is more reasonable. Such "codon-position" models have been shown to explain the evolution of codon data better than homogenous models in previous studies. This paper examines the ways in which codon-position models outperform homogeneous models and characterizes the differences in estimates of model parameters across CPs. Using the PANDIT database of multiple species DNA sequence alignments, we quantify the differences in the evolutionary processes at the 3 CPs in a systematic and comprehensive manner, characterizing previously undescribed features of protein evolution. We relate our findings to the functional constraints imposed by the genetic code, protein function, and the types of mutation that cause synonymous and nonsynonymous codon changes. The results increase our understanding of selective constraints and could be incorporated into phylogenetic analyses or gene-finding techniques in the future. The methods used are extended to an overlapping reading frame data set, and we discover that overlapping reading frames do not necessarily cause more stringent evolutionary constraints.     

Analyses of evolutionary dynamics in viruses are hindered by a time-dependent bias in rate estimates

Time-scales of viral evolution and emergence have been studied widely, but are often poorly understood. Molecular analyses of viral evolutionary time-scales generally rely on estimates of rates of nucleotide substitution, which vary by several orders of magnitude depending on the timeframe of measurement. We analysed data from all major groups of viruses and found a strong negative relationship between estimates of nucleotide substitution rate and evolutionary timescale. Strikingly, this relationship was upheld both within and among diverse groups of viruses. A detailed case study of primate lentiviruses revealed that the combined effects of sequence saturation and purifying selection can explain this time-dependent pattern of rate variation. Therefore, our analyses show that studies of evolutionary time-scales in viruses require a reconsideration of substitution rates as a dynamic, rather than as a static, feature of molecular evolution. Improved modelling of viral evolutionary rates has the potential to change our understanding of virus origins.     

A thermodynamic theory of codon bias in viral genes

The relationship between degeneracy in the genetic code and the occurrence of a strong codon bias is examined, with particular reference to a group of viral genomes. The present paper shows how codon bias may have been imposed by thermodynamic considerations at the time the primitive DNA first formed in the primordial soup. Using a four-state Ising-like model with stacking interactions between successive base pairs, we show how primeval periodic DNA polymers could have arisen the remnants of which are still observed in codon biases today.     

ADHD symptoms in healthy adults are associated with stressful life events and negative memory bias

Stressful life events, especially Childhood Trauma, predict ADHD symptoms. Childhood Trauma and negatively biased memory are risk factors for affective disorders. The association of life events and bias with ADHD symptoms may inform about the etiology of ADHD. Memory bias was tested using a computer task in N = 675 healthy adults. Life events and ADHD symptoms were assessed using questionnaires. The mediation of the association between life events and ADHD symptoms by memory bias was examined. We explored the roles of different types of life events and of ADHD symptom clusters. Life events and memory bias were associated with overall ADHD symptoms as well as inattention and hyperactivity/impulsivity symptom clusters. Memory bias mediated the association of Lifetime Life Events, specifically Childhood Trauma, with ADHD symptoms. Negatively biased memory may be a cognitive marker of the effects of Childhood Trauma on the development and/or persistence of ADHD symptoms.     

Homologous recombination proteins are associated with centrosomes and are required for mitotic stability

In response to DNA damage, cells need robust repair mechanisms to complete the cell cycle successfully. Severe forms of DNA damage are repaired by homologous recombination (HR), in which the XRCC2 protein plays a vital role. Cells deficient in XRCC2 also show disruption of the centrosome, a key component of the mitotic apparatus. We find that this centrosome disruption is dynamic and when it occurs during mitosis it is linked directly to the onset of mitotic catastrophe in a significant fraction of the XRCC2-deficient cells. However, we also show for the first time that XRCC2 and other HR proteins, including the key recombinase RAD51, co-localize with the centrosome. Co-localization is maintained throughout the cell cycle, except when cells are finishing mitosis when RAD51 accumulates in the midbody between the separating cells. Taken together, these data suggest a tight functional linkage between the centrosome and HR proteins, potentially to coordinate the deployment of a DNA damage response at vulnerable phases of the cell cycle.     

Eukaryotic mRNAs encoding abundant and scarce proteins are statistically dissimilar in many structural features

It is well known that non-coding mRNA sequences are dissimilar in many structural features. For individual mRNAs correlations were found for some of these features and their translational efficiency. However, no systematic statistical analysis was undertaken to relate protein abundance and structural characteristics of mRNA encoding the given protein. We have demonstrated that structural and contextual features of eukaryotic mRNAs encoding high- and low-abundant proteins differ in the 5′ untranslated regions (UTR). Statistically, 5′ UTRs of low-expression mRNAs are longer, their guanine plus cytosine content is higher, they have a less optimal context of the translation initiation codons of the main open reading frames and contain more frequently upstream AUG than 5′ UTRs of high-expression mRNAs. Apart from the differences in 5′ UTRs, high-expression mRNAs contain stronger termination signals. Structural features of low- and high-expression mRNAs are likely to contribute to the yield of their protein products.     

Heat-shock proteins are associated with hnRNA in Drosophila melanogaster tissue culture cells

Ribonucleoprotein complexes of Drosophila melanogaster Kc tissue culture cells grown at 24°C or heat-shocked at 37°C were cross-linked in vivo by u.v. irradiation. Cross-linked heterogeneous nuclear ribonucleoprotein (hnRNP) complexes were fractionated by oligo(dT)-cellulose chromatography and CsCI density centrifugation. The hnRNP complexes of both 24°C and 37°C culture cells possess buoyant densities in CsCI between = 1.38 g/cm^-3 and 1.43 g/cm^-3. The ³⁵S-labelled proteins bound to the hnRNA of 37°C culture cells correspond in mol. wt. to the so-called heat-shock proteins of 70 K, 68 K, 27 K, 26 K, 23 K and 22 K. The 70 K and 68 K proteins are also present in hnRNP complexes of 24°C culture cells. In addition, several other Drosophila hnRNPs of 140 K, 56 K, 45 K, 43 K, 38 K, 37 K and 34 K, whose synthesis is strongly repressed under heat-shock conditions, could be identified. The results demonstrate that the so-called heat-shock proteins possess a function as RNPs.     

Eukaryotic precursor proteins are processed by Escherichia coli outer membrane protein OmpP.

A new specific endopeptidase that cleaves eukaryotic precursor proteins has been found in Escherichia coli K but not in E. coli B strains. After purification, protein sequencing and Western blotting, the endopeptidase was shown to be identical with E. coli outer membrane protein OmpP [Kaufmann, A., Stierhof, Y.-D. & Henning, U. (1994) J. Bacteriol. 176, 359-367]. Further characterization of enzymatic properties of the new peptidase was performed. Comparison of the cleavage specificities of the newly found endopeptidase and that of rat mitochondrial processing peptidase (MPP) showed that patterns of proteolytic cleavage on the investigated precursor proteins by both enzymes are similar. By using three mitochondrial precursor proteins, the specificity assigned to OmpP previously, a cleavage position between two basic amino-acid residues, was extended to a three amino-acid recognition sequence. Positions +1 to +3 of this extended recognition site consist of an amino-acid residue with a small aliphatic side chain such as alanine or serine, a large hydrophobic residue such as leucine or valine followed by an arginine residue. Additionally, structural motifs of the substrate seem to be required for OmpP cleavage.     

Active zone proteins are dynamically associated with synaptic ribbons in rat pinealocytes

Synaptic ribbons (SRs) are prominent organelles that are abundant in the ribbon synapses of sensory neurons where they represent a specialization of the cytomatrix at the active zone (CAZ). SRs occur not only in neurons, but also in neuroendocrine pinealocytes where their function is still obscure. In this study, we report that pinealocyte SRs are associated with CAZ proteins such as Bassoon, Piccolo, CtBP1, Munc13-1, and the motorprotein KIF3A and, therefore, consist of a protein complex that resembles the ribbon complex of retinal and other sensory ribbon synapses. The pinealocyte ribbon complex is biochemically dynamic. Its protein composition changes in favor of Bassoon, Piccolo, and Munc13-1 at night and in favor of KIF3A during the day, whereas CtBP1 is equally present during the night and day. The diurnal dynamics of the ribbon complex persist under constant darkness and decrease after stimulus deprivation of the pineal gland by constant light. Our findings indicate that neuroendocrine pinealocytes possess a protein complex that resembles the CAZ of ribbon synapses in sensory organs and whose dynamics are under circadian regulation.     

The myc proteins are not associated with chromatin in mitotic cells.

The protein products of cellular and viral myc oncogenes are detected in nuclei by immunofluorescence. No myc fluorescence is found in nucleoli. In mitotic cells the myc antigens are not found associated with metaphase chromosomes, but are diffusely distributed throughout the cytoplasm. Cytoplasmic myc fluorescence is first observed when chromatin begins to condense in early prophase. Granular nuclear myc fluorescence is again discerned in telophase cells, when the nuclear envelope is formed and becomes more prominent upon cytokinesis; concomitantly the diffuse cytoplasmic myc staining is lost. These results suggest that myc proteins not only bind to DNA or chromatin, but are also associated with other structural systems in the nuclei.     

Ancestral inference and the study of codon bias evolution: implications for molecular evolutionary analyses of the Drosophila melanogaster subgroup

Reliable inference of ancestral sequences can be critical to identifying both patterns and causes of molecular evolution. Robustness of ancestral inference is often assumed among closely related species, but tests of this assumption have been limited. Here, we examine the performance of inference methods for data simulated under scenarios of codon bias evolution within the Drosophila melanogaster subgroup. Genome sequence data for multiple, closely related species within this subgroup make it an important system for studying molecular evolutionary genetics. The effects of asymmetric and lineage-specific substitution rates (i.e., varying levels of codon usage bias and departures from equilibrium) on the reliability of ancestral codon usage was investigated. Maximum parsimony inference, which has been widely employed in analyses of Drosophila codon bias evolution, was compared to an approach that attempts to account for uncertainty in ancestral inference by weighting ancestral reconstructions by their posterior probabilities. The latter approach employs maximum likelihood estimation of rate and base composition parameters. For equilibrium and most non-equilibrium scenarios that were investigated, the probabilistic method appears to generate reliable ancestral codon bias inferences for molecular evolutionary studies within the D. melanogaster subgroup. These reconstructions are more reliable than parsimony inference, especially when codon usage is strongly skewed. However, inference biases are considerable for both methods under particular departures from stationarity (i.e., when adaptive evolution is prevalent). Reliability of inference can be sensitive to branch lengths, asymmetry in substitution rates, and the locations and nature of lineage-specific processes within a gene tree. Inference reliability, even among closely related species, can be strongly affected by (potentially unknown) patterns of molecular evolution in lineages ancestral to those of interest.     

Unusual codon bias occurring within insertion sequences in Escherichia coli

The large open reading frames of insertion sequences from Escherichia coli were examined for their spatial pattern of codon usage bias and distribution of rarely used codons. There is a bias in codon usage that is generally lower toward the terminal ends of the coding regions, which is reflected in the occurrence of an excess of nonpreferred codons in the 3 portions of the coding regions as compared with the 5 portions. In contrast, typical chromosomal genes have a lower codon usage bias toward the 5 ends of the coding regions. These results imply that the selective forces reflected in codon usage bias may differ according to position within the coding sequence. In addition, these constraints apparently differ in important ways between genes contained in insertion sequences and those in the chromosome.