Human evolution and the mitochondrial genome

The use of mitochondrial DNA (mtDNA) continues to dominate studies of human genetic variation and evolution. Recent work has re-affirmed the strict maternal inheritance of mtDNA, yielded new insights into the extent and nature of intra-individual variation, supported a recent African origin of human mtDNA, and amply demonstrated the utility of mtDNA in tracing population history and in analyses of ancient remains.     

Adaptive evolution of the Homo mitochondrial genome

Adaptive evolution of 12 protein-coding mitochondrial genes in members of genus Homo (Denisova hominin (H. sp. Altai), Neandertals (H. neanderthalensis) and modern humans (H. sapiens)) has been evaluated by assessing the pattern of changes in the physicochemical properties of amino acid replacements during the primate evolution. It has been found that in the Homo molecular adaptation (positive destabilizing selection) become apparent in the form of 12 radical amino acid replacements accompanied by statistically significant (P < 0.001) changes of physicochemical properties that probably had the functional consequences. These replacements have occurred on the stage of a common ancestor of the Homo (in CO2 and CytB genes) as well as with the appearance of the common ancestor of Neandertals and modern humans (in CO1 and ND5 genes). Radical amino acid replacements were mainly revealed in the cytochrome c oxidase complex IV and cytochrome bc1 complex III, thus coinciding with general trend of increasing of non-synonymous changes in mtDNA genes coding subunits of complexes III and IV proteins in anthropoid primates.     

Nonadaptive evolution of mitochondrial genome size

Genomes vary greatly in size and complexity, and identifying the evolutionary forces that have generated this variation remains a major goal in biology. A controversial proposal is that most changes in genome size are initially deleterious and therefore are linked to episodes of decrease in effective population sizes. Support for this hypothesis comes from large-scale comparative analyses, but vanishes when phylogenetic nonindependence is taken into account. Another approach to test this hypothesis involves analyzing sequence evolution among clades where duplications have recently fixed. Here we show that episodes of fixation of duplications in mitochondrial genomes of the gecko Heteronotia binoei (two independent clades) and of mantellid frogs (five distinct branches) coincide with reductions in the ability of selection to purge slightly deleterious mutations. Our results support the idea that genome complexity can arise through nonadaptive processes in tetrapods.     

Sequential adaptive introgression of the mitochondrial genome in Drosophila yakuba and Drosophila santomea

Interspecific hybridization provides the unique opportunity for species to tap into genetic variation present in a closely related species and potentially take advantage of beneficial alleles. It has become increasingly clear that when hybridization occurs, mitochondrial DNA (mtDNA) often crosses species boundaries, raising the possibility that it could serve as a recurrent target of natural selection and source of species' adaptations. Here we report the sequences of 46 complete mitochondrial genomes of Drosophila yakuba and Drosophila santomea, two sister species known to produce hybrids in nature (~3%). At least two independent events of mtDNA introgression are uncovered in this study, including an early invasion of the D. yakuba mitochondrial genome that fully replaced the D. santomea mtDNA native haplotypes and a more recent, ongoing event centred in the hybrid zone. Interestingly, this recent introgression event bears the signature of Darwinian natural selection, and the selective haplotype can be found at low frequency in Africa mainland populations of D. yakuba. We put forward the possibility that, because the effective population size of D. santomea is smaller than that of D. yakuba, the faster accumulation of mildly deleterious mutations associated with Muller's ratchet in the former species may have facilitated the replacement of the mutationally loaded mitochondrial genome of D. santomea by that of D. yakuba.     

Principles of genome evolution in the Drosophila melanogaster species group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

Modeling the evolution of the human mitochondrial genome.

Mitochondrial DNA data have been used extensively to study evolution and early human origins. These applications require estimates of the rate at which nucleotide substitutions occur in the DNA sequence. We consider the problem of estimating substitution rates in the presence of site-to-site rate variation. A coalescent model is presented that allows for different substitution rates for purines and pyrimidines, as well as more detailed models that allow fast and slow rates within each of the purine and pyrimidine classes. A method for estimating such rates is presented. Even for these simple models of site heterogeneity, there are, typically, insufficient data to obtain reliable estimates of site-specific substitution rates. However, estimates of the average rate across all sites appear to be relatively stable even in the presence of site heterogeneity. Simulations of models with site-to-site variation in mutation rate show that hypervariable sites can produce peaks in the pairwise difference curves that have previously been attributed to population dynamics.     

Organization and evolution of the mitochondrial genome of yeast

Summary The mitochondrial genome of yeast (S. cerevisiae orS. carlsbergensis) appears to be formed by 60–70 genetic units, each one of which is formed by (1) a GC-rich sequence, possibly having a regulatory role; (2) a gene, and (3) an AT-rich spacer, which probably is not transcribed. Recombination in this genome appears to underlie a number of important phenomena. The organization of the mitochondrial genome of yeast and these recombinational events are discussed in relationship with the organization and evolution of the nuclear genome of eukaryotes.     

线粒体基因组的遗传与进化研究进展

线粒体基因组是当前生命科学的热门话题之一。文章根据国内外有关线粒体DNA的结构、表达过程和遗传特征方面的最新研究成果,重点介绍线粒体基因组与核基因组关系、线粒体遗传密码及基因组的进货线索等问题,并简要说明了线粒体基因组遗传分析的要点。     

Variation in synonymous codon use and DNA polymorphism within the Drosophila genome

A strong negative correlation between the rate of amino-acid substitution and codon usage bias in Drosophila has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. To further explore this possibility we have investigated polymorphism and divergence at three kinds of sites: synonymous, nonsynonymous and intronic in relation to codon bias in D. melanogaster and D. simulans. We confirmed that protein evolution is one of the main explicative parameters for interlocus codon bias variation (r(2) approximately 40%). However, intron or synonymous diversities, which could have been expected to be good indicators of local interference [here defined as the additional increase of drift due to selection on tightly linked sites, also called 'genetic draft' by Gillespie (2000)] did not covary significantly with codon bias or with protein evolution. Concurrently, levels of polymorphism were reduced in regions of low recombination rates whereas codon bias was not. Finally, while nonsynonymous diversities were very well correlated between species, neither synonymous nor intron diversities observed in D. melanogaster were correlated with those observed in D. simulans. All together, our results suggest that the selective constraint on the protein is a stable component of gene evolution while local interference is not. The pattern of variation in genetic draft along the genome therefore seems to be instable through evolutionary times and should therefore be considered as a minor determinant of codon bias variance. We argue that selective constraints for optimal codon usage are likely to be correlated with selective constraints on the protein, both between codons within a gene, as previously suggested, and also between genes within a genome.     

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

The complete mitochondrial genome of Fusarium oxysporum: insights into fungal mitochondrial evolution

The complete mitochondrial DNA (mtDNA) sequence was determined for the phytopathogenic fungus Fusarium oxysporum. It is 34,477 bp long, maps circularly, and encodes for 14 protein-coding, 25 tRNA and 2 rRNA genes. The nucleotide and amino acid data sets from its 14 concatenated protein-coding mitochondrial (mt) genes were used along with gene order comparisons for an extensive phylogenetic study of the Subphylum Pezizomycotina. Our results are in agreement with current taxonomic treatments and additionally provide better statistical support for all relationships within Pezizomycotina when compared to analyses based on single or few gene data sets. The gene order of F. oxysporum was consistent with that established in the order Hypocreales (Class: Sordariomycetes) and enhanced previous suppositions on the ancestral state of Sordariomycetes. In comparison with mt genomes of the other orders it added further insights to the evolution of Pezizomycotina.     

Accelerated evolution of the mitochondrial genome in an alloplasmic line of durum wheat

Background

Wheat is an excellent plant species for nuclear mitochondrial interaction studies due to availability of large collection of alloplasmic lines. These lines exhibit different vegetative and physiological properties than their parents. To investigate the level of sequence changes introduced into the mitochondrial genome under the alloplasmic condition, three mitochondrial genomes of the Triticum-Aegilops species were sequenced: 1) durum alloplasmic line with the Ae. longissima cytoplasm that carries the T. turgidum nucleus designated as (lo) durum, 2) the cytoplasmic donor line, and 3) the nuclear donor line.

Results

The mitochondrial genome of the T. turgidum was 451,678 bp in length with high structural and nucleotide identity to the previously characterized T. aestivum genome. The assembled mitochondrial genome of the (lo) durum and the Ae. longissima were 431,959 bp and 399,005 bp in size, respectively. The high sequence coverage for all three genomes allowed analysis of heteroplasmy within each genome. The mitochondrial genome structure in the alloplasmic line was genetically distant from both maternal and paternal genomes. The alloplasmic durum and the Ae. longissima carry the same versions of atp6, nad6, rps19-p, cob and cox2 exon 2 which are different from the T. turgidum parent. Evidence of paternal leakage was also observed by analyzing nad9 and orf359 among all three lines. Nucleotide search identified a number of open reading frames, of which 27 were specific to the (lo) durum line.

Conclusions

Several heteroplasmic regions were observed within genes and intergenic regions of the mitochondrial genomes of all three lines. The number of rearrangements and nucleotide changes in the mitochondrial genome of the alloplasmic line that have occurred in less than half a century was significant considering the high sequence conservation between the T. turgidum and the T. aestivum that diverged from each other 10,000 years ago. We showed that the changes in genes were not limited to paternal leakage but were sufficiently significant to suggest that other mechanisms, such as recombination and mutation, were responsible. The newly formed ORFs, differences in gene sequences and copy numbers, heteroplasmy, and substoichiometric changes show the potential of the alloplasmic condition to accelerate evolution towards forming new mitochondrial genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-67) contains supplementary material, which is available to authorized users.     

The mitochondrial genome of the moss Physcomitrella patens sheds new light on mitochondrial evolution in land plants

The phylogenetic positions of bryophytes and charophytes, together with their genome features, are important for understanding early land plant evolution. Here we report the complete nucleotide sequence (105,340 bp) of the circular-mapping mitochondrial DNA of the moss Physcomitrella patens. Available evidence suggests that the multipartite structure of the mitochondrial genome in flowering plants does not occur in Physcomitrella. It contains genes for 3 rRNAs (rnl, rns, and rrn5), 24 tRNAs, and 42 conserved mitochondrial proteins (14 ribosomal proteins, 4 ccm proteins, 9 nicotinamide adenine dinucleotide dehydrogenase subunits, 5 ATPase subunits, 2 succinate dehydrogenase subunits, apocytochrome b, 3 cytochrome oxidase subunits, and 4 other proteins). We estimate that 5 tRNA genes are missing that might be encoded by the nuclear genome. The overall mitochondrial genome structure is similar in Physcomitrella, Chara vulgaris, Chaetosphaeridium globosum, and Marchantia polymorpha, with easily identifiable inversions and translocations. Significant synteny with angiosperm and chlorophyte mitochondrial genomes was not detected. Phylogenetic analysis of 18 conserved proteins suggests that the moss-liverwort clade is sister to angiosperms, which is consistent with a previous analysis of chloroplast genes but is not consistent with some analyses using mitochondrial sequences. In Physcomitrella, 27 introns are present within 16 genes. Nine of its intron positions are shared with angiosperms and 4 with Marchantia, which in turn shares only one intron position with angiosperms. The phylogenetic analysis as well as the syntenic structure suggest that the mitochondrial genomes of Physcomitrella and Marchantia retain prototype features among land plant mitochondrial genomes.