Accumulation of slightly deleterious mutations in mitochondrial protein-coding genes of large versus small mammals

Nonlinear dynamics of nonsynonymous (dN) and synonymous (dS) substitution rates affects inference of selection

Selection modulates gene sequence evolution in different ways by constraining potential changes of amino acid sequences (purifying selection) or by favoring new and adaptive genetic variants (positive selection). The number of nonsynonymous differences in a pair of protein-coding sequences can be used to quantify the mode and strength of selection. To control for regional variation in substitution rates, the proportionate number of nonsynonymous differences (d(N)) is divided by the proportionate number of synonymous differences (d(S)). The resulting ratio (d(N)/d(S)) is a widely used indicator for functional divergence to identify particular genes that underwent positive selection. With the ever-growing amount of genome data, summary statistics like mean d(N)/d(S) allow gathering information on the mode of evolution for entire species. Both applications hinge on the assumption that d(S) and mean d(S) (approximately branch length) are neutral and adequately control for variation in substitution rates across genes and across organisms, respectively. We here explore the validity of this assumption using empirical data based on whole-genome protein sequence alignments between human and 15 other vertebrate species and several simulation approaches. We find that d(N)/d(S) does not appropriately reflect the action of selection as it is strongly influenced by its denominator (d(S)). Particularly for closely related taxa, such as human and chimpanzee, d(N)/d(S) can be misleading and is not an unadulterated indicator of selection. Instead, we suggest that inconsistencies in the behavior of d(N)/d(S) are to be expected and highlight the idea that this behavior may be inherent to taking the ratio of two randomly distributed variables that are nonlinearly correlated. New null hypotheses will be needed to adequately handle these nonlinear dynamics.

The smallest avian genomes are found in hummingbirds

It has often been suggested that the genome sizes of birds are constrained relative to other tetrapods owing to the high metabolic demands of powered flight and the link between nuclear DNA content and red blood cell size. This hypothesis predicts that hummingbirds, which engage in energy-intensive hovering flight, will display especially constrained genomes even relative to other birds. We report genome size measurements for 37 species of hummingbirds that confirm this prediction. Our results suggest that genome size was reduced before the divergence of extant hummingbird lineages, and that only minimal additional reduction occurred during hummingbird diversification. Unlike in some other avian taxa, the small amount of variation observed within hummingbirds is not explained by variation in respiratory and flight-related parameters. Unexpectedly, genome size appears to have increased in four unrelated hummingbird species whose distributions are centred on humid forests of the upper-tropical elevational zone on the eastern slope of the Andes. This suggests that the secondary expansion of the genome may have been mediated by biogeographical and demographic effects.

Signature of diversifying selection on members of the pentatricopeptide repeat protein family in Arabidopsis lyrata

Pentatricopeptide repeat (PPR) proteins compose a family of nuclear-encoded transcriptional regulators of cytoplasmic genes. They have shown dramatic expansion in copy number in plants, and although the functional importance of many remains unclear, a subset has been repeatedly implicated as nuclear restorers for cytoplasmic male sterility. Here we investigate the molecular population genetics and molecular evolution of seven single-copy PPR genes in the outcrossing model plant Arabidopsis lyrata. In comparison with neutral reference loci, we find, on average, elevated levels of polymorphism and an excess of high-frequency variants at these PPR genes, suggesting that natural selection is maintaining polymorphism at some of these loci. This elevation in diversity persists when we control for divergence and generally decreases in the flanking regions, suggesting that these genes are themselves the targets of selection. Some of the PPR genes also demonstrate elevated population differentiation, which is consistent with spatially varying selection. In contrast, no comparable patterns are observed at these loci in A. thaliana, providing no evidence for the action of balancing selection in this selfing species. Taken together, these results suggest that a subset of PPR genes may be subject to balancing selection associated with ongoing cytonuclear coevolution in the outcrossing A. lyrata, which is possibly mediated either by intergenomic conflict or by compensatory evolution.

Divergent evolution in the cytoplasmic domains of PRLR and GHR genes in Artiodactyla

Background

Prolactin receptor (PRLR) and growth hormone receptor (GHR) belong to the large superfamily of class 1 cytokine receptors. Both of them have been identified as candidate genes affecting key quantitative traits, like growth and reproduction in livestock. We have previously studied the molecular anatomy of the cytoplasmic domain of GHR in different cattle breeds and artiodactyl species. In this study we have analysed the corresponding cytoplasmic signalling region of PRLR.

Results

We sequenced PRLR gene exon 10, coding for the major part of the cytoplasmic domain, from cattle, American bison, European bison, yak, sheep, pig and wild boar individuals. We found different patterns of variation in the two receptors within and between ruminants and pigs. Pigs and bison species have no variation within GHR exon 10, but show high haplotype diversity for the PRLR exon 10. In cattle, PRLR shows lower diversity than GHR. The Bovinae PRLR haplotype network fits better the known phylogenetic relationships between the species than that of the GHR, where differences within cattle breeds are larger than between the different species in the subfamily. By comparison with the wild boar haplotypes, a high number of subsequent nonsynonymous substitutions seem to have accumulated in the pig PRLR exon 10 after domestication.

Conclusion

Both genes affect a multitude of traits that have been targets of selection after domestication. The genes seem to have responded differently to different selection pressures imposed by human artificial selection. The results suggest possible effects of selective sweeps in GHR before domestication in the pig lineage or species divergence in the Bison lineage. The PRLR results may be explained by strong directional selection in pigs or functional switching.

Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability

The evolution of flight is the most important feature of birds, and this ability has helped them become one of the most successful groups of vertebrates. However, some species have independently lost their ability to fly. The degeneration of flight ability is a long process, and some species remain transitional locomotive models. Most of the energy required for locomotion is supplied by mitochondria via oxidative phosphorylation. Thus, rapidly flying birds should require a more energy efficient metabolism than weakly flying or flightless species. Therefore, we speculated that evolutionary constraints acted on the mtDNA of birds with different locomotive abilities. To test this hypothesis, we compared 76 complete avian mitochondrial genomes. Weakly flying and flightless birds, compared to rapidly flying birds, accumulated more nonsynonymous nucleotide substitutions relative to synonymous substitutions. Even after controlling for mutation rate, this trend remained significant. This finding was further tested for its generality by examining 214 complete mammalian mitochondrial genomes. The same as birds, a negative correlation was also found for the K(a)/K(s) ratio and locomotive speed. Our results demonstrated that, in addition to the previously described role for effective population size, functional constraints due to locomotion play an important role in the evolution of mtDNA.

Comparative genomics of Drosophila mtDNA: Novel features of conservation and change across functional domains and lineages

To gain insight on mitochondrial DNA (mtDNA) evolution, we assembled and analyzed the mitochondrial genomes of Drosophila erecta, D. ananassae, D. persimilis, D. willistoni, D. mojavensis, D. virilis and D. grimshawi together with the sequenced mtDNAs of the melanogaster subgroup. Genomic comparisons across the well-defined Drosophila phylogeny impart power for detecting conserved mtDNA regions that maintain metabolic function and regions that evolve uniquely on lineages. Evolutionary rate varies across intergenic regions of the mtDNA. Rapidly evolving intergenic regions harbor the majority of mitochondrial indel divergence. In contrast, patterns of nearly perfect conservation within intergenic regions reveal a refined set of nucleotides underlying the binding of transcription termination factors. Sequencing of 5' cDNA ends indicates that cytochrome C oxidase I (CoI) has a novel (T/C)CG start codon and that perfectly conserved regions upstream of two NADH dehydrogenase (ND) genes are transcribed and likely extend these protein sequences. Substitutions at synonymous sites in the Drosophila mitochondrial proteomes reflect a mutation process that is biased toward A and T nucleotides and differs between mtDNA strands. Differences in codon usage bias across genes reveal that weak selection at silent sites may offset the mutation bias. The mutation-selection balance at synonymous sites has also diverged between the Drosophila and Sophophora lineages. Rates of evolution are highly heterogeneous across the mitochondrial proteome, with ND accumulating many more amino acid substitutions than CO. These oxidative phosphorylation complex-specific rates of evolution vary across lineages and may reflect physiological and ecological change across the Drosophila phylogeny.

The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals

BACKGROUND

During the last ten years, major advances have been made in characterizing and understanding the evolution of mitochondrial DNA, the most popular marker of molecular biodiversity. Several important results were recently reported using mammals as model organisms, including (i) the absence of relationship between mitochondrial DNA diversity and life-history or ecological variables, (ii) the absence of prominent adaptive selection, contrary to what was found in invertebrates, and (iii) the unexpectedly large variation in neutral substitution rate among lineages, revealing a possible link with species maximal longevity. We propose to challenge these results thanks to the bird/mammal comparison. Direct estimates of population size are available in birds, and this group presents striking life-history trait differences with mammals (higher mass-specific metabolic rate and longevity). These properties make birds the ideal model to directly test for population size effects, and to discriminate between competing hypotheses about the causes of substitution rate variation.

RESULTS

A phylogenetic analysis of cytochrome b third-codon position confirms that the mitochondrial DNA mutation rate is quite variable in birds, passerines being the fastest evolving order. On average, mitochondrial DNA evolves slower in birds than in mammals of similar body size. This result is in agreement with the longevity hypothesis, and contradicts the hypothesis of a metabolic rate-dependent mutation rate. Birds show no footprint of adaptive selection on cytochrome b evolutionary patterns, but no link between direct estimates of population size and cytochrome b diversity. The mutation rate is the best predictor we have of within-species mitochondrial diversity in birds. It partly explains the differences in mitochondrial DNA diversity patterns observed between mammals and birds, previously interpreted as reflecting Hill-Robertson interferences with the W chromosome.

CONCLUSION

Mitochondrial DNA diversity patterns in birds are strongly influenced by the wide, unexpected variation of mutation rate across species. From a fundamental point of view, these results are strongly consistent with a relationship between species maximal longevity and mitochondrial mutation rate, in agreement with the mitochondrial theory of ageing. Form an applied point of view, this study reinforces and extends the message of caution previously expressed for mammals: mitochondrial data tell nothing about species population sizes, and strongly depart the molecular clock assumption.

Evidence for variation in the effective population size of animal mitochondrial DNA

Background

It has recently been shown that levels of diversity in mitochondrial DNA are remarkably constant across animals of diverse census population sizes and ecologies, which has led to the suggestion that the effective population of mitochondrial DNA may be relatively constant.

Results

Here we present several lines of evidence that suggest, to the contrary, that the effective population size of mtDNA does vary, and that the variation can be substantial. First, we show that levels of mitochondrial and nuclear diversity are correlated within all groups of animals we surveyed. Second, we show that the effectiveness of selection on non-synonymous mutations, as measured by the ratio of the numbers of non-synonymous and synonymous polymorphisms, is negatively correlated to levels of mitochondrial diversity. Finally, we estimate the effective population size of mitochondrial DNA in selected mammalian groups and show that it varies by at least an order of magnitude.

Conclusions

We conclude that there is variation in the effective population size of mitochondria. Furthermore we suggest that the relative constancy of DNA diversity may be due to a negative correlation between the effective population size and the mutation rate per generation.

Inverse relationship between longevity and evolutionary rate of mitochondrial proteins in mammals and birds

Recently, an unexpected, positive correlation between the rate of evolution of mitochondrial proteins and longevity was reported. Here we re-analyze this relationship in various mammalian lineages using a bayesian phylogenetic analysis of amino-acid sequences, allowing for variable evolutionary rates across sites and species. A negative relationship between protein evolutionary rate and species longevity is reported for all oxidative phosphorylation complexes. A detailed analysis of the cytochrome b in 528 mammals reinforced this result, which contradicts previous publications. Reconducting the analysis in birds yielded similar results. We explain the discrepancy between this and previous reports by our improved taxon sampling and more appropriate methodology: unlike distance-based methods, the tree-based bayesian approach can take into account the high variation of substitution rate across amino-acid sites, and the resulting multiple substitution events. We discuss how our analysis contradicts Rottenberg's rationale, but does not dismiss his proposal of a longevity-dependent selective pressure on mitochondrial mutation rate in mammals and birds. This is because his interpretation invokes adaptation as the single evolutionary force at work, disregarding the effects of mutation, genetic drift, and purifying selection.

Near neutrality, rate heterogeneity, and linkage govern mitochondrial genome evolution in Atlantic cod (Gadus morhua) and other gadine fish

The mitochondrial DNA (mtDNA) genome figures prominently in evolutionary investigations of vertebrate animals due to a suite of characteristics that include absence of Darwinian selection, high mutation rate, and inheritance as a single linkage group. Given complete linkage and selective neutrality, mtDNA gene trees are expected to correspond to intraspecific phylogenies, and mtDNA diversity will reflect population size. The validity of these assumptions is, however, rarely tested on a genome-wide scale. Here, we analyze rates and patterns of molecular evolution among 32 whole mitochondrial genomes of Atlantic Cod (Gadus morhua) as compared with its sister taxon, the walleye pollock (Gadus [Theragra] chalcogrammus), and genomes of seven other gadine codfish. We evaluate selection within G. morhua, between sister species, and among species and intraspecific measures of linkage disequilibrium and recombination within G. morhua. Strong rate heterogeneity occurs among sites and genes at all levels of hierarchical comparison, consistent with variation in mutation rates across the genome. Neutrality indices (dN/dS) are significantly greater than unity among G. morhua genomes and between sister species, which suggests that polymorphisms within species are slightly deleterious, as expected under the nearly neutral theory of molecular evolution. Among species of gadines, dN/dS ratios are heterogeneous among genes, consistent with purifying selection and variation in functional constraint among genes rather than positive selection. The dN/dS ratio for ND4L is anomalously high across all hierarchical levels. There is no evidence for recombination within G. morhua. These patterns contrast strongly with those reported for humans: genome-wide patterns in other vertebrates should be investigated to elucidate the complex patterns of mtDNA molecular evolution.

An evolutionary footprint of age-related natural selection in mitochondrial DNA

By comparing mtDNA sequences between different orders of mammals, we show that both longevity and generation time are significantly correlated with the nucleotide content of the mtDNA. Specifically, there is a positive correlation between generation time and mt GC content. This correlation is repeated, at a finer evolutionary scale, within the primates. Moreover, a comparison of human and chimpanzee mtDNAs shows that the effect has been very pronounced during the short evolutionary period since the divergence of these two species, with human mtDNA showing a GC-biased pattern of substitution at the variable sites. In addition to these DNA sequence patterns, comparisons between the human and the chimp mt protein sequences also revealed a surprisingly high substitution rate for threonine residues, resulting in a reduction of threonine in the human mt proteome. These patterns of both DNA and protein evolution can be explained by a balance between AT-biased mutational pressure and age-related purifying selection.

Mitigating mutational meltdown in mammalian mitochondria

Animal mitochondrial genomes have high rates of sequence evolution, and should decay from the accumulation of deleterious mutations. But the purging of mutant mtDNAs in a pedigree of "mutator mice" reveals the speed and power of purifying selection to maintain mitochondrial function.

Patterns of molecular evolution in Caenorhabditis preclude ancient origins of selfing

The evolution of self-fertilization can mediate pronounced changes in genomes as a by-product of a drastic reduction in effective population size and the concomitant accumulation of slightly deleterious mutations by genetic drift. In the nematode genus Caenorhabditis, a highly selfing lifestyle has evolved twice independently, thus permitting an opportunity to test for the effects of mode of reproduction on patterns of molecular evolution on a genomic scale. Here we contrast rates of nucleotide substitution and codon usage bias among thousands of orthologous groups of genes in six species of Caenorhabditis, including the classic model organism Caenorhabditis elegans. Despite evidence that weak selection on synonymous codon usage is pervasive in the history of all species in this genus, we find little difference among species in the patterns of codon usage bias and in replacement-site substitution. Applying a model of relaxed selection on codon usage to the C. elegans and C. briggsae lineages suggests that self-fertilization is unlikely to have evolved more than approximately 4 million years ago, which is less than a quarter of the time since they shared a common ancestor with outcrossing species. We conclude that the profound changes in mating behavior, physiology, and developmental mechanisms that accompanied the transition from an obligately outcrossing to a primarily selfing mode of reproduction evolved in the not-too-distant past.

Higher origination and extinction rates in larger mammals

Do large mammals evolve faster than small mammals or vice versa? Because the answer to this question contributes to our understanding of how life-history affects long-term and large-scale evolutionary patterns, and how microevolutionary rates scale-up to macroevolutionary rates, it has received much attention. A satisfactory or consistent answer to this question is lacking, however. Here, we take a fresh look at this problem using a large fossil dataset of mammals from the Neogene of the Old World (NOW). Controlling for sampling biases, calculating per capita origination and extinction rates of boundary-crossers and estimating survival probabilities using capture-mark-recapture (CMR) methods, we found the recurring pattern that large mammal genera and species have higher origination and extinction rates, and therefore shorter durations. This pattern is surprising in the light of molecular studies, which show that smaller animals, with their shorter generation times and higher metabolic rates, have greater absolute rates of evolution. However, higher molecular rates do not necessarily translate to higher taxon rates because both the biotic and physical environments interact with phenotypic variation, in part fueled by mutations, to affect origination and extinction rates. To explain the observed pattern, we propose that the ability to evolve and maintain behavior such as hibernation, torpor and burrowing, collectively termed "sleep-or-hide" (SLOH) behavior, serves as a means of environmental buffering during expected and unexpected environmental change. SLOH behavior is more common in some small mammals, and, as a result, SLOH small mammals contribute to higher average survivorship and lower origination probabilities among small mammals.

Molecular evolution of the cytochrome c oxidase subunit 5A gene in primates

BACKGROUND

Many electron transport chain (ETC) genes show accelerated rates of nonsynonymous nucleotide substitutions in anthropoid primate lineages, yet in non-anthropoid lineages the ETC proteins are typically highly conserved. Here, we test the hypothesis that COX5A, the ETC gene that encodes cytochrome c oxidase subunit 5A, shows a pattern of anthropoid-specific adaptive evolution, and investigate the distribution of this protein in catarrhine brains.

RESULTS

In a dataset comprising 29 vertebrate taxa, including representatives from all major groups of primates, there is nearly 100% conservation of the COX5A amino acid sequence among extant, non-anthropoid placental mammals. The most recent common ancestor of these species lived about 100 million years (MY) ago. In contrast, anthropoid primates show markedly elevated rates of nonsynonymous evolution. In particular, branch site tests identify five positively selected codons in anthropoids, and ancestral reconstructions infer that substitutions in these codons occurred predominantly on stem lineages (anthropoid, ape and New World monkey) and on the human terminal branch. Examination of catarrhine brain samples by immunohistochemistry characterizes for the first time COX5A protein distribution in the primate neocortex, and suggests that the protein is most abundant in the mitochondria of large-size projection neurons. Real time quantitative PCR supports previous microarray results showing COX5A is expressed in cerebral cortical tissue at a higher level in human than in chimpanzee or gorilla.

CONCLUSION

Taken together, these results suggest that both protein structural and gene regulatory changes contributed to COX5A evolution during humankind's ancestry. Furthermore, these findings are consistent with the hypothesis that adaptations in ETC genes contributed to the emergence of the energetically expensive anthropoid neocortex.

Life-history traits drive the evolutionary rates of mammalian coding and noncoding genomic elements

A comprehensive phylogenetic framework is indispensable for investigating the evolution of genomic features in mammals as a whole, and particularly in humans. Using the ENCODE sequence data, we estimated mammalian neutral evolutionary rates and selective pressures acting on conserved coding and noncoding elements. We show that neutral evolutionary rates can be explained by the generation time (GT) hypothesis. Accordingly, primates (especially humans), having longer GTs than other mammals, display slower rates of neutral evolution. The evolution of constrained elements, particularly of nonsynonymous sites, is in agreement with the expectations of the nearly neutral theory of molecular evolution. We show that rates of nonsynonymous substitutions (dN) depend on the population size of a species. The results are robust to the exclusion of hypermutable CpG prone sites. The average rate of evolution in conserved noncoding sequences (CNCs) is 1.7 times higher than in nonsynonymous sites. Despite this, CNCs evolve at similar or even lower rates than nonsynonymous sites in the majority of basal branches of the eutherian tree. This observation could be the result of an overall gradual or, alternatively, lineage-specific relaxation of CNCs. The latter hypothesis was supported by the finding that 3 of the 20 longest CNCs displayed significant relaxation of individual branches. This observation may explain why the evolution of CNCs fits the expectations of the nearly neutral theory less well than the evolution of nonsynonymous sites.