Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Similar Efficacies of Selection Shape Mitochondrial and Nuclear Genes in Both Drosophila melanogaster and Homo sapiens

Deleterious mutations contribute to polymorphism even when selection effectively prevents their fixation. The efficacy of selection in removing deleterious mitochondrial mutations from populations depends on the effective population size (N_e) of the mitochondrial DNA and the degree to which a lack of recombination magnifies the effects of linked selection. Using complete mitochondrial genomes from Drosophila melanogaster and nuclear data available from the same samples, we reexamine the hypothesis that nonrecombining animal mitochondrial DNA harbor an excess of deleterious polymorphisms relative to the nuclear genome. We find no evidence of recombination in the mitochondrial genome, and the much-reduced level of mitochondrial synonymous polymorphism relative to nuclear genes is consistent with a reduction in N_e. Nevertheless, we find that the neutrality index, a measure of the excess of nonsynonymous polymorphism relative to the neutral expectation, is only weakly significantly different between mitochondrial and nuclear loci. This difference is likely the result of the larger proportion of beneficial mutations in X-linked relative to autosomal loci, and we find little to no difference between mitochondrial and autosomal neutrality indices. Reanalysis of published data from Homo sapiens reveals a similar lack of a difference between the two genomes, although previous studies have suggested a strong difference in both species. Thus, despite a smaller N_e, mitochondrial loci of both flies and humans appear to experience similar efficacies of purifying selection as do loci in the recombining nuclear genome.

Time dependency of foamy virus evolutionary rate estimates

Background

It appears that substitution rate estimates co-vary very strongly with their timescale of measurement; the shorter the timescale, the higher the estimated value. Foamy viruses have a long history of co-speciation with their hosts, and one of the lowest estimated rates of evolution among viruses. However, when their rate of evolution is estimated over short timescales, it is more reminiscent of the rapid rates seen in other RNA viruses. This discrepancy between their short-term and long-term rates could be explained by the time-dependency of substitution rate estimates. Several empirical models have been proposed and used to correct for the time-dependent rate phenomenon (TDRP), such as a vertically-translated exponential rate decay model and a power-law rate decay model. Nevertheless, at present, it is still unclear which model best describes the rate dynamics. Here, we use foamy viruses as a case study to empirically describe the phenomenon and to determine how to correct rate estimates for its effects. Four empirical models were investigated: (i) a vertically-translated exponential rate decay model, (ii) a simple exponential rate decay model, (iii) a vertically-translated power-law rate decay model, and (iv) a simple power-law rate decay model.

Results

Our results suggest that the TDRP is likely responsible for the large discrepancy observed in foamy virus short-term and long-term rate estimates, and the simple power-law rate decay model is the best model for inferring evolutionary timescales. Furthermore, we demonstrated that, within the Bayesian phylogenetic framework, currently available molecular clocks can severely bias evolutionary date estimates, indicating that they are inadequate for correcting for the TDRP. Our analyses also suggest that different viral lineages may have different TDRP dynamics, and this may bias date estimates if it is unaccounted for.

Conclusions

As evolutionary rate estimates are dependent on their measurement timescales, their values must be used and interpreted under the context of the timescale of rate estimation. Extrapolating rate estimates across large timescales for evolutionary inferences can severely bias the outcomes. Given that the TDRP is widespread in nature but has been noted only recently the estimated timescales of many viruses may need to be reconsidered and re-estimated. Our models could be used as a guideline to further improve current phylogenetic inference tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0408-z) contains supplementary material, which is available to authorized users.

Signs of Selection in Synonymous Sites of the Mitochondrial Cytochrome b Gene of Baikal Oilfish (Comephoridae) by mRNA Secondary Structure Alterations

Studies over the past decade have shown a significant role of synonymous mutations in posttranscriptional regulation of gene expression, which is particularly associated with messenger RNA (mRNA) secondary structure alterations. Most studies focused on prokaryote genomes and the nuclear genomes of eukaryotes while little is known about the regulation of mitochondrial DNA (mtDNA) gene expression. This paper reveals signs of selection in synonymous sites of the mitochondrial cytochrome b gene (Cytb) of Baikal oilfish or golomyankas (Comephoridae) directed towards altering the secondary structure of the mRNA and probably altering the character of mtDNA gene expression. Our findings are based on comparisons of intraspecific genetic variation patterns of small golomyanka (Comephorus dybowski) and two genetic groups of big golomyanka (Comephorus dybowskii). Two approaches were used: (i) analysis of the distribution of synonymous mutations between weak-AT (W) and strong-GC (S) nucleotides within species and groups in accordance with mutation directions from central to peripheral haplotypes and (ii) approaches based on the predicted mRNA secondary structure.

Selection against heteroplasmy explains the evolution of uniparental inheritance of mitochondria

Why are mitochondria almost always inherited from one parent during sexual reproduction? Current explanations for this evolutionary mystery include conflict avoidance between the nuclear and mitochondrial genomes, clearing of deleterious mutations, and optimization of mitochondrial-nuclear coadaptation. Mathematical models, however, fail to show that uniparental inheritance can replace biparental inheritance under any existing hypothesis. Recent empirical evidence indicates that mixing two different but normal mitochondrial haplotypes within a cell (heteroplasmy) can cause cell and organism dysfunction. Using a mathematical model, we test if selection against heteroplasmy can lead to the evolution of uniparental inheritance. When we assume selection against heteroplasmy and mutations are neither advantageous nor deleterious (neutral mutations), uniparental inheritance replaces biparental inheritance for all tested parameter values. When heteroplasmy involves mutations that are advantageous or deleterious (non-neutral mutations), uniparental inheritance can still replace biparental inheritance. We show that uniparental inheritance can evolve with or without pre-existing mating types. Finally, we show that selection against heteroplasmy can explain why some organisms deviate from strict uniparental inheritance. Thus, we suggest that selection against heteroplasmy explains the evolution of uniparental inheritance.

Mitochondria contain genes that encode the machinery needed to power cells. Unlike the nuclear genome, the mitochondrial genome is typically inherited from one parent only (uniparental inheritance). The most common explanation for uniparental inheritance is the genomic conflict theory, which states that uniparental inheritance evolved to prevent the spread of ‘selfish’ mitochondria that replicate quickly but produce energy inefficiently. Current explanations have a major problem: when using realistic parameters, mathematical models cannot show that uniparental inheritance can replace biparental inheritance. Clearly, we need a new explanation that fits with standard population-genetic theory. Recent evidence suggests cells may incur a cost when they carry multiple types of mitochondria. Here we show mathematically that uniparental inheritance could have evolved to avoid the costs of maintaining multiple mitochondrial lineages within a cell. Our results explain the long-standing evolutionary mystery of uniparental inheritance and provide insight into the evolution of mating types and binary sexes. Selection against heteroplasmy also has implications for the evolution of the mitochondrial genome because new mitochondrial haplotypes always lead to heteroplasmy before becoming fixed in the population. Thus, selection against heteroplasmy may explain why mtDNA coding-genes have slower substitution rates than analogous genes within the nucleus.

Maternal ancestry and population history from whole mitochondrial genomes

MtDNA has been a widely used tool in human evolutionary and population genetic studies over the past three decades. Its maternal inheritance and lack of recombination have offered the opportunity to explore genealogical relationships among individuals and to study the frequency differences of matrilineal clades among human populations at continental and regional scales. The whole mtDNA genome sequencing delivers molecular resolution that is sufficient to distinguish patterns that have arisen over thousands of years. However, mutation rate is highly variable among the functional and non-coding domains of mtDNA which makes it challenging to obtain accurate split dates of the mitochondrial clades. Due to the shallow coalescent time of mitochondrial TMRCA at approximately 100 to 200 thousand years (ky), mtDNA data have only limited power to inform us about the more distant past and the early stages of human evolutionary history. The variation shared by mitochondrial genomes of individuals drawn from different continents outside Africa has been used to illuminate the details of the colonization process of the Old World, whereas regional patterns of variation have been at the focus of studies addressing questions of a more recent time scale. In the era of whole nuclear genome sequencing, mitochondrial genomes are continuing to be informative as a unique tool for the assessment of female-specific aspects of the demographic history of human populations.

Asymmetric Context-Dependent Mutation Patterns Revealed through Mutation-Accumulation Experiments

Despite the general assumption that site-specific mutation rates are independent of the local sequence context, a growing body of evidence suggests otherwise. To further examine context-dependent patterns of mutation, we amassed 5,645 spontaneous mutations in wild- type (WT) and mismatch-repair deficient (MMR(-)) mutation-accumulation (MA) lines of the gram-positive model organism Bacillus subtilis. We then analyzed>7,500 spontaneous base-substitution mutations across B. subtilis, Escherichia coli, and Mesoplasma florum WT and MMR(-) MA lines, finding a context-dependent mutation pattern that is asymmetric around the origin of replication. Different neighboring nucleotides can alter site-specific mutation rates by as much as 75-fold, with sites neighboring G:C base pairs or dimers involving alternating pyrimidine-purine and purine-pyrimidine nucleotides having significantly elevated mutation rates. The influence of context-dependent mutation on genome architecture is strongest in M. florum, consistent with the reduced efficiency of selection in organisms with low effective population size. If not properly accounted for, the disparities arising from patterns of context-dependent mutation can significantly influence interpretations of positive and purifying selection.

Maternal age effect and severe germ-line bottleneck in the inheritance of human mitochondrial DNA

The manifestation of mitochondrial DNA (mtDNA) diseases depends on the frequency of heteroplasmy (the presence of several alleles in an individual), yet its transmission across generations cannot be readily predicted owing to a lack of data on the size of the mtDNA bottleneck during oogenesis. For deleterious heteroplasmies, a severe bottleneck may abruptly transform a benign (low) frequency in a mother into a disease-causing (high) frequency in her child. Here we present a high-resolution study of heteroplasmy transmission conducted on blood and buccal mtDNA of 39 healthy mother-child pairs of European ancestry (a total of 156 samples, each sequenced at ∼20,000× per site). On average, each individual carried one heteroplasmy, and one in eight individuals carried a disease-associated heteroplasmy, with minor allele frequency ≥1%. We observed frequent drastic heteroplasmy frequency shifts between generations and estimated the effective size of the germ-line mtDNA bottleneck at only ∼30-35 (interquartile range from 9 to 141). Accounting for heteroplasmies, we estimated the mtDNA germ-line mutation rate at 1.3 × 10(-8) (interquartile range from 4.2 × 10(-9) to 4.1 × 10(-8)) mutations per site per year, an order of magnitude higher than for nuclear DNA. Notably, we found a positive association between the number of heteroplasmies in a child and maternal age at fertilization, likely attributable to oocyte aging. This study also took advantage of droplet digital PCR (ddPCR) to validate heteroplasmies and confirm a de novo mutation. Our results can be used to predict the transmission of disease-causing mtDNA variants and illuminate evolutionary dynamics of the mitochondrial genome.

Reproductive clonality in protozoan pathogens--truth or artefact?

The debate around the frequency and importance of genetic exchange in parasitic protozoa is now several decades old. Recently, fresh assertions have been made that predominant clonal evolution explains the population structures of several key protozoan pathogens. Here, we present an alternative perspective. On the assumption that much apparent clonality may be an artefact of inadequate sampling and study design, we review current research to define why sex might be so difficult to detect in protozoan parasite populations. In doing so, we contrast laboratory models of genetic exchange in parasitic protozoa with natural patterns of genetic diversity and consider the fitness advantage of sex at different evolutionary scales. We discuss approaches to improve the accuracy of efforts to characterize genetic exchange in the field. We also examine the implications of the first population genomic studies for the debate around sex and clonality in parasitic protozoa and discuss caveats for the future.

Keeping mtDNA in shape between generations

Since the unexpected discovery that mitochondria contain their own distinct DNA molecules, studies of the mitochondrial DNA (mtDNA) have yielded many surprises. In animals, transmission of the mtDNA genome is explicitly non-Mendelian, with a very high number of genome copies being inherited from the mother after a drastic bottleneck. Recent work has begun to uncover the molecular details of this unusual mode of transmission. Many surprising variations in animal mitochondrial biology are known; however, a series of recent studies have identified a core of evolutionarily conserved mechanisms relating to mtDNA inheritance, e.g., mtDNA bottlenecks during germ cell development, selection against specific mtDNA mutation types during maternal transmission, and targeted destruction of sperm mitochondria. In this review, we outline recent literature on the transmission of mtDNA in animals and highlight the implications for human health and ageing.

On the adaptive value of cytoplasmic genomes in plants

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Population genomic analysis reveals highly conserved mitochondrial genomes in the yeast species Lachancea thermotolerans

The increasing availability of mitochondrial (mt) sequence data from various yeasts provides a tool to study genomic evolution within and between different species. While the genomes from a range of lineages are available, there is a lack of information concerning intraspecific mtDNA diversity. Here, we analyzed the mt genomes of 50 strains from Lachancea thermotolerans, a protoploid yeast species that has been isolated from several locations (Europe, Asia, Australia, South Africa, and North / South America) and ecological sources (fruit, tree exudate, plant material, and grape and agave fermentations). Protein-coding genes from the mtDNA were used to construct a phylogeny, which reflected a similar, yet less resolved topology than the phylogenetic tree of 50 nuclear genes. In comparison to its sister species Lachancea kluyveri, L. thermotolerans has a smaller mt genome. This is due to shorter intergenic regions and fewer introns, of which the latter are only found in COX1. We revealed that L. kluyveri and L. thermotolerans share similar levels of intraspecific divergence concerning the nuclear genomes. However, L. thermotolerans has a more highly conserved mt genome with the coding regions characterized by low rates of nonsynonymous substitution. Thus, in the mt genomes of L. thermotolerans, stronger purifying selection and lower mutation rates potentially shape genome diversity in contract to what was found for L. kluyveri, demonstrating that the factors driving mt genome evolution are different even between closely related species.

The coevolutionary period of Wolbachia pipientis infecting Drosophila ananassae and its impact on the evolution of the host germline stem cell regulating genes

The endosymbiotic bacteria Wolbachia pipientis is known to infect a wide range of arthropod species yet less is known about the coevolutionary history it has with its hosts. Evidence of highly identical W. pipientis strains in evolutionary divergent hosts suggests horizontal transfer between hosts. For example, Drosophila ananassae is infected with a W. pipientis strain that is nearly identical in sequence to a strain that infects both D. simulans and D. suzukii, suggesting recent horizontal transfer among these three species. However, it is unknown whether the W. pipientis strain had recently invaded all three species or a more complex infectious dynamic underlies the horizontal transfers. Here, we have examined the coevolutionary history of D. ananassae and its resident W. pipientis to infer its period of infection. Phylogenetic analysis of D. ananassae mitochondrial DNA and W. pipientis DNA sequence diversity revealed the current W. pipientis infection is not recent. In addition, we examined the population genetics and molecular evolution of several germline stem cell (GSC) regulating genes of D. ananassae. These studies reveal significant evidence of recent and long-term positive selection at stonewall in D. ananassae, whereas pumillio showed patterns of variation consistent with only recent positive selection. Previous studies had found evidence for adaptive evolution of two key germline differentiation genes, bag of marbles (bam) and benign gonial cell neoplasm (bgcn), in D. melanogaster and D. simulans and proposed that the adaptive evolution at these two genes was driven by arms race between the host GSC and W. pipientis. However, we did not find any statistical departures from a neutral model of evolution for bam and bgcn in D. ananassae despite our new evidence that this species has been infected with W. pipientis for a period longer than the most recent infection in D. melanogaster. In the end, analyzing the GSC regulating genes individually showed two of the seven genes to have evidence of selection. However, combining the data set and fitting a specific population genetic model significant proportion of the nonsynonymous sites across the GSC regulating genes were driven to fixation by positive selection. Clearly the GSC system is under rapid evolution and potentially multiple drivers are causing the rapid evolution.

Regulation of mitochondrial genome inheritance by autophagy and ubiquitin-proteasome system: implications for health, fitness, and fertility

Mitochondria, the energy-generating organelles, play a role in numerous cellular functions including adenosine triphosphate (ATP) production, cellular homeostasis, and apoptosis. Maternal inheritance of mitochondria and mitochondrial DNA (mtDNA) is universally observed in humans and most animals. In general, high levels of mitochondrial heteroplasmy might contribute to a detrimental effect on fitness and disease resistance. Therefore, a disposal of the sperm-derived mitochondria inside fertilized oocytes assures normal preimplantation embryo development. Here we summarize the current research and knowledge concerning the role of autophagic pathway and ubiquitin-proteasome-dependent proteolysis in sperm mitophagy in mammals, including humans. Current data indicate that sperm mitophagy inside the fertilized oocyte could occur along multiple degradation routes converging on autophagic clearance of paternal mitochondria. The influence of assisted reproductive therapies (ART) such as intracytoplasmic sperm injection (ICSI), mitochondrial replacement (MR), and assisted fertilization of oocytes from patients of advanced reproductive age on mitochondrial function, inheritance, and fitness and for the development and health of ART babies will be of particular interest to clinical audiences. Altogether, the study of sperm mitophagy after fertilization has implications in the timing of evolution and developmental and reproductive biology and in human health, fitness, and management of mitochondrial disease.

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.

Emergence of resistance to atovaquone-proguanil in malaria parasites: insights from computational modeling and clinical case reports

The usefulness of atovaquone-proguanil (AP) as an antimalarial treatment is compromised by the emergence of atovaquone resistance during therapy. However, the origin of the parasite mitochondrial DNA (mtDNA) mutation conferring atovaquone resistance remains elusive. Here, we report a patient-based stochastic model that tracks the intrahost emergence of mutations in the multicopy mtDNA during the first erythrocytic parasite cycles leading to the malaria febrile episode. The effect of mtDNA copy number, mutation rate, mutation cost, and total parasite load on the mutant parasite load per patient was evaluated. Computer simulations showed that almost any infected patient carried, after four to seven erythrocytic cycles, de novo mutant parasites at low frequency, with varied frequencies of parasites carrying varied numbers of mutant mtDNA copies. A large interpatient variability in the size of this mutant reservoir was found; this variability was due to the different parameters tested but also to the relaxed replication and partitioning of mtDNA copies during mitosis. We also report seven clinical cases in which AP-resistant infections were treated by AP. These provided evidence that parasiticidal drug concentrations against AP-resistant parasites were transiently obtained within days after treatment initiation. Altogether, these results suggest that each patient carries new mtDNA mutant parasites that emerge before treatment but are killed by high starting drug concentrations. However, because the size of this mutant reservoir is highly variable from patient to patient, we propose that some patients fail to eliminate all of the mutant parasites, repeatedly producing de novo AP treatment failures.

Parallel responses of bees to Pleistocene climate change in three isolated archipelagos of the southwestern Pacific

The impacts of glacial cycles on the geographical distribution and size of populations have been explored for numerous terrestrial and marine taxa. However, most studies have focused on high latitudes, with only a few focused on the response of biota to the last glacial maximum (LGM) in equatorial regions. Here, we examine how population sizes of key bee fauna in the southwest Pacific archipelagos of Fiji, Vanuatu and Samoa have fluctuated over the Quaternary. We show that all three island faunas suffered massive population declines, roughly corresponding in time to the LGM, followed by rapid expansion post-LGM. Our data therefore suggest that Pleistocene climate change has had major impacts across a very broad tropical region. While other studies indicate widespread Holarctic effects of the LGM, our data suggest a much wider range of latitudes, extending to the tropics, where these climate change repercussions were important. As key pollinators, the inferred changes in these bee faunas may have been critical in the development of the diverse Pacific island flora. The magnitude of these responses indicates future climate change scenarios may have alarming consequences for Pacific island systems involving pollinator-dependent plant communities and agricultural crops.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

Quantitative assessment of heteroplasmy of mitochondrial genome: perspectives in diagnostics and methodological pitfalls

The role of alterations of mitochondrial DNA (mtDNA) in the development of human pathologies is not understood well. Most of mitochondrial mutations are characterized by the phenomenon of heteroplasmy which is defined as the presence of a mixture of more than one type of an organellar genome within a cell or tissue. The level of heteroplasmy varies in wide range, and the expression of disease is dependent on the percent of alleles bearing mutations, thus allowing consumption that an upper threshold level may exist beyond which the mitochondrial function collapses. Recent findings have demonstrated that some mtDNA heteroplasmic mutations are associated with widely spread chronic diseases, including atherosclerosis and cancer. Actually, each etiological mtDNA mutation has its own heteroplasmy threshold that needs to be measured. Therefore, quantitative evaluation of a mutant allele of mitochondrial genome is an obvious methodological challenge, since it may be a keystone for diagnostics of individual genetic predisposition to the disease. This review provides a comprehensive comparison of methods applicable to the measurement of heteroplasmy level of mitochondrial mutations associated with the development of pathology, in particular, in atherosclerosis and its clinical manifestations.

Variegated expression of Hsp22 transgenic reporters indicates cell-specific patterns of aging in Drosophila oenocytes

The cytoplasmic chaperone gene Hsp70 and the mitochondrial chaperone gene Hsp22 are upregulated during normal aging in Drosophila in tissue-general patterns. In addition, Hsp22 reporters are dramatically upregulated during aging in a subset of the oenocytes (liver-like cells). Hsp22 reporter expression varied dramatically between individual oenocytes and between groups of oenocytes located in adjacent body segments, and was negatively correlated with accumulation of age pigment, indicating cell-specific and cell-lineage-specific patterns of oenocyte aging. Conditional transgenic systems were used to express 88 transgenes to search for trans-regulators of the Hsp70 and Hsp22 reporters during aging. The wingless gene increased tissue-general upregulation of both Hsp70 and Hsp22 reporters. In contrast, the mitochondrial genes MnSOD and Hsp22 increased expression of Hsp22 reporters in the oenocytes and decreased accumulation of age pigment in these cells. The data suggest that cell-specific and cell lineage-specific patterns of mitochondrial malfunction contribute to oenocyte aging.

Oxidative stress is not a major contributor to somatic mitochondrial DNA mutations

The accumulation of somatic mitochondrial DNA (mtDNA) mutations is implicated in aging and common diseases of the elderly, including cancer and neurodegenerative disease. However, the mechanisms that influence the frequency of somatic mtDNA mutations are poorly understood. To develop a simple invertebrate model system to address this matter, we used the Random Mutation Capture (RMC) assay to characterize the age-dependent frequency and distribution of mtDNA mutations in the fruit fly Drosophila melanogaster. Because oxidative stress is a major suspect in the age-dependent accumulation of somatic mtDNA mutations, we also used the RMC assay to explore the influence of oxidative stress on the somatic mtDNA mutation frequency. We found that many of the features associated with mtDNA mutations in vertebrates are conserved in Drosophila, including a comparable somatic mtDNA mutation frequency (∼10(-5)), an increased frequency of mtDNA mutations with age, and a prevalence of transition mutations. Only a small fraction of the mtDNA mutations detected in young or old animals were G∶C to T∶A transversions, a signature of oxidative damage, and loss-of-function mutations in the mitochondrial superoxide dismutase, Sod2, had no detectable influence on the somatic mtDNA mutation frequency. Moreover, a loss-of-function mutation in Ogg1, which encodes a DNA repair enzyme that removes oxidatively damaged deoxyguanosine residues (8-hydroxy-2'-deoxyguanosine), did not significantly influence the somatic mtDNA mutation frequency of Sod2 mutants. Together, these findings indicate that oxidative stress is not a major cause of somatic mtDNA mutations. Our data instead suggests that somatic mtDNA mutations arise primarily from errors that occur during mtDNA replication. Further studies using Drosophila should aid in the identification of factors that influence the frequency of somatic mtDNA mutations.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Wolbachia variants induce differential protection to viruses in Drosophila melanogaster: a phenotypic and phylogenomic analysis

Wolbachia are intracellular bacterial symbionts that are able to protect various insect hosts from viral infections. This tripartite interaction was initially described in Drosophila melanogaster carrying wMel, its natural Wolbachia strain. wMel has been shown to be genetically polymorphic and there has been a recent change in variant frequencies in natural populations. We have compared the antiviral protection conferred by different wMel variants, their titres and influence on host longevity, in a genetically identical D. melanogaster host. The phenotypes cluster the variants into two groups--wMelCS-like and wMel-like. wMelCS-like variants give stronger protection against Drosophila C virus and Flock House virus, reach higher titres and often shorten the host lifespan. We have sequenced and assembled the genomes of these Wolbachia, and shown that the two phenotypic groups are two monophyletic groups. We have also analysed a virulent and over-replicating variant, wMelPop, which protects D. melanogaster even better than the closely related wMelCS. We have found that a ~21 kb region of the genome, encoding eight genes, is amplified seven times in wMelPop and may be the cause of its phenotypes. Our results indicate that the more protective wMelCS-like variants, which sometimes have a cost, were replaced by the less protective but more benign wMel-like variants. This has resulted in a recent reduction in virus resistance in D. melanogaster in natural populations worldwide. Our work helps to understand the natural variation in wMel and its evolutionary dynamics, and inform the use of Wolbachia in arthropod-borne disease control.

Monophyly of Wolbachia pipientis genomes within Drosophila melanogaster: geographic structuring, titre variation and host effects across five populations

Wolbachia pipientis is one of the most widely studied endosymbionts today, yet we know little about its short-term adaptation and evolution. Here, using a set of 91 inbred Drosophila melanogaster lines from five populations, we explore patterns of diversity and recent evolution in the Wolbachia strain wMel. Within the D. melanogaster lines, we identify six major mitochondrial clades and four wMel clades. Concordant with past studies, the Wolbachia haplotypes contain an overall low level of nucleotide diversity, yet they still display geographic structuring. Using Bayesian analysis informed with demographic estimates of colonization times, we estimate that all extant D. melanogaster mitochondrial haplotypes coalesce to a Wolbachia-infected ancestor approximately 2200 years ago. Finally, we measure wMel titre within the infected flies and find that titre varies across populations, an effect attributable to host genetic factors. This demonstration of local phenotypic divergence suggests that intraspecific host genetic variation plays a key role in shaping this model symbiotic system.

Evolutionary perspectives on the links between mitochondrial genotype and disease phenotype

BACKGROUND

Disorders of the mitochondrial respiratory chain are heterogeneous in their symptoms and underlying genetics. Simple links between candidate mutations and expression of disease phenotype typically do not exist. It thus remains unclear how the genetic variation in the mitochondrial genome contributes to the phenotypic expression of complex traits and disease phenotypes.

SCOPE OF REVIEW

I summarize the basic genetic processes known to underpin mitochondrial disease. I highlight other plausible processes, drawn from the evolutionary biological literature, whose contribution to mitochondrial disease expression remains largely empirically unexplored. I highlight recent advances to the field, and discuss common-ground and -goals shared by researchers across medical and evolutionary domains.

MAJOR CONCLUSIONS

Mitochondrial genetic variance is linked to phenotypic variance across a variety of traits (e.g. reproductive function, life expectancy) fundamental to the upkeep of good health. Evolutionary theory predicts that mitochondrial genomes are destined to accumulate male-harming (but female-friendly) mutations, and this prediction has received proof-of-principle support. Furthermore, mitochondrial effects on the phenotype are typically manifested via interactions between mitochondrial and nuclear genes. Thus, whether a mitochondrial mutation is pathogenic in effect can depend on the nuclear genotype in which is it expressed.

GENERAL SIGNIFICANCE

Many disease phenotypes associated with OXPHOS malfunction might be determined by the outcomes of mitochondrial-nuclear interactions, and by the evolutionary forces that historically shaped mitochondrial DNA (mtDNA) sequences. Concepts and results drawn from the evolutionary sciences can have broad, but currently under-utilized, applicability to the medical sciences and provide new insights into understanding the complex genetics of mitochondrial disease. This article is part of a Special Issue entitled Frontiers of Mitochondrial Research.

Introgression in the Drosophila subobscura--D. Madeirensis sister species: evidence of gene flow in nuclear genes despite mitochondrial differentiation

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.

Mitochondrial cox sequences of Nilaparvata lugens and Sogatella furcifera (Hemiptera, Delphacidae): low specificity among Asian planthopper populations

The brown planthoppers (BPH) Nilaparvata lugens (Stål) and the white-backed planthoppers (WBPH) Sogatella furcifera (Horváth) annually migrate from tropical and subtropical regions to temperate regions in Asia, including Japan, Korea and northern China. To elucidate the genetic divergence based on geography of planthoppers and to estimate their migration route on the basis of molecular data, we analysed a part of their mitochondrial genome sequences. Sequences of cytochrome oxidase subunit I (cox1) - transfer RNA for Leu (trnL2) - cox2 were determined for 579 BPH (1,928 bp) and 464 WBPH (1,927 bp) individuals collected from 31 and 25 locations, respectively, in East and Southeast Asia. Thirty and 20 mitochondrial haplotypes were detected for BPH and WBPH, respectively. Single populations of both planthoppers included multiple haplotypes, and many haplotypes were shared in some populations and areas. The most frequently detected haplotypes accounted for approximately 50% of all BPH and WBPH individuals. To evaluate gene flow among planthoppers in different regions in Asia, pairwise fixation index (Fst) values were calculated. For BPH, high Fst values (0.580-0.926) were shown between planthoppers in Papua New Guinea (PNG) and the other areas and moderate Fst values (0.176-0.362) were observed between those in southern Philippines and other areas. For WBPH, the Fst value was the highest between Taiwan and southern Vietnam (0.236), and low among the other areas. AMOVA indicated no genetic structure among eight areas, excluding southern Philippines and PNG, for BPH, and among ten areas for WBPH. These data indicate that both planthoppers do not show much differentiation of local populations and/or have genetically intermixed Asian populations. These data also indicate that it may be difficult to distinguish regional planthopper populations on the basis of differences in mitochondrial sequences.

Somatic deleterious mutation rate in a woody plant: estimation from phenotypic data

We conducted controlled crosses in populations of the long-lived clonal shrub, Vaccinium angustifolium (lowbush blueberry) to estimate inbreeding depression and mutation parameters associated with somatic deleterious mutation. Inbreeding depression level was high, with many plants failing to set fruit after self-pollination. We also compared fruit set from autogamous pollinations (pollen collected from within the same inflorescence) with fruit set from geitonogamous pollinations (pollen collected from the same plant but from inflorescences separated by several meters of branch growth). The difference between geitonogamous versus autogamous fitness within single plants is referred to as 'autogamy depression' (AD). AD can be caused by somatic deleterious mutation. AD was significantly different from zero for fruit set. We developed a maximum-likelihood procedure to estimate somatic mutation parameters from AD, and applied it to geitonogamous and autogamous fruit set data from this experiment. We infer that, on average, approximately three sublethal, partially dominant somatic mutations exist within the crowns of the plants studied. We conclude that somatic mutation in this woody plant results in an overall genomic deleterious mutation rate that exceeds the rate measured to date for annual plants. Some implications of this result for evolutionary biology and agriculture are discussed.

An Incompatibility between a mitochondrial tRNA and its nuclear-encoded tRNA synthetase compromises development and fitness in Drosophila

Mitochondrial transcription, translation, and respiration require interactions between genes encoded in two distinct genomes, generating the potential for mutations in nuclear and mitochondrial genomes to interact epistatically and cause incompatibilities that decrease fitness. Mitochondrial-nuclear epistasis for fitness has been documented within and between populations and species of diverse taxa, but rarely has the genetic or mechanistic basis of these mitochondrial–nuclear interactions been elucidated, limiting our understanding of which genes harbor variants causing mitochondrial–nuclear disruption and of the pathways and processes that are impacted by mitochondrial–nuclear coevolution. Here we identify an amino acid polymorphism in the Drosophila melanogaster nuclear-encoded mitochondrial tyrosyl–tRNA synthetase that interacts epistatically with a polymorphism in the D. simulans mitochondrial-encoded tRNA^Tyr to significantly delay development, compromise bristle formation, and decrease fecundity. The incompatible genotype specifically decreases the activities of oxidative phosphorylation complexes I, III, and IV that contain mitochondrial-encoded subunits. Combined with the identity of the interacting alleles, this pattern indicates that mitochondrial protein translation is affected by this interaction. Our findings suggest that interactions between mitochondrial tRNAs and their nuclear-encoded tRNA synthetases may be targets of compensatory molecular evolution. Human mitochondrial diseases are often genetically complex and variable in penetrance, and the mitochondrial–nuclear interaction we document provides a plausible mechanism to explain this complexity.

The ancient symbiosis between two prokaryotes that gave rise to the eukaryotic cell has required genomic cooperation for at least a billion years. Eukaryotic cells respire through the coordinated expression of their nuclear and mitochondrial genomes, both of which encode the proteins and RNAs required for mitochondrial transcription, translation, and aerobic respiration. Genetic interactions between these genomes are hypothesized to influence the effects of mitochondrial mutations on disease and drive mitochondrial–nuclear coevolution. Here we characterize the molecular cause and the cellular and organismal consequences of a mitochondrial–nuclear interaction in Drosophila between naturally occurring mutations in a mitochondrial tRNA and a nuclear-encoded tRNA synthetase. These mutations have little effect on their own; but, when combined, they severely compromise development and reproduction. tRNA synthetases attach the appropriate amino acid onto their cognate tRNA, and this reaction is required for efficient and accurate protein synthesis. We show that disruption of this interaction compromises mitochondrial function, providing hypotheses for the variable penetrance of diseases associated with mitochondrial tRNAs and for which pathways and processes are likely to be affected by mitochondrial–nuclear interactions.

Coevolution of Drosophila melanogaster mtDNA and Wolbachia genotypes

Maternally inherited microorganisms can influence the mtDNA pattern of variation in hosts. This influence is driven by selection among symbionts and can cause the frequency of mitochondrial variants in the population to eventually increase or decrease. Wolbachia infection is common and widespread in Drosophila melanogaster populations. We compared genetic variability of D. melanogaster mitotypes with Wolbachia genotypes among isofemale lines associated with different geographic locations and time intervals to study coevolution of the mtDNA and Wolbachia. Phylogenetic analysis of D. melanogaster mtDNA revealed two clades diverged in Africa, each associated with one of the two Wolbachia genotype groups. No evidence of horizontal transmission of Wolbachia between maternal lineages has been found. All the mtDNA variants that occur in infected isofemale lines are found in uninfected isofemale lines and vice versa, which is indicative of a recent loss of infection from some maternal fly lineages and confirms a significant role of Wolbachia in the D. melanogaster mtDNA pattern of variation. Finally, we present a comparative analysis of biogeographic distribution of D. melanogaster mitotypes all over the world.

Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster

Wolbachia are maternally inherited symbiotic bacteria, commonly found in arthropods, which are able to manipulate the reproduction of their host in order to maximise their transmission. The evolutionary history of endosymbionts like Wolbachia can be revealed by integrating information on infection status in natural populations with patterns of sequence variation in Wolbachia and host mitochondrial genomes. Here we use whole-genome resequencing data from 290 lines of Drosophila melanogaster from North America, Europe, and Africa to predict Wolbachia infection status, estimate relative cytoplasmic genome copy number, and reconstruct Wolbachia and mitochondrial genome sequences. Overall, 63% of Drosophila strains were predicted to be infected with Wolbachia by our in silico analysis pipeline, which shows 99% concordance with infection status determined by diagnostic PCR. Complete Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with strict vertical transmission through the maternal cytoplasm and imperfect transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most recent common ancestor of all Wolbachia and mitochondrial genomes in D. melanogaster dates to around 8,000 years ago. We find evidence for a recent global replacement of ancestral Wolbachia and mtDNA lineages, but our data suggest that the derived wMel lineage arose several thousand years ago, not in the 20th century as previously proposed. Our data also provide evidence that this global replacement event is incomplete and is likely to be one of several similar incomplete replacement events that have occurred since the out-of-Africa migration that allowed D. melanogaster to colonize worldwide habitats. This study provides a complete genomic analysis of the evolutionary mode and temporal dynamics of the D. melanogaster–Wolbachia symbiosis, as well as important resources for further analyses of the impact of Wolbachia on host biology.

Host–microbe interactions play important roles in the physiology, development, and ecology of many organisms. Studying how hosts and their microbial symbionts evolve together over time is crucial for understanding the impact that microbes have on host biology. With the advent of high-throughput sequencing technologies, it is now possible to obtain complete genomic information for hosts and their associated microbes. Here we use whole-genome sequences from ∼300 strains of the fruitfly Drosophila melanogaster to reveal the evolutionary history of this model species and its intracellular bacterial symbiont Wolbachia. The major findings of this study are that Wolbachia in D. melanogaster is inherited strictly through the egg with no evidence of horizontal transfer from other species, that the genealogies of Wolbachia and mitochondrial genomes are virtually the same, and that both Wolbachia and mitochondrial genomes show evidence for a recent incomplete global replacement event, which has left remnant lineages in North America, Europe, and Africa. We also use the fact that Wolbachia and mitochondrial genomes have the same genealogy to estimate the rate of molecular evolution for Wolbachia, which allows us to put dates on key events in the history of this important host–microbe model system.

Properties of spontaneous mutational variance and covariance for wing size and shape in Drosophila melanogaster

We estimated mutational variance-covariance matrices, M, for wing shape and size in two genotypes of Drosophila melanogaster after 192 generations of mutation accumulation. We characterized 21 potentially independent aspects of wing shape and size using geometric morphometrics, and analyzed the data using a likelihood-based factor-analytic approach. We implement a previously unused analysis that describes those directions with the greatest difference in evolvability between pairs of matrices. There are significant mutational effects on 19 of 21 possible aspects of wing form, consistent with the high dimensionality of standing genetic variation for wing shape previously identified in D. melanogaster. Mutations have partially recessive effects, consistent with average dominance around 0.25. Sex-specific matrices are relatively similar, although male-specific matrices are slightly larger, as expected due to dosage compensation on the X chromosome. Genotype-specific matrices are quite different. Matrices may differ both because of sampling error based on small samples of mutations with large phenotypic effects, and because of the mutational properties of the genotypes. Genotypic differences are likely to be involved, as the two genotypes have different molecular mutation rates and properties.

Extraordinary genome stability in the ciliate Paramecium tetraurelia

Mutation plays a central role in all evolutionary processes and is also the basis of genetic disorders. Established base-substitution mutation rates in eukaryotes range between ∼5 × 10(-10) and 5 × 10(-8) per site per generation, but here we report a genome-wide estimate for Paramecium tetraurelia that is more than an order of magnitude lower than any previous eukaryotic estimate. Nevertheless, when the mutation rate per cell division is extrapolated to the length of the sexual cycle for this protist, the measure obtained is comparable to that for multicellular species with similar genome sizes. Because Paramecium has a transcriptionally silent germ-line nucleus, these results are consistent with the hypothesis that natural selection operates on the cumulative germ-line replication fidelity per episode of somatic gene expression, with the germ-line mutation rate per cell division evolving downward to the lower barrier imposed by random genetic drift. We observe ciliate-specific modifications of widely conserved amino acid sites in DNA polymerases as one potential explanation for unusually high levels of replication fidelity.

Male-killing Wolbachia and mitochondrial selective sweep in a migratory African insect

Background

Numerous recent studies have shown that resident symbiotic microorganisms of insects play a fundamental role in host ecology and evolution. The lepidopteran pest, African armyworm (Spodoptera exempta), is a highly migratory and destructive species found throughout sub-Saharan Africa, that can experience eruptive outbreaks within the space of a single generation, making predicting population dynamics and pest control forecasting extremely difficult. Three strains of Wolbachia have recently been identified infecting this species in populations sampled from Tanzania. In this study, we examined the interaction between Wolbachia pipiensis infections and the co-inherited marker, mtDNA, within populations of armyworm, as a means to investigate the population biology and evolutionary history of Wolbachia and its host.

Results

A Wolbachia-infected isofemale line was established in the laboratory. Phenotypic studies confirmed the strain wExe1 as a male-killer. Partial sequencing of the mitochondrial COI gene from 164 individual field-collected armyworm of known infection status revealed 17 different haplotypes. There was a strong association between Wolbachia infection status and mtDNA haplotype, with a single dominant haplotype, haplo1 (90.2% prevalence), harbouring the endosymbiont. All three Wolbachia strains were associated with this haplotype. This indicates that Wolbachia may be driving a selective sweep on armyworm haplotype diversity. Despite very strong biological and molecular evidence that the samples represent a single species (including from nuclear 28S gene markers), the 17 haplotypes did not fall into a monophyletic clade within the Spodoptera genus; with six haplotypes (2 each from 3 geographically separate populations) differing by >11% in their nucleotide sequence to the other eleven.

Conclusions

This study suggests that three strains of Wolbachia may be driving a selective sweep on armyworm haplotype diversity, and that based on COI sequence data, S. exempta is not a monophyletic group within the Spodoptera genus. This has clear implications for the use of mtDNA as neutral genetic markers in insects, and also demonstrates the impact of Wolbachia infections on host evolutionary genetics.

Factors affecting the relative abundance of nuclear copies of mitochondrial DNA (numts) in hominoids

Although nuclear copies of mitochondrial DNA (numts) can originate from any portion of the mitochondrial genome, evidence from humans suggests that more variable parts of the mitochondrial genome, such as the mitochondrial control region (MCR), are under-represented in the nucleus. This apparent deficit might arise from the erosion of sequence identity in numts originating from rapidly evolving mitochondrial sequences. However, the extent to which mitochondrial sequence properties impacts the number of numts detected in genomic surveys has not been evaluated. In order to address this question, we: (1) conducted exhaustive BLAST searches of MCR numts in three hominoid genomes; (2) assessed numt prevalence across the four MCR sub-domains (HV1, CCD, HV2, and MCR(F)); (3) estimated their insertion rates in great apes (Hominoidea); and (4) examined the relationship between mitochondrial DNA variability and numt prevalence in sequences originating from MCR and coding regions of the mitochondrial genome. Results indicate a marked deficit of numts from HV2 and MCR(F) MCR sub-domains in all three species. These MCR sub-domains exhibited the highest proportion of variable sites and the lowest number of detected numts per mitochondrial site. Variation in MCR insertion rate between lineages was also observed with a pronounced burst in recent integrations within chimpanzees and orangutans. A deficit of numts from HV2/MCR(F) was observed regardless of age, whereas HV1 is under-represented only in older numts (>25 million years). Finally, more variable mitochondrial genes also exhibit a lower identity with nuclear copies and because of this, appear to be under-represented in human numt databases.

Whole genome sequencing of mutation accumulation lines reveals a low mutation rate in the social amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10(-9), with a Poisson confidence interval of 4.1×10(-9) - 9.5×10(-9), per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10(-11), with a Poisson confidence interval ranging from 7.4×10(-13) to 1.6×10(-10), is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.

Mitochondrial genome evolution in a single protoploid yeast species

Mitochondria are organelles, which play a key role in some essential functions, including respiration, metabolite biosynthesis, ion homeostasis, and apoptosis. The vast numbers of mitochondrial DNA (mtDNA) sequences of various yeast species, which have recently been published, have also helped to elucidate the structural diversity of these genomes. Although a large corpus of data are now available on the diversity of yeast species, little is known so far about the mtDNA diversity in single yeast species. To study the genetic variations occurring in the mtDNA of wild yeast isolates, we performed a genome-wide polymorphism survey on the mtDNA of 18 Lachancea kluyveri (formerly Saccharomyces kluyveri) strains. We determined the complete mt genome sequences of strains isolated from various geographical locations (in North America, Asia, and Europe) and ecological niches (Drosophila, tree exudates, soil). The mt genome of the NCYC 543 reference strain is 51,525 bp long. It contains the same core of genes as Lachancea thermotolerans, the nearest relative to L. kluyveri. To explore the mt genome variations in a single yeast species, we compared the mtDNAs of the 18 isolates. The phylogeny and population structure of L. kluyveri provide clear-cut evidence for the existence of well-defined geographically isolated lineages. Although these genomes are completely syntenic, their size and the intron content were found to vary among the isolates studied. These genomes are highly polymorphic, showing an average diversity of 28.5 SNPs/kb and 6.6 indels/kb. Analysis of the SNP and indel patterns showed the existence of a particularly high overall level of polymorphism in the intergenic regions. The dN/dS ratios obtained are consistent with purifying selection in all these genes, with the noteworthy exception of the VAR1 gene, which gave a very high ratio. These data suggest that the intergenic regions have evolved very fast in yeast mitochondrial genomes.

Detection of ultra-rare mutations by next-generation sequencing

Next-generation DNA sequencing promises to revolutionize clinical medicine and basic research. However, while this technology has the capacity to generate hundreds of billions of nucleotides of DNA sequence in a single experiment, the error rate of ~1% results in hundreds of millions of sequencing mistakes. These scattered errors can be tolerated in some applications but become extremely problematic when "deep sequencing" genetically heterogeneous mixtures, such as tumors or mixed microbial populations. To overcome limitations in sequencing accuracy, we have developed a method termed Duplex Sequencing. This approach greatly reduces errors by independently tagging and sequencing each of the two strands of a DNA duplex. As the two strands are complementary, true mutations are found at the same position in both strands. In contrast, PCR or sequencing errors result in mutations in only one strand and can thus be discounted as technical error. We determine that Duplex Sequencing has a theoretical background error rate of less than one artifactual mutation per billion nucleotides sequenced. In addition, we establish that detection of mutations present in only one of the two strands of duplex DNA can be used to identify sites of DNA damage. We apply the method to directly assess the frequency and pattern of random mutations in mitochondrial DNA from human cells.

Phylogenetic analysis of four nuclear protein-encoding genes largely corroborates the traditional classification of Bivalvia (Mollusca)

Revived interest in molluscan phylogeny has resulted in a torrent of molecular sequence data from phylogenetic, mitogenomic, and phylogenomic studies. Despite recent progress, basal relationships of the class Bivalvia remain contentious, owing to conflicting morphological and molecular hypotheses. Marked incongruity of phylogenetic signal in datasets heavily represented by nuclear ribosomal genes versus mitochondrial genes has also impeded consensus on the type of molecular data best suited for investigating bivalve relationships. To arbitrate conflicting phylogenetic hypotheses, we evaluated the utility of four nuclear protein-encoding genes-ATP synthase β, elongation factor-1α, myosin heavy chain type II, and RNA polymerase II-for resolving the basal relationships of Bivalvia. We sampled all five major lineages of bivalves (Archiheterodonta, Euheterodonta [including Anomalodesmata], Palaeoheterodonta, Protobranchia, and Pteriomorphia) and inferred relationships using maximum likelihood and Bayesian approaches. To investigate the robustness of the phylogenetic signal embedded in the data, we implemented additional datasets wherein length variability and/or third codon positions were eliminated. Results obtained include (a) the clade (Nuculanida+Opponobranchia), i.e., the traditionally defined Protobranchia; (b) the monophyly of Pteriomorphia; (c) the clade (Archiheterodonta+Palaeoheterodonta); (d) the monophyly of the traditionally defined Euheterodonta (including Anomalodesmata); and (e) the monophyly of Heteroconchia, i.e., (Palaeoheterodonta+Archiheterodonta+Euheterodonta). The stability of the basal tree topology to dataset manipulation is indicative of signal robustness in these four genes. The inferred tree topology corresponds closely to those obtained by datasets dominated by nuclear ribosomal genes (18S rRNA and 28S rRNA), controverting recent taxonomic actions based solely upon mitochondrial gene phylogenies.

The relation between the neutrality index for mitochondrial genes and the distribution of mutational effects on fitness

We explore factors affecting patterns of polymorphism and divergence (as captured by the neutrality index) at mammalian mitochondrial loci. To do this, we develop a population genetic model that incorporates a fraction of neutral amino acid sites, mutational bias, and a probability distribution of selection coefficients against new nonsynonymous mutations. We confirm, by reanalyzing publicly available datasets, that the mitochondrial cyt-b gene shows a broad range of neutrality indices across mammalian taxa, and explore the biological factors that can explain this observation. We find that observed patterns of differences in the neutrality index, polymorphism, and divergence are not caused by differences in mutational bias. They can, however, be explained by a combination of a small fraction of neutral amino acid sites, weak selection acting on most amino acid mutations, and differences in effective population size among taxa.

How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.

Fitness effects of mutations in bacteria

Mutation is the primary source of variation in any organism. Without it, natural selection cannot operate and organisms cannot adapt to novel environments. Mutation is also generally a source of defect: many mutations are not neutral but cause fitness decreases in the organisms where they arise. In bacteria, another important source of variation is horizontal gene transfer. This source of variation can also cause beneficial or deleterious effects. Determining the distribution of fitness effects of mutations in different environments and genetic backgrounds is an active research field. In bacteria, knowledge of these distributions is key for understanding important traits. For example, for determining the dynamics of microorganisms with a high genomic mutation rate (mutators), and for understanding the evolution of antibiotic resistance, and the emergence of pathogenic traits. All of these characteristics are extremely relevant for human health both at the individual and population levels. Experimental evolution has been a valuable tool to address these questions. Here, we review some of the important findings of mutation effects in bacteria revealed through laboratory experiments.