Positive Selection Scanning Reveals Decoupling of Enzymatic Activities of Carbamoyl Phosphate Synthetase in Helicobacter pylori

In an effort to detect factors which may be under positive selection, a survey for such genes in two pathogenic strains of Helicobacter pylori (J99 and 26695) was performed. Based on an analysis of synonymous and nonsynonymous substitutions, we identified 19 candidate genes under positive selection. A search for homologues with known crystallographic structures revealed Escherichia coli carbomoyl phosphate synthetase as a homologue of H. pylori carbamoyl phosphate synthetase. Carbamoyl phosphate synthetase as isolated from E. coli is a heterodimeric enzyme that possesses two different but coupled functionalities and is involved in the first committed step in the separate biosynthetic pathways for arginine and pyrimidine nucleotides. In this study, we provide evidence indicating that one of these functionalities appears to be under selective pressure. Reports from previously published site-directed mutagenesis studies point to a decoupling of amidotransferase and synthetase activities. Implications of these findings for a metabolic enzyme under positive selection are discussed in terms of the mechanisms of H. pylori pathogenesis. Received: 11 June 2001 / Accepted: 12 September 2001     

Evolutionary Lability of Context-Dependent Codon Bias in Bacteria

In bacteria, synonymous codon usage can be considerably affected by base composition at neighboring sites. Such context-dependent biases may be caused by either selection against specific nucleotide motifs or context-dependent mutation biases. Here we consider the evolutionary conservation of context-dependent codon bias across 11 completely sequenced bacterial genomes. In particular, we focus on two contextual biases previously identified in Escherichia coli; the avoidance of out-of-frame stop codons and AGG motifs. By identifying homologues of E. coli genes, we also investigate the effect of gene expression level in Haemophilus influenzae and Mycoplasma genitalium. We find that while context-dependent codon biases are widespread in bacteria, few are conserved across all species considered. Avoidance of out-of-frame stop codons does not apply to all stop codons or amino acids in E. coli, does not hold for different species, does not increase with gene expression level, and is not relaxed in Mycoplasma spp., in which the canonical stop codon, TGA, is recognized as tryptophan. Avoidance of AGG motifs shows some evolutionary conservation and increases with gene expression level in E. coli, suggestive of the action of selection, but the cause of the bias differs between species. These results demonstrate that strong context-dependent forces, both selective and mutational, operate on synonymous codon usage but that these differ considerably between genomes. Received: 6 May 1999 / Accepted: 29 October 1999     

Limitations of compositional approach to identifying horizontally transferred genes

Genes with atypical G+C content and pattern of codon usage in a certain genome are possibly of exotic origin, and this idea has been applied to identify horizontal events. In this way, it was postulated that a total of 755 genes in the E. coli genome are relics of horizontal events after the divergence of E. coli from the Salmonella lineage 100 million years ago (Lawrence and Ochman, 1998). In this paper we propose a new way to study sequence composition more thoroughly. We found that although the 755 genes differ in composition from other genes in the E. coli genome, the difference is minor. If we accepted that these genes are horizontally transferred, then (1) it would be more likely that they were transferred from genomes evolutionarily closely related to E. coli; but (2) the dating method used by Lawrence and Ochman (1997, 1998) largely underestimated the average age of introduced sequences in the E. coli genome, in particular, most of the 755 genes should be introduced into E. coli before, instead of after, the divergence of E. coli from the Salmonella lineage. Our study reveals that atypical G+C content and pattern of codon usage are not reliable indicators of horizontal gene transfer events. Received: 27 September 2000 / Accepted: 9 April 2001     

Mutants of the Lactose Carrier of Escherichia coli Which Show Altered Sugar Recognition Plus a Severe Defect in Sugar Accumulation

Lactose and melibiose are actively accumulated by the wild-type Escherichia coli lactose carrier, which is an integral membrane protein energized by the proton motive force. Mutants of the E. coli lactose carrier were isolated by their ability to grow on minimal plates with succinate plus IPTG in the presence of the toxic lactose analog β-thio-o-nitrophenylgalactoside (TONPG). TONPG-resistant mutants were streaked on melibiose MacConkey indicator plates, and red clones were picked. These melibiose positive mutants were then streaked on lactose MacConkey plates, and white clones were picked. Transport assays indicated that the mutants had altered sugar recognition and a defect in sugar accumulation. The mutants had a poor apparent K _m for both lactose and melibiose in transport. One mutant had almost no ability to take up lactose, but melibiose downhill transport was 58% (V _max ) of normal. All of the mutants accumulated methyl-α-d-galactopyranoside (TMG) to only 8% or less of normal, and two failed to accumulate. Immunoblot analysis of the mutant lactose carrier proteins indicated that loss of sugar transport activity was not due to loss of expression in the membrane. Nucleotide sequencing of the lacY gene from the mutants revealed changes in the following amino acids of the lactose carrier: M23I, W151L, G257D, A295D and G377V. Two of the mutants (G257D and G377V) are novel in that they represent the first amino acids in periplasmic loops to be implicated with changes in sugar recognition. We conclude that the amino acids M23, W151, G257, A295 and G377 of the E. coli lactose carrier play either a direct or an indirect role in sugar recognition and accumulation. Received: 12 October 1999/Revised: 21 December 1999     

Physical and Genetic Map of the Genome of Buchnera, the Primary Endosymbiont of the Pea Aphid Acyrthosiphon pisum

Abstract The genome of Buchnera, an endosymbiotic bacterium of the pea aphid Acyrthosiphon pisum, was characterized by pulse-field gel electrophoresis (PFGE) as a circular DNA molecule of 657 kb. The enzymes I-CeuI, CpoI, ApaI, SmaI, NaeI, SacII, MluI, FspI, and NruI were used to cleave the DNA of Buchnera into fragments of suitable size for PFGE analysis. A physical map of the Buchnera genome, including restriction fragments from seven of these enzymes, was constructed using double cutting, partial digestion, and hybridization with linking fragments, and 29 genes and operons were localized on the map. In addition, the genomic map of Buchnera was compared with those of Escherichia coli and Haemophilus influenzae. The gene order in Buchnera is more similar to that of E. coli than to H. influenzae. The dramatic shrinkage of the Buchnera genome compared with those of other members of the closely related Enterobacteriaceae family is discussed in terms of evolution under the influence of the intracellular symbiotic association. Received: 20 May 1998 / Accepted: 1 July 1998     

Analysis of the cellular functions of Escherichia coli operons and their conservation in Bacillus subtilis

The common assumption of operons as composed of genes that cooperate in a biological process is confirmed here by showing that Escherichia coli operons tend to be composed of genes that belong to the same general class of cellular function. Furthermore, the comparison between the genomic organization of E. coli and that of Bacillus subtilis shows that the genes that are homologous to genes that belong to experimentally characterized E. coli operons tend to cluster in neighboring regions of the genome. This tendency is greater for the subset of E. coli operons whose genes belong to a single functional class. These observations indicate strong evolutionary pressure that, translated into functional constraints, leads to the inclusion of many essential functions in conserved operons and clusters in these two distant species.     

Identification of Sequences Encoding the Detoxification Metalloisomerase Glyoxalase I in Microbial Genomes from Several Pathogenic Organisms

The ubiquitous glyoxalase system, which is composed of two enzymes, removes cellular cytotoxic methylglyoxal (MG). In an effort to identify critical residues conserved in the evolution of the first enzyme in this system, glyoxalase I (GlxI), as well as the structural implications of sequence alterations in this enzyme, a search of the National Center for Biotechnology Information (NCBI) database of unfinished genomes was undertaken. Eleven putative GlxI sequences from pathogenic organisms were identified and analyses of these sequences in relation to the known and previously identified GlxI enzymes were performed. Several of these sequences show a very high similarity to the Escherichia coli GlxI sequence, most notably the 79% identity of the sequence identified from Yersinia pestis, the causative agent of bubonic plague. In addition to the conservation of residues critical to binding the catalytic metal in all of the proposed GlxI enzymes, four regions in the Homo sapiens GlxI enzyme are absent in all of the bacterial GlxI sequences, with the exception of Pseudomonas putida. Removal of these regions may alter the active-site conformation of the bacterial enzymes in relation to that of the H. sapiens. These differences may be targeted for the development of inhibitors selective to the bacterial enzymes. Received: 13 October 1999 / Accepted: 17 January 2000     

Correlation Between Transepithelial Na+ Transport and Transepithelial Water Movement Across Isolated Frog Skin (Rana esculenta)

In the present work the coupling under short-circuited conditions between the net Na⁺-influx across isolated frog skin and the transepithelial transport of water was examined i.e., the short-circuit current (I _sc ) and the transepithelial water movement (TEWM) were measured simultaneously. It has been shown repeatedly that the I _sc across isolated frog skin is equal to the net transepithelial Na⁺ transport. Furthermore the coupling between transepithelial uptake of NaCl under open-circuit conditions and TEWM was also measured. The addition of antidiuretic hormone (AVT) to skins incubated under short-circuited conditions resulted in an increase in the I _sc and TEWM. Under control conditions I _sc was 9.14 ± 2.43 and in the presence of AVT 45.9 ± 7.3 neq cm⁻² min⁻¹ (n= 9) and TEWM changed from 12.45 ± 4.46 to 132.8 ± 15.8 nL cm⁻² min⁻¹. The addition of the Na⁺ channel blocking agent amiloride resulted in a reduction both in I _sc and TEWM, and a linear correlation between I _sc and TEWM was found. The correlation corresponds to that 160 ± 15 (n= 7) molecules of water follow each Na⁺ across the skin. In another series of experiments it was found that there was a linear correlation between I _sc and the increase in apical osmolarity needed to stop the TEWM. The data presented indicate that the observed coupling between the net transepithelial Na⁺ transport and TEWM is caused by local osmosis. Received: 16 October 1996/Revised: 6 March 1997     

The Interaction of Protein Structure, Selection, and Recombination on the Evolution of the Type-1 Fimbrial Major Subunit (fimA) from Escherichia coli

Fimbrial adhesins allow bacteria to interact with and attach to their environment. The bacteria possibly benefit from these interactions, but all external structures including adhesins also allow bacteria to be identified by other organisms. Thus adhesion molecules might be under multiple forms of selection including selection to constrain functional interactions or evolve novel epitopes to avoid recognition. We address these issues by studying genetic diversity in the Escherichia coli type-1 fimbrial major subunit, fimA. Overall, sequence diversity in fimA is high (π= 0.07) relative to that in other E. coli genes. High diversity is a function of positive diversifying selection, as detected by d _N/d _S ratios higher than 1.0, and amino acid residuces subject to diversifying selection are nonrandomly clustered on the exterior surface of the peptide. In addition, McDonald and Kreitman tests suggest that there has been historical but not current directional selection at fimA between E. coli and Salmonella. Finally, some regions of the fimA peptide appear to be under strong structural constraint within E. coli, particularly the interior regions of the molecule that is involved in subunit to subunit interaction. Recombination also plays a major role contributing to E. coli fimA allelic variation and estimates of recombination (2N _e c) and mutation (2N _eμ) are about the same. Recombination may act to separate the diverse evolutionary forces in different regions of the fimA peptide. Received: 13 April 2000 / Accepted: 28 October 2000     

Multiple Mechanosensitive Ion Channels from Escherichia coli,Activated at Different Thresholds of Applied Pressure

Mechanosensitive ion channels from Escherichia coli were studied in giant proteoliposomes reconstituted from an inner membrane fraction, or in giant round cells in which the outer membrane and the cell wall had been disrupted by a lysozyme-EDTA treatment and a mild osmotic shock. Patch-clamp experiments revealed the presence in these two preparations of an array of different conductances (100 to 2,300 pS in 0.1 m KCl) activated by stretch. The electrical activity induced by stretch in the native membrane was complex, due to the activation of several different conductances. In contrast, patches of proteoliposomes generally contained clusters of identical conductances, which differed from patch to patch. These experiments are consistent with the notion that these different conductances correspond to different proteins in the plasma membrane of E. coli, which segregate into clusters of identical channels on dilution involved in reconstitution in proteoliposomes. These conductances could be grouped into three subfamilies of poorly selective channels. In both preparations, the higher the conductance, the higher was the negative pressure needed for activation. We discuss the putative role of these channels as parts of a multicomponent osmoregulatory system. Received: 23 May 1995/Revised: 31 January 1996     

Evolution of Orthologous Intronless and Intron-Bearing Globin Genes in Two Insect Species

While globin genes ctt-2β and ctt-9.1 in Chironomus thummi thummi each have a single intron, all of the other insect globin genes reported so far are intronless. We analyzed four globin genes linked to the two intron-bearing genes in C. th. thummi. Three have a single intron at the same position as ctt-2β and ctt-9.1; the fourth is intronless and lies between intron bearing genes. Finally, in addition to its intron, one gene (ctt-13RT) was recently interrupted by retrotransposition. Phylogenetic analyses show that the six genes in C. th. thummi share common ancestry with five globin genes in the distantly related species C. tentans, and that a 5-gene ancestral cluster predates the divergence of the two species. One gene in the ancestral cluster gave rise to ctn-ORFB in C. tentans, and duplicated in C. th. thummi to create ctt-11 and ctt-12. From parsimonious calculations of evolutionary distances since speciation, ctt-11, ctt-12, and ctn-ORFB evolved rapidly, while ctn-ORFE in C. tentans evolved slowly compared to other globin genes in the clusters. While these four globins are under selective pressure, we suggest that most chironomid globin genes were not selected for their unique function. Instead, we propose that high gene copy number itself was selected because conditions favored organisms that could synthesize more hemoglobin. High gene copy number selection to produce more of a useful product may be the basis of forming multigene families, all of whose members initially accumulate neutral substitutions while retaining essential function. Maintenance of a large family of globin genes not only ensured high levels of hemoglobin production, but may have facilitated the extensive divergence of chironomids into as many as 5000 species. Received: 31 December 1996 / Accepted: 16 May 1997     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998     

Reactivity of Cysteines in the Transmembrane Region of the Na,K-ATPase α Subunit Probed with Hg2+

To gain insight into the structure and conformational coupling in the Na,K-ATPase, this study characterized the reaction of the α1 subunit transmembrane cysteines with a small probe. Intact HeLa cells expressing heterologous Na,K-ATPase were treated with (μm) HgCl₂ after placing the enzyme predominantly in either of two conformations, phosphorylated E2P.Na/E2P or dephosphorylated ATP.E1.K/ATP.E1. Under both conditions the treatment led to enzyme inactivation following a double exponential kinetic as determined by ouabain-sensitive K⁺ uptake measurements. However, the rate constant of the slow reacting component was ten times larger when the protein was probed in a medium that would favor enzyme phosphorylation. Enzymes carrying mutations of cysteines located in the α1 subunit transmembrane region were used to identify the reacting–SH groups. Replacement Cys104Ser reduced enzyme inactivation by removing the slow reacting component under both treatment conditions. Replacement of Cys964 reduced the inactivation rate constant of the fast reacting component (79%) and removed the slow reacting component when the dephosphorylated enzyme was treated with Hg²⁺. Moreover, Cys964Ser substituted enzyme was insensitive to Hg²⁺ when treated under phosphorylation conditions. These results indicate that Cys964 is involved in the fast inactivation by Hg²⁺. Although the double mutant Cys964, 104Ser was still partially inactivated by treatment under nonphosphorylating conditions, an enzyme devoid of transmembrane cysteines was insensitive to Hg²⁺ under all treatment conditions. Thus, this enzyme provides a background where accessibility of engineered transmembrane cysteines can be tested. Received: 13 March 2000/Revised: 23 June 2000     

Detection of Genes with Atypical Nucleotide Sequence in Microbial Genomes

Along the gene, nucleotides in various codon positions tend to exert a slight but observable influence on the nucleotide choice at neighboring positions. Such context biases are different in different organisms and can be used as genomic signatures. In this paper, we will focus specifically on the dinucleotide composed of a third codon position nucleotide and its succeeding first position nucleotide. Using the 16 possible dinucleotide combinations, we calculate how well individual genes conform to the observed mean dinucleotide frequencies of an entire genome, forming a distance measure for each gene. It is found that genes from different genomes can be separated with a high degree of accuracy, according to these distance values. In particular, we address the problem of recent horizontal gene transfer, and how imported genes may be evaluated by their poor assimilation to the host's context biases. By concentrating on the third- and succeeding first position nucleotides, we eliminate most spurious contributions from codon usage and amino-acid requirements, focusing mainly on mutational effects. Since imported genes are expected to converge only gradually to genomic signatures, it is possible to question whether a gene present in only one of two closely related organisms has been imported into one organism or deleted in the other. Striking correlations between the proposed distance measure and poor homology are observed when Escherichia coli genes are compared to Salmonella typhi, indicating that sets of outlier genes in E. coli may contain a high number of genes that have been imported into E. coli, and not deleted in S. typhi. Received: 16 January 2001 / Accepted: 30 August 2001     

Effect of Mutation of Potassium-Efflux System,KefA, on Mechanosensitive Channels in the Cytoplasmic Membrane of Escherichia coli

The effect of a kefA mutation on the mechanosensitive channels in the cytoplasmic membrane of Escherichia coli was established by introducing a mutation of the kefA gene into wild-type E. coli by P1 transduction. The mutation of the kefA gene not only made the cells sensitive to K⁺ in the medium but also changed the mechanosensitive channel activity. The kefA mutation did not change the conductances of the two mechanosensitive channels in the cytoplasmic membrane of E. coli, but it prolonged the channel open time. Also, the kefA mutation made the cells more sensitive to pressure in comparison to wild-type cells. The high sensitivity to pressure of the kefA mutant was not modulated by betaine or by the potassium gradient across the membrane. The effect of the kefA mutation on mechanosensitive channels was not due to a membrane fluidity change. KefA might be a regulator for mechanosensitive channels. Received: 6 September 1995/Revised: 13 December 1995     

Evidence for horizontal transfer from streptococcus to escherichia coli of the kfid gene encoding the K5-specific UDP-glucose dehydrogenase

Capsular polysaccharides are important virulence factors both in Gram-positive and Gram-negative bacteria. A similar cluster organization of the genes involved in the synthesis of bacterial exopolysaccharides has been postulated in both cases, suggesting that these clusters evolved by module assembly. Horizontal gene transfer has been postulated to explain the polymorphism found in these cellular polymers. The cap1K and cap3A genes coding for the pneumococcal type 1 and type 3 UDP-glucose dehydrogenases, respectively, have been compared with other UDP-sugar dehydrogenases. We have observed that the evolutionary distance between Cap1K and Cap3A is approximately equal to that found between Cap1K (or Cap3A) and other UDP-GlcDH of families evolutionarily distant like KfiD, the dehydrogenase from Escherichia coli K5. On the basis of comparisons of G + C content, patterns of synonymous and nonsynonymous substitutions, dinucleotide frequencies, and codon usage bias, we conclude that the kfiD gene has been introduced into E. coli from an exogenous source, probably from a streptococcal species. Received: 26 May 1997 / Accepted: 30 July 1997     

Effects of Density and Gating of Delayed-Rectifier Potassium Channels on Resting Membrane Potential and its Fluctuations

The aim of this study is to evaluate directly, using a reduced experimental system, the nature of interactions between voltage-gated potassium channels and the resting membrane potential. Xenopus oocytes were injected with various concentrations of cRNA coding for a delayed-rectifier potassium channel Shaker-IR. The effects of the density and kinetics of the expressed channels on resting membrane potential is explored in isolated (``inside-out') patches. The channel density is given in terms of maximal conductance (G _max), measured from the maximal slope of the I-V curve under voltage clamp conditions. The capacitance of the experimental setup is approximately 1 pF. At high channel densities (G _max > 10 pA/mV) the mean membrane potential is stabilized at approximately −60 mV. This resting membrane potential is more than 35 mV positive to the reversal potential for potassium ions under the same experimental conditions. Analyses of voltage clamp experiments indicate that at high channel densities the mean membrane potential is determined by the rates of channel activation and deactivation, but is not affected by the rates involved in the process of slow (C-type) inactivation. In contrast, at lower channel densities membrane potential is very unstable, and its mean value and amplitude of fluctuations are strongly affected by the process of slow (C-type) inactivation. Received: 21 March 1996/Revised: 6 August 1996     

N-Carbamoyl Amino Acid Solid–Gas Nitrosation by NO/NO x : A New Route to Oligopeptides via α-Amino Acid N-Carboxyanhydride. Prebiotic Implications

α-N-Carbamoyl amino acid (CAA), whose conditions of formation in a prebiotic hydrosphere have been described previously (Taillades et al. 1998), could have been an important intermediate in prebiotic peptide synthesis through reaction with atmospheric NO _x . Nitrosation of solid CAA (glycine or valine derivative) by a 4/1 NO/O₂ gaseous mixture (1 atm) yields N-carboxyanhydride (NCA) quantitatively in less than 1 h at room temperature. The crude solid NCA undergoes quantitative oligomerization (from trimer to nonamer under the conditions we used) when treated with a (bi)carbonate aqueous buffer at pH 9. We therefore suggest that part of the prebiotic amino acid activation/polymerization process may have taken place in a dry phase (``drying-lagoon' scenario). Received: 23 June 1998 / Accepted: 7 December 1998     

Genetic Characterization of Plasmids Containing Genes Encoding Enzymes of Leucine Biosynthesis in Endosymbionts (Buchnera) of Aphids

The prokaryotic endosymbionts (Buchnera) of aphids are known to provision their hosts with amino acids that are limiting in the aphid diet. Buchnera from the aphids Schizaphis graminum and Diuraphis noxia have plasmids containing leuABCD, genes that encode enzymes of the leucine biosynthetic pathway, as well as genes encoding proteins probably involved in plasmid replication (repA1 and repA2) and an open reading frame (ORF1) of unknown function. The newly reported plasmids closely resemble a plasmid previously described in Buchnera of the aphid Rhopalosiphum padi [Bracho AM, Martínez-Torres D, Moya A, Latorre A (1995) J Mol Evol 41:67–73]. Nucleotide sequence comparisons indicate conserved regions which may correspond to an origin of replication and two promoters, as well as inverted repeats, one of which resembles a rho-independent terminator. Phylogenetic analyses based on amino acid sequences of leu gene products and ORF1 resulted in trees identical to those obtained from endosymbiont chromosomal genes and the plasmid-borne trpEG. These results are consistent with a single evolutionary origin of the leuABCD-containing plasmid in a common ancestor of Aphididae and the lack of plasmid exchange between endosymbionts of different aphid species. Trees for ORF1 and repA (based on both nucleotides and amino acids) are used to examine the basis for leu plasmid differences between Buchnera of Thelaxes suberi and Aphididae. The most plausible explanation is that a single transfer of the leu genes to an ancestral replicon was followed by rearrangements. The related replicon in Buchnera of Pemphigidae, which lacks leuABCD, appears to represent the ancestral condition, implying that the plasmid location of the leu genes arose after the Pemphigidae diverged from other aphid families. This conclusion parallels previously published observations for the unrelated trpEG plasmid, which is present in Aphididae and absent in Pemphigidae. Recruitment of amino acid biosynthetic genes to plasmids has been ongoing in Buchnera lineages after the infection of aphid hosts. Received: 9 March 1998 / Accepted: 18 May 1998.     

Coding sequence divergence between two closely related plant species: Arabidopsis thaliana and Brassica rapa ssp. pekinensis

To characterize the coding-sequence divergence of closely related genomes, we compared DNA sequence divergence between sequences from a Brassica rapa ssp. pekinensis EST library isolated from flower buds and genomic sequences from Arabidopsis thaliana. The specific objectives were (i) to determine the distribution of and relationship between K _a and K _s, (ii) to identify genes with the lowest and highest K _a:K _s values, and (iii) to evaluate how codon usage has diverged between two closely related species. We found that the distribution of K _a:K _s was unimodal, and that substitution rates were more variable at nonsynonymous than synonymous sites, and detected no evidence that K _a and K _s were positively correlated. Several genes had K _a:K _s values equal to or near zero, as expected for genes that have evolved under strong selective constraint. In contrast, there were no genes with K _a:K _s >1 and thus we found no strong evidence that any of the 218 sequences we analyzed have evolved in response to positive selection. We detected a stronger codon bias but a lower frequency of GC at synonymous sites in A. thaliana than B. rapa. Moreover, there has been a shift in the profile of most commonly used synonymous codons since these two species diverged from one another. This shift in codon usage may have been caused by stronger selection acting on codon usage or by a shift in the direction of mutational bias in the B. rapa phylogenetic lineage.