
==== Front
Genome BiolGenome Biology1465-69061465-6914BioMed Central London gb-2005-6-2-r141569394310.1186/gb-2005-6-2-r14ResearchAccelerated evolution associated with genome reduction in a free-living prokaryote Dufresne Alexis 1dufresne@sb-roscoff.frGarczarek Laurence 1garczare@sb-roscoff.frPartensky Frdric 1partensky@sb-roscoff.fr1 Station Biologique, UMR 7127 CNRS et Universit Paris 6, BP74, 29682 Roscoff Cedex, France2005 14 1 2005 6 2 R14 R14 5 10 2004 2 12 2004 7 12 2004 Copyright  2005 Dufresne et al.; licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Prochlorococcus sp. are marine bacteria with very small genomes. The mechanisms by which these reduced genomes have evolved appears, however, to be distinct from those that have led to small genome size in intracellular bacteria.

Background
Three complete genomes of Prochlorococcus species, the smallest and most abundant photosynthetic organism in the ocean, have recently been published. Comparative genome analyses reveal that genome shrinkage has occurred within this genus, associated with a sharp reduction in G+C content. As all examples of genome reduction characterized so far have been restricted to endosymbionts or pathogens, with a host-dependent lifestyle, the observed genome reduction in Prochlorococcus is the first documented example of such a process in a free-living organism.

Results
Our results clearly indicate that genome reduction has been accompanied by an increased rate of protein evolution in P. marinus SS120 that is even more pronounced in P. marinus MED4. This acceleration has affected every functional category of protein-coding genes. In contrast, the 16S rRNA gene seems to have evolved clock-like in this genus. We observed that MED4 and SS120 have lost several DNA-repair genes, the absence of which could be related to the mutational bias and the acceleration of amino-acid substitution.

Conclusions
We have examined the evolutionary mechanisms involved in this process, which are different from those known from host-dependent organisms. Indeed, most substitutions that have occurred in Prochlorococcus have to be selectively neutral, as the large size of populations imposes low genetic drift and strong purifying selection. We assume that the major driving force behind genome reduction within the Prochlorococcus radiation has been a selective process favoring the adaptation of this organism to its environment. A scenario is proposed for genome evolution in this genus.
==== Body
Background
The size of bacterial genomes is primarily the result of two counteracting processes: the acquisition of new genes by gene duplication or by horizontal gene transfer; and the deletion of non-essential genes. Genomic flux created by these gains and losses of genetic information can substantially alter gene content. This process drives divergence of bacterial species and eventually adaptation to new ecological niches [1]. In some cases, gene deletion may prevail over gene acquisition, leading to genome reduction. This process has occurred several times during evolution and has been well documented for cellular organelles [2,3], obligate pathogens such as Mycoplasma genitalium [4] or phytoplasmas [5] and symbionts such as the insect endosymbiont Buchnera [6-8] or the hyperthermophile Nanoarchaeum equitans [9]. In the case of organelles, the degree of genome reduction can be extensive as a result of massive gene transfer into the host nucleus, allowing maintenance of the corresponding functions in the resulting composite organism. Mitochondrial or chloroplast genomes, for instance, can be as small as 6 kilobases (kb) [10] and 35 kb [11], respectively. In the case of obligate host-dependent bacteria, the reduction is more limited because the relationships with their hosts are less intimate than for organelles in eukaryotic cells. Thus, obligatory pathogens need to retain a minimum of functions that allow them to infect new hosts and to avoid host defenses, and obligate endosymbionts carry genes which are absolutely necessary for host survival. For instance, a substantial part (approximately 10 %) of the Buchnera genome is devoted to biosynthesis of amino acids which are essential to its host [6].

So far, all characterized examples of genome reduction have been associated with a change from a free-living to a host-dependent lifestyle [12]. It is therefore intriguing that a similar phenomenon of genome reduction has occurred within the free-living marine cyanobacterial genus Prochlorococcus [13-15]. The latter is present at high abundance (often over 105 cells/ml) in all nutrient-poor areas of the world's oceans between 40N and 40S and is probably the most abundant photosynthetic organism on Earth [16,17]. It has been shown that two major ecotypes exist within this genus [18]. The first is adapted to grow at the base of the illuminated layer and displays a high divinyl-chlorophyll b to a ratio; the second inhabits the upper layer of the ocean and has a low divinyl-chlorophyll b to a ratio [19]. The genome of one high-light-adapted (HL) strain, Prochlorococcus marinus MED4 [14], and of two low-light-adapted (LL) strains, P. marinus SS120 [13] and Prochlorococcus species MIT9313 [14], have recently been sequenced and annotated.

Phylogenetic trees based on 16S rRNA sequences [18] or 16S-23S ribosomal internal transcribed spacer sequences [20] show that Prochlorococcus sp. MIT9313 branches at the base of the Prochlorococcus radiation, close to the Synechococcus group [21]. In contrast, the Prochlorococcus HL clade, encompassing the MED4 strain, appears to be the most recently evolved Prochlorococcus group, consistent with the fact that this clade is much less diversified than are the LL clades.

Despite the close relatedness of these strains, their genomes vary widely in terms of size, G+C content and the number of protein-coding genes (Table 1). While the general characteristics of the MIT9313 genome are very similar to those of the Synechococcus sp. WH8102 genome [22], MED4 has the smallest genome for a photosynthetic organism known to date and the SS120 genome is only 90 kb larger. Furthermore, this genome reduction is clearly accompanied by a drift in G+C content, a phenomenon that commonly occurs during the evolution of host-dependent genomes [23]. However, the evolutionary mechanisms involved in the genome reductive process are most probably different from those that have occurred in host-dependent organisms. Using comparative sequence analyses of the four genomes of marine picocyanobacteria published to date, we have attempted to better understand the causes and consequences of this phenomenon and to address the relationships between genome reduction and niche adaptation in marine picocyanobacteria.

Results
Synteny and genome stability
Alignments of whole genomes show a strong conservation of the gene order between MED4 and SS120 (Figure 1a). There are only five inversions larger than 20 kb between these two genomes. In contrast, the large number of inversions and translocations and the shorter size of the colinear segments between SS120 and MIT9313 on the one hand and MIT9313 and WH8102 on the other hand (Figure 1b,c) indicate that extensive genome rearrangements have occurred not only between Synechococcus and Prochlorococcus but also between MIT9313 and the two other Prochlorococcus strains (see also Figure 2 in [14]). The degree of synteny observed between the four marine picocyanobacteria genomes strengthens the hypothesis of a more recent divergence of the clades containing MED4 and SS120 than of the clade containing MIT9313.

Overall genome composition
The downsizing of MED4 and SS120 genomes during evolution is associated with a genome-wide adenine (A) and thymine (T) enrichment (Table 1). The bias is most pronounced at neutral sites such as intergenic regions (MED4, 76.6% A+T; SS120, 69.3% A+T) and third-codon positions of protein-coding genes (MED4, 79.7% A+T; SS120, 73.85% A+T). This bias has little effect on ribosomal RNA genes (5S, 16S and 23S) which have a G+C content greater than 50% in all four picocyanobacterial genomes. In both MED4 and SS120, the single rRNA gene cluster can easily be spotted as a G+C-rich anomaly compared to the rest of the genome (see for example, Figure 1 in [15]). In direct contrast, protein-coding genes are strongly affected by the extreme base composition of these genomes. First, the bias influences codon usage since, for a given amino acid, AT-rich codons are preferentially used (Figure 3a). Second, the amino-acid composition of the proteins themselves is affected (Figure 3b). Indeed, when compared to Prochlorococcus sp. MIT9313 and Synechococcus sp. WH8102, the genes of P. marinus MED4 and SS120 contain fewer amino acids encoded by G+C-rich codons (for example, alanine or arginine) and more amino acids encoded by A+T-rich codons (for example, isoleucine or lysine).

Orthologous gene pool size
A total of 1,306 orthologs belonging to all major functional categories are common to the four genomes (see Additional data file 1) and probably constitute an estimate of the core of genes conserved in all marine picocyanobacteria. This is sensibly more than the pool of around 1,000 orthologs identified by W.R. Hess [15]. The difference certainly results from the use by the latter author of a low E-value threshold (10e-12) for BLAST comparisons. In contrast, our analysis is based on identification of reciprocal best hits without the use of any particular threshold (apart from the default BLAST threshold) and consequently allows the detection of orthologous relationships whatever the gene lengths or the level of similarity. Still, our ortholog identification process is rather strict and the set of orthologs identified in this study probably corresponds to a lower estimate of the actual number of orthologs shared by the four genomes. This set of genes represents a substantial percentage of the total pool of all protein-coding genes in P. marinus MED4 (73.2%) and SS120 (69.2%) and about half of the gene set in Prochlorococcus sp. MIT9313 (56.2%) and Synechococcus sp. WH8102 (51.1%). These percentages are consistent with the differences in the respective number of genes within these genomes (Table 1) and are compatible with the assumption that a massive gene loss has occurred in MED4 and SS120 during their evolution from a Prochlorococcus ancestor with a larger genome [13-15].

Accelerated rate of evolution of protein-coding genes in Prochlorococcus
Because biased base composition seems to constrain amino-acid usage in the Prochlorococcus genomes, we have investigated whether it also affects the rate of protein sequence evolution in these genomes. We used the 1,306 orthologs common to the four genomes to estimate the amino-acid substitution rate in each genome. Branch lengths calculated for a given tree topology (the same topology as for the 16S rRNA gene tree; see Figure 2) are 0.46, 0.22, 0.16 and 0.14 amino acid substitutions per site for branches dA, dB, dC and dX, respectively. Using Synechococcus sp. WH8102 as the outgroup, we tested the rate-constancy hypothesis and computed the ratios of branch lengths. Relative rate tests (two-cluster and branch length tests) indicate that protein sequences evolved at significantly different rates (P < 0.001) between MED4, SS120 and MIT9313. Therefore the hypothesis of a constant evolutionary rate between these strains can be rejected for protein-coding genes. The calculation of branch-length ratios reveals that the amino-acid substitution rate is 2.04-fold higher inMED4 than in SS120 (dA/dB) and 3.81-fold higher in MED4 than in MIT9313 ((dA+dX)/dC). This rate is also 2.31-fold higher for SS120 than for MIT9313 ((dB+dX)/dC). Computation of branch lengths for each functional category shows that the increased rate of amino-acid replacement in protein sequences concerns every category (Figure 4 and Table 2). These results imply that the rate of amino-acid substitution increased during evolution of the Prochlorococcus genus concomitantly with genome reduction and increase in A+T content.

Synonymous and nonsynonymous substitutions
The ratio of the rate of nonsynonymous substitutions (dN) to the rate of synonymous substitutions (dS) is commonly used to measure the relative rate of purifying selection acting at the protein level. We determined dS and dN for each gene pair of every group of orthologs and their values were averaged for each genome. Surprisingly, we observed saturation at synonymous sites for all genome pairs (dS > 2) and the calculation of the dN/dS ratio was thus impossible. Still, the average dN was higher between MED4 and SS120 (0.36) than between SS120 and MIT9313 (0.32). The lowest dN was observed between MIT9313 and WH8102 (0.24), a finding which is consistent with the relative acceleration of amino-acid substitutions in MED4 and in SS120.

DNA-repair systems
A shift in base composition may reflect the loss of DNA-repair genes and we therefore determined the presence or absence of genes involved in these mechanisms. As the mutational pressure is toward a high A+T content in both MED4 and SS120, we looked more closely at those genes whose absence could increase the frequency of G:C to A:T mutations. Among the genes putatively encoding DNA-repair enzymes identified in MIT9313 and WH8102, a few are missing in SS120 and/or MED4 (Table 3). Both MED4 and SS120 lack the ada gene, which encodes 6-O-methylguanine-DNA methyltransferase, which repairs alkylated forms of guanine and thymine in DNA. Such alkylations generate lesions that can lead to G:C to A:T transversions [24]. Interestingly, the MED4 genome is the only one among the four picocyanobacteria not to encode the A/G-specific DNA glycosylase MutY, as previously noted by Rocap and co-workers [14]. This enzyme acts with MutT (NTP pyrophosphohydrolase) and MutM (formamido-pyrimidine-DNA glycosylase) in the GO system to avoid misincorporation of oxidized guanine (8-oxoG) in DNA and to repair the base mismatches A:8-oxoG [25]. In Escherichia coli, knocking out both mutM and mutY translates into a 1,000-fold increase of G:C to A:T transversions in comparison to the wild-type strain [26]. In addition to MutT and MutY, MIT9313 and WH8102 encode a third enzyme of the NUDIX hydrolase family that is missing in MED4 and SS120. This hydrolase could act to prevent mutations. However because of the broad substrate specificity of this family, one cannot know with certainty the function of this protein. Likewise, two genes coding for enzymes of the RecF pathway have been lost either by both MED4 and SS120 (DNA helicase RecQ) or only by MED4 (exonuclease RecJ).

Discussion
The process of genome reduction which has occurred within the Prochlorococcus radiation has to our knowledge never been observed so far in any other free-living prokaryote. Since Prochlorococcus sp. MIT9313 has a genome size very similar to that of Synechococcus sp. WH8102 (2.4 megabase-pair (Mbp)), as well as several other marine Synechococcus spp. (M. Ostrowski and D. Scanlan, personal communication), it is reasonable to assume that the common ancestor of all Prochlorococcus species also had a genome size around 2.4 Mbp. Under this hypothesis, the genome reduction which has occurred in MED4 would correspond to around 31%. By comparison, the extent of genome reduction in the insect endosymbiont Buchnera, as compared to a reconstructed ancestral genome, is around 77% [27]. The genome of P. marinus SS120 - and a fortiori the MED4 genome - is considered to be near minimal for a free-living oxyphototrophic organism [13]. It would seem that genome reduction in these organisms probably cannot proceed below a certain limit, corresponding to a gene pool containing all the essential genes of biosynthetic pathways and housekeeping functions (probably including most of the 1,306 four-way orthologous genes identified in this study) plus a number of other genes, including genus-specific as well as niche-specific genes. For instance, MED4 encodes a number of photolyase-related proteins, a few specific ABC transporters (for cyanate, for example; [14] and data not shown). These specific compounds might be critical for survival in the upper water layer, which receives high photon fluxes, UV light and is nutrient-depleted, but less so for life deeper in the water column.

If both Prochlorococcus lineages and host-dependent organisms have undergone genome reduction associated with accelerated substitution rates, these phenomena must have arisen from very different causes as the resulting gene repertoires of the two types of organisms differ tremendously. Indeed, the genome evolution of endosymbionts and obligatory pathogens is driven by two main processes which have mutually reinforcing effects on genome size and evolutionary rates. Being confined inside their host, these bacteria have tiny population sizes and are regularly bottlenecked at each host generation or at each new host infection. Consequently, they experience a strong genetic drift [28] involving an increase in substitution rate. This acceleration results in the accumulation at random of slightly deleterious mutations in protein-coding genes [8,29] as well as in rRNA genes [29,30]. This genetic drift enhances the downsizing of the genome through inactivation and then elimination of potentially beneficial but dispensable genes. Among these, there have been a number of DNA-repair genes, the disappearance of which could have further increased the mutation rate [6,31-33]. Furthermore, a number of genes may be subject to a relaxation of purifying selection which is therefore rendered less effective in maintaining gene function. This relaxation particularly affects genes which have become useless because they are redundant in their host genome, such as genes involved in the biosynthesis of amino acids, nucleotides, fatty acids and even ATP [4-6,8,9,32]. Selection pressure is also reduced for genes involved in environmental sensing and regulatory systems, such as two-component systems, because of the much buffered environment offered by the host [6].

In the free-living genus Prochlorococcus, the very large size of field populations [34] means that these populations are subject to much lower genetic drift and their genomes are subject to much stronger purifying selection than are those of endosymbionts and pathogens [35]. Consequently, the observed accelerated rate of evolution probably results merely from the increase in the mutation rate, which in turn is probably due to the loss of DNA-repair genes, even if one should note that, in P. marinus SS120 only two such genes are missing (Table 3). We observed a similar acceleration of amino-acid substitutions for all functional categories (Figure 4). This finding is more consistent with a global increase in the mutation rate than with relaxed selection, the latter being unlikely to occur to the same extent at all loci. We also assume that most amino-acid substitutions that have occurred in Prochlorococcus proteins are neutral; that is, they have not altered protein function. Indeed, populations of the HL clade which, like MED4, have the most derived protein sequences of all Prochlorococcus species, appear to be the most abundant photosynthetic organisms in the upper layer of the temperate and inter-tropical oceans [16]. Such an ecological success would hardly be possible for organisms handicapped by a large number of slightly deleterious mutations, especially given the fact that most genes are single copy, and so compensation of gene function is generally not possible. The effect of the maintenance of a high level of purifying selection on counteracting deleterious substitutions is particularly obvious in the rRNA genes. Contrary to the protein-coding genes, relative rate tests did not show any significant differences in the rates of evolution of the 16S rRNA genes in the four marine picocyanobacterial genomes, and thus there is no evidence that either SS120 or MED4 could have accumulated mutations destabilizing the secondary structure of their 16S rRNA molecule. One noteworthy consequence of the acceleration in the rates of evolution of protein-coding genes in Prochlorococcus is that phylogenetic reconstructions based on protein sequences are biased. Indeed, this leads to much longer branches for these two strains than for MIT9313. The resulting tree topology most often does not support that obtained with the 16S rRNA gene, for which the molecular clock hypothesis holds true according to our analyses. Thus, rRNA genes are likely to be among the few genes that will give reliable estimates of the phylogenetic distances between Prochlorococcus strains.

If it is neither the relaxation of purifying selection nor an increase in genetic drift that has been the main factor causing Prochlorococcus genome reduction, an alternative possibility is that the latter could be the result of a selective process favoring the adaptation of Prochlorococcus to its environment. The apparently better ecological success in oligotrophic areas of Prochlorococcus species compared to their close relative Synechococcus [16,34], strongly suggests that the reduction of Prochlorococcus genome size could provide a competitive advantage to the former. Indeed, extensive comparisons of the gene complements of these two organisms show very few examples - at least among genes for which function is known - of the occurrence of specific genes in MED4 which could explain its better adaptation (data not shown). One noteworthy exception is the presence in Prochlorococcus, but not Synechococcus, of flavodoxin and ferritin, two proteins that possibly give Prochlorococcus a better resistance to iron stress. Apart from that, Synechococcus appears more like a generalist, in particular with regard to nitrogen or phosphorus uptake and assimilation [22], and should a priori be more suited to sustain competition. Hence, we assume that the key to the success of Prochlorococcus resides less in the development of a specific complex or pathway to cope better with unfavorable conditions than in the simplification of its genome and cell organization, which can allow this organism to make substantial economies in energy and material for cell maintenance.

The mere reduction in genome size per se is a potential source of substantial economies for the cell, as it reduces the amount of nitrogen and phosphorus, two particularly limiting elements in the upper part of the ocean, which are necessary, for instance, in DNA synthesis. Another advantage is that it allows a concomitant reduction in cell volume. It has been previously suggested (see, for example [36]) that, for a phytoplanktonic organism, a small cell volume confers two selective advantages by reducing self-shading (the package effect) and by increasing the cell surface-to-volume ratio, which can improve nutrient uptake. The first advantage would improve the fitness of the LL strains, whereas the second would offer an advantage to the HL strains living in nutrient-depleted surface waters. Finally, cell division is less costly for a small than for a large cell. On the basis of these observations, we assume that the major driving force for genome reduction within the Prochlorococcus radiation has been the selection for a more economical lifestyle. The bias toward an A+T-rich genome in MED4 and SS120 is also consistent with this hypothesis, as it can be seen as a way to economize on nitrogen. Indeed, an AT base-pair contains seven atoms of nitrogen, one less than a GC base-pair.

With this hypothesis in mind, we propose a possible scenario for the evolution of Prochlorococcus genomes. Using a rate of 16S rRNA divergence of 1% per 50 million years [37], one can estimate that the differentiation of these two genera is as recent as 150 million years, as the molecular clock hypothesis holds for this gene in Prochlorococcus and Synechococcus. The ancestral Prochlorococcus cells must have developed in the LL niche, a niche probably left free by other picocyanobacteria. Given the considerable difference in genome size between the LL strains MIT9313 and SS120, it appears that genome reduction itself must have started in one (or possibly several) lineage(s) within the LL niche some time after Prochlorococcus differentiation from its common ancestor with marine Synechococcus species. Why the selection has affected only one (or some?) and not all Prochlorococcus lineages remains unclear. Examination of the gene repertoire of P. marinus SS120 [13] suggests that this genome reduction must have concerned the random loss of dispensable genes from many different pathways. At some point during evolution, some genes involved in DNA repair have been affected; these would include the ada gene, which may be responsible for the shift in base composition, but also possibly several others, not necessarily involved in GC to AT mutation repair (see Table 3). Loss of these genes may have led to an increase in the mutation rate and therefore in the rate of evolution of protein-coding genes, accompanied by a more rapid genome shrinkage and a shift of base composition toward AT. It is worth noting that one likely consequence of this genome-wide compositional shift is the absence of the adaptive codon bias in the genomes of Prochlorococcus species MED4 and SS120. AT-rich codons are preferentially used whatever the amino acid (Figure 3a). Thus, codon usage in these genomes appears to reflect more the local base-composition bias than the selection for a more efficient translation through the use of optimal codons. The same conclusion has been drawn for other small genomes with high A+T content [28,38].

Later during evolution (around 80 million years ago, according to the degree of 16S rRNA sequence divergence between MED4 and SS120) one LL population which probably already had a significantly reduced cell and genome size must have progressively adapted to the HL niche and eventually recolonized the upper layer. How this change in ecological niche was possible is still hard to define. Comparison of the gene set that differs between the LL-adapted SS120 and the HL-adapted MED4 shows that very few genes might be sufficient to shift from one to the other niche, including a multiplication of hli genes [39] and the differential retention of genes which were present in the common ancestor of Prochlorococcus and Synechococcus, (such as the photolyases and cyanate transporters mentioned above) and were secondarily lost in the LL-adapted lineages.

Conclusions
Genome evolution in the free-living genus Prochlorococcus has similar features to that in host-dependent prokaryotes: genome reduction, bias toward a low G+C content, acceleration in the evolution rate of protein-coding genes, and loss of DNA-repair genes. In contrast to the latter organisms, however, in Prochlorococcus this evolution does not appear to be the result of genetic drift or relaxed selection being exerted on some gene categories. Indeed, purifying selection is very efficient in Prochlorococcus, as rRNA genes have evolved at a similar rate in all genomes. Despite the decrease in G+C content and an accelerated rate of evolution of protein-coding genes, purifying selection must also act on these genes and avoid potentially deleterious mutations. We hypothesize that a reduction in genome size (which allows a concomitant reduction in cell size and substantial economies in energy and nutrients) can constitute a selective advantage for life in the open ocean, both at depths where photon energy is low and in surface waters where nutrients are scarce.

Genome shrinkage in Prochlorococcus has led to populations highly specialized to narrow ecological niches, at the expense of versatility and competitiveness in changing conditions. Indeed, not only is the distribution of the Prochlorococcus genus limited to low latitudes (40N and 40S, see [34]) but the different ecotypes are themselves more or less confined to a restricted part of the euphotic layer [40]; for example, they experience only limited changes in temperature and salinity. Paradoxically, because warm oligotrophic areas constitute a very large part of the world's oceans, the ecological niches (both LL and HL) occupied by Prochlorococcus species are huge, and thus this organism appears globally, despite its specialization, as one of the most successful oxyphototrophs on Earth.

Materials and methods
Genome sequence data
The complete genome sequences and annotations of Prochlorococcus marinus MED4, P. marinus SS120, Prochlorococcus sp. MIT9313 and Synechococcus sp. WH8102 (accession numbers: NC_005071, NC_005072, NC_005042 and NC_005070 respectively) were downloaded from the Genome division of the NCBI Entrez system. A few additional genes which were modeled in at least one genome and were present in the other genomes but not modeled (because of their small size, for example) were included in our dataset (see Additional data file 2).

Alignment of whole genomes
Genome sequences translated in their six reading frames were aligned with the Promer program of the MUMmer 3.0 system [41].

Codon and amino-acid usage
Codon usage was computed for every open reading frame (ORF) of each genome with the EMBOSS program cusp. Amino-acid usage was derived from the results produced by cusp.

Identification of orthologous proteins
We used a sequence-similarity based approach which is similar to the procedure used for the cluster of orthologous groups (COGs [42]). For each genome pair, all-against-all BLAST [43] comparisons were performed using protein sequences and reciprocal genome-specific best hits were identified. We considered genes as being probable orthologs when they were included in groups of size four in which each gene was the best hit of the three others. From similarity searches against the COG database, orthologs were assigned to functional categories according to those defined for the COG system. Because of the lack of a particular category for photosynthesis genes, the latter were assigned to the 'energy production and conversion' COG category. Other genes which fell into more than one of the 19 COG categories have been assigned to a supplementary category called 'miscellaneous'.

Phylogenetic branch length estimations
Protein sequences from each of the groups of four orthologous genes were aligned using ClustalW [44] with default parameters. After exclusion of all gap sites, individual alignments were concatenated in one super-alignment of 388,120 sites. Gamma distances [45] with an alpha parameter of 1 were estimated between each pair of sequences of the super-alignment. Phylogenetic branch lengths were calculated from distances with the ordinary least-squares method [45]. Relative rate tests (two-cluster test and Branch length test) were applied in order to test the constancy of amino-acid substitution rates between the three Prochlorococcus genomes (hypothesis of the molecular clock). The same analysis was applied to orthologs of each functional category.

Estimate of synonymous and nonsynonymous substitution rates
Nucleotide sequences of each group of orthologs were aligned with Protal2dna according to alignments of their corresponding amino-acid sequences [46]. Pairwise estimates of the synonymous (dS) and non-synonymous (dN) substitution rates were obtained from the Yn00 program of the PAML 3.13 package [47].

Additional data files
The following additional data are available with the online version of this article. Additional data file 1 lists the orthologous genes classified by functional category. Orthologous genes were assigned to the functional categories of COG system. Photosynthesis genes were assigned to the 'energy production and conversion' COG category. Genes falling in more than one of the 19 COG categories have been assigned to a supplementary category called 'miscellaneous'. Additional data file 2 is a fasta file of orthologous genes which were modeled in at least one genome and present but not modeled in the other genomes.

Supplementary Material
Additional data file 1
The orthologous genes classified by functional category

Click here for additional data file

 Additional data file 2
A a fasta file of orthologous genes which were modeled in at least one genome and present but not modeled in the other genomes

Click here for additional data file

 Acknowledgements
We are very grateful to Martin Ostrowski and Dave Scanlan for their critical reading of the manuscript. This work was supported by the European Union Program MARGENES (QLRT-2001-01226), the EU FP6 Network of Excellence 'Marine Genomics Europe' and by the French programs Genomer (Rgion Bretagne) and Ouest-Genopole. AD is supported by a doctoral fellowship from Rgion Bretagne.

Figures and Tables
Figure 1 Alignments of complete genome sequences of marine picocyanobacteria. Genome sequences are translated in their six reading frames. (a) Comparison of the MED4 and SS120 genomes; (b) comparison of the SS120 and MIT9313 genomes; (c) comparison of the MIT9313 and WH8102 genomes. Colinear segments are shown in red and inversions in green. Translocated segments are above or below the diagonal.

Figure 2 Phylogenetic tree of 16S rRNA genes from the four marine picocyanobacteria. Neighbor-joining tree with Kimura 2-parameter correction. The bootstrap value (1,000 replications) is shown in boldface. Lengths of the branches dA, dB, dC and dX (see text) are given below the branches. N1, node 1, branchpoint between MED4 and SS120; N2, node 2, branchpoint between MIT9313 and Node 1.

Figure 3 Influence of mutational bias in codon usage and amino-acid usage. (a) Percentage use of AT-rich codons in the four marine picocyanobacteria. Amino acids are ranked according to AT content of their respective codons. Methionine and tryptophan, which are both encoded by only one codon, have been discarded from the analysis. (b) Percentage use of amino acids in marine picocyanobacteria.

Figure 4 Amino-acid substitution rate per functional category. Branch lengths computed for each functional category between (a) MED4 and SS120, (b) SS120 and MIT9313 and (c) MED4 and MIT9313. In the three comparisons, branch-length values are aligned along a line with a slope much greater than 1, indicating that acceleration of the substitution rates occurs in every functional category. Axes represent the number of amino-acid substitutions per site. Red circle, amino-acid transport and metabolism; green circle, carbohydrate transport and metabolism; yellow circle, cell-cycle control; blue triangle, cell wall/membrane biogenesis; pink circle, coenzyme metabolism; red square, defense mechanisms; green square, energy production and conversion; yellow square, function unknown; blue square, general function prediction only; pink square, inorganic ion transport and metabolism; red triangle, intracellular trafficking; green triangle, lipid transport and metabolism; yellow triangle, nucleotide transport and metabolism; blue circle, posttranslational modification, protein turnover; pink triangle, replication, recombination and repair; red diamond, secondary metabolite biosynthesis, transport and catabolism; green diamond, signal transduction mechanisms; yellow diamond, transcription; blue diamond, translation; black circle, miscellaneous.

Table 1 General features of the genomes of the four marine picocyanobacteria used in this study

Genome	Size (Mbp)	GC%	Number of protein-coding genes	
P. marinus MED4	1.66	30.8	1,716	
P. marinus SS120	1.75	36.4	1,882	
Prochlorococcus sp. MIT9313	2.41	50.7	2,273	
Synechococcus sp. WH8102	2.43	59.4	2,525	
Table 2 Number and percentage of orthologous genes per functional category

Category	Number of genes	% of orthologs	
Amino-acid transport and metabolism	94	7.2	
Carbohydrate transport and metabolism	50	3.8	
Cell-cycle control	17	1.3	
Cell wall/membrane biogenesis	55	4.2	
Coenzyme metabolism	99	7.6	
Defense mechanisms	14	1.1	
Energy production and conversion	106	8.1	
Function unknown	269	20.6	
General function prediction only	116	8.9	
Inorganic ion transport and metabolism	47	3.6	
Intracellular trafficking	13	1.0	
Lipid transport and metabolism	25	1.9	
Nucleotide transport and metabolism	39	3.0	
Posttranslational modification, protein turnover	59	4.5	
Replication, recombination and repair	51	3.9	
Secondary metabolite biosynthesis, transport and catabolism	6	0.5	
Signal transduction mechanisms	11	0.8	
Transcription	26	2.0	
Translation	127	9.7	
Miscellaneous	82	6.3	
Table 3 DNA-repair genes missing only in P. marinus MED4 or in both MED4 and SS120

Gene	COG	Product	MED4	SS120	MIT9313	WH8102	

                              ada/ogt
                           	0350	6-O-methylguanine-DNA methyltransferase	-	-	PMT0269	SYNW1680	

                              mutY
                           	1194	A/G-specific DNA glycosylase	-	Pro1789	PMT0135	SYNW0115	
recQ	0514	Superfamily II DNA helicase	-	-	PMT0189	SYNW1958	
recJ	0608	Single-stranded DNA-specific exonuclease	-	Pro0984	PMT0761	SYNW1206	
exoI/xseA	1570	Exonuclease VII large subunit	-	Pro0111	PMT1641	SYNW2181	
xseB	1722	Exonuclease VII small subunit	-	Pro0112	PMT1642	SYNW2182	
-	0494	NUDIX hydrolase family	-	-	PMT1026	SYNW1334	
Genes in bold are involved in repair of G:C to A:T mutations.
==== Refs
Lawrence JG Roth JR  Charlebois RL  Genomic flux: genome evolution by gene loss and acquisition. Organization of the Prokaryotic Genome 1999 Washington, DC: American Society for Microbiology 263 289 
Andersson SG Kurland CG  Reductive evolution of resident genomes. Trends Microbiol 1998 6 263 268 9717214 10.1016/S0966-842X(98)01312-2 
Martin W  Gene transfer from organelles to the nucleus: frequent and in big chunks. Proc Natl Acad Sci USA 2003 100 8612 8614 12861078 10.1073/pnas.1633606100 
Fraser CM Gocayne JD White O Adams MD Clayton RA Fleischmann RD Bult CJ Kerlavage AR Sutton G Kelley JM  The minimal gene complement of Mycoplasma genitalium. Science 1995 270 397 403 7569993 
Oshima K Kakizawa S Nishigawa H Jung HY Wei W Suzuki S Arashida R Nakata D Miyata S Ugaki M Namba S  Reductive evolution suggested from the complete genome sequence of a plant-pathogenic phytoplasma. Nat Genet 2004 36 27 29 14661021 10.1038/ng1277 
Shigenobu S Watanabe H Hattori M Sakaki Y Ishikawa H  Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature 2000 407 81 86 10993077 10.1038/35024074 
Tamas I Klasson L Canback B Naslund AK Eriksson AS Wernegreen JJ Sandstrom JP Moran NA Andersson SG  50 million years of genomic stasis in endosymbiotic bacteria. Science 2002 296 2376 2379 12089438 10.1126/science.1071278 
van Ham RC Kamerbeek J Palacios C Rausell C Abascal F Bastolla U Fernandez JM Jimenez L Postigo M Silva FJ  Reductive genome evolution in Buchnera aphidicola. Proc Natl Acad Sci USA 2003 100 581 586 12522265 10.1073/pnas.0235981100 
Waters E Hohn MJ Ahel I Graham DE Adams MD Barnstead M Beeson KY Bibbs L Bolanos R Keller M  The genome of Nanoarchaeum equitans: insights into early archaeal evolution and derived parasitism. Proc Natl Acad Sci USA 2003 100 12984 12988 14566062 10.1073/pnas.1735403100 
Conway DJ Fanello C Lloyd JM Al-Joubori BM Baloch AH Somanath SD Roper C Oduola AM Mulder B Povoa MM  Origin of Plasmodium falciparum malaria is traced by mitochondrial DNA. Mol Biochem Parasitol 2000 111 163 171 11087926 10.1016/S0166-6851(00)00313-3 
Kohler S Delwiche CF Denny PW Tilney LG Webster P Wilson RJ Palmer JD Roos DS  A plastid of probable green algal origin in Apicomplexan parasites. Science 1997 275 1485 1489 9045615 10.1126/science.275.5305.1485 
Moran NA  Microbial minimalism: genome reduction in bacterial pathogens. Cell 2002 108 583 586 11893328 10.1016/S0092-8674(02)00665-7 
Dufresne A Salanoubat M Partensky F Artiguenave F Axmann IM Barbe V Duprat S Galperin MY Koonin EV Le Gall F  Genome sequence of the cyanobacterium Prochlorococcus marinus SS120, a nearly minimal oxyphototrophic genome. Proc Natl Acad Sci USA 2003 100 10020 10025 12917486 10.1073/pnas.1733211100 
Rocap G Larimer FW Lamerdin J Malfatti S Chain P Ahlgren NA Arellano A Coleman M Hauser L Hess WR  Genome divergence in two Prochlorococcus ecotypes reflects oceanic niche differentiation. Nature 2003 424 1042 1047 12917642 10.1038/nature01947 
Hess WR  Genome analysis of marine photosynthetic microbes and their global role. Curr Opin Biotechnol 2004 15 191 198 15193326 10.1016/j.copbio.2004.03.007 
Partensky F Blanchot J Vaulot D  Charpy L, Larkum AWD  Differential distribution and ecology of Prochlorococcus and Synechococcus in oceanic waters: a review. Marine Cyanobacteria 1999 Monaco: Muse Ocanographique 457 475 
Garcia-Pichel F Belnap J Neuer S Schanz F  Estimates of cyanobacterial biomass and its distribution. Archiv Hydrobiol 2003 109 Suppl 148 213 228 
Moore LR Rocap G Chisholm SW  Physiology and molecular phylogeny of coexisting Prochlorococcus ecotypes. Nature 1998 393 464 467 9624000 10.1038/30965 
Moore LR Chisholm SW  Photophysiology of the marine cyanobacterium Prochlorococcus: ecotypic differences among cultured isolates. Limnol Oceanogr 1999 44 628 638 
Rocap G Distel DL Waterbury JB Chisholm SW  Resolution of Prochlorococcus and Synechococcus ecotypes by using 16S-23S ribosomal DNA internal transcribed spacer sequences. Appl Environ Microbiol 2002 68 1180 1191 11872466 10.1128/AEM.68.3.1180-1191.2002 
Fuller NJ Marie D Partensky F Vaulot D Post AF Scanlan DJ  Clade-specific 16S ribosomal DNA oligonucleotides reveal the predominance of a single marine Synechococcus clade throughout a stratified water column in the Red Sea. Appl Environ Microbiol 2003 69 2430 2443 12732508 10.1128/AEM.69.5.2430-2443.2003 
Palenik B Brahamsha B Larimer FW Land M Hauser L Chain P Lamerdin J Regala W Allen EE McCarren J  The genome of a motile marine Synechococcus. Nature 2003 424 1037 1042 12917641 10.1038/nature01943 
Moran NA  Tracing the evolution of gene loss in obligate bacterial symbionts. Curr Opin Microbiol 2003 6 512 518 14572545 10.1016/j.mib.2003.08.001 
Mackay WJ Han S Samson LD  DNA alkylation repair limits spontaneous base substitution mutations in Escherichia coli. J Bacteriol 1994 176 3224 3230 8195077 
Michaels ML Cruz C Grollman AP Miller JH  Evidence that MutY and MutM combine to prevent mutations by an oxidatively damaged form of guanine in DNA. Proc Natl Acad Sci USA 1992 89 7022 7025 1495996 
Horst JP Wu TH Marinus MG  Escherichia coli mutator genes. Trends Microbiol 1999 7 29 36 10068995 10.1016/S0966-842X(98)01424-3 
Moran NA Mira A  The process of genome shrinkage in the obligate symbiont Buchnera aphidicola. Genome Biol 2001 2 research0054.1 0054.12 11790257 10.1186/gb-2001-2-12-research0054 
Wernegreen JJ Moran NA  Evidence for genetic drift in endosymbionts (Buchnera): analyses of protein-coding genes. Mol Biol Evol 1999 16 83 97 10331254 
Moran NA  Accelerated evolution and Muller's rachet in endosymbiotic bacteria. Proc Natl Acad Sci USA 1996 93 2873 2878 8610134 10.1073/pnas.93.7.2873 
Lambert JD Moran NA  Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria. Proc Natl Acad Sci USA 1998 95 4458 4462 9539759 10.1073/pnas.95.8.4458 
Koonin EV Mushegian AR Rudd KE  Sequencing and analysis of bacterial genomes. Curr Biol 1996 6 404 416 8723345 10.1016/S0960-9822(02)00508-0 
Glass JI Lefkowitz EJ Glass JS Heiner CR Chen EY Cassell GH  The complete sequence of the mucosal pathogen Ureaplasma urealyticum. Nature 2000 407 757 762 11048724 10.1038/35037619 
Akman L Yamashita A Watanabe H Oshima K Shiba T Hattori M Aksoy S  Genome sequence of the endocellular obligate symbiont of tsetse flies, Wigglesworthia glossinidia. Nat Genet 2002 32 402 407 12219091 10.1038/ng986 
Partensky F Hess WR Vaulot D  Prochlorococcus, a marine photosynthetic prokaryote of global significance. Microbiol Mol Biol Rev 1999 63 106 127 10066832 
Ohta T  The nearly neutral theory of molecular evolution. Annu Rev Ecol Syst 1992 23 263 286 10.1146/annurev.es.23.110192.001403 
Chisholm SW  Falkowski PG, Woodhead AD  Phytoplankton size. Primary Productivity and Biogeochemical Cycles in the Sea 1992 New York: Plenum Press 213 237 
Ochman H Wilson AC  Evolution in bacteria: evidence for a universal substitution rate in cellular genomes. J Mol Evol 1987 26 74 86 3125340 
Andersson SG Sharp PM  Codon usage and base composition in Rickettsia prowazekii. J Mol Evol 1996 42 525 536 8662004 
Bhaya D Dufresne A Vaulot D Grossman A  Analysis of the hli gene family in marine and freshwater cyanobacteria. FEMS Microbiol Lett 2002 215 209 219 12399037 10.1016/S0378-1097(02)00934-5 
West NJ Scanlan DJ  Niche-partitioning of Prochlorococcus populations in a stratified water column in the eastern North Atlantic Ocean. Appl Environ Microbiol 1999 65 2585 2591 10347047 
Kurtz S Phillippy A Delcher AL Smoot M Shumway M Antonescu C Salzberg SL  Versatile and open software for comparing large genomes. Genome Biol 2004 5 R12 14759262 10.1186/gb-2004-5-2-r12 
Tatusov RL Koonin EV Lipman DJ  A genomic perspective on protein families. Science 1997 278 631 637 9381173 10.1126/science.278.5338.631 
Altschul SF Madden TL Schaffer AA Zhang J Zhang Z Miller W Lipman DJ  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997 25 3389 3402 9254694 10.1093/nar/25.17.3389 
Thompson JD Higgins DG Gibson TJ  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994 22 4673 4680 7984417 
Nei M Kumar S  Molecular Evolution and Phylogenetic 2000 Oxford: Oxford University Press 
protal2dna 
Yang Z  PAML: A program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 1997 13 555 556 9367129

