A family exhibiting heteroplasmy in the human mitochondrial DNA control region reveals both somatic mosaicism and pronounced segregation of mitotypes

A family exhibiting heteroplasmy at position 16 355 in hypervariable region I of the human mtDNA control region has been identified. This family consists of a mother, daughter, and son. DNA samples extracted from blood stains, buccal swabs, and hairs from these individuals were amplified by PCR and sequenced utilizing fluoresence-labeled dye terminator chemistry in an automated DNA sequencer. In both the daughter and mother, heteroplasmy was observed in DNA extracted from blood stains, buccal swabs, and hairs. In the blood stains, the proportion of cytosine was greater than thymine in both individuals. Buccal swab extracts showed a more balanced contribution from the two nucleotides. Telogenic hair root and hair shaft samples exhibited a wide range of nucleotide contributions at this position, from predominately cytosine in some samples to predominately thymine in others. The apparent stochastic segregation of mitotypes in hair samples is discussed from a forensic viewpoint, and the mechanism of mtDNA heteroplasmy is considered. Received: 6 November 1996 / Accepted: 13 February 1997     

Screening for mitochondrial DNA heteroplasmy in children at risk for mitochondrial disease

Temporal temperature gradient gel electrophoresis was used to screen 70% of the mtDNA, including all 22 tRNA genes, for heteroplasmy in 75 children with neuromuscular and/or multi-system dysfunction and elevated lactate levels, and in 95 controls. Standard PCR/ASO (allele specific oligonucleotide) and Southern analyses were also employed. Excluding common length variants, heteroplasmy was found in 22 patients and two controls (P < 0.001), with four patients demonstrating heteroplasmy in two locations each. Of the 23 heteroplasmic variants sequenced among the patients, 17 were novel point variants in the control region (CR) and only two involved tRNA genes. Heteroplasmy is highly associated with the disease group, and is predominately found in the CR, an area rarely studied in patient populations. These variants may be pathological mutations or disease markers.     

No relationship found between point heteroplasmy in mitochondrial DNA control region and age range,sex and haplogroup in human hairs

The analysis of heteroplasmy (presence of more than one type of mitochondrial DNA in an individual) is used as a tool in human identification studies, anthropology, and most currently in studies that relate heteroplasmy with longevity. The frequency of heteroplasmy and its correlation with age has been analyzed using different tissues such as blood, muscle, heart, bone and brain and in different regions of mitochondrial DNA, but this analysis had never been performed using hair samples. In this study, samples of hair were sequenced in order to ascertain whether the presence or not of heteroplasmy varied according to age, sex and origin of haplogroup individuals. The samples were grouped by age (3 groups), gender (male and female) and haplogroup of origin (European, African and Native American), and analyzed using the chi-square statistical test (χ²). Based in statistical results obtained, we conclude that there is no relationship between heteroplasmy and sex, age and haplogroup origin using hair samples.     

Nonsynonymous site heteroplasmy in fish mitochondrial DNA

Heteroplasmic nucleotide polymorphisms are rarely observed in wild animal mitochondrial DNA. The occurrence of such site heteroplasmy is expected to be extremely rare at nonsynonymous sites where the number of nucleotide substitutions per site is low due to functional constraints. This report deals with nonsynonymous mitochondrial heteroplasmy from two wild fish species, chum salmon and Japanese flounder. We detected an A/C nonsynonymous heteroplasmic site corresponding to putative amino acids, Ile or Met, in NADH dehydrogenase subunit-5 (ND5) region of chum salmon. The heteroplasmic site was at the 3rd position of 58th codon. As for Japanese flounder we detected a C/T nonsynonymous heteroplasmic site corresponding to putative amino acids, Leu or Pro, in ND4 region. The heteroplasmic site was at the 2nd position of 450th codon. We also verified heteroplasmy at these sites by sequencing cloned fragments.     

Nuclear mitochondrial interplay in the modulation of the homopolymeric tract length heteroplasmy in the control (D-loop) region of the mitochondrial DNA

We have studied the genetic characteristics of a homopolymeric tract length heteroplasmy associated with the 16189C variant in the mtDNA D-loop control region to identify the factor(s) involved in the generation of the length heteroplasmy. The relative proportion of the various lengths of the polycytosines (i.e., the pattern of the length heteroplasmy) is maintained in an individual, and previous evidence shows that it is regenerated de novo following cell divisions. The pattern is maintained in maternally related individuals, suggestive of mtDNA determinants. Of the 38 individuals with the 16189C variant studied, 39% were found to exhibit the (16180)AAACCCCCCCCCCC(16193) variant associated with A16183C polymorphism [(11C)-group], while 53% showed the (16180)AACCCCCCCCCCCC(16193) variant associated with a further A16182C polymorphism [(12C)-group]. Haplotype analysis of the mtDNA revealed a specific association of the longer mean length of the poly[C] in the (12C)-group with haplogroup B. A similar association was also observed in the (11C)-group, but with a novel haplogroup. Cybrid constructions revealed that the involvement of nuclear factor(s) in the generation of the length heteroplasmy is prominent in homopolymeric tract of eight cytosines. The nuclearly coded factor(s) is/are presumably related to the fidelity of the nuclearly coded components of the mitochondrial DNA replication machinery.     

Transmission of mitochondrial DNA heteroplasmy in normal pedigrees

The presence of multiple mitochondrial genotypes (heteroplasmy) has been studied in normal individuals. Six multigenerational normal families were screened for heteroplasmy by PCR of the mitochondrial control region and the cytochrome c oxidase intergenic regions. Two individuals from different families exhibited multiple length polymorphisms in a homopolymeric tract at positions 16 184–16 193 and a grandmother in a third family was heteroplasmic for both cytosine and thymidine at position 15 945. Although the 15 945 T variant comprised 28% of the grandmother’s mitochondrial DNA, this sequence was not present in any of her descendants. Heteroplasmy was detected in 2.5% of the 96 mother-offspring pairs, consistent with the possibility that it may not be rare. Received: 18 August 1997 / Accepted: 10 November 1997     

A new mitochondrial disease associated with mitochondrial DNA heteroplasmy.

A variable combination of developmental delay, retinitis pigmentosa, dementia, seizures, ataxia, proximal neurogenic muscle weakness, and sensory neuropathy occurred in four members of a family and was maternally transmitted. There was no histochemical evidence of mitochondrial myopathy. Blood and muscle from the patients contained two populations of mitochondrial DNA, one of which had a previously unreported restriction site for AvaI. Sequence analysis showed that this was due to a point mutation at nucleotide 8993, resulting in an amino acid change from a highly conserved leucine to arginine in subunit 6 of mitochondrial H(+)-ATPase. There was some correlation between clinical severity and the amount of mutant mitochondrial DNA in the patients; this was present in only small quantities in the blood of healthy elderly relatives in the same maternal line.     

Analysis of phylogenetically reconstructed mutational spectra in human mitochondrial DNA control region

Analysis of mutations in mitochondrial DNA is an important issue in population and evolutionary genetics. To study spontaneous base substitutions in human mitochondrial DNA we reconstructed the mutational spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. An excess of pyrimidine transitions was found both in HVS I and II regions. By means of classification analysis numerous mutational hotspots were revealed in these spectra. Context analysis of hotspots revealed a complex influence of neighboring bases on mutagenesis in the HVS I region. Further statistical analysis suggested that a transient misalignment dislocation mutagenesis operating in monotonous runs of nucleotides play an important role for generating base substitutions in mitochondrial DNA and define context properties of mtDNA. Our results suggest that dislocation mutagenesis in HVS I and II is a fingerprint of errors produced by DNA polymerase gamma in the course of human mitochondrial DNA replication     

Naturally occurring mitochondrial DNA heteroplasmy in the MRL mouse

The MRL/MpJ mouse is an inbred laboratory strain of Mus musculus, known to exhibit enhanced autoimmunity, increased wound healing, and increased regeneration properties. We report the full-length mitochondrial DNA (mtDNA) sequence of the MRL mouse (Accession # EU450583), and characterize the discovery of two naturally occurring heteroplasmic sites. The first is a T3900C substitution in the TPsiC loop of the tRNA methionine gene (tRNA-Met; mt-Tm). The second is a heteroplasmic insertion of 1-6 adenine nucleotides in the A-tract of the tRNA arginine gene (tRNA-Arg; mt-Tr) at positions 9821-9826. The level of heteroplasmy varied independently at these two sites in MRL individuals. The length of the tRNA-Arg A-tract increased with age, but heteroplasmy at the tRNA-Met site did not change with age. The finding of naturally occurring mtDNA heteroplasmy in an inbred strain of mouse makes the MRL mouse a powerful new experimental model for studies designed to explore therapeutic measures to alter the cellular burden of heteroplasmy.     

MitoBamAnnotator: A web-based tool for detecting and annotating heteroplasmy in human mitochondrial DNA sequences

The use of Next-Generation Sequencing of mitochondrial DNA is becoming widespread in biological and clinical research. This, in turn, creates a need for a convenient tool that detects and analyzes heteroplasmy. Here we present MitoBamAnnotator, a user friendly web-based tool that allows maximum flexibility and control in heteroplasmy research. MitoBamAnnotator provides the user with a comprehensively annotated overview of mitochondrial genetic variation, allowing for an in-depth analysis with no prior knowledge in programming.     

Variability of the human mitochondrial DNA control region sequences in the Lithuanian population

The Lithuanians and Latvians are the only two Baltic cultures that survived until today. Since the Neolithic period the native inhabitants of the present-day Lithuanian territory have not been replaced by any other ethnic group. Therefore the genetic characterization of the present-day Lithuanians may shed some light on the early history of the Balts. We have analysed 120 DNA samples from two Lithuanian ethnolinguistic groups (Aukstaiciai and Zemaiciai) by direct sequencing of the first hypervariable segment (HVI) of the control region of mitochondrial DNA (mtDNA) and restriction enzyme digestion for polymorphic site 00073. On the basis of specific nucleotide substitutions the obtained sequences were classified to mtDNA haplogroups. This revealed the presence of almost all European haplogroups (except X) in the Lithuanian sample, including those that expanded through Europe in the Palaeolithic and those whose expansion occurred during the Neolithic. Molecular diversity indices (gene diversity 0.97, nucleotide diversity 0.012 and mean number of pairwise differences 4.5) were within the range usually reported in European populations. No significant differences between Aukstaiciai and Zemaiciai subgroups were found, but some slight differences need further investigation.     

Size polymorphism and heteroplasmy in the mitochondrial DNA of lower vertebrates

The mitochondrial DNA of the bowfin fish and each of two species of treefrogs displays large-scale size variation. Within each species, mitochondrial genomes span more than a 700 base pair range, and the size polymorphism is localized to one portion of the genome. In addition, about 5 percent of the total 357 individuals surveyed were observed to carry two size classes of mtDNA. These findings are among the few documented instances of extensive within-species mtDNA size polymorphism and individual heteroplasmy, and constitute exceptions to previously reached generalizations about the molecular basis of mtDNA variation.     

Generation of VNTRs and heteroplasmy by sequence turnover in the mitochondrial control region of two elephant seal species

We describe an unusual repetitive DNA region located in the 3′ end of the light (L)-strand in the mitochondrial control region of two elephant seal species. The array of tandem repeats shows both VNTR (variable-number tandem repeat) and sequence variation and is absent from 12 compared mammalian species, except for the occurrence in the same location of a distinct repetitive region in rabbit mtDNA and a similar repeat in the harbor seal. The sequence composition and arrangement of the repeats differ considerably between the northern elephant seal (Mirounga angustirostris) and the southern species (M. leonina) despite an estimated divergence time of 1 MY (based on an mtDNA-RNA gene and the nonrepetitive control region). Analysis of repeat sequence relationships within and between species indicate that divergence in sequence and structure of repeats has involved both slippagelike and unequal crossingover processes of turnover, generating very high levels of divergence and heteroplasmy. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

Length heteroplasmy in the first hypervariable segment of the human mtDNA control region.

The first hypervariable segment of the human mtDNA control region contains a homopolymeric tract of cytosines between nt 16184 and 16193, interrupted at position 16189 by a thymine, according to the Cambridge reference sequence. A variant commonly found in population screening is a T-to-C transition at nt 16189, resulting in an uninterrupted homopolymeric tract. Direct sequencing of individuals with this variant produces a characteristic blurred sequence in nucleotides beyond the tract. Sequencing clones from these individuals revealed that this is caused by high levels of length heteroplasmy in the homopolymeric tract and low levels of length heteroplasmy in the four adenines following the tract. We have developed a rapid method involving densitometry of sequencing gels to quantify the relative proportions of different length variants present in an individual. We have used this to study the proportions of length variants in individuals from three twin pairs and two maternal lineages. While unrelated individuals usually have different proportions of length variants, all maternally related individuals studied have the same proportions, even if they are only distantly related. It is not obvious how identical heteroplasmic profiles are maintained in maternally related individuals, but some possible mechanisms are suggested.     

Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.     

Analysis of heteroplasmy in the major noncoding region of mitochondrial DNA in the blood and atherosclerotic plaques of carotid arteries

For identification of somatic mitochondrial DNA (mtDNA) mutations, the mtDNA major noncoding region (D-loop) sequence in blood samples and carotid atherosclerosis plaques from patients with atherosclerosis was analyzed. Five point heteroplasmic positions were observed in 4 of 23 individuals (17%). Only in two cases could heteroplasmy have resulted from somatic mutation, whereas three heteroplasmic positions were found in both vascular tissue and blood. In addition, length heteroplasmy in a polycytosine stretches was registered at nucleotide positions 303–315 in 16 individuals, and also in the 16184–16193 region in four patients. The results suggest that somatic mtDNA mutations can occur during atherosclerosis, but some heteroplasmic mutations may appear in all tissues, possibly being inherited.     

RecA-dependent DNA repair results in increased heteroplasmy of the Arabidopsis mitochondrial genome

Plant mitochondria have very active DNA recombination activities that are responsible for its plastic structures and that should be involved in the repair of double-strand breaks in the mitochondrial genome. Little is still known on plant mitochondrial DNA repair, but repair by recombination is believed to be a major determinant in the rapid evolution of plant mitochondrial genomes. In flowering plants, mitochondria possess at least two eubacteria-type RecA proteins that should be core components of the mitochondrial repair mechanisms. We have performed functional analyses of the two Arabidopsis (Arabidopsis thaliana) mitochondrial RecAs (RECA2 and RECA3) to assess their potential roles in recombination-dependent repair. Heterologous expression in Escherichia coli revealed that RECA2 and RECA3 have overlapping as well as specific activities that allow them to partially complement bacterial repair pathways. RECA2 and RECA3 have similar patterns of expression, and mutants of either display the same molecular phenotypes of increased recombination between intermediate-size repeats, thus suggesting that they act in the same recombination pathways. However, RECA2 is essential past the seedling stage and should have additional important functions. Treatment of plants with several DNA-damaging drugs further showed that RECA3 is required for different recombination-dependent repair pathways that significantly contribute to plant fitness under stress. Replication repair of double-strand breaks results in the accumulation of crossovers that increase the heteroplasmic state of the mitochondrial DNA. It was shown that these are transmitted to the plant progeny, enhancing the potential for mitochondrial genome evolution.     

The distribution of mitochondrial DNA heteroplasmy due to random genetic drift

Cells containing pathogenic mutations in mitochondrial DNA (mtDNA) generally also contain the wild-type mtDNA, a condition called heteroplasmy. The amount of mutant mtDNA in a cell, called the heteroplasmy level, is an important factor in determining the amount of mitochondrial dysfunction and therefore the disease severity. mtDNA is inherited maternally, and there are large random shifts in heteroplasmy level between mother and offspring. Understanding the distribution in heteroplasmy levels across a group of offspring is an important step in understanding the inheritance of diseases caused by mtDNA mutations. Previously, our understanding of the heteroplasmy distribution has been limited to just the mean and variance of the distribution. Here we give equations, adapted from the work of Kimura on random genetic drift, for the full mtDNA heteroplasmy distribution. We describe how to use the Kimura distribution in mitochondrial genetics, and we test the Kimura distribution against human, mouse, and Drosophila data sets.     

An unusual source of apparent mitochondrial heteroplasmy: duplicate mitochondrial control regions in Thalassarche albatrosses

Molecular ecologists, in search of suitable molecular markers, frequently PCR-amplify regions of mitochondrial DNA from total DNA extracts. This approach, although common, is prone to the co-amplification of nuclear copies of transposed DNA sequences (numts), which can then generate apparent mitochondrial sequence heteroplasmy. In this study we describe the discovery of apparent mitochondrial sequence heteroplasmy in Thalassarche albatrosses but eliminate the possibility of true sequence heteroplasmy and numts and instead reveal the source of the apparent heteroplasmy to be a duplicated control region. The two control regions align easily but are not identical in sequence or in length. Comparisons of functionally significant conserved sequence blocks do not provide evidence of degeneration in either duplicate. Phylogenetic analyses of domain I of both control region copies in five Thalassarche species indicate that they are largely evolving in concert; however, a short section within them is clearly evolving independently. To our knowledge this is the first time contrasting evolutionary patterns have been reported for duplicate control regions. Available evidence suggests that this duplication may be taxonomically widespread, so the results presented here should be considered in future evolutionary studies targeting the control region of all Procellariiformes and potentially other closely related avian groups.