Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

Single-cell mtDNA dynamics in tumors is driven by coregulation of nuclear and mitochondrial genomes

The extent of cell-to-cell variation in tumor mitochondrial DNA (mtDNA) copy number and genotype, and the phenotypic and evolutionary consequences of such variation, are poorly characterized. Here we use amplification-free single-cell whole-genome sequencing (Direct Library Prep (DLP+)) to simultaneously assay mtDNA copy number and nuclear DNA (nuDNA) in 72,275 single cells derived from immortalized cell lines, patient-derived xenografts and primary human tumors. Cells typically contained thousands of mtDNA copies, but variation in mtDNA copy number was extensive and strongly associated with cell size. Pervasive whole-genome doubling events in nuDNA associated with stoichiometrically balanced adaptations in mtDNA copy number, implying that mtDNA-to-nuDNA ratio, rather than mtDNA copy number itself, mediated downstream phenotypes. Finally, multimodal analysis of DLP+ and single-cell RNA sequencing identified both somatic loss-of-function and germline noncoding variants in mtDNA linked to heteroplasmy-dependent changes in mtDNA copy number and mitochondrial transcription, revealing phenotypic adaptations to disrupted nuclear/mitochondrial balance.

Single-cell analysis reveals context-dependent, cell-level selection of mtDNA

Heteroplasmy occurs when wild-type and mutant mitochondrial DNA (mtDNA) molecules co-exist in single cells1. Heteroplasmy levels change dynamically in development, disease and ageing2,3, but it is unclear whether these shifts are caused by selection or drift, and whether they occur at the level of cells or intracellularly. Here we investigate heteroplasmy dynamics in dividing cells by combining precise mtDNA base editing (DdCBE)4 with a new method, SCI-LITE (single-cell combinatorial indexing leveraged to interrogate targeted expression), which tracks single-cell heteroplasmy with ultra-high throughput. We engineered cells to have synonymous or nonsynonymous complex I mtDNA mutations and found that cell populations in standard culture conditions purge nonsynonymous mtDNA variants, whereas synonymous variants are maintained. This suggests that selection dominates over simple drift in shaping population heteroplasmy. We simultaneously tracked single-cell mtDNA heteroplasmy and ancestry, and found that, although the population heteroplasmy shifts, the heteroplasmy of individual cell lineages remains stable, arguing that selection acts at the level of cell fitness in dividing cells. Using these insights, we show that we can force cells to accumulate high levels of truncating complex I mtDNA heteroplasmy by placing them in environments where loss of biochemical complex I activity has been reported to benefit cell fitness. We conclude that in dividing cells, a given nonsynonymous mtDNA heteroplasmy can be harmful, neutral or even beneficial to cell fitness, but that the 'sign' of the effect is wholly dependent on the environment.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

Next-Generation Sequencing to Characterize Mitochondrial Genomic DNA Heteroplasmy

Mitochondria play a very important role in many crucial cellular functions. Each eukaryotic cell contains hundreds of mitochondria with hundreds of mitochondrial genomes. Mutant and wild-type mitochondrial DNA (mtDNA) may co-exist as heteroplasmy and cause human disease. The purpose of the protocols in this article is to simultaneously determine the mtDNA sequence and quantify the heteroplasmy level using parallel sequencing. The protocols include mitochondrial genomic DNA PCR amplification of two full-length products using two distinct sets of PCR primers. The PCR products are mixed at an equimolar ratio, and the samples are then barcoded and sequenced with high-throughput next-generation sequencing technology. This technology is highly sensitive, specific, and accurate in determining mtDNA mutations and the degree/level of heteroplasmy. © 2022 Wiley Periodicals LLC. Basic Protocol 1: PCR amplification of mitochondrial DNA Basic Protocol 2: Analysis of next-generation sequencing of mitochondrial DNA Basic Protocol 3: Mutect2 pipeline for automated sample processing and large-scale data analysis.

Sonlicromanol improves neuronal network dysfunction and transcriptome changes linked to m.3243A>G heteroplasmy in iPSC-derived neurons

Mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes (MELAS) is often caused by an adenine to guanine variant at m.3243 (m.3243A>G) of the MT-TL1 gene. To understand how this pathogenic variant affects the nervous system, we differentiated human induced pluripotent stem cells (iPSCs) into excitatory neurons with normal (low heteroplasmy) and impaired (high heteroplasmy) mitochondrial function from MELAS patients with the m.3243A>G pathogenic variant. We combined micro-electrode array (MEA) measurements with RNA sequencing (MEA-seq) and found reduced expression of genes involved in mitochondrial respiration and presynaptic function, as well as non-cell autonomous processes in co-cultured astrocytes. Finally, we show that the clinical phase II drug sonlicromanol can improve neuronal network activity when treatment is initiated early in development. This was intricately linked with changes in the neuronal transcriptome. Overall, we provide insight in transcriptomic changes in iPSC-derived neurons with high m.3243A>G heteroplasmy, and show the pathology is partially reversible by sonlicromanol.

Mitochondrial DNA polymorphisms and biogenesis genes in primary and metastatic uveal melanoma cell lines

PURPOSE

This study was designed to identify mitochondrial (mt) DNA variations in primary and metastatic uveal melanoma (UM) cell lines and their relation with cell metabolism to gain insight into metastatic progression.

METHOD

The entire mtDNA genomes were sequenced using Sanger sequencing from two primary UM cell lines (92.1 and MEL270) and two cell lines (OMM2.3 and OMM2.5) derived from liver metastases of the MEL270 patient. The mtDNA copy numbers determined by the ratio of nDNA versus mtDNA. qRT-PCR was used to evaluate expression levels of mitochondrial biogenesis genes.

RESULTS

Sequencing showed that cell line MEL270 and metastases-derived OMM2.3 and OMM2.5 cell lines had homoplasmic single nucleotide polymorphisms (SNPs) representing J1c7a haplogroup, whereas 92.1 cells had mtDNA H31a haplogroup. mtDNA copy numbers were significantly higher in primary cell lines. The metastatic UM cells showed down-regulation of POLG, TFAM, NRF-1 and SIRT1 compared to their primary MEL270 cells. PGC-1α was downregulated in 92.1 and upregulated in MEL270, OMM2.3 and OMM2.5.

CONCLUSIONS

Our finding suggests that within metastatic cells, the heteroplasmic SNPs, copy numbers and mitochondrial biogenesis genes are modulated differentially compared to their primary UM cells. Therefore, investigating pathogenic mtDNA variants associated with cancer metabolic susceptibility may provide future therapeutic strategies in metastatic UM.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.

Deep sequencing shows that accumulation of potentially pathogenic mtDNA mutations rather than mtDNA copy numbers may be associated with early embryonic loss

PURPOSE

To explore the relationship between mitochondrial DNA quantity and heteroplasmy and early embryonic loss.

METHODS

A total of 150 villous samples from patients with spontaneous abortion (SA, n = 75) or induced abortion (IA, n = 75) were collected. qPCR and next-generation sequencing (NGS) were used to test mitochondrial DNA quantity and heteroplasmy. Missense mutations with a CADD score > 15 and heteroplasmy ≥ 70% were defined as potentially pathogenic mutations.

RESULTS

With respect to mitochondrial DNA copy numbers, there was no significant difference between the SA and IA groups (median (IQR), 566 (397-791) vs. 614 (457-739); P = 0.768) or between the euploid and aneuploid groups (median (IQR), 516 (345-730) vs. 599 (423-839); P = 0.107). mtDNA copy numbers were not associated with spontaneous abortion using logistic regression analysis (P = 0.196, 95% CI 1.000-1.001). In addition, more patients harbored possibly pathogenic mtDNA mutations in their chorionic villi in the SA group (70.7%, 53/75) compared with the IA group (54.7%, 41/75; P < 0.05). However, there was no statistical difference between the euploid (80%, 24/30) and aneuploid groups (64.4%, 29/45; p = 0.147).

CONCLUSION

Early embryonic loss and the formation of aneuploidy were not related to mtDNA copy number. Patients with spontaneous abortion were more likely to have possibly pathogenic mutations in their mtDNA, and this may assist in purifying pathogenic mtDNA. However, whether the accumulation of these potentially morbific mtDNA mutations caused early embryonic loss requires further investigation.

Maternally inherited mitochondrial respiratory disorders: from pathogenetic principles to therapeutic implications

Maternally inherited mitochondrial respiratory disorders are rare, progressive, and multi-systemic diseases that remain intractable, with no effective therapeutic interventions. Patients share a defective oxidative phosphorylation pathway responsible for mitochondrial ATP synthesis, in most cases due to pathogenic mitochondrial variants transmitted from mother to child or to a rare de novo mutation or large-scale deletion of the mitochondrial genome. The clinical diagnosis of these mitochondrial diseases is difficult due to exceptionally high clinical variability, while their genetic diagnosis has improved with the advent of next-generation sequencing. The mechanisms regulating the penetrance of the mitochondrial variants remain unresolved with the patient's nuclear background, epigenomic regulation, heteroplasmy, mitochondrial haplogroups, and environmental factors thought to act as rheostats. The lack of animal models mimicking the phenotypic manifestations of these disorders has hampered efforts toward curative therapies. Patient-derived cellular paradigms provide alternative models for elucidating the pathogenic mechanisms and screening pharmacological small molecules to enhance mitochondrial function. Recent progress has been made in designing promising approaches to curtail the negative impact of dysfunctional mitochondria and alleviate clinical symptoms: 1) boosting mitochondrial biogenesis; 2) shifting heteroplasmy; 3) reprogramming metabolism; and 4) administering hypoxia-based treatment. Here, we discuss their varying efficacies and limitations and provide an outlook on their therapeutic potential and clinical application.

The Mitochondrial Dysfunction Hypothesis in Autism Spectrum Disorders: Current Status and Future Perspectives

Autism spectrum disorders (ASDs) constitute a set of heterogeneous neurodevelopmental conditions, characterized by a wide genetic variability that has led to hypothesize a polygenic origin. The metabolic profiles of patients with ASD suggest a possible implication of mitochondrial pathways. Although different physiological and biochemical studies reported deficits in mitochondrial oxidative phosphorylation in subjects with ASD, the role of mitochondrial DNA variations has remained relatively unexplored. In this review, we report and discuss very recent evidence to demonstrate the key role of mitochondrial disorders in the development of ASD.

A novel heteroplasmic mutation in mitochondrial tRNAArg gene associated with non-dystrophic myotonias

Non-dystrophic myotonias (NDM) are rare diseases caused by defects in skeletal muscle chloride and sodium ion channels. It is well established that high-energy consuming tissues such as muscular and nervous systems are exclusively dependent on the ATP generation by mitochondria. The mitochondrial dysfunction, which is caused by mitochondrial DNA mutations, played an important role in the pathogenesis of non-dystrophic myotonias. The purpose of this study is to identify mitochondrial tRNA mutations in non-dystrophic myotonias patients. In this study, 45 Iranian patients with non-dystrophic myotonia were investigated for intracellular ATP content and the mutation screening in all the mitochondrial tRNA genes by DNA sequencing. Our findings showed that lymphocyte intracellular ATP is significantly decreased in NDM patients compared with control subjects (p = 0.001). We found nine mutations in mitochondrial tRNA genes, including m.4454 T > C (in the TψC loop of tRNA^Met), m.5568 A > G (tRNA^Trp), m.5794 T > C (in the anticodon loop of tRNA^Cys), novel m.10438 A > T, and m.10462 T > C (in anticodon loop and ACC stem of tRNA^Arg), m.12308 A > G (tRNA^Leu(CUN)) and m.15907 A > G, m.15924 A > G, and m.15928 G > A (in the anticodon stem of tRNA^Thr) in 31 NDM patients. These results suggest that novel m.10438 A > T mutation is involved in NDM patients and reinforces the significant association between this mutation in mitochondrial tRNA^Arg Gene and NDM patients (p = 0.008).