Mitochondrial DNA Heteroplasmy Research Articles

High-Throughput Single-Cell, Single-Mitochondrial DNA Assay Using Hydrogel Droplet Microfluidics.

There is growing interest in understanding the biological implications of single cell heterogeneity and heteroplasmy of mitochondrial DNA (mtDNA), but current methodologies for single-cell mtDNA analysis limit the scale of analysis to small cell populations. Although droplet microfluidics have increased the throughput of single-cell genomic, RNA, and protein analysis, their application to sub-cellular organelle analysis has remained a largely unsolved challenge. Here, we introduce an agarose-based droplet microfluidic approach for single-cell, single-mtDNA analysis, which allows simultaneous processing of hundreds of individual mtDNA molecules within >10,000 individual cells. Our microfluidic chip encapsulates individual cells in agarose beads, designed to have a sufficiently dense hydrogel network to retain mtDNA after lysis and provide a robust scaffold for subsequent multi-step processing and analysis. To mitigate the impact of the high viscosity of agarose required for mtDNA retention on the throughput of microfluidics, we developed a parallelized device, successfully achieving ~95 % mtDNA retention from single cells within our microbeads at >700,000 drops/minute. To demonstrate utility, we analyzed specific regions of the single-mtDNA using a multiplexed rolling circle amplification (RCA) assay. We demonstrated compatibility with both microscopy, for digital counting of individual RCA products, and flow cytometry for higher throughput analysis.

High‐Throughput Single‐Cell, Single‐Mitochondrial DNA Assay Using Hydrogel Droplet Microfluidics

AbstractThere is growing interest in understanding the biological implications of single cell heterogeneity and heteroplasmy of mitochondrial DNA (mtDNA), but current methodologies for single‐cell mtDNA analysis limit the scale of analysis to small cell populations. Although droplet microfluidics have increased the throughput of single‐cell genomic, RNA, and protein analysis, their application to sub‐cellular organelle analysis has remained a largely unsolved challenge. Here, we introduce an agarose‐based droplet microfluidic approach for single‐cell, single‐mtDNA analysis, which allows simultaneous processing of hundreds of individual mtDNA molecules within >10,000 individual cells. Our microfluidic chip encapsulates individual cells in agarose beads, designed to have a sufficiently dense hydrogel network to retain mtDNA after lysis and provide a robust scaffold for subsequent multi‐step processing and analysis. To mitigate the impact of the high viscosity of agarose required for mtDNA retention on the throughput of microfluidics, we developed a parallelized device, successfully achieving ~95 % mtDNA retention from single cells within our microbeads at >700,000 drops/minute. To demonstrate utility, we analyzed specific regions of the single‐mtDNA using a multiplexed rolling circle amplification (RCA) assay. We demonstrated compatibility with both microscopy, for digital counting of individual RCA products, and flow cytometry for higher throughput analysis.

Vitamin B12 status and folic acid supplementation influence mitochondrial heteroplasmy levels in mice.

One-carbon metabolism is a complex network of metabolic reactions that are essential for cellular function including DNA synthesis. Vitamin B12 and folate are micronutrients that are utilized in this pathway and their deficiency can result in the perturbation of one-carbon metabolism and subsequent perturbations in DNA replication and repair. This effect has been well characterized in nuclear DNA but to date, mitochondrial DNA (mtDNA) has not been investigated extensively. Mitochondrial variants have been associated with several inherited and age-related disease states; therefore, the study of factors that impact heteroplasmy are important for advancing our understanding of the mitochondrial genome's impact on human health. Heteroplasmy studies require robust and efficient mitochondrial DNA enrichment to carry out in-depth mtDNA sequencing. Many of the current methods for mtDNA enrichment can introduce biases and false-positive results. Here, we use a method that overcomes these limitations and have applied it to assess mitochondrial heteroplasmy in mouse models of altered one-carbon metabolism. Vitamin B12 deficiency was found to cause increased levels of mitochondrial DNA heteroplasmy across all tissues that were investigated. Folic acid supplementation also contributed to elevated mitochondrial DNA heteroplasmy across all mouse tissues investigated. Heteroplasmy analysis of human data from the Framingham Heart Study suggested a potential sex-specific effect of folate and vitamin B12 status on mitochondrial heteroplasmy. This is a novel relationship that may have broader consequences for our understanding of one-carbon metabolism, mitochondrial-related disease and the influence of nutrients on DNA mutation rates.

The Relationship between Mitochondrial Genome Mutations in Monocytes and the Development of Obesity and Coronary Heart Disease.

Metabolic disorders, including obesity, are often accompanied by an increased risk of cardiovascular complications. Monocytes are the common link between obesity and cardiovascular diseases (CVDs). The bias of innate cellular immunity towards pro-inflammatory activation stimulates the development of diseases associated with chronic inflammation, in particular metabolic disorders, including obesity, as well as CVDs. Disorders in the functional state of monocytes and activation of inflammation may be associated with mitochondrial dysfunction. Mutations accumulating in mitochondrial DNA with age may lead to mitochondrial dysfunction and may be considered a potential marker for developing chronic inflammatory diseases. The present study aimed to study the relationship between mitochondrial heteroplasmy in CD14+ monocytes and cardiovascular risk factors in 22 patients with obesity and coronary heart disease (CHD) by comparing them to 22 healthy subjects. It was found that single-nucleotide variations (SNV) A11467G have a negative correlation with total cholesterol (r = -0.82, p < 0.05), low density lipoproteins (LDL) (r = -0.82, p < 0.05), with age (r = -0.57, p < 0.05) and with mean carotid intima-media thickness (cIMT) (r = -0.43, p < 0.05) and a positive correlation with HDL level (r = 0.71, p < 0.05). SNV 576insC positively correlated with body mass index (BMI) (r = 0.60, p < 0.001) and LDL level (r = 0.43, p < 0.05). SNV A1811G positively correlated with mean cIMT (r = 0.60, p < 0.05). It was revealed that some variants of mitochondrial DNA (mtDNA) heteroplasmy are associated with CVD risk factors. The results demonstrate the potential for using these molecular genetic markers to develop personalized CVD and metabolic disorder treatments.

Twists in genetic mitochondria: Unraveling a case of mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes overlapping with myoclonic epilepsy with ragged red fibers

This case report presents a 23-year-old female with a complex medical history, initially diagnosed with juvenile myoclonic epilepsy, later discovered to have the m.3243A&gt;G variant associated with mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes (MELAS). Despite a typical myoclonic epilepsy with ragged red fibers (MERRF) phenotype, genetic testing confirmed a singular MELAS mutation. The patient exhibited diverse symptoms, including seizures, chronic migraines, myoclonus, and visual disturbances. Diagnostic studies revealed basal ganglia calcifications and progressive brain MRI abnormalities consistent with mitochondrial dysfunction. Additionally, ophthalmologic exams identified features indicative of MELAS pigmentary retinopathy. The discussion highlights the clinical and genetic diversity of mitochondrial disorders, emphasizing the challenges in distinguishing overlapping phenotypes. The patient’s mitochondrial DNA heteroplasmy influenced the variable onset and severity of symptoms. The case underscores the importance of comprehensive genetic investigations, as the m.3243A&gt;G variant is associated with both MELAS and MERRF, leading to diagnostic complexities. Management strategies, primarily symptom-focused due to the lack of a standardized treatment approach, include prophylactic arginine, taurine, and coenzyme Q10 supplementation. The patient’s response to various medications, including vagal nerve stimulator placement, Botox injections, and novel treatments like ASP0367, demonstrates the ongoing efforts to address symptoms and improve quality of life. This report contributes to the understanding of mitochondrial overlap syndromes, offering insights into the diagnostic challenges and management complexities associated with MELAS and MERRF. The case underscores the need for a multidisciplinary approach, combining clinical, genetic, and therapeutic considerations, to optimize care for individuals with mitochondrial encephalopathies.

The metabolic control theory: Its development and its application to mitochondrial oxidative phosphorylation

Metabolic Control Theory (MCT) and Metabolic Control Analysis (MCA) are the two sides, theoretical and experimental, of the measurement of the sensitivity of metabolic networks in the vicinity of a steady state. We will describe the birth and the development of this theory from the first analyses of linear pathways up to a global mathematical theory applicable to any metabolic network. We will describe how the theory, given the global nature of mitochondrial oxidative phosphorylation, solved the problem of what controls mitochondrial ATP synthesis and then how it led to a better understanding of the differential tissue expression of human mitochondrial pathologies and of the heteroplasmy of mitochondrial DNA, leading to the concept of the threshold effect.

Vanadium Pentoxide Exposure Causes Strain-Dependent Changes in Mitochondrial DNA Heteroplasmy, Copy Number, and Lesions, but Not Nuclear DNA Lesions.

Interstitial lung diseases (ILDs) are lethal lung diseases characterized by pulmonary inflammation and progressive lung interstitial scarring. We previously developed a mouse model of ILD using vanadium pentoxide (V2O5) and identified several gene candidates on chromosome 4 associated with pulmonary fibrosis. While these data indicated a significant genetic contribution to ILD susceptibility, they did not include any potential associations and interactions with the mitochondrial genome that might influence disease risk. To conduct this pilot work, we selected the two divergent strains we previously categorized as V2O5-resistant C57BL6J (B6) and -responsive DBA/2J (D2) and compared their mitochondrial genome characteristics, including DNA variants, heteroplasmy, lesions, and copy numbers at 14- and 112-days post-exposure. While we did not find changes in the mitochondrial genome at 14 days post-exposure, at 112 days, we found that the responsive D2 strain exhibited significantly fewer mtDNA copies and more lesions than control animals. Alongside these findings, mtDNA heteroplasmy frequency decreased. These data suggest that mice previously shown to exhibit increased susceptibility to pulmonary fibrosis and inflammation sustain damage to the mitochondrial genome that is evident at 112 days post-V2O5 exposure.

Quantitative Analysis of Mitochondrial DNA Heteroplasmy in Urinary Podocytes of Myoclonus Epilepsy With Ragged-Red Fibers Syndrome

Quantitative Analysis of Mitochondrial DNA Heteroplasmy in Urinary Podocytes of Myoclonus Epilepsy With Ragged-Red Fibers Syndrome

CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome

Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former’s inefficiency and the latter’s inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional ‘Correction’ mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in ‘Correction’ mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.

Neuropathological hallmarks in autopsied cases with mitochondrial diseases caused by the mitochondrial 3243A>G mutation.

The mitochondrial (m.) 3243A>G mutation is known to be associated with various mitochondrial diseases including mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS). Their clinical symptoms have been estimated to occur with an increased mitochondrial DNA (mtDNA) heteroplasmy and reduced activity of oxidative phosphorylation (OXPHOS) complexes, but their trends in the central nervous system remain unknown. Six autopsied mutant cases and three disease control cases without the mutation were enrolled in this study. The mutant cases had a disease duration of 1-27 years. Five of six mutant cases were compatible with MELAS. In the mutant cases, cortical lesions including a laminar necrosis were frequently observed in the parietal, lateral temporal, and occipital lobes; less frequently in the frontal lobe including precentral gyrus; and not at all in the medial temporal lobe. The mtDNA heteroplasmy in brain tissue samples of the mutant cases was strikingly high, ranging from 53.8% to 85.2%. The medial temporal lobe was preserved despite an inhospitable environment having high levels of mtDNA heteroplasmy and lactic acid. OXPHOS complex I was widely decreased in the mutant cases. The swelling of smooth muscle cells in the vessels on the leptomeninges, with immunoreactivity (IR) against mitochondria antibody, and a decreased nuclear/cytoplasmic ratio of choroidal epithelial cells were observed in all mutant cases but in none without the mutation. Common neuropathological findings such as cortical laminar necrosis and basal ganglia calcification were not always observed in the mutant cases. A high level of mtDNA heteroplasmy was observed throughout the brain in spite of heterogeneous cortical lesions. A lack of medial temporal lesion, mitochondrial vasculopathy in vessels on the leptomeninges, and an increased cytoplasmic size of epithelial cells in the choroid plexus could be neuropathological hallmarks helpful in the diagnosis of mitochondrial diseases.

Novel economical, accurate, sensitive, single-cell analytical method for mitochondrial DNA quantification in mtDNA mutation carriers.

Although a variety of analytical methods have been developed to detect mitochondrial DNA (mtDNA) heteroplasmy, there are special requirements of mtDNA heteroplasmy quantification for women carrying mtDNA mutations receiving the preimplantation genetic diagnosis (PGD) and prenatal diagnosis (PD) in clinic. These special requirements include various sample types, large samplenumber, long-term follow-up, and the need for detection of single-cell from biopsied embryos. Therefore, developing an economical, accurate, high-sensitive, and single-cell analytical method for mtDNA heteroplasmy is necessary. In this study, we developed the Sanger sequencing combined droplet digital polymerase chain reaction (ddPCR) method for mtDNA quantification and compared the results to next-generation sequencing (NGS). A total of seventeen families with twelve mtDNA mutations were recruited in this study. The results showed that both Sanger sequencing and ddPCR could be used to analyze the mtDNA heteroplasmy in single-cell samples. There was no statistically significant difference in heteroplasmy levels in common samples with high heteroplasmy (≥ 5%), low heteroplasmy (< 5%), and single-cell samples, either between Sanger sequencing and NGS methods, or between ddPCR and NGS methods (P > 0.05). However, Sanger sequencing was unable to detect extremely low heteroplasmy accurately. But even in samples with extremely low heteroplasmy (0.40% and 0.92%), ddPCR was always able to quantify them. Compared to NGS, Sanger sequencing combined ddPCR analytical methods greatly reduced the cost of sequencing. In conclusion, this study successfully established an economical, accurate, sensitive, single-cell analytical method based on the Sanger sequencing combined ddPCR methods for mtDNA heteroplasmy quantification in aclinical setting.

Analysis of MT-ATP8 gene variants reported in patients by modeling in silico and in yeast model organism

Defects in ATP synthase functioning due to the substitutions in its two mitochondrially encoded subunits a and 8 lead to untreatable mitochondrial diseases. Defining the character of variants in genes encoding these subunits is challenging due to their low frequency, heteroplasmy of mitochondrial DNA in patients’ cells and polymorphisms of mitochondrial genome. We successfully used yeast S. cerevisiae as a model to study the effects of variants in MT-ATP6 gene and our research led to understand how eight amino acid residues substitutions impact the proton translocation through the channel formed by subunit a and c-ring of ATP synthase at the molecular level. Here we applied this approach to study the effects of the m.8403T>C variant in MT-ATP8 gene. The biochemical data from yeast mitochondria indicate that equivalent mutation is not detrimental for the yeast enzyme functioning. The structural analysis of substitutions in subunit 8 introduced by m.8403T>C and five other variants in MT-ATP8 provides indications about the role of subunit 8 in the membrane domain of ATP synthase and potential structural consequences of substitutions in this subunit.

Quantitative haplotype-resolved analysis of mitochondrial DNA heteroplasmy in Human single oocytes, blastoids, and pluripotent stem cells.

Maternal mitochondria are the sole source of mtDNA for every cell of the offspring. Heteroplasmic mtDNA mutations inherited from the oocyte are a common cause of metabolic diseases and associated with late-onset diseases. However, the origin and dynamics of mtDNA heteroplasmy remain unclear. We used our individual Mitochondrial Genome sequencing (iMiGseq) technology to study mtDNA heterogeneity, quantitate single nucleotide variants (SNVs) and large structural variants (SVs), track heteroplasmy dynamics, and analyze genetic linkage between variants at the individual mtDNA molecule level in single oocytes and human blastoids. Our study presented the first single-mtDNA analysis of the comprehensive heteroplasmy landscape in single human oocytes. Unappreciated levels of rare heteroplasmic variants well below the detection limit of conventional methods were identified in healthy human oocytes, of which many are reported to be deleterious and associated with mitochondrial disease and cancer. Quantitative genetic linkage analysis revealed dramatic shifts of variant frequency and clonal expansions of large SVs during oogenesis in single-donor oocytes. iMiGseq of a single human blastoid suggested stable heteroplasmy levels during early lineage differentiation of naïve pluripotent stem cells. Therefore, our data provided new insights of mtDNA genetics and laid a foundation for understanding mtDNA heteroplasmy at early stages of life.

Avoiding misleading estimates using mtDNA heteroplasmy statistics to study bottleneck size and selection.

Mitochondrial DNA heteroplasmy samples can shed light on vital developmental and genetic processes shaping mitochondrial DNA populations. The sample means and sample variance of a set of heteroplasmy observations are typically used both to estimate bottleneck sizes and to perform fits to the theoretical "Kimura" distribution in seeking evidence for mitochondrial DNA selection. However, each of these applications raises problems. Sample statistics do not generally provide optimal fits to the Kimura distribution and so can give misleading results in hypothesis testing, including false positive signals of selection. Using sample variance can give misleading results for bottleneck size estimates, particularly for small samples. These issues can and do lead to false positive results for mitochondrial DNA mechanisms-all published experimental datasets we re-analyzed, reported as displaying departures from the Kimura model, do not in fact give evidence for such departures. Here we outline a maximum likelihood approach that is simple to implement computationally and addresses all of these issues. We advocate the use of maximum likelihood fits and explicit hypothesis tests, not fits and Kolmogorov-Smirnov tests via summary statistics, for ongoing work with mitochondrial DNA heteroplasmy.

Mitochondrial DNA Sequencing and Heteroplasmy Quantification by Next Generation Sequencing.

Over the last 10years, next generation sequencing (NGS) became the gold standard for both diagnosis and discovery of new disease genes responsible for heterogeneous disorders, such as mitochondrial encephalomyopathies. The application of this technology to mtDNA mutations poses extra challenges compared to other genetic conditions because of the peculiarities of mitochondrial genetics and the requirement for proper NGS data management and analysis. Here, we describe a detailed, clinically relevant protocol to sequence the whole mtDNA and quantify heteroplasmy levels of mtDNA variants, starting from total DNA through the generation of a single PCR amplicon.

Author Correction: Segregation of mitochondrial DNA heteroplasmy through a developmental genetic bottleneck in human embryos

Author Correction: Segregation of mitochondrial DNA heteroplasmy through a developmental genetic bottleneck in human embryos

Mitochondrial DNA heteroplasmy distinguishes disease manifestation in PINK1/PRKN-linked Parkinson's disease.

Biallelic mutations in PINK1/PRKN cause recessive Parkinson's disease. Given the established role of PINK1/Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA integrity and inflammation as disease modifiers in carriers of mutations in these genes. Mitochondrial DNA integrity was investigated in a large collection of biallelic (n = 84) and monoallelic (n = 170) carriers of PINK1/PRKN mutations, idiopathic Parkinson's disease patients (n = 67) and controls (n = 90). In addition, we studied global gene expression and serum cytokine levels in a subset. Affected and unaffected PINK1/PRKN monoallelic mutation carriers can be distinguished by heteroplasmic mitochondrial DNA variant load (area under the curve = 0.83, CI 0.74-0.93). Biallelic PINK1/PRKN mutation carriers harbour more heteroplasmic mitochondrial DNA variants in blood (P = 0.0006, Z = 3.63) compared to monoallelic mutation carriers. This enrichment was confirmed in induced pluripotent stem cell-derived (controls, n = 3; biallelic PRKN mutation carriers, n = 4) and post-mortem (control, n = 1; biallelic PRKN mutation carrier, n = 1) midbrain neurons. Last, the heteroplasmic mitochondrial DNA variant load correlated with IL6 levels in PINK1/PRKN mutation carriers (r = 0.57, P = 0.0074). PINK1/PRKN mutations predispose individuals to mitochondrial DNA variant accumulation in a dose- and disease-dependent manner.

Mitochondrial DNA copy number changes, heteroplasmy, and mutations in plasma-derived exosomes and brain tissue of glioblastoma patients

Glioblastoma is the most common malignant tumor of the central nervous system (CNS) in adults. Glioblastoma cells show increased glucose consumption associated with poor prognosis. Since mitochondria play a crucial role in energy metabolism, mutations and copy number changes of mitochondrial DNA may serve as biomarkers. As the brain is difficult to access, analysis of mitochondria directly from the brain tissue represents a challenge. Exosome analysis is an alternative (still poorly explored) approach to investigate molecular changes in CNS tumors.We analyzed brain tissue DNA and plasma-derived exosomal DNA (exoDNA) of 44 glioblastoma patients and 40 control individuals. Quantitative real-time PCR was performed to determine mtDNA copy numbers and the Kruskal-Wallis and Mann-Whitney U test were used for statistical analysis of data. Subsequently, sequencing libraries were prepared and sequenced on the MiSeq platform to identify mtDNA point mutations.Tissue mtDNA copy number was different among controls and patients in multiple comparisons. A similar tendency was detected in exosomes. Based on NGS analysis, several mtDNA point mutations showed slightly different frequencies between cases and controls, but the clinical relevance of these observations is difficult to assess and likely less than that of overall mtDNA copy number changes. Allele frequencies of variants were used to determine the level of heteroplasmy (found to be higher in exo-mtDNA of control individuals).Despite the suggested potential, the use of such biomarkers for the screening and/or diagnosis of glioblastomas is still limited, thus further studies are needed.

Precise and simultaneous quantification of mitochondrial DNA heteroplasmy and copy number by digital PCR

Mitochondrial DNA (mtDNA) is present in multiple copies and phenotypic consequences of mtDNA mutations depend on the mutant load surpassing a specific threshold. Additionally, changes in mtDNA copy number can impact mitochondrial ATP production, resulting in disease. Therefore, the precise determination of mtDNA heteroplasmy and copy number is crucial to the study of mitochondrial diseases. However, current methods can be imprecise, and quantifying small changes in either heteroplasmy or copy number is challenging. We developed a new approach to measure mtDNA heteroplasmy using a single digital PCR (dPCR) probe. This method is based on the observation that fluorescent-labeled probes in dPCR exhibit different intensities depending on the presence of a single nucleotide change in the sequence bound by the probe. This finding allowed us to precisely and simultaneously determine mtDNA copy number and heteroplasmy levels using duplex dPCR. We tested this approach in two different models (human and mouse), which proved faster and more internally controlled when compared to other published methods routinely used in the mitochondrial genetics field. We believe this approach could be broadly applicable to the detection and quantification of other mixed genetic variations.

UNCOVERING THE MISSING PIECE: COMPREHENSIVE REVIEW OF THE TRANSMISSION OF MITOCHONDRIAL DNA HETEROPLASMY TO THE PREIMPLANTATION EMBRYO

The objectives for this study were to: 1) Summarize the methods and key findings from published literature regarding mitochondrial DNA (mtDNA) heteroplasmy transmission and rate of de novo mutations seen in offspring; 2) Evaluate the utility of genome-wide mtDNA heteroplasmy analysis in preimplantation embryos; and 3) Discuss technical and ethical challenges of whole-genome mtDNA analysis in preimplantation embryos.