Recurrent Tissue-Specific mtDNA Mutations Are Common in Humans

Mitochondrial DNA (mtDNA) variation can affect phenotypic variation; therefore, knowing its distribution within and among individuals is of importance to understanding many human diseases. Intra-individual mtDNA variation (heteroplasmy) has been generally assumed to be random. We used massively parallel sequencing to assess heteroplasmy across ten tissues and demonstrate that in unrelated individuals there are tissue-specific, recurrent mutations. Certain tissues, notably kidney, liver and skeletal muscle, displayed the identical recurrent mutations that were undetectable in other tissues in the same individuals. Using RFLP analyses we validated one of the tissue-specific mutations in the two sequenced individuals and replicated the patterns in two additional individuals. These recurrent mutations all occur within or in very close proximity to sites that regulate mtDNA replication, strongly implying that these variations alter the replication dynamics of the mutated mtDNA genome. These recurrent variants are all independent of each other and do not occur in the mtDNA coding regions. The most parsimonious explanation of the data is that these frequently repeated mutations experience tissue-specific positive selection, probably through replication advantage.     

A sensitive denaturing gradient-Gel electrophoresis assay reveals a high frequency of heteroplasmy in hypervariable region 1 of the human mtDNA control region

A population study of heteroplasmy in the hypervariable region 1 (HV1) portion of the human mtDNA control region was performed. Blood samples from 253 randomly chosen individuals were examined using a sensitive denaturing gradient-gel electrophoresis (DGGE) system. This method is capable of detecting heteroplasmic proportions as low as 1% and virtually all heteroplasmy where the minor component is > or = 5%. Heteroplasmy was observed in 35 individuals (13.8%; 95% confidence interval [CI] 9.6-18.0). Of these individuals, 33 were heteroplasmic at one nucleotide position, whereas 2 were heteroplasmic at two different positions (a condition known as "triplasmy"). Although heteroplasmy occurred at a total of 16 different positions throughout HV1, it was most frequently observed at positions 16093 (n=13) and 16129 (n=6). In addition, the majority of heteroplasmic variants occurred at low proportions and could not be detected by direct sequencing of PCR products. This study indicates that low-level heteroplasmy in HV1 is relatively common and that it occurs at a broad spectrum of sites. Our results corroborate those of other recent reports indicating that heteroplasmy in the control region is more common than was previously believed-a finding that is of potential importance to evolutionary studies and forensic applications that are based on mtDNA variation.     

Detecting Heteroplasmy from High-Throughput Sequencing of Complete Human Mitochondrial DNA Genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.     

Dynamics of mitochondrial heteroplasmy in three families investigated via a repeatable re-sequencing study

Background

Originally believed to be a rare phenomenon, heteroplasmy - the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual - is emerging as an important component of eukaryotic genetic diversity. Heteroplasmies can be used as genetic markers in applications ranging from forensics to cancer diagnostics. Yet the frequency of heteroplasmic alleles may vary from generation to generation due to the bottleneck occurring during oogenesis. Therefore, to understand the alterations in allele frequencies at heteroplasmic sites, it is of critical importance to investigate the dynamics of maternal mtDNA transmission.

Results

Here we sequenced, at high coverage, mtDNA from blood and buccal tissues of nine individuals from three families with a total of six maternal transmission events. Using simulations and re-sequencing of clonal DNA, we devised a set of criteria for detecting polymorphic sites in heterogeneous genetic samples that is resistant to the noise originating from massively parallel sequencing technologies. Application of these criteria to nine human mtDNA samples revealed four heteroplasmic sites.

Conclusions

Our results suggest that the incidence of heteroplasmy may be lower than estimated in some other recent re-sequencing studies, and that mtDNA allelic frequencies differ significantly both between tissues of the same individual and between a mother and her offspring. We designed our study in such a way that the complete analysis described here can be repeated by anyone either at our site or directly on the Amazon Cloud. Our computational pipeline can be easily modified to accommodate other applications, such as viral re-sequencing.     

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.     

Mitochondrial Mutations in Subjects with Psychiatric Disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.     

Homopolymeric tract heteroplasmy in mtDNA from tissues and single oocytes: support for a genetic bottleneck.

While mtDNA polymorphisms at single base positions are common, the overwhelming majority of the mitochondrial genomes within a single individual are usually identical. When there is a point-mutation difference between a mother and her offspring, there may be a complete switching of mtDNA type within a single generation. It is generally assumed that there is a genetic bottleneck whereby a single or small number of founder mtDNA(s) populate the organism, but it is not known at which stages the restriction/amplification of mtDNA subtype(s) occur, and this uncertainty impedes antenatal diagnosis for mtDNA disorders. Length polymorphisms in homopolymeric tracts have been demonstrated in the large noncoding region of mtDNA. We have developed a new method, T-PCR (trimmed PCR), to quantitate heteroplasmy for two of these tracts (D310 and D16189). D310 variation is sufficient to indicate clonal origins of tissues and single oocytes. Tissues from normal individuals often possessed more than one length variant (heteroplasmy). However, there was no difference in the pattern of the length variants between somatic tissues in any control individual when bulk samples were taken. Oocytes from normal women undergoing in vitro fertilization were frequently heteroplasmic for length variants, and in two cases the modal length of the D310 tract differed in individual oocytes from the same woman. These data suggest that a restriction/amplification event, which we attribute to clonal expansion of founder mtDNA(s), has occurred by the time oocytes are mature, although further segregation may occur at a later stage. In contrast to controls, the length distribution of the D310 tract varied between tissues in a patient with heteroplasmic mtDNA rearrangements, suggesting that these mutants influence segregation. These findings have important implications for the genetic counselling of patients with pathogenic mtDNA mutations.     

Base composition analysis of human mitochondrial DNA using electrospray ionization mass spectrometry: a novel tool for the identification and differentiation of humans

In traditional approaches, mitochondrial DNA (mtDNA) variation is exploited for forensic identity testing by sequencing the two hypervariable regions of the human mtDNA control region. To reduce time and labor, single nucleotide polymorphism (SNP) assays are being sought to possibly replace sequencing. However, most SNP assays capture only a portion of the total variation within the desired regions, require a priori knowledge of the position of the SNP in the genome, and are generally not quantitative. Furthermore, with mtDNA, the clustering of SNPs complicates the design of SNP extension primers or hybridization probes. This article describes an automated electrospray ionization mass spectrometry method that can detect a number of clustered SNPs within an amplicon without a priori knowledge of specific SNP positions and can do so quantitatively. With this technique, the base composition of a PCR amplicon, less than 140 nucleotides in length, can be calculated. The difference in base composition between two samples indicates the presence of an SNP. Therefore, no post-PCR analytical construct needs to be developed to assess variation within a fragment. Of the 2754 different mtDNA sequences in the public forensic mtDNA database, nearly 90% could be resolved by the assay. The mass spectrometer is well suited to characterize and quantitate heteroplasmic samples or those containing mixtures. This makes possible the interpretation of mtDNA mixtures (as well as mixtures when assaying other SNPs). This assay can be expanded to assess genetic variation in the coding region of the mtDNA genome and can be automated to facilitate analysis of a large number of samples such as those encountered after a mass disaster.     

Inheritance of mitochondrial DNA recombinants in double-heteroplasmic families: potential implications for phylogenetic analysis

Recently, somatic recombination of human mitochondrial DNA (mtDNA) was discovered in skeletal muscle. To determine whether recombinant mtDNA molecules can be transmitted through the germ line, we investigated two families, each harboring two inherited heteroplasmic mtDNA mutations. Using allele-specific polymerase chain reaction and single-cell and single-molecule mutational analyses, we discovered, in both families, all four possible allelic combinations of the two heteroplasmic mutations (tetraplasmy), the hallmark of mtDNA recombination. We strongly suggest that these recombinant mtDNA molecules were inherited rather than de novo generated somatically, because they (1) are highly abundant and (2) are present in different tissues of maternally related family members, including young individuals. Moreover, the comparison of the complete mtDNA sequence of one of the families with database sequences revealed an irregular, nontreelike pattern of mutations, reminiscent of a reticulation. We therefore propose that certain reticulations of the human mtDNA phylogenetic tree might be explained by recombination of coexisting mtDNA molecules harboring multiple mutations.     

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.     

一种新的线粒体基因组DNA捕获探针的制备及初步应用

目的:建立新的线粒体基因组DNA杂交捕获探针制备方法并用进行初步应用。方法:通过PCR技术扩增特异线粒体DNA片段,并与生物素偶联,最后与标记磁珠的亲和素混合获得捕获探针。并自行制备的线粒体基因组DNA文库捕获探针与肝癌全基因组测序文库进行液相杂交。分离捕获产物后PCR扩增并进行测序分析。结果:成功建立了线粒体基因组杂交捕获探针制备方法并成功分离线粒体基因组DNA;对测序数据的分析显示:90%以上测序数据来自线粒体基因组DNA,且覆盖率达到100%,且均一性良好。检测到的同质性变异位点数量和异质性变异位点数量与全基因组测序数据产生的结果接近(P=0.9152,P=0.8409)。结论:新方法制备的线粒体基因组DNA杂交捕获探针可以从全基因组文库中高效捕获线粒体基因组DNA测序文库。     

Heteroplasmy of the human mtDNA control region remains constant during life

In a longitudinal, retrospective study, we monitored the level of heteroplasmy at nucleotide position (nt) 309 and nt 16189 of the control region of human mtDNA. As a unique source of DNA, we analyzed multiple cervical-cell samples collected, during 1 or 2 decades, from four women with heteroplasmy at either nt 309 or nt 16189. According to accurate, quantitative analysis by solid-phase minisequencing, the level of heteroplasmy remained stable in the cervical-cell samples from all four women during the time studied. We also analyzed autopsy samples from several different tissues, all containing nt 309 in heteroplasmic form, of one of the women, who was deceased. On the basis of our results, heteroplasmy in the control region of mtDNA seems to be inherited and is not the result of somatic age-related 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.     

Human aging and somatic point mutations in mtDNA: A comparative study of generational differences (grandparents and grandchildren)

The accumulation of somatic mutations in mtDNA is correlated with aging. In this work, we sought to identify somatic mutations in the HVS-1 region (D-loop) of mtDNA that might be associated with aging. For this, we compared 31 grandmothers (mean age: 63 ± 2.3 years) and their 62 grandchildren (mean age: 15 ± 4.1 years), the offspring of their daughters. Direct DNA sequencing showed that mutations absent in the grandchildren were detected in a presumably homoplasmic state in three grandmothers and in a heteroplasmic state in an additional 13 grandmothers; no mutations were detected in the remaining 15 grandmothers. However, cloning followed by DNA sequencing in 12 grandmothers confirmed homoplasia in only one of the three mutations previously considered to be homoplasmic and did not confirm heteroplasmy in three out of nine grandmothers found to be heteroplasmic by direct sequencing. Thus, of 12 grandmothers in whom mtDNA was analyzed by cloning, eight were heteroplasmic for mutations not detected in their grandchildren. In this study, the use of genetically related subjects allowed us to demonstrate the occurrence of age-related (> 60 years old) mutations (homoplasia and heteroplasmy). It is possible that both of these situations (homoplasia and heteroplasmy) were a long-term consequence of mitochondrial oxidative phosphorylation that can lead to the accumulation of mtDNA mutations throughout life.     

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5–10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10–95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.     

Whole mitochondrial DNA variations in hippocampal surgical specimens and blood samples with high-throughput sequencing: A case of mesial temporal lobe epilepsy with hippocampal sclerosis

Introduction

Hippocampal sclerosis is the most common lesion in patients with mesial temporal lobe epilepsy. Recently, there has been growing evidence on the involvement of mitochondria also in sporadic forms of epilepsy. In addition, it has been increasingly argued that mitochondrial dysfunction has an important role in epileptogenesis and seizure generation in temporal lobe epilepsy. Although mtDNA polymorphisms have been identified as potential risk factors for neurological diseases, the link between homoplasmy and heteroplasmy within tissues is not clear. We investigated whether mitochondrial DNA (mtDNA) polymorphisms are involved in a case report of a patient with mesial temporal lobe epilepsy-hippocampal sclerosis (MTLE-HS).

Design

We report the whole genome mtDNA deep sequencing results and clinical features of a 36-year-old woman with MTLE-HS. We used pyrosequencing technology to sequence a whole mitochondrial genome isolated from six different regions of her brain and blood. To assess the possible role of mitochondrial DNA variations in affected tissues, we compared all specimens from different regions of the hippocampus and blood.

Results

In total, 35 homoplasmic and 18 heteroplasmic variations have been detected in 6 different regions of the hippocampus and in blood samples. While the samples did not display any difference in homoplasmic variations, it has been shown that hippocampus regions contain more heteroplasmic variations than blood. The number of heteroplasmic variations was highest in the CA2 region of the brain and accumulated in ND2, ND4 and ND5 genes. Also, dentate and subiculum regions of the hippocampus had similar heteroplasmic variation profiles.

Discussion

We present a new rare example of parallel mutation at 16223 position. Our case suggests that defects in mitochondrial function might be underlying the pathogenesis of seizures in temporal lobe epilepsy.     

A new method for analysis of mitochondrial DNA point mutations and assess levels of heteroplasmy

Determination of mitochondrial DNA (mtDNA) heteroplasmy for the diagnosis of patients with mitochondrial disorders is a difficult task due to the coexistence of wild-type and mutant genomes. We have developed a new method for genotyping and quantification of heteroplasmic point mutations in mtDNA based on the SNaPshot technology. We compared the data of this method with the widely used "last hot-cycle" PCR-RFLP method by studying 15 patients carrying mtDNA mutations. We showed that SNaPshot is an accurate, reproducible, and sensitive technique for the determination of heteroplasmic mtDNA mutations in different tissues from patients, and it is a promising system to be used in prenatal and postnatal diagnosis of mtDNA-associated disorders.     

Detection of mosaic pattern of mitochondrial DNA alterations in different populations of cells from the same endometrial tumor

Somatic mitochondrial DNA (mtDNA) alterations including point mutations and microsatellite instability (MSI) have been frequently detected in human cancers. To further explore the extensiveness of mtDNA alterations, we have analyzed the occurrence of somatic mtDNA mutations in different populations of endometrial cancer cells from the same tumor tissues as compared with adjacent non-tumor cells. Laser-captured micro-dissection was used to harvest endometrial cancer cells from separated areas of the same tumor and adjacent normal cells. Total DNA isolated from micro-dissected cells was PCR amplified and analyzed for mtDNA alterations by polyacrylamide gel electrophoresis and DNA sequencing. Multiple mtDNA alterations were detected in different portions of the same tumor. Different populations of endometrial cancer cells carried different patterns of mtDNA mutations. Interestingly, unlike previous reports, most mutations were found to be heteroplasmic. We have demonstrated the occurrence of hyper-variability of mtDNA alterations in a single piece of tumor tissue. Our observations support the hypothesis that the accumulation of mtDNA alterations is random and expands independently. The data presented here showed the heterogeneity of cancer cells in terms of mtDNA alterations in endometrial cancer.     

The inheritance of heteroplasmy in guppies

Guppies Poecilia reticulata from the Rio Grande, Trinidad are heteroplasmic; individuals possess up to nine different-sized mtDNA haplotypes. A PCR survey of mtDNA length variation that included mothers and embryos suggests that a large number of mitochondrial genomes (possibly within a much smaller number of organelles) pass from one generation to the next.