Streamlining and simplification of microbial genome architecture

Intergenic transposable elements are not randomly distributed in bacteria

Insertion sequences (ISs) are mobile genetic elements in bacterial genomes. In general, intergenic IS elements are probably less deleterious for their hosts than intragenic ISs, simply because they have a lower likelihood of disrupting native genes. However, since promoters, Shine-Dalgarno sequences, and transcription factor binding sites are intergenic and upstream of genes, I hypothesized that not all neighboring gene orientations (NGOs) are selectively equivalent for IS insertion. To test this, I analyzed the NGOs of all intergenic ISs in 326 fully sequenced bacterial chromosomes. Of the 116 genomes with enough IS elements for statistical analysis, 68 have significantly more ISs between convergently oriented genes than expected, and 46 have significantly fewer ISs between divergently oriented genes. This suggests that natural selection molds intergenic IS distributions because they are least intrusive between convergent gene pairs and most intrusive between divergent gene pairs.

Constraints and plasticity in genome and molecular-phenome evolution

Multiple constraints variously affect different parts of the genomes of diverse life forms. The selective pressures that shape the evolution of viral, archaeal, bacterial and eukaryotic genomes differ markedly, even among relatively closely related animal and bacterial lineages; by contrast, constraints affecting protein evolution seem to be more universal. The constraints that shape the evolution of genomes and phenomes are complemented by the plasticity and robustness of genome architecture, expression and regulation. Taken together, these findings are starting to reveal complex networks of evolutionary processes that must be integrated to attain a new synthesis of evolutionary biology.

Extensive loss of RNA editing sites in rapidly evolving Silene mitochondrial genomes: selection vs. retroprocessing as the driving force

Theoretical arguments suggest that mutation rates influence the proliferation and maintenance of RNA editing. We identified RNA editing sites in five species within the angiosperm genus Silene that exhibit highly divergent mitochondrial mutation rates. We found that mutational acceleration has been associated with rapid loss of mitochondrial editing sites. In contrast, we did not find a significant difference in the frequency of editing in chloroplast genes, which lack the mutation rate variation observed in the mitochondrial genome. As found in other angiosperms, the rate of substitution at RNA editing sites in Silene greatly exceeds the rate at synonymous sites, a pattern that has previously been interpreted as evidence for selection against RNA editing. Alternatively, we suggest that editing sites may experience higher rates of C-to-T mutation than other portions of the genome. Such a pattern could be caused by gene conversion with reverse-transcribed mRNA (i.e., retroprocessing). If so, the genomic distribution of RNA editing site losses in Silene suggests that such conversions must be occurring at a local scale such that only one or two editing sites are affected at a time. Because preferential substitution at editing sites appears to occur in angiosperms regardless of the mutation rate, we conclude that mitochondrial rate accelerations within Silene have "fast-forwarded" a preexisting pattern but have not fundamentally changed the evolutionary forces acting on RNA editing sites.

Low nucleotide diversity for the expanded organelle and nuclear genomes of Volvox carteri supports the mutational-hazard hypothesis

The noncoding-DNA content of organelle and nuclear genomes can vary immensely. Both adaptive and nonadaptive explanations for this variation have been proposed. This study addresses a nonadaptive explanation called the mutational-hazard hypothesis and applies it to the mitochondrial, plastid, and nuclear genomes of the multicellular green alga Volvox carteri. Given the expanded architecture of the V. carteri organelle and nuclear genomes (60-85% noncoding DNA), the mutational-hazard hypothesis would predict them to have less silent-site nucleotide diversity (π(silent)) than their more compact counterparts from other eukaryotes-ultimately reflecting differences in 2N(g)μ (twice the effective number of genes per locus in the population times the mutation rate). The data presented here support this prediction: Analyses of mitochondrial, plastid, and nuclear DNAs from seven V. carteri forma nagariensis geographical isolates reveal low values of π(silent) (0.00038, 0.00065, and 0.00528, respectively), much lower values than those previously observed for the more compact organelle and nuclear DNAs of Chlamydomonas reinhardtii (a close relative of V. carteri). We conclude that the large noncoding-DNA content of the V. carteri genomes is best explained by the mutational-hazard hypothesis and speculate that the shift from unicellular to multicellular life in the ancestor that gave rise to V. carteri contributed to a low V. carteri population size and thus a reduced 2N(g)μ. Complete mitochondrial and plastid genome maps for V. carteri are also presented and compared with those of C. reinhardtii.

Conditions for the evolution of gene clusters in bacterial genomes

Genes encoding proteins in a common pathway are often found near each other along bacterial chromosomes. Several explanations have been proposed to account for the evolution of these structures. For instance, natural selection may directly favour gene clusters through a variety of mechanisms, such as increased efficiency of coregulation. An alternative and controversial hypothesis is the selfish operon model, which asserts that clustered arrangements of genes are more easily transferred to other species, thus improving the prospects for survival of the cluster. According to another hypothesis (the persistence model), genes that are in close proximity are less likely to be disrupted by deletions. Here we develop computational models to study the conditions under which gene clusters can evolve and persist. First, we examine the selfish operon model by re-implementing the simulation and running it under a wide range of conditions. Second, we introduce and study a Moran process in which there is natural selection for gene clustering and rearrangement occurs by genome inversion events. Finally, we develop and study a model that includes selection and inversion, which tracks the occurrence and fixation of rearrangements. Surprisingly, gene clusters fail to evolve under a wide range of conditions. Factors that promote the evolution of gene clusters include a low number of genes in the pathway, a high population size, and in the case of the selfish operon model, a high horizontal transfer rate. The computational analysis here has shown that the evolution of gene clusters can occur under both direct and indirect selection as long as certain conditions hold. Under these conditions the selfish operon model is still viable as an explanation for the evolution of gene clusters.

Genes involved in a common pathway or function are frequently found near each other on bacterial chromosomes. A number of hypotheses have been previously presented to explain this observation. A particularly influential theory is the selfish operon model, which posits that horizontal transfer could promote gene clustering by favouring transfer of arrangements of genes that are close together. Subsequent theoretical development and analysis of genomic data have contributed to the debate about the plausibility of this model. Here, by re-examining the evolutionary dynamics of gene clusters, we provide and discuss conditions under which gene clusters can evolve. We find that first, some form of bias for clustering is required for clusters to evolve. This bias can be in the form of bias in horizontal transfer towards genes that are close together, or direct natural selection for gene proximity. Our computational work does not present a theoretical obstacle to the selfish operon model as a possible explanation for the evolution of gene clusters.

Localized plasticity in the streamlined genomes of vinyl chloride respiring Dehalococcoides

Vinyl chloride (VC) is a human carcinogen and widespread priority pollutant. Here we report the first, to our knowledge, complete genome sequences of microorganisms able to respire VC, Dehalococcoides sp. strains VS and BAV1. Notably, the respective VC reductase encoding genes, vcrAB and bvcAB, were found embedded in distinct genomic islands (GEIs) with different predicted integration sites, suggesting that these genes were acquired horizontally and independently by distinct mechanisms. A comparative analysis that included two previously sequenced Dehalococcoides genomes revealed a contextually conserved core that is interrupted by two high plasticity regions (HPRs) near the Ori. These HPRs contain the majority of GEIs and strain-specific genes identified in the four Dehalococcoides genomes, an elevated number of repeated elements including insertion sequences (IS), as well as 91 of 96 rdhAB, genes that putatively encode terminal reductases in organohalide respiration. Only three core rdhA orthologous groups were identified, and only one of these groups is supported by synteny. The low number of core rdhAB, contrasted with the high rdhAB numbers per genome (up to 36 in strain VS), as well as their colocalization with GEIs and other signatures for horizontal transfer, suggests that niche adaptation via organohalide respiration is a fundamental ecological strategy in Dehalococccoides. This adaptation has been exacted through multiple mechanisms of recombination that are mainly confined within HPRs of an otherwise remarkably stable, syntenic, streamlined genome among the smallest of any free-living microorganism.

On the need for widespread horizontal gene transfers under genome size constraint

Background

While eukaryotes primarily evolve by duplication-divergence expansion (and reduction) of their own gene repertoire with only rare horizontal gene transfers, prokaryotes appear to evolve under both gene duplications and widespread horizontal gene transfers over long evolutionary time scales. But, the evolutionary origin of this striking difference in the importance of horizontal gene transfers remains by and large a mystery.

Hypothesis

We propose that the abundance of horizontal gene transfers in free-living prokaryotes is a simple but necessary consequence of two opposite effects: i) their apparent genome size constraint compared to typical eukaryote genomes and ii) their underlying genome expansion dynamics through gene duplication-divergence evolution, as demonstrated by the presence of many tandem and block repeated genes. In principle, this combination of genome size constraint and underlying duplication expansion should lead to a coalescent-like process with extensive turnover of functional genes. This would, however, imply the unlikely, systematic reinvention of functions from discarded genes within independent phylogenetic lineages. Instead, we propose that the long-term evolutionary adaptation of free-living prokaryotes must have resulted in the emergence of efficient non-phylogenetic pathways to circumvent gene loss.

Implications

This need for widespread horizontal gene transfers due to genome size constraint implies, in particular, that prokaryotes must remain under strong selection pressure in order to maintain the long-term evolutionary adaptation of their "mutualized" gene pool, beyond the inevitable turnover of individual prokaryote species. By contrast, the absence of genome size constraint for typical eukaryotes has presumably relaxed their need for widespread horizontal gene transfers and strong selection pressure. Yet, the resulting loss of genetic functions, due to weak selection pressure and inefficient gene recovery mechanisms, must have ultimately favored the emergence of more complex life styles and ecological integration of many eukaryotes.

Reviewers

This article was reviewed by Pierre Pontarotti, Eugene V Koonin and Sergei Maslov.

Streamlining and large ancestral genomes in Archaea inferred with a phylogenetic birth-and-death model

Homologous genes originate from a common ancestor through vertical inheritance, duplication, or horizontal gene transfer. Entire homolog families spawned by a single ancestral gene can be identified across multiple genomes based on protein sequence similarity. The sequences, however, do not always reveal conclusively the history of large families. To study the evolution of complete gene repertoires, we propose here a mathematical framework that does not rely on resolved gene family histories. We show that so-called phylogenetic profiles, formed by family sizes across multiple genomes, are sufficient to infer principal evolutionary trends. The main novelty in our approach is an efficient algorithm to compute the likelihood of a phylogenetic profile in a model of birth-and-death processes acting on a phylogeny.

We examine known gene families in 28 archaeal genomes using a probabilistic model that involves lineage- and family-specific components of gene acquisition, duplication, and loss. The model enables us to consider all possible histories when inferring statistics about archaeal evolution. According to our reconstruction, most lineages are characterized by a net loss of gene families. Major increases in gene repertoire have occurred only a few times. Our reconstruction underlines the importance of persistent streamlining processes in shaping genome composition in Archaea. It also suggests that early archaeal genomes were as complex as typical modern ones, and even show signs, in the case of the methanogenic ancestor, of an extremely large gene repertoire.

Deletional bias across the three domains of life

Elevated levels of genetic drift are hypothesized to be a dominant factor that influences genome size evolution across all life-forms. However, increased levels of drift appear to be correlated with genome expansion in eukaryotes but with genome contraction in bacteria, suggesting that these two groups of organisms experience vastly different mutational inputs and selective constraints. To determine the contribution of small insertion and deletion events to the differences in genome organization between eukaryotes and prokaryotes, we systematically surveyed 17 taxonomic groups across the three domains of life. Based on over 5,000 indel events in noncoding regions, we found that deletional events outnumbered insertions in all groups examined. The extent of deletional bias, when measured by the total length of insertions to deletions, revealed a marked disparity between eukaryotes and prokaryotes, whereas the ratio was close to one in the three eukaryotic groups examined, deletions outweighed insertions by at least a factor of 10 in most prokaryotes. Moreover, the strength of deletional bias is associated with the proportion of coding regions in prokaryotic genomes. Considering that genetic drift is a stochastic process and does not discriminate the exact nature of mutations, the degree of bias toward deletions provides an explanation to the differential responses of eukaryotes and prokaryotes to elevated levels of drift. Furthermore, deletional bias, rather than natural selection, is the primary mechanism by which the compact gene packing within most prokaryotic genomes is maintained.

The consequences of genetic drift for bacterial genome complexity

Genetic drift, which is particularly effective within small populations, can shape the size and complexity of genomes by affecting the fixation of deleterious mutations. In Bacteria, assessing the contribution of genetic drift to genome evolution is problematic because the usual methods, based on intraspecific polymorphisms, can be thwarted by difficulties in delineating species' boundaries. The increased availability of sequenced bacterial genomes allows application of an alternative estimator of drift, the genome-wide ratio of replacement to silent substitutions in protein-coding sequences. This ratio, which reflects the action of purifying selection across the entire genome, shows a strong inverse relationship with genome size, indicating that drift promotes genome reduction in bacteria.

Deletion rate evolution and its effect on genome size and coding density

Deletion rates are thought to be important factors in determining the genome size of organisms in nature. Although it is indisputable that deletions, and thus deletion rates, affect genome size, it is unclear how, or indeed if, genome size is regulated via the deletion rate. Here, we employ a mathematical model to determine the evolutionary fate of deletion rate mutants. Simulations are employed to explore the interactions between deletions, deletion rate mutants, and genome size. The results show that, in this model, the fate of deletion rate mutants will depend on the fraction of essential genomic material, on the frequency of sexual recombination, as well as on the population size of the organism. We find that there is no optimal deletion rate in any state. However, at one critical coding density, all changes in deletion rate are neutral and the rate may drift either up or down. As a consequence, the coding density of the genome is expected to fluctuate around this critical density. Characteristic differences in the impact of deletion rate mutations on prokaryote and eukaryote genomes are described.

Phenotypic and genotypic characterization of multiple strains of the diazotrophic cyanobacterium, Crocosphaera watsonii, isolated from the open ocean

Diazotrophic cyanobacteria have long been recognized as important sources of reduced nitrogen (N) and therefore are important ecosystem components. Until recently, species of the filamentous cyanobacterium Trichodesmium were thought to be the primary sources of fixed N to the open ocean euphotic zone. It is now recognized that unicellular cyanobacteria are also important contributors, with members of the oligotrophic genus Crocosphaera being the only cultured examples. Herein we genetically and phenotypically characterize 10 strains isolated from the tropical Atlantic and North Pacific Oceans, and show that although all of the strains are highly similar at the genetic level, with the internal transcribed sequence (ITS) region sequence varying by approximately 2 bp on average, there are many unexpected phenotypic differences between the isolates (e.g. cell size, temperature optima and range, extracellular material excretion and variability in rates of nitrogen fixation). However based on the observed sequence similarity, we propose that all of these isolates are members of the genus Crocosphaera (type strain Crocosphaera watsonii WH8501), and that the phenotypic diversity we see may reflect ecologically important variation relevant for modelling N(2) fixation in the oligotrophic ocean.

Ecological specialization in a spatially structured population of the thermophilic cyanobacterium Mastigocladus laminosus

Laboratory evolution experiments suggest the potential for microbial populations to contribute significant ecological variation to ecosystems, yet the functional importance of genetic diversity within natural populations of microorganisms is largely unknown. Here, we investigated the distribution of genetic and phenotypic variation for a population of the cyanobacterium Mastigocladus laminosus distributed along the temperature gradient of White Creek, Yellowstone NP. A total of 153 laboratory strains were directly isolated from five sites with mean annual temperatures ranging between 39 and 54 degrees C. Genetic characterization at four nitrogen metabolism genes identified 15 closely related lineages in the population sample. These lineages were distributed nonrandomly along White Creek, but the observed geographic structure could not be explained by limited dispersal capabilities. Temperature performance experiments with six M. laminosus lineages that maximized their respective relative abundances at different positions along the gradient provided evidence for niche differentiation within the population. Niche differentiation included a tradeoff in performance at high and low temperatures, respectively. The physiological variation of these lineages in laboratory culture was generally well matched to the prevailing temperature conditions experienced by these organisms in situ. These results suggest that sympatric diversification along an ecological selection gradient can be a potent source of evolutionary innovation in microbial populations.

Genetic diversity in the Paramecium aurelia species complex

Current understanding of the population genetics of free-living unicellular eukaryotes is limited, and the amount of genetic variability in these organisms is still a matter of debate. We characterized-reproductively and genetically-worldwide samples of multiple Paramecium species belonging to a cryptic species complex, Paramecium aurelia, whose species have been shown to be reproductively isolated. We found that levels of genetic diversity both in the nucleus and in the mitochondrion are substantial within groups of reproductively compatible P. aurelia strains but drop considerably when strains are partitioned according to their phylogenetic groupings. Our study reveals the existence of discrepancies between the mating behavior of a number of P. aurelia strains and their multilocus genetic profile, a controversial finding that has major consequences for both the current methods of species assignment and the species problem in the P. aurelia complex.

Trends in prokaryotic evolution revealed by comparison of closely related bacterial and archaeal genomes

In order to explore microevolutionary trends in bacteria and archaea, we constructed a data set of 41 alignable tight genome clusters (ATGCs). We show that the ratio of the medians of nonsynonymous to synonymous substitution rates (dN/dS) that is used as a measure of the purifying selection pressure on protein sequences is a stable characteristic of the ATGCs. In agreement with previous findings, parasitic bacteria, notwithstanding the sometimes dramatic genome shrinkage caused by gene loss, are typically subjected to relatively weak purifying selection, presumably owing to relatively small effective population sizes and frequent bottlenecks. However, no evidence of genome streamlining caused by strong selective pressure was found in any of the ATGCs. On the contrary, a significant positive correlation between the genome size, as well as gene size, and selective pressure was observed, although a variety of free-living prokaryotes with very close selective pressures span nearly the entire range of genome sizes. In addition, we examined the connections between the sequence evolution rate and other genomic features. Although gene order changes much faster than protein sequences during the evolution of prokaryotes, a strong positive correlation was observed between the "rearrangement distance" and the amino acid distance, suggesting that at least some of the events leading to genome rearrangement are subjected to the same type of selective constraints as the evolution of amino acid sequences.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

Population size and genome size in fishes: a closer look

The several thousand-fold range in genome size among animals has remained a subject of active research and debate for more than half a century, but no satisfactory explanation has yet been provided. Many one-dimensional models have been postulated, but so far none has been successful in accounting for observed patterns in genome size diversity. The recent model based on differences in effective population size appeared to gain empirical support with a study of genome size and inferred effective population size in fishes, but there were several questionable aspects of the analysis. First, it was based on an assumption that microsatellite heterozygosity indicates long-term effective population size, whereas in actuality these markers evolve quickly and are sensitive to demographic events. Second, it included both ancient polyploids and non-polyploids, the former of which did not gain their current genome sizes through the accumulation of slightly deleterious mutations as required in the model. Third, the analysis neglected the tremendous influence that Pleistocene glaciation bottlenecks had on heterozygosities in freshwater (and far less so, marine) fishes. In sum, it is apparent that genomes reached their current sizes in most fishes long before contemporary microsatellite heterozygosities were shaped, and that ancient polyploidy rather than the accumulation of mildly deleterious transposon insertions in small populations is the dominant factor that has influenced the large end of the range of genome sizes among fishes.

Nucleotide diversity in the mitochondrial and nuclear compartments of Chlamydomonas reinhardtii: investigating the origins of genome architecture

Background

The magnitude of intronic and intergenic DNA can vary substantially both within and among evolutionary lineages; however, the forces responsible for this disparity in genome compactness are conjectural. One explanation, termed the mutational-burden hypothesis, posits that genome compactness is primarily driven by two nonadaptive processes: mutation and random genetic drift – the effects of which can be discerned by measuring the nucleotide diversity at silent sites (π_silent), defined as noncoding sites and the synonymous sites of protein-coding regions. The mutational-burden hypothesis holds that π_silent is negatively correlated to genome compactness. We used the model organism Chlamydomonas reinhardtii, which has a streamlined, coding-dense mitochondrial genome and an noncompact, intron-rich nuclear genome, to investigate the mutational-burden hypothesis. For measuring π_silent we sequenced the complete mitochondrial genome and portions of 7 nuclear genes from 7 geographical isolates of C. reinhardtii.

Results

We found significantly more nucleotide diversity in the nuclear compartment of C. reinhardtii than in the mitochondrial compartment: net values of π_silent for the nuclear and mitochondrial genomes were 32 × 10^-3 and 8.5 × 10^-3, respectively; and when insertions and deletions (indels) are factored in, these values become 49 × 10^-3 for the nuclear DNA and 11 × 10^-3 for the mitochondrial DNA (mtDNA). Furthermore, our investigations of C. reinhardtii revealed 4 previously undiscovered mitochondrial introns, one of which contains a fragment of the large-subunit (LSU) rRNA gene and another of which is found in a region of the LSU-rRNA gene not previously reported (for any taxon) to contain introns.

Conclusion

At first glance our results are in opposition to the mutational-burden hypothesis: π_silent was approximately 4 times greater in the nuclear compartment of C. reinhardtii relative to the mitochondrial compartment. However, when we consider the encumbrance of noncoding DNA in each of these C. reinhardtii compartments, we conclude that introns in the mtDNA impose a greater burden than those in the nuclear DNA and suggest that the same may be true for the intergenic regions. Overall, we cannot reject the mutational-burden hypothesis and feel that more data on nucleotide diversity from green algae and other protists are needed.

The impact of genome analyses on our understanding of ammonia-oxidizing bacteria

The availability of whole-genome sequences for ammonia-oxidizing bacteria (AOB) has led to dramatic increases in our understanding of these environmentally important microorganisms. Their genomes are smaller than many other members of the proteobacteria and may indicate genome reductions consistent with their limited lifestyle. The genomes have a surprising level of gene repetition including genes for ammonia catabolism, iron acquisition, and insertion sequences. The gene profiles reveal limited genes for catabolism and transport of complex organic compounds, but complete pathways for some other compounds. This led to the observation of chemolithoheterotrophic growth of Nitrosomonas europaea. Genes for sucrose synthesis/degradation were identified. The core metabolic module of aerobic ammonia oxidation, the extraction of electrons from hydroxylamine to generate proton-motive force and reductant, has evolutionary roots in the denitrification inventory of anaerobic sulfur-dependent bacteria. The extension by ammonia monooxygenase provides a mechanism to feed this module using ammonia and O(2).

The rise and fall of deleterious mutation

It is well established that selection is less efficient in small populations than in large ones. Here we review the impact of this effect by considering the gradual selective purging of deleterious mutation over time. We outline an approach to explore the dynamics of this process, and highlight its profound implications.

Bacterial group II introns: not just splicing

Group II introns are both catalytic RNAs (ribozymes) and mobile retroelements that were discovered almost 14 years ago. It has been suggested that eukaryotic mRNA introns might have originated from the group II introns present in the alphaproteobacterial progenitor of the mitochondria. Bacterial group II introns are of considerable interest not only because of their evolutionary significance, but also because they could potentially be used as tools for genetic manipulation in biotechnology and for gene therapy. This review summarizes what is known about the splicing mechanisms and mobility of bacterial group II introns, and describes the recent development of group II intron-based gene-targetting methods. Bacterial group II intron diversity, evolutionary relationships, and behaviour in bacteria are also discussed.

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.