Molecular Biology and Evolution 18:1789-1799 (2001)
© 2001 Society for Molecular Biology and Evolution
Ongoing Evolution of Strand Composition in Bacterial Genomes
Atelier de BioInformatique, Université Paris VI, Paris, France
Unité REG, URA 2171, Institut Pasteur, Paris, France;
HKU-Pasteur Research Centre, Pokfulam, Hong Kong
 |
Abstract |
|---|
 TOP Abstract IntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
We tried to identify the substitutions involved in the establishment of replication strand bias, which has been recognized as an important evolutionary factor in the evolution of bacterial genomes. First, we analyzed the composition asymmetry of 28 complete bacterial genomes and used it to test the possibility that asymmetric deamination of cytosine might be at the origin of the bias. The model showed significant correlation to the data but left unexplained a significant portion of the variance and indicated a systematic underestimation of GC skews in comparison with TA skews. Second, we analyzed the substitutions acting on the genes from five fully sequenced Chlamydia genomes that had not suffered strand switch since speciation. This analysis showed that substitutions were not at equilibrium in Chlamydia trachomatis or in C. muridarum and that strand bias is still an on-going process in these genes. Third, we identified substitutions involved in the adaptation of genes that had switched strands after speciation. These genes adapted quickly to the skewed composition of the new strand, mostly due to C
T, AG, and CG asymmetric substitutions. This observation was reinforced by the analysis of genes that switched strands after divergence between Bacillus subtilis and B. halodurans. Finally, we propose a more extended model based on the analysis of the substitution asymmetries of Chlamydia. This model fits well with the data provided by bacterial genomes presenting strong strand bias. |
Introduction |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
Mutational pressures leading to dramatic differences between the nucleotide compositions of genomes have long been recognized among bacteria (Sueoka 1962
). These patterns produce heterogeneities in the chromosome, either because there is horizontal transfer from genomes with very different compositions (Ochman, Lawrence, and Groisman 2000
) or because the mechanisms causing the bias act differentially in certain regions of the chromosome (Gautier 2000
). Such is the case for replicating strand-associated biases. Under no-strand bias conditions, one would expect to find the equalities A = T and C = G in each strand of DNA (Lobry 1995
). However, the analysis of complete bacterial genomes reveals an important asymmetry between the leading and the lagging replicating strands (Lobry 1996
). These observations have been confirmed in many bacteria and have allowed for the determination of the putative origins of replication in a substantial number of bacterial genomes (Grigoriev 1998
; Salzberg et al. 1998
; Lopez et al. 1999
; Rocha, Danchin, and Viari 1999a
). Interestingly, all genomes present the same asymmetry, G is relatively more abundant than C in the genes coded in the leading strand and less abundant in the genes coded in the lagging strand, which is frequently accompanied by a larger abundance of T over A in the leading strand. These biases propagate into higher-order biases in a correlated way, thereby changing the relative frequencies of codons and amino acids of genes and corresponding proteins in each of the replicating strands (Perrière, Lobry, and Thioulouse 1996
; McInerney 1998
; Rocha, Danchin, and Viari 1999a
).
Bacterial genes transcribed at high levels have a preference for positioning into the leading strand, presumably to minimize collisions between the replication fork and the transcription bubble (McLean, Wolfe, and Devine 1998
). As a consequence, more genes are coded from the leading strand than from the lagging strand in most bacteria (in Gram-positive bacteria, this may go up to 80% of the genes in the leading strand) (Rocha et al. 2000
). Since protein-coding sequences do not contain the same abundance of G and C, this bias accumulates with the replication bias (Sueoka 1999
). Methods designed to deal with this problem revealed that replication strand bias was not due to the asymmetric distribution of genes in the bacterial chromosomes (Rocha, Danchin, and Viari 1999a
; Mackiewicz et al. 1999
).
Explanations for strand biases as a by-product of mechanisms other than replication also include transcription-coupled repair, codon usage bias, and oligonucleotide bias (Francino et al. 1996
; Mrázek and Karlin 1998
; Salzberg et al. 1998
). The former are related to the asymmetrical distribution of highly expressed genes along the two DNA strands. However, the data on codon usage bias (Moszer, Rocha, and Danchin 1999
) and on expression arrays (Tao et al. 1999
) seem to indicate that only a reduced number of genes are highly expressed at exponential growth. Indeed, differences in transcription levels have not been found to constitute a major cause of replication-linked bias (Tillier and Collins 2000a
). Finally, the contribution of signals such as the 
sequence to strand bias was found to be very small due to the fraction of the genome they occupy (Tillier and Collins 2000a
).
Among the theories aimed at explaining strand bias based on the asymmetry of the replication bubble, the cytosine deamination theory enjoys the most attention (Frank and Lobry 1999
). The deamination of cytosine in DNA occurs at significant rates in vivo and leads to the formation of uracil, which is excised by the action of uracil-DNA glycosylase (Lindahl 1993
). The rate of cytosine deamination increases by a factor of 140 when the DNA is single-stranded (Beletskii and Bhagwat 1996
). Methylation of cytosine is known to increase the rate of deamination by a factor of 4. In this case, the deamination produces a T, which cannot be corrected by the glycosylase (Coulondre et al. 1978
). Since the leading strand is more exposed in the single-stranded state (in order to serve as template for the synthesis of the lagging strand) (Marians 1992
), and CT mutations would induce the formation of GC skews, cytosine deamination has been proposed to be at the basis of strand bias (Frank and Lobry 1999
). This hypothesis has the advantage of explaining GC and TA skews (and larger GC skews in G+C-poor genomes) within a known mechanistic model and based on a well-established mutational hot spot.
Recently, a large number of studies have accounted for strand asymmetries (for reviews, see Francino and Ochman 1997
; Frank and Lobry 1999
). However, several important questions still remain unanswered: (1) Are the genomes that present strong strand bias at compositional equilibrium? (2) What are the major substitutions associated with the establishment of the bias? (3) Is there a simple specific function/mutation responsible for the bias? To tackle these questions, we benefited from the existence of the complete genome sequences of several very closely related bacteria, in particular five Chlamydia genomes (Stephens et al. 1998
; Kalman et al. 1999
; Read et al. 2000
; Shirai et al. 2000
) and two species of Bacillus (Kunst et al. 1997
; Takami et al. 2000
).
Asymmetrical changes are usually studied either by phylogenetic reconstruction of homologous sequences or by detection of deviations from the parity of bases in the genome. In the present work, we explored both types of methodologies, taking advantage of the very strong strand bias of Chlamydia and Bacillus genomes and their extensive homology.
|
Materials and Methods |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
Data
Sequence data for all complete bacterial genomes were retrieved from GenBank (http://www.ncbi.nlm.nih.gov). We analyzed the following complete genomes, using the annotations contained in their respective GenBank files: Aquifex aeolicus (Aae), Bacillus halodurans (Bha), B. subtilis (Bsu), Borrelia burgdorferi (Bbu), Buchnera sp. (Bsp), Caulobacter crescentus (Ccr), Campylobacter jejuni (Cje), Chlamydia pneumoniae CWL029 (CpnC), C. pneumoniae AR39 (CpnA), C. pneumoniae J130 (CpnJ), C. trachomatis serovar D (Ctr), C. muridarum (Cmu), Deinococcus radiodurans (two chromosomes) (Dra), Escherichia coli (Eco), Haemophilus influenzae (Hin), Helicobacter pylori 26695 (Hpy), Mycoplasma genitalium (Mge), M. pneumoniae (Mpn), Mycobacterium tuberculosis (Mtu), Neisseria meningitidis MC58 (Nme), Pseudomonas aeruginosa (Pae), Rickettsia prowazekii (Rpr), Synechocystis spp. C125 (Ssp), Treponema pallidum (Tpa), Thermotoga maritima (Tma), Ureaplasma urealyticum (Uur), Vibrio cholerae (two chromosomes) (Vch), and Xylella fastidiosa (Xfa).
Statistical Analysis of Skews
Identification of Biased Genomes Through Linear Discriminant Analysis
Simple cumulative GC and TA skews are sensitive to different populations of genes in the two replicating strands. Therefore, we used linear discriminant analysis to identify genomes with significant strand bias (Rocha, Danchin, and Viari 1999a
). We considered a genome to contain a significant strand bias when the maximal accuracy of our method (percentage of true positives in the classification of leading- and lagging-strand genes) was better than that of the best of 10 random genomes of the same size and composition. In practice, this implies that a genome displays significant strand bias if the accuracy of the discrimination between the replicating strands is larger than 60%–75% (depending on genome length). Because different bacterial strains within a species are often very similar in sequence, and in order not to bias the results, we analyzed only one representative strain for each bacterial species.
Identifying Origins of Replication by GC Skew
The skews for a given sequence were defined as (Lobry 1996
)
|
|
).
Genes' GC Skews
We computed GC and TA skews for the gene sequences using equations (1) and (2) . To compare the differences in skews between the genes present in the different strands, we computed 
GC and TA skews. These quantities are defined as the difference between the average skews of the genes in the leading strand and the ones in the lagging strand. Considering Nleading and Nlagging, the numbers of genes in the leading and lagging strands, respectively, one obtains
|
|
Analysis of Similarity
Definition of Homologous Genes
Two genes were considered homologous if they coded for proteins similar both in sequence and in size. To identify homologous genes, we performed pairwise comparisons of all proteins of all proteome pairs, filtering potential hits with P < 10-5 in BlastP and a maximal difference of protein lengths of 20%. Subsequently, we aligned the sequences using a variant of the classical dynamic programming algorithm for global alignment, where one counts 0-weight for gaps at both ends of the largest sequence using the BLOSUM62 matrix (Erickson and Sellers 1983
). Finally, we retained pairs of proteins with more than 40% similarity.
Classification of Orthologous Genes
Two homologous genes were considered to be orthologous if they were each other's best matches in the respective genomes. We obtained 687 sets of five orthologous for the Chlamydia set and 2,123 sets of two orthologous for the Bacilli set. Using these sets, we analyzed the conservation of the gene organization between genomes by displaying a scatter-plot of the positions of orthologous genes in the different genomes. We further defined three classes of genes for Chlamydia and for Bacillus: genes present in all genomes in the same replicating strand (SS), genes present in different replicating strands (DS), and genes not present in other genomes according to our stringent orthology criteria (NS). For Chlamydia, this resulted in 49 DS and 638 SS genes, whereas for Bacillus we obtained 372 DS and 1751 SS genes. Naturally, the number of NS genes changed from species to species. For example, we obtained 1,977 NS genes for B. subtilis and 131 NS genes for C. muridarum.
Characterization of Orthologous Genes
The Chlamydia set was extremely interesting due to the small divergence between the five genomes. The similarity of the ribosomal 16S subunit was 93.5% between C. trachomatis and C. pneumoniae, 97.4% among C. trachomatis and C. muridarum genomes, and >99% among C. pneumoniae strains. This is in accordance with the phylogeny of Chlamydia that proposes a first speciation event between C. pneumoniae and the pair C. trachomatis/C. muridarum (Everett, Bush, and Andersen 1999
). The divergence between the two Bacillus 16S subunits (94.4%) was intermediate to that between different species of Chlamydia. A further advantage of the Chlamydia set was the intermediate level of synteny between the elements: within C. pneumoniae and within the pair C. trachomatis/C. muridarum, there was complete conservation of gene order, whereas between these two sets there was less conservation. As a result, the classification of orthologous genes was valid for the genomes of C. pneumoniae on one side and for C. trachomatis and C. muridarum on the other. Also, the vast majority of NS genes were present either in all C. pneumoniae strains or in both C. trachomatis and C. muridarum. Since Chlamydia and Bacillus contain strand biases, one expects a strong and significant signal in these genomes. Finally, we can neglect the effects of differences in G+C content, because they were similar within the two groups of genomes (table 1
).
|
Characterization of Differences in Alignments
Directed Changes
We made multiple alignments of protein and DNA sequences of SS and DS genes using CLUSTAL W with default parameters (Thompson, Higgins, and Gibson 1994
). Multiple alignments of proteins were back-translated into DNA in order to determine rates of synonymous (Ks) and nonsynonymous (Ka) substitutions (Li 1997
). Since these values are sensitive to low levels of sequence identity, we took into consideration only pairs of genes with Ks < 2.0 for the analysis of SS genes. We used this set of DNA alignments to analyze the spectrum of changes observed among the orthologous genes. For this analysis, we invoked parsimony by counting a change only when a column of the multiple alignments presented a nucleotide consensus minus one (e.g., one column with four G's and one C is a GC change in the sequence with C). Since each alignment corresponded to a set of SS genes whose positions in the chromosome were known, we could separate the alignments into two subsets; one corresponding to all genes in the leading strand (346 genes) and the other corresponding to genes in the lagging strand (292 genes). Consequently, we analyzed the changes occurring in each strand separately. For these analyses, we used JaDis, a publicly available program designed to compute distances between nucleic acid sequences (http://pbil.univ-lyon1.fr/) (Gonçalves et al. 1999
). We performed a similar analysis for DS genes. In both cases, the absolute values of substitutions from i to j were converted to relative substitution frequencies fij by dividing by the average number of nucleotides i in the sequence and expressing the result as a percentage among all types of nucleotide substitution (Gojobori, Li, and Graur 1982
).
Undirected Changes
Since the switch of replicating strand took place before the divergence C. trachomatis/C. muridarum or of C. pneumoniae strains, we used a complementary approach to analyze DS genes, using pairwise alignments to count mismatches. This analysis should account for the adaptation of genes since switching strands, whereas the analysis of DS genes with multiple alignments only accounts for adaptation since speciation. To limit our analysis to well-conserved proteins, we imposed a stricter threshold of similarity (>60% in protein sequence), which resulted in a set of proteins exhibiting an average similarity of 75%. These proteins were more constrained at the amino acid level, which renders pairwise alignments more reliable. For the two bacilli, we lacked a sufficiently close outgroup to determine the direction of mutations, and the same method was used: (1) We aligned the sequences and identified the mismatches. (2) Having arbitrarily chosen one of the two species as the reference species (RS), we cataloged all mismatches in classes XRS : YnonRS (X 
Y).(3) We separated all DS genes into two categories: the genes that were in the leading strand in the reference species (lagging in the other genome), and those that were in the lagging strand in the reference species (leading in the other genome). (4) For each of these two categories, we computed the difference between XRS : YnonRS and YRS : XnonRS, which indicates the asymmetry in terms of mismatches. The comparison of the asymmetries between the two categories indicates the asymmetry at the basis of the adaptation to the new strand.
It is important to emphasize the differences between this and the preceding analysis. Suppose that a gene is in the leading strand of the reference species. In the analysis of the multiple alignment of SS genes, a mismatch C4,lead : T1,lead indicates a CleadTlead substitution. In the pairwise alignment of DS genes, a mismatch Tlead : Clag cannot be interpreted as a CleadTlead change. This mismatch may have been caused by CT in the leading strand of the reference genome or by TC in the lagging strand of the other genome (in the previous analysis, the probability of four TleadClead changes was neglected, referring to parsimony). Therefore, C : T mismatches in DS genes indicate either a CleadTlead change or an AleadGlead change. We indicate this ambiguity by C(A)T(G).
Further Analysis of Alignments
To analyze the importance of multiple substitutions in our data set, we looked for changes in the SS multiple alignments between C. trachomatis and C. muridarum (which corresponded to a larger evolutionary distance than the one between C. pneumoniae strains). The average Ks was below one substitution per synonymous site, and the upper Ks value for the 95% confidence interval was 0.84 (for nonsynonymous sites [Ka] it was 0.07).
|
Results |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
Strand Biases Among Bacteria
Prevalence of Strand Biases
Twenty-one out of 28 chromosomes (representing 26 species) of bacteria had significant strand biases (table 1 ). Exceptions were the Mycoplasmas, D. radiodurans, A. aeolicus, and Synechocystis sp., many of which have been described previously (Rocha, Danchin, and Viari 1999a
). The largest strand bias (as measured by maximal discrimination) was observed in the genome of Borrelia burgdorferi, followed by C. trachomatis, T. pallidum, C. pneumoniae, and chromosome 2 of V. cholerae. Although many of these highly biased genomes presented low G+C contents, the correlation between accuracy and G+C content was not significantly different from 0 (Spearman's rho = -0.18, P > 0.3). In fact, the above-mentioned genomes lacking strand bias (with the exception of D. radiodurans) also had G+C contents well below 50%.
Correlation Between Discrimination and Skews
GC skews of the genomes were closely correlated with the accuracy of discrimination between strands. This was less evident for the TA skews. In fact, the Spearman's rank correlation between maximal discrimination and GC skew was 0.77 (P < 0.001), but it was only 0.40 between maximal discrimination and TA skew (P < 0.1). Because the cytosine deamination theory predicts proportionality between the two skews, we built a model for the theory and tested it with the data of biased genomes.
Testing the Cytosine Deamination Theory
The composition of a hypothetical unbiased gene in terms of the four bases is NA,0, NC,0, NG,0, NT,0, taken as the mean of the average composition of the leading- and lagging-strand genes. Suppose that the gene is in the leading strand and will suffer the corresponding asymmetry. The deamination theory predicts that strand bias will induce a CT change with probability z for each C during a period t. Therefore, the composition of the gene remains unchanged for G and A, but not for C and T:
|
|
|
|
T asymmetries to predict coherent values for GCt and TAt within each genome. Because the deamination theory is a mutational theory, the best fit of the model should correspond to the analysis of third positions of codons. Note that this model refers to CT asymmetries, of which cytosine deamination is only one particular case.
The 21 bacterial chromosomes with strand bias showed a significant correlation between the two values of zt predicted by the two expressions in equation (6)
using third positions of codons. In fact, after excluding Buchnera sp. (see below), the Spearman's rank correlation between the two terms was 0.83 (P < 0.001), suggesting an important contribution of CT asymmetries to the establishment of the strand bias. However, a systematic underevaluation of GC skew was also apparent (P < 0.05, Wilcoxon test). In 17 out of 20 chromosomes, the number of changes required to explain the TA skews was larger than the one required to explain the GC skew. The systematic underestimation of GC skews suggests the existence of other sources of bias.
Preliminary regression analysis of the zt data indicated heteroscedasticity, with variance increasing with the values of zt. Standard procedures of regression analysis recommend a logarithmic transformation of the data in such cases (Zar 1996
). The linear regression on the transformed data (excluding Buchnera sp. and H. pylori) resulted in a fit presenting a coefficient of determination of 0.68 (P < 0.001). Although the regression line fits the data well, it does not fit the expected result (fig. 1
). This confirms the results of the Wilcoxon test, suggesting that a model based solely on CT asymmetries underestimates GC skews.
|
Strand Bias Oscillations in Chlamydia
Chlamydia SS Genes Are Not at Equilibrium
We identified the recent substitution asymmetries in the SS genes of Chlamydia by analyzing the changes in C. trachomatis, C. muridarum, and C. pneumoniae. For C. pneumoniae, we considered all changes observed in the genomes of the three strains to obtain larger counts. Because these strains diverged rather recently, the parameters of the substitution matrix had very large confidence intervals. Indeed, our analysis used only 0.1% of the positions in the multiple alignments for C. pneumoniae, compared with 6.0% for C. trachomatis and C. muridarum. The frequencies of substitutions for the other Chlamydia are presented in table 2 . Most changes are not significantly unbalanced between replicating strands in C. trachomatis and C. muridarum, and none are statistically significant in C. pneumoniae. In C. muridarum, we observed that SS genes from both strands were getting richer in G and T and poorer in A and C (P < 0.001; Wilcoxon tests). In C. trachomatis, SS genes were getting richer in G and C and poorer in A and T (P < 0.001; Wilcoxon tests). When we compared the evolutions of the relative compositions in the two strands, we observed that the C. muridarum genes in the leading strand were getting richer in G and C and poorer in A and T (P < 0.02; Wilcoxon tests). In C. trachomatis, the leading strand genes were getting richer in T and G (P < 0.01; Wilcoxon tests) when compared with the lagging strand (differences for A and C are not statistically significant). As a consequence, there was no significant evolution of the GC and TA skews in the genes of the two strands for C. muridarum, and there was an increase in GC skew in C. trachomatis (P < 0.01; Wilcoxon test), with no significant differences for TA skews. For C. pneumoniae, the differences were not statistically significant. The most important asymmetries in substitution frequencies in C. muridarum and C. trachomatis were A
G–GA (2.64% and 0.86%, respectively) and CT–TC (1.45% and -0.87%, respectively), but only the former were statistically significant (P < 0.01; Wilcoxon test).
|
The Fate of Inverted Genes and the Evolution of
GC SkewsInversions in Chlamydia
Inverted genes in genomes with strand bias are expected to adapt fast to the composition of the new strand (Rocha, Danchin, and Viari 1999a
; Tillier and Collins 2000b
). Hence, we computed GC skews in third codon positions of DS genes in C. pneumoniae AR39 and C. muridarum (orthologous in different replicating strands) to test if they had evolved toward the bias of the new strand. We observed that all but two DS genes had GC skews typical of the new strand (fig. 2
). We then proceeded to test if DS genes were distinguishable from SS and NS genes of the same current strand. Taking C. muridarum as a reference, we observed that GC skew at position 3 was significantly different between the three types of genes (P < 0.001; Kruskal-Wallis test), but, surprisingly, DS genes possessed a larger GC skew (0.567) than SS (0.492) or NS (0.54) genes; i.e., the order of the bias was DS 
NS > SS (P < 0.05; Tukey-Kramer test). The same analysis applied to C. pneumoniae AR39 genes revealed similar results, but in this case the order was NS 
DS > SS (P < 0.01; Tukey-Kramer test).
|
Inversions in Bacillus
A similar pattern was found among DS genes in Bacillus. Taking B. subtilis as a reference, we observed
GC skews in the order NS (0.171) DS (0.167) > SS (0.124) (P < 0.05; Tukey Kramer test), and similar results were obtained using B. halodurans as the reference. As in Chlamydia, the genes in Bacillus that have suffered an inversion have acquired the composition corresponding to the new host strand. Also, in both genera, DS genes exhibited biases larger than expected for their new strand.
Characterization of Mutations in Inverted Genes
Within each of the two Chlamydia monophyletic groups the genomes were collinear. Because strand switch took place before the speciation events, the analysis of 34 DS genes was first performed using only a pair of genomes. We chose C. muridarum and C. pneumoniae AR39 for simplicity (both sequences start at the origin of replication). DS genes suffered a strand switch since speciation of these two species, and this induced opposite replication biases (Tillier and Collins 2000b
). Hence, the substitutions that took place in these genes provide clues for the establishment of the bias. We observed significant asymmetries in A(G)C(T)–C(T)A(G), in CG, and especially in C(A)T(G)–T(G)C(A) (table 3
). In this analysis, the directionality of changes could not be determined (see Materials and Methods).
|
These changes led to an increase in G and T in the leading strand and an increase in C and A in the lagging strand. Similar results were obtained for the DS genes in the two bacilli. For both Bacillus and Chlamydia, C and G increased in the respective strands at a faster rate than A and T. This is because C
G–GC is asymmetric, but AT–TA is not. We analyzed the changes in the third codon positions for both sets in order to check if selective constraints at the amino acid level could be responsible for part of the signal. The differences were not statistically significant (P > 0.1; Wilcoxon test).
Although the strand switch took place before speciation events within C. pneumoniae or between C. trachomatis and C. muridarum, one may suppose that some of the change occurred after speciation. In this case, the analysis of DS genes could be done using the multiple alignments, as for SS genes, with the advantage that the direction of the substitutions can be determined. We built the table of relative substitution frequencies for the DS genes (table 4
), which revealed three significantly different frequencies of substitution: CT–TC (11.4%), AG–GA (8.4%). and CG–GC (2.2%) (P < 0.01; Wilcoxon tests). The AC–CA difference (2.4%) had the same magnitude as that of CG-GC, but it was not statistically significant (P > 0.1).
|
Refining the Model
We started by proposing a model based on C
T asymmetries, which were clearly insufficient to explain the data. The analysis of pairwise alignments of DS genes suggested that CG had about half the importance of C(A)T(G)–T(G)C(A) in the establishment of the asymmetry. The analysis of multiple alignments indicated that CT–TC was only slightly more important than AG–GA. Since these values were for Chlamydia and they may change for different species, we make the simplification that the three asymmetries all have the same relative importance. Hence, we can rewrite equation (5)
as
|
|
Genes Evolve at Different Rates Depending on Position and Type
Amounts of Changes in Genes in the Different Replicating Strands
The lagging-strand genes presented more changes among SS genes both for Chlamydia (6.6%, P < 0.001; Wilcoxon test) and for Bacillus (6.4%, P < 0.001).
Different Evolution Rates for NS, SS, and DS Genes
One would expect SS genes to exhibit higher GC skews, but we observed the opposite. Hence, we tested to determine if genes in different strands and belonging to different types (i.e., SS, DS) evolved at similar rates. Both in Chlamydia and in Bacillus, DS genes evolved significantly faster than SS genes (P < 0.001; Wilcoxon test). Between C. trachomatis and C. muridarum, similarity scores increase as follows: NS < DS < SS (P < 0.01; Tukey-Kramer test).
Ka/Ks Ratios
We investigated the Ka/Ks ratios of orthologs of C. trachomatis and C. muridarum and of orthologs of C. muridarum and C. pneumoniae AR39 (note that we removed all pairs for which Ks > 2.0; see Materials and Methods). The median Ka/Ks value for orthologs in C. trachomatis/C. muridarum was 0.07 (0.069 in the lagging strand and 0.073 in the leading strand) and 0.12 for orthologs in C. muridarum/C. pneumoniae AR39 (0.120 in the lagging strand and 0.117 in the leading strand). These differences between leading- and lagging-strand genes are not statistically significant. The comparison of SS and DS genes among bacilli revealed similar values of Ka/Ks for both sets (respectively, 0.12 and 0.14).
|
Discussion |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
Ongoing Strand Bias in Bacteria
Following previous observations (Mackiewicz et al. 1999
; Rocha, Danchin, and Viari 1999a
), we observed pervasiveness of strand bias in bacterial genomes, with GC skews correlating better with strand discrimination than TA skews. This suggests that strand bias should not have its origin only on CT asymmetries, a hypothesis which is reinforced by the analyses of the models and of the asymmetric substitution rates. The extended model for the impact of asymmetry in gene evolution revealed a good fit to the data and one single outlier: M. tuberculosis. Nevertheless, Buchnera sp. is a borderline case (fig. 1
), which may not be surprising given that it has the characteristics of an endocellular symbiont suffering a genome reduction that involves the lost of a considerable amount of repair mechanisms (Shigenobu et al. 2000
). One should note that the species under comparison diverged a long time ago, have very different lifestyles, and contain proteomes heterogeneous in size and composition, different G+C contents, and very different sets of replication and repair machineries. All of these effects result in the considerable amount of variance that remains unexplained by the linear regression analyses.
One might expect that substitution frequencies in SS genes in stable genomes with strong GC skews such as C. muridarum and C. trachomatis might be at equilibrium. Instead, they were found to entail variations in the nucleotide compositions of the two strands. As a result, in C. trachomatis, there is a net increase in GC skew in SS genes (nonsignificant for C. muridarum). This means that, at least in this species, strand bias is still an ongoing process even in SS genes. Hence, it is not surprising that in both Chlamydia and Bacillus, DS genes have adapted to the composition of the new strand. It was previously reported that paralogs in different replicating strands of B. burgdorferi presented signs of adaptation to the respective strands (Lafay et al. 1999
; Rocha, Danchin, and Viari 1999a
), and a recent study of DS genes in Chlamydia demonstrated that they adapt fast to the new strand (Tillier and Collins 2000b
). Our results confirm these previous works and provide some clues for the substitutions at the basis of the adaptation.
Our analysis of SS and DS genes used closely related sequences, but not enough to assure complete avoidance of multiple substitutions. To avoid a large number of multiple substitutions and to produce faithful alignments, we restricted our study to the most conserved proteins. These proteins have evolved less due to functional constraints. We are then observing substitutions for which the weight of selection may be important, which incorporates a bias, particularly due to the codon usage bias and transcription-coupled repair of highly expressed genes. Chlamydia and Bacillus have different codon usage biases (Moszer, Rocha, and Danchin 1999
; Romero, Zavala, and Musto 2000
), but the analysis of the substitutions in the two groups is concordant. If codon usage bias was biasing the results in a very important way, this should not happen. Also, it has been found that biases associated with highly expressed genes do not significantly interfere with replication bias (Mackiewicz et al. 1999
; Tillier and Collins 2000a
), and one may expect this to not significantly change the results. The interference of some biases in the analyses of other biases is a very important question that unfortunately remains unsolved. Analyses of other close genomes with important strand bias may shed further light on this question.
Different Causes for a Simple Bias
The two major asymmetries in DS genes of C. muridarum and C. trachomatis since speciation are AG–GA and CT–TC. This is consistent with the observations of SS genes that also indicate that these two substitutions are the most asymmetric (although the latter is not statistically significant). It is also consistent with the analysis of DS genes using pairwise alignments in order to capture the substitutions during all the processes of adaptation. Unfortunately, the latter analysis does not allow one distinguish between CT and AG. Both asymmetries induce increases in GC and TA skews and lead to proportionality between the increases. CT asymmetries may be assigned to preferential cytosine deamination in single-stranded DNA, although other hypotheses based on CT asymmetries may also be compatible with the data. G : CA : T transitions dominate the spectra of mutations of E. coli; however, most studies have focused on CT mutations, not on AG (Frank and Lobry 1999
), and few data are available on AG mutations.
Two other different substitution frequencies are asymmetric in the adaptation process of DS genes: CG–GC and AC–CA. Both were also observed in the analyses of pairwise alignments, but in this case AC–CA was indistinguishable from GT (although GT was not significant in the multiple-alignment approaches). Interestingly, CG and GC are among the most rare mutations observed in E. coli (Hutchinson 1996
). However, the asymmetry is not necessarily correlated with the absolute number of substitutions. Indeed, also in our data set, CG and GC are systematically the most rare substitutions (tables 2 and 4
). Nevertheless, it is puzzling to observe that the pairwise alignments indicate a stronger role for CG–GC in the adaptation to the strand than do the multiple alignments (7% and 2.2%, respectively). One may consider two explanations for this observation: (1) CG asymmetries are more important in earlier phases of the adaptation, or (2) multiple substitutions in the pairwise alignments are biasing our results concerning these rare mutations (pairwise alignments represent a longer period of evolution). Although the AC–CA asymmetry plays no direct role in the establishment of the bias (it just converts TA skew in GC skew), one may suppose that part of the CG substitutions in fact correspond to CAG multiple substitutions. Independent of its origin, the CG asymmetry results in an increase in GC skew without an increase in TA skew. This may explain the systematic underevaluation of GC skews by the model based only on CT asymmetry: the contribution of CG asymmetry would "correct" this effect (note that AT is not asymmetric).
We have suggested elsewhere that genome shuffling would disturb strand bias (Rocha, Danchin, and Viari 1999b
). However, one might also suppose mutation rates and repair efficiency to play an important role in the tuning of the process. As for mutation, using the example of cytosine deamination, the methylation of cytosine increases CT mutations. In bacteria, such methylation is provided by restriction modification systems, which are constantly being acquired and lost by horizontal transfer (Jeltsch and Pingoud 1996
). Therefore, mutation could be a cause of the change in bias with time. As for repair, several reports indicate that the mismatch repair system compensates for CT substitution asymmetries, even when the cytosine is methylated (Jones, Wagner, and Radman 1987
). A small change in the repair machinery could originate an increased capacity for repair and a consequent reduction in the asymmetry. The efficiency of mismatch repair does change along the evolutionary history of bacteria (Taddei et al. 1997
; Sniegowski et al. 2000
). It has also been shown that replicating strands exhibit different replication accuracy rates (Izuta, Roberts, and Kunkel 1995
; Iwaki et al. 1996
; Fijalkowska et al. 1998
) and that repair might be involved in it (Radman 1998
).
Since the analyses of DS genes were performed between genomes that diverged some time ago, we cannot assume that the substitutions we observe are devoid of selective constraints. Nevertheless, we observed that switched genes adapted fast to the new strand and that this adaptation resulted in a small Ka/Ks ratio. This suggests that amino acid bias is not at the origin of such biases, even if the adaptation process involves changing the amino acid content of the coded proteins (Perrière, Lobry, and Thioulouse 1996
; Lafay et al. 1999
). Effects of amino acid selection on strand bias have also been discarded by other analyses (Mackiewicz et al. 1999
; Tillier and Collins 2000b
). However, selective effects at the DNA level cannot be ruled out through this approach. In particular, one cannot exclude the possibility that strand bias is the result of selection for more efficient replication by the asymmetric replication bubble.
|
Conclusions |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
The existence of strand bias has important consequences for the study of bacterial molecular evolution. First, it indicates that in many bacteria the use of substitution matrices that do not take into account strand asymmetries provides a poor approximation of real data. Second, it indicates that a gene may suffer a process of accelerated evolution just through a change of replicating strand. This may also make the discrimination of paralogy from orthology difficult, especially in large functional families. Third, it provides an additional signal to use in the detection of horizontal transfer and genome rearrangements, but only for recent events. Fourth, it may elucidate still unknown particularities in the processes underlying DNA replication and repair.
|
Acknowledgements |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
Alain Viari played a very important role in earlier discussions of this work. Isabelle Gonçalves programmed in JaDis the function that allows the determination of mutations of the type consensus - 1 in the multiple alignments. We are grateful for comments and suggestions from Elisabeth Tillier, Carmen Gomes, and Isabelle Gonçalves on previous versions of the manuscript. The criticisms and suggestions of two anonymous referees constituted important contributions to this work.
|
Footnotes |
|---|
Howard Ochman, Reviewing Editor
1 Abbreviations: DS, orthologous genes present in different replicating strands (i.e., leading versus lagging); Ka, nonsynonymous substitution rate; Ks, synonymous substitution rate; NS, genes without orthologs in the other species; SS, orthologous genes present in the same replicating strand. 
2 Keywords: replication
strand bias
mutation
genome analysis
sequence evolution
3 Address for correspondence and reprints: Eduardo P. C. Rocha, Atelier de BioInformatique, Université Paris VI, 12 Rue Cuvier, 75005 Paris, France. erocha{at}abi.snv.jussieu.fr
|
References |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionConclusionsAcknowledgementsReferences |
|---|
Beletskii A., A. S. Bhagwat, 1996 Transcription-induced mutations: increase in C to T mutations in the non-transcribed strand during transcription in Escherichia coli Proc. Natl. Acad. Sci. USA 93:13919-13924
Coulondre C., J. H. Miller, P. J. Farabaugh, W. Gilbert, 1978 Molecular basis of base substitution hotspots in Escherichia coli Nature 274:775-780[Medline]
Erickson B. W., P. H. Sellers, 1983 Recognition of patterns in genetic sequences Pp. 55–91 in D. Sankoff and J. B. Kruskal, eds. Time warps, string edits, and macromolecules: the theory and practice of sequence comparison. Addison-Wesley, Reading, Mass
Everett K. D., R. M. Bush, A. A. Andersen, 1999 Emended description of the order Chlamydiales, proposal of Parachlamydiaceae fam. nov. and Simkaniaceae fam. nov., each containing one monotypic genus, revised taxonomy of the family Chlamydiaceae, including a new genus and five new species, and standards for the identification of organisms Int. J. Syst. Bacteriol 49:415-440
Fijalkowska I. J., P. Jonczyk, M. M. Tkaczyk, M. Bialokorska, R. M. Schaaper, 1998 Unequal fidelity of leading strand and lagging strand DNA replication on the Escherichia coli genome Proc. Natl. Acad. Sci. USA 95:10020-10025
Francino M. P., L. Chao, M. A. Riley, H. Ochman, 1996 Asymmetries generated by transcription-coupled repair in enterobacterial genes Science 272:107-109[Abstract]
Francino M. P., H. Ochman, 1997 Strand asymmetries in DNA evolution Trends Genet 13:240-245[Web of Science][Medline]
Frank A. C., J. R. Lobry, 1999 Asymmetric patterns: a review of possible underlying mutational or selective mechanisms Gene 238:65-77[Web of Science][Medline]
Gautier C., 2000 Compositional bias in DNA Curr. Opin. Genet. Dev 10:656-661[Web of Science][Medline]
Gojobori T., W.-H. Li, D. Graur, 1982 Patterns of nucleotide substitution in pseudogenes and functional genes J. Mol. Evol 18:360-369[Web of Science][Medline]
Gonçalves I., M. Robinson, G. Perriere, D. Mouchiroud, 1999 JaDis: computing distances between nucleic acid sequences Bioinformatics 15:424-425
Grigoriev A., 1998 Analyzing genomes with cumulative skew diagrams Nucleic Acids Res 26:2286-2290
Hutchinson F., 1996 Mutagenesis Pp. 2218–2235 in F. Neidhardt, R. Curtiss, J. L. Ingraham, E. C. Lin, K. B. Low, B. Magasanik, W. S. Reznikoff, M. Riley, M. Schaechter, and H. E. Umbarger, eds. Escherichia coli and Salmonella: cellular and molecular biology. ASM Press, Washington, D.C
Iwaki T., A. Kawamura, Y. Ishino, K. Kohno, Y. Kano, N. Goshima, M. Yara, M. Furusawa, H. Doi, F. Imamoto, 1996 Preferential replication-dependent mutagenesis in the lagging DNA strand in Escherichia coli Mol. Gen. Genet 251:657-664[Web of Science][Medline]
Izuta S., J. D. Roberts, T. A. Kunkel, 1995 Replication error rates for G. dGTP, T.dGTP, and A.dGTP mispairs and evidence for differential proofreading by leading and lagging strand DNA replication complexes in human cells J. Biol. Chem 270:2595-2600
Jeltsch A., A. Pingoud, 1996 Horizontal gene transfer contributes to the wide distribution and evolution of type II restriction-modification systems J. Mol. Evol 42:91-96[Web of Science][Medline]
Jones M., R. Wagner, M. Radman, 1987 Mismatch repair of deaminated 5-methyl-cytosine J. Mol. Biol 194:155-159[Web of Science][Medline]
Kalman S., W. Mitchell, R. Marathe, C. Lammel, J. Fan, R. W. Hyman, L. Olinger, J. Grimwood, R. W. Davis, R. S. Stephens, 1999 Comparative genomes of Chlamydia pneumoniae and Chlamydia trachomatis Nat. Genet 21:385-389[Web of Science][Medline]
Kunst F., N. Ogasawara, I. Moszer, et al. (151 co-authors) 1997 The complete genome sequence of the Gram-positive bacterium Bacillus subtilis Nature 390:249-256[Medline]
Lafay B., A. T. Lloyd, M. J. McLean, K. M. Devine, P. M. Sharp, K. H. Wolfe, 1999 Proteome composition and codon usage in spirochaetes: species-specific and DNA strand-specific mutational biases Nucleic Acids Res 27:1642-1649
Li W.-H., 1997 Molecular evolution Sinauer, Sunderland, Mass
Lindahl T., 1993 Instability and decay of the primary structure of DNA Nature 362:709-715[Medline]
Lobry J. R., 1995 Properties of a general model of DNA evolution under no-strand bias conditions J. Mol. Evol 40:326-330[Web of Science][Medline]
———. 1996 Asymetric substitution patterns in the two DNA strands of bacteria Mol. Biol. Evol 13:660-665[Abstract]
Lopez P., H. Philippe, H. Myllykallio, P. Forterre, 1999 Identification of putative chromosomal origins of replication in Archaea Mol. Microbiol 32:883-886[Web of Science][Medline]
McInerney J. O., 1998 Replicational and transcriptional selection on codon usage in Borrelia burgdorferi Proc. Natl. Acad. Sci. USA 95:10698-10703
Mackiewicz P., A. Gierlik, M. Kowalczuk, M. R. Dudek, S. Cebrat, 1999 How does replication-associated mutational pressure influence amino acid composition of proteins? Genome Res 9:409-416
McLean M. J., K. H. Wolfe, K. M. Devine, 1998 Base composition skews, replication orientation and gene orientation in 12 prokaryote genomes J. Mol. Evol 47:691-696[Web of Science][Medline]
Marians K. J., 1992 Prokaryotic DNA replication Annu. Rev. Biochem 61:673-719[Web of Science][Medline]
Moszer I., E. P. C. Rocha, A. Danchin, 1999 Codon usage and lateral gene transfer in Bacillus subtilis Curr. Opin. Microbiol 2:524-528[Web of Science][Medline]
Mrázek J., S. Karlin, 1998 Strand compositional asymmetry in bacterial and large viral genomes Proc. Natl. Acad. Sci. USA 95:3720-3725
Ochman H., J. G. Lawrence, E. A. Groisman, 2000 Lateral gene transfer and the nature of bacterial innovation Nature 405:299-304[Medline]
Perrière G., J. R. Lobry, J. Thioulouse, 1996 Correspondence discriminant analysis: a multivariate method for comparing classes of protein and nucleic acid sequences CABIOS 12:519-524
Radman M., 1998 DNA replication: one strand may be more equal Proc. Natl. Acad. Sci. USA 95:9718-9719
Read T. D., R. C. Brunham, C. Shen, et al. (25 co-authors) 2000 Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39 Nucleic Acids Res 28:1397-1406
Rocha E. P. C., A. Danchin, A. Viari, 1999a Universal replication bias in bacteria Mol. Microbiol 32:11-16[Web of Science][Medline]
———. 1999b Functional and evolutionary roles of long repeats in prokaryotes Res. Microbiol 150:725-733[Medline]
Rocha E. P. C., P. Guerdoux-Jamet, I. Moszer, A. Viari, A. Danchin, 2000 Implication of gene distribution in the bacterial chromosome for the bacterial cell factory J. Biotechnol 78:209-219[Web of Science][Medline]
Romero H., A. Zavala, H. Musto, 2000 Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces Nucleic Acids Res 28:2084-2090
Salzberg S. L., A. J. Salzberg, A. R. Kerlavage, J.-F. Tomb, 1998 Skewed oligomers and origins of replication Gene 217:57-67[Web of Science][Medline]
Shigenobu S., H. Watanabe, M. Hattori, Y. Sakaki, H. Ishikawa, 2000 Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS Nature 407:81-86[Medline]
Shirai M., H. Hirakawa, M. Kimoto, et al. (77 co-authors) 2000 Comparison of whole genome sequences of Chlamydia pneumoniae J138 from Japan and CWL029 from USA Nucleic Acids Res 28:2311-2314
Sniegowski P. D., P. J. Gerrish, T. Johnson, A. Shaver, 2000 The evolution of mutation rates: separating causes from consequences Bioessays 22:1057-1066[Web of Science][Medline]
Stephens R. S., S. Kalman, C. Lammel, et al. (12 co-authors) 1998 Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis Science 282:754-759
Sueoka N., 1962 On the genetic basis of variation and heterogeneity of DNA base composition Proc. Natl. Acad. Sci. USA 48:582-591
———. 1999 Two aspects of DNA base composition: G+C content and translation-coupled deviation from intra-strand rule of A=T and G=C J. Mol. Evol 49:49-62[Web of Science][Medline]
Taddei F., I. Matic, B. Godelle, M. Radman, 1997 To be a mutator, or how pathogenic and commensal bacteria can evolve rapidly Trends Microbiol 5:427-429[Web of Science][Medline]
Takami H., K. Nakasone, Y. Takaki, et al. (12 co-authors) 2000 Complete genome sequence of the alkaliphilic bacterium Bacillus halodurans and genomic sequence comparison with Bacillus subtilis Nucleic Acids Res 28:4317-4331
Tao H., C. Bausch, C. Richmond, F. R. Blattner, T. Conway, 1999 Functional genomics: expression analysis of Escherichia coli growing on minimal and rich media J. Bacteriol 181:6425-6440
Thompson J. D., D. G. Higgins, T. J. Gibson, 1994 CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice Nucleic Acids Res 22:4673-4680
Tillier E. R., R. A. Collins, 2000a The contributions of replication orientation, gene direction, and signal sequences to base-composition asymmetries in bacterial genomes J. Mol. Evol 50:249-257[Web of Science][Medline]
———. 2000b Replication orientation affects the rate and direction of bacterial gene evolution J. Mol. Evol 51:459-463[Web of Science][Medline]
Zar J. H., 1996 Biostatistical analysis. 3rd edition Prentice Hall, Upper Saddle River, N.J
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  W. Hao and G. B. Golding Does Gene Translocation Accelerate the Evolution of Laterally Transferred Genes? Genetics, August 1, 2009; 182(4): 1365 - 1375. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 N. V. Sernova and M. S. Gelfand Identification of replication origins in prokaryotic genomes Brief Bioinform, September 1, 2008; 9(5): 376 - 391. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 X. J. Min and D. A. Hickey DNA Asymmetric Strand Bias Affects the Amino Acid Composition of Mitochondrial Proteins DNA Res, November 1, 2007; (2007) dsm019v1. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. P.C. Rocha, M. Touchon, and E. J. Feil Similar compositional biases are caused by very different mutational effects Genome Res., December 1, 2006; 16(12): 1537 - 1547. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 K. R. Sakharkar, M. K. Sakharkar, C. Verma, and V. T. K. Chow Comparative study of overlapping genes in bacteria, with special reference to Rickettsia prowazekii and Rickettsia conorii Int J Syst Evol Microbiol, May 1, 2005; 55(3): 1205 - 1209. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. P. C. Rocha The replication-related organization of bacterial genomes Microbiology, June 1, 2004; 150(6): 1609 - 1627. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 A. Fuglsang Evolution of Prokaryotic DNA: Intragenic and Extragenic Divergences Observed with Orthologs from Three Related Species Mol. Biol. Evol., June 1, 2004; 21(6): 1152 - 1159. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. P. C. Rocha and A. Danchin Gene essentiality determines chromosome organisation in bacteria Nucleic Acids Res., November 15, 2003; 31(22): 6570 - 6577. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. T. Herbeck, D. P. Wall, and J. J. Wernegreen Gene expression level influences amino acid usage, but not codon usage, in the tsetse fly endosymbiont Wigglesworthia Microbiology, September 1, 2003; 149(9): 2585 - 2596. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
G. Achaz, E. Coissac, P. Netter, and E. P. C. Rocha Associations Between Inverted Repeats and the Structural Evolution of Bacterial Genomes Genetics, August 1, 2003; 164(4): 1279 - 1289. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Y. Zivanovic, P. Lopez, H. Philippe, and P. Forterre Pyrococcus genome comparison evidences chromosome shuffling-driven evolution Nucleic Acids Res., May 1, 2002; 30(9): 1902 - 1910. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||