Mitochondrial genome copy number measured by DNA sequencing in human blood is strongly associated with metabolic traits via cell-type composition differences

Biological and translational attributes of mitochondrial DNA copy number: Laboratory perspective to clinical relevance

The mitochondrial DNA copy number (mtDNAcn) plays a vital role in cellular energy metabolism and mitochondrial health. As mitochondria are responsible for adenosine triphosphate production through oxidative phosphorylation, maintaining an appropriate mtDNAcn level is vital for the overall cellular function. Alterations in mtDNAcn have been linked to various diseases, including neurodegenerative disorders, metabolic conditions, and cancers, making it an important biomarker for understanding the disease pathogenesis. The accurate estimation of mtDNAcn is essential for clinical applications. Quantitative polymerase chain reaction and next-generation sequencing are commonly employed techniques with distinct advantages and limitations. Clinically, mtDNAcn serves as a valuable indicator for early diagnosis, disease progression, and treatment response. For instance, in oncology, elevated mtDNAcn levels in blood samples are associated with tumor aggressiveness and can aid in monitoring treatment efficacy. In neurodegenerative diseases such as Alzheimer’s and Parkinson’s, altered mtDNAcn patterns provide insights into disease mechanisms and progression. Understanding and estimating mtDNAcn are critical for advancing diagnostic and therapeutic strategies in various medical fields. As research continues to uncover the implications of mtDNAcn alterations, its potential as a clinical biomarker is likely to expand, thereby enhancing our ability to diagnose and manage complex diseases.

Exploring the interplay between mitochondrial dysfunction, early life adversity and bipolar disorder

OBJECTIVE

Mitochondria are essential for energy production and reactive oxygen species (ROS) generation, with changes in ROS levels or energy demands affecting mitochondrial DNA (mtDNA) copy numbers, indicating mitochondrial function. Early life adversity (ELA) affects mitochondrial dynamics, influencing long-term health. Both ELA and mitochondrial abnormalities have been independently associated with bipolar disorder (BD). This study aims to explore the complex interplay between mitochondrial dysfunction, ELA, and BD.

METHODS

The study included 60 participants diagnosed with BD and 66 healthy controls (HCs). Data were collected using the Childhood Trauma Questionnaire (CTQ), and leukocyte mtDNA copy number (MCN) was determined from blood samples.

RESULTS

The results indicated the CTQ sum scores were significantly higher in the BD group, reflecting greater exposure to ELA. In HCs, a marginally significant nonlinear relationship between the square of the CTQ sum score and MCN was found. Further analysis demonstrated a significant interaction between ELA and BD on MCN (p = 0.023), highlighting a critical connection between ELA and mitochondrial dysfunction in BD and reinforcing its biological underpinnings.

CONCLUSIONS

Future treatments for BD might target mitochondrial dysfunctions related to chronic stress, with potential pharmaceuticals designed to address these issues and mitigate the negative effects of chronic stress.

Engineering mtDNA deletions by reconstituting end joining in human mitochondria

Recent breakthroughs in the genetic manipulation of mitochondrial DNA (mtDNA) have enabled precise base substitutions and the efficient elimination of genomes carrying pathogenic mutations. However, reconstituting mtDNA deletions linked to mitochondrial myopathies remains challenging. Here, we engineered mtDNA deletions in human cells by co-expressing end-joining (EJ) machinery and targeted endonucleases. Using mitochondrial EJ (mito-EJ) and mito-ScaI, we generated a panel of clonal cell lines harboring a ∼3.5 kb mtDNA deletion across the full spectrum of heteroplasmy. Investigating these cells revealed a critical threshold of ∼75% deleted genomes, beyond which oxidative phosphorylation (OXPHOS) protein depletion, metabolic disruption, and impaired growth in galactose-containing media were observed. Single-cell multiomic profiling identified two distinct nuclear gene deregulation responses: one triggered at the deletion threshold and another progressively responding to heteroplasmy. Ultimately, we show that our method enables the modeling of disease-associated mtDNA deletions across cell types and could inform the development of targeted therapies.

Mitochondrial DNA Copy Number as a Potential Biomarker for the Severity of Motor Symptoms and Prognosis in Parkinson's Disease

BACKGROUND

Mitochondrial function influences Parkinson's disease (PD) through the accumulation of pathogenic alpha-synuclein, oxidative stress, impaired autophagy, and neuroinflammation. The mitochondrial DNA copy number (mtDNA-CN), representing the number of mitochondrial DNA copies within a cell, serves as an easily assessable proxy for mitochondrial function.

OBJECTIVE

This study aimed to assess the diagnostic and prognostic capabilities of mtDNA-CN in PD.

METHODS

We assessed mtDNA-CN in blood samples using whole genome sequencing from 405 patients with PD and 200 healthy controls (HC). We examined the relationship between mtDNA-CN levels and motor symptom severity in PD, as well as their association with dementia development in patients with early-PD (within 3 years of diagnosis).

RESULTS

mtDNA-CN levels were significantly lower in patients with PD compared with HC (P = 1.1 × 10-5). A negative correlation was discovered between mtDNA-CN level and motor severity in PD (correlation coefficient = -0.20; P = 0.008). Among 210 patients with early-PD, Cox regression analysis demonstrated an association between lower mtDNA-CN levels and a higher risk of developing dementia (hazard ratio [HR] = 0.41, 95% confidence interval: 0.20-0.86, P = 0.02), even after adjusting for age and blood cell count (HR = 0.41, 95% confidence interval: 0.18-0.92, P = 0.03). However, mtDNA-CN levels did not significantly correlate with motor progression in PD.

CONCLUSION

Our findings suggest that blood mtDNA-CN may function as a diagnostic biomarker for PD and a prognostic marker for dementia in patients with PD. © 2025 International Parkinson and Movement Disorder Society.

Blood-derived mitochondrial DNA copy number is associated with Alzheimer disease, Alzheimer-related biomarkers and serum metabolites

BACKGROUND

Blood-derived mitochondrial DNA copy number (mtDNA-CN) is a proxy measurement of mitochondrial function in the peripheral and central systems. Abnormal mtDNA-CN not only indicates impaired mtDNA replication and transcription machinery but also dysregulated biological processes such as energy and lipid metabolism. However, the relationship between mtDNA-CN and Alzheimer disease (AD) is unclear.

METHODS

We performed two-sample Mendelian randomization (MR) using publicly available summary statistics from GWAS for mtDNA-CN and AD to investigate the causal relationship between mtDNA-CN and AD. We estimated mtDNA-CN using whole-genome sequence data from blood and brain samples of 13,799 individuals from the Alzheimer's Disease Sequencing Project. Linear and Cox proportional hazards models adjusting for age, sex, and study phase were used to assess the association of mtDNA-CN with AD. The association of AD biomarkers and serum metabolites with mtDNA-CN in blood was evaluated in Alzheimer's Disease Neuroimaging Initiative using linear regression. We conducted a causal mediation analysis to test the natural indirect effects of mtDNA-CN change on AD risk through the significantly associated biomarkers and metabolites.

RESULTS

MR analysis suggested a causal relationship between decreased blood-derived mtDNA-CN and increased risk of AD (OR = 0.68; P = 0.013). Survival analysis showed that decreased mtDNA-CN was significantly associated with higher risk of conversion from mild cognitive impairment to AD (HR = 0.80; P = 0.002). We also identified significant associations of mtDNA-CN with brain FDG-PET (β = 0.103; P = 0.022), amyloid-PET (β = 0.117; P = 0.034), CSF amyloid-β (Aβ) 42/40 (β=-0.124; P = 0.017), CSF t-Tau (β = 0.128; P = 0.015), p-Tau (β = 0.140; P = 0.008), and plasma NFL (β=-0.124; P = 0.004) in females. Several lipid species, amino acids, biogenic amines in serum were also significantly associated with mtDNA-CN. Causal mediation analyses showed that about a third of the effect of mtDNA-CN on AD risk was mediated by plasma NFL (P = 0.009), and this effect was more significant in females (P < 0.005).

CONCLUSIONS

Our study indicates that mtDNA-CN measured in blood is predictive of AD and is associated with AD biomarkers including plasma NFL particularly in females. Further, we illustrate that decreased mtDNA-CN possibly increases AD risk through dysregulation of mitochondrial lipid metabolism and inflammation.

Engineering mtDNA Deletions by Reconstituting End-Joining in Human Mitochondria

Recent breakthroughs in the genetic manipulation of mitochondrial DNA (mtDNA) have enabled the precise introduction of base substitutions and the effective removal of genomes carrying harmful mutations. However, the reconstitution of mtDNA deletions responsible for severe mitochondrial myopathies and age-related diseases has not yet been achieved in human cells. Here, we developed a method to engineer specific mtDNA deletions in human cells by co-expressing end-joining (EJ) machinery and targeted endonucleases. As a proof-of-concept, we used mito-EJ and mito-ScaI to generate a panel of clonal cell lines harboring a ∼3.5 kb mtDNA deletion with the full spectrum of heteroplasmy. Investigating these isogenic cells revealed a critical threshold of ∼75% deleted genomes, beyond which cells exhibited depletion of OXPHOS proteins, severe metabolic disruption, and impaired growth in galactose-containing media. Single-cell multiomic analysis revealed two distinct patterns of nuclear gene deregulation in response to mtDNA deletion accumulation; one triggered at the deletion threshold and another progressively responding to increasing heteroplasmy. In summary, the co-expression of mito-EJ and programable nucleases provides a powerful tool to model disease-associated mtDNA deletions in different cell types. Establishing a panel of cell lines with a large-scale deletion at varying levels of heteroplasmy is a valuable resource for understanding the impact of mtDNA deletions on diseases and guiding the development of potential therapeutic strategies.

Highlights

Combining prokaryotic end-joining with targeted endonucleases generates specific mtDNA deletions in human cellsEngineering a panel of cell lines with a large-scale deletion that spans the full spectrum of heteroplasmy75% heteroplasmy is the threshold that triggers mitochondrial and cellular dysfunctionTwo distinct nuclear transcriptional programs in response to mtDNA deletions: threshold-triggered and heteroplasmy-sensing.

Early life exposure to mercury and relationships with telomere length and mitochondrial DNA content in European children

BACKGROUND

Telomere length (TL) and mitochondrial function expressed as mitochondrial DNA copy number (mtDNAcn) are biomarkers of aging and oxidative stress and inflammation, respectively. Methylmercury (MeHg), a common pollutant in fish, induces oxidative stress. We hypothesized that elevated oxidative stress from exposure to MeHg decreases mtDNAcn and shortens TL.

METHODS

Study participants are 6-11-year-old children from the HELIX multi-center birth cohort study, comprising six European countries. Prenatal and postnatal total mercury (THg) concentrations were measured in blood samples, TL and mtDNAcn were determined in child DNA. Covariates and confounders were obtained by questionnaires. Robust regression models were run, considering sociodemographic and lifestyle covariates, as well as fish consumption. Sex, ethnicity, and fish consumption interaction models were also run.

RESULTS

We found longer TL with higher pre- and postnatal THg blood concentrations, even at low-level THg exposure according to the RfD proposed by the US EPA. The prenatal association showed a significant linear relationship with a 3.46 % increase in TL for each unit increased THg. The postnatal association followed an inverted U-shaped marginal non-linear relationship with 1.38 % an increase in TL for each unit increased THg until reaching a cut-point at 0.96 μg/L blood THg, from which TL attrition was observed. Higher pre- and postnatal blood THg concentrations were consistently related to longer TL among cohorts and no modification effect of fish consumption nor children's sex was observed. No association between THg exposure and mtDNAcn was found.

DISCUSSION

We found evidence that THg is associated with TL but the associations seem to be time- and concentration-dependent. Further studies are needed to clarify the mechanism behind the telomere changes of THg and related health effects.

Fetal cord plasma herpesviruses and preeclampsia: an observational cohort study

A previous study suggested that fetal inheritance of chromosomally integrated human herpesvirus 6 (ici-HHV6) is associated with the hypertensive pregnancy disorder preeclampsia (PE). We aimed to study this question utilizing cord plasma samples (n = 1276) of the Finnish Genetics of Preeclampsia Consortium (FINNPEC) cohort: 539 from a pregnancy with PE and 737 without. We studied these samples and 30 placentas from PE pregnancies by a multiplex qPCR for the DNAs of all nine human herpesviruses. To assess the population prevalence of iciHHV-6, we studied whole-genome sequencing data from blood-derived DNA of 3421 biobank subjects. Any herpes viral DNA was detected in only two (0.37%) PE and one (0.14%) control sample (OR 2.74, 95% CI 0.25-30.4). One PE sample contained iciHHV-6B and another HHV-7 DNA. The control's DNA was of iciHHV-6B; the fetus having growth restriction and preterm birth without PE diagnosis. Placentas showed no herpesviruses. In the biobank data, 3 of 3421 subjects (0.08%) had low level HHV-6B but no iciHHV-6. While iciHHV-6 proved extremely rare, both fetuses with iciHHV-6B were growth-restricted, preterm, and from a pregnancy with maternal hypertension. Our findings suggest that human herpesviruses are not a significant cause of PE, whereas iciHHV-6 may pose some fetal risk.

Plausible Role of Mitochondrial DNA Copy Number in Neurodegeneration-a Need for Therapeutic Approach in Parkinson's Disease (PD)

Parkinson's disease (PD) is an advancing age-associated progressive brain disorder which has various diverse factors, among them mitochondrial dysfunction involves in dopaminergic (DA) degeneration. Aging causes a rise in mitochondrial abnormalities which leads to structural and functional modifications in neuronal activity and cell death in PD. This ends in deterioration of mitochondrial function, mitochondrial alterations, mitochondrial DNA copy number (mtDNA CN) and oxidative phosphorylation (OXPHOS) capacity. mtDNA levels or mtDNA CN in PD have reported that mtDNA depletion would be a predisposing factor in PD pathogenesis. To maintain the mtDNA levels, therapeutic approaches have been focused on mitochondrial biogenesis in PD. The depletion of mtDNA levels in PD can be influenced by autophagic dysregulation, apoptosis, neuroinflammation, oxidative stress, sirtuins, and calcium homeostasis. The current review describes the regulation of mtDNA levels and discusses the plausible molecular pathways in mtDNA CN depletion in PD pathogenesis. We conclude by suggesting further research on mtDNA depletion which might show a promising effect in predicting and diagnosing PD.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

ANKLE1 cleaves mitochondrial DNA and contributes to cancer risk by promoting apoptosis resistance and metabolic dysregulation

Alleles within the chr19p13.1 locus are associated with increased risk of both ovarian and breast cancer and increased expression of the ANKLE1 gene. ANKLE1 is molecularly characterized as an endonuclease that efficiently cuts branched DNA and shuttles between the nucleus and cytoplasm. However, the role of ANKLE1 in mammalian development and homeostasis remains unknown. In normal development ANKLE1 expression is limited to the erythroblast lineage and we found that ANKLE1's role is to cleave the mitochondrial genome during erythropoiesis. We show that ectopic expression of ANKLE1 in breast epithelial-derived cells leads to genome instability and mitochondrial DNA (mtDNA) cleavage. mtDNA degradation then leads to mitophagy and causes a shift from oxidative phosphorylation to glycolysis (Warburg effect). Moreover, mtDNA degradation activates STAT1 and expression of epithelial-mesenchymal transition (EMT) genes. Reduction in mitochondrial content contributes to apoptosis resistance, which may allow precancerous cells to avoid apoptotic checkpoints and proliferate. These findings provide evidence that ANKLE1 is the causal cancer susceptibility gene in the chr19p13.1 locus and describe mechanisms by which higher ANKLE1 expression promotes cancer risk.

Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

Influence of environmental factors and genetic variation on mitochondrial DNA copy number

Mitochondrial DNA copy number (mtDNA CN) has been shown to be highly heritable and associated with traits of interest in humans. However, studies are lacking in the literature for livestock species such as beef cattle. In this study, 2,371 individuals from a crossbred beef population comprising the Germplasm Evaluation program from the U.S. Meat Animal Research Center had samples of blood, leucocyte, or semen collected for low-pass sequencing (LPS) that resulted in both nuclear DNA (nuDNA) and mitochondrial DNA (mtDNA) sequence reads. Mitochondrial DNA CN was estimated based on the ratio of mtDNA to nuDNA coverages. Genetic parameters for mtDNA CN were estimated from an animal model based on a genomic relationship matrix (~87K SNP from the nuDNA). Different models were used to test the effects of tissue, sex, age at sample collection, heterosis, and breed composition. Maternal effects, assessed by fitting a maternal additive component and by fitting eleven SNP on the mtDNA, were also obtained. As previously reported, mtDNA haplotypes were used to classify individuals into Taurine haplogroups (T1, T2, T3/T4, and T5). Estimates of heritability when fitting fixed effects in addition to the intercept were moderate, ranging from 0.11 to 0.31 depending on the model. From a model ignoring contemporary group, semen samples had the lowest mtDNA CN, as expected, followed by blood and leucocyte samples (P ≤ 0.001). The effect of sex and the linear and quadratic effects of age were significant (P ≤ 0.02) depending on the model. When significant, females had greater mtDNA CN than males. The effects of heterosis and maternal heterosis were not significant (P ≥ 0.47). The estimates of maternal and mtDNA heritability were near zero (≤0.03). Most of the samples (98%) were classified as haplogroup T3. Variation was observed in the mtDNA within Taurine haplogroups, which enabled the identification of 24 haplotypes. These results suggest that mtDNA CN is under nuclear genetic control and would respond favorably to selection.

Genome-wide association studies of metabolites in Finnish men identify disease-relevant loci

Few studies have explored the impact of rare variants (minor allele frequency < 1%) on highly heritable plasma metabolites identified in metabolomic screens. The Finnish population provides an ideal opportunity for such explorations, given the multiple bottlenecks and expansions that have shaped its history, and the enrichment for many otherwise rare alleles that has resulted. Here, we report genetic associations for 1391 plasma metabolites in 6136 men from the late-settlement region of Finland. We identify 303 novel association signals, more than one third at variants rare or enriched in Finns. Many of these signals identify genes not previously implicated in metabolite genome-wide association studies and suggest mechanisms for diseases and disease-related traits.