The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.

Replicational and transcriptional selection on codon usage in Borrelia burgdorferi

With more than 10 fully sequenced, publicly available prokaryotic genomes, it is now becoming possible to gain useful insights into genome evolution. Before the genome era, many evolutionary processes were evaluated from limited data sets and evolutionary models were constructed on the basis of small amounts of evidence. In this paper, I show that genes on the Borrelia burgdorferi genome have two separate, distinct, and significantly different codon usages, depending on whether the gene is transcribed on the leading or lagging strand of replication. Asymmetrical replication is the major source of codon usage variation. Replicational selection is responsible for the higher number of genes on the leading strands, and transcriptional selection appears to be responsible for the enrichment of highly expressed genes on these strands. Replicational–transcriptional selection, therefore, has an influence on the codon usage of a gene. This is a new paradigm of codon selection in prokaryotes.

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.

Evolution of codon usage bias in Drosophila

We first review what is known about patterns of codon usage bias in Drosophila and make the following points: (i) Drosophila genes are as biased or more biased than those in microorganisms. (ii) The level of bias of genes and even the particular pattern of codon bias can remain phylogenetically invariant for very long periods of evolution. (iii) However, some genes, even very tightly linked genes, can change very greatly in codon bias across species. (iv) Generally G and especially C are favored at synonymous sites in biased genes. (v) With the exception of aspartic acid, all amino acids contribute significantly and about equally to the codon usage bias of a gene. (vi) While most individual amino acids that can use G or C at synonymous sites display a preference for C, there are exceptions: valine and leucine, which prefer G. (vii) Finally, smaller genes tend to be more biased than longer genes. We then examine possible causes of these patterns and discount mutation bias on three bases: there is little evidence of regional mutation bias in Drosophila, mutation bias is likely toward A+T (the opposite of codon usage bias), and not all amino acids display the preference for the same nucleotide in the wobble position. Two lines of evidence support a selection hypothesis based on tRNA pools: highly biased genes tend to be highly and/or rapidly expressed, and the preferred codons in highly biased genes optimally bind the most abundant isoaccepting tRNAs. Finally, we examine the effect of bias on DNA evolution and confirm that genes with high codon usage bias have lower rates of synonymous substitution between species than do genes with low codon usage bias. Surprisingly, we find that genes with higher codon usage bias display higher levels of intraspecific synonymous polymorphism. This may be due to opposing effects of recombination.

RNA editing in Arabidopsis mitochondria effects 441 C to U changes in ORFs

On the basis of the sequence of the mitochondrial genome in the flowering plant Arabidopsis thaliana, RNA editing events were systematically investigated in the respective RNA population. A total of 456 C to U, but no U to C, conversions were identified exclusively in mRNAs, 441 in ORFs, 8 in introns, and 7 in leader and trailer sequences. No RNA editing was seen in any of the rRNAs or in several tRNAs investigated for potential mismatch corrections. RNA editing affects individual coding regions with frequencies varying between 0 and 18.9% of the codons. The predominance of RNA editing events in the first two codon positions is not related to translational decoding, because it is not correlated with codon usage. As a general effect, RNA editing increases the hydrophobicity of the coded mitochondrial proteins. Concerning the selection of RNA editing sites, little significant nucleotide preference is observed in their vicinity in comparison to unedited C residues. This sequence bias is, per se, not sufficient to specify individual C nucleotides in the total RNA population in Arabidopsis mitochondria.

Highly expressed and alien genes of the Synechocystis genome

Comparisons of codon frequencies of genes to several gene classes are used to characterize highly expressed and alien genes on the Synechocystis PCC6803 genome. The primary gene classes include the ensemble of all genes (average gene), ribosomal protein (RP) genes, translation processing factors (TF) and genes encoding chaperone/degradation proteins (CH). A gene is predicted highly expressed (PHX) if its codon usage is close to that of the RP/TF/CH standards but strongly deviant from the average gene. Putative alien (PA) genes are those for which codon usage is significantly different from all four classes of gene standards. In Synechocystis, 380 genes were identified as PHX. The genes with the highest predicted expression levels include many that encode proteins vital for photosynthesis. Nearly all of the genes of the RP/TF/CH gene classes are PHX. The principal glycolysis enzymes, which may also function in CO2 fixation, are PHX, while none of the genes encoding TCA cycle enzymes are PHX. The PA genes are mostly of unknown function or encode transposases. Several PA genes encode polypeptides that function in lipopolysaccharide biosynthesis. Both PHX and PA genes often form significant clusters (operons). The proteins encoded by PHX and PA genes are described with respect to functional classifications, their organization in the genome and their stoichiometry in multi-subunit complexes.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.

Mercuric ion reduction and resistance in transgenic Arabidopsis thaliana plants expressing a modified bacterial merA gene.

With global heavy metal contamination increasing, plants that can process heavy metals might provide efficient and ecologically sound approaches to sequestration and removal. Mercuric ion reductase, MerA, converts toxic Hg2+ to the less toxic, relatively inert metallic mercury (Hg0) The bacterial merA sequence is rich in CpG dinucleotides and has a highly skewed codon usage, both of which are particularly unfavorable to efficient expression in plants. We constructed a mutagenized merA sequence, merApe9, modifying the flanking region and 9% of the coding region and placing this sequence under control of plant regulatory elements. Transgenic Arabidopsis thaliana seeds expressing merApe9 germinated, and these seedlings grew, flowered, and set seed on medium containing HgCl2 concentrations of 25-100 microM (5-20 ppm), levels toxic to several controls. Transgenic merApe9 seedlings evolved considerable amounts of Hg0 relative to control plants. The rate of mercury evolution and the level of resistance were proportional to the steady-state mRNA level, confirming that resistance was due to expression of the MerApe9 enzyme. Plants and bacteria expressing merApe9 were also resistant to toxic levels of Au3+. These and other data suggest that there are potentially viable molecular genetic approaches to the phytoremediation of metal ion pollution.

Transgenic barley expressing a protein-engineered, thermostable (1,3-1,4)-beta-glucanase during germination.

The codon usage of a hybrid bacterial gene encoding a thermostable (1,3-1,4)-beta-glucanase was modified to match that of the barley (1,3-1,4)-beta-glucanase isoenzyme EII gene. Both the modified and unmodified bacterial genes were fused to a DNA segment encoding the barley high-pI alpha-amylase signal peptide downstream of the barley (1,3-1,4)-beta-glucanase isoenzyme EII gene promoter. When introduced into barley aleurone protoplasts, the bacterial gene with adapted codon usage directed synthesis of heat stable (1,3-1,4)-beta-glucanase, whereas activity of the heterologous enzyme was not detectable when protoplasts were transfected with the unmodified gene. In a different expression plasmid, the codon modified bacterial gene was cloned downstream of the barley high-pI alpha-amylase gene promoter and signal peptide coding region. This expression cassette was introduced into immature barley embryos together with plasmids carrying the bar and the uidA genes. Green, fertile plants were regenerated and approximately 75% of grains harvested from primary transformants synthesized thermostable (1,3-1,4)-beta-glucanase during germination. All three trans genes were detected in 17 progenies from a homozygous T1 plant.

Alternative translation initiation site usage results in two functionally distinct forms of the GATA-1 transcription factor.

GATA-1 is a zinc-finger transcription factor that plays a critical role in the normal development of hematopoietic cell lineages. In human and murine erythroid cells a previously undescribed 40-kDa protein is detected with GATA-1-specific antibodies. We show that the 40-kDa GATA-1 (GATA-1s) is produced by the use of an internal AUG initiation codon in the GATA-1 transcript. The GATA-1 proteins share identical binding activity and form heterodimers in erythroleukemic cells but differ in their transactivation potential and in their expression in developing mouse embryos.

Generation of secretable and nonsecretable interleukin 15 isoforms through alternate usage of signal peptides

Two isoforms of human interleukin 15 (IL-15) exist. One isoform has a shorter putative signal peptide (21 amino acids) and its transcript shows a tissue distribution pattern that is distinct from that of the alternative IL-15 isoform with a 48-aa signal peptide. The 21-aa signal isoform is preferentially expressed in tissues such as testis and thymus. Experiments using different combinations of signal peptides and mature proteins (IL-2, IL-15, and green fluorescent protein) showed that the short signal peptide regulates the fate of the mature protein by controlling the intracellular trafficking to nonendoplasmic reticulum sites, whereas the long signal peptide both regulates the rate of protein translation and functions as a secretory signal peptide. As a consequence, the IL-15 associated with the short signal peptide is not secreted, but rather is stored intracellularly, appearing in the nucleus and cytoplasmic components. Such production of an intracellular lymphokine is not typical of other soluble interleukin systems, suggesting a biological function for IL-15 as an intracellular molecule.

Infrequent Translation of a Nonsense Codon Is Sufficient to Decrease mRNA Level

In many organisms nonsense mutations decrease the level of mRNA. In the case of mammalian cells, it is still controversial whether translation is required for this nonsense-mediated RNA decrease (NMD). Although previous analyzes have shown that conditions that impede translation termination at nonsense codons also prevent NMD, the residual level of termination was unknown in these experiments. Moreover, the conditions used to impede termination might also have interfered with NMD in other ways. Because of these uncertainties, we have tested the effects of limiting translation of a nonsense codon in a different way, using two mutations in the immunoglobulin μ heavy chain gene. For this purpose we exploited an exceptional nonsense mutation at codon 3, which efficiently terminates translation but nonetheless maintains a high level of μ mRNA. We have shown 1) that translation of Ter462 in the double mutant occurs at only ∼4% the normal frequency, and 2) that Ter462 in cis with Ter3 can induce NMD. That is, translation of Ter462 at this low (4%) frequency is sufficient to induce NMD.

Autoinhibition of serotonin cells: An intrinsic regulatory mechanism sensitive to the pattern of usage of the cells

After periods of high-frequency firing, the normal rhythmically active serotonin (5HT)-containing neurosecretory neurons of the lobster ventral nerve cord display a period of suppressed spike generation and reduced synaptic input that we refer to as “autoinhibition.” The duration of this autoinhibition is directly related to the magnitude and duration of the current injection triggering the high-frequency firing. More interesting, however, is that the autoinhibition is inversely related to the initial firing frequency of these cells within their normal range of firing (0.5–3 Hz). This allows more active 5HT neurons to resume firing after shorter durations of inhibition than cells that initially fired at slower rates. Although superfused 5HT inhibits the spontaneous firing of these cells, the persistence of autoinhibition in saline with no added calcium, in cadmium-containing saline, and in lobsters depleted of serotonin suggests that intrinsic membrane properties account for the autoinhibition. A similar autoinhibition is seen in spontaneously active octopamine neurons but is absent from spontaneously active γ-aminobutyric acid cells. Thus, this might be a characteristic feature of amine-containing neurosecretory neurons. The 5HT cells of vertebrate brain nuclei share similarities in firing frequencies, spike shapes, and inhibition by 5HT with the lobster cells that were the focus of this study. However, the mechanism suggested to underlie autoinhibition in vertebrate neurons is that 5HT released from activated or neighboring cells acts back on inhibitory autoreceptors that are found on the dendrites and cell bodies of these neurons.