An Approach to Population and Evolutionary Genetic Theory for Genes in Mitochondria and Chloroplasts, and Some Results

We developed population genetic theory for organelle genes, using an infinite alleles model appropriate for molecular genetic data, and considering the effects of mutation and random drift on the frequencies of selectively neutral alleles. The effects of maternal inheritance and vegetative segregation of organelle genes are dealt with by defining new effective gene numbers, and substituting these for 2N(e) in classical theory of nuclear genes for diploid organisms. We define three different effective gene numbers. The most general is N(lambda), defined as a function of population size, number of organelle genomes per cell, and proportions of genes contributed by male and female gametes to the zygote. In many organisms, vegetative segregation of organelle genomes and intracellular random drift of organelle gene frequencies combine to produce a predominance of homoplasmic cells within individuals in the population. Then, the effective number of organelle genes is N(eo), a simple function of the numbers of males and females and of the maternal and paternal contributions to the zygote. Finally, when the paternal contribution is very small, N( eo) is closely approximated by the number of females, N( f). Then if the sex ratio is 1, the mean time to fixation or loss of new mutations is approximately two times longer for nuclear genes than for organelle genes, and gene diversity is approximately four times greater. The difference between nuclear and organelle genes disappears or is reversed in animals in which males have large harems. The differences between nuclear and organelle gene behavior caused by maternal inheritance and vegetative segregation are generally small and may be overshadowed by differences in mutation rates to neutral alleles. For monoecious organisms, the effective number of organelle genes is approximately equal to the total population size N. We also show that a population can be effectively subdivided for organelle genes at migration rates which result in panmixis for nuclear genes, especially if males migrate more than females.     

Mutation accumulation in transfer RNAs: molecular evidence for Muller's ratchet in mitochondrial genomes

The accumulation of deleterious mutations is thought to be a major factor preventing the long-term persistence of obligately asexual lineages relative to their sexual ancestors. This phenomenon is also of potential relevance to sexual species that harbor asexually propagating organelle genomes. A comparative study of the transfer RNA genes in animal mitochondrial and nuclear genomes demonstrates that the former accumulate nucleotide substitutions much more rapidly than do the latter, and several lines of evidence are consistent with the idea that the excess substitutions are mildly deleterious. First, the average binding stability between complementary strands in the stems of mitochondrial tRNAs is less than half that in nuclear tRNAs. Second, most loop sizes in the mitochondrial tRNAs have experienced a net reduction in size over evolutionary time, and they are nearly 50 times more variable in the mitochondrial than in the nuclear genome. Third, although nearly 20% of the nucleotides in nuclear tRNA genes (particularly those involved in tertiary interactions) are invariant across all animal taxa and all tRNA species, there are no invariant sites in the mitochondrial tRNAs. These observations, as well as results from recent laboratory experiments, are consistent with the hypothesis that nonrecombining organelle genomes are subject to gradual loss of fitness due to the cumulative chance fixation of mildly deleterious mutations.      

The mutational hazard hypothesis of organelle genome evolution: 10 years on

Why is there such a large variation in size and noncoding DNA content among organelle genomes? One explanation is that this genomic variation results from differences in the rates of organelle mutation and random genetic drift, as opposed to being the direct product of natural selection. Along these lines, the mutational hazard hypothesis (MHH) holds that ‘excess’ DNA is a mutational liability (because it increases the potential for harmful mutations) and, thus, has a greater tendency to accumulate in an organelle system with a low mutation rate as opposed to one with a high rate of mutation. Various studies have explored this hypothesis and, more generally, the relationship between organelle genome architecture and the mode and efficiency of organelle DNA repair. Although some of these investigations are in agreement with the MHH, others have contradicted it; nevertheless, they support a central role of mutation, DNA maintenance pathways and random genetic drift in fashioning organelle chromosomes. Arguably, one of the most important contributions of the MHH is that it has sparked crucial, widespread discussions about the importance of nonadaptive processes in genome evolution.     

Mutational decay and age of chloroplast and mitochondrial genomes transferred recently to angiosperm nuclear chromosomes

Transfers of organelle DNA to the nucleus established several thousand functional genes in eukaryotic chromosomes over evolutionary time. Recent transfers have also contributed nonfunctional plastid (pt)- and mitochondrion (mt)-derived DNA (termed nupts and numts, respectively) to plant nuclear genomes. The two largest transferred organelle genome copies are 131-kb nuptDNA in rice (Oryza sativa) and 262-kb numtDNA in Arabidopsis (Arabidopsis thaliana). These transferred copies were compared in detail with their bona fide organelle counterparts, to which they are 99.77% and 99.91% identical, respectively. No evidence for purifying selection was found in either nuclear integrant, indicating that they are nonfunctional. Mutations attributable to 5-methylcytosine hypermutation have occurred at a 6- to 10-fold higher rate than other point mutations in Arabidopsis numtDNA and rice nuptDNA, respectively, revealing this as a major mechanism of mutational decay for these transferred organelle sequences. Short indels occurred preferentially within homopolymeric stretches but were less frequent than point mutations. The 131-kb nuptDNA is absent in the O. sativa subsp. indica or Oryza rufipogon nuclear genome, suggesting that it was transferred within the O. sativa subsp. japonica lineage and, as revealed by sequence comparisons, after its divergence from the indica chloroplast lineage. The time of the transfer for the rice nupt was estimated as 148,000 (74,000--296,000) years ago and that for the Arabidopsis numtDNA as 88,000 (44,000--176,000) years ago. The results reveal transfer and integration of entire organelle genomes into the nucleus as an ongoing evolutionary process and uncover mutational mechanisms affecting organelle genomes recently transferred into a new mutational environment.     

Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates

Genome size and complexity vary tremendously among eukaryotic species and their organelles. Comparisons across deeply divergent eukaryotic lineages have suggested that variation in mutation rates may explain this diversity, with increased mutational burdens favoring reduced genome size and complexity. The discovery that mitochondrial mutation rates can differ by orders of magnitude among closely related angiosperm species presents a unique opportunity to test this hypothesis. We sequenced the mitochondrial genomes from two species in the angiosperm genus Silene with recent and dramatic accelerations in their mitochondrial mutation rates. Contrary to theoretical predictions, these genomes have experienced a massive proliferation of noncoding content. At 6.7 and 11.3 Mb, they are by far the largest known mitochondrial genomes, larger than most bacterial genomes and even some nuclear genomes. In contrast, two slowly evolving Silene mitochondrial genomes are smaller than average for angiosperms. Consequently, this genus captures approximately 98% of known variation in organelle genome size. The expanded genomes reveal several architectural changes, including the evolution of complex multichromosomal structures (with 59 and 128 circular-mapping chromosomes, ranging in size from 44 to 192 kb). They also exhibit a substantial reduction in recombination and gene conversion activity as measured by the relative frequency of alternative genome conformations and the level of sequence divergence between repeat copies. The evolution of mutation rate, genome size, and chromosome structure can therefore be extremely rapid and interrelated in ways not predicted by current evolutionary theories. Our results raise the hypothesis that changes in recombinational processes, including gene conversion, may be a central force driving the evolution of both mutation rate and genome structure.     

Rate accelerations in nuclear 18S rDNA of mycoheterotrophic and parasitic angiosperms

Rate variation in genes from all three genomes has been observed frequently in plant lineages with a parasitic and mycoheterotrophic mode of life. While the loss of photosynthetic ability leads to a relaxation of evolutionary constraints in genes involved in the photosynthetic apparatus, it remains to be determined how prevalent increased substitution rates are in nuclear DNA of non-photosynthetic angiosperms. In this study we infer rates of molecular evolution of 18S rDNA of all parasitic and mycoheterotorphic plant families (except Lauraceae and Polygalaceae) using relative rate tests. In several holoparasitic and mycoheterotrophic plant lineages extremely high substitution rates are observed compared to other photosynthetic angiosperms. The position and frequency of these substitutions have been identified to understand the mutation dynamics of 18S rRNA in achlorophyllous plants. Despite the presence of significantly elevated substitution rates, very few mutations occur in major functional and structural regions of the small ribosomal molecule, providing evidence that the efficiency of the translational apparatus in non-photosynthetic plants has not been affected.     

Reducing the genome size of organelles favours gene transfer to the nucleus

Endosymbiotic organelles exhibit strong genetic erosion during their evolution as a result of the loss of unnecessary genes and of gene transfer to the nucleus. The reasons for this erosion are much debated. Unidirectionality of DNA exchange between cell compartments could favour biased gene transfer, but selection might also act to favour nuclear localization of genes, for example, because organelles accumulate more mutations than do nuclei. Selection for rapid replication might be a general cause of organelle genome reduction. This selection also accounts for the compactness of organelle genomes.     

Organelle Gene Diversity under Migration, Mutation, and Drift: Equilibrium Expectations, Approach to Equilibrium, Effects of Heteroplasmic Cells, and Comparison to Nuclear Genes

We developed stochastic population genetic theory for mitochondrial and chloroplast genes, using an infinite alleles model appropriate for molecular genetic data. We considered the effects of mutation, random drift, and migration in a finite island model on selectively neutral alleles. Recurrence equations were obtained for the expectation of gene diversities within zygotes, within colonies, and between colonies. The variables are number and sizes of colonies, migration rates, sex ratios, degree of paternal transmission, number of germ line cell divisions, effective number of segregating organelle genomes, and mutation rate. Computer solutions of the recurrence equations were used to study the approach to equilibrium. Gene diversities equilibrate slowly, while GST, used to measure population subdivision, equilibrates rapidly. Approximate equilibrium equations for gene diversities and GST can be obtained by substituting Neo and me, simple functions of the numbers of breeding or migrating males and females and of the degree of paternal transmission, for the effective numbers of genes and migration rates in the corresponding equations for nuclear genes. The approximate equations are not valid when the diversity within individuals is large compared to that between individuals, as is often true for the D-loop of animal mtDNA. We used the exact equations to verify that organelle genes often show more subdivision than nuclear genes; however, we also identified the range of breeding and migrating sex ratios for which population subdivision is greater for nuclear genes. Finally, we show that gene diversities are higher for nuclei than for organelles over a larger range of sex ratios in a subdivided population than in a panmictic population.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.     

Origin, evolution and genetic effects of nuclear insertions of organelle DNA

In eukaryotes, nuclear genomes are subject to an influx of DNA from mitochondria and plastids. The nuclear insertion of organellar sequences can occur during the illegitimate repair of double-stranded breaks. After integration, nuclear organelle DNA is modified by point mutations, and by deletions. Insertion of organelle DNA into nuclear genes is not rare and can potentially have harmful effects. In humans, some insertions of nuclear mitochondrial DNA are associated with heritable diseases. It remains to be determined whether nuclear organelle DNA can contribute beneficially to gene evolution.     

Natural selection on functional modules, a genome-wide analysis

Classically, the functional consequences of natural selection over genomes have been analyzed as the compound effects of individual genes. The current paradigm for large-scale analysis of adaptation is based on the observed significant deviations of rates of individual genes from neutral evolutionary expectation. This approach, which assumed independence among genes, has not been able to identify biological functions significantly enriched in positively selected genes in individual species. Alternatively, pooling related species has enhanced the search for signatures of selection. However, grouping signatures does not allow testing for adaptive differences between species. Here we introduce the Gene-Set Selection Analysis (GSSA), a new genome-wide approach to test for evidences of natural selection on functional modules. GSSA is able to detect lineage specific evolutionary rate changes in a notable number of functional modules. For example, in nine mammal and Drosophilae genomes GSSA identifies hundreds of functional modules with significant associations to high and low rates of evolution. Many of the detected functional modules with high evolutionary rates have been previously identified as biological functions under positive selection. Notably, GSSA identifies conserved functional modules with many positively selected genes, which questions whether they are exclusively selected for fitting genomes to environmental changes. Our results agree with previous studies suggesting that adaptation requires positive selection, but not every mutation under positive selection contributes to the adaptive dynamical process of the evolution of species.     

Why mitochondrial genes are most often found in nuclei

A very small fraction of the proteins required for the propagation and function of mitochondria are coded by their genomes, while nuclear genes code the vast majority. We studied the migration of genes between the two genomes when transfer mechanisms mediate this exchange. We could calculate the influence of differential mutation rates, as well as that of biased transfer rates, on the partitioning of genes between the two genomes. We observe no significant difference in partitioning for haploid and diploid cell populations, but the effective size of cell populations is important. For infinitely large effective populations, higher mutation rates in mitochondria than in nuclear genomes are required to drive mitochondrial genes to the nuclear genome. In the more realistic case of finite populations, gene transfer favoring the nucleus and/or higher mutation rates in the mitochondrion will drive mitochondrial genes to the nucleus. We summarize experimental data that identify a gene transfer process mediated by vacuoles that favors the accumulation of mitochondrial genes in the nuclei of modern cells. Finally, we compare the behavior of mitochondrial genes for which transfer to the nucleus is neutral or influenced by purifying selection.     

Construction of a highly error-prone DNA polymerase for developing organelle mutation systems

A novel family of DNA polymerases replicates organelle genomes in a wide distribution of taxa encompassing plants and protozoans. Making error-prone mutator versions of gamma DNA polymerases revolutionised our understanding of animal mitochondrial genomes but similar advances have not been made for the organelle DNA polymerases present in plant mitochondria and chloroplasts. We tested the fidelities of error prone tobacco organelle DNA polymerases using a novel positive selection method involving replication of the phage lambda cI repressor gene. Unlike gamma DNA polymerases, ablation of 3′–5′ exonuclease function resulted in a modest 5–8-fold error rate increase. Combining exonuclease deficiency with a polymerisation domain substitution raised the organelle DNA polymerase error rate by 140-fold relative to the wild type enzyme. This high error rate compares favourably with error-rates of mutator versions of animal gamma DNA polymerases. The error prone organelle DNA polymerase introduced mutations at multiple locations ranging from two to seven sites in half of the mutant cI genes studied. Single base substitutions predominated including frequent A:A (template: dNMP) mispairings. High error rate and semi-dominance to the wild type enzyme in vitro make the error prone organelle DNA polymerase suitable for elevating mutation rates in chloroplasts and mitochondria.     

Reductive genome evolution from the mother of Rickettsia

The Rickettsia genus is a group of obligate intracellular α-proteobacteria representing a paradigm of reductive evolution. Here, we investigate the evolutionary processes that shaped the genomes of the genus. The reconstruction of ancestral genomes indicates that their last common ancestor contained more genes, but already possessed most traits associated with cellular parasitism. The differences in gene repertoires across modern Rickettsia are mainly the result of differential gene losses from the ancestor. We demonstrate using computer simulation that the propensity of loss was variable across genes during this process. We also analyzed the ratio of nonsynonymous to synonymous changes (Ka/Ks) calculated as an average over large sets of genes to assay the strength of selection acting on the genomes of Rickettsia, Anaplasmataceae, and free-living γ-proteobacteria. As a general trend, Ka/Ks were found to decrease with increasing divergence between genomes. The high Ka/Ks for closely related genomes are probably due to a lag in the removal of slightly deleterious nonsynonymous mutations by natural selection. Interestingly, we also observed a decrease of the rate of gene loss with increasing divergence, suggesting a similar lag in the removal of slightly deleterious pseudogene alleles. For larger divergence (Ks > 0.2), Ka/Ks converge toward similar values indicating that the levels of selection are roughly equivalent between intracellular α-proteobacteria and their free-living relatives. This contrasts with the view that obligate endocellular microorganisms tend to evolve faster as a consequence of reduced effectiveness of selection, and suggests a major role of enhanced background mutation rates on the fast protein divergence in the obligate intracellular α-proteobacteria.     

The influence of selection on the evolutionary distance estimated from the base changes observed between homologous nucleotide sequences

In most studies of molecular evolution, the nucleotide base at a site is assumed to change with the apparent rate under functional constraint, and the comparison of base changes between homologous genes is thought to yield the evolutionary distance corresponding to the site-average change rate multiplied by the divergence time. However, this view is not sufficiently successful in estimating the divergence time of species, but mostly results in the construction of tree topology without a time-scale. In the present paper, this problem is investigated theoretically by considering that observed base changes are the results of comparing the survivals through selection of mutated bases. In the case of weak selection, the time course of base changes due to mutation and selection can be obtained analytically, leading to a theoretical equation showing how the selection has influence on the evolutionary distance estimated from the enumeration of base changes. This result provides a new method for estimating the divergence time more accurately from the observed base changes by evaluating both the strength of selection and the mutation rate. The validity of this method is verified by analysing the base changes observed at the third codon positions of amino acid residues with four-fold codon degeneracy in the protein genes of mammalian mitochondria; i.e. the ratios of estimated divergence times are fairly well consistent with a series of fossil records of mammals. Throughout this analysis, it is also suggested that the mutation rates in mitochondrial genomes are almost the same in different lineages of mammals and that the lineage-specific base-change rates indicated previously are due to the selection probably arising from the preference of transfer RNAs to codons.     

Conservative and compensatory evolution in oxidative phosphorylation complexes of angiosperms with highly divergent rates of mitochondrial genome evolution

Interactions between nuclear and mitochondrial gene products are critical for eukaryotic cell function. Nuclear genes encoding mitochondrial‐targeted proteins (N‐mt genes) experience elevated rates of evolution, which has often been interpreted as evidence of nuclear compensation in response to elevated mitochondrial mutation rates. However, N‐mt genes may be under relaxed functional constraints, which could also explain observed increases in their evolutionary rate. To disentangle these hypotheses, we examined patterns of sequence and structural evolution in nuclear‐ and mitochondrial‐encoded oxidative phosphorylation proteins from species in the angiosperm genus Silene with vastly different mitochondrial mutation rates. We found correlated increases in N‐mt gene evolution in species with fast‐evolving mitochondrial DNA. Structural modeling revealed an overrepresentation of N‐mt substitutions at positions that directly contact mutated residues in mitochondrial‐encoded proteins, despite overall patterns of conservative structural evolution. These findings support the hypothesis that selection for compensatory changes in response to mitochondrial mutations contributes to the elevated rate of evolution in N‐mt genes. We discuss these results in light of theories implicating mitochondrial mutation rates and mitonuclear coevolution as drivers of speciation and suggest comparative and experimental approaches that could take advantage of heterogeneity in rates of mtDNA evolution across eukaryotes to evaluate such theories.     

Similar relative mutation rates in the three genetic compartments of Mesostigma and Chlamydomonas

Levels of nucleotide substitution at silent sites in organelle versus nuclear DNAs have been used to estimate relative mutation rates among these compartments and explain lineage-specific features of genome evolution. Synonymous substitution divergence values in animals suggest that the rate of mutation in the mitochondrial DNA is 10-50 times higher than that of the nuclear DNA, whereas overall data for most seed plants support relative mutation rates in mitochondrial, plastid, and nuclear DNAs of 1:3:10. Little is known about relative mutation rates in green algae, as substitution rate data is limited to only the mitochondrial and nuclear genomes of the chlorophyte Chlamydomonas. Here, we measure silent-site substitution rates in the plastid DNA of Chlamydomonas and the three genetic compartments of the streptophyte green alga Mesostigma. In contrast to the situation in animals and land plants, our results support similar relative mutation rates among the three genetic compartments of both Chlamydomonas and Mesostigma. These data are discussed in relation to published intra-species genetic diversity data for the three genetic compartments of Chlamydomonas and are ultimately used to address contemporary hypotheses on the organelle genome evolution. To guide future work, we describe evolutionary divergence data of all publically available Mesostigma viride strains and identify, for the first time, three distinct lineages of Mesostigma.