Identification of cytoplasmically transferred mitochondrial DNA in female germlines of Drosophila and its propagation in the progeny

Summary The composition of mitochondrial DNA (mtDNA) was analyzed in single female flies that developed from fertilized Drosophila melanogaster eggs, into which germ plasm of D. simulans had been introduced. HpaII cleavage patterns showed that all 12 individual female flies examined had developed from eggs in which 37%–71% of the total mtDNA was D. simulans mtDNA (Ds mtDNA) and the rest was D. melanogaster mtDNA (Dm mtDNA). The stability of this heteroplasmic state in these isofemale lines was monitored for seven generations at both individual and population levels. Results showed that the heteroplasmy of Dm and Ds mtDNAs was stably transmitted for at least three generations at the population level, but showed stochastic segregation at the individual level. After 4–6 generations, all individuals lost Ds mtDNA. The mechanisms of preferential loss of Ds mtDNA and of transmission of heteroplasmic mtDNA to descendants are discussed.     

The on‐again,off‐again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co‐evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.     

Molecular Population Genetics of Mtdna Size Variation in Crickets

Nucleotide sequence analysis of a region of cricket (Gryllus firmus) mtDNA showing discrete length variation revealed tandemly repeated sequences 220 base pairs (bp) in length. The repeats consist of 206 bp sequences bounded by the dyad symmetric sequence 5'GGGGGCATGCCCCC3'. The sequence data showed that mtDNA size variation in this species is due to variation in the number of copies of tandem repeats. Southern blot analysis was used to document the frequency of crickets heteroplasmic for two or more different-sized mtDNAs. In New England populations of G. firmus and a close relative Gryllus pennsylvanicus approximately 60% of the former and 45% of the latter were heteroplasmic. From densitometry of autoradiographs the frequencies of mtDNA size classes were determined for the population samples and are shown to very different in the two species. However, in populations where hybridization between the two species has occurred, the frequencies of size classes and cytoplasmic genotypes in each species' distinct mtDNA lineage were shifted in a manner suggesting nuclear-cytoplasmic interactions. The data were applied to reported diversity indices and hierarchical statistics. The hierarchical statistics indicated that the greatest proportion of variation for mtDNA size was due to variation among individuals in their cytoplasmic genotypes (heteroplasmic or homoplasmic state). The diversity indices were used to estimate a per-generation mutation rate for size variants of 10(-4). The data are discussed in light of the relationship between genetic drift and mutation in maintaining variation for mtDNA size.     

Heteroplasmy of Short Tandem Repeats in Mitochondrial DNA of Atlantic Cod, Gadus Morhua

The mitochondrial DNA of the Atlantic cod (Gadus morhua) contains a tandem array of 40-bp repeats in the D-loop region of the molecule. Variation among molecules in the copy number of these repeats results in mtDNA length variation and heteroplasmy (the presence of more than one form of mtDNA in an individual). In a sample of fish collected from different localities around Iceland and off George's Bank, each individual was heteroplasmic for two or more mtDNAs ranging in repeat copy number from two (common) to six (rare). An earlier report on mtDNA heteroplasmy in sturgeon (Acipenser transmontanus) presented a competitive displacement model for length mutations in mtDNAs containing tandem arrays and the cod data deviate from this model. Depending on the nature of putative secondary structures and the location of D-loop strand termination, additional mechanisms of length mutation may be needed to explain the range of mtDNA length variants maintained in these populations. The balance between genetic drift and mutation in maintaining this length polymorphism is estimated through a hierarchical analysis of diversity of mtDNA length variation in the Iceland samples. Eighty percent of the diversity lies within individuals, 8% among individuals and 12% among localities. An estimate of theta = 2N(eo) mu greater than 1 indicates that this system is characterized by a high mutation rate and is governed primarily by deterministic dynamics. The sequences of repeat arrays from fish collected in Norway, Iceland and George's Bank show no nucleotide variation suggesting that there is very little substructuring to the North Atlantic cod population.     

Replicative advantage and tissue-specific segregation of RR mitochondrial DNA between C57BL/6 and RR heteroplasmic mice

To investigate the interactions between mtDNA and nuclear genomes, we produced heteroplasmic maternal lineages by transferring the cytoplasts between the embryos of two mouse strains, C57BL/6 (B6) and RR. A total of 43 different nucleotides exist in the displacement-loop (D-loop) region of mtDNA between B6 and RR. Heteroplasmic embryos were reconstructed by electrofusion using a blastomere from a two-cell stage embryo of one strain and an enucleated blastomere from a two-cell stage embryo of the other strain. Equivalent volumes of both types of mtDNAs were detected in blastocyst stage embryos. However, the mtDNA from the RR strain became biased in the progeny, regardless of the source of the nuclear genome. The RR mtDNA population was very high in most of the tissues examined but was relatively low in the brain and the heart. An age-related increase of RR mtDNA was also observed in the blood. The RR mtDNAs in the reconstructed embryos and in the embryos collected from heteroplasmic mice showed a different segregation pattern during early embryonic development. These results suggest that the RR mtDNA has a replicative advantage over B6 mtDNA during embryonic development and differentiation, regardless of the type of nuclear genome.     

Non-Random mtDNA Segregation Patterns Indicate a Metastable Heteroplasmic Segregation Unit in m.3243A>G Cybrid Cells

Many pathogenic mitochondrial DNA mutations are heteroplasmic, with a mixture of mutated and wild-type mtDNA present within individual cells. The severity and extent of the clinical phenotype is largely due to the distribution of mutated molecules between cells in different tissues, but mechanisms underpinning segregation are not fully understood. To facilitate mtDNA segregation studies we developed assays that measure m.3243A>G point mutation loads directly in hundreds of individual cells to determine the mechanisms of segregation over time. In the first study of this size, we observed a number of discrete shifts in cellular heteroplasmy between periods of stable heteroplasmy. The observed patterns could not be parsimoniously explained by random mitotic drift of individual mtDNAs. Instead, a genetically metastable, heteroplasmic mtDNA segregation unit provides the likely explanation, where stable heteroplasmy is maintained through the faithful replication of segregating units with a fixed wild-type/m.3243A>G mutant ratio, and shifts occur through the temporary disruption and re-organization of the segregation units. While the nature of the physical equivalent of the segregation unit remains uncertain, the factors regulating its organization are of major importance for the pathogenesis of mtDNA diseases.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

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

Frequency and Pattern of Heteroplasmy in the Complete Human Mitochondrial Genome

Determining the levels of human mitochondrial heteroplasmy is of utmost importance in several fields. In spite of this, there are currently few published works that have focused on this issue. In order to increase the knowledge of mitochondrial DNA (mtDNA) heteroplasmy, the main goal of this work is to investigate the frequency and the mutational spectrum of heteroplasmy in the human mtDNA genome. To address this, a set of nine primer pairs designed to avoid co-amplification of nuclear DNA (nDNA) sequences of mitochondrial origin (NUMTs) was used to amplify the mitochondrial genome in 101 individuals. The analysed individuals represent a collection with a balanced representation of genders and mtDNA haplogroup distribution, similar to that of a Western European population. The results show that the frequency of heteroplasmic individuals exceeds 61%. The frequency of point heteroplasmy is 28.7%, with a widespread distribution across the entire mtDNA. In addition, an excess of transitions in heteroplasmy were detected, suggesting that genetic drift and/or selection may be acting to reduce its frequency at population level. In fact, heteroplasmy at highly stable positions might have a greater impact on the viability of mitochondria, suggesting that purifying selection must be operating to prevent their fixation within individuals. This study analyses the frequency of heteroplasmy in a healthy population, carrying out an evolutionary analysis of the detected changes and providing a new perspective with important consequences in medical, evolutionary and forensic fields.     

Mitochondrial DNA Transmission Genetics in Crickets

This paper presents the results of a single generation study of the transmission genetics of mitochondrial DNA in the field cricket Gryllus firmus. In this species, individuals heteroplasmic for at least two different-sized mitochondrial genomes can be collected easily from natural populations. The frequencies of mtDNA size variants in heteroplasmic females and samples of their offspring were estimated by densitometry of autoradiographs. The variance in mitochondrial genotype frequencies among the offspring of heteroplasmic females indicates that, through genetic drift, fixation would take several hundred animal generations. Differences between the observations and data on mtDNA transmission in yeast and cows are discussed in light of the differences in organelle sampling regime and early developmental events in these species. Our data also show shifts in genotype frequencies in the transmission from mother to offspring that suggest a bias in favor of smaller genomes. The nature of mtDNA size variation in natural populations of crickets is discussed in reference to a mutation-selection balance.     

Selfish Little Circles: Transmission Bias and Evolution of Large Deletion-Bearing Mitochondrial DNA in Caenorhabditis briggsae Nematodes

Selfish DNA poses a significant challenge to genome stability and organismal fitness in diverse eukaryotic lineages. Although selfish mitochondrial DNA (mtDNA) has known associations with cytoplasmic male sterility in numerous gynodioecious plant species and is manifested as petite mutants in experimental yeast lab populations, examples of selfish mtDNA in animals are less common. We analyzed the inheritance and evolution of mitochondrial DNA bearing large heteroplasmic deletions including nad5 gene sequences (nad5Δ mtDNA), in the nematode Caenorhabditis briggsae. The deletion is widespread in C. briggsae natural populations and is associated with deleterious organismal effects. We studied the inheritance patterns of nad5Δ mtDNA using eight sets of C. briggsae mutation-accumulation (MA) lines, each initiated from a different natural strain progenitor and bottlenecked as single hermaphrodites across generations. We observed a consistent and strong drive toward higher levels of deletion-bearing molecules in the heteroplasmic pool of mtDNA after ten generations of bottlenecking. Our results demonstrate a uniform transmission bias whereby nad5Δ mtDNA accumulates to higher levels relative to intact mtDNA in multiple genetically diverse natural strains of C. briggsae. We calculated an average 1% per-generation transmission bias for deletion-bearing mtDNA relative to intact genomes. Our study, coupled with known deleterious phenotypes associated with high deletion levels, shows that nad5Δ mtDNA are selfish genetic elements that have evolved in natural populations of C. briggsae, offering a powerful new system to study selfish mtDNA dynamics in metazoans.     

Detection of mutations in mitochondrial DNA by droplet digital PCR

Mutations in mitochondrial DNA (mtDNA) may result in various pathological processes. Detection of mutant mtDNAs is a problem for diagnostic practice that is complicated by heteroplasmy – a phenomenon of the inferring presence of at least two allelic variants of the mitochondrial genome. Also, the level of heteroplasmy largely determines the profile and severity of clinical manifestations. Here we discuss detection of mutations in heteroplasmic mtDNA using up-todate methods that have not yet been introduced as routine clinical assays. These methods can be used for detecting mutations in mtDNA to verify diagnosis of “mitochondrial disease”, studying dynamics of mutant mtDNA in body tissues of patients, as well as investigating structural features of mtDNAs. Original data on allele-specific discrimination of m.11778G>A mutation by droplet digital PCR are presented, which demonstrate an opportunity for simultaneous detection and quantitative assessment of mutations in mtDNAs.     

Distribution and threshold expression of the tRNA(Lys) mutation in skeletal muscle of patients with myoclonic epilepsy and ragged-red fibers (MERRF).

We investigated the distribution and expression of mutant mtDNAs carrying the A-to-G mutation at position 8344 in the tRNA(Lys) gene in the skeletal muscle of four patients with myoclonus epilepsy and ragged-red fibers (MERRF). The proportion of mutant genomes was greater than 80% of total mtDNAs in muscle samples of all patients and was associated with a decrease in the activity of cytochrome c oxidase (COX). The vast majority of myoblasts, cloned from the satellite-cell population in the same muscles, were homoplasmic for the mutation. The overall proportion of mutant mtDNAs in this population was similar to that in differentiated muscle, suggesting that the ratio of mutant to wild-type mtDNAs in skeletal muscle is determined either in the ovum or during early development and changes little with age. Translation of all mtDNA-encoded genes was severely depressed in homoplasmic mutant myoblast clones but not in heteroplasmic or wild-type clones. The threshold for biochemical expression of the mutation was determined in heteroplasmic myotubes formed by fusion of different proportions of mutant and wild-type myoblasts. The magnitude of the decrease in translation in myotubes containing mutant mtDNAs was protein specific. Complex I and IV subunits were more affected than complex V subunits, and there was a rough correlation with both protein size and number of lysine residues. Approximately 15% wild-type mtDNAs restored translation and COX activity to near normal levels. These results show that the A-to-G substitution in tRNA(Lys) is a functionally recessive mutation that can be rescued by intraorganellar complementation with a small proportion of wild-type mtDNAs and explain the steep threshold for expression of the MERRF clinical phenotype.     

Multiple Origins of a Mitochondrial Mutation Conferring Deafness

A point mutation (1555G) in the smaller ribosomal subunit of the mitochondrial DNA (mtDNA) has been associated with maternally inherited traits of hypersensitivity to streptomycin and sensorineural deafness in a number of families from China, Japan, Israel, and Africa. To determine whether this distribution was the result of a single or multiple mutational events, we carried out genetic distance analysis and phylogenetic analysis of 10 independent mtDNA D-loop sequences from Africa and Asia. The mtDNA sequence diversity was high (2.21%). Phylogenetic analysis assigned 1555G-bearing haplotypes at very divergent points in the human mtDNA evolutionary tree, and the 1555G mutations occur in many cases on race-specific mtDNA haplotypes, both facts are inconsistent with a recent introgression of the mutation into these races. The simplest interpretation of the available data is that there have been multiple origins of the 1555G mutation. The genetic distance among mtDNAs bearing the pathogenic 1555G mutation is much larger than among mtDNAs bearing either evolutionarily neutral or weakly deleterious nucleotide substitutions (such as the 4336G mutation). These results are consistent with the view that pathogenic mtDNA haplotypes such as 1555G arise on disparate mtDNA lineages which because of negative natural selection leave relatively few related descendants. The co-existence of the same mutation with deafness in individuals with very different nuclear and mitochondrial genetic backgrounds confirms the pathogenicity of the 1555G mutation.     

Selective sweeps of mitochondrial DNA can drive the evolution of uniparental inheritance

Although the uniparental (or maternal) inheritance of mitochondrial DNA (mtDNA) is widespread, the reasons for its evolution remain unclear. Two main hypotheses have been proposed: selection against individuals containing different mtDNAs (heteroplasmy) and selection against “selfish” mtDNA mutations. Recently, uniparental inheritance was shown to promote adaptive evolution in mtDNA, potentially providing a third hypothesis for its evolution. Here, we explore this hypothesis theoretically and ask if the accumulation of beneficial mutations provides a sufficient fitness advantage for uniparental inheritance to invade a population in which mtDNA is inherited biparentally. In a deterministic model, uniparental inheritance increases in frequency but cannot replace biparental inheritance if only a single beneficial mtDNA mutation sweeps through the population. When we allow successive selective sweeps of mtDNA, however, uniparental inheritance can replace biparental inheritance. Using a stochastic model, we show that a combination of selection and drift facilitates the fixation of uniparental inheritance (compared to a neutral trait) when there is only a single selective mtDNA sweep. When we consider multiple mtDNA sweeps in a stochastic model, uniparental inheritance becomes even more likely to replace biparental inheritance. Our findings thus suggest that selective sweeps of beneficial mtDNA haplotypes can drive the evolution of uniparental inheritance.     

Characteristics and Distribution of Large Tandem Duplications in Brook Stickleback (Culaea Inconstans) Mitochondrial DNA

Most animal mitochondrial DNAs (mtDNAs) range in size from 15 to 18 kb, but increased sizes up to ~40 kb are occasionally found. We investigated large size variation in mtDNA of the brook stickleback fish, Culaea inconstans, and characterized four large (2.7-5.8 kb) tandem duplications. Duplications differ in size, frequency of occurrence, and degree of associated heteroplasmy, but each includes the control region and one or more adjacent genes. Duplications are correlated with two mtDNA lineages sampled from 31 populations. L(1) duplications (3.2-4.8 kb) were present in all lineage I individuals (n = 121, 19 populations); 53 fish were heteroplasmic due to variation in the copy number of a tandemly repeated 270-bp sequence within the duplicated region. In contrast, duplications L(2), L(3), and L(4) (2.7-5.8 kb) occurred in only 117 of 174 lineage II fish, in eight of 14 populations. Nine fish with L(3) or L(4) duplications were heteroplasmic, possessing some mtDNAs that lacked duplications (normal-length mtDNAs). Heteroplasmy in L(2) was associated with a small variable region near the ND5 gene. Phylogenetic analysis of restriction sites in Culaea mtDNAs and haplotype-defining sequence differences present in both copies argue for multiple independent events that gave rise to three of the four duplications.     

Evidence for Non-Neutrality of Mitochondrial DNA Haplotypes in Drosophila Pseudoobscura

Mitochondrial DNA (mtDNA) haplotypes usually are assumed to be neutral, unselected markers of evolving female lineages. This assumption was tested by monitoring haplotype frequencies in 12 experimental populations of Drosophila pseudoobscura which were polymorphic for mtDNA haplotypes. Populations were maintained for at least 10 generations, and in one case for 32 generations, while tests of mtDNA selective neutrality were conducted. In an initial population, formed from a mixture of two strains with different mitochondrial haplotypes, the frequency of the Bogota haplotype increased 46% in 3 generations, reaching an apparent equilibrium frequency of 82% after 32 generations. Perturbation of this equilibrium by addition of the less common haplotype resulted in a rapid, dramatic increase in frequency of the second haplotype, and a return to essentially the same equilibrium frequency as before perturbation. This behavior is not consistent with mtDNA neutrality, nor is the equilibrium consistent with a simple model of constant selection on the haploid mtDNAs. Replicate cage experiments with mtDNA haplotypes did not always generate the same result as the initial cage. Several lines of evidence, including manipulations of the nuclear genome, support the idea that both nuclear and mitochondrial genomes are involved in the dramatic mtDNA frequency changes. In another experiment, strong female viability selection was implicated via mtDNA frequency changes. Although the causes of the dramatic mtDNA frequency changes in our populations are not obvious, it is clear that Drosophila mitochondrial haplotypes are not always simply neutral markers. Our findings are relevant to the introduction of a novel mtDNA variant from one species or one population into another. Such introductions could be strongly favored by selection, even if it is sporadic.     

Microsatellite Evolution in the Mitochondrial Genome of Bechstein’s Bat (Myotis bechsteinii)

Being highly polymorphic, microsatellites are widely used genetic markers. They are abundant throughout the nuclear genomes of eukaryotes but rare in the mitochondrial genomes (mtDNA) of animals. We describe a short but highly polymorphic AT microsatellite in the mtDNA control region of Bechstein’s bat and discuss the role of mutation, genetic drift, and selection in maintaining its variability. As heteroplasmy and hence mutation rate were positively correlated with repeat number, a simple mutation model cannot explain the observed frequency distribution of AT copy numbers. Because of the unimodal distribution of repeat numbers found in heteroplasmic individuals, single step mutations are likely to be the predominant mechanism of copy number alternations. Above a certain copy number (seven repeats), deletions of single dinucleotide repeats seem to be more common than additions, which results in a decrease in frequency of long alleles. Heteroplasmy was inherited from mothers to their offspring and no evidence of paternal inheritance of mitochondria was found. Genetic differences accumulated with more distant ancestry, which suggests that microsatellites can be useful genetic markers in population genetics.[Reviewing Editor: Dr. Rafael Zardoya]     

mtDNA and the origin of Caucasians: identification of ancient Caucasian-specific haplogroups, one of which is prone to a recurrent somatic duplication in the D-loop region.

mtDNA sequence variation was examined in 175 Caucasians from the United States and Canada by PCR amplification and high-resolution restriction-endonuclease analysis. The majority of the Caucasian mtDNAs were subsumed within four mtDNA lineages (haplogroups) defined by mutations that are rarely seen in Africans and Mongoloids. The sequence divergence of these haplogroups indicates that they arose early in Caucasian radiation and gave raise to modern European mtDNAs. Although ancient, none of these haplogroups is old enough to be compatible with a Neanderthal origin, suggesting that Homo sapiens sapiens displaced H. s. neanderthaliensis, rather than mixed with it. The mtDNAs of one of these haplogroups have a unique homoplasmic insertion between nucleotide pair (np) 573 and np 574, within the D-loop control region. This insertion makes these mtDNAs prone to a somatic mutation that duplicates a 270-bp portion of the D-loop region between np 309 and np 572. This finding suggests that certain nonpathogenic mtDNA mutations could predispose individuals to mtDNA rearrangements.     

Mitochondria,maternal inheritance,and asymmetric fitness: Why males die younger

Mitochondrial function is achieved through the cooperative interaction of two genomes: one nuclear (nuDNA) and the other mitochondrial (mtDNA). The unusual transmission of mtDNA, predominantly maternal without recombination is predicted to affect the fitness of male offspring. Recent research suggests the strong sexual dimorphism in aging is one such fitness consequence. The uniparental inheritance of mtDNA results in a selection asymmetry; mutations that affect only males will not respond to natural selection, imposing a male‐specific mitochondrial mutation load. Prior work has implicated this male‐specific mutation load in disease and infertility, but new data from fruit flies suggests a prominent role for mtDNA in aging; across many taxa males almost invariably live shorter lives than females. Here we discuss this new work and identify some areas of future research that might now be encouraged to explore what may be the underpinning cause of the strong sexual dimorphism in aging. Editor's suggested further reading in BioEssays: Mitonuclear match: Optimizing fitness and fertility over generations drives ageing within generations Abstract Mitochondrial manoeuvres: Latest insights and hypotheses on mitochondrial partitioning during mitosis in Saccharomyces cerevisiae Abstract Mitochondria and the culture of the Borg Abstract