Mitochondrion

For many years, the genetic manipulation of mitochondrial DNA was largely hampered by inefficient delivery of nucleic acids to mitochondria. However, the development of mitoCBEs, such as mitochondrial cytosine base editors (DdCBEs), which catalyse C*G-to-T*A conversions, and more recently, mitoABEs, such as transcription-activator-like effector (TALE)-linked deaminases (TALEDs) enabling A*T-to-G*C conversion, has transformed this field. Generally, mitochondrial base editors exhibit high on-target efficiency and are straightforward to design and use. Nonetheless, unintended off-target effects cannot be overlooked and should be assessed consistently with each experiment, which can be challenging without specialised bioinformatic expertise. Here, we introduce Mitochondrial Base Editor Analysis Package (MitoBEAP), which, to our knowledge, is the first R package specifically designed to analyse next-generation sequencing data from base-edited mtDNA samples. The package facilitates the analysis of potential off-target effects, offers multiple visualisation options, and allows customisation of graphics and thresholds for calculations. As a proof of concept, this study demonstrates how MitoBEAP can be utilised to measure the efficiency of DdCBE treatment targeting human 12S rRNA, as well as to identify potentially harmful off-target conversions across the mtDNA.

The existence and functional relevance of mitochondrial DNA methylation remain controversial. Here, we systematically profiled cytosine methylation and hydroxymethylation across human brain and blood tissues spanning healthy and malignant states using orthogonal sequencing approaches that avoid chemical conversion during library preparation. While nuclear DNA exhibited canonical methylation patterns, mitochondrial DNA consistently showed negligible signal, indistinguishable from background technical noise. By mapping cytosine-guanine sites between mitochondrial DNA and nuclear-embedded mitochondrial sequences, we demonstrate the potential of these nuclear counterparts to confound not only cytosine methylation but also hydroxymethylation measurements, corroborating and extending prior findings implicating nuclear contamination as a potential source of apparent mitochondrial epigenetic signals. Additional technical factors that inflate apparent mtDNA methylation signals were identified, including sequence context biases, flow cell chemistries, and coverage-dependent discrepancies between the heavy and light strands. Collectively, these results provide convergent evidence against the presence of biologically meaningful cytosine methylation or hydroxymethylation in mitochondrial DNA. These findings caution against interpreting apparent mtDNA methylation signals in human adult tissues as meaningful without rigorous orthogonal validation and comprehensive consideration of technical and analytical confounding factors.

Mitochondrial DNA (mtDNA) maintenance disorders arise from defects in mtDNA replication or repair, frequently resulting in extensive deletions or depletion of mtDNA. Mitochondrial genome maintenance exonuclease 1 (MGME1) is a nuclear-encoded nuclease essential for mtDNA replication and genome stability, and biallelic pathogenic variants in MGME1 cause mitochondrial DNA depletion syndrome 11. Here, we report a novel homozygous MGME1 missense variant c.820G>A (p. Ala274Thr) in five affected individuals from unrelated South Indian families presenting with proximal myopathy, chronic progressive external ophthalmoplegia, and cardiac and renal involvement. Patient-derived cells exhibited a significant reduction in mtDNA copy number, consistent with impaired mtDNA maintenance. Mechanistic analyses combining imaging-based and biochemical approaches demonstrated that the MGME1 variant disrupts both mtDNA replication and repair. Functional characterization further revealed defective oxidative phosphorylation and reduced mitochondrial membrane potential, confirming mitochondrial dysfunction. Collectively, our findings establish the pathogenicity of this novel MGME1 variant and expand the clinical and molecular spectrum of MGME1-associated mitochondrial disease, linking impaired mtDNA replication to multisystemic mitochondrial dysfunction. Graphical abstract summarySchematic overview of mitochondrial DNA (mtDNA) replication in wild-type and MGME1 mutant conditions. In wild-type cells, the coordinated activity of the mitochondrial replication machinery (TWNK, POLG, LIG3, TFAM, and MGME1) maintains mtDNA integrity and supports normal oxidative phosphorylation and ATP production. In contrast, the MGME1 mutation disrupts mtDNA replication and processing, leading to replication defects, mtDNA depletion, and mitochondrial dysfunction. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=116 SRC="FIGDIR/small/725852v1_ufig1.gif" ALT="Figure 1"> View larger version (42K): org.highwire.dtl.DTLVardef@50187forg.highwire.dtl.DTLVardef@78799borg.highwire.dtl.DTLVardef@177e7c6org.highwire.dtl.DTLVardef@1118c68_HPS_FORMAT_FIGEXP M_FIG C_FIG

Mitochondrial Encephalopathy, Lactic Acidosis, and Stroke-like episodes (MELAS) syndrome is primarily caused by the heteroplasmic m.3243A>G/MT-TL1 pathogenic variant. Patients exhibit elevated circulating cell-free mtDNA (cf-mtDNA) in plasma, which acts as a damage-associated molecular pattern. Using patient-derived fibroblasts and neuronal progenitors, as well as transmitochondrial cytoplasmic hybrids (cybrids), we show that mutant cells release higher levels of cf-mtDNA than wild-type controls, demonstrating that the m.3243A>G pathogenic variant drives mtDNA release. Mechanistically, increased mitochondrial oxidative stress promotes mtDNA oxidation and fragmentation, leading to Ca2+ overload and subsequent mtDNA extrusion. This, in turn, triggers inflammasome activation and pyroptosis, resulting in the secretion of pro-inflammatory cytokines and the activation of innate immune pathways. Pharmacological inhibition of the Mitochondrial Calcium Uniporter (MCU) or Voltage-Dependent Anion Channel (VDAC) reduced mtDNA release, confirming their involvement. Overall, our findings reveal a previously unrecognized mechanism in MELAS linking mitochondrial dysfunction to innate immune activation, with potential implications for therapeutic intervention. TeaserMELAS mutation drives mtDNA release, triggering inflammation via oxidative stress, calcium imbalance, and inflammasome activation

Background: Biallelic variants in GFM2, encoding mitochondrial elongation factor G2 (mtEFG2), a GTPase involved in the termination stage of mitochondrial translation, cause autosomal recessive combined oxidative phosphorylation deficiency. Noncoding structural variants may be missed by exome sequencing but can disrupt splicing and provide opportunities for variant-specific therapeutic rescue. We investigated the molecular mechanism underlying suspected Leigh syndrome in an infant with mitochondrial disease and evaluated whether splice-switching oligonucleotide (SSO) treatment could correct the pathogenic splicing defect. Methods: The proband underwent exome sequencing followed by short-read and long-read whole genome sequencing. RNA sequencing, reverse-transcription PCR, quantitative PCR, and cycloheximide treatment were used to characterize the effect of the identified intronic duplication on GFM2 splicing and transcript stability. Patient-derived fibroblasts were treated with SSOs targeting the aberrant splice junction. Rescue was assessed by RNA studies, western blotting, and spectrophotometric measurement of cytochrome c oxidase (COX). Results: Whole genome sequencing identified a paternally-inherited GFM2 missense variant, NM_032380.5:c.2195C>T p.(Pro732Leu), in trans to a maternally-inherited 221-nucleotide intronic duplication, NM_032380.5:c.2029-741_2029-521dup. RNA studies revealed a 87-nucleotide pseudoexon, generated by activation of a cryptic acceptor splice site within the duplicated sequence. The resulting transcript harbored a premature termination codon (PTC) and underwent nonsense-mediated decay, as confirmed by cycloheximide rescue. Together with reduced mtEFG2 protein levels on western blot, the findings supported a loss-of-function mechanism. Enzymatic analysis of affected fibroblasts showed reduced activity of the mtDNA-dependent complex IV subunit COX, with preservation of the nuclear-encoded complex II enzyme succinate dehydrogenase and the control enzyme citrate synthase, consistent with impaired mitochondrial translation. A SSO targeting the aberrant intron-pseudoexon junction nearly abolished pseudoexon inclusion, restored correctly spliced GFM2 transcript from the duplication-containing allele, increased mtEFG2 protein levels, and significantly improved COX activity. Conclusions: This study identifies a pathogenic intronic GFM2 duplication that causes mitochondrial disease through pseudoexon activation and nonsense-mediated decay. The findings demonstrate the value of integrated genome and transcriptome analysis for exome-negative mitochondrial disease and provide in-vitro proof of concept that SSOs can restore transcript processing, protein expression, and mitochondrial respiratory-chain function in patient-derived cells.

Mitochondrial DNA (mtDNA) heteroplasmy, which is the coexistence of wild-type and mutant mtDNA variants within the same cell, plays a critical role in modulating cellular phenotype as well as disease severity and penetrance. Bulk RNA sequencing is not able to detect cell-to-cell variability in heteroplasmy, limiting our understanding of mitochondrial pathological mechanisms. In this study, we leverage single-cell RNA sequencing (scRNA-seq) combined with a robust bioinformatics pipeline to characterize mtDNA heteroplasmy. We employed four fibroblast lines from patients harboring heteroplasmic mtDNA pathogenic variants in genes encoding respiratory complex I subunits. While RNA heteroplasmy corresponds to DNA-based measurements at the bulk-level, single-cell analysis uncovers a diverged distribution: most cells have near-homoplasmic (wild-type or mutant) mtDNA, with few cells showing intermediate levels. Furthermore, we find that high mutation levels correlate with transcriptional profile changes, though these responses are highly sample-specific, suggesting that nuclear background and cellular context critically influence mitochondrial dysfunction and compensatory mechanisms. Our findings highlight the power of single-cell technologies to better understand the complex link between mtDNA genetic diversity and mitochondrial phenotypic variability, and to study crucial aspects in mitochondrial biology and pathology, such as clonal dynamics, at single-cell resolution.

Mitochondrial diseases (MDs) are genetically and phenotypically diverse and can be difficult to diagnose. Prevalence estimates derive largely from diagnosed cases and may underestimate population MD risk. Population-based studies are limited in scope and number but indicate MD variants are common. As genomic sequencing advances have made comprehensive population-based evaluation feasible, we sought to evaluate nuclear MD variation in a population cohort to understand variant prevalence and differences in MD risk estimates We identified disease-associated nuclear gene variants in 270 nuclear MD genes across 2,845 healthy older individuals in the Medical Genome Reference Bank. From Pathogenic or Likely Pathogenic Variants (PLPVs) we estimated autosomal recessive (AR) and autosomal dominant (AD) MD risk for individual genes and all nuclear variant-associated MDs. We identified 554 PLPV alleles representing 357 unique variants in 145 genes. Combined AR MD risk was estimated at 25.8 per 100,000 (95% CI 18.7 to 32.9), or 1 in 3,880 individuals. SPG7 (12.65 per 100,000; 95% CI 7.52-20.6) and POLG (4.23 per 100,000; 95% CI 2.10-8.01) contributed the greatest single gene AR MD risks and OPA1 variants posed the greatest AD MD risk. We observed a high rate of MD-associated nuclear gene variation in this healthy older cohort. The estimated lifetime AR MD risk was higher than commonly quoted prevalence estimates for all MDs, and the presence of common AD variants suggests variant penetrance may be lower than previously understood. These data help contextualise population MD risk and may inform clinical counselling and care.

Genome sequencing of the heterogeneous primary mitochondrial disorders (PMD) frequently reveals variants of uncertain significance that require functional tests for diagnosis, and does not identify variants in all patients. We analyzed mitochondrial enzyme assays, blue native polyacrylamide gel electrophoresis (BN-PAGE) with in-gel activity staining, complex I assembly blot, and select protein abundances in fibroblasts of a case series of 204 PMD patients divided into functional classes, in comparison to 51 controls and 53 differential diagnostic conditions. Overall, sensitivity and specificity for respiratory chain enzyme assays were 46% and 93% respectively, for BN-PAGE 40% and 98%, for complex I assembly assay 49% and 99%. The overall sensitivity of all tests was 76%, specificity 93%, with positive predictive value 96% and negative predictive value 67%. Categories with high sensitivity were isolated complex deficiencies, nuclear DNA-encoded mitochondrial protein synthesis defects, co-factor defects, and mitochondrial amino-acyl-tRNA synthetase conditions when aided by protein abundance. Mitochondrial DNA mutations and maintenance disorders showed poor sensitivities. Secondary dysfunctions were rare. A complete battery of functional tests showed strong diagnostic clinical utility in fibroblasts.

-Lipoic acid (LA) is widely included in "mitochondrial cocktails" recommended to patients with primary mitochondrial disorders, yet its mechanism of action remains unclear. Here, we define the intracellular availability and functional utilization of LA in mammalian cells. We show that LA exists in two functionally distinct cellular pools: a low-abundance free pool and a protein-bound pool generated through mitochondrial fatty acid synthesis (mtFAS). Disruption of the mtFAS pathway abolishes protein lipoylation and impairs oxidative phosphorylation without altering free LA levels. Conversely, supplementation with exogenous LA markedly increases free intracellular LA without restoring protein lipoylation, mitochondrial respiration, or cell proliferation. Instead, the cellular effects of LA supplementation resemble those of the antioxidant N-acetylcysteine. These findings clarify the mechanism of action of a widely used mitochondrial supplement and identify a fundamental disconnect between cellular LA abundance and mitochondrial utilization, challenging the rationale for using LA supplementation to restore mitochondrial function.

Mitochondrial DNA (mtDNA) copy number is widely used as a biomarker for mitochondrial function and disease risk, yet its relationship to mtDNA gene expression - one important functional output - remains poorly understood. Prior studies examining this relationship have largely relied on heterogeneous tissue samples, where confounding by cell-type composition obscures the underlying biology. We rigorously tested this relationship in lymphoblastoid cell lines (LCLs), where we find no correlation between mtDNA copy number and gene expression across 731 individuals, and minimal association across 49 GTEx tissues except whole blood. Using population genetic modeling of heteroplasmy drift between DNA and RNA, we estimate that effectively {approx}50 out of 813 mtDNA templates are transcriptionally active in LCLs, indicating low mtDNA accessibility. This confirms, using an independent method and a different cell type, previous observations in HeLa cells, where mtDNA is largely compacted into nucleoids. Together, our results demonstrate that mtDNA copy number and expression are largely decoupled, and that above a certain rate-limiting threshold, mtDNA accessibility -- rather than absolute copy number -- is the more relevant quantity for explaining inter-individual differences in gene expression. This challenges the interpretation of mtDNA copy number as a proxy for mitochondrial transcriptional output and highlights the need for a more mechanistic understanding of mtDNA copy number associations with disease-relevant traits. Cis- and trans-eQTL mapping further reveals that genetic regulation of mtDNA gene expression operates primarily through post-transcriptional mechanisms rather than transcription initiation, yet despite the high inter-individual variance in mtDNA gene expression, genetic variation underlying mtDNA regulation appears to be under strong selective constraint.

We show that chronic impairment of mitochondrial respiration is associated with marked accumulation of cytochrome c (Cytc) protein. Using SCO2-deficient HCT116 cells lacking functional cytochrome c oxidase and wild-type cells exposed to sustained hypoxia, we found that substantial mitochondrial Cytc accumulation parallels reduced electron flux through Cytc. SCO2-deficient cells exhibited equally elevated Cytc levels under normoxia (19% O2) and hypoxia (0.1-3% O2). Wild-type cells under sustained hypoxia accumulated Cytc, reaching levels comparable to those in SCO2-deficient cells. This effect was reversible upon reoxygenation. Increased Cytc protein levels were also observed in other cell models, including primary cortical neurons cultured under chronic hypoxia and in cerebral cortex tissue from hypoxia-exposed mice. Cytc accumulation occurred independently of CYCS transcription, mRNA translation, HIF activation, ROS production and changes in mitochondrial network. Pharmacological inhibition of complex III was likewise accompanied by increased Cytc levels, whereas mitochondrial uncoupling had no effect, suggesting that impaired electron transfer rather than membrane depolarisation per se underlies this association. Raman spectroscopy revealed enrichment of reduced Cytc and an increased Cytc-to-cytochrome b ratio in respiration-deficient cells. Further supporting a stabilisation-based mechanism, the fraction of membrane-unbound ferro-Cytc was decreased in SCO2-deficient cells, consistent with moderate cardiolipin enrichment, which is known to enhance retention of Cytc at the inner mitochondrial membrane. Despite elevated mitochondrial Cytc content, SCO2-deficient cells were less susceptible to apoptosis induced by intermittent hypoxia or dichloroacetate. Together, these findings indicate that reduced electron flux through complex IV is associated with Cytc accumulation through increased protein stability and membrane retention without enhancing apoptotic sensitivity.

Disorders of mitochondrial oxidative phosphorylation affecting multiple respiratory chain complexes are among the most common causes of mitochondrial disease in humans. However, impaired energy metabolism alone does not fully account for tissue-specific vulnerability and disease progression, suggesting that additional molecular mechanisms contribute to disease pathology. We previously identified PTCD1 variants in a child with infantile cardiomyopathy associated with a combined respiratory chain deficiency. Here, we establish the pathogenicity of three PTCD1 (NM_015545.4) variants in which p.(Arg113Trp) and p.(Gly184Arg) segregate in cis whereas p.(Arg130*) is present in trans, demonstrating that disrupted mitochondrial proteostasis contributes to tissue damage in PTCD1 deficiency. PTCD1 patient cardiac tissue characterisation revealed impaired mitoribosome biogenesis, alongside increased aggregation of selective mitochondrial matrix proteins. Cell models expressing individual and combined PTCD1 missense variants, coupled with proteomics, recapitulated the protein aggregation, with the cis p.(Arg113Trp);p.(Gly184Arg) combination showing the most severe effect. Protein aggregation was accompanied by altered OPA1 processing and mitochondrial network remodelling. Our findings establish accumulating proteotoxic stress arising from impaired mitoribosome assembly as a pathogenic mechanism in post-mitotic tissues, driving PTCD1 cardiomyopathy.

Codon usage bias is a fundamental genomic characteristic that prefers non-random preferential use of synonymous codons. It is a major determinant of translational efficiency, gene regulation, and molecular evolution. However, the evolutionary bias and functional relevance of codon usage bias across the plant lineage is poorly defined and yet to understand what are the major factors responsible for relative synonymous codon usage (RSCU) in genomes and how codon usage bias influences the gene regulation, molecular evolution genomes. A genome-wide codon usage bias study of coding DNA sequences of 262 plant genome was conducted. It encompassed more than 4.6 billion codons from > 11 million coding sequences. Relative synonymous codon usage, codon adaptation index, codon-anticodon mapping, effective number of codon (ENC)-GC3, GC1,2-GC3, parity rule 2 (PR2-bias), molecular economy, and machine learning approaches were used for the study. It was found that codon usage bias was strongly non-random and exhibited a clear phylogenetic structuring. The higher plants favoured A/T-ending, whereas early-diverging lineages were enriched in G/C-ending codons. Analysis of RSCU, codon adaptation index, and codon-anticodon pairing indicated that translational selection is mediated by tRNA availability, contributing sustainability to these molecular patterns. Machine-learning approaches identified a small subset of codons having outsized influence on genome-wide codon usage landscapes. Further studies revealed the presence of robust inverse relationships between the effective number of codons and GC content at synonymous third positions. Neutrality analysis revealed approximately 61% of variation was driven by mutational pressure, tempered by selective constraints. Phylogenetic reconstruction showed a progressive relaxation of codon bias from algae to angiosperms while maintaining a conserved molecular economy cost of ~ 30 ATP per codon across the lineages. The study revealed codon usage bias is lineage-specific evolutionary conserved trait governed by mutation, selection, and translational optimization.

Multiple acyl-CoA dehydrogenase deficiency (MADD) is a mitochondrial lipid storage myopathy characterized by impaired fatty acid {beta}-oxidation, mitochondrial dysfunction, and progressive neuromuscular and cardiac disease. MADD is most commonly caused by pathogenic variants in electron transfer flavoprotein dehydrogenase (ETFDH), which encodes electron transfer flavoprotein-ubiquinone oxidoreductase (Etf-QO), a critical redox enzyme that transfers electrons from acyl-CoA dehydrogenases to the mitochondrial electron transport chain. Defective Etf-QO activity disrupts electron flow, promotes reactive oxygen species (ROS) production, and impairs cellular energy metabolism, linking abnormal lipid oxidation to oxidative stress-mediated tissue damage. To investigate the role of redox imbalance in MADD pathogenesis, we generated CRISPR/Cas9 knock-in Drosophila melanogaster models carrying patient-relevant Etf-QO missense mutations (L127R, S296C, and L399F; corresponding to human L138R, S307C, and L409F) within conserved FAD- and ubiquinone-binding domains. Mutant flies developed progressive locomotor impairment, reduced muscle performance, and marked lipid droplet accumulation in skeletal muscle, cardiac tissue, and fat bodies, indicating systemic defects in mitochondrial lipid utilization. Cardiac analyses demonstrated reduced fractional shortening, prolonged heart period, and increased arrhythmia index, consistent with metabolic cardiomyopathy associated with mitochondrial oxidative stress. In vivo respirometry revealed significantly decreased oxygen consumption, reflecting impaired oxidative phosphorylation. At the molecular level, mutant flies exhibited elevated ROS levels and ATP depletion, accompanied by increased expression of AMPK, PGC-1, and Tfam, suggesting activation of energy stress signaling and compensatory mitochondrial biogenesis. Importantly, endurance exercise significantly improved locomotor and cardiac function while reducing lipid accumulation and oxidative stress. Together, these findings establish a redox-centered in vivo model of MADD and identify oxidative stress as a major driver of disease pathology and a potential therapeutic target.

Many flowering plant species have adopted an aquatic lifestyle, contrasting with their terrestrial ancestors. Adapting to an aquatic environment required numerous evolutionary changes, including gene expansion and contraction. One of the most striking contractions has been observed in the genomes of seagrasses, where the ACO and ACS genes, involved in ethylene biosynthesis, are very few in number or even completely absent. To confirm this adaptation, we identified traces of gene loss in the genomes of four seagrass species, in the form of pseudogenes. Surprisingly, no gene loss was found in the species that had completely lost the function of ethylene synthesis, likely indicating an ancient loss of these genes. Conversely, several pseudogenes were found in the species where the ACO and ACS genes are contracting, indicating a recent and potentially ongoing process. We used the same approach on Utricularia gibba, a submerged freshwater plant, and also found a reduced number of ACO and ACS genes. In contrast, two terrestrial species closely related to seagrasses and U. gibba found a higher number of ACO and ACS genes, with no definitive evidence of gene loss. These results confirm that the loss of ethylene biosynthesis function in seagrasses is indeed linked to gene loss and suggests that it is an adaptation to a submerged rather than a marine lifestyle.

Malonate is often described as an endogenous inhibitor of complex II of the electron transport chain. However, the cellular source of malonate is unclear, and current knowledge concerning its metabolism is limited to the action of a single enzyme, Acyl-CoA Synthetase Family Member 3 (ACSF3), which converts malonate to malonyl-CoA in the mitochondrial matrix. One potential route of malonate metabolism downstream of ACSF3 is its consumption by the mitochondrial fatty acid synthesis (mtFAS) pathway. However, studies examining the link between ACSF3 and mtFAS have yielded conflicting results. We developed a novel mass spectrometry approach to perform stable isotope tracing into products of mtFAS, and found that while malonate is in fact a carbon source for mtFAS, ACSF3 is not required for malonate incorporation into mtFAS products. Using this method to trace other nutrients into mtFAS, we also found evidence of acetyl-CoA carboxylase 1 (ACC1)-dependent malonate synthesis from glucose. We further show that ACC1 is required for optimal mtFAS activity, with downstream effects on oxidative phosphorylation. Together these findings establish the malonate as a regulated endogenous intermediate that supports mtFAS activity and mitochondrial oxidative function.

Expression of the chloroplast AAA+ chaperone CLPD gene increases during senescence and drought, but its functional role in chloroplast proteostasis is poorly understood. This study provides a comprehensive analysis of Arabidopsis CLPD protein accumulation across development from early seedlings to senescence, and compares results to its homologs CLPC1,2, as well as CLPB3 and cpHSP90. The developmental consequences of complete loss of CLPD expression (clpd-1), as well as overexpression of functional CLPD or CLPD impaired in ATP hydrolysis (CLPD-TRAP), were determined in Arabidopsis. clpd-1 has accelerated seedling development while functional CLPD overexpression lines, but not CLPD-TRAP, have delayed development. To determine if CLPD is a bona fide CLP chaperone associating with the CLPPRT protease and to identify in vivo candidate substrates, we employed the CLPD-TRAP line during the vegetative and flowering (senescent) growth stages. Affinity purification of CLPD-TRAP followed by mass spectrometry showed high enrichment of the CLP protease complex, thus providing direct support for the role of CLPD in substrate delivery to the CLP protease. CLPC1,2 were also highly enriched in the CLPD-TRAP interactome, suggesting hetero-oligomerization and cooperation between the three chaperones is likely. Nine chloroplast candidate substrates were identified in the CLPD-interactomes, including: FHY2 involved in riboflavin synthesis, THI1 and THIC involved in thiamin metabolism, and four proteins of unknown function. Several of these have been previously identified as potential CLPC1 substrates; however, others appear to be specific to CLPD. CLPD acts in substrate selection within a heteromeric CLPC-CLPD hexamer, likely to make unique contributions through its divergent N-terminus.

Significance Missense mutations can lead to pathological effects in human cells. Predictive methods that account for structural context, such as AlphaMissense, can provide pathogenicity scores. The accumulation of pathogenicity hotspots can reveal important structural features within individual proteins of protein families, such as GLUT transporters. Mapping pathogenicity scores onto the structure can thus provide a mechanistic explanation of the protein function necessary for its role in the cell. Abstract Non-synonymous amino acid substitutions (missense mutations) are common in the general population; some are causative of serious disease. Depending on their structural context, they can disrupt protein function, folding, or dynamics. Computational predictive methods developed in recent years, such as AlphaMissense, provide new insights into how missense mutations affect protein structure by predicting and mapping their pathogenicity across each amino acid in the human proteome. In this study, we identify recurring patterns of pathogenicity prediction across the GLUT family membrane transporters encoded by genes slc2a1-14. Within the GLUT transporter family, we observe higher pathogenicity profiles in the transmembrane domains, particularly in pore-lining and binding-site residues. Predicted missense pathogenicity is elevated throughout residues assigned to the central cavity, suggesting sensitivity of the transport pathway. Another finding shows higher pathogenicity in specific transmembrane helices of the protein, with the same pattern across all proteins. On the other hand, we observed lower pathogenicity values in some representatives of the GLUT family. These findings show that the pathogenicity of glucose transport within the GLUT family may be shaped by functional redundancy and physiological essentiality across GLUT groups.

Differentiation of skeletal muscle is associated with increased mitochondrial biogenesis and reliance of oxidative phosphorylation (OXPHOS). The terminal enzyme complex in the electron transport chain, cytochrome c oxidase (COX), requires copper for its assembly and activity, and copper delivery to mitochondria is essential for OXPHOS. However, when mitochondrial copper becomes essential during skeletal myoblast differentiation is not known. Here, we show that genetic deficiency of the mitochondrial copper and phosphate carrier SLC25A3 induced prior to myoblast differentiation leads to the formation of smaller myotubes, but SLC25A3 deficiency induced in mature myotubes leads to cell death and detachment. Both phenotypes are recapitulated upon genetic knockdown of COX17, a critical assembly protein for both COX copper cofactors, or by chemical inhibition of COX. Importantly, myotube death caused by SLC25A3 deficiency is rescued by copper supplementation or expression of an SLC25A3 variant that transports copper but not phosphate. Taken together these data support a model wherein copper transport by SLC25A3 and copper delivery to COX is critical for survival in mature myotubes.

Ancient environmental DNA (aeDNA) from permafrost, lake, cave, and marine sediments provides a rich source of genetic data that captures broad perspectives of past biodiversity. Accurate dating is crucial for discovering ecologically relevant patterns from aeDNA, and molecular clock dating would allow for sample ages to be estimated from the recovered genetic material itself instead of the geological components. However, the fragmented and damaged nature of short-read ancient DNA (aDNA) from multiple taxonomic sources poses significant challenges and has limited this dating approach for aeDNA. Here we developed ratePlacer, a phylogeny-based method for analyzing aeDNA that can combine information from many short reads in a sample while accounting for DNA damage to provide maximum likelihood estimates of sample ages. Simulations demonstrate that ratePlacer accurately dates samples even under the fragmented, damaged conditions characteristic of aeDNA and outperforms Bayesian tip-dating approaches for taxonomically mixed samples commonly found in aeDNA. Yet age estimates from re-dating Kap Kobenhavn varied across taxa, highlighting the difficulty of molecular clock dating in aeDNA. This dating also revealed elevated G[->]T and C[->]A mismatches consistent with oxidative damage. These patterns reveal aDNA damage beyond deamination and that remains understudied, suggesting that aeDNA should be carefully evaluated in genomic and evolutionary analyses. The new dating method, ratePlacer, extends molecular clock dating of aDNA from single-specimen to pooled environmental DNA data, where traditional methods struggle.