Trends in Codon and Amino Acid Usage in Thermotoga maritima

The usage of synonymous codons and the frequencies of amino acids were investigated in the complete genome of the bacterium Thermotoga maritima using a multivariate statistical approach. The GC3 content of each gene was the most prominent source of variation of codon usage. Surprisingly the usage of UGU and UGC (synonymous triplets coding for Cys, the least frequent amino acid in this species) was detected as the second most prominent source of variation. However, this result is probably an artifact due to the very low frequency of Cys together with the nonbiased composition of this genome. The third trend was related to the preferential usage of a subset of codons among highly expressed genes, and these triplets are presumed to be translationally optimal. Concerning the amino acid usage, the hydropathy level of each protein (and therefore the frequency of charged residues) was the main trend, while the second factor was related to the frequency of usage of the smaller residues, suggesting that the cell economy strongly influences the architecture of the proteins. The third axis of the analysis discriminated the usage of Phe, Tyr, Trp (aromatic residues) plus Cys, Met, and His. These six residues have in common the property of being the preferential targets of reactive oxygen species, and therefore the anaerobic condition of T. maritima is an important factor for the amino acid frequencies. Finally, the Cys content of each protein was the fourth trend. Received: 22 June 2001 / Accepted: 1 October 2001     

A Test of Amino Acid Reversibility

In studies of molecular evolution, the assumption that protein evolution is reversible has often been made, but rarely tested. Here we use a large set of orthologous murid protein coding sequences to perform a simple test of reversibility, and find no evidence to reject the assumption of reversibility in protein evolution. Received: 10 October 2000 / Accepted: 18 January 2001     

Volume-regulatory Amino Acid Release from the Protozoan Parasite Crithidia luciliae

The unicellular protozoan parasite, Crithidia luciliae, responded to osmotic swelling by undergoing a regulatory volume decrease. This process was accompanied by the efflux of amino acids (predominantly alanine, proline and glycine). The relative loss of the electroneutral amino acids proline, valine, alanine and glycine was greater than that for the anionic amino acid, glutamate; there was negligible loss of the cationic amino acids, lysine, arginine and ornithine. The characteristics of amino acid release were investigated using a radiolabeled form of the nonmetabolized alanine analogue α-aminoisobutyrate. α-Aminoisobutyrate efflux was activated within a few seconds of a reduction of the osmolality, and inactivated rapidly (again within a few seconds) on restoration of isotonicity. The initial rate of efflux of α-aminoisobutyrate from cells in hypotonic medium was unaffected by the extracellular amino acid concentration. Hypotonically activated α-aminoisobutyrate efflux (as well as the associated regulatory volume decrease) was inhibited by the sulfhydryl reagent N-ethylmaleimide but was not inhibited by a range of anion transport blockers. As in the efflux experiments, unidirectional influx rates for α-aminoisobutyrate increased markedly following reduction of the osmolality, consistent with the swelling-activated amino acid release mechanism allowing the flux of solutes in both directions. Hypotonically activated α-aminoisobutyrate influx showed no tendency to saturate up to an extracellular concentration of 50 mm. The functional characteristics of the amino acid release mechanism are those of a channel, with a preference for electroneutral and anionic amino acids over cationic amino acids. However, the pharmacology of the system differs from that of the anion-selective channels that are thought to mediate the volume-regulatory efflux of organic osmolytes from vertebrate cells. Received: 13 May 1996/Revised: 9 July 1996     

Nucleotide Composition Bias Affects Amino Acid Content in Proteins Coded by Animal Mitochondria

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content. Received: 3 July 1996 / Accepted: 6 November 1996     

Evolution of Chordate Actin Genes: Evidence from Genomic Organization and Amino Acid Sequences

The origin and evolutionary relationship of actin isoforms was investigated in chordates by isolating and characterizing two new ascidian cytoplasmic and muscle actin genes. The exon–intron organization and sequences of these genes were compared with those of other invertebrate and vertebrate actin genes. The gene HrCA1 encodes a cytoplasmic (nonmuscle)-type actin, whereas the MocuMA2 gene encodes an adult muscle-type actin. Our analysis of these genes showed that intron positions are conserved among the deuterostome actin genes. This suggests that actin gene families evolved from a single actin gene in the ancestral deuterostome. Sequence comparisons and molecular phylogenetic analyses also suggested a close relationship between the ascidian and vertebrate actin isoforms. It was also found that there are two distinct lineages of muscle actin isoforms in ascidians: the larval muscle and adult body-wall isoforms. The four muscle isoforms in vertebrates show a closer relationship to each other than to the ascidian muscle isoforms. Similarly, the two cytoplasmic isoforms in vertebrates show a closer relationship to each other than to the ascidian and echinoderm cytoplasmic isoforms. In contrast, the two types of ascidian muscle actin diverge from each other. The close relationship between the ascidian larval muscle actin and the vertebrate muscle isoforms was supported by both neighbor-joining and maximum parsimony analyses. These results suggest that the chordate ancestor had at least two muscle actin isoforms and that the vertebrate actin isoforms evolved after the separation of the vertebrates and urochordates. Received: 20 June 1996 / Accepted: 16 October 1996     

Synonymous and Nonsynonymous Substitutions in Mammalian Genes: Intragenic Correlations

Previous investigations indicated that synonymous and nonsynonymous substitution rates are correlated in mammalian genes. In the present work, this correlation has been studied at the intragenic level using a dataset of 48 orthologous genes from species belonging to at least four different mammalian orders. The results obtained show that the intragenic variability in synonymous rates is correlated with that of nonsynonymous rates. Moreover, the variation in GC level (and especially of C level) of silent positions along each gene is correlated with the variation in synonymous rate. These results reinforce the previous conclusions that synonymous and nonsynonymous rates as well as GC levels of silent positions are to some extent under common selective constraints. Received: 10 July 1997 / Accepted: 13 August 1997     

Plastid Genome Phylogeny and a Model of Amino Acid Substitution for Proteins Encoded by Chloroplast DNA

Maximum likelihood (ML) phylogenies based on 9,957 amino acid (AA) sites of 45 proteins encoded in the plastid genomes of Cyanophora, a diatom, a rhodophyte (red algae), a euglenophyte, and five land plants are compared with respect to several properties of the data, including between-site rate variation and aberrant amino acid composition in individual species. Neighbor-joining trees from AA LogDet distances and ML analyses are seen to be congruent when site rate variability was taken into account. Four feasible trees are identified in these analyses, one of which is preferred, and one of which is almost excluded by statistical criteria. A transition probability matrix for the general reversible Markov model of amino acid substitutions is estimated from the data, assuming each of these four trees. In all cases, the tree with diatom and rhodophyte as sister taxa was clearly favored. The new transition matrix based on the best tree, called cpREV, takes into account distinct substitution patterns in plastid-encoded proteins and should be useful in future ML inferences using such data. A second rate matrix, called cpREV*, based on a weighted sum of rate matrices from different trees, is also considered. Received: 3 June 1999 / Accepted: 26 November 1999     

Gap Junctional Channels Regulate Acid Secretion in the Mammalian Gastric Gland

Gap junction channels are regarded as a primary pathway for intercellular message transfer, including calcium wave propagation. Our study identified two gap junctional proteins, connexin26 and connexin32, in rat gastric glands by RT-PCR, Western blot analysis, and immunofluorescence. We demonstrated a potential physiological role of the gap junctional channels in the acid secretory process using the calcium indicator fluo-3, and microinjection of Lucifer Yellow. Application of gastrin (10⁻⁷ m) to the basolateral membrane resulted in the induction of uniphasic calcium signals in adjacent parietal cells. In addition, single parietal cell microinjections in intact glands with the cell-impermeant dye Lucifer Yellow resulted in a transfer of dye from the injected cell to the adjacent parietal cell following gastrin stimulation, demonstrating gastrin-induced cell-to-cell communication. Both calcium wave propagation and Lucifer Yellow transfer were blocked by the gap junction inhibitor 18α-glycyrrhetinic acid. Our studies demonstrate that functional gap junction channels in gastric glands provide an effective means for rapid cell-to-cell communication and allow for the rapid onset of acid secretion. Received: 4 December 2000/Revised: 5 June 2001     

New Glycoprotein-Associated Amino Acid Transporters

The L-type amino acid transporter LAT1 has recently been identified as being a disulfide-linked ``light chain' of the ubiquitously expressed glycoprotein 4F2hc/CD98. Several LAT1-related transporters have been identified, which share the same putative 12-transmembrane segment topology and also associate with the single transmembrane domain 4F2hc protein. They display differing amino acid substrate specificities, transport kinetics and localizations such as, for instance, y⁺LAT1 which is localized at the basolateral membrane of transporting epithelia, and the defect of which causes lysinuric protein intolerance. The b^0,+AT transporter which associates with the 4F2hc-related rBAT protein to form the luminal high-affinity diamino acid transporter defective in cystinuria, belongs to the same family of glycoprotein-associated amino acid transporters (gpaATs). These glycoprotein-associated transporters function as amino acid exchangers. They extend the specificity range of vectorial amino acid transport when located in the same membrane as carriers that unidirectionally transport one of the exchanged substrates. gpaATs belong to a phylogenetic cluster within the amino acid/polyamine/choline (APC) superfamily of transporters. This cluster, which we designate the LAT family (named after its first vertebrate member), includes some members from nematodes, yeast and bacteria. The latter of these proteins presumably lack association with a second subunit. In this review, we focus on the animal members of the LAT cluster that form, together with some of the nematode members, the family of glycoprotein-associated amino acid transporters (gpaAT family). Received: 20 July 1999/Revised: 7 September 1999     

The Binding Specificity of Amino Acid Transport System y+L in Human Erythrocytes is Altered by Monovalent Cations

System y⁺L is a broad-scope amino acid transporter which binds and translocates cationic and neutral amino acids. Na⁺ replacement with K⁺ does not affect lysine transport, but markedly decreases the affinity of the transporter for l-leucine and l-glutamine. This observation suggests that the specificity of system y⁺L varies depending on the ionic composition of the medium. Here we have studied the interaction of the carrier with various amino acids in the presence of Na⁺, K⁺, Li⁺ and guanidinium ion. In agreement with the prediction, the specificity of system y⁺L was altered by the monovalent cations. In the presence of Na⁺, l-leucine was the neutral amino acid that interacted more powerfully. Elongation of the side chain (glycine - l-norleucine) strengthened binding. In contrast, bulkiness at the level of the β carbon was detrimental. In K⁺, the carrier behaved as a cationic amino acid specific carrier, interacting weakly with neutral amino acids. Li⁺ was found to potentiate neutral amino acid binding and in general the apparent affinities were higher than in Na⁺; elongation of the nonpolar side chain made a more important contribution to binding and the carrier was more tolerant towards β carbon substitution. Guanidinium stimulated the interaction of the carrier with neutral amino acids, but the effect was restricted to certain analogues (e.g., l-leucine, l-glutamine, l-methionine). Thus, in the presence of guanidinium, the carrier discriminates sharply among different neutral amino acids. The results suggest that the monovalent cations stabilize different carrier conformations. Received: 22 January 1996/Revised: 26 April 1996     

On the Evolution of Redundancy in Genetic Codes

We simulate a deterministic population genetic model for the coevolution of genetic codes and protein-coding genes. We use very simple assumptions about translation, mutation, and protein fitness to calculate mutation-selection equilibria of codon frequencies and fitness in a large asexual population with a given genetic code. We then compute the fitnesses of altered genetic codes that compete to invade the population by translating its genes with higher fitness. Codes and genes coevolve in a succession of stages, alternating between genetic equilibration and code invasion, from an initial wholly ambiguous coding state to a diversified frozen coding state. Our simulations almost always resulted in partially redundant frozen genetic codes. Also, the range of simulated physicochemical properties among encoded amino acids in frozen codes was always less than maximal. These results did not require the assumption of historical constraints on the number and type of amino acids available to codes nor on the complexity of proteins, stereochemical constraints on the translational apparatus, nor mechanistic constraints on genetic code change. Both the extent and timing of amino-acid diversification in genetic codes were strongly affected by the message mutation rate and strength of missense selection. Our results suggest that various omnipresent phenomena that distribute codons over sites with different selective requirements—such as the persistence of nonsynonymous mutations at equilibrium, the positive selection of the same codon in different types of sites, and translational ambiguity—predispose the evolution of redundancy and of reduced amino acid diversity in genetic codes. Received: 21 December 2000 / Accepted: 12 March 2001     

rBAT is an Amino Acid Exchanger with Variable Stoichiometry

The rBAT protein, when expressed in Xenopus oocytes, was previously shown to reproduce the selectivity of the Na⁺-independent neutral and basic amino acid transport system called b^o,+. More recently, the capacity of rBAT to generate a transmembrane current was demonstrated when addition of neutral amino acids stimulated the efflux of cations (presumably basic amino acids) in rBAT-injected oocytes. In the present paper, aminoisobutyric acid (AIB), a neutral amino acid analogue, was shown to induce outward currents (efflux of basic amino acids) through rBAT similar to those caused by alanine in terms of affinity, maximal currents and I-V curves. Despite generating similar currents, the AIB transport rate was more than 30 times lower than that of alanine, thus challenging the assumption that rBAT functions as a classical exchanger. Experiments using a cut-open oocyte voltage clamp demonstrated that AIB was capable of stimulating rBAT-mediated currents from either side of the membrane. AIB, like alanine, was able to stimulate the efflux of radiolabeled alanine and arginine while no rBAT-mediated efflux was measurable in the absence of external rBAT substrates. These results demonstrate that (i) the presence of amino acids is required on both sides of the membrane for rBAT to mediate amino acid flux and thus rBAT must be some type of exchanger but (ii) rBAT-mediated amino acid influx is not stoichiometrically related to the efflux. A model of a ``double gated pore' is proposed to account for these properties of rBAT, which contravene standard models of exchangers and other transporters. Received: 15 June 1995/Revised: 21 September 1995     

Rates of Conservative and Radical Nonsynonymous Nucleotide Substitutions in Mammalian Nuclear Genes

To understand the process and mechanism of protein evolution, it is important to know what types of amino acid substitutions are more likely to be under selection and what types are mostly neutral. An amino acid substitution can be classified as either conservative or radical, depending on whether it involves a change in a certain physicochemical property of the amino acid. Assuming Kimura's two-parameter model of nucleotide substitution, I present a method for computing the numbers of conservative and radical nonsynonymous (amino acid altering) nucleotide substitutions per site and estimate these rates for 47 nuclear genes from mammals. The results are as follows. (1) The average radical/conservative rate ratio is 0.81 for charge changes, 0.85 for polarity changes, and 0.49 when both polarity and volume changes are considered. (2) The radical/conservative rate ratio is positively correlated with the nonsynonymous/synonymous rate ratio for charge changes or when both polarity and volume changes are considered. (3) Both the conservative/synonymous rate ratio and the radical/synonymous rate ratio are lower in the rodent lineage than in the primate or artiodactyl lineage, suggesting more intense purifying selection in the rodent lineage, for both conservative and radical nonsynonymous substitutions. (4) Neglecting transition/transversion bias would cause an underestimation of both radical and conservative rates and the ratio thereof. (5) Transversions induce more dramatic genetic alternations than transitions in that transversions produce more amino acid altering changes and among which, more radical changes. Received: 6 April 1999 / Accepted: 16 August 1999     

Molecular Evolution Modeled as a Fractal Renewal Point Process in Agreement with the Dispersion of Substitutions in Mammalian Genes

A fractal renewal point process (FRPP) is used to model molecular evolution in agreement with the relationship between the variance and the mean numbers of nonsynonymous and synonymous substitutions in mammals. Like other episodic models such as the doubly stochastic Poisson process, this model accounts for the large variances observed in amino acid substitution rates, but unlike certain other episodic models, it also accounts for the increase in the index of dispersion with the mean number of substitutions in Ohta's (1995) data. We find that this correlation is significant for nonsynonymous substitutions at the 1% level and for synonymous substitutions at the 10% level, even after removing lineage effects and when using Bulmer's (1989) unbiased estimator of the index of dispersion. This model is simpler than most other overdispersed models of evolution in the sense that it is fully specified by a single interevent probability distribution. Interpretations in terms of chaotic dynamics and in terms of chance and selection are discussed. Received: 12 January 1998 / Accepted: 19 May 1998     

Phenylalanine-Binding RNAs and Genetic Code Evolution

We isolated RNAs by selection–amplification, selecting for affinity to Phe–Sepharose and elution with free l-phenylalanine. Constant sequences did not contain Phe condons or anticodons, to avoid any possible confounding influence on initially randomized sequences. We examined the eight most frequent Phe-binding RNAs for inclusion of coding triplets. Binding sites were defined by nucleotide conservation, protection, and interference data. Together these RNAs comprise 70% of the 105 sequenced RNAs. The K _D for the strongest sites is ≈50 μM free amino acid, with strong stereoselectivity. One site strongly distinguishes free Phe from Trp and Tyr, a specificity not observed previously. In these eight Phe-binding RNAs, Phe codons are not significantly associated with Phe binding sites. However, among 21 characterized RNAs binding Phe, Tyr, Arg, and Ile, containing 1342 total nucleotides, codons are 2.7-fold more frequent within binding sites than in surrounding sequences in the same molecules. If triplets were not specifically related to binding sites, the probability of this distribution would be 4.8 × 10⁻¹¹. Therefore, triplet concentration within amino acid binding sites taken together is highly likely. In binding sites for Arg, Tyr, and Ile cognate codons are overrepresented. Thus Arg, Tyr, and Ile may be amino acids whose codons were assigned during an era of direct RNA–amino acid affinity. In contrast, Phe codons arguably were assigned by another criterion, perhaps during later code evolution.     

Circular Permutations in the Molecular Evolution of DNA Methyltransferases

Circular permutations of genes during molecular evolution often are regarded as elusive, although a simple model can explain these rearrangements. The model assumes that first a gene duplication of the precursor gene occurs in such a way that both genes become fused in frame, leading to a tandem protein. After generation of a new start codon within the 5′ part of the tandem gene and a stop at an equivalent position in the 3′ part of the gene, a protein is encoded that represents a perfect circular permutation of the precursor gene product. The model is illustrated here by the molecular evolution of adenine-N⁶ DNA methyltransferases. β- and γ-type enzymes of this family can be interconverted by a single circular permutation event. Interestingly, tandem proteins, proposed as evolutionary intermediates during circular permutation, can be directly observed in the case of adenine methyltransferases, because some enzymes belonging to type IIS, like the FokI methyltransferase, are built up by two fused enzymes, both of which are active independently of each other. The mechanism for circular permutation illustrated here is very easy and applicable to every protein. Thus, circular permutation can be regarded as a normal process in molecular evolution and a changed order of conserved amino acid motifs should not be interpreted to argue against divergent evolution. Received: 17 November 1998 / Accepted: 19 February 1999     

Sequences and Evolution of Human and Squirrel Monkey Blue Opsin Genes

The sequences of the entire blue opsin gene in the squirrel monkey (Saimiri boliviensis) and the five introns of the human blue opsin gene were obtained. Intron 3 of these genes contains an Alu sequence and intron 4 contains a partial mer13 sequence. A comparison of the squirrel monkey opsin sequence with published mammalian opsin sequences shows that features believed to be functionally critical are all conserved. However, the blue opsin has evolved twice as fast as rhodopsin and is only as conservative as the β globin, which has evolved at the average rate of mammalian proteins. Interestingly, the interhelical loops are, on average, actually more conservative than the transmembrane α helical regions. The introns of the blue opsin gene have evolved at the average rate of introns in primate genes. Received: 5 August 1996 / Accepted: 2 October 1996