fastMitoCalc: an ultra-fast program to estimate mitochondrial DNA copy number from whole-genome sequences

Mitochondrial DNA copy number associated dementia risk by somatic mutations and frailty

Mitochondrial dysfunction is linked to physical impairment and dementia. Mitochondrial DNA copy number (mtDNAcn) from blood may predict cognitive decline and dementia risk, but the effect of somatic mutations or frailty is unknown. We estimated mtDNAcn using fastMitoCalc and microheteroplasmies using mitoCaller, from Whole Genome Sequencing (WGS) data. In 189,566 participants free of dementia at study entry (mean age = 56 ± 8), we examined the association between mtDNAcn and subsequent dementia diagnosis using Cox regression. Cognition was assessed in a subset on average 8.9 years later. We examined the associations between mtDNAcn and cognitive measures using multivariable linear regression, adjusted for demographic factors, mtDNAcn-related parameters, and apolipoprotein E ε4 status. We further stratified by frailty and microheteroplasmies. Over an average follow-up of 13.2 years, 3533 participants developed dementia. Each SD higher mtDNAcn (16) was associated with 4.2% lower all-cause dementia hazard (HR = 0.958, p = 0.030), 6% lower non-AD dementia hazard (HR = 0.933, p = 0.022), and not-AD dementia hazard. The associations between mtDNAcn and all-cause dementia and non-AD dementia were stronger among those who were pre-frail or frail or with higher microheteroplasmies. Higher mtDNAcn was associated with higher DSST scores (p = 0.036) and significant only among those with higher microheteroplasmies or frailty (p = 0.029 and 0.048, respectively). mtDNAcn was also associated with delta TMT and paired associate learning only in pre-frail/frail participants (p = 0.007 and 0.045, respectively). Higher WGS-based mtDNAcn in human blood is associated with lower dementia risk, specifically non-AD dementia, and specific cognitive function. The relationships appear stronger in high somatic mutations or frailty. Future studies are warranted to investigate biological underpinnings.

Peripheral immune cell abundance differences link blood mitochondrial DNA copy number and Parkinson's disease

Mitochondrial dysfunction plays an important role in Parkinson's disease (PD), with mitochondrial DNA copy number (mtDNA-CN) emerging as a potential marker for mitochondrial health. We investigated the links between blood mtDNA-CN and PD severity and risk using the Accelerating Medicines Partnership program for Parkinson's Disease dataset, replicating our results in the UK Biobank. Our findings reveal that reduced blood mtDNA-CN levels are associated with heightened PD risk and increased severity of motor symptoms and olfactory dysfunction. We estimated blood cell composition using complete blood cell profile when available or RNA-sequencing data as a surrogate. After adjusting for blood cell composition, the associations between mtDNA-CN and PD risk and clinical symptoms became non-significant. Bidirectional Mendelian randomization analysis also found no evidence of a direct causal relationship between blood mtDNA-CN and PD susceptibility. Hence peripheral inflammatory immune responses rather than mitochondrial dysfunction underpin these previously identified associations in PD.

Inference of forensic body fluids/tissues based on mitochondrial DNA copy number: a preliminary study

The inference of body fluids and tissues is critical in reconstructing crime scenes and inferring criminal behaviors. Nevertheless, present methods are incompatible with conventional DNA genotyping, and additional testing might result in excessive consumption of forensic scene materials. This study aims to investigate the feasibility of distinguishing common body fluids/tissues through the difference in mitochondrial DNA copy number (mtDNAcn). Four types of body fluids/tissues were analyzed in this study - hair, saliva, semen, and skeletal muscle. MtDNAcn was estimated by dividing the read counts of mitochondrial DNA to that of nuclear DNA (RRmt/nu). Results indicated that there were significant differences in RRmt/nu between different body fluids/tissues. Specifically, hair samples exhibited the highest RRmt/nu (log10RRmt/nu: 4.3 ± 0.28), while semen samples showed the lowest RRmt/nu (log10RRmt/nu: -0.1 ± 0.28). RRmt/nu values for DNA samples without extraction were notably higher (approximately 2.9 times) than those obtained after extraction. However, no significant difference in RRmt/nu was observed between various age and gender groups. Hierarchical clustering and Kmeans clustering analyses showed that body fluids/tissues of the same type clustered closely to each other and could be inferred with high accuracy. In conclusion, this study demonstrated that the simultaneous detection of nuclear and mitochondrial DNA made it possible to perform conventional DNA analyses and body fluid/tissue inference at the same time, thus killing two birds with one stone. Furthermore, mtDNAcn has the potential to serve as a novel and promising biomarker for the identification of body fluids/tissues.

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.

A method for measuring mitochondrial DNA copy number in pediatric populations

The mitochondrion is a multifunctional organelle that modulates multiple systems critical for homeostasis during pathophysiological stress. Variation in mitochondrial DNA (mtDNA) copy number (mtDNAcn), a key mitochondrial change associated with chronic stress, is an emerging biomarker for disease pathology and progression. mtDNAcn can be quantified from whole blood samples using qPCR to determine the ratio of mtDNA to nuclear DNA. However, the collection of blood samples in pediatric populations, particularly in infants and young children, can be technically challenging, yield much smaller volume samples, and can be distressing for the patients and their caregivers. Therefore, we have validated a mtDNAcn assay utilizing DNA from simple buccal swabs (Isohelix SK-2S) and report here it's performance in specimens from infants (age = <12 months). Utilizing qPCR to amplify ∼200 bp regions from two mitochondrial (ND1, ND6) and two nuclear (BECN1, NEB) genes, we demonstrated absolute (100%) concordance with results from low-pass whole genome sequencing (lpWGS). We believe that this method overcomes key obstacles to measuring mtDNAcn in pediatric populations and creates the possibility for development of clinical assays to measure mitochondrial change during pathophysiological stress.

Approaches and challenges in identifying, quantifying, and manipulating dynamic mitochondrial genome variations

Mitochondria, the cellular powerhouses, possess their own unique genetic system, including replication, transcription, and translation. Studying these processes is crucial for comprehending mitochondrial disorders, energy production, and their related diseases. Over the past decades, various approaches have been applied in detecting and quantifying mitochondrial genome variations with also the purpose of manipulation of mitochondria or mitochondrial genome for therapeutics. Understanding the scope and limitations of above strategies is not only fundamental to the understanding of basic biology but also critical for exploring disease-related novel target(s), as well to develop innovative therapies. Here, this review provides an overview of different tools and techniques for accurate mitochondrial genome variations identification, quantification, and discuss novel strategies for the manipulation of mitochondria to develop innovative therapeutic interventions, through combining the insights gained from the study of mitochondrial genetics with ongoing single cell omics combined with advanced single molecular tools.

Subtype-Specific Association of Mitochondrial DNA Copy Number With Poststroke/TIA Outcomes in 10 241 Patients in China

BACKGROUND

Mitochondrial DNA copy number (mtDNA-CN) is associated with the severity and mortality in patients with stroke, but the associations in different stroke subtypes remain unexplored.

METHODS

We conducted an observational prospective cohort analysis on patients with ischemic stroke or transient ischemic attack enrolled in the Third China National Stroke Registry. We applied logistic models to assess the association of mtDNA-CN with functional outcome (modified Rankin Scale score, 3-6 versus 0-2) and Cox proportional hazard models to assess the association with stroke recurrence (treating mortality as a competing risk) and mortality during a 12-month follow-up, adjusting for sex, age, physical activity, National Institutes of Health Stroke Scale at admission, history of stroke and peripheral artery disease, small artery occlusion, and interleukin-6. Subgroup analyses stratified by age and stroke subtypes were conducted.

RESULTS

The Third China National Stroke Registry enrolled 15 166 patients, of which 10 241 with whole-genome sequencing data were retained (mean age, 62.2 [SD, 11.2] years; 68.8% men). The associations between mtDNA-CN and poststroke/transient ischemic attack outcomes were specific to patients aged ≤65 years, with lower mtDNA-CN significantly associated with stroke recurrence in 12 months (subdistribution hazard ratio, 1.15 per SD lower mtDNA-CN [95% CI, 1.04-1.27]; P=5.2×10-3) and higher all-cause mortality in 3 months (hazard ratio, 2.19 [95% CI, 1.41-3.39]; P=5.0×10-4). Across subtypes, the associations of mtDNA-CN with stroke recurrence were specific to stroke of undetermined cause (subdistribution hazard ratio, 1.28 [95% CI, 1.11-1.48]; P=6.6×10-4). In particular, lower mtDNA-CN was associated with poorer functional outcomes in stroke of undetermined cause patients diagnosed with embolic stroke of undetermined source (odds ratio, 1.53 [95% CI, 1.20-1.94]; P=5.4×10-4), which remained significant after excluding patients with recurrent stroke (odds ratio, 1.49 [95% CI, 1.14-1.94]; P=3.0×10-3).

CONCLUSIONS

Lower mtDNA-CN is associated with higher stroke recurrence rate and all-cause mortality, as well as poorer functional outcome at follow-up, among stroke of undetermined cause, embolic stroke of undetermined source, and younger patients.

Unsupervised machine learning identifies distinct ALS molecular subtypes in post-mortem motor cortex and blood expression data

Amyotrophic lateral sclerosis (ALS) displays considerable clinical and genetic heterogeneity. Machine learning approaches have previously been utilised for patient stratification in ALS as they can disentangle complex disease landscapes. However, lack of independent validation in different populations and tissue samples have greatly limited their use in clinical and research settings. We overcame these issues by performing hierarchical clustering on the 5000 most variably expressed autosomal genes from motor cortex expression data of people with sporadic ALS from the KCL BrainBank (N = 112). Three molecular phenotypes linked to ALS pathogenesis were identified: synaptic and neuropeptide signalling, oxidative stress and apoptosis, and neuroinflammation. Cluster validation was achieved by applying linear discriminant analysis models to cases from TargetALS US motor cortex (N = 93), as well as Italian (N = 15) and Dutch (N = 397) blood expression datasets, for which there was a high assignment probability (80-90%) for each molecular subtype. The ALS and motor cortex specificity of the expression signatures were tested by mapping KCL BrainBank controls (N = 59), and occipital cortex (N = 45) and cerebellum (N = 123) samples from TargetALS to each cluster, before constructing case-control and motor cortex-region logistic regression classifiers. We found that the signatures were not only able to distinguish people with ALS from controls (AUC 0.88 ± 0.10), but also reflect the motor cortex-based disease process, as there was perfect discrimination between motor cortex and the other brain regions. Cell types known to be involved in the biological processes of each molecular phenotype were found in higher proportions, reinforcing their biological interpretation. Phenotype analysis revealed distinct cluster-related outcomes in both motor cortex datasets, relating to disease onset and progression-related measures. Our results support the hypothesis that different mechanisms underpin ALS pathogenesis in subgroups of patients and demonstrate potential for the development of personalised treatment approaches. Our method is available for the scientific and clinical community at https://alsgeclustering.er.kcl.ac.uk .

Mitochondrial DNA Copy Number and Heteroplasmy in Monozygotic Twins Discordant for Schizophrenia

Schizophrenia (SZ) is a severe, complex, and common mental disorder with high heritability (80%), an adult age of onset, and high discordance (∼50%) in monozygotic twins (MZ). Extensive studies on familial and non-familial cases have implicated a number of segregating mutations and de novo changes in SZ that may include changes to the mitochondrial genome. Yet, no single universally causal variant has been identified, highlighting its extensive genetic heterogeneity. This report specifically focuses on the assessment of changes in the mitochondrial genome in a unique set of monozygotic twins discordant (MZD) for SZ using blood. Genomic DNA from six pairs of MZD twins and two sets of parents (N = 16) was hybridized to the Affymetrix Human SNP Array 6.0 to assess mitochondrial DNA copy number (mtDNA-CN). Whole genome sequencing (WGS) and quantitative polymerase chain reaction (qPCR) was performed for a subset of MZD pairs and their parents and was also used to derive mtDNA-CN estimates. The WGS data were further analyzed to generate heteroplasmy (HP) estimates. Our results show that mtDNA-CN estimates for within-pair and mother-child differences were smaller than comparisons involving unrelated individuals, as expected. MZD twins showed discordance in mtDNA-CN estimates and displayed concordance in directionality of differences for mtDNA-CN across all technologies. Further, qPCR performed better than Affymetrix in estimating mtDNA-CN based on relatedness. No reliable differences in HP were detected between MZD twins. The within-MZD differences in mtDNA-CN observed represent postzygotic somatic changes that may contribute to discordance of MZ twins for diseases, including SZ.

The importance of genome sequencing: unraveling SSBP1 variant missed by exome sequencing

BACKGROUND

Single-stranded DNA-binding protein 1 (SSBP1) plays an essential role in mitochondrial DNA (mtDNA) replication and maintenance, as well as development of retina. Here, we describe the clinical findings and genetic basis of a family with two members affected with bilateral optic atrophy.

MATERIALS AND METHODS

Clinical data were retrospectively collected from an electronic medical record system. Genetic results were obtained using exome sequencing (ES) and genome sequencing (GS).

RESULTS

A 36-year-old man presented with low vision in both eyes since early childhood, with a best-corrected visual acuity of 20/500 in both eyes. He exhibited generalized optic atrophy and diffuse retinal nerve fiber layer thinning without retinal degeneration in both eyes. The family history was consistent with autosomal dominant traits. ES was performed; however, we did not identify any pathogenic variants in the known dominant optic atrophy genes. Subsequently, GS was performed, and it revealed a novel heterozygous c.364A>G p.(Lys122Glu) variant in SSBP1. In silico prediction supported it as deleterious, while segregation analysis detected it in his affected mother and his unaffected sister. No foveopathy or retinal degeneration was observed in the patient's family members.

CONCLUSIONS

We report a novel pathogenic heterozygous SSBP1 variant in a family with autosomal dominant optic atrophy and incomplete penetrance. Furthermore, we demonstrated that GS is advantageous over ES even for the discovery of coding variants, providing uniform coverage. Therefore, GS should be emphasized to improve the molecular diagnostic rate of inherited optic neuropathy.

A bioinformatics pipeline for estimating mitochondrial DNA copy number and heteroplasmy levels from whole genome sequencing data

Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

Single-cell mtDNA heteroplasmy in colorectal cancer

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

Personality traits are consistently associated with blood mitochondrial DNA copy number estimated from genome sequences in two genetic cohort studies

Background

Mitochondrial DNA copy number (mtDNAcn) in tissues and blood can be altered in conditions like diabetes and major depression and may play a role in aging and longevity. However, little is known about the association between mtDNAcn and personality traits linked to emotional states, metabolic health, and longevity. This study tests the hypothesis that blood mtDNAcn is related to personality traits and mediates the association between personality and mortality.

Methods

We assessed the big five personality domains and facets using the Revised NEO Personality Inventory (NEO-PI-R), assessed depressive symptoms with the Center for Epidemiologic Studies Depression Scale (CES-D), estimated mtDNAcn levels from whole-genome sequencing, and tracked mortality in participants from the Baltimore Longitudinal Study of Aging. Results were replicated in the SardiNIA Project.

Results

We found that mtDNAcn was negatively associated with the Neuroticism domain and its facets and positively associated with facets from the other four domains. The direction and size of the effects were replicated in the SardiNIA cohort and were robust to adjustment for potential confounders in both samples. Consistent with the Neuroticism finding, higher depressive symptoms were associated with lower mtDNAcn. Finally, mtDNAcn mediated the association between personality and mortality risk.

Conclusions

To our knowledge, this is the first study to show a replicable association between mtDNAcn and personality. Furthermore, the results support our hypothesis that mtDNAcn is a biomarker of the biological process that explains part of the association between personality and mortality.

Funding

Support for this work was provided by the Intramural Research Program of the National Institute on Aging (Z01-AG000693, Z01-AG000970, and Z01-AG000949) and the National Institute of Neurological Disorders and Stroke of the National Institutes of Health. AT was also supported by the National Institute on Aging of the National Institutes of Health Grant R01AG068093.

Mitochondrial DNA copy number and heteroplasmy load correlate with skeletal muscle oxidative capacity by P31 MR spectroscopy

The association between blood‐based estimates of mitochondrial DNA parameters, mitochondrial DNA copy number (mtDNA‐CN) and heteroplasmy load, with skeletal muscle bioenergetic capacity was evaluated in 230 participants of the Baltimore Longitudinal Study of Aging (mean age:74.7 years, 53% women). Participants in the study sample had concurrent data on muscle oxidative capacity (τ_PCr) assessed by ³¹P magnetic resonance spectroscopy, and mitochondrial DNA parameters estimated from whole‐genome sequencing data. In multivariable linear regression models, adjusted for age, sex, extent of phosphocreatine (PCr) depletion, autosomal sequencing coverage, white blood cell total, and differential count, as well as platelet count, mtDNA‐CN and heteroplasmy load were not significantly associated with τ_PCr (both p > 0.05). However, in models evaluating whether the association between mtDNA‐CN and τ_PCr varied by heteroplasmy load, there was a significant interaction between mtDNA‐CN and heteroplasmy load (p = 0.037). In stratified analysis, higher mtDNA‐CN was significantly associated with lower τ_PCr among participants with high heteroplasmy load (n = 84, β (SE) = −0.236 (0.115), p‐value = 0.044), but not in those with low heteroplasmy load (n = 146, β (SE) = 0.046 (0.119), p‐value = 0.702). Taken together, mtDNA‐CN and heteroplasmy load provide information on muscle bioenergetics. Thus, mitochondrial DNA parameters may be considered proxy measures of mitochondrial function that can be used in large epidemiological studies, especially when comparing subgroups.

Emerging methods for and novel insights gained by absolute quantification of mitochondrial DNA copy number and its clinical applications

The past thirty years have seen a surge in interest in pathophysiological roles of mitochondria, and the accurate quantification of mitochondrial DNA copy number (mCN) in cells and tissue samples is a fundamental aspect of assessing changes in mitochondrial health and biogenesis. Quantification of mCN between studies is surprisingly variable due to a combination of physiological variability and diverse protocols being used to measure this endpoint. The advent of novel methods to quantify nucleic acids like digital polymerase chain reaction (dPCR) and high throughput sequencing offer the ability to measure absolute values of mCN. We conducted an in-depth survey of articles published between 1969 -- 2020 to create an overview of mCN values, to assess consensus values of tissue-specific mCN, and to evaluate consistency between methods of assessing mCN. We identify best practices for methods used to assess mCN, and we address the impact of using specific loci on the mitochondrial genome to determine mCN. Current data suggest that clinical measurement of mCN can provide diagnostic and prognostic value in a range of diseases and health conditions, with emphasis on cancer and cardiovascular disease, and the advent of means to measure absolute mCN should improve future clinical applications of mCN measurements.

Pathogenetic and Prognostic Implications of Increased Mitochondrial Content in Multiple Myeloma

Many studies over the last 20 years have investigated the role of mitochondrial DNA (mtDNA) alterations in carcinogenesis. However, the status of the mtDNACN in MM and its implication in the pathogenesis of the disease remains unclear. We examined changes in plasma cell mtDNACN across different stages of MM by applying RT-PCR and high-throughput sequencing analysis. We observed a significant increase in the average mtDNACN in myeloma cells compared with healthy plasma cells (157 vs. 40 copies; p = 0.02). We also found an increase in mtDNACN in SMM and newly diagnosed MM (NDMM) paired samples and in consecutive relapses in the same patient. Survival analysis revealed the negative impact of a high mtDNACN in progression-free survival in NDMM (p = 0.005). Additionally, we confirmed the higher expression of mitochondrial biogenesis regulator genes in myeloma cells than in healthy plasma cells and we detected single nucleotide variants in several genes involved in mtDNA replication. Finally, we found that there was molecular similarity between "rapidly-progressing SMM" and MM regarding mtDNACN. Our data provide evidence that malignant transformation of myeloma cells involves the activation of mitochondrial biogenesis, resulting in increased mtDNA levels, and highlights vulnerabilities and potential therapeutic targets in the treatment of MM. Accordingly, mtDNACN tracking might guide clinical decision-making and management of complex entities such as high-risk SMM.

Mitochondrial DNA copy number in human disease: the more the better?

Most of the genetic information has been lost or transferred to the nucleus during the evolution of mitochondria. Nevertheless, mitochondria have retained their own genome that is essential for oxidative phosphorylation (OXPHOS). In mammals, a gene‐dense circular mitochondrial DNA (mtDNA) of about 16.5 kb encodes 13 proteins, which constitute only 1% of the mitochondrial proteome. Mammalian mtDNA is present in thousands of copies per cell and mutations often affect only a fraction of them. Most pathogenic human mtDNA mutations are recessive and only cause OXPHOS defects if present above a certain critical threshold. However, emerging evidence strongly suggests that the proportion of mutated mtDNA copies is not the only determinant of disease but that also the absolute copy number matters. In this review, we critically discuss current knowledge of the role of mtDNA copy number regulation in various types of human diseases, including mitochondrial disorders, neurodegenerative disorders and cancer, and during ageing. We also provide an overview of new exciting therapeutic strategies to directly manipulate mtDNA to restore OXPHOS in mitochondrial diseases.

Asthma and its relationship to mitochondrial copy number: Results from the Asthma Translational Genomics Collaborative (ATGC) of the Trans-Omics for Precision Medicine (TOPMed) program

Background

Mitochondria support critical cellular functions, such as energy production through oxidative phosphorylation, regulation of reactive oxygen species, apoptosis, and calcium homeostasis.

Objective

Given the heightened level of cellular activity in patients with asthma, we sought to determine whether mitochondrial DNA (mtDNA) copy number measured in peripheral blood differed between individuals with and without asthma.

Methods

Whole genome sequence data was generated as part of the Trans-Omics for Precision Medicine (TOPMed) Program on participants from the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-ethnicity (SAPPHIRE) and the Study of African Americans, Asthma, Genes, & Environment II (SAGE II). We restricted our analysis to individuals who self-identified as African American (3,651 asthma cases and 1,344 controls). Mitochondrial copy number was estimated using the sequencing read depth ratio for the mitochondrial and nuclear genomes. Respiratory complex expression was assessed using RNA-sequencing.

Results

Average mitochondrial copy number was significantly higher among individuals with asthma when compared with controls (SAPPHIRE: 218.60 vs. 200.47, P<0.001; SAGE II: 235.99 vs. 223.07, P<0.001). Asthma status was significantly associated with mitochondrial copy number after accounting for potential explanatory variables, such as participant age, sex, leukocyte counts, and mitochondrial haplogroup. Despite the consistent relationship between asthma status and mitochondrial copy number, the latter was not associated with time-to-exacerbation or patient-reported asthma control. Mitochondrial respiratory complex gene expression was disproportionately lower in individuals with asthma when compared with individuals without asthma and other protein-encoding genes.

Conclusions

We observed a robust association between asthma and higher mitochondrial copy number. Asthma having an effect on mitochondria function was also supported by lower respiratory complex gene expression in this group.

Thinking outside the nucleus: Mitochondrial DNA copy number in health and disease

Mitochondrial DNA copy number (mtDNA-CN) is a biomarker of mitochondrial function and levels of mtDNA-CN have been reproducibly associated with overall mortality and a number of age-related diseases, including cardiovascular disease, chronic kidney disease, and cancer. Recent advancements in techniques for estimating mtDNA-CN, in particular the use of DNA microarrays and next-generation sequencing data, have led to the comprehensive assessment of mtDNA-CN across these and other diseases and traits. The importance of mtDNA-CN measures to disease and these advancing technologies suggest the potential for mtDNA-CN to be a useful biomarker in the clinic. While the exact mechanism(s) underlying the association of mtDNA-CN with disease remain to be elucidated, we review the existing literature which supports roles for inflammatory dynamics, immune function and alterations to cell signaling as consequences of variation in mtDNA-CN. We propose that future studies should focus on characterizing longitudinal, cell-type and cross-tissue profiles of mtDNA-CN as well as improving methods for measuring mtDNA-CN which will expand the potential for its use as a clinical biomarker.

Investigating Human Mitochondrial Genomes in Single Cells

Mitochondria host multiple copies of their own small circular genome that has been extensively studied to trace the evolution of the modern eukaryotic cell and discover important mutations linked to inherited diseases. Whole genome and exome sequencing have enabled the study of mtDNA in a large number of samples and experimental conditions at single nucleotide resolution, allowing the deciphering of the relationship between inherited mutations and phenotypes and the identification of acquired mtDNA mutations in classical mitochondrial diseases as well as in chronic disorders, ageing and cancer. By applying an ad hoc computational pipeline based on our MToolBox software, we reconstructed mtDNA genomes in single cells using whole genome and exome sequencing data obtained by different amplification methodologies (eWGA, DOP-PCR, MALBAC, MDA) as well as data from single cell Assay for Transposase Accessible Chromatin with high-throughput sequencing (scATAC-seq) in which mtDNA sequences are expected as a byproduct of the technology. We show that assembled mtDNAs, with the exception of those reconstructed by MALBAC and DOP-PCR methods, are quite uniform and suitable for genomic investigations, enabling the study of various biological processes related to cellular heterogeneity such as tumor evolution, neural somatic mosaicism and embryonic development.

Measuring biological aging in humans: A quest

The global population of individuals over the age of 65 is growing at an unprecedented rate and is expected to reach 1.6 billion by 2050. Most older individuals are affected by multiple chronic diseases, leading to complex drug treatments and increased risk of physical and cognitive disability. Improving or preserving the health and quality of life of these individuals is challenging due to a lack of well-established clinical guidelines. Physicians are often forced to engage in cycles of "trial and error" that are centered on palliative treatment of symptoms rather than the root cause, often resulting in dubious outcomes. Recently, geroscience challenged this view, proposing that the underlying biological mechanisms of aging are central to the global increase in susceptibility to disease and disability that occurs with aging. In fact, strong correlations have recently been revealed between health dimensions and phenotypes that are typical of aging, especially with autophagy, mitochondrial function, cellular senescence, and DNA methylation. Current research focuses on measuring the pace of aging to identify individuals who are "aging faster" to test and develop interventions that could prevent or delay the progression of multimorbidity and disability with aging. Understanding how the underlying biological mechanisms of aging connect to and impact longitudinal changes in health trajectories offers a unique opportunity to identify resilience mechanisms, their dynamic changes, and their impact on stress responses. Harnessing how to evoke and control resilience mechanisms in individuals with successful aging could lead to writing a new chapter in human medicine.

Evaluation of mitochondrial DNA copy number estimation techniques

Mitochondrial DNA copy number (mtDNA-CN), a measure of the number of mitochondrial genomes per cell, is a minimally invasive proxy measure for mitochondrial function and has been associated with several aging-related diseases. Although quantitative real-time PCR (qPCR) is the current gold standard method for measuring mtDNA-CN, mtDNA-CN can also be measured from genotyping microarray probe intensities and DNA sequencing read counts. To conduct a comprehensive examination on the performance of these methods, we use known mtDNA-CN correlates (age, sex, white blood cell count, Duffy locus genotype, incident cardiovascular disease) to evaluate mtDNA-CN calculated from qPCR, two microarray platforms, as well as whole genome (WGS) and whole exome sequence (WES) data across 1,085 participants from the Atherosclerosis Risk in Communities (ARIC) study and 3,489 participants from the Multi-Ethnic Study of Atherosclerosis (MESA). We observe mtDNA-CN derived from WGS data is significantly more associated with known correlates compared to all other methods (p < 0.001). Additionally, mtDNA-CN measured from WGS is on average more significantly associated with traits by 5.6 orders of magnitude and has effect size estimates 5.8 times more extreme than the current gold standard of qPCR. We further investigated the role of DNA extraction method on mtDNA-CN estimate reproducibility and found mtDNA-CN estimated from cell lysate is significantly less variable than traditional phenol-chloroform-isoamyl alcohol (p = 5.44x10-4) and silica-based column selection (p = 2.82x10-7). In conclusion, we recommend the field moves towards more accurate methods for mtDNA-CN, as well as re-analyze trait associations as more WGS data becomes available from larger initiatives such as TOPMed.

Quantitative mitochondrial DNA copy number determination using droplet digital PCR with single-cell resolution

Mitochondria are involved in a number of diverse cellular functions, including energy production, metabolic regulation, apoptosis, calcium homeostasis, cell proliferation, and motility, as well as free radical generation. Mitochondrial DNA (mtDNA) is present at hundreds to thousands of copies per cell in a tissue-specific manner. mtDNA copy number also varies during aging and disease progression and therefore might be considered as a biomarker that mirrors alterations within the human body. Here, we present a new quantitative, highly sensitive droplet digital PCR (ddPCR) method, droplet digital mitochondrial DNA measurement (ddMDM), to measure mtDNA copy number not only from cell populations but also from single cells. Our developed assay can generate data in as little as 3 h, is optimized for 96-well plates, and also allows the direct use of cell lysates without the need for DNA purification or nuclear reference genes. We show that ddMDM is able to detect differences between samples whose mtDNA copy number was close enough as to be indistinguishable by other commonly used mtDNA quantitation methods. By utilizing ddMDM, we show quantitative changes in mtDNA content per cell across a wide variety of physiological contexts including cancer progression, cell cycle progression, human T cell activation, and human aging.

Influence of cell distribution and diabetes status on the association between mitochondrial DNA copy number and aging phenotypes in the InCHIANTI study

Mitochondrial DNA copy number (mtDNA-CN) estimated in whole blood is a novel marker of mitochondrial mass and function that can be used in large population-based studies. Analyses that attempt to relate mtDNA-CN to specific aging phenotypes may be confounded by differences in the distribution of blood cell types across samples. Also, low or high mtDNA-CN may have a different meaning given the presence of diseases associated with mitochondrial damage. We evaluated the impact of blood cell type distribution and diabetes status on the association between mtDNA-CN and aging phenotypes, namely chronologic age, interleukin-6, hemoglobin, and all-cause mortality, among 672 participants of the InCHIANTI study. After accounting for white blood cell count, platelet count, and white blood cell proportions in multivariate models, associations of mtDNA-CN with age and interleukin-6 were no longer statistically significant. Evaluation of a statistical interaction by diabetes status suggested heterogeneity of effects in the analysis of mortality (P < 0.01). The magnitude and direction of associations between mtDNA-CN estimated from blood samples and aging phenotypes are influenced by the sample cell type distribution and disease status. Therefore, accounting for these factors may aid understanding of the relevance of mitochondrial DNA copy number to health and aging.