Purifying selection in mitochondria, free-living and obligate intracellular proteobacteria
© Mamirova et al; licensee BioMed Central Ltd. 2007
Received: 14 September 2006
Accepted: 12 February 2007
Published: 12 February 2007
The effectiveness of elimination of slightly deleterious mutations depends mainly on drift and recombination frequency. Here we analyze the influence of these two factors on the strength of the purifying selection in mitochondrial and proteobacterial orthologous genes taking into account the differences in the organism lifestyles.
(I) We found that the probability of fixation of nonsynonymous substitutions (K n /K s ) in mitochondria is significantly lower compared to obligate intracellular bacteria and even marginally significantly lower compared to free-living bacteria. The comparison of bacteria of different lifestyles demonstrates more effective elimination of slightly deleterious mutations in (II) free-living bacteria as compared to obligate intracellular species and in (III) obligate intracellular parasites as compared to obligate intracellular symbionts. (IV) Finally, we observed that the level of the purifying selection (i.e. 1-K n /K s ) increases with the density of mobile elements in bacterial genomes.
This study shows that the comparison of patterns of molecular evolution of orthologous genes between ecologically different groups of organisms allow to elucidate the genetic consequences of their various lifestyles. Comparing the strength of the purifying selection among proteobacteria with different lifestyles we obtained results, which are in concordance with theoretical expectations: (II) low effective population size and level of recombination in obligate intracellular proteobacteria lead to less effective elimination of mutations compared to free-living relatives; (III) rare horizontal transmissions, i.e. effectively zero recombination level in symbiotic obligate intracellular bacteria leads to less effective purifying selection than in parasitic obligate intracellular bacteria; (IV) the increased frequency of recombination in bacterial genomes with high mobile element density leads to a more effective elimination of slightly deleterious mutations. At the same time, (I) more effective purifying selection in relatively small populations of nonrecombining mitochondria as compared to large populations of recombining proteobacteria was unexpected. We hypothesize that additional features such as the high number of protein-protein interactions or female germ-cell atresia increase evolutionary constraints and maintain the effective purifying selection in mitochondria, but more work is needed to definitely establish these additional features.
Theoretically, elimination of recent slightly deleterious mutations must be more effective in species with high population size [1, 2], since low stochastic sampling variance decreases the power of random genetic drift. Also, according to the theory, meiotic segregation and recombination allow for independent mixing of alleles, thus leading to more effective elimination of slightly deleterious mutations [3–6] and halting of the Muller's ratchet [7, 8]. Thus both the population size and the recombination level determine the effectiveness of the purifying selection.
Although the influences of population size and recombination level on the effectiveness of the purifying selection are interconnected and often mutually reinforcing, it is possible to distinguish between them by comparing species with different population sizes, levels of recombination, or both. For example, the influence of the population size on probability of fixation of slightly deleterious mutations was observed in comparison of primates and rodents [9–11], mammals, birds, and drosophilids , island- versus continent-inhabiting populations of same species , and large- versus small-bodied mammals. The role of the recombination level was demonstrated in asexual versus sexual Daphnia . The influence of both population size and recombination level, most likely, leads to the differences in the molecular evolution patterns of endosymbiotic and free-living bacteria [16–18] and to degeneration of the neo-Y chromosome of Drosophila miranda .
The strength of the purifying selection, acting on different genes, is mainly determined by structural or functional constraints on the encoded protein. Since orthologous proteins in closely related species most likely operate under similar constraints, it is possible to investigate more intimate relationships between the rates of molecular evolution, the population sizes of the analyzed species and the level of recombination in their genomes. Here we perform a comprehensive analysis of purifying selection in various groups of proteobacteria and mammalian mitochondria.
The relative effectiveness of the purifying selection was estimated as the rate of fixation of nonsynonymous substitutions (K n /K s ) and radical amino acid substitutions (K r /K c ) in seven orthologous genes, encoding subunits of a respiratory chain enzyme (NADH:ubiquinone-oxidoreductase or complex I) from bacterial and mitochondrial genomes. These genes encode components of one multi-subunit complex and represent a large fraction of all mitochondrial protein-coding genes. Orthologous genes, encoding subunits of complex I, exist in all aerobic bacteria.
Since the respiratory chain proteins are highly conserved, we assume that the role of the positive selection is negligible, and thus these genes are appropriate for the analysis of the purifying selection. Throughout the work we empirically corroborate this assumption.
Firstly, we compare mammalian mitochondria with bacteria. Since mitochondria are descendants of ancient endosymbiotic proteobacteria, it is interesting to compare the purifying selection in orthologous mitochondrial and proteobacterial genes with regards to the differences of their ecology. Mitochondrial genes of vertebrates are strictly maternally inherited, and thus asexual, and have low population size . Conversely, the bacteria (especially free-living) are characterized by nonzero recombination levels and high population sizes. Therefore, we can expect that mitochondria are subject to less effective purifying selection than bacteria.
Secondly, we consider the purifying selection in various groups of proteobacteria, taking into account their ecology. From the ecological point of view, the analyzed sample of bacterial species consists of two main groups: obligate intracellular species (parasites or mutualists) and free-living species able to exist outside of the host cell. The free-living protebacteria are characterized by higher population size and higher recombination level, and thus one may expect that they are subject to more effective purifying selection as compared to obligate intracellular species [21, 18].
Thirdly, we compare the purifying selection in obligate intracellular symbionts (gamma-proteobacteria) and obligate intracellular parasites (alpha-proteobacteria). It is generally assumed that among obligate intracellular bacteria, symbionts with strictly vertical inheritance (such as obligate dietary endosymbiotic gamma-proteobacteria: Buchnera and Blochmannia, that are required for the survival and reproduction of their insect host) possess lower level of recombination as compared with parasites (such as some alpha-proteobacteria: human pathogen Rickettsia and reproductive parasite of arthropods Wolbachia) which are able to switch from one host to another and thus contact and recombine with a novel gene pool . Thus we test the hypothesis that the transmission differences among obligate intracellular species, determined by the differences in their ecology, might influence the effectiveness of purifying selection through recombination frequency.
Fourthly, we analyze the correlation of the strength of purifying selection with the mobile elements density in bacterial genomes. We expect that the mobile element density in bacteria should positively correlate with the recombination level, as it seems likely that these are two reciprocally reinforcing factors. The inter-genomic recombination level in bacteria is mainly determined by the activity of the mobile elements which able to perform conjugative transfer (conjugative plasmids, conjugative elements) or transduction (prophages of temperate phages) [23, 24]. Number of other mobile elements, without ability to induce recombination, is evolutionary maintained in concordance with the recombination level in bacterial genomes  as it is in eukaryotes [25–33]. Thus we use the density of all mobile elements as a proxy measure of the recombination activity in bacterial genomes.
Unexpectedly, the level of purifying selection in mitochondria is significantly higher as compared with obligate intracellular bacteria and even marginally significantly higher as compared with free-living bacteria. This observation contradicts the widely accepted opinion about degradation of mitochondrial genomes and we discuss some factors, which can increase selection constraints of mitochondrial genes. The comparison of bacteria demonstrates a more effective elimination of mutations in free-living species as compared to obligate intracellular ones with low effective population size and level of recombination. The purifying selection in parasitic obligate intracellular bacteria is more effective as compared to symbiotic ones, possibly due to the higher recombination level in the former. To elucidate the purifying effect of recombination, we demonstrate that the effectiveness of purifying selection (1-K n /K s ) positively correlates with the density of mobile genes, i.e. with the frequency of recombination.
(I) Evolutionary models
Further the bacterial tree was divided into 3 monophyletic sub-trees (alpha-, beta and gamma-proteobacteria), and the mammalian tree was divided into glires and primates sub-trees, and each sub-tree was analyzed separately. For each sub-tree three models were implemented: models 0 and 2a defined as above, and model 2c with K n /K s estimated separately for each external branch and a single K n /K s for all internal branches.
Among any two compared nested models, the one with more degrees of freedom had a significantly higher maximum likelihood (ML) value (i.e. 2ΔlnL was larger with statistical significance of at least a 5% for a χ2 distribution with the number of the degrees of freedom dependent on the model), thus providing a better fit to the analyzed data. For the complete bacterial and mammalian trees: ML (models 2b) > ML (model 2a) > ML (model 0); for sub-trees: ML (models 2c) > ML (models 2a) > ML (model 0).
Since in model 2c the external branches of seven species had no nonsynonymous substitutions (K n = 0 and thus K n /K s = 0), these species were not considered further. Therefore, our final dataset for this model consisted of 110 species, including 36 alpha-proteobacteria (20 free-living and 16 obligate intracellular parasites), 20 free-living beta-proteobacteria, 22 gamma-proteobacteria (16 free-living and 6 obligate intracellular symbionts), 13 glires and 19 primates.
Means and standard errors of the K n /K s values obtained in models 2a, 2b and 2c.
glires 0.043 ± 0.0021, 0.049 ± 0.0042
0.070 ± 0.0017, 0.073 ± 0.0065
primates 0.088 ± 0.0025, 0.089 ± 0.0087
free-living 0.102 ± 0.0022, 0.148 ± 0.0303
beta 0.081 ± 0.0038, 0.048 ± 0.0095
0.148 ± 0.0026), 0.192 ± 0.0300
alpha 0.121 ± 0.0036, 0.243 ± 0.0610
gamma 0.089 ± 0.0038, 0.153 ± 0.0643
obligate 0.305 ± 0.0087, 0.304 ± 0.0689
alpha 0.244 ± 0.0088, 0.224 ± 0.0503
gamma 0.445 ± 0.0203, 0.519 ± 0.2004
(II) Lack of positive selection
For all phylogenetic groups (all mammals, primates, glires, all proteobacteria, α-, β-, γ-proteobacteria) model 7 (K n /K s values assumed to be beta-distributed on interval (0, 1)) had significantly lower maximum likelihood (ML) value than model 8 (one extra class of sites which may have K n /K s values larger than one). However, the K n /K s values of the extra class sites were less then one in all mammals (0.298) and in all proteobacteria (0.384), as well as in each group separately: 0.319 in primates, 0.223 in glires, 0.513 in α-, 0.239 in β- and 0.462 in γ-proteobacteria. Thus we did not observe positive selection in the analyzed genes.
(III) Comparison of mitochondria with bacteria
The K n /K s values of mitochondrial genes are significantly lower than those for all proteobacteria (P = 0.013), as well as for all obligate proteobacteria (P < 0.001) and obligate alpha-proteobacteria (P < 0.001), but only marginally significantly lower than those for all free-living proteobacteria (P = 0.068; Fig. 2). The primate K n /K s values are larger than the glires ones (P < 0.001) and lower than K n /K s for obligate alpha-proteobacteria (P = 0.007). Average values of K n and K s estimated separately for each group are given in Table in additional file 1.
(IV) Bacterial ecology
In the analysis of all proteobacteria, the obligate intracellular proteobacteria have higher K n /K s , as compared with the free-living ones (P = 0.018; Fig. 2). For gamma-proteobacteria, the obligate intracellular ones have higher K n /K s than the free-living ones (P = 0.032), while for alpha-proteobacteria the obligate intracellular have K n /K s similar to that of free-living ones (P = 0.815; Fig. 2).
The K n /K s values of obligate intracellular alpha-proteobacteria (parasitic) are smaller than those of obligate intracellular gamma-proteobacteria (symbiotic) (Fig. 2) at a marginally significant level (P = 0.054). The free-living alpha-proteobacteria have higher K n /K s values as compared to free-living beta-proteobacteria (P = 0.003), and similar to those of free-living gamma-proteobacteria (P = 0.319; Fig. 2). The K n /K s values for free-living gamma-proteobacteria are marginally higher than for beta-proteobacteria (P = 0.080; Fig. 2). Averages values of K n and K s estimated separately for each group are given in additional file 1.
(V) Fixation of radical nonsynonymous substitutions (Kr/Kc)
Two of the four K r /K c ratios (charge-based K r /K c and volume-based K r /K c ) demonstrate significantly higher values in all bacteria compared with all mitochondria (P < 0.001 and P = 0.049, respectively). None of the K r /K c ratios demonstrate significant difference between free-living and obligate intracellular bacteria.
(VI) Influence of branch lengths
Pairwise group comparisons of K n /K s values (estimated in model 2c) for all analyzed species (column "all branches") and for species with external branches of the optimal length (column "optimal branches"). In each comparison, the first mentioned group (before slash) has lower K n /K s average value compared to the second group (after slash). The significance of differences of the averages was estimated using the t-test. The number of species from both groups is in the brackets, respectively. See the text for details.
P = 0.013 (32/80)
P = 0.032 (28/59)
P < 0.001 (32/22)
P < 0.001 (28/15)
P < 0.001 (32/16)
P < 0.001 (28/9)
P = 0.068 (32/56)
P = 0.171 (28/44)
P < 0.001 (13/19)
P < 0.001 (12/16)
free-living proteobacteria/obligate proteobacteria
P = 0.018 (56/22)
P = 0.002 (44/15)
free-living gamma-proteobacteria/obligate gama-proteobacteria
P = 0.032 (16/6)
P = 0.059 (13/6)
free-living alpha-proteobacteria/obligate alpha-proteobacteria
P = 0.815 (20/16)
P = 0.163 (17/9)
obligate alpha-proteobacteria/obligate gamma-proteobacteria
P = 0.054 (16/6)
P = 0.165 (9/6)
(VII) Density of mobile genes
The K n /K s values were negatively correlated with the mobile elements density (hereafter MED: plasmid, prophage and transposable element genes per one megabase of genome size; Kendall tau = -0.305, P = 0.0014, N = 46) (Fig. 3). After exclusion of two species with no mobile elements (Buchnera aphidicola APS and Buchnera aphidicola Sg), we performed linear regression of Ln(MED) on Ln(K n /K s ) that corroborated the observed trend (Ln(K n /K s ) = -1.297 - 0.363×Ln(MED), r2 = 0.138, P = 0.013). Additionally, we took into account nonindependence of data using linear mixed-effects model and still observed significant negative regression (Ln(K n /K s ) = -0.670 - 0.619×Ln(MED), P = 0.0012).
Unexpectedly, mitochondrial genes, despite low effective population sizes and no recombination, seem to eliminate slightly deleterious mutations more effectively than orthologous genes in obligate intracellular and even free-living bacteria which have huge population sizes and nonzero recombination level.
Moreover, we demonstrated that this result cannot be explained by the positive selection events in bacteria. The K n /K s values were < 1 in all branch-specific models (2a, 2b, 2c), where the K n /K s values were estimated separately for some branches or sub-trees (but were averaged across all sites). Further, the K n /K s values were < 1 for all classes of sites analyzed under site-specific models (7, 8) which account for heterogeneous selective pressure on codons (but are averaged across all branches of the phylogeny). However as the adaptive evolution could affect only few sites at several time points, both approaches of averaging rates over sites or over time (over branches) could fail to detect positive selection. Here we argue that the design of our experiment, especially (I) analysis of conserved genes on (II) external branches of the tree (see model 2c), allows us to disregard the effect of possible positive selection. Indeed, since we analyzed only external branches of the tree, possible events of positive selection in the deep-branch phylogeny did not influence our results. On the other hand it is unlikely that events of positive selection on external branches (for example in bacterial species) could be so numerous and so universal as to fully determine the observed trend.
Our result is unlikely to be biased due to the saturation of synonymous substitutions on long branches or to stochastic errors on too short branches, the since analysis of species with optimal branch length demonstrated similar trends in majority of analyses. Besides, the K r /K c analysis that is less sensitive to the saturation effect, demonstrated similar trends.
Although we used the models that incorporate the possibility of different codon usage between the compared groups of species (see Methods, par. (b)), it is still could be possible, that the codon bias effect was not fully accounted for. We can suppose that free-living bacteria with high Ne possess higher K n /K s values (as compared to mitochondria) due to elimination of some synonymous substitutions through selection for more preferable codons. Such selection would lead to underestimation of K s . To test this possibility we estimated the codon usage as a Kullback-Leibler distance or relative entropy , where f is the observed probability distribution of m codons in an analyzed genome (downloaded from the Codon Usage Database ) and f' is the codon frequency under the assumption of uniform synonymous codon usage. If the compared distributions are similar, the relative entropy value should be small. We obtained that the average relative entropy is maximal for genomes of free-living bacteria (mean K = 257, N = 51), intermediate for genomes of obligate bacteria (mean K = 225, N = 19) and minimal for mitochondrial genes (mean K = 191, N = 33). However the comparison of averages did not reveal any significant differences between the groups (the Mann-Whitney U-test) and therefore we concluded that the codon bias should not strongly influence our results.
As we look at only a single sequence from each species and ignore polymorphism within the species, some polymorphic mutations may contribute to K n /K s estimated in our work. For a synonymous substitutions (assuming their neutrality), the expected polymorphism at mutation-drift equilibrium is proportional to the effective population size . At the same time polymorphic slightly-deleterious nonsynonymous substitutions should segregate at a higher frequency in populations with low N e as compared to populations with high N e . Thus the both effects should lead to underestimation of K n /K s values of bacterial genes and overestimation of K n /K s values of mitochondrial genes. Since bacterial K n /K s values, obtained in our work are higher than mitochondrial ones, elimination of the influence of polymorphic substitutions (i.e. increase in bacterial K n /K s values and decrease in mitochondrial K n /K s values) should only strengthen our results.
Thus positive selection, codon usage, stochastic errors or saturation effects and influence of polymorphic substitutions can not provide alternative explanations to the conjecture that purifying selection is more efficient in mammalian mitochondria compared to obligate bacterial species and is at least the same compared to free-living bacterial species.
Most likely, although the mitochondrial genes have the same functions as their bacterial orthologs, they experience stronger evolutionary constraints. There may be several possible causes for that. Firstly, the increased number of subunits in mammalian Complex I (45  versus 14 ) can lead to a larger number of protein-protein interactions and consequently to a slower rate of the protein evolution . Secondly, it has been suggested that effective elimination of deleterious mutations from mtDNA during the female germ-cell atresia is caused by the preferential apoptosis of egg cells with defect mtDNA [42–44]. Thirdly, additional factors such as the protein essentiality or the expression level could increase evolutionary constraints on the mitochondrial genes [45, 46]. Although it seems likely that both protein essentiality as well as the expression level of the mitochondrial genes should be higher compared to their bacterial counterparts, we are not aware of any experimental studies corroborating these suggestions. Thus future experiments in this area may help to identify the main cause of the increased evolutionary stability of the mitochondrial protein-coding genes as compared to their bacterial orthologs.
The widely accepted opinion about degradation of mitochondrial genome is based on theoretical  and empirical backgrounds [48–50]. Since mitochondria possess low effective population size and no recombination, the molecular evolution of mtDNA genes should be associated with a high rate of accumulation of slightly deleterious mutations. A number of studies compared mitochondrial and nuclear genes and corroborated these predictions [48–50]. However, because the genes in nucleus and in mitochondria are not orthologous, this approach is not fully adequate. Here, we compared orthologous genes between mitochondria and various proteobacteria, and observed no degradation of mitochondrial genes. Therefore, we suggest that some other above-mentioned factors have to influence on molecular evolution of mtDNA and maintain highly-effective purifying selection in mitochondria despite their low effective population size and absence of recombination.
Although this study does not refute completely the irreversible accumulation of slightly deleterious mutations in nonrecombining mitochondrial genes, we argue that the rate of the Muller ratchet in modern mammalian mitochondrial genes is at least slower than that in the ecological equivalents of mitochondria, obligate intracellular endosymbionts.
The rate of the Muller ratchet depends on the number of optimal individuals n o with fewest mutations. If n o is small, there is a chance that all no individuals would die without offspring and the ratchet will click round one notch, i.e. the new optimal class of individuals will became slightly worse since it will contain more mutations. Haigh demonstrated that under the additive fitness effect of deleterious mutations, the number of individuals in the optimal class is n o = , where N e is the effective population size, U is the mutational rate per genome per generation, and s is the selection coefficient of deleterious mutations . We argue that no in mitochondria is significantly higher than in obligate intracellular symbionts due to higher s (because of the lower values of K n /K s in mitochondria observed here) and lower or similar U (because of the smaller mitochondrial genome size (~16500 bp), but higher mutational rate per nucleotide per generation). Thus, assuming all other parameters equal (N e and mutational rate per nucleotide per generation), the rate of the Muller ratchet is expected to be lower in mitochondria.
Although two of four K r /K c values demonstrate more effective purifying selection in mammalian mitochondria as compared to bacteria, it seems likely, that the K r /K c values in general are less sensitive to the variations in the purifying selection between species.
Higher K n /K s values in obligate intracellular versus free-living bacteria are consistent with the theory as well as with previous observations. This difference is most likely due to low N e and/or recombination level in the obligate intracellular bacteria.
We observed increased effectiveness of purifying selection in obligate intracellular alpha-proteobacterial parasites as compared to obligate intracellular gamma-proteobacterial symbionts. Since both groups pass through the bottlenecks due to their obligate intracellular lifestyle, differences in effective population sizes are unlikely to cause the observed differences in effectiveness of selection. It has been suggested that, since obligate intracellular parasitic microbes change their hosts more frequently than strictly vertically transmitted symbiotic microbes, parasites enjoy recombination among different colonies during horizontal transmissions . Therefore, the most likely explanation for the observed difference in the effectiveness of selection is the higher recombination level in the obligate intracellular parasitic alpha-proteobacteria.
One recently described genome-level difference between obligate intracellular parasites and endosymbionts is the higher number of mobile elements in the former . The negative regression between the mobile element density and the K n /K s values provides empirical evidence that high level of recombination leads to more effective elimination of slightly deleterious mutations. However, causal relationships of various types of mobile elements and purifying selection seem to be different. (I) The negative regression of transposon density on K n /K s is most likely an effect of recombination on both these values, since recombination increases both transposon density and effectiveness of purifying selection (i.e. 1-K n /K s ). That would mean that there may be no direct causal links between transposon density and purifying selection. (II) However, other mobile elements, such as plasmids and temperate phages, determine recombination events , and therefore their high density causes the increase in the effectiveness of purifying selection. Therefore, recombination induced by mobile elements has an additional role, increasing effectiveness of the purifying selection in bacterial genomes, which is beneficial for the host genome evolution [52, 53].
This study shows that the comparison of patterns of molecular evolution between ecologically different groups elucidates the genetic consequences of various ecological lifestyles. Comparing the strength of the purifying selection among proteobacteria with different lifestyles we obtained results, which are in concordance with theoretical expectations: low effective population size and level of recombination in obligate intracellular proteobacteria lead to less effective elimination of mutations compared to free-living relatives; rare horizontal transmission, i.e. effectively zero recombination level in symbiotic obligate intracellular bacteria leads to less effective purifying selection than in parasitic obligate intracellular bacteria; the high frequency of recombination in bacterial genomes with high mobile element density lead to a more effective elimination of deleterious mutations. At the same time, more effective purifying selection in relatively small populations of asexual mitochondria as compared to large populations of sexual proteobacteria was unexpected since it contradicts the common theory. It seems that additional features such as the high number of protein-protein interactions or female germ-cell atresia maintain the effective purifying selection in mitochondria.
(a) Data and preprocessing
Sequences of seven orthologous genes (mitochondrial: ND1-6 and ND4L; bacterial: NuoA, NuoH, NuoJ, NuoK, NuoL, NuoM, NuoN) of alpha-, beta-, gamma-proteobacteria and mitochondria of primates and glires (Rodentia and Lagomorpha), were downloaded from the National Center for Biotechnology Information database (NCBI)  using blastp  with genes of E. coli K12 as a query. Thus we obtained the data for 83 bacterial and 34 mammalian species. These genes were then translated into amino acid sequences in silico. The obtained amino acid sequences for each gene were aligned using E-INS-i MAFFT model [56, 57] and then reverse translated in silico to obtain the nucleotide sequence alignments. The aligned amino acid and corresponding nucleotide sequences were concatenated into a single amino acid (nucleotide) sequence for each species (3353 codons in length). The concatenated amino acid sequences were used to reconstruct the phylogenetic trees for each monophyletic group (alpha-, beta-, gamma-proteobacteria, rodents and primates) using PHYML [58, 59]. All bacterial trees (for alpha-, beta- and gamma-proteobacteria separately) were rooted using the Homo sapiens mitochondrial sequence as an outgroup. Similarly, the primates and glires trees were rooted with Rickettsia rickettsia as an outgroup. The constructed mammalian and microbial trees are shown in Figure 1. The homologous codon positions that exist only in bacteria (mammals) and consequently are absent in all mammals (respectively, all bacteria) were removed from all sequences. The resulting sequence was 2126 codons in length. The concatenated nucleotide sequences were used to analyze the pattern of nucleotide substitutions and to reconstruct the nucleotide sequences in the internal nodes of the trees.
(b) Ratio of nonsynonymous to synonymous nucleotide substitutions (Kn/Ks)
The codon-based likelihood model suggested by Goldman and Yang  and implemented in the program codeml of PAML package  provides a useful framework for estimation of synonymous and nonsynonymous substitution rates (K n /K s ). The transition/transversion rate bias and the codon frequency bias are found to have significant effects on the estimation of synonymous and nonsynonymous rates, and approximate methods do not adequately account for these factors . Thus, the codon-based likelihood approach is preferable, because this model can easily incorporate these parameters. Because of that, the ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions (K n /K s ) in this study was estimated using the program codeml from the PAML package .
To take into account possible species differences in transition/transversion ratio, this parameter was set as free and thus it was estimated from the sequence data in each model. To take into account possible differences in the codon usage bias, all codon frequencies also were set as free parameters.
To estimate the K n /K s values several models were implemented for the mammalian and bacterial trees: model 0 with a single K n /K s value for all branches of the tree; model 2a with two K n /K s values estimated separately for all external and all internal branches; and model 2b with K n /K s values estimated separately for groups of external branches and a single K n /K s value for all internal branches. In model 2b, bacterial external branches were grouped by ecology (two groups: free-living and obligate intracellular organisms) or by taxonomy and ecology (five groups: free-living and obligate intracellular alpha-proteobacteria (parasitic); free-living beta-proteobacteria; free-living and obligate intracellular gamma-proteobacteria (symbiotic)). Since mammalian mitochondria are obligate intracellular symbionts with strictly vertical inheritance, in model 2b the external branches of the mammalian tree were grouped only by taxonomy (primates and glires).
Since an independent estimate of K n /K s for each external branch of the whole tree would have taken too much processor time, the bacterial tree was divided into 3 monophyletic sub-trees (alpha-, beta and gamma-proteobacteria), and the mammalian tree was divided into sub-trees corresponding to glires and primates, and each sub-tree was analyzed separately. For each sub-tree three models were implemented: models 0 and 2a defined as above, and model 2c with K n /K s estimated separately for each external branch and a single K n /K s for all internal branches. The resulting maximum likelihood (ML) values of the models were compared using the likelihood-ratio test .
(c) Positive selection
As only few codons in a gene may be under positive selection, while the rest of the gene is under purifying selection, the average K n /K s value is likely to be below one despite the action of positive selection on the gene. All above mentioned models (see the previous paragraph) assume a single K n /K s ratio for all codons and thus represent a conservative test of positive selection. Adding another class of codons with different K n /K s values may better describe the distribution of K n /K s across codons and reveal codons under positive selection [63, 64]. For this goal we used models 7 and 8 in the codeml program. A neutral model (7) with K n /Ks values assumed to be beta-distributed on interval (0, 1) was compared with a selection model (8) with one extra class of sites which may have K n /K s values exceeding 1. Positive selection is inferred if the estimate of the K n /K s value of this extra class indeed is larger then one and if the Likelihood Ratio Test (LRT) is significant. The LRT is performed by taking the negative of twice the log-likelihood difference between the nested models (7 and 8) and comparing this to the χ2 distribution with two degrees of freedom.
(d) Ratio of radical to conservative substitutions (Kr/Kc)
Based on the algorithm by Hughes and coauthors , which was modified by to take into account the transition/transversion rate bias , we computed the rates of conservative (K c ) and radical (K r ) nonsynonymous substitutions for nucleotide sequences of modern species and their most recent reconstructed ancestors (the ancestral nucleotide sequences were reconstructed using the method by Yang et al.  implemented in PAML under model 2c). Since the K c and K r rates were small (< 0.3), the Jukes-Cantor formula was used to correct for multiple hits; that is, our K r /K c ratio is identical to the d R /d C ratio in Zhang's notation .
The twenty amino acids were classified into groups in four different ways according to their physicochemical properties. We relied largely on Zhang's classification, were the amino acids differ by charge, polarity and both polarity and volume . Additionally, we used Taylor's classification, that classifies the amino acids by volume alone . Amino acid substitutions within groups (i.e. when ancestral and modern amino acids in homologous sites belong to the same group) were regarded as conserved, while those between groups as radical.
(e) Mobile element genes
For the bacterial species with complete genome sequences available, we determined the number of genes with mobile DNA-related functions. For this we used the automated annotation of the Comprehensive Microbial Resource of the Institute of Genomic Research (TIGR). This annotation classifies genes into nineteen functional-role categories, among which the mobile-DNA category specifies plasmid, prophage and transposable-element genes . We used the number of genes from this category divided by the genome size in megabases as an estimate of the mobile element density in each species.
(f) Statistical analysis
All statistical analyses were done in the R language (R Development Core Team 2005) . The K n /K s values obtained from models 2a and 2b were compared using the standard error (SE) values from PAML . The K n /K s distributions, obtained using model 2c for pairs of groups with different ecology, were compared using the parametric t-test. The relationship of K n /K s and the mobile elements density was analyzed using three approaches: the Kendall rank correlation which is robust to the nonlinearity of a relationship; the ordinal linear regression for loge transformed data; and the regression analysis which takes into account possible non-independence of data. Indeed, since the number of mobile genes may possess some phylogenetic inertia and thus may not be independent among species due to shared ancestry, we performed regression analysis based on the linear mixed-effects model (the function lme in the nlme package  of R language (R Development Core Team 2005)), which explicitly takes into account the hierarchical (nested) structure of comparative-species data . Here, the data are nested into the following two levels: species within families (Rhizobiales, Rickettsiales; Burkholderiales, Neisseriales; Enterobacteriales, Pseudomonadales) and within orders (alpha-, beta-, and gamma-proteobacteria, correspondently).
We are grateful to George Bazykin for useful discussions and editing of manuscript, Sergei V. Shestakov for stimulating discussion of roles of mobile elements in bacteria. The manuscript has greatly benefited from stimulating comments and suggestions of the two anonymous reviewers. This study was partially supported by grants from the Howard Hughes Medical Institute (55005610), INTAS (05-1000008-8028) and the Russian Academy of Science (Program "Molecular and Cellular Biology"). K.P. was supported by the Russian Fund of Basic Research (grants 07-04-00521 and 07-04-01756).
- Wright S: Evolution in Mendelian populations. Genetics. 1931, 16: 97-159.PubMed CentralPubMedGoogle Scholar
- Kimura M: The neutral theory of molecular evolution. 1983, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Kimura M, Maruyama T: The mutation load with epistatic gene interactions in fitness. Genetics. 1966, 54: 1337-1351.PubMed CentralPubMedGoogle Scholar
- Kondrashov AS: Selection against harmful mutations in large sexual and asexual populations. Genet Res Camb. 1982, 40: 325-332.View ArticleGoogle Scholar
- Kondrashov AS: Deleterious mutations and the evolution of sexual reproduction. Nature. 1988, 336: 435-440. 10.1038/336435a0.View ArticlePubMedGoogle Scholar
- Charlesworth B: Mutation-selection balance and the evolutionary advantage of sex and recombination. Genet Res Camb. 1990, 55: 199-221.View ArticleGoogle Scholar
- Muller HJ: The relation of recombination to mutational advance. Mutat Res. 1964, 1: 2-9.View ArticleGoogle Scholar
- Felsenstein J: The evolutionary advantage of recombination. Genetics. 1974, 78: 737-756.PubMed CentralPubMedGoogle Scholar
- Li W-H, Tanimura M, Sharp PM: An evaluation of the molecular clock hypothesis using mammalian DNA sequences. J Mol Evol. 1987, 25: 330-342.View ArticlePubMedGoogle Scholar
- Ohta T: Synonymous and nonsynonymous substitutions in mammalian genes and the nearly neutral theory. J Mol Evol. 1995, 40: 56-63. 10.1007/BF00166595.View ArticlePubMedGoogle Scholar
- Eyre-Walker A, Keightley PD, Smith NGS, Gaffney D: Quantifying the slightly deleterious mutation model of molecular evolution. Mol Biol Evol. 2002, 19: 2142-2149.View ArticlePubMedGoogle Scholar
- Keightley PD, Eyre-Walker A: Deleterious Mutations and the Evolution of Sex. Science. 2000, 290: 331-333. 10.1126/science.290.5490.331.View ArticlePubMedGoogle Scholar
- Woolfit M, Bromham L: Population size and molecular evolution on islands. Proc R Soc B. 2005, 272: 2277-2282. 10.1098/rspb.2005.3217.PubMed CentralView ArticlePubMedGoogle Scholar
- Popadin K, Polishchuk LV, Mamirova L, Knorre D, Gunbin K: Accumulation of slightly deleterious mutations in mitochondrial protein-coding genes: the larger the mammal, the worse the genes. BMC Evol Biol. 2007, 7 (1): 17-10.1186/1471-2148-7-17.PubMed CentralView ArticlePubMedGoogle Scholar
- Paland S, Lynch M: Transitions to asexuality result in excess amino acid substitutions. Science. 2006, 311: 990-992. 10.1126/science.1118152.View ArticlePubMedGoogle Scholar
- Wernegreen JJ, Moran NA: Evidence for genetic drift in endosymbionts (Buchnera): analyses of protein-coding genes. Mol Biol Evol. 1999, 16: 83-97.View ArticlePubMedGoogle Scholar
- Wernegreen JJ: Genome evolution in bacterial endosymbionts of insects. Nature Reviews Genetics. 2002, 3: 850-861. 10.1038/nrg931.View ArticlePubMedGoogle Scholar
- Woolfit M, Bromham L: Increased rates of sequence evolution in endosymbiotic bacteria and fungi with small effective population sizes. Mol Biol Evol. 2003, 20: 1545-1555. 10.1093/molbev/msg167.View ArticlePubMedGoogle Scholar
- Bachtrog D: Sex chromosome evolution: molecular aspects of Y-chromosome degeneration in Drosophila. Genome Res. 2005, 15: 1393-1401. 10.1101/gr.3543605.PubMed CentralView ArticlePubMedGoogle Scholar
- Birky CW: Uniparental inheritance of mitochondrial and chloroplast genes: Mechanisms and evolution. Proc Nat Acad Sci USA. 1995, 92: 11331-11338. 10.1073/pnas.92.25.11331.PubMed CentralView ArticlePubMedGoogle Scholar
- Moran NA: Accelerated evolution and Muller's ratchet in endosymbiotic bacteria. Proc Nat Acad Sci USA. 1996, 93: 2873-2878. 10.1073/pnas.93.7.2873.PubMed CentralView ArticlePubMedGoogle Scholar
- Bordenstein SR, Reznikoff WS: Mobile DNA in obligate intracellular bacteria. Nature Reviews/Microbiology. 2005, 3: 688-699.PubMedGoogle Scholar
- Levin BR: The evolution of sex in bacteria. Evolution of sex: An Examination of Current Ideas. 1988, Michod, Levin, Sinauer Associates Inc., U.S, 194-211.Google Scholar
- Thomas CM, Nielsen KM: Mechanisms of, and barriers to, horizontal gene transfer between bacteria. Nature reviews microbiology. 2005, 3: 711-721. 10.1038/nrmicro1234.View ArticlePubMedGoogle Scholar
- Hickey DA: Selfish DNA: a sexually – transmitted nuclear parasite. Genetics. 1982, 101: 519-531.PubMed CentralPubMedGoogle Scholar
- Zeyl C, Bell G: Symbiotic DNA in eukaryotic genomes. Trends in Ecology and Evolution. 1996, 11: 10-15. 10.1016/0169-5347(96)81058-5.View ArticlePubMedGoogle Scholar
- Burt A, Trivers R: Selfish DNA and breeding system in flowering plants. Proc R Soc Lond B. 1998, 265: 141-146. 10.1098/rspb.1998.0275.View ArticleGoogle Scholar
- Duret L, Marais G, Biemont C: Transposons, but not retrotransposons are located preferentially in regions of high recombination rate in Caenorhabiditis elegans. Genetics. 2000, 156: 1661-1669.PubMed CentralPubMedGoogle Scholar
- Arkhipova I, Meselson M: Transposable elements in sexual and ancient asexual taxa. Proc Natl Acad Sci U S A. 2000, 98: 14473-14477. 10.1073/pnas.97.26.14473.View ArticleGoogle Scholar
- Schön I, Martens K: Transposable elements and asexual reproduction. TREE. 2000, 15: 287-288.PubMedGoogle Scholar
- Schoen I, Martens K: Are ancient asexuals less burdened? Selfish DNA, transposons and reproductive mode. Journal of Natural History. 2002, 36: 379-390. 10.1080/00222930110089148.View ArticleGoogle Scholar
- Sullender BW, Crease TJ: The behavior of Daphnia pulex transposable elements in cyclically and obligately parthenogenetic populations. J Mol Evol. 2001, 53: 63-69.PubMedGoogle Scholar
- Nuzhdin SV, Petrov DA: Transposable elements in clonal lineages: lethal hangover from sex. Biological Journal of the Linnean Society. 2003, 79: 33-41. 10.1046/j.1095-8312.2003.00188.x.View ArticleGoogle Scholar
- Li W-H: Molecular evolution. 1997, University of Texas, Health Science Center at Houston: Sinauer Associates Inc., U.SGoogle Scholar
- Wyckoff GJ, Malcom CM, Vallender EJ, Lahn BT: A highly unexpected strong correlation between fixation probability of nonsynonymous mutations and mutation rate. Trends in Genetics. 2005, 21: 381-390. 10.1016/j.tig.2005.05.005.View ArticlePubMedGoogle Scholar
- Swanson WJ, Nielsen R, Yang Q: Pervasive Adaptive Evolution in Mammalian Fertilization Proteins. Mol Biol Evol. 2003, 20: 18-20.View ArticlePubMedGoogle Scholar
- Lemay P: Relative entropy (Kullback-Leibler distance). Download date is 05. 11. 2006, [http://tecfa.unige.ch/~lemay/thesis/THX-Doctorat/node159.html]
- Nakamura Y, Gojobori T, Ikemura T: Codon usage tabulated from the international DNA sequence databases: status for the year 2000. Nucl Acids Res. 2000, 28: 292-10.1093/nar/28.1.292.PubMed CentralView ArticlePubMedGoogle Scholar
- Carroll J, Fearnley IM, Skehel JM, Shannon RJ, Hirst J, Walker JE: Bovine complex I is a complex of forty-five different subunits. J Biol Chem.Google Scholar
- Yagi T: The bacterial energy-transducing NADH-quinone oxidoreductases. Biochim Biophys Acta. 1993, 1141: 1-17. 10.1016/0005-2728(93)90182-F.View ArticlePubMedGoogle Scholar
- Fraser HB, Wall DP, Hirsh AE: A simple dependence between protein evolution rate and the number of protein-protein interactions. BMC Evolutionary Biology. 2003, 3: 11-16. 10.1186/1471-2148-3-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Krakauer DC, Mira A: Mitochondria and germ-cell death. Nature. 1999, 400: 125-126. 10.1038/22026.View ArticlePubMedGoogle Scholar
- Perez GI, Trbovich AM, Gosden RG, Tilly JL: Mitochondria and the death of oocytes. Nature. 2000, 403: 500-501. 10.1038/35000651.View ArticlePubMedGoogle Scholar
- Bergstrom CT, Pritchard J: Germline bottlenecks and evolutionary maintenance og mitochondrial genomes. Genetics. 1998, 149: 2135-2146.PubMed CentralPubMedGoogle Scholar
- Hirsh AE, Fraser HB: Protein dispensability and rate of evolution. Nature. 2001, 411: 1046-1049. 10.1038/35082561.View ArticlePubMedGoogle Scholar
- Pal C, Papp B, Hurst LD: Highly Expressed Genes in Yeast Evolve Slowly. Genetics. 2000, 158: 927-931.Google Scholar
- Hoekstra RF: Evolutionary origin and consequences of uniparental mitochondrial inheritance. Human Reproduction. 2000, 15: 102-111.View ArticlePubMedGoogle Scholar
- Lynch M: Mutation accumulation in transfer RNAs: molecular evidence for Muller's ratchet in mitochondrial genomes. Evolution. 1996, 13: 209-220.Google Scholar
- Lynch M: Mutation accumulation in nuclear, organelle, and prokaryotic transfer RNA genes. Evolution. 1997, 14: 914-925.Google Scholar
- Lynch M, Blanchard JL: Deleterious mutation accumulation in organelle genomes. Genetica. 1998, 102–103: 29-39. 10.1023/A:1017022522486.View ArticlePubMedGoogle Scholar
- Haigh J: The accumulation of deleterious genes in a population – Muller's ratchet. Theor Popul Biol. 1978, 14 (2): 251-67. 10.1016/0040-5809(78)90027-8.View ArticlePubMedGoogle Scholar
- Frost LS, Leplae R, Summers AO, Toussaint A: Mobile genetic elements: the agents of open source evolution. Nature Reviews Microbiology. 2005, 3: 722-732. 10.1038/nrmicro1235.View ArticlePubMedGoogle Scholar
- Kidwell MG, Lisch DR: Transposable elements and host genome evolution. TREE. 2000, 15: 95-99.PubMedGoogle Scholar
- NCBI. Download date is January 2006, [http://ncbi.nlm.nih.gov/genomes/ORGANELLES/animalabout.html]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticlePubMedGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Research. 2002, 30: 3059-3066. 10.1093/nar/gkf436.PubMed CentralView ArticlePubMedGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Research. 2005, 33: 511-518. 10.1093/nar/gki198.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML online – a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Research. 2005, 33: 557-559. 10.1093/nar/gki352.View ArticleGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Molecular Biology and Evolution. 1994, 11: 725-736.PubMedGoogle Scholar
- Yang ZH: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13: 555-556.PubMedGoogle Scholar
- Yang Z, Nielsen R: Synonymous and nonsynonymous rate variation in nuclear genes of mammals. J Mol Evol. 1998, 46: 409-418. 10.1007/PL00006320.View ArticlePubMedGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AK: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 431-449.PubMed CentralPubMedGoogle Scholar
- Yang Z, Swanson WJ: Codon-substitution models to detect adaptive evolution that account for heterogeneous selective pressures among site classes. Mol Biol Evol. 2002, 19: 49-57.View ArticlePubMedGoogle Scholar
- Hughes AL, Ota T, Nei M: Positive Darwinian selection promotes charge profile diversity in the antigen-binding cleft of class I major-histocompatibility-complex molecules. Mol Biol Evol. 1990, 7: 515-524.PubMedGoogle Scholar
- Zhang J: Rates of conservative and radical nonsynonymous nucleotide substitutions in mammalian nuclear genes. J Mol Evol. 2000, 50: 56-68.PubMedGoogle Scholar
- Yang Z, Goldman N, Friday AE: Maximum likelihood trees from DNA sequences: a peculiar statistical estimation problem. Syst Biol. 1995, 44: 384-399. 10.2307/2413599.View ArticleGoogle Scholar
- Taylor WR: The classification of amino acid conservation. J Theor Biol. 1986, 119: 205-218. 10.1016/S0022-5193(86)80075-3.View ArticlePubMedGoogle Scholar
- Peterson JD, Umayam LA, Dickinson T, Hickey EK, White O: The comprehensive microbial resource. Nucleic Acids Res. 2001, 29: 123-125. 10.1093/nar/29.1.123.PubMed CentralView ArticlePubMedGoogle Scholar
- R Development Core Team 2005 R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, [http://www.R-project.org]
- Pinheiro J, Bates D, DebRoy S, Sarkar D: nlme: Linear and nonlinear mixed effects models. 2005, R package version 3.1-60Google Scholar
- Pinheiro JC, Bates DM: Mixed Effects Models in S and S-PLUS (Statistics & Computing S). 2000, New York: Springer VerlagView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.