ChREBP is activated by reductive stress and mediates GCKR-associated metabolic traits.

Common genetic variants in glucokinase regulator (GCKR), which encodes GKRP, a regulator of hepatic glucokinase (GCK), influence multiple metabolic traits in genome-wide association studies (GWASs), making GCKR one of the most pleiotropic GWAS loci in the genome. It is unclear why. Prior work has demonstrated that GCKR influences the hepatic cytosolic NADH/NAD+ ratio, also referred to as reductive stress. Here, we demonstrate that reductive stress is sufficient to activate the transcription factor ChREBP and necessary for its activation by the GKRP-GCK interaction, glucose, and ethanol. We show that hepatic reductive stress induces GCKR GWAS traits such as increased hepatic fat, circulating FGF21, and circulating acylglycerol species, which are also influenced by ChREBP. We define the transcriptional signature of hepatic reductive stress and show its upregulation in fatty liver disease and downregulation after bariatric surgery in humans. These findings highlight how a GCKR-reductive stress-ChREBP axis influences multiple human metabolic traits.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans.

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

Effectors Enabling Adaptation to Mitochondrial Complex I Loss in Hürthle Cell Carcinoma.

Oncocytic (Hürthle cell) carcinoma of the thyroid (HCC) is genetically characterized by complex I mitochondrial DNA mutations and widespread chromosomal losses. Here, we utilize RNA sequencing and metabolomics to identify candidate molecular effectors activated by these genetic drivers. We find glutathione biosynthesis, amino acid metabolism, mitochondrial unfolded protein response, and lipid peroxide scavenging to be increased in HCC. A CRISPR-Cas9 knockout screen in a new HCC model reveals which pathways are key for fitness, and highlights loss of GPX4, a defense against lipid peroxides and ferroptosis, as a strong liability. Rescuing complex I redox activity with the yeast NADH dehydrogenase (NDI1) in HCC cells diminishes ferroptosis sensitivity, while inhibiting complex I in normal thyroid cells augments ferroptosis induction. Our work demonstrates unmitigated lipid peroxide stress to be an HCC vulnerability that is mechanistically coupled to the genetic loss of mitochondrial complex I activity. SIGNIFICANCE: HCC harbors abundant mitochondria, mitochondrial DNA mutations, and chromosomal losses. Using a CRISPR-Cas9 screen inspired by transcriptomic and metabolomic profiling, we identify molecular effectors essential for cell fitness. We uncover lipid peroxide stress as a vulnerability coupled to mitochondrial complex I loss in HCC. See related article by Frank et al., p. 1884. This article is highlighted in the In This Issue feature, p. 1749.

Multi-omics identifies large mitoribosomal subunit instability caused by pathogenic MRPL39 variants as a cause of pediatric onset mitochondrial disease.

MRPL39 encodes one of 52 proteins comprising the large subunit of the mitochondrial ribosome (mitoribosome). In conjunction with 30 proteins in the small subunit, the mitoribosome synthesizes the 13 subunits of the mitochondrial oxidative phosphorylation (OXPHOS) system encoded by mitochondrial Deoxyribonucleic acid (DNA). We used multi-omics and gene matching to identify three unrelated individuals with biallelic variants in MRPL39 presenting with multisystem diseases with severity ranging from lethal, infantile-onset (Leigh syndrome spectrum) to milder with survival into adulthood. Clinical exome sequencing of known disease genes failed to diagnose these patients; however quantitative proteomics identified a specific decrease in the abundance of large but not small mitoribosomal subunits in fibroblasts from the two patients with severe phenotype. Re-analysis of exome sequencing led to the identification of candidate single heterozygous variants in mitoribosomal genes MRPL39 (both patients) and MRPL15. Genome sequencing identified a shared deep intronic MRPL39 variant predicted to generate a cryptic exon, with transcriptomics and targeted studies providing further functional evidence for causation. The patient with the milder disease was homozygous for a missense variant identified through trio exome sequencing. Our study highlights the utility of quantitative proteomics in detecting protein signatures and in characterizing gene-disease associations in exome-unsolved patients. We describe Relative Complex Abundance analysis of proteomics data, a sensitive method that can identify defects in OXPHOS disorders to a similar or greater sensitivity to the traditional enzymology. Relative Complex Abundance has potential utility for functional validation or prioritization in many hundreds of inherited rare diseases where protein complex assembly is disrupted.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans.

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.

Combinatorial GxGxE CRISPR screen identifies SLC25A39 in mitochondrial glutathione transport linking iron homeostasis to OXPHOS.

The SLC25 carrier family consists of 53 transporters that shuttle nutrients and co-factors across mitochondrial membranes. The family is highly redundant and their transport activities coupled to metabolic state. Here, we use a pooled, dual CRISPR screening strategy that knocks out pairs of transporters in four metabolic states - glucose, galactose, OXPHOS inhibition, and absence of pyruvate - designed to unmask the inter-dependence of these genes. In total, we screen 63 genes in four metabolic states, corresponding to 2016 single and pair-wise genetic perturbations. We recover 19 gene-by-environment (GxE) interactions and 9 gene-by-gene (GxG) interactions. One GxE interaction hit illustrates that the fitness defect in the mitochondrial folate carrier (SLC25A32) KO cells is genetically buffered in galactose due to a lack of substrate in de novo purine biosynthesis. GxG analysis highlights a buffering interaction between the iron transporter SLC25A37 (A37) and the poorly characterized SLC25A39 (A39). Mitochondrial metabolite profiling, organelle transport assays, and structure-guided mutagenesis identify A39 as critical for mitochondrial glutathione (GSH) import. Functional studies reveal that A39-mediated glutathione homeostasis and A37-mediated mitochondrial iron uptake operate jointly to support mitochondrial OXPHOS. Our work underscores the value of studying family-wide genetic interactions across different metabolic environments.

Mitochondrial DNA variation across 56,434 individuals in gnomAD.

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

Loss of LUC7L2 and U1 snRNP subunits shifts energy metabolism from glycolysis to OXPHOS.

Oxidative phosphorylation (OXPHOS) and glycolysis are the two major pathways for ATP production. The reliance on each varies across tissues and cell states, and can influence susceptibility to disease. At present, the full set of molecular mechanisms governing the relative expression and balance of these two pathways is unknown. Here, we focus on genes whose loss leads to an increase in OXPHOS activity. Unexpectedly, this class of genes is enriched for components of the pre-mRNA splicing machinery, in particular for subunits of the U1 snRNP. Among them, we show that LUC7L2 represses OXPHOS and promotes glycolysis by multiple mechanisms, including (1) splicing of the glycolytic enzyme PFKM to suppress glycogen synthesis, (2) splicing of the cystine/glutamate antiporter SLC7A11 (xCT) to suppress glutamate oxidation, and (3) secondary repression of mitochondrial respiratory supercomplex formation. Our results connect LUC7L2 expression and, more generally, the U1 snRNP to cellular energy metabolism.

Fatal perinatal mitochondrial cardiac failure caused by recurrent de novo duplications in the ATAD3 locus.

BACKGROUND: In about half of all patients with a suspected monogenic disease, genomic investigations fail to identify the diagnosis. A contributing factor is the difficulty with repetitive regions of the genome, such as those generated by segmental duplications. The ATAD3 locus is one such region, in which recessive deletions and dominant duplications have recently been reported to cause lethal perinatal mitochondrial diseases characterized by pontocerebellar hypoplasia or cardiomyopathy, respectively. METHODS: Whole exome, whole genome and long-read DNA sequencing techniques combined with studies of RNA and quantitative proteomics were used to investigate 17 subjects from 16 unrelated families with suspected mitochondrial disease. FINDINGS: We report six different de novo duplications in the ATAD3 gene locus causing a distinctive presentation including lethal perinatal cardiomyopathy, persistent hyperlactacidemia, and frequently corneal clouding or cataracts and encephalopathy. The recurrent 68 Kb ATAD3 duplications are identifiable from genome and exome sequencing but usually missed by microarrays. The ATAD3 duplications result in the formation of identical chimeric ATAD3A/ATAD3C proteins, altered ATAD3 complexes and a striking reduction in mitochondrial oxidative phosphorylation complex I and its activity in heart tissue. CONCLUSIONS: ATAD3 duplications appear to act in a dominant-negative manner and the de novo inheritance infers a low recurrence risk for families, unlike most pediatric mitochondrial diseases. More than 350 genes underlie mitochondrial diseases. In our experience the ATAD3 locus is now one of the five most common causes of nuclear-encoded pediatric mitochondrial disease but the repetitive nature of the locus means ATAD3 diagnoses may be frequently missed by current genomic strategies. FUNDING: Australian NHMRC, US Department of Defense, Japanese AMED and JSPS agencies, Australian Genomics Health Alliance and Australian Mito Foundation.

MitoCarta3.0: an updated mitochondrial proteome now with sub-organelle localization and pathway annotations.

The mammalian mitochondrial proteome is under dual genomic control, with 99% of proteins encoded by the nuclear genome and 13 originating from the mitochondrial DNA (mtDNA). We previously developed MitoCarta, a catalogue of over 1000 genes encoding the mammalian mitochondrial proteome. This catalogue was compiled using a Bayesian integration of multiple sequence features and experimental datasets, notably protein mass spectrometry of mitochondria isolated from fourteen murine tissues. Here, we introduce MitoCarta3.0. Beginning with the MitoCarta2.0 inventory, we performed manual review to remove 100 genes and introduce 78 additional genes, arriving at an updated inventory of 1136 human genes. We now include manually curated annotations of sub-mitochondrial localization (matrix, inner membrane, intermembrane space, outer membrane) as well as assignment to 149 hierarchical 'MitoPathways' spanning seven broad functional categories relevant to mitochondria. MitoCarta3.0, including sub-mitochondrial localization and MitoPathway annotations, is freely available at http://www.broadinstitute.org/mitocarta and should serve as a continued community resource for mitochondrial biology and medicine.

Genetic Screen for Cell Fitness in High or Low Oxygen Highlights Mitochondrial and Lipid Metabolism.

Human cells are able to sense and adapt to variations in oxygen levels. Historically, much research in this field has focused on hypoxia-inducible factor (HIF) signaling and reactive oxygen species (ROS). Here, we perform genome-wide CRISPR growth screens at 21%, 5%, and 1% oxygen to systematically identify gene knockouts with relative fitness defects in high oxygen (213 genes) or low oxygen (109 genes), most without known connection to HIF or ROS. Knockouts of many mitochondrial pathways thought to be essential, including complex I and enzymes in Fe-S biosynthesis, grow relatively well at low oxygen and thus are buffered by hypoxia. In contrast, in certain cell types, knockout of lipid biosynthetic and peroxisomal genes causes fitness defects only in low oxygen. Our resource nominates genetic diseases whose severity may be modulated by oxygen and links hundreds of genes to oxygen homeostasis.

A patient with homozygous nonsense variants in two Leigh syndrome disease genes: Distinguishing a dual diagnosis from a hypomorphic protein-truncating variant.

Leigh syndrome is a mitochondrial disease caused by pathogenic variants in over 85 genes. Whole exome sequencing of a patient with Leigh-like syndrome identified homozygous protein-truncating variants in two genes associated with Leigh syndrome; a reported pathogenic variant in PDHX (NP_003468.2:p.(Arg446*)), and an uncharacterized variant in complex I (CI) assembly factor TIMMDC1 (NP_057673.2:p.(Arg225*)). The TIMMDC1 variant was predicted to truncate 61 amino acids at the C-terminus and functional studies demonstrated a hypomorphic impact of the variant on CI assembly. However, the mutant protein could still rescue CI assembly in TIMMDC1 knockout cells and the patient's clinical phenotype was not clearly distinct from that of other patients with the same PDHX defect. Our data suggest that the hypomorphic effect of the TIMMDC1 protein-truncating variant does not constitute a dual diagnosis in this individual. We recommend cautious assessment of variants in the C-terminus of TIMMDC1 and emphasize the need to consider the caveats detailed within the American College of Medical Genetics and Genomics (ACMG) criteria when assessing variants.

Hypoxia Rescues Frataxin Loss by Restoring Iron Sulfur Cluster Biogenesis.

Friedreich's ataxia (FRDA) is a devastating, multisystemic disorder caused by recessive mutations in the mitochondrial protein frataxin (FXN). FXN participates in the biosynthesis of Fe-S clusters and is considered to be essential for viability. Here we report that when grown in 1% ambient O2, FXN null yeast, human cells, and nematodes are fully viable. In human cells, hypoxia restores steady-state levels of Fe-S clusters and normalizes ATF4, NRF2, and IRP2 signaling events associated with FRDA. Cellular studies and in vitro reconstitution indicate that hypoxia acts through HIF-independent mechanisms that increase bioavailable iron as well as directly activate Fe-S synthesis. In a mouse model of FRDA, breathing 11% O2 attenuates the progression of ataxia, whereas breathing 55% O2 hastens it. Our work identifies oxygen as a key environmental variable in the pathogenesis associated with FXN depletion, with important mechanistic and therapeutic implications.

Widespread Chromosomal Losses and Mitochondrial DNA Alterations as Genetic Drivers in Hürthle Cell Carcinoma.

Hürthle cell carcinoma of the thyroid (HCC) is a form of thyroid cancer recalcitrant to radioiodine therapy that exhibits an accumulation of mitochondria. We performed whole-exome sequencing on a cohort of primary, recurrent, and metastatic tumors, and identified recurrent mutations in DAXX, TP53, NRAS, NF1, CDKN1A, ARHGAP35, and the TERT promoter. Parallel analysis of mtDNA revealed recurrent homoplasmic mutations in subunits of complex I of the electron transport chain. Analysis of DNA copy-number alterations uncovered widespread loss of chromosomes culminating in near-haploid chromosomal content in a large fraction of HCC, which was maintained during metastatic spread. This work uncovers a distinct molecular origin of HCC compared with other thyroid malignancies.

Early loss of mitochondrial complex I and rewiring of glutathione metabolism in renal oncocytoma.

Renal oncocytomas are benign tumors characterized by a marked accumulation of mitochondria. We report a combined exome, transcriptome, and metabolome analysis of these tumors. Joint analysis of the nuclear and mitochondrial (mtDNA) genomes reveals loss-of-function mtDNA mutations occurring at high variant allele fractions, consistent with positive selection, in genes encoding complex I as the most frequent genetic events. A subset of these tumors also exhibits chromosome 1 loss and/or cyclin D1 overexpression, suggesting they follow complex I loss. Transcriptome data revealed that many pathways previously reported to be altered in renal oncocytoma were simply differentially expressed in the tumor's cell of origin, the distal nephron, compared with other nephron segments. Using a heuristic approach to account for cell-of-origin bias we uncovered strong expression alterations in the gamma-glutamyl cycle, including glutathione synthesis (increased GCLC) and glutathione degradation. Moreover, the most striking changes in metabolite profiling were elevations in oxidized and reduced glutathione as well as γ-glutamyl-cysteine and cysteinyl-glycine, dipeptide intermediates in glutathione biosynthesis, and recycling, respectively. Biosynthesis of glutathione appears adaptive as blockade of GCLC impairs viability in cells cultured with a complex I inhibitor. Our data suggest that loss-of-function mutations in complex I are a candidate driver event in renal oncocytoma that is followed by frequent loss of chromosome 1, cyclin D1 overexpression, and adaptive up-regulation of glutathione biosynthesis.

Spatiotemporal compartmentalization of hepatic NADH and NADPH metabolism.

Compartmentalization is a fundamental design principle of eukaryotic metabolism. Here, we review the compartmentalization of NAD+/NADH and NADP+/NADPH with a focus on the liver, an organ that experiences the extremes of biochemical physiology each day. Historical studies of the liver, using classical biochemical fractionation and measurements of redox-coupled metabolites, have given rise to the prevailing view that mitochondrial NAD(H) pools tend to be oxidized and important for energy homeostasis, whereas cytosolic NADP(H) pools tend to be highly reduced for reductive biosynthesis. Despite this textbook view, many questions still remain as to the relative size of these subcellular pools and their redox ratios in different physiological states, and to what extent such redox ratios are simply indicators versus drivers of metabolism. By performing a bioinformatic survey, we find that the liver expresses 352 known or predicted enzymes composing the hepatic NAD(P)ome, i.e. the union of all predicted enzymes producing or consuming NADP(H) or NAD(H) or using them as a redox co-factor. Notably, less than half are predicted to be localized within the cytosol or mitochondria, and a very large fraction of these genes exhibit gene expression patterns that vary during the time of day or in response to fasting or feeding. A future challenge lies in applying emerging new genetic tools to measure and manipulate in vivo hepatic NADP(H) and NAD(H) with subcellular and temporal resolution. Insights from such fundamental studies will be crucial in deciphering the pathogenesis of very common diseases known to involve alterations in hepatic NAD(P)H, such as diabetes and fatty liver disease.

Defective mitochondrial rRNA methyltransferase MRM2 causes MELAS-like clinical syndrome.

Defects in nuclear-encoded proteins of the mitochondrial translation machinery cause early-onset and tissue-specific deficiency of one or more OXPHOS complexes. Here, we report a 7-year-old Italian boy with childhood-onset rapidly progressive encephalomyopathy and stroke-like episodes. Multiple OXPHOS defects and decreased mtDNA copy number (40%) were detected in muscle homogenate. Clinical features combined with low level of plasma citrulline were highly suggestive of mitochondrial encephalopathy, lactic acidosis and stroke-like episodes (MELAS) syndrome, however, the common m.3243 A > G mutation was excluded. Targeted exome sequencing of genes encoding the mitochondrial proteome identified a damaging mutation, c.567 G > A, affecting a highly conserved amino acid residue (p.Gly189Arg) of the MRM2 protein. MRM2 has never before been linked to a human disease and encodes an enzyme responsible for 2'-O-methyl modification at position U1369 in the human mitochondrial 16S rRNA. We generated a knockout yeast model for the orthologous gene that showed a defect in respiration and the reduction of the 2'-O-methyl modification at the equivalent position (U2791) in the yeast mitochondrial 21S rRNA. Complementation with the mrm2 allele carrying the equivalent yeast mutation failed to rescue the respiratory phenotype, which was instead completely rescued by expressing the wild-type allele. Our findings establish that defective MRM2 causes a MELAS-like phenotype, and suggests the genetic screening of the MRM2 gene in patients with a m.3243 A > G negative MELAS-like presentation.

Biallelic C1QBP Mutations Cause Severe Neonatal-, Childhood-, or Later-Onset Cardiomyopathy Associated with Combined Respiratory-Chain Deficiencies.

Complement component 1 Q subcomponent-binding protein (C1QBP; also known as p32) is a multi-compartmental protein whose precise function remains unknown. It is an evolutionary conserved multifunctional protein localized primarily in the mitochondrial matrix and has roles in inflammation and infection processes, mitochondrial ribosome biogenesis, and regulation of apoptosis and nuclear transcription. It has an N-terminal mitochondrial targeting peptide that is proteolytically processed after import into the mitochondrial matrix, where it forms a homotrimeric complex organized in a doughnut-shaped structure. Although C1QBP has been reported to exert pleiotropic effects on many cellular processes, we report here four individuals from unrelated families where biallelic mutations in C1QBP cause a defect in mitochondrial energy metabolism. Infants presented with cardiomyopathy accompanied by multisystemic involvement (liver, kidney, and brain), and children and adults presented with myopathy and progressive external ophthalmoplegia. Multiple mitochondrial respiratory-chain defects, associated with the accumulation of multiple deletions of mitochondrial DNA in the later-onset myopathic cases, were identified in all affected individuals. Steady-state C1QBP levels were decreased in all individuals' samples, leading to combined respiratory-chain enzyme deficiency of complexes I, III, and IV. C1qbp-/- mouse embryonic fibroblasts (MEFs) resembled the human disease phenotype by showing multiple defects in oxidative phosphorylation (OXPHOS). Complementation with wild-type, but not mutagenized, C1qbp restored OXPHOS protein levels and mitochondrial enzyme activities in C1qbp-/- MEFs. C1QBP deficiency represents an important mitochondrial disorder associated with a clinical spectrum ranging from infantile lactic acidosis to childhood (cardio)myopathy and late-onset progressive external ophthalmoplegia.