An atlas of mitochondrial DNA genotype-phenotype associations in the UK Biobank

Nanobiopsy investigation of the subcellular mtDNA heteroplasmy in human tissues

Mitochondrial function is critical to continued cellular vitality and is an important contributor to a growing number of human diseases. Mitochondrial dysfunction is typically heterogeneous, mediated through the clonal expansion of mitochondrial DNA (mtDNA) variants in a subset of cells in a given tissue. To date, our understanding of the dynamics of clonal expansion of mtDNA variants has been technically limited to the single cell-level. Here, we report the use of nanobiopsy for subcellular sampling from human tissues, combined with next-generation sequencing to assess subcellular mtDNA mutation load in human tissue from mitochondrial disease patients. The ability to map mitochondrial mutation loads within individual cells of diseased tissue samples will further our understanding of mitochondrial genetic diseases.

Genetic variants affecting mitochondrial function provide further insights for kidney disease

BACKGROUND

Chronic kidney disease (CKD) is a complex disorder that has become a high prevalence global health problem, with diabetes being its predominant pathophysiologic driver. Autosomal genetic variation only explains some of the predisposition to kidney disease. Variations in the mitochondrial genome (mtDNA) and nuclear-encoded mitochondrial genes (NEMG) are implicated in susceptibility to kidney disease and CKD progression, but they have not been thoroughly explored. Our aim was to investigate the association of variation in both mtDNA and NEMG with CKD (and related traits), with a particular focus on diabetes.

METHODS

We used the UK Biobank (UKB) and UK-ROI, an independent collection of individuals with type 1 diabetes mellitus (T1DM) patients.

RESULTS

Fourteen mitochondrial variants were associated with estimated glomerular filtration rate (eGFR) in UKB. Mitochondrial variants and haplogroups U, H and J were associated with eGFR and serum variables. Mitochondrial haplogroup H was associated with all the serum variables regardless of the presence of diabetes. Mitochondrial haplogroup X was associated with end-stage kidney disease (ESKD) in UKB. We confirmed the influence of several known NEMG on kidney disease and function and found novel associations for SLC39A13, CFL1, ACP2 or ATP5G1 with serum variables and kidney damage, and for SLC4A1, NUP210 and MYH14 with ESKD. The G allele of TBC1D32-rs113987180 was associated with higher risk of ESKD in patients with diabetes (OR:9.879; CI95%:4.440-21.980; P = 2.0E-08). In UK-ROI, AGXT2-rs71615838 and SURF1-rs183853102 were associated with diabetic nephropathies, and TFB1M-rs869120 with eGFR.

CONCLUSIONS

We identified novel variants both in mtDNA and NEMG which may explain some of the missing heritability for CKD and kidney phenotypes. We confirmed the role of MT-ND5 and mitochondrial haplogroup H on renal disease (serum variables), and identified the MT-ND5-rs41535848G variant, along with mitochondrial haplogroup X, associated with higher risk of ESKD. Despite most of the associations were independent of diabetes, we also showed potential roles for NEMG in T1DM.

Mitochondrial DNA variants and microbiota: An experimental strategy to identify novel therapeutic potential in chronic inflammatory diseases

We previously demonstrated that mice carrying natural mtDNA variants of the FVB/NJ strain (m.7778 G>T in the mt-Atp8 gene in mitochondrial complex V), namely C57BL/6 J-mtFVB/NJ (B6-mtFVB), exhibited (i) partial protection from experimental skin inflammatory diseases in an anti-murine type VII collagen antibody-induced skin inflammation model and psoriasiform dermatitis model; (ii) significantly altered metabolites, including short-chain fatty acids, according to targeted metabolomics of liver, skin and lymph node samples; and (iii) a differential composition of the gut microbiota according to bacterial 16 S rRNA gene sequencing of stool samples compared to wild-type C57BL/6 J (B6) mice. To further dissect these disease-contributing factors, we induced an experimental antibody-induced skin inflammatory disease in gnotobiotic mice. We performed shotgun metagenomic sequencing of caecum contents and untargeted metabolomics of liver, CD4+ T cell, and caecum content samples from conventional B6-mtFVB and B6 mice. We identified D-glucosamine as a candidate mediator that ameliorated disease severity in experimental antibody-induced skin inflammation by modulating immune cell function in T cells, neutrophils and macrophages. Because mice carrying mtDNA variants of the FVB/NJ strain show differential disease susceptibility to a wide range of experimental diseases, including diet-induced atherosclerosis in low-density lipoprotein receptor knockout mice and collagen antibody-induced arthritis in DBA/1 J mice, this experimental approach is valuable for identifying novel therapeutic options for skin inflammatory conditions and other chronic inflammatory diseases to which mice carrying specific mtDNA variants show differential susceptibility.

Assessing the joint effects of mitochondrial genes and physical activity on the psychiatric phenotype of subjective well-being based on the UK Biobank data

Subjective well-being (SWB) is an important measure for mental health status. Previous research has shown that physical activity can affect an individual's well-being, yet the underlying molecular mechanism remains to be clarified. In this study, we aim to evaluate the potential interactions between mitochondrial genes and physical activity (PA) as well as their combined effects on individual well-being. SWB phenotype data in UK Biobank were enrolled for this study including nine aspects such as work/job satisfaction, health satisfaction, family relationship satisfaction, friendships satisfaction, financial situation satisfaction, ever depressed for a whole week, general happiness, general happiness with own health and belief that own life is meaningful. We made analysis for each aspects separately. Firstly, mitochondria-wide association studies (MiWAS) was conducted to assess the association of mitochondrial Single Nucleotide Polymorphisms SNP with each aspect of SWB. Then an interaction analysis of mitochondrial DNA (mtDNA) mutation and PA was performed to evaluate their joint effect on SWB status. Meanwhile, these two analysis were made for female and male group separately as well as the total samples, all under the control of possible confounding factors including gender, age, Townsend Deprivation Index (TDI), education, alcohol consumption, smoking habits, and 10 principal components. MiWAS analysis identified 45 mtSNPs associated with 9 phenotypes of SWB. For example, m.15218A > G on MT-CYB in the health satisfaction phenotype of the total subjects. Gender-specific analyses found 30 mtSNPs in females and 58 in males, involving 13 mtGenes. In mtDNA-PA interaction analysis, we also identified 10 significant mtDNA-PA interaction sets for SWB. For instance, m.13020 T > C (MT-ND5) was associated with the SWB financial situation satisfaction phenotype in all subjects (P = 0.00577). In addition, MiWAS analysis identified 12 mtGene variants associated with SWB, as MT-ND1 and MT-ND2. However, in mtDNA-PA interactions we detected 7 mtDNA affecting psychiatric disorders occurring, as in the friendships satisfaction phenotype (m.3394 T > C on MT-ND1). Our study results suggest an implication of the interaction between mitochondrial function and physical activity in the risk of psychiatric disorder development.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses

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

CD4+ T cell mitochondrial genotype in Multiple Sclerosis: a cross-sectional and longitudinal analysis

Multiple Sclerosis (MS) is a chronic autoimmune demyelinating disease of the central nervous system (CNS), with a largely unknown etiology, where mitochondrial dysfunction likely contributes to neuroaxonal loss and brain atrophy. Mirroring the CNS, peripheral immune cells from patients with MS, particularly CD4+ T cells, show inappropriate mitochondrial phenotypes and/or oxidative phosphorylation (OxPhos) insufficiency, with a still unknown contribution of mitochondrial DNA (mtDNA). We hypothesized that mitochondrial genotype in CD4+ T cells might influence MS disease activity and progression. Thus, we performed a retrospective cross-sectional and longitudinal study on patients with a recent diagnosis of either Clinically Isolated Syndrome (CIS) or Relapsing-Remitting MS (RRMS) at two timepoints: 6 months (VIS1) and 36 months (VIS2) after disease onset. Our primary outcomes were the differences in mtDNA extracted from CD4+ T cells between: (I) patients with CIS/RRMS (PwMS) at VIS1 and age- and sex-matched healthy controls (HC), in the cross-sectional analysis, and (II) different diagnostic evolutions in PwMS from VIS1 to VIS2, in the longitudinal analysis. We successfully performed mtDNA whole genome sequencing (mean coverage: 2055.77 reads/base pair) in 183 samples (61 triplets). Nonetheless, mitochondrial genotype was not associated with a diagnosis of CIS/RRMS, nor with longitudinal diagnostic evolution.

Cytoplasmic genome contributions to domestication and improvement of modern maize

BACKGROUND

Studies on maize evolution and domestication are largely limited to the nuclear genomes, and the contribution of cytoplasmic genomes to selection and domestication of modern maize remains elusive. Maize cytoplasmic genomes have been classified into fertile (NA and NB) and cytoplasmic-nuclear male-sterility (CMS-S, CMS-C, and CMS-T) groups, but their contributions to modern maize breeding have not been systematically investigated.

RESULTS

Here we report co-selection and convergent evolution between nuclear and cytoplasmic genomes by analyzing whole genome sequencing data of 630 maize accessions modern maize and its relatives, including 24 fully assembled mitochondrial and chloroplast genomes. We show that the NB cytotype is associated with the expansion of modern maize to North America, gradually replaces the fertile NA cytotype probably through unequal division, and predominates in over 90% of modern elite inbred lines. The mode of cytoplasmic evolution is increased nucleotypic diversity among the genes involved in photosynthesis and energy metabolism, which are driven by selection and domestication. Furthermore, genome-wide association study reveals correlation of cytoplasmic nucleotypic variation with key agronomic and reproductive traits accompanied with the diversification of the nuclear genomes.

CONCLUSIONS

Our results indicate convergent evolution between cytoplasmic and nuclear genomes during maize domestication and breeding. These new insights into the important roles of mitochondrial and chloroplast genomes in maize domestication and improvement should help select elite inbred lines to improve yield stability and crop resilience of maize hybrids.

Penetrance and expressivity of mitochondrial variants in a large clinically unselected population

Whole genome sequencing (WGS) from large clinically unselected cohorts provides a unique opportunity to assess the penetrance and expressivity of rare and/or known pathogenic mitochondrial variants in population. Using WGS from 179 862 clinically unselected individuals from the UK Biobank, we performed extensive single and rare variant aggregation association analyses of 15 881 mtDNA variants and 73 known pathogenic variants with 15 mitochondrial disease-relevant phenotypes. We identified 12 homoplasmic and one heteroplasmic variant (m.3243A>G) with genome-wide significant associations in our clinically unselected cohort. Heteroplasmic m.3243A>G (MAF = 0.0002, a known pathogenic variant) was associated with diabetes, deafness and heart failure and 12 homoplasmic variants increased aspartate aminotransferase levels including three low-frequency variants (MAF ~0.002 and beta~0.3 SD). Most pathogenic mitochondrial disease variants (n = 66/74) were rare in the population (<1:9000). Aggregated or single variant analysis of pathogenic variants showed low penetrance in unselected settings for the relevant phenotypes, except m.3243A>G. Multi-system disease risk and penetrance of diabetes, deafness and heart failure greatly increased with m.3243A>G level ≥ 10%. The odds ratio of these traits increased from 5.61, 12.3 and 10.1 to 25.1, 55.0 and 39.5, respectively. Diabetes risk with m.3243A>G was further influenced by type 2 diabetes genetic risk. Our study of mitochondrial variation in a large-unselected population identified novel associations and demonstrated that pathogenic mitochondrial variants have lower penetrance in clinically unselected settings. m.3243A>G was an exception at higher heteroplasmy showing a significant impact on health making it a good candidate for incidental reporting.

Mitochondrial genetic variation and risk of chronic kidney disease and acute kidney injury in UK Biobank participants

Experimental models suggest an important role for mitochondrial dysfunction in the pathogenesis of chronic kidney disease (CKD) and acute kidney injury (AKI), but little is known regarding the impact of common mitochondrial genetic variation on kidney health. We sought to evaluate associations of inherited mitochondrial DNA (mtDNA) variation with risk of CKD and AKI in a large population-based cohort. We categorized UK Biobank participants who self-identified as white into eight distinct mtDNA haplotypes, which were previously identified based on their associations with phenotypes associated with mitochondrial DNA copy number, a measure of mitochondrial function. We used linear and logistic regression models to evaluate associations of these mtDNA haplotypes with estimated glomerular filtration rate by serum creatinine and cystatin C (eGFRCr-CysC, N = 362,802), prevalent (N = 416 cases) and incident (N = 405 cases) end-stage kidney disease (ESKD), AKI defined by diagnostic codes (N = 14,170 cases), and urine albumin/creatinine ratio (ACR, N = 114,662). The mean age was 57 ± 8 years and the mean eGFR was 90 ± 14 ml/min/1.73 m2. MtDNA haplotype was significantly associated with eGFR (p = 2.8E-12), but not with prevalent ESKD (p = 5.9E-2), incident ESKD (p = 0.93), AKI (p = 0.26), or urine ACR (p = 0.54). The association of mtDNA haplotype with eGFR remained significant after adjustment for diabetes mellitus and hypertension (p = 1.2E-10). When compared to the reference haplotype, mtDNA haplotypes I (β = 0.402, standard error (SE) = 0.111; p = 2.7E-4), IV (β = 0.430, SE = 0.073; p = 4.2E-9), and V (β = 0.233, SE = 0.050; p = 2.7E-6) were each associated with higher eGFR. Among self-identified white UK Biobank participants, mtDNA haplotype was associated with eGFR, but not with ESKD, AKI or albuminuria.

An individual with Sarmatian-related ancestry in Roman Britain

In the second century CE the Roman Empire had increasing contact with Sarmatians, nomadic Iranian speakers occupying an area stretching from the Pontic-Caspian steppe to the Carpathian mountains, both in the Caucasus and in the Danubian borders of the empire.1,2,3 In 175 CE, following their defeat in the Marcomannic Wars, emperor Marcus Aurelius drafted Sarmatian cavalry into Roman legions and deployed 5,500 Sarmatian soldiers to Britain, as recorded by contemporary historian Cassius Dio.4,5 Little is known about where the Sarmatian cavalry were stationed, and no individuals connected with this historically attested event have been identified to date, leaving its impact on Britain largely unknown. Here we document Caucasus- and Sarmatian-related ancestry in the whole genome of a Roman-period individual (126-228 calibrated [cal.] CE)-an outlier without traceable ancestry related to local populations in Britain-recovered from a farmstead site in present-day Cambridgeshire, UK. Stable isotopes support a life history of mobility during childhood. Although several scenarios are possible, the historical deployment of Sarmatians to Britain provides a parsimonious explanation for this individual's extraordinary life history. Regardless of the factors behind his migrations, these results highlight how long-range mobility facilitated by the Roman Empire impacted provincial locations outside of urban centers.

The genetic puzzle of a SOD1-patient with ocular ptosis and a motor neuron disease: a case report

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease with a complex genetic architecture, showing monogenic, oligogenic, and polygenic inheritance. In this study, we describe the case of a 71 years-old man diagnosed with ALS with atypical clinical features consisting in progressive ocular ptosis and sensorineural deafness. Genetic analyses revealed two heterozygous variants, in the SOD1 (OMIM*147450) and the TBK1 (OMIM*604834) genes respectively, and furthermore mitochondrial DNA (mtDNA) sequencing identified the homoplasmic m.14484T>C variant usually associated with Leber's Hereditary Optic Neuropathy (LHON). We discuss how all these variants may synergically impinge on mitochondrial function, possibly contributing to the pathogenic mechanisms which might ultimately lead to the neurodegenerative process, shaping the clinical ALS phenotype enriched by adjunctive clinical features.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine

INTRODUCTION

Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied.

OBJECTIVES

To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region.

METHODS

We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples.

RESULTS

Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1.

CONCLUSION

Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

Mitochondria-wide association study observed significant interactions of mitochondrial respiratory and the inflammatory in the development of anxiety and depression

The aim of this study was to investigate the possible interaction of mitochondrial dysfunction and inflammatory cytokines in the risk of anxiety and depression. We utilized the UK Biobank for the sample of this study. A mitochondria-wide association(MiWAS) and interaction analysis was performed to investigate the interaction effects of mitochondrial DNA (mtDNA)×C-reactive protein (CRP) on the risks of self-reported anxiety (N = 72,476), general anxiety disorder (GAD-7) scores (N = 80,853), self-reported depression (N = 80,778), Patient Health Questionnaire (PHQ-9) scores (N = 80,520) in total samples, females and males, respectively, adjusting for sex, age, Townsend deprivation index (TDI), education score, alcohol intake, smoking and 10 principal components. In all, 25 mtSNPs and 10 mtSNPs showed significant level of association with self-reported anxiety and GAD-7 score respectively. A total of seven significant mtDNA × CRP interactions were found for anxiety, such as m.3915G>A(MT-ND1) for self-reported anxiety in total subjects (P = 6.59 × 10^-3), m.4561T>C(MT-ND2) (P = 3.04 × 10^-3) for GAD-7 score in total subjects. For depression, MiWAS identified 17 significant mtSNPs for self-reported depression and 14 significant mtSNPs for PHQ-9 scores. 17 significant mtDNA associations (2 for self-reported depression and 15 for PHQ-9 score) was identified, such as m.14869G>A(MT-CYB; P = 2.22 × 10^-3) associated with self-reported depression and m.4561T>C (MT-ND2; P value = 3.02 × 10^-8) associated with PHQ-9 score in all subjects. In addition, 5 common mtDNA shared with anxiety and depression were found in MiWAS, and 4 common mtDNA variants were detected to interact with CRP for anxiety and depression, such as m.9899T>C(MT-CO3). Our study suggests the important interaction effects of mitochondrial function and CRP on the risks of anxiety and depression.

Association of mitochondrial DNA variation with high myopia in a Han Chinese population

High myopia (HM), which is characterized by oxidative stress, is one of the leading causes of visual impairment and blindness across the world. Family and population genetic studies have uncovered nuclear-genome variants in proteins functioned in the mitochondria. However, whether mitochondrial DNA mutations are involved in HM remains unexplored. Here, we performed the first large-scale whole-mitochondrial genome study in 9613 HM cases and 9606 control subjects of Han Chinese ancestry for identifying HM-associated mitochondrial variants. The single-variant association analysis identified nine novel genetic variants associated with HM reaching the entire mitochondrial wide significance level, including rs370378529 in ND2 with an odds ratio (OR) of 5.25. Interestingly, eight out of nine variants were predominantly located in related sub-haplogroups, i.e. m.5261G > A in B4b1c, m.12280A > G in G2a4, m.7912G > A in D4a3b, m.94G > A in D4e1, m.14857 T > C in D4e3, m.14280A > G in D5a2, m.16272A > G in G2a4, m.8718A > G in M71 and F1a3, indicating that the sub-haplogroup background can increase the susceptible risk for high myopia. The polygenic risk score analysis of the target and validation cohorts indicated a high accuracy for predicting HM with mtDNA variants (AUC = 0.641). Cumulatively, our findings highlight the critical roles of mitochondrial variants in untangling the genetic etiology of HM.

Secondary structure of the human mitochondrial genome affects formation of deletions

BACKGROUND

Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions.

RESULTS

By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6-9 kb and another within 13-16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6-9 kb and 13-16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470-8482 bp (base pair) and a second arm at 13,447-13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging.

CONCLUSIONS

Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species.

Mitochondrial DNA variation in Alzheimer's disease reveals a unique microprotein called SHMOOSE

Mitochondrial DNA variants have previously associated with disease, but the underlying mechanisms have been largely elusive. Here, we report that mitochondrial SNP rs2853499 associated with Alzheimer's disease (AD), neuroimaging, and transcriptomics. We mapped rs2853499 to a novel mitochondrial small open reading frame called SHMOOSE with microprotein encoding potential. Indeed, we detected two unique SHMOOSE-derived peptide fragments in mitochondria by using mass spectrometry-the first unique mass spectrometry-based detection of a mitochondrial-encoded microprotein to date. Furthermore, cerebrospinal fluid (CSF) SHMOOSE levels in humans correlated with age, CSF tau, and brain white matter volume. We followed up on these genetic and biochemical findings by carrying out a series of functional experiments. SHMOOSE acted on the brain following intracerebroventricular administration, differentiated mitochondrial gene expression in multiple models, localized to mitochondria, bound the inner mitochondrial membrane protein mitofilin, and boosted mitochondrial oxygen consumption. Altogether, SHMOOSE has vast implications for the fields of neurobiology, Alzheimer's disease, and microproteins.

Mitochondrial Base Editing: Recent Advances towards Therapeutic Opportunities

Mitochondria are critical organelles that form networks within our cells, generate energy dynamically, contribute to diverse cell and organ function, and produce a variety of critical signaling molecules, such as cortisol. This intracellular microbiome can differ between cells, tissues, and organs. Mitochondria can change with disease, age, and in response to the environment. Single nucleotide variants in the circular genomes of human mitochondrial DNA are associated with many different life-threatening diseases. Mitochondrial DNA base editing tools have established novel disease models and represent a new possibility toward personalized gene therapies for the treatment of mtDNA-based disorders.

Oxidized mitochondrial DNA induces gasdermin D oligomerization in systemic lupus erythematosus

Although extracellular DNA is known to form immune complexes (ICs) with autoantibodies in systemic lupus erythematosus (SLE), the mechanisms leading to the release of DNA from cells remain poorly characterized. Here, we show that the pore-forming protein, gasdermin D (GSDMD), is required for nuclear DNA and mitochondrial DNA (mtDNA) release from neutrophils and lytic cell death following ex vivo stimulation with serum from patients with SLE and IFN-γ. Mechanistically, the activation of FcγR downregulated Serpinb1 following ex vivo stimulation with serum from patients with SLE, leading to spontaneous activation of both caspase-1/caspase-11 and cleavage of GSDMD into GSDMD-N. Furthermore, mtDNA oxidization promoted GSDMD-N oligomerization and cell death. In addition, GSDMD, but not peptidyl arginine deiminase 4 is necessary for extracellular mtDNA release from low-density granulocytes from SLE patients or healthy human neutrophils following incubation with ICs. Using the pristane-induced lupus model, we show that disease severity is significantly reduced in mice with neutrophil-specific Gsdmd deficiency or following treatment with the GSDMD inhibitor, disulfiram. Altogether, our study highlights an important role for oxidized mtDNA in inducing GSDMD oligomerization and pore formation. These findings also suggest that GSDMD might represent a possible therapeutic target in SLE.

Insulin Resistance in Mitochondrial Diabetes

Mitochondrial diabetes (MD) is generally classified as a genetic defect of β-cells. The main pathophysiology is insulin secretion failure in pancreatic β-cells due to impaired mitochondrial ATP production. However, several reports have mentioned the presence of insulin resistance (IR) as a clinical feature of MD. As mitochondrial dysfunction is one of the important factors causing IR, we need to focus on IR as another pathophysiology of MD. In this special issue, we first briefly summarized the insulin signaling and molecular mechanisms of IR. Second, we overviewed currently confirmed pathogenic mitochondrial DNA (mtDNA) mutations from the MITOMAP database. The variants causing diabetes were mostly point mutations in the transfer RNA (tRNA) of the mitochondrial genome. Third, we focused on these variants leading to the recently described "tRNA modopathies" and reviewed the clinical features of patients with diabetes. Finally, we discussed the pathophysiology of MD caused by mtDNA mutations and explored the possible mechanism underlying the development of IR. This review should be beneficial to all clinicians involved in diagnostics and therapeutics related to diabetes and mitochondrial diseases.

Mitochondrial haplogroups and cognitive progression in Parkinson's disease

Mitochondria are a culprit in the onset of Parkinson's disease, but their role during disease progression is unclear. Here we used Cox proportional hazards models to exam the effect of variation in the mitochondrial genome on longitudinal cognitive and motor progression over time in 4064 patients with Parkinson's disease. Mitochondrial macro-haplogroup was associated with reduced risk of cognitive disease progression in the discovery and replication population. In the combined analysis, patients with the super macro-haplogroup J, T, U# had a 41% lower risk of cognitive progression with P = 2.42 × 10-6 compared to those with macro-haplogroup H. Exploratory analysis indicated that the common mitochondrial DNA variant, m.2706A>G, was associated with slower cognitive decline with a hazard ratio of 0.68 (95% confidence interval 0.56-0.81) and P = 2.46 × 10-5. Mitochondrial haplogroups were not appreciably linked to motor progression. This initial genetic survival study of the mitochondrial genome suggests that mitochondrial haplogroups may be associated with the pace of cognitive progression in Parkinson's disease over time.

Is the disease risk and penetrance in Leber hereditary optic neuropathy actually low?

Pedigree analysis showed that a large proportion of Leber hereditary optic neuropathy (LHON) family members who carry a mitochondrial risk variant never lose vision. Mitochondrial haplotype appears to be a major factor influencing the risk of vision loss from LHON. Mitochondrial variants, including m.14484T>C and m.11778G>A, have been added to gene arrays, and thus many patients and research participants are tested for LHON mutations. Analysis of the UK Biobank and Australian cohort studies found more than 1 in 1,000 people in the general population carry either the m.14484T>C or the m.11778G>A LHON variant. None of the subset of carriers examined had visual acuity at 20/200 or worse, suggesting a very low penetrance of LHON. Haplogroup analysis of m.14484T>C carriers showed a high rate of haplogroup U subclades, previously shown to have low penetrance in pedigrees. Penetrance calculations of the general population are lower than pedigree calculations, most likely because of modifier genetic factors. This Matters Arising Response paper addresses the Watson et al. (2022) Matters Arising paper, published concurrently in The American Journal of Human Genetics.

The mitochondrial challenge: Disorders and prevention strategies

This short review provides basic knowledge on mitochondrial inheritance, its disorders, and potential ways to overcome them in human reproductive medicine. The latter are currently mostly associated with the so-called mitochondrial replacement (nuclear transfer) procedures, performed at different stages and with slight technical differences. Being promising but obviously highly invasive, these procedures require detailed investigation of their delayed effects on embryogenesis, pregnancy and future health. A special attention is paid to the newest available data on these issues, as well as to their limitations and possible further research directions.

Mitochondrial Unfolded Protein Response and Integrated Stress Response as Promising Therapeutic Targets for Mitochondrial Diseases

The development and application of high-throughput omics technologies have enabled a more in-depth understanding of mitochondrial biosynthesis metabolism and the pathogenesis of mitochondrial diseases. In accordance with this, a host of new treatments for mitochondrial disease are emerging. As an essential pathway in maintaining mitochondrial proteostasis, the mitochondrial unfolded protein response (UPR^mt) is not only of considerable significance for mitochondrial substance metabolism but also plays a fundamental role in the development of mitochondrial diseases. Furthermore, in mammals, the integrated stress response (ISR) and UPR^mt are strongly coupled, functioning together to maintain mitochondrial function. Therefore, ISR and UPR^mt show great application prospects in the treatment of mitochondrial diseases. In this review, we provide an overview of the molecular mechanisms of ISR and UPR^mt and focus on them as potential targets for mitochondrial disease therapy.

A library of base editors for the precise ablation of all protein-coding genes in the mouse mitochondrial genome

The development of curative treatments for mitochondrial diseases, which are often caused by mutations in mitochondrial DNA (mtDNA) that impair energy metabolism and other aspects of cellular homoeostasis, is hindered by an incomplete understanding of the underlying biology and a scarcity of cellular and animal models. Here we report the design and application of a library of double-stranded-DNA deaminase-derived cytosine base editors optimized for the precise ablation of every mtDNA protein-coding gene in the mouse mitochondrial genome. We used the library, which we named MitoKO, to produce near-homoplasmic knockout cells in vitro and to generate a mouse knockout with high heteroplasmy levels and no off-target edits. MitoKO should facilitate systematic and comprehensive investigations of mtDNA-related pathways and their impact on organismal homoeostasis, and aid the generation of clinically meaningful in vivo models of mtDNA dysfunction.

Modulating mitochondrial DNA mutations: factors shaping heteroplasmy in the germ line and somatic cells

Until recently it was thought that most humans only harbor one type of mitochondrial DNA (mtDNA), however, deep sequencing and single-cell analysis has shown the converse - that mixed populations of mtDNA (heteroplasmy) are the norm. This is important because heteroplasmy levels can change dramatically during transmission in the female germ line, leading to high levels causing severe mitochondrial diseases. There is also emerging evidence that low level mtDNA mutations contribute to common late onset diseases such as neurodegenerative disorders and cardiometabolic diseases because the inherited mutation levels can change within developing organs and non-dividing cells over time. Initial predictions suggested that the segregation of mtDNA heteroplasmy was largely stochastic, with an equal tendency for levels to increase or decrease. However, transgenic animal work and single-cell analysis have shown this not to be the case during germ-line transmission and in somatic tissues during life. Mutation levels in specific mtDNA regions can increase or decrease in different contexts and the underlying molecular mechanisms are starting to be unraveled. In this review we provide a synthesis of recent literature on the mechanisms of selection for and against mtDNA variants. We identify the most pertinent gaps in our understanding and suggest ways these could be addressed using state of the art techniques.

Mitochondrial signal transduction

The analogy of mitochondria as powerhouses has expired. Mitochondria are living, dynamic, maternally inherited, energy-transforming, biosynthetic, and signaling organelles that actively transduce biological information. We argue that mitochondria are the processor of the cell, and together with the nucleus and other organelles they constitute the mitochondrial information processing system (MIPS). In a three-step process, mitochondria (1) sense and respond to both endogenous and environmental inputs through morphological and functional remodeling; (2) integrate information through dynamic, network-based physical interactions and diffusion mechanisms; and (3) produce output signals that tune the functions of other organelles and systemically regulate physiology. This input-to-output transformation allows mitochondria to transduce metabolic, biochemical, neuroendocrine, and other local or systemic signals that enhance organismal adaptation. An explicit focus on mitochondrial signal transduction emphasizes the role of communication in mitochondrial biology. This framework also opens new avenues to understand how mitochondria mediate inter-organ processes underlying human health.

Common mitochondrial haplogroups as modifiers of the onset-age for critical COVID-19

As a key regulator of innate immunity, mitochondrial function is essential to maintain antiviral activities. Common mitochondrial DNA variants (haplogroups) have been associated with different physiological capacities and the nrisk of developing several diseases. Haplogroup H was associated with increased survival among sepsis patients, and lower risk of progression toward AIDS in HIV infected and lower manifestation of severe manifestation of herpex virus disease. We studied 316 Spanish with critical COVID-19, and found that the 7028C (haplogroup H) was protective among patients with early-onset disease (≤65 vs > 65 years, p = 0.01), while the ancestral 16223T was a risk factor for early-onset critical COVID-19 (OR = 3.36, 95 %CI = 1.49-7.54). Our work suggested that common mitochondrial variants may serve as predictors of COVID-19 severity. Additional studies to confirm this effect from other populations are of special interest.

Single-cell multiomics analyses of spindle-transferred human embryos suggest a mostly normal embryonic development

Mitochondrial DNA (mtDNA) mutations are often associated with incurable diseases and lead to detectable pathogenic variants in 1 out of 200 babies. Uncoupling of the inheritance of mtDNA and the nuclear genome by spindle transfer (ST) can potentially prevent the transmission of mtDNA mutations from mother to offspring. However, no well-established studies have critically assessed the safety of this technique. Here, using single-cell triple omics sequencing method, we systematically analyzed the genome (copy number variation), DNA methylome, and transcriptome of ST and control blastocysts. The results showed that, compared to that in control embryos, the percentage of aneuploid cells in ST embryos did not significantly change. The epiblast, primitive endoderm, and trophectoderm (TE) of ST blastocysts presented RNA expression profiles that were comparable to those of control blastocysts. However, the DNA demethylation process in TE cells of ST blastocysts was slightly slower than that in the control blastocysts. Collectively, our results suggest that ST seems generally safe for embryonic development, with a relatively minor delay in the DNA demethylation process at the blastocyst stage.

Uncoupling of the inheritance of mtDNA and the nuclear genome by spindle transfer could prevent the transmission of mtDNA mutations. Systematic single-cell multiomic analyses of spindle transferred human embryos suggest this technique seems generally safe for human embryonic development and deserves further scientific evaluation and clinical testing.

Exploring the Effects of Mitonuclear Interactions on Mitochondrial DNA Gene Expression in Humans

Most mitochondrial protein complexes include both nuclear and mitochondrial gene products, which coevolved to work together. This coevolution can be disrupted due to disparity in genetic ancestry between the nuclear and mitochondrial genomes in recently admixed populations. Such mitonuclear DNA discordance might result in phenotypic effects. Several nuclear-encoded proteins regulate expression of mitochondrial DNA (mtDNA) genes. We hypothesized that mitonuclear DNA discordance affects expression of genes encoded by mtDNA. To test this, we utilized the data from the GTEx project, which contains expression levels for ∼100 African Americans and >600 European Americans. The varying proportion of African and European ancestry in recently admixed African Americans provides a range of mitonuclear discordance values, which can be correlated with mtDNA gene expression levels (adjusted for age and ischemic time). In contrast, European Americans did not undergo recent admixture. We demonstrated that, for most mtDNA protein-coding genes, expression levels in energetically-demanding tissues were lower in African Americans than in European Americans. Furthermore, gene expression levels were lower in individuals with higher mitonuclear discordance, independent of population. Moreover, we found a negative correlation between mtDNA gene expression and mitonuclear discordance. In African Americans, the average value of African ancestry was higher for nuclear-encoded mitochondrial than non-mitochondrial genes, facilitating a match in ancestry with the mtDNA and more optimal interactions. These results represent an example of a phenotypic effect of mitonuclear discordance on human admixed populations, and have potential biomedical applications.

MitoVisualize: a resource for analysis of variants in human mitochondrial RNAs and DNA

SUMMARY

We present MitoVisualize, a new tool for analysis of the human mitochondrial DNA (mtDNA). MitoVisualize enables visualization of: (i) the position and effect of variants in mitochondrial transfer RNA and ribosomal RNA secondary structures alongside curated variant annotations, (ii) data across RNA structures, such as to show all positions with disease-associated variants or with post-transcriptional modifications and (iii) the position of a base, gene or region in the circular mtDNA map, such as to show the location of a large deletion. All visualizations can be easily downloaded as figures for reuse. MitoVisualize can be useful for anyone interested in exploring mtDNA variation, though is designed to facilitate mtDNA variant interpretation in particular.

AVAILABILITY AND IMPLEMENTATION

MitoVisualize can be accessed via https://www.mitovisualize.org/. The source code is available at https://github.com/leklab/mito_visualize/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Germline mutations in mitochondrial complex I reveal genetic and targetable vulnerability in IDH1-mutant acute myeloid leukaemia

The interaction of germline variation and somatic cancer driver mutations is under-investigated. Here we describe the genomic mitochondrial landscape in adult acute myeloid leukaemia (AML) and show that rare variants affecting the nuclear- and mitochondrially-encoded complex I genes show near-mutual exclusivity with somatic driver mutations affecting isocitrate dehydrogenase 1 (IDH1), but not IDH2 suggesting a unique epistatic relationship. Whereas AML cells with rare complex I variants or mutations in IDH1 or IDH2 all display attenuated mitochondrial respiration, heightened sensitivity to complex I inhibitors including the clinical-grade inhibitor, IACS-010759, is observed only for IDH1-mutant AML. Furthermore, IDH1 mutant blasts that are resistant to the IDH1-mutant inhibitor, ivosidenib, retain sensitivity to complex I inhibition. We propose that the IDH1 mutation limits the flexibility for citrate utilization in the presence of impaired complex I activity to a degree that is not apparent in IDH2 mutant cells, exposing a mutation-specific metabolic vulnerability. This reduced metabolic plasticity explains the epistatic relationship between the germline complex I variants and oncogenic IDH1 mutation underscoring the utility of genomic data in revealing metabolic vulnerabilities with implications for therapy.

Mitochondrial metabolism has been associated with tumourigenesis in acute myeloid leukaemia (AML) and currently considered as a potential therapeutic target. Here, the authors show, in patients with AML, that germline mutations in mitochondrial complex I are mutually exclusive with somatic mutations in the metabolic enzyme IDH1, and find IDH1 mutant cells have increased sensitivity to complex I inhibitors.

Mitochondria-derived peptides in aging and healthspan

The mechanisms that explain mitochondrial dysfunction in aging and healthspan continue to be studied, but one element has been unexplored: microproteins. Small open reading frames in circular mitochondria DNA can encode multiple microproteins, called mitochondria-derived peptides (MDPs). Currently, eight MDPs have been published: humanin, MOTS-c, and SHLPs 1–6. This Review describes recent advances in microprotein discovery with a focus on MDPs. It discusses what is currently known about MDPs in aging and how this new understanding could add to the way we understand age-related diseases including type 2 diabetes, cancer, and neurodegenerative diseases at the genomic, proteomic, and drug-development levels.

Mitochondrial 1555 G&gt;A variant as a potential risk factor for childhood glioblastoma

Background

Childhood glioblastoma multiforme (GBM) is a highly aggressive disease with low survival, and its etiology, especially concerning germline genetic risk, is poorly understood. Mitochondria play a key role in putative tumorigenic processes relating to cellular oxidative metabolism, and mitochondrial DNA variants were not previously assessed for association with pediatric brain tumor risk.

Methods

We conducted an analysis of 675 mitochondrial DNA variants in 90 childhood GBM cases and 2789 controls to identify enrichment of mitochondrial variant associated with GBM risk. We also performed this analysis for other glioma subtypes including pilocytic astrocytoma. Nuclear-encoded mitochondrial gene variants were also analyzed.

Results

We identified m1555 A>G was significantly associated with GBM risk (adjusted OR 29.30, 95% CI 5.25–163.4, P-value 9.5 X 10^–4). No association was detected for other subtypes. Haplotype analysis further supported the independent risk contributed by m1555 G>A, instead of a haplogroup joint effect. Nuclear-encoded mitochondrial gene variants identified significant associations in European (rs62036057 in WWOX, adjusted OR = 2.99, 95% CI 1.88–4.75, P-value = 3.42 X 10^–6) and Hispanic (rs111709726 in EFHD1, adjusted OR = 3.57, 95% CI 1.99–6.40, P-value = 1.41 X 10^–6) populations in ethnicity-stratified analyses.

Conclusion

We report for the first time a potential role played by a functional mitochondrial ribosomal RNA variant in childhood GBM risk, and a potential role for both mitochondrial and nuclear-mitochondrial DNA polymorphisms in GBM tumorigenesis. These data implicate cellular oxidative metabolic capacity as a contributor to the etiology of pediatric glioblastoma.

The potential of mitochondrial genome engineering

Mitochondria are subject to unique genetic control by both nuclear DNA and their own genome, mitochondrial DNA (mtDNA), of which each mitochondrion contains multiple copies. In humans, mutations in mtDNA can lead to devastating, heritable, multi-system diseases that display different tissue-specific presentation at any stage of life. Despite rapid advances in nuclear genome engineering, for years, mammalian mtDNA has remained resistant to genetic manipulation, hampering our ability to understand the mechanisms that underpin mitochondrial disease. Recent developments in the genetic modification of mammalian mtDNA raise the possibility of using genome editing technologies, such as programmable nucleases and base editors, for the treatment of hereditary mitochondrial disease.

Examining the Association between Mitochondrial Genome Variation and Coronary Artery Disease

Large-scale genome-wide association studies have identified hundreds of single-nucleotide variants (SNVs) significantly associated with coronary artery disease (CAD). However, collectively, these explain <20% of the heritability. Hypothesis: Here, we hypothesize that mitochondrial (MT)-SNVs might present one potential source of this “missing heritability”. Methods: We analyzed 265 MT-SNVs in ~500,000 UK Biobank individuals, exploring two different CAD definitions: a more stringent (myocardial infarction and/or revascularization; HARD = 20,405), and a more inclusive (angina and chronic ischemic heart disease; SOFT = 34,782). Results: In HARD cases, the most significant (p < 0.05) associations were for m.295C>T (control region) and m.12612A>G (ND5), found more frequently in cases (OR = 1.05), potentially related to reduced cardiorespiratory fitness in response to exercise, as well as for m.12372G>A (ND5) and m.11467A>G (ND4), present more frequently in controls (OR = 0.97), previously associated with lower ROS production rate. In SOFT cases, four MT-SNVs survived multiple testing corrections (at FDR < 5%), all potentially conferring increased CAD risk. Of those, m.11251A>G (ND4) and m.15452C>A (CYB) have previously shown significant associations with body height. In line with this, we observed that CAD cases were slightly less physically active, and their average body height was ~2.00 cm lower compared to controls; both traits are known to be related to increased CAD risk. Gene-based tests identified CO2 associated with HARD/SOFT CAD, whereas ND3 and CYB associated with SOFT cases (p < 0.05), dysfunction of which has been related to MT oxidative stress, obesity/T2D (CO2), BMI (ND3), and angina/exercise intolerance (CYB). Finally, we observed that macro-haplogroup I was significantly (p < 0.05) more frequent in HARD cases vs. controls (3.35% vs. 3.08%), potentially associated with response to exercise. Conclusions: We found only spurious associations between MT genome variation and HARD/SOFT CAD and conclude that more MT-SNV data in even larger study cohorts may be needed to conclusively determine the role of MT DNA in CAD.

Heteroplasmy of Wild Type Mitochondrial DNA Variants in Mice Causes Metabolic Heart Disease With Pulmonary Hypertension and Frailty

Background: In most eukaryotic cells, the mitochondrial DNA (mtDNA) is uniparentally transmitted and present in multiple copies derived from the clonal expansion of maternally inherited mtDNA. All copies are therefore near-identical, or homoplasmic. The presence of more than one mtDNA variant in the same cytoplasm can arise naturally or result from new medical technologies aimed at preventing mitochondrial genetic diseases and improving fertility. The latter is called divergent non-pathological mtDNAs heteroplasmy (DNPH). We hypothesized that DNPH is maladaptive and usually prevented by the cell. Methods: We engineered and characterized DNPH mice throughout their lifespan using transcriptomic, metabolomic, biochemical, physiological and phenotyping techniques. We focused on in vivo imaging techniques for non-invasive assessment of cardiac and pulmonary energy metabolism. Results: We show that DNPH impairs mitochondrial function, with profound consequences in critical tissues that cannot resolve heteroplasmy, particularly cardiac and skeletal muscle. Progressive metabolic stress in these tissues leads to severe pathology in adulthood, including pulmonary hypertension and heart failure, skeletal muscle wasting, frailty, and premature death. Symptom severity is strongly modulated by the nuclear context. Conclusions: Medical interventions that may generate DNPH should address potential incompatibilities between donor and recipient mtDNA.

Interactions between nuclear and mitochondrial SNPs and Parkinson's disease risk

Interactions between the products of the nuclear and mitochondrial genomes are critical for the function of most eukaryotic cells. Recently the introduction of mitochondrial replacement therapy has raised the question of incompatibilities between mitochondrial and nuclear variants, and their potential influence on the genetic makeup of human populations. Such interactions could also contribute to the variability of the penetrance of pathogenic DNA variants. This led us to investigate the frequencies of combinations of nuclear and mitochondrial SNP alleles (mitonuclear combinations) in healthy individuals (n = 5375) and in a cohort of patients with Parkinson's disease (PD, n = 2210). In the unaffected population, we were not able to find associations between nuclear and mitochondrial variants with a false discovery rate below 0.05 after accounting for multiple testing (i.e., the number of combinations examined). However, in the PD cohort, five combinations surpassed this threshold. Next, after combining both cohorts, we investigated whether these associations were modulated by disease status. All five combinations were significant (p < 10-3 for all tests). These combinations also showed significant evidence for an effect of the interaction between the mitochondrial and nuclear variants on disease risk. Their nuclear components mapped to TBCA, NIBAN3, and GLT25D1 and an uncharacterised intergenic region. In summary, starting from a single cohort design we identified combinations of nuclear and mitochondrial variants affecting PD disease risk.

In vivo mitochondrial base editing via adeno-associated viral delivery to mouse post-mitotic tissue

Mitochondria host key metabolic processes vital for cellular energy provision and are central to cell fate decisions. They are subjected to unique genetic control by both nuclear DNA and their own multi-copy genome - mitochondrial DNA (mtDNA). Mutations in mtDNA often lead to clinically heterogeneous, maternally inherited diseases that display different organ-specific presentation at any stage of life. For a long time, genetic manipulation of mammalian mtDNA has posed a major challenge, impeding our ability to understand the basic mitochondrial biology and mechanisms underpinning mitochondrial disease. However, an important new tool for mtDNA mutagenesis has emerged recently, namely double-stranded DNA deaminase (DddA)-derived cytosine base editor (DdCBE). Here, we test this emerging tool for in vivo use, by delivering DdCBEs into mouse heart using adeno-associated virus (AAV) vectors and show that it can install desired mtDNA edits in adult and neonatal mice. This work provides proof-of-concept for use of DdCBEs to mutagenize mtDNA in vivo in post-mitotic tissues and provides crucial insights into potential translation to human somatic gene correction therapies to treat primary mitochondrial disease phenotypes.

Mutations in mitochondrial DNA can lead to clinically heterogeneous disease. Here the authors demonstrate in vivo base editing of mouse mitochondrial DNA in a post-mitotic tissue by AAV delivery of DddA-derived cytosine base editor (DdCBE).

Comprehensive summary of mitochondrial DNA alterations in the postmortem human brain: A systematic review

BACKGROUND

Mitochondrial DNA (mtDNA) encodes 37 genes necessary for synthesizing 13 essential subunits of the oxidative phosphorylation system. mtDNA alterations are known to cause mitochondrial disease (MitD), a clinically heterogeneous group of disorders that often present with neuropsychiatric symptoms. Understanding the nature and frequency of mtDNA alterations in health and disease could be a cornerstone in disentangling the relationship between biochemical findings and clinical symptoms of brain disorders. This systematic review aimed to summarize the mtDNA alterations in human brain tissue reported to date that have implications for further research on the pathophysiological significance of mtDNA alterations in brain functioning.

METHODS

We searched the PubMed and Embase databases using distinct terms related to postmortem human brain and mtDNA up to June 10, 2021. Reports were eligible if they were empirical studies analysing mtDNA in postmortem human brains.

FINDINGS

A total of 158 of 637 studies fulfilled the inclusion criteria and were clustered into the following groups: MitD (48 entries), neurological diseases (NeuD, 55 entries), psychiatric diseases (PsyD, 15 entries), a miscellaneous group with controls and other clinical diseases (5 entries), ageing (20 entries), and technical issues (5 entries). Ten entries were ascribed to more than one group. Pathogenic single nucleotide variants (pSNVs), both homo- or heteroplasmic variants, have been widely reported in MitD, with heteroplasmy levels varying among brain regions; however, pSNVs are rarer in NeuD, PsyD and ageing. A lower mtDNA copy number (CN) in disease was described in most, but not all, of the identified studies. mtDNA deletions were identified in individuals in the four clinical categories and ageing. Notably, brain samples showed significantly more mtDNA deletions and at higher heteroplasmy percentages than blood samples, and several of the deletions present in the brain were not detected in the blood. Finally, mtDNA heteroplasmy, mtDNA CN and the deletion levels varied depending on the brain region studied.

INTERPRETATION

mtDNA alterations are well known to affect human tissues, including the brain. In general, we found that studies of MitD, NeuD, PsyD, and ageing were highly variable in terms of the type of disease or ageing process investigated, number of screened individuals, studied brain regions and technology used. In NeuD and PsyD, no particular type of mtDNA alteration could be unequivocally assigned to any specific disease or diagnostic group. However, the presence of mtDNA deletions and mtDNA CN variation imply a role for mtDNA in NeuD and PsyD. Heteroplasmy levels and threshold effects, affected brain regions, and mitotic segregation patterns of mtDNA alterations may be involved in the complex inheritance of NeuD and PsyD and in the ageing process. Therefore, more information is needed regarding the type of mtDNA alteration, the affected brain regions, the heteroplasmy levels, and their relationship with clinical phenotypes and the ageing process.

FUNDING

Hospital Universitari Institut Pere Mata; Institut d'Investigació Sanitària Pere Virgili; Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación (PI18/00514).

Human studies of mitochondrial biology demonstrate an overall lack of binary sex differences: A multivariate meta-analysis

Mitochondria are maternally inherited organelles that play critical tissue-specific roles, including hormone synthesis and energy production, that influence human development, health, and aging. However, whether mitochondria from women and men exhibit consistent biological differences remains unclear, representing a major gap in knowledge. This meta-analysis systematically examined four domains and six subdomains of mitochondrial biology (total 39 measures), including mitochondrial content, respiratory capacity, reactive oxygen species (ROS) production, morphometry, and mitochondrial DNA copy number. Standardized effect sizes (Hedge's g) of sex differences were computed for each measure using data in 2258 participants (51.5% women) from 50 studies. Only two measures demonstrated aggregate binary sex differences: higher mitochondrial content in women's WAT and isolated leukocyte subpopulations (g = 0.20, χ² p = .01), and higher ROS production in men's skeletal muscle (g = 0.49, χ² p < .0001). Sex differences showed weak to no correlation with age or BMI. Studies with small sample sizes tended to overestimate effect sizes (r = -.17, p < .001), and sex differences varied by tissue examined. Our findings point to a wide variability of findings in the literature concerning possible binary sex differences in mitochondrial biology. Studies specifically designed to capture sex- and gender-related differences in mitochondrial biology are needed, including detailed considerations of physical activity and sex hormones.

Mitochondrial DNA variation across 56,434 individuals in gnomAD

Genomic databases of allele frequency are extremely helpful for evaluating clinical variants of unknown significance; however, until now, databases such as the Genome Aggregation Database (gnomAD) have focused on nuclear DNA and have ignored the mitochondrial genome (mtDNA). Here we present a pipeline to call mtDNA variants that addresses three technical challenges: (i) detecting homoplasmic and heteroplasmic variants, present respectively in all or a fraction of mtDNA molecules, (ii) circular mtDNA genome, and (iii) misalignment of nuclear sequences of mitochondrial origin (NUMTs). We observed that mtDNA copy number per cell varied across gnomAD cohorts and influenced the fraction of NUMT-derived false-positive variant calls, which can account for the majority of putative heteroplasmies. To avoid false positives, we excluded contaminated samples, cell lines, and samples prone to NUMT misalignment due to few mtDNA copies. Furthermore, we report variants with heteroplasmy greater than 10%. We applied this pipeline to 56,434 whole genome sequences in the gnomAD v3.1 database that includes individuals of European (58%), African (25%), Latino (10%), and Asian (5%) ancestry. Our gnomAD v3.1 release contains population frequencies for 10,850 unique mtDNA variants at more than half of all mtDNA bases. We report frequencies within each nuclear ancestral population and mitochondrial haplogroup. Homoplasmic variants account for most variant calls (98%) and unique variants (85%). We observed that 1/250 individuals carry a pathogenic mtDNA variant with heteroplasmy above 10%. These mtDNA population allele frequencies are freely accessible and will aid in diagnostic interpretation and research studies.

Changes of Gut Microbiota by Natural mtDNA Variant Differences Augment Susceptibility to Metabolic Disease and Ageing

We recently reported on two mouse strains carrying different single nucleotide variations in the mitochondrial complex I gene, i.e., B6-mtBPL mice carrying m.11902T>C and B6-mtALR carrying m.4738C>A. B6-mtBPL mice exhibited a longer lifespan and a lower metabolic disease susceptibility despite mild mitochondrial functional differences in steady-state. As natural polymorphisms in the mitochondrial DNA (mtDNA) are known to be associated with distinct patterns of gut microbial composition, we further investigated the gut microbiota composition in these mice strains. In line with mouse phenotypes, we found a significantly lower abundance of Proteobacteria, which is positively associated with pathological conditions, in B6-mtBPL compared to B6-mtALR mice. A prediction of functional profile of significantly differential bacterial genera between these strains revealed an involvement of glucose metabolism pathways. Whole transcriptome analysis of liver samples from B6-mtBPL and B6-mtALR mice confirmed these findings. Thus, both host gene expression and gut microbial changes caused by the mtDNA variant differences may contribute to the ageing and metabolic phenotypes observed in these mice strains. Since gut microbiota are easier to modulate, compared with mtDNA variants, identification of such mtDNA variants, specific gut bacterial species and bacterial metabolites may be a potential intervention to modulate common diseases, which are differentially susceptible to individuals with different mtDNA variants.

From Forensics to Clinical Research: Expanding the Variant Calling Pipeline for the Precision ID mtDNA Whole Genome Panel

Despite a multitude of methods for the sample preparation, sequencing, and data analysis of mitochondrial DNA (mtDNA), the demand for innovation remains, particularly in comparison with nuclear DNA (nDNA) research. The Applied Biosystems™ Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific, USA) is an innovative library preparation kit suitable for degraded samples and low DNA input. However, its bioinformatic processing occurs in the enterprise Ion Torrent Suite™ Software (TSS), yielding BAM files aligned to an unorthodox version of the revised Cambridge Reference Sequence (rCRS), with a heteroplasmy threshold level of 10%. Here, we present an alternative customizable pipeline, the PrecisionCallerPipeline (PCP), for processing samples with the correct rCRS output after Ion Torrent sequencing with the Precision ID library kit. Using 18 samples (3 original samples and 15 mixtures) derived from the 1000 Genomes Project, we achieved overall improved performance metrics in comparison with the proprietary TSS, with optimal performance at a 2.5% heteroplasmy threshold. We further validated our findings with 50 samples from an ongoing independent cohort of stroke patients, with PCP finding 98.31% of TSS's variants (TSS found 57.92% of PCP's variants), with a significant correlation between the variant levels of variants found with both pipelines.

Gaining Insight into Mitochondrial Genetic Variation and Downstream Pathophysiology: What Can i(PSCs) Do?

Mitochondria are specialized organelles involved in energy production that have retained their own genome throughout evolutionary history. The mitochondrial genome (mtDNA) is maternally inherited and requires coordinated regulation with nuclear genes to produce functional enzyme complexes that drive energy production. Each mitochondrion contains 5-10 copies of mtDNA and consequently, each cell has several hundreds to thousands of mtDNAs. Due to the presence of multiple copies of mtDNA in a mitochondrion, mtDNAs with different variants may co-exist, a condition called heteroplasmy. Heteroplasmic variants can be clonally expanded, even in post-mitotic cells, as replication of mtDNA is not tied to the cell-division cycle. Heteroplasmic variants can also segregate during germ cell formation, underlying the inheritance of some mitochondrial mutations. Moreover, the uneven segregation of heteroplasmic variants is thought to underlie the heterogeneity of mitochondrial variation across adult tissues and resultant differences in the clinical presentation of mitochondrial disease. Until recently, however, the mechanisms mediating the relation between mitochondrial genetic variation and disease remained a mystery, largely due to difficulties in modeling human mitochondrial genetic variation and diseases. The advent of induced pluripotent stem cells (iPSCs) and targeted gene editing of the nuclear, and more recently mitochondrial, genomes now provides the ability to dissect how genetic variation in mitochondrial genes alter cellular function across a variety of human tissue types. This review will examine the origins of mitochondrial heteroplasmic variation and propagation, and the tools used to model mitochondrial genetic diseases. Additionally, we discuss how iPSC technologies represent an opportunity to advance our understanding of human mitochondrial genetics in disease.

Blood mitochondrial DNA copy number: What are we counting?

There is growing scientific interest to develop scalable biological measures that capture mitochondrial (dys)function. Mitochondria have their own genome, the mitochondrial DNA (mtDNA). It has been proposed that the number of mtDNA copies per cell (mtDNA copy number; mtDNAcn) reflects mitochondrial health. The common availability of stored DNA material or existing DNA sequencing data, especially from blood and other easy-to-collect samples, has made its quantification a popular approach in clinical and epidemiological studies. However, the interpretation of mtDNAcn is not univocal, and either a reduction or elevation in mtDNAcn can indicate dysfunction. The major determinants of blood-derived mtDNAcn are the heterogeneous cell type composition of leukocytes and platelet abundance, which can change with time of day, aging, and with disease. Hematopoiesis is a likely driver of blood mtDNAcn. Here we discuss the rationale and available methods to quantify mtDNAcn, the influence of blood cell type variations, and consider important gaps in knowledge that need to be resolved to maximize the scientific value around the investigation of blood mtDNAcn.

A globally diverse reference alignment and panel for imputation of mitochondrial DNA variants

BACKGROUND

Variation in mitochondrial DNA (mtDNA) identified by genotyping microarrays or by sequencing only the hypervariable regions of the genome may be insufficient to reliably assign mitochondrial genomes to phylogenetic lineages or haplogroups. This lack of resolution can limit functional and clinical interpretation of a substantial body of existing mtDNA data. To address this limitation, we developed and evaluated a large, curated reference alignment of complete mtDNA sequences as part of a pipeline for imputing missing mtDNA single nucleotide variants (mtSNVs). We call our reference alignment and pipeline MitoImpute.

RESULTS

We aligned the sequences of 36,960 complete human mitochondrial genomes downloaded from GenBank, filtered and controlled for quality. These sequences were reformatted for use in imputation software, IMPUTE2. We assessed the imputation accuracy of MitoImpute by measuring haplogroup and genotype concordance in data from the 1000 Genomes Project and the Alzheimer's Disease Neuroimaging Initiative (ADNI). The mean improvement of haplogroup assignment in the 1000 Genomes samples was 42.7% (Matthew's correlation coefficient = 0.64). In the ADNI cohort, we imputed missing single nucleotide variants.

CONCLUSION

These results show that our reference alignment and panel can be used to impute missing mtSNVs in existing data obtained from using microarrays, thereby broadening the scope of functional and clinical investigation of mtDNA. This improvement may be particularly useful in studies where participants have been recruited over time and mtDNA data obtained using different methods, enabling better integration of early data collected using less accurate methods with more recent sequence data.