mtDNA Variation and Analysis Using Mitomap and Mitomaster

mtDNA analysis using Mitopore

Mitochondrial DNA (mtDNA) analysis is crucial for the diagnosis of mitochondrial disorders, forensic investigations, and basic research. Existing pipelines are complex, expensive, and require specialized personnel. In many cases, including the diagnosis of detrimental single nucleotide variants (SNVs), mtDNA analysis is still carried out using Sanger sequencing. Here, we developed a simple workflow and a publicly available webserver named Mitopore that allows the detection of mtDNA SNVs, indels, and haplogroups. To simplify mtDNA analysis, we tailored our workflow to process noisy long-read sequencing data for mtDNA analysis, focusing on sequence alignment and parameter optimization. We implemented Mitopore with eliBQ (eliminate bad quality reads), an innovative quality enhancement that permits the increase of per-base quality of over 20% for low-quality data. The whole Mitopore workflow and webserver were validated using patient-derived and induced pluripotent stem cells harboring mtDNA mutations. Mitopore streamlines mtDNA analysis as an easy-to-use fast, reliable, and cost-effective analysis method for both long- and short-read sequencing data. This significantly enhances the accessibility of mtDNA analysis and reduces the cost per sample, contributing to the progress of mtDNA-related research and diagnosis.

Genetic susceptibility to optic neuropathy in patients with alcohol use disorder

BACKGROUND

The pathophysiology of toxico-nutritional optic neuropathies remains debated, with no clear understanding of the respective roles played by the direct alcohol toxicity, smoking and the often associated vitamin deficiencies, which are risk factors for optic neuropathy. Our aim was to investigate genetic susceptibility in patients with bilateral infraclinical optic neuropathy associated with chronic alcohol use disorder.

METHODS

This retrospective cohort study included 102 visually asymptomatic patients with documented alcohol use disorder from a French reference center. Optic neuropathy was identified with optical coherence tomography (OCT), after which genetic susceptibility in the group of affected patients was investigated. Genetic testing was performed using panel sequencing of 87 nuclear genes and complete mitochondrial DNA sequencing.

RESULTS

Optic neuropathy was detected in 36% (37/102) of the included patients. Genetic testing of affected patients disclosed two patients (2/30, 6.7%) with optic neuropathy associated with pathogenic variants affecting the SPG7 gene and five patients (5/30, 16.7%) who harbored variants of uncertain significance close to probable pathogenicity in the genes WFS1, LOXL1, MMP19, NR2F1 and PMPCA. No pathogenic mitochondrial DNA variants were found in this group.

CONCLUSIONS

OCT can detect presence of asymptomatic optic neuropathy in patients with chronic alcohol use disorder. Furthermore, genetic susceptibility to optic neuropathy in this setting is found in almost a quarter of affected patients. Further studies may clarify the role of preventative measures in patients who might be predisposed to avoidable visual loss and blindness.

Identification of a novel MT-ND3 variant and restoring mitochondrial function by allotopic expression of MT-ND3 gene

Mitochondrial diseases are caused by nuclear, or mitochondrial DNA (mtDNA) variants and related co-factors. Here, we report a novel m.10197G > C variant in MT-ND3 in a patient, and two other patients with m.10191 T > C. MT-ND3 variants are known to cause Leigh syndrome or mitochondrial complex I deficiency. We performed the functional analyses of the novel m.10197G > C variant that significantly lowered MT-ND3 protein levels, causing complex I assembly and activity deficiency, and reduction of ATP synthesis. We adapted a previously described re-engineering technique of delivering mitochondrial genes into mitochondria through codon optimization for nuclear expression and translation by cytoplasmic ribosomes to rescue defects arising from the MT-ND3 variants. We constructed mitochondrial targeting sequences along with the codon-optimized MT-ND3 and imported them into the mitochondria. To achieve the goal, we imported codon-optimized MT-ND3 into mitochondria in three patients with m.10197G > C and m.10191 T > C missense variants in the MT-ND3. Nuclear expression of the MT-ND3 gene partially restored protein levels, complex I deficiency, and significant improvement of ATP production indicating a functional rescue of the mutant phenotype. The codon-optimized nuclear expression of mitochondrial protein and import inside the mitochondria can supplement the requirements for ATP in energy-deficient mitochondrial disease patients.

A systematic review on the biochemical threshold of mitochondrial genetic variants

Mitochondrial DNA (mtDNA) variants cause a range of diseases from severe pediatric syndromes to aging-related conditions. The percentage of mtDNA copies carrying a pathogenic variant, variant allele frequency (VAF), must reach a threshold before a biochemical defect occurs, termed the biochemical threshold. Whether the often-cited biochemical threshold of >60% VAF is similar across mtDNA variants and cell types is unclear. In our systematic review, we sought to identify the biochemical threshold of mtDNA variants in relation to VAF by human tissue/cell type. We used controlled vocabulary terms to identify articles measuring oxidative phosphorylation (OXPHOS) complex activities in relation to VAF. We identified 76 eligible publications, describing 69, 12, 16, and 49 cases for complexes I, III, IV, and V, respectively. Few studies evaluated OXPHOS activities in diverse tissue types, likely reflective of clinical access. A number of cases with similar VAFs for the same pathogenic variant had varying degrees of residual activity of the affected complex, alluding to the presence of modifying variants. Tissues and cells with VAFs <60% associated with low complex activities were described, suggesting the possibility of a biochemical threshold of <60%. Using Kendall rank correlation tests, the VAF of the m.8993T > G variant correlated with complex V activity in skeletal muscle (τ = -0.58, P = 0.01, n = 13); however, no correlation was observed in fibroblasts (P = 0.7, n = 9). Our systematic review highlights the need to investigate the biochemical threshold over a wider range of VAFs in disease-relevant cell types to better define the biochemical threshold for specific mtDNA variants.

Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

Screening for Rare Mitochondrial Genome Variants Reveals a Potentially Novel Association between MT-CO1 and MT-TL2 Genes and Diabetes Phenotype

Variations in several nuclear genes predisposing humans to the development of MODY diabetes have been very well characterized by modern genetic diagnostics. However, recent reports indicate that variants in the mtDNA genome may also be associated with the diabetic phenotype. As relatively little research has addressed the entire mitochondrial genome in this regard, the aim of the present study is to evaluate the genetic variations present in mtDNA among individuals susceptible to MODY diabetes. In total, 193 patients with a MODY phenotype were tested with a custom panel with mtDNA enrichment. Heteroplasmic variants were selected for further analysis via further sequencing based on long-range PCR to evaluate the potential contribution of frequent NUMTs (acronym for nuclear mitochondrial DNA) insertions. Twelve extremely rare variants with a potential damaging character were selected, three of which were likely to be the result of NUMTs from the nuclear genome. The variant m.3243A>G in MT-TL1 was responsible for 3.5% of MODY cases in our study group. In addition, a novel, rare, and possibly pathogenic leucine variant m.12278T>C was found in MT-TL2. Our findings also found the MT-CO1 gene to be over-represented in the study group, with a clear phenotype-genotype correlation observed in one family. Our data suggest that heteroplasmic variants in MT-COI and MT-TL2 genes may play a role in the pathophysiology of glucose metabolism in humans.

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Light-strand bias and enriched zones of embedded ribonucleotides are associated with DNA replication and transcription in the human-mitochondrial genome

Abundant ribonucleoside-triphosphate (rNTP) incorporation into DNA by DNA polymerases in the form of ribonucleoside monophosphates (rNMPs) is a widespread phenomenon in nature, resulting in DNA-structural change and genome instability. The rNMP distribution, characteristics, hotspots and association with DNA metabolic processes in human mitochondrial DNA (hmtDNA) remain mostly unknown. Here, we utilize the ribose-seq technique to capture embedded rNMPs in hmtDNA of six different cell types. In most cell types, the rNMPs are preferentially embedded on the light strand of hmtDNA with a strong bias towards rCMPs; while in the liver-tissue cells, the rNMPs are predominately found on the heavy strand. We uncover common rNMP hotspots and conserved rNMP-enriched zones across the entire hmtDNA, including in the control region, which links the rNMP presence to the frequent hmtDNA replication-failure events. We show a strong correlation between coding-sequence size and rNMP-embedment frequency per nucleotide on the non-template, light strand in all cell types, supporting the presence of transient RNA-DNA hybrids preceding light-strand replication. Moreover, we detect rNMP-embedment patterns that are only partly conserved across the different cell types and are distinct from those found in yeast mtDNA. The study opens new research directions to understand the biology of hmtDNA and genomic rNMPs.

Children born after assisted reproduction more commonly carry a mitochondrial genotype associating with low birthweight

Children conceived through assisted reproductive technologies (ART) have an elevated risk of lower birthweight, yet the underlying cause remains unclear. Our study explores mitochondrial DNA (mtDNA) variants as contributors to birthweight differences by impacting mitochondrial function during prenatal development. We deep-sequenced the mtDNA of 451 ART and spontaneously conceived (SC) individuals, 157 mother-child pairs and 113 individual oocytes from either natural menstrual cycles or after ovarian stimulation (OS) and find that ART individuals carried a different mtDNA genotype than SC individuals, with more de novo non-synonymous variants. These variants, along with rRNA variants, correlate with lower birthweight percentiles, independent of conception mode. Their higher occurrence in ART individuals stems from de novo mutagenesis associated with maternal aging and OS-induced oocyte cohort size. Future research will establish the long-term health consequences of these changes and how these findings will impact the clinical practice and patient counselling in the future.

Bioanthropological analysis of human remains from the archaic and classic period discovered in Puyil cave, Mexico

OBJECTIVES

Determine the geographic place of origin and maternal lineage of prehistoric human skeletal remains discovered in Puyil Cave, Tabasco State, Mexico, located in a region currently populated by Olmec, Zoque and Maya populations.

MATERIALS AND METHODS

All specimens were radiocarbon (14 C) dated (beta analytic), had dental modifications classified, and had an analysis of 13 homologous reference points conducted to evaluate artificial cranial deformation (ACD). Following DNA purification, hypervariable region I (HVR-1) of the mitogenome was amplified and Sanger sequenced. Finally, Next Generation Sequencing (NGS) was performed for total DNA. Mitochondrial DNA (mtDNA) variants and haplogroups were determined using BioEdit 7.2 and IGV software and confirmed with MITOMASTER and WebHome softwares.

RESULTS

Radiocarbon dating (14 C) demonstrated that the inhabitants of Puyil Cave lived during the Archaic and Classic Periods and displayed tabular oblique and tabular mimetic ACD. These pre-Hispanic remains exhibited five mtDNA lineages: A, A2, C1, C1c and D4. Network analysis revealed a close genetic affinity between pre-Hispanic Puyil Cave inhabitants and contemporary Maya subpopulations from Mexico and Guatemala, as well as individuals from Bolivia, Brazil, the Dominican Republic, and China.

CONCLUSIONS

Our results elucidate the dispersal of pre-Hispanic Olmec and Maya ancestors and suggest that ACD practices are closely related to Olmec and Maya practices. Additionally, we conclude that ACD has likely been practiced in the region since the Middle-Archaic Period.

SkQ3 Exhibits the Most Pronounced Antioxidant Effect on Isolated Rat Liver Mitochondria and Yeast Cells

Oxidative stress is involved in a wide range of age-related diseases. A critical role has been proposed for mitochondrial oxidative stress in initiating or promoting these pathologies and the potential for mitochondria-targeted antioxidants to fight them, making their search and testing a very urgent task. In this study, the mitochondria-targeted antioxidants SkQ1, SkQ3 and MitoQ were examined as they affected isolated rat liver mitochondria and yeast cells, comparing SkQ3 with clinically tested SkQ1 and MitoQ. At low concentrations, all three substances stimulated the oxidation of respiratory substrates in state 4 respiration (no ADP addition); at higher concentrations, they inhibited the ADP-triggered state 3 respiration and the uncoupled state, depolarized the inner mitochondrial membrane, contributed to the opening of the mPTP (mitochondrial permeability transition pore), did not specifically affect ATP synthase, and had a pronounced antioxidant effect. SkQ3 was the most active antioxidant, not possessing, unlike SkQ1 or MitoQ, prooxidant activity with increasing concentrations. In yeast cells, all three substances reduced prooxidant-induced intracellular oxidative stress and cell death and prevented and reversed mitochondrial fragmentation, with SkQ3 being the most efficient. These data allow us to consider SkQ3 as a promising potential therapeutic agent to mitigate pathologies associated with oxidative stress.

VarCards2: an integrated genetic and clinical database for ACMG-AMP variant-interpretation guidelines in the human whole genome

VarCards, an online database, combines comprehensive variant- and gene-level annotation data to streamline genetic counselling for coding variants. Recognising the increasing clinical relevance of non-coding variations, there has been an accelerated development of bioinformatics tools dedicated to interpreting non-coding variations, including single-nucleotide variants and copy number variations. Regrettably, most tools remain as either locally installed databases or command-line tools dispersed across diverse online platforms. Such a landscape poses inconveniences and challenges for genetic counsellors seeking to utilise these resources without advanced bioinformatics expertise. Consequently, we developed VarCards2, which incorporates nearly nine billion artificially generated single-nucleotide variants (including those from mitochondrial DNA) and compiles vital annotation information for genetic counselling based on ACMG-AMP variant-interpretation guidelines. These annotations include (I) functional effects; (II) minor allele frequencies; (III) comprehensive function and pathogenicity predictions covering all potential variants, such as non-synonymous substitutions, non-canonical splicing variants, and non-coding variations and (IV) gene-level information. Furthermore, VarCards2 incorporates 368 820 266 documented short insertions and deletions and 2 773 555 documented copy number variations, complemented by their corresponding annotation and prediction tools. In conclusion, VarCards2, by integrating over 150 variant- and gene-level annotation sources, significantly enhances the efficiency of genetic counselling and can be freely accessed at http://www.genemed.tech/varcards2/.

The clinical, myopathological, and genetic analysis of 155 Chinese mitochondrial ophthalmoplegia patients with mitochondrial DNA single large deletions

BACKGROUND

Progressive external ophthalmoplegia (PEO) is a common subtype of mitochondrial encephalomyopathy.

OBJECTIVE

The study aimed to investigate the relationship between mitochondrial DNA (mtDNA) abnormalities, muscle pathology, and clinical manifestations in Chinese patients with single large-scale mtDNA deletion presenting with PEO.

METHODS

This is a retrospective single-center study. Patients with PEO who had a single large deletion in mitochondrial DNA were included in this study. The associations were analyzed between mtDNA deletion patterns, myopathological changes, and clinical characteristics.

RESULTS

In total, 155 patients with mitochondrial PEO carrying single large-scale mtDNA mutations were enrolled, including 137 chronic progressive external ophthalmoplegia (CPEO) and 18 Kearns-Sayre syndrome (KSS) patients. The onset ages were 9.61 ± 4.12 in KSS and 20.15 ± 9.06 in CPEO. The mtDNA deletions ranged from 2225 bp to 9131 bp, with m.8470_13446del being the most common. The KSS group showed longer deletions than the CPEO group (p = 0.004). Additionally, a higher number of deleted genes encoding respiratory chain complex subunits (p = 0.001) and tRNA genes (p = 0.009) were also observed in the KSS group. A weak negative correlation between the mtDNA deletion size and ages of onset (p < 0.001, r = -0.369) was observed. The proportion of ragged red fibers, ragged blue fibers, and cytochrome c negative fibers did not correlate significantly with onset ages (p > 0.05). However, a higher percentage of abnormal muscle fibers corresponds to an increased prevalence of exercise intolerance, limb muscle weakness, dysphagia, and cerebellar ataxia.

CONCLUSION

We reported a large Chinese cohort consisting of mitochondrial PEO patients with single large-scale mtDNA deletions. Our results demonstrated that the length and locations of mtDNA deletions may influence onset ages and clinical phenotypes. The severity of muscle pathology could not only indicate diagnosis but also may be associated with clinical manifestations beyond the extraocular muscles.

Post-Transcriptional Methylation of Mitochondrial-tRNA Differentially Contributes to Mitochondrial Pathology

Human mitochondrial tRNAs (mt-tRNAs), critical for mitochondrial biogenesis, are frequently associated with pathogenic mutations. These mt-tRNAs have unusual sequence motifs and require post-transcriptional modifications to stabilize their fragile structures. However, whether a modification that stabilizes a wild-type (WT) mt-tRNA structure would also stabilize its pathogenic variants is unknown. Here we show that the N 1 -methylation of guanosine at position 9 (m 1 G9) of mt-Leu(UAA), while stabilizing the WT tRNA, has an opposite and destabilizing effect on variants associated with MELAS (mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes). This differential effect is further demonstrated by the observation that demethylation of m 1 G9, while damaging to the WT tRNA, is beneficial to the major pathogenic variant, improving its structure and activity. These results have new therapeutic implications, suggesting that the N 1 -methylation of mt-tRNAs at position 9 is a determinant of pathogenicity and that controlling the methylation level is an important modulator of mt-tRNA-associated diseases.

mtDNA Single-Nucleotide Variants Associated with Type 2 Diabetes

Type 2 diabetes (T2D) is a chronic systemic disease with a complex etiology, characterized by insulin resistance and mitochondrial dysfunction in various cell tissues. To explore this relationship, we conducted a secondary analysis of complete mtDNA sequences from 1261 T2D patients and 1105 control individuals. Our findings revealed significant associations between certain single-nucleotide polymorphisms (SNPs) and T2D. Notably, the variants m.1438A>G (rs2001030) (controls: 32 [27.6%], T2D: 84 [72.4%]; OR: 2.46; 95%CI: 1.64-3.78; p < 0.001), m.14766C>T (rs193302980) (controls: 498 [36.9%], T2D: 853 [63.1%]; OR: 2.57, 95%CI: 2.18-3.04, p < 0.001), and m.16519T>C (rs3937033) (controls: 363 [43.4%], T2D: 474 [56.6%]; OR: 1.24, 95%CI: 1.05-1.47, p = 0.012) were significantly associated with the likelihood of developing diabetes. The variant m.16189T>C (rs28693675), which has been previously documented in several studies across diverse populations, showed no association with T2D in our analysis (controls: 148 [13.39] T2D: 171 [13.56%]; OR: 1.03; 95%CI: 0.815-1.31; p = 0.83). These results provide evidence suggesting a link between specific mtDNA polymorphisms and T2D, possibly related to association rules, topological patterns, and three-dimensional conformations associated with regions where changes occur, rather than specific point mutations in the sequence.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia

The analysis of mitochondrial DNA (mtDNA) hypervariable region (HVR) sequence data from ancient human remains provides valuable insights into the genetic structure and population dynamics of ancient populations. mtDNA is particularly useful in studying ancient populations, because it is maternally inherited and has a higher mutation rate compared to nuclear DNA. To determine the genetic structure of three Colombian pre-Hispanic populations and compare them with current populations, we determined the haplotypes from human bone remains by sequencing several mitochondrial DNA segments. A wide variety of mitochondrial polymorphisms were obtained from 33 samples. Our results support a high population heterogeneity among pre-Hispanic populations in Colombia.

Autosomal recessive Leber's hereditary optic neuropathy caused by a homozygous variant in DNAJC30 gene

Recently, Stenton et al. (2021) described a new, autosomal recessive inheritance pattern of Leber's hereditary optic neuropathy (LHON) caused by missense variants in the DNAJC30 gene. The DNAJC30 c.152A > G, p.(Tyr51Cys) variant was by far the most common variant reported in patients originating from Eastern Europe, therefore, it is believed to be a founder variant in these populations. We report the first two cases of DNAJC30-linked autosomal recessive LHON in a young male and a female originating from Estonia. The patients presented severe loss of central vision and clinical features indistinguishable from mitochondrial LHON. The whole exome sequencing carried out in the male patient and the next-generation sequencing panel in the young female patient identified the same homozygous missense variant in the DNAJC30 gene. Our cases further reinforce the pathogenicity of c.152A > G, p.(Tyr51Cys) DNAJC30 variant causing autosomal recessive LHON. According to the gnomAD database, the allele frequency of this variant in the Estonian population is 0.8%, translating into a prevalence of carriers of 1:60. It is the highest among different gnomAD populations. Applying the Hardy-Weinberg equation, an estimated 92 persons in the Estonian population carry the homozygous variant c.152A > G, p.(Tyr51Cys) in DNAJC30. In patients with LHON, we advise sequencing both the DNAJC30 gene and mitochondrial DNA simultaneously.

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

Mitochondrial Genome Variation in Polish Elite Athletes

Energy efficiency is one of the fundamental athletic performance-affecting features of the cell and the organism as a whole. Mitochondrial DNA (mtDNA) variants and haplogroups have been linked to the successful practice of various sports, but despite numerous studies, understanding of the correlation is far from being comprehensive. In this study, the mtDNA sequence and copy number were determined for 99 outstanding Polish male athletes performing in power (n = 52) or endurance sports (n = 47) and 100 controls. The distribution of haplogroups, single nucleotide variant association, heteroplasmy, and mtDNA copy number were analyzed in the blood and saliva. We found no correlation between any haplogroup, single nucleotide variant, especially rare or non-synonymous ones, and athletic performance. Interestingly, heteroplasmy was less frequent in the study group, especially in endurance athletes. We observed a lower mtDNA copy number in both power and endurance athletes compared to controls. This could result from an inactivity of compensatory mechanisms activated by disadvantageous variants present in the general population and indicates a favorable genetic makeup of the athletes. The results emphasize a need for a more comprehensive analysis of the involvement of the mitochondrial genome in physical performance, combining nucleotide and copy number analysis in the context of nuclear gene variants.

APOGEE 2: multi-layer machine-learning model for the interpretable prediction of mitochondrial missense variants

Mitochondrial dysfunction has pleiotropic effects and is frequently caused by mitochondrial DNA mutations. However, factors such as significant variability in clinical manifestations make interpreting the pathogenicity of variants in the mitochondrial genome challenging. Here, we present APOGEE 2, a mitochondrially-centered ensemble method designed to improve the accuracy of pathogenicity predictions for interpreting missense mitochondrial variants. Built on the joint consensus recommendations by the American College of Medical Genetics and Genomics/Association for Molecular Pathology, APOGEE 2 features an improved machine learning method and a curated training set for enhanced performance metrics. It offers region-wise assessments of genome fragility and mechanistic analyses of specific amino acids that cause perceptible long-range effects on protein structure. With clinical and research use in mind, APOGEE 2 scores and pathogenicity probabilities are precompiled and available in MitImpact. APOGEE 2's ability to address challenges in interpreting mitochondrial missense variants makes it an essential tool in the field of mitochondrial genetics.

Restoration of mitochondrial function through activation of hypomodified tRNAs with pathogenic mutations associated with mitochondrial diseases

Mutations in mitochondrial (mt-)tRNAs frequently cause mitochondrial dysfunction. Mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS), and myoclonus epilepsy associated with ragged red fibers (MERRF) are major clinical subgroups of mitochondrial diseases caused by pathogenic point mutations in tRNA genes encoded in mtDNA. We previously reported a severe reduction in the frequency of 5-taurinomethyluridine (τm5U) and its 2-thiouridine derivative (τm5s2U) in the anticodons of mutant mt-tRNAs isolated from the cells of patients with MELAS and MERRF, respectively. The hypomodified tRNAs fail to decode cognate codons efficiently, resulting in defective translation of respiratory chain proteins in mitochondria. To restore the mitochondrial activity of MELAS patient cells, we overexpressed MTO1, a τm5U-modifying enzyme, in patient-derived myoblasts. We used a newly developed primer extension method and showed that MTO1 overexpression almost completely restored the τm5U modification of the MELAS mutant mt-tRNALeu(UUR). An increase in mitochondrial protein synthesis and oxygen consumption rate suggested that the mitochondrial function of MELAS patient cells can be activated by restoring the τm5U of the mutant tRNA. In addition, we confirmed that MTO1 expression restored the τm5s2U of the mutant mt-tRNALys in MERRF patient cells. These findings pave the way for epitranscriptomic therapies for mitochondrial diseases.

Mitochondrial DNA mutations can influence the post-implantation development of human mosaic embryos

Introduction: Several healthy euploid births have been reported following the transfer of mosaic embryos, including both euploid and aneuploid blastomeres. This has been attributed to a reduced number of aneuploid cells, as previously reported in mice, but remains poorly explored in humans. We hypothesized that mitochondrial function, one of the most critical factors for embryonic development, can influence human post-implantation embryonic development, including a decrease of aneuploid cells in mosaic embryos. Methods: To clarify the role of mitochondrial function, we biopsied multiple parts of each human embryo and observed the remaining embryos under in vitro culture as a model of post-implantation development (n = 27 embryos). Karyotyping, whole mitochondrial DNA (mtDNA) sequencing, and mtDNA copy number assays were performed on all pre- and post-culture samples. Results: The ratio of euploid embryos was significantly enhanced during in vitro culture, whereas the ratio of mosaic embryos was significantly reduced. Furthermore, post-culture euploid and culturable embryos had significantly few mtDNA mutations, although mtDNA copy numbers did not differ. Discussion: Our results indicate that aneuploid cells decrease in human embryos post-implantation, and mtDNA mutations might induce low mitochondrial function and influence the development of post-implantation embryos with not only aneuploidy but also euploidy. Analyzing the whole mtDNA mutation number may be a novel method for selecting a better mosaic embryo for transfer.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

Multimodal single-cell analysis of nonrandom heteroplasmy distribution in human retinal mitochondrial disease

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multisystem disease. The wide range of manifestations observed in patients with mitochondrial disease results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify nonrandom distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA-Seq, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profiled the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and from healthy control donors. Utilizing the retina as a model for complex multilineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm-derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm-derived lineage, namely the vasculature of the choroid, was near homoplasmic for the WT allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the nonrandom nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

Single-cell multi-omics of mitochondrial DNA disorders reveals dynamics of purifying selection across human immune cells

Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.

Modeling of mitochondrial genetic polymorphisms reveals induction of heteroplasmy by pleiotropic disease locus 10398A>G

Mitochondrial (MT) dysfunction has been associated with several neurodegenerative diseases including Alzheimer's disease (AD). While MT-copy number differences have been implicated in AD, the effect of MT heteroplasmy on AD has not been well characterized. Here, we analyzed over 1800 whole genome sequencing data from four AD cohorts in seven different tissue types to determine the extent of MT heteroplasmy present. While MT heteroplasmy was present throughout the entire MT genome for blood samples, we detected MT heteroplasmy only within the MT control region for brain samples. We observed that an MT variant 10398A>G (rs2853826) was significantly associated with overall MT heteroplasmy in brain tissue while also being linked with the largest number of distinct disease phenotypes of all annotated MT variants in MitoMap. Using gene-expression data from our brain samples, our modeling discovered several gene networks involved in mitochondrial respiratory chain and Complex I function associated with 10398A>G. The variant was also found to be an expression quantitative trait loci (eQTL) for the gene MT-ND3. We further characterized the effect of 10398A>G by phenotyping a population of lymphoblastoid cell-lines (LCLs) with and without the variant allele. Examination of RNA sequence data from these LCLs reveal that 10398A>G was an eQTL for MT-ND4. We also observed in LCLs that 10398A>G was significantly associated with overall MT heteroplasmy within the MT control region, confirming the initial findings observed in post-mortem brain tissue. These results provide novel evidence linking MT SNPs with MT heteroplasmy and open novel avenues for the investigation of pathomechanisms that are driven by this pleiotropic disease associated loci.

Clinical and bi-genomic DNA findings of patients suspected to have mitochondrial diseases

Background: Mitochondrial diseases are the most common group of inherited metabolic disorders, causing difficulties in definite diagnosis due to clinical and genetic heterogeneity. Clinical components are predominantly associated with pathogenic variants shown in nuclear or mitochondrial genomes that affect vital respiratory chain function. The development of high-throughput sequencing technologies has accelerated the elucidation of the genetic etiology of many genetic diseases that previously remained undiagnosed. Methods: Thirty affected patients from 24 unrelated families with clinical, radiological, biochemical, and histopathological evaluations considered for mitochondrial diseases were investigated. DNA isolated from the peripheral blood samples of probands was sequenced for nuclear exome and mitochondrial DNA (mtDNA) analyses. MtDNA sequencing was also performed from the muscle biopsy material in one patient. For segregation, Sanger sequencing is performed for pathogenic alterations in five other affected family members and healthy parents. Results: Exome sequencing revealed 14 different pathogenic variants in nine genes encoding mitochondrial function peptides (AARS2, EARS2, ECHS1, FBXL4, MICOS13, NDUFAF6, OXCT1, POLG, and TK2) in 12 patients from nine families and four variants in genes encoding important for muscle structure (CAPN3, DYSF, and TCAP) in six patients from four families. Three probands carried pathogenic mtDNA variations in two genes (MT-ATP6 and MT-TL1). Nine variants in five genes are reported for the first time with disease association: (AARS2: c.277C>T/p.(R93*), c.845C>G/p.(S282C); EARS2: c.319C>T/p.(R107C), c.1283delC/p.(P428Lfs*); ECHS1: c.161G>A/p.(R54His); c.202G>A/p.(E68Lys); NDUFAF6: c.479delA/p.(N162Ifs*27); and OXCT1: c.1370C>T/p.(T457I), c.1173-139G>T/p.(?). Conclusion: Bi-genomic DNA sequencing clarified genetic etiology in 67% (16/24) of the families. Diagnostic utility by mtDNA sequencing in 13% (3/24) and exome sequencing in 54% (13/24) of the families prioritized searching for nuclear genome pathologies for the first-tier test. Weakness and muscle wasting observed in 17% (4/24) of the families underlined that limb-girdle muscular dystrophy, similar to mitochondrial myopathy, is an essential point for differential diagnosis. The correct diagnosis is crucial for comprehensive genetic counseling of families. Also, it contributes to making treatment-helpful referrals, such as ensuring early access to medication for patients with mutations in the TK2 gene.

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.

Nanopore sequencing identifies a higher frequency and expanded spectrum of mitochondrial DNA deletion mutations in human aging

Mitochondrial DNA (mtDNA) deletion mutations cause many human diseases and are linked to age-induced mitochondrial dysfunction. Mapping the mutation spectrum and quantifying mtDNA deletion mutation frequency is challenging with next-generation sequencing methods. We hypothesized that long-read sequencing of human mtDNA across the lifespan would detect a broader spectrum of mtDNA rearrangements and provide a more accurate measurement of their frequency. We employed nanopore Cas9-targeted sequencing (nCATS) to map and quantitate mtDNA deletion mutations and develop analyses that are fit-for-purpose. We analyzed total DNA from vastus lateralis muscle in 15 males ranging from 20 to 81 years of age and substantia nigra from three 20-year-old and three 79-year-old men. We found that mtDNA deletion mutations detected by nCATS increased exponentially with age and mapped to a wider region of the mitochondrial genome than previously reported. Using simulated data, we observed that large deletions are often reported as chimeric alignments. To address this, we developed two algorithms for deletion identification which yield consistent deletion mapping and identify both previously reported and novel mtDNA deletion breakpoints. The identified mtDNA deletion frequency measured by nCATS correlates strongly with chronological age and predicts the deletion frequency as measured by digital PCR approaches. In substantia nigra, we observed a similar frequency of age-related mtDNA deletions to those observed in muscle samples, but noted a distinct spectrum of deletion breakpoints. NCATS-mtDNA sequencing allows the identification of mtDNA deletions on a single-molecule level, characterizing the strong relationship between mtDNA deletion frequency and chronological aging.

Case report: A rare variant m.4135T>C in the MT-ND1 gene leads to Leber hereditary optic neuropathy and altered respiratory chain supercomplexes

Leber hereditary optic neuropathy is a primary mitochondrial disease characterized by acute visual loss due to the degeneration of retinal ganglion cells. In this study, we describe a patient carrying a rare missense heteroplasmic variant in MT-ND1, NC_012920.1:m.4135T>C (p.Tyr277His) manifesting with a typical bilateral painless decrease of the visual function, triggered by physical exercise or higher ambient temperature. Functional studies in muscle and fibroblasts show that amino acid substitution Tyr277 with His leads to only a negligibly decreased level of respiratory chain complex I (CI), but the formation of supercomplexes and the activity of the enzyme are disturbed noticeably. Our data indicate that although CI is successfully assembled in the patient's mitochondria, its function is hampered by the m.4135T>C variant, probably by stabilizing CI in its inactive form. We conclude that the m.4135T>C variant together with a combination of external factors is necessary to manifest the phenotype.

Fish Models for Exploring Mitochondrial Dysfunction Affecting Neurodegenerative Disorders

Neurodegenerative disorders are characterized by the progressive loss of neuronal structure or function, resulting in memory loss and movement disorders. Although the detailed pathogenic mechanism has not been elucidated, it is thought to be related to the loss of mitochondrial function in the process of aging. Animal models that mimic the pathology of a disease are essential for understanding human diseases. In recent years, small fish have become ideal vertebrate models for human disease due to their high genetic and histological homology to humans, ease of in vivo imaging, and ease of genetic manipulation. In this review, we first outline the impact of mitochondrial dysfunction on the progression of neurodegenerative diseases. Then, we highlight the advantages of small fish as model organisms, and present examples of previous studies regarding mitochondria-related neuronal disorders. Lastly, we discuss the applicability of the turquoise killifish, a unique model for aging research, as a model for neurodegenerative diseases. Small fish models are expected to advance our understanding of the mitochondrial function in vivo, the pathogenesis of neurodegenerative diseases, and be important tools for developing therapies to treat diseases.

Role of SIRT3 in mitochondrial biology and its therapeutic implications in neurodegenerative disorders

Sirtuin 3 (SIRT3), a mitochondrial deacetylase expressed preferentially in high-metabolic-demand tissues including the brain, requires NAD⁺ as a cofactor for catalytic activity. It regulates various processes such as energy homeostasis, redox balance, mitochondrial quality control, mitochondrial unfolded protein response (UPRmt), biogenesis, dynamics and mitophagy by altering protein acetylation status. Reduced SIRT3 expression or activity causes hyperacetylation of hundreds of mitochondrial proteins, which has been linked with neurological abnormalities, neuro-excitotoxicity and neuronal cell death. A body of evidence has suggested, SIRT3 activation as a potential therapeutic modality for age-related brain abnormalities and neurodegenerative disorders.

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.

The emergent role of mitochondrial RNA modifications in metabolic alterations

Mitochondrial epitranscriptomics refers to the modifications occurring in all the different RNA types of mitochondria. Although the number of mitochondrial RNA modifications is less than those in cytoplasm, substantial evidence indicates that they play a critical role in accurate protein synthesis. Recent evidence supported those modifications in mitochondrial RNAs also have crucial implications in mitochondrial-related diseases. In the light of current knowledge about the involvement, the association between mitochondrial RNA modifications and diseases arises from studies focusing on mutations in both mitochondrial and nuclear DNA genes encoding enzymes involved in such modifications. Here, we review the current evidence available for mitochondrial RNA modifications and their role in metabolic disorders, and we also explore the possibility of using them as promising targets for prevention and early detection. Finally, we discuss future directions of mitochondrial epitranscriptomics in these metabolic alterations, and how these RNA modifications may offer a new diagnostic and theragnostic avenue for preventive purposes. This article is categorized under: RNA Processing > RNA Editing and Modification.

Detection of UV-Induced Deletions in Mitochondrial DNA

Mitochondrial DNA (mtDNA) mutations are found in several human pathologies and are associated with aging. Deletion mutations in mtDNA result in the loss of essential genes for mitochondrial function. Over 250 deletion mutations have been reported and the common deletion is the most frequent mtDNA deletion linked to disease. This deletion removes 4977 base pairs of mtDNA. It has previously been shown that exposure to UVA radiation can promote the formation of the common deletion. Furthermore, aberrations in mtDNA replication and repair are associated with formation of the common deletion. However, molecular mechanisms describing the formation of this deletion are poorly characterized. This chapter describes a method to irradiate human skin fibroblasts with physiological doses of UVA and the subsequent detection of the common deletion by quantitative PCR analysis.

The tRNA identity landscape for aminoacylation and beyond

tRNAs are key partners in ribosome-dependent protein synthesis. This process is highly dependent on the fidelity of tRNA aminoacylation by aminoacyl-tRNA synthetases and relies primarily on sets of identities within tRNA molecules composed of determinants and antideterminants preventing mischarging by non-cognate synthetases. Such identity sets were discovered in the tRNAs of a few model organisms, and their properties were generalized as universal identity rules. Since then, the panel of identity elements governing the accuracy of tRNA aminoacylation has expanded considerably, but the increasing number of reported functional idiosyncrasies has led to some confusion. In parallel, the description of other processes involving tRNAs, often well beyond aminoacylation, has progressed considerably, greatly expanding their interactome and uncovering multiple novel identities on the same tRNA molecule. This review highlights key findings on the mechanistics and evolution of tRNA and tRNA-like identities. In addition, new methods and their results for searching sets of multiple identities on a single tRNA are discussed. Taken together, this knowledge shows that a comprehensive understanding of the functional role of individual and collective nucleotide identity sets in tRNA molecules is needed for medical, biotechnological and other applications.

Spectrum of germline and somatic mitochondrial DNA variants in Tuberous Sclerosis Complex

Tuberous Sclerosis Complex (TSC) is caused by loss of function variants in either TSC1 or TSC2 and is characterized by broad phenotypic heterogeneity. Currently, there is limited knowledge regarding the role of the mitochondrial genome (mtDNA) in TSC pathogenesis. In this study, we aimed to determine the prevalence and spectrum of germline and somatic mtDNA variants in TSC and identify potential disease modifiers. Analysis of mtDNA amplicon massively parallel sequencing (aMPS) data, off-target mtDNA from whole-exome sequencing (WES), and/or qPCR, revealed mtDNA alterations in 270 diverse tissues (139 TSC-associated tumors and 131 normal tissue samples) from 199 patients and six healthy individuals. Correlation of clinical features to mtDNA variants and haplogroup analysis was done in 102 buccal swabs (age: 20-71 years). No correlation was found between clinical features and either mtDNA variants or haplogroups. No pathogenic variants were identified in the buccal swab samples. Using in silico analysis, we identified three predicted pathogenic variants in tumor samples: MT-ND4 (m.11742G>A, p. Cys328Tyr, VAF: 43%, kidney angiomyolipoma), MT-CYB (m.14775T>C, p. Leu10Pro, VAF: 43%, LAM abdominal tumor) and MT-CYB (m.15555C>T, p. Pro270Leu, VAF: 7%, renal cell carcinoma). Large deletions of the mitochondrial genome were not detected. Analysis of tumors from 23 patients with corresponding normal tissue did not reveal any recurrent tumor-associated somatic variants. The mtDNA/gDNA ratio between tumors and corresponding normal tissue was also unchanged. Overall, our findings demonstrate that the mitochondrial genome is highly stable across tissues and within TSC-associated tumors.

Modopathies Caused by Mutations in Genes Encoding for Mitochondrial RNA Modifying Enzymes: Molecular Mechanisms and Yeast Disease Models

In eukaryotes, mitochondrial RNAs (mt-tRNAs and mt-rRNAs) are subject to specific nucleotide modifications, which are critical for distinct functions linked to the synthesis of mitochondrial proteins encoded by mitochondrial genes, and thus for oxidative phosphorylation. In recent years, mutations in genes encoding for mt-RNAs modifying enzymes have been identified as being causative of primary mitochondrial diseases, which have been called modopathies. These latter pathologies can be caused by mutations in genes involved in the modification either of tRNAs or of rRNAs, resulting in the absence of/decrease in a specific nucleotide modification and thus on the impairment of the efficiency or the accuracy of the mitochondrial protein synthesis. Most of these mutations are sporadic or private, thus it is fundamental that their pathogenicity is confirmed through the use of a model system. This review will focus on the activity of genes that, when mutated, are associated with modopathies, on the molecular mechanisms through which the enzymes introduce the nucleotide modifications, on the pathological phenotypes associated with mutations in these genes and on the contribution of the yeast Saccharomyces cerevisiae to confirming the pathogenicity of novel mutations and, in some cases, for defining the molecular defects.

Neurogenetic motor disorders

Advances in the field of neurogenetics have practical applications in rapid diagnosis on blood and body fluids to extract DNA, obviating the need for invasive investigations. The ability to obtain a presymptomatic diagnosis through genetic screening and biomarkers can be a guide to life-saving disease-modifying therapy or enzyme replacement therapy to compensate for the deficient disease-causing enzyme. The benefits of a comprehensive neurogenetic evaluation extend to family members in whom identification of the causal gene defect ensures carrier detection and at-risk counseling for future generations. This chapter explores the many facets of the neurogenetic evaluation in adult and pediatric motor disorders as a primer for later chapters in this volume and a roadmap for the future applications of genetics in neurology.

Compact zinc finger base editors that edit mitochondrial or nuclear DNA in vitro and in vivo

DddA-derived cytosine base editors (DdCBEs) use programmable DNA-binding TALE repeat arrays, rather than CRISPR proteins, a split double-stranded DNA cytidine deaminase (DddA), and a uracil glycosylase inhibitor to mediate C•G-to-T•A editing in nuclear and organelle DNA. Here we report the development of zinc finger DdCBEs (ZF-DdCBEs) and the improvement of their editing performance through engineering their architectures, defining improved ZF scaffolds, and installing DddA activity-enhancing mutations. We engineer variants with improved DNA specificity by integrating four strategies to reduce off-target editing. We use optimized ZF-DdCBEs to install or correct disease-associated mutations in mitochondria and in the nucleus. Leveraging their small size, we use a single AAV9 to deliver into heart, liver, and skeletal muscle in post-natal mice ZF-DdCBEs that efficiently install disease-associated mutations. While off-target editing of ZF-DdCBEs is likely too high for therapeutic applications, these findings demonstrate a compact, all-protein base editing research tool for precise editing of organelle or nuclear DNA without double-strand DNA breaks.

Analysis of Human Clinical Mutations of Mitochondrial ND1 in a Bacterial Model System for Complex I

The most common causes of mitochondrial dysfunction and disease include mutations in subunits and assembly factors of Complex I. Numerous mutations in the mitochondrial gene ND1 have been identified in humans. Currently, a bacterial model system provides the only method for rapid construction and analysis of mutations in homologs of human ND1. In this report, we have identified nine mutations in human ND1 that are reported to be pathogenic and are located at subunit interfaces. Our hypothesis was that these mutations would disrupt Complex I assembly. Seventeen mutations were constructed in the homologous nuoH gene in an E. coli model system. In addition to the clinical mutations, alanine substitutions were constructed in order to distinguish between a deleterious effect from the introduction of the mutant residue and the loss of the original residue. The mutations were moved to an expression vector containing all thirteen genes of the E. coli nuo operon coding for Complex I. Membrane vesicles were prepared and rates of deamino-NADH oxidase activity and proton translocation were measured. Samples were also tested for assembly by native gel electrophoresis and for expression of NuoH by immunoblotting. A range of outcomes was observed: Mutations at four of the sites allow normal assembly with moderate activity (50−76% of wild type). Mutations at the other sites disrupt assembly and/or activity, and in some cases the outcomes depend upon the amino acid introduced. In general, the outcomes are consistent with the proposed pathogenicity in humans.

Genomic architecture of autism from comprehensive whole-genome sequence annotation

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Combination of common mtDNA variants results in mitochondrial dysfunction and a connective tissue dysregulation

Mitochondrial dysfunction can be associated with a range of clinical manifestations. Here, we report a family with a complex phenotype including combinations of connective tissue, neurological, and metabolic symptoms that were passed on to all surviving children. Analysis of the maternally inherited mtDNA revealed a novel genotype encompassing the haplogroup J - defining mitochondrial DNA (mtDNA) ND5 m.13708G>A (A458T) variant arising on the mtDNA haplogroup H7A background, an extremely rare combination. Analysis of transmitochondrial cybrids with the 13708A-H7 mtDNA revealed a lower mitochondrial respiration, increased reactive oxygen species production (mROS), and dysregulation of connective tissue gene expression. The mitochondrial dysfunction was exacerbated by histamine, explaining why all eight surviving children inherited the dysfunctional histidine decarboxylase allele (W327X) from the father. Thus, certain combinations of common mtDNA variants can cause mitochondrial dysfunction, mitochondrial dysfunction can affect extracellular matrix gene expression, and histamine-activated mROS production can augment the severity of mitochondrial dysfunction. Most important, we have identified a previously unreported genetic cause of mitochondrial disorder arising from the incompatibility of common, nonpathogenic mtDNA variants.

Dissecting the multifaceted contribution of the mitochondrial genome to autism spectrum disorder

Autism spectrum disorder (ASD) is a clinically heterogeneous class of neurodevelopmental conditions with a strong, albeit complex, genetic basis. The genetic architecture of ASD includes different genetic models, from monogenic transmission at one end, to polygenic risk given by thousands of common variants with small effects at the other end. The mitochondrial DNA (mtDNA) was also proposed as a genetic modifier for ASD, mostly focusing on maternal mtDNA, since the paternal mitogenome is not transmitted to offspring. We extensively studied the potential contribution of mtDNA in ASD pathogenesis and risk through deep next generation sequencing and quantitative PCR in a cohort of 98 families. While the maternally-inherited mtDNA did not seem to predispose to ASD, neither for haplogroups nor for the presence of pathogenic mutations, an unexpected influence of paternal mtDNA, apparently centered on haplogroup U, came from the Italian families extrapolated from the test cohort (n = 74) when compared to the control population. However, this result was not replicated in an independent Italian cohort of 127 families and it is likely due to the elevated paternal age at time of conception. In addition, ASD probands showed a reduced mtDNA content when compared to their unaffected siblings. Multivariable regression analyses indicated that variants with 15%-5% heteroplasmy in probands are associated to a greater severity of ASD based on ADOS-2 criteria, whereas paternal super-haplogroups H and JT were associated with milder phenotypes. In conclusion, our results suggest that the mtDNA impacts on ASD, significantly modifying the phenotypic expression in the Italian population. The unexpected finding of protection induced by paternal mitogenome in term of severity may derive from a role of mtDNA in influencing the accumulation of nuclear de novo mutations or epigenetic alterations in fathers' germinal cells, affecting the neurodevelopment in the offspring. This result remains preliminary and needs further confirmation in independent cohorts of larger size. If confirmed, it potentially opens a different perspective on how paternal non-inherited mtDNA may predispose or modulate other complex diseases.

Loss of Drosophila Clueless differentially affects the mitochondrial proteome compared to loss of Sod2 and Pink1

Mitochondria contain their own DNA, mitochondrial DNA, which encodes thirteen proteins. However, mitochondria require thousands of proteins encoded in the nucleus to carry out their many functions. Identifying the definitive mitochondrial proteome has been challenging as methods isolating mitochondrial proteins differ and different tissues and organisms may have specialized proteomes. Mitochondrial diseases arising from single gene mutations in nucleus encoded genes could affect the mitochondrial proteome, but deciphering which effects are due to loss of specific pathways or to accumulated general mitochondrial damage is difficult. To identify specific versus general effects, we have taken advantage of mutations in three Drosophila genes, clueless, Sod2, and Pink1, which are required for mitochondrial function through different pathways. We measured changes in each mutant’s mitochondrial proteome using quantitative tandem mass tag mass spectrometry. Our analysis identified protein classes that are unique to each mutant and those shared between them, suggesting that some changes in the mitochondrial proteome are due to general mitochondrial damage whereas others are gene specific. For example, clueless mutants had the greatest number of less and more abundant mitochondrial proteins whereas loss of all three genes increased stress and metabolism proteins. This study is the first to directly compare in vivo steady state levels of mitochondrial proteins by examining loss of three pathways critical for mitochondrial function. These data could be useful to understand disease etiology, and how mutations in genes critical for mitochondrial function cause specific mitochondrial proteomic changes as opposed to changes due to generalized mitochondrial damage.

Mitochondrial DNA population variation is not associated with Alzheimer's in the Japanese population: A consistent finding across global populations

Several mitochondrial DNA (mtDNA) haplogroup association studies have suggested that common mtDNA variants are associated with multifactorial diseases, including Alzheimer’s disease (AD). However, such studies have also produced conflicting results. A new mtDNA association model, the ‘variant load model’ (VLM), has been applied to multiple disease phenotypes. Application of the VLM in a 2017 study failed to find different variant loads in AD patients compared to controls, in two cohorts of European origin. The study also suggested a lower variant load in healthy elderly individuals, but could offer no replicate cohort to support this observation. Here, the VLM is applied to Japanese mtDNA sequences; in doing so, we explored the role of mtDNA variation in AD and ageing in a different global population. Consistent with the previous findings using the VLM in two populations of European origin, we found no evidence for an association between rarer, non-haplogroup associated variation and the development of AD. However, the result in the context of ageing that suggested those with fewer mildly deleterious mutations might undergo healthier ageing, was not replicated. In contrast to our previous study, our present results suggest that those living to advanced old age may harbour more mildly deleterious mtDNA variations. Importantly our analysis showed this finding is not primarily being driven by many rare population variants dispersed across the mtDNA, but by a few more frequent variants with high MutPred scores. It is suggested the variants in question do not exert a mildly deleterious effect in their most frequent haplogroup context.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

Accurate analysis of mitochondrial DNA is important for mitochondrial disease clinical research and diagnostics. Here, authors present a method using Cas9 cleavage, nanopore sequencing and a custom pipeline to identify pathogenic variants, deletions and accurately quantify heteroplasmy to below 1%.

Whole exome and genome sequencing in mendelian disorders: a diagnostic and health economic analysis

Whole genome sequencing (WGS) improves Mendelian disorder diagnosis over whole exome sequencing (WES); however, additional diagnostic yields and costs remain undefined. We investigated differences between diagnostic and cost outcomes of WGS and WES in a cohort with suspected Mendelian disorders. WGS was performed in 38 WES-negative families derived from a 64 family Mendelian cohort that previously underwent WES. For new WGS diagnoses, contemporary WES reanalysis determined whether variants were diagnosable by original WES or unique to WGS. Diagnostic rates were estimated for WES and WGS to simulate outcomes if both had been applied to the 64 families. Diagnostic costs were calculated for various genomic testing scenarios. WGS diagnosed 34% (13/38) of WES-negative families. However, contemporary WES reanalysis on average 2 years later would have diagnosed 18% (7/38 families) resulting in a WGS-specific diagnostic yield of 19% (6/31 remaining families). In WES-negative families, the incremental cost per additional diagnosis using WGS following WES reanalysis was AU$36,710 (£19,407;US$23,727) and WGS alone was AU$41,916 (£22,159;US$27,093) compared to WES-reanalysis. When we simulated the use of WGS alone as an initial genomic test, the incremental cost for each additional diagnosis was AU$29,708 (£15,705;US$19,201) whereas contemporary WES followed by WGS was AU$36,710 (£19,407;US$23,727) compared to contemporary WES. Our findings confirm that WGS is the optimal genomic test choice for maximal diagnosis in Mendelian disorders. However, accepting a small reduction in diagnostic yield, WES with subsequent reanalysis confers the lowest costs. Whether WES or WGS is utilised will depend on clinical scenario and local resourcing and availability.