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.

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.

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.

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: (i) detecting homoplasmic and heteroplasmic variants, present respectively in all or a fraction of mtDNA molecules, (ii) circular mtDNA genome, and (iii) misalignment of nuclear sequences of mitochondrial origin (NUMTs). We observed that mtDNA copy number per cell varied across gnomAD cohorts and influenced the fraction of NUMT-derived false-positive variant calls, which can account for the majority of putative heteroplasmies. To avoid false positives, we excluded contaminated samples, cell lines, and samples prone to NUMT misalignment due to few mtDNA copies. Furthermore, we report variants with heteroplasmy greater than 10%. We applied this pipeline to 56,434 whole genome sequences in the gnomAD v3.1 database that includes individuals of European (58%), African (25%), Latino (10%), and Asian (5%) ancestry. Our gnomAD v3.1 release contains population frequencies for 10,850 unique mtDNA variants at more than half of all mtDNA bases. We report frequencies within each nuclear ancestral population and mitochondrial haplogroup. Homoplasmic variants account for most variant calls (98%) and unique variants (85%). We observed that 1/250 individuals carry a pathogenic mtDNA variant with heteroplasmy above 10%. These mtDNA population allele frequencies are freely accessible and will aid in diagnostic interpretation and research studies.

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.

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.

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.

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.

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.

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.

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.

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.

Biallelic Mutations in MRPS34 Lead to Instability of the Small Mitoribosomal Subunit and Leigh Syndrome

The synthesis of all 13 mitochondrial DNA (mtDNA)-encoded protein subunits of the human oxidative phosphorylation (OXPHOS) system is carried out by mitochondrial ribosomes (mitoribosomes). Defects in the stability of mitoribosomal proteins or mitoribosome assembly impair mitochondrial protein translation, causing combined OXPHOS enzyme deficiency and clinical disease. Here we report four autosomal-recessive pathogenic mutations in the gene encoding the small mitoribosomal subunit protein, MRPS34, in six subjects from four unrelated families with Leigh syndrome and combined OXPHOS defects. Whole-exome sequencing was used to independently identify all variants. Two splice-site mutations were identified, including homozygous c.321+1G>T in a subject of Italian ancestry and homozygous c.322-10G>A in affected sibling pairs from two unrelated families of Puerto Rican descent. In addition, compound heterozygous MRPS34 mutations were identified in a proband of French ancestry; a missense (c.37G>A [p.Glu13Lys]) and a nonsense (c.94C>T [p.Gln32^∗]) variant. We demonstrated that these mutations reduce MRPS34 protein levels and the synthesis of OXPHOS subunits encoded by mtDNA. Examination of the mitoribosome profile and quantitative proteomics showed that the mitochondrial translation defect was caused by destabilization of the small mitoribosomal subunit and impaired monosome assembly. Lentiviral-mediated expression of wild-type MRPS34 rescued the defect in mitochondrial translation observed in skin fibroblasts from affected subjects, confirming the pathogenicity of MRPS34 mutations. Our data establish that MRPS34 is required for normal function of the mitoribosome in humans and furthermore demonstrate the power of quantitative proteomic analysis to identify signatures of defects in specific cellular pathways in fibroblasts from subjects with inherited disease.

CLIC, a tool for expanding biological pathways based on co-expression across thousands of datasets

In recent years, there has been a huge rise in the number of publicly available transcriptional profiling datasets. These massive compendia comprise billions of measurements and provide a special opportunity to predict the function of unstudied genes based on co-expression to well-studied pathways. Such analyses can be very challenging, however, since biological pathways are modular and may exhibit co-expression only in specific contexts. To overcome these challenges we introduce CLIC, CLustering by Inferred Co-expression. CLIC accepts as input a pathway consisting of two or more genes. It then uses a Bayesian partition model to simultaneously partition the input gene set into coherent co-expressed modules (CEMs), while assigning the posterior probability for each dataset in support of each CEM. CLIC then expands each CEM by scanning the transcriptome for additional co-expressed genes, quantified by an integrated log-likelihood ratio (LLR) score weighted for each dataset. As a byproduct, CLIC automatically learns the conditions (datasets) within which a CEM is operative. We implemented CLIC using a compendium of 1774 mouse microarray datasets (28628 microarrays) or 1887 human microarray datasets (45158 microarrays). CLIC analysis reveals that of 910 canonical biological pathways, 30% consist of strongly co-expressed gene modules for which new members are predicted. For example, CLIC predicts a functional connection between protein C7orf55 (FMC1) and the mitochondrial ATP synthase complex that we have experimentally validated. CLIC is freely available at www.gene-clic.org. We anticipate that CLIC will be valuable both for revealing new components of biological pathways as well as the conditions in which they are active.

Comparative Analysis of Mitochondrial N-Termini from Mouse, Human, and Yeast

The majority of mitochondrial proteins are encoded in the nuclear genome, translated in the cytoplasm, and directed to the mitochondria by an N-terminal presequence that is cleaved upon import. Recently, N-proteome catalogs have been generated for mitochondria from yeast and from human U937 cells. Here, we applied the subtiligase method to determine N-termini for 327 proteins in mitochondria isolated from mouse liver and kidney. Comparative analysis between mitochondrial N-termini from mouse, human, and yeast proteins shows that whereas presequences are poorly conserved at the sequence level, other presequence properties are extremely conserved, including a length of ∼20-60 amino acids, a net charge between +3 to +6, and the presence of stabilizing amino acids at the N-terminus of mature proteins that follow the N-end rule from bacteria. As in yeast, ∼80% of mouse presequence cleavage sites match canonical motifs for three mitochondrial peptidases (MPP, Icp55, and Oct1), whereas the remainder do not match any known peptidase motifs. We show that mature mitochondrial proteins often exist with a spectrum of N-termini, consistent with a model of multiple cleavage events by MPP and Icp55. In addition to analysis of canonical targeting presequences, our N-terminal dataset allows the exploration of other cleavage events and provides support for polypeptide cleavage into two distinct enzymes (Hsd17b4), protein cleavages key for signaling (Oma1, Opa1, Htra2, Mavs, and Bcs2l13), and in several cases suggests novel protein isoforms (Scp2, Acadm, Adck3, Hsdl2, Dlst, and Ogdh). We present an integrated catalog of mammalian mitochondrial N-termini that can be used as a community resource to investigate individual proteins, to elucidate mechanisms of mammalian mitochondrial processing, and to allow researchers to engineer tags distally to the presequence cleavage.

A Genome-wide CRISPR Death Screen Identifies Genes Essential for Oxidative Phosphorylation

Oxidative phosphorylation (OXPHOS) is the major pathway for ATP production in humans. Deficiencies in OXPHOS can arise from mutations in either mitochondrial or nuclear genomes and comprise the largest collection of inborn errors of metabolism. At present we lack a complete catalog of human genes and pathways essential for OXPHOS. Here we introduce a genome-wide CRISPR "death screen" that actively selects dying cells to reveal human genes required for OXPHOS, inspired by the classic observation that human cells deficient in OXPHOS survive in glucose but die in galactose. We report 191 high-confidence hits essential for OXPHOS, including 72 underlying known OXPHOS diseases. Our screen reveals a functional module consisting of NGRN, WBSCR16, RPUSD3, RPUSD4, TRUB2, and FASTKD2 that regulates the mitochondrial 16S rRNA and intra-mitochondrial translation. Our work yields a rich catalog of genes required for OXPHOS and, more generally, demonstrates the power of death screening for functional genomic analysis.

MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins

Mitochondria are complex organelles that house essential pathways involved in energy metabolism, ion homeostasis, signalling and apoptosis. To understand mitochondrial pathways in health and disease, it is crucial to have an accurate inventory of the organelle's protein components. In 2008, we made substantial progress toward this goal by performing in-depth mass spectrometry of mitochondria from 14 organs, epitope tagging/microscopy and Bayesian integration to assemble MitoCarta (www.broadinstitute.org/pubs/MitoCarta): an inventory of genes encoding mitochondrial-localized proteins and their expression across 14 mouse tissues. Using the same strategy we have now reconstructed this inventory separately for human and for mouse based on (i) improved gene transcript models, (ii) updated literature curation, including results from proteomic analyses of mitochondrial sub-compartments, (iii) improved homology mapping and (iv) updated versions of all seven original data sets. The updated human MitoCarta2.0 consists of 1158 human genes, including 918 genes in the original inventory as well as 240 additional genes. The updated mouse MitoCarta2.0 consists of 1158 genes, including 967 genes in the original inventory plus 191 additional genes. The improved MitoCarta 2.0 inventory provides a molecular framework for system-level analysis of mammalian mitochondria.

Expansion of biological pathways based on evolutionary inference

The availability of diverse genomes makes it possible to predict gene function based on shared evolutionary history. This approach can be challenging, however, for pathways whose components do not exhibit a shared history but rather consist of distinct "evolutionary modules." We introduce a computational algorithm, clustering by inferred models of evolution (CLIME), which inputs a eukaryotic species tree, homology matrix, and pathway (gene set) of interest. CLIME partitions the gene set into disjoint evolutionary modules, simultaneously learning the number of modules and a tree-based evolutionary history that defines each module. CLIME then expands each module by scanning the genome for new components that likely arose under the inferred evolutionary model. Application of CLIME to ∼1,000 annotated human pathways and to the proteomes of yeast, red algae, and malaria reveals unanticipated evolutionary modularity and coevolving components. CLIME is freely available and should become increasingly powerful with the growing wealth of eukaryotic genomes.

Mutation in the novel nuclear-encoded mitochondrial protein CHCHD10 in a family with autosomal dominant mitochondrial myopathy

Mitochondrial myopathies belong to a larger group of systemic diseases caused by morphological or biochemical abnormalities of mitochondria. Mitochondrial disorders can be caused by mutations in either the mitochondrial or nuclear genome. Only 5% of all mitochondrial disorders are autosomal dominant. We analyzed DNA from members of the previously reported Puerto Rican kindred with an autosomal dominant mitochondrial myopathy (Heimann-Patterson et al. 1997). Linkage analysis suggested a putative locus on the pericentric region of the long arm of chromosome 22 (22q11). Using the tools of integrative genomics, we established chromosome 22 open reading frame 16 (C22orf16) (later designated as CHCHD10) as the only high-scoring mitochondrial candidate gene in our minimal candidate region. Sequence analysis revealed a double-missense mutation (R15S and G58R) in cis in CHCHD10 which encodes a coiled coil-helix-coiled coil-helix protein of unknown function. These two mutations completely co-segregated with the disease phenotype and were absent in 1,481 Caucasian and 80 Hispanic (including 32 Puerto Rican) controls. Expression profiling showed that CHCHD10 is enriched in skeletal muscle. Mitochondrial localization of the CHCHD10 protein was confirmed using immunofluorescence in cells expressing either wild-type or mutant CHCHD10. We found that the expression of the G58R, but not the R15S, mutation induced mitochondrial fragmentation. Our findings identify a novel gene causing mitochondrial myopathy, thereby expanding the spectrum of mitochondrial myopathies caused by nuclear genes. Our findings also suggest a role for CHCHD10 in the morphologic remodeling of the mitochondria.

Next generation sequencing with copy number variant detection expands the phenotypic spectrum of HSD17B4-deficiency

Background

D-bifunctional protein deficiency, caused by recessive mutations in HSD17B4, is a severe, infantile-onset disorder of peroxisomal fatty acid oxidation. Few affected patients survive past two years of age. Compound heterozygous mutations in HSD17B4 have also been reported in two sisters diagnosed with Perrault syndrome (MIM # 233400), who presented in adolescence with ovarian dysgenesis, hearing loss, and ataxia.

Case presentation

An adult male presented with cerebellar ataxia, peripheral neuropathy, hearing loss, and azoospermia. The clinical presentation, in combination with biochemical findings in serum, urine, and muscle biopsy, suggested a mitochondrial disorder. Commercial genetic testing of 18 ataxia and mitochondrial disease genes was negative. Targeted exome sequencing followed by analysis of single nucleotide variants and small insertions/deletions failed to reveal a genetic basis of disease. Application of a computational algorithm to infer copy number variants (CNVs) from exome data revealed a heterozygous 12 kb deletion of exons 10–13 of HSD17B4 that was compounded with a rare missense variant (p.A196V) at a highly conserved residue. Retrospective review of patient records revealed mildly elevated ratios of pristanic:phytanic acid and arachidonic:docosahexaenoic acid, consistent with dysfunctional peroxisomal fatty acid oxidation.

Conclusion

Our case expands the phenotypic spectrum of HSD17B4-deficiency, representing the first male case reported with infertility. Furthermore, it points to crosstalk between mitochondria and peroxisomes in HSD17B4-deficiency and Perrault syndrome.

CLYBL is a polymorphic human enzyme with malate synthase and β-methylmalate synthase activity

CLYBL is a human mitochondrial enzyme of unknown function that is found in multiple eukaryotic taxa and conserved to bacteria. The protein is expressed in the mitochondria of all mammalian organs, with highest expression in brown fat and kidney. Approximately 5% of all humans harbor a premature stop polymorphism in CLYBL that has been associated with reduced levels of circulating vitamin B12. Using comparative genomics, we now show that CLYBL is strongly co-expressed with and co-evolved specifically with other components of the mitochondrial B12 pathway. We confirm that the premature stop polymorphism in CLYBL leads to a loss of protein expression. To elucidate the molecular function of CLYBL, we used comparative operon analysis, structural modeling and enzyme kinetics. We report that CLYBL encodes a malate/β-methylmalate synthase, converting glyoxylate and acetyl-CoA to malate, or glyoxylate and propionyl-CoA to β-methylmalate. Malate synthases are best known for their established role in the glyoxylate shunt of plants and lower organisms and are traditionally described as not occurring in humans. The broader role of a malate/β-methylmalate synthase in human physiology and its mechanistic link to vitamin B12 metabolism remain unknown.

EMRE is an essential component of the mitochondrial calcium uniporter complex

The mitochondrial uniporter is a highly selective calcium channel in the organelle's inner membrane. Its molecular components include the EF-hand-containing calcium-binding proteins mitochondrial calcium uptake 1 (MICU1) and MICU2 and the pore-forming subunit mitochondrial calcium uniporter (MCU). We sought to achieve a full molecular characterization of the uniporter holocomplex (uniplex). Quantitative mass spectrometry of affinity-purified uniplex recovered MICU1 and MICU2, MCU and its paralog MCUb, and essential MCU regulator (EMRE), a previously uncharacterized protein. EMRE is a 10-kilodalton, metazoan-specific protein with a single transmembrane domain. In its absence, uniporter channel activity was lost despite intact MCU expression and oligomerization. EMRE was required for the interaction of MCU with MICU1 and MICU2. Hence, EMRE is essential for in vivo uniporter current and additionally bridges the calcium-sensing role of MICU1 and MICU2 with the calcium-conducting role of MCU.

Macrocytic anemia and mitochondriopathy resulting from a defect in sideroflexin 4

We used exome sequencing to identify mutations in sideroflexin 4 (SFXN4) in two children with mitochondrial disease (the more severe case also presented with macrocytic anemia). SFXN4 is an uncharacterized mitochondrial protein that localizes to the mitochondrial inner membrane. sfxn4 knockdown in zebrafish recapitulated the mitochondrial respiratory defect observed in both individuals and the macrocytic anemia with megaloblastic features of the more severe case. In vitro and in vivo complementation studies with fibroblasts from the affected individuals and zebrafish demonstrated the requirement of SFXN4 for mitochondrial respiratory homeostasis and erythropoiesis. Our findings establish mutations in SFXN4 as a cause of mitochondriopathy and macrocytic anemia.

Comparative RNA editing in autistic and neurotypical cerebella

Adenosine-to-inosine (A-to-I) RNA editing is a neurodevelopmentally regulated epigenetic modification shown to modulate complex behavior in animals. Little is known about human A-to-I editing, but it is thought to constitute one of many molecular mechanisms connecting environmental stimuli and behavioral outputs. Thus, comprehensive exploration of A-to-I RNA editing in human brains may shed light on gene-environment interactions underlying complex behavior in health and disease. Synaptic function is a main target of A-to-I editing, which can selectively recode key amino acids in synaptic genes, directly altering synaptic strength and duration in response to environmental signals. Here, we performed a high-resolution survey of synaptic A-to-I RNA editing in a human population, and examined how it varies in autism, a neurodevelopmental disorder in which synaptic abnormalities are a common finding. Using ultra-deep (>1000 × ) sequencing, we quantified the levels of A-to-I editing of 10 synaptic genes in postmortem cerebella from 14 neurotypical and 11 autistic individuals. A high dynamic range of editing levels was detected across individuals and editing sites, from 99.6% to below detection limits. In most sites, the extreme ends of the population editing distributions were individuals with autism. Editing was correlated with isoform usage, clusters of correlated sites were identified, and differential editing patterns examined. Finally, a dysfunctional form of the editing enzyme adenosine deaminase acting on RNA B1 was found more commonly in postmortem cerebella from individuals with autism. These results provide a population-level, high-resolution view of A-to-I RNA editing in human cerebella and suggest that A-to-I editing of synaptic genes may be informative for assessing the epigenetic risk for autism.

Targeted exome sequencing of suspected mitochondrial disorders

OBJECTIVE

To evaluate the utility of targeted exome sequencing for the molecular diagnosis of mitochondrial disorders, which exhibit marked phenotypic and genetic heterogeneity.

METHODS

We considered a diverse set of 102 patients with suspected mitochondrial disorders based on clinical, biochemical, and/or molecular findings, and whose disease ranged from mild to severe, with varying age at onset. We sequenced the mitochondrial genome (mtDNA) and the exons of 1,598 nuclear-encoded genes implicated in mitochondrial biology, mitochondrial disease, or monogenic disorders with phenotypic overlap. We prioritized variants likely to underlie disease and established molecular diagnoses in accordance with current clinical genetic guidelines.

RESULTS

Targeted exome sequencing yielded molecular diagnoses in established disease loci in 22% of cases, including 17 of 18 (94%) with prior molecular diagnoses and 5 of 84 (6%) without. The 5 new diagnoses implicated 2 genes associated with canonical mitochondrial disorders (NDUFV1, POLG2), and 3 genes known to underlie other neurologic disorders (DPYD, KARS, WFS1), underscoring the phenotypic and biochemical overlap with other inborn errors. We prioritized variants in an additional 26 patients, including recessive, X-linked, and mtDNA variants that were enriched 2-fold over background and await further support of pathogenicity. In one case, we modeled patient mutations in yeast to provide evidence that recessive mutations in ATP5A1 can underlie combined respiratory chain deficiency.

CONCLUSION

The results demonstrate that targeted exome sequencing is an effective alternative to the sequential testing of mtDNA and individual nuclear genes as part of the investigation of mitochondrial disease. Our study underscores the ongoing challenge of variant interpretation in the clinical setting.

MPV17 Mutations Causing Adult-Onset Multisystemic Disorder With Multiple Mitochondrial DNA Deletions

OBJECTIVE To identify the cause of an adult-onset multisystemic disease with multiple deletions of mitochondrial DNA (mtDNA). DESIGN Case report. SETTING University hospitals. PATIENT A 65-year-old man with axonal sensorimotor peripheral neuropathy, ptosis, ophthalmoparesis, diabetes mellitus, exercise intolerance, steatohepatopathy, depression, parkinsonism, and gastrointestinal dysmotility. RESULTS Skeletal muscle biopsy revealed ragged-red and cytochrome- c oxidase-deficient fibers, and Southern blot analysis showed multiple mtDNA deletions. No deletions were detected in fibroblasts, and the results of quantitative polymerase chain reaction showed that the amount of mtDNA was normal in both muscle and fibroblasts. Exome sequencing using a mitochondrial library revealed compound heterozygous MPV17 mutations (p.LysMet88-89MetLeu and p.Leu143*), a novel cause of mtDNA multiple deletions. CONCLUSIONS In addition to causing juvenile-onset disorders with mtDNA depletion, MPV17 mutations can cause adult-onset multisystemic disease with multiple mtDNA deletions.

Next-generation sequencing reveals DGUOK mutations in adult patients with mitochondrial DNA multiple deletions

The molecular diagnosis of mitochondrial disorders still remains elusive in a large proportion of patients, but advances in next generation sequencing are significantly improving our chances to detect mutations even in sporadic patients. Syndromes associated with mitochondrial DNA multiple deletions are caused by different molecular defects resulting in a wide spectrum of predominantly adult-onset clinical presentations, ranging from progressive external ophthalmoplegia to multi-systemic disorders of variable severity. The mutations underlying these conditions remain undisclosed in half of the affected subjects. We applied next-generation sequencing of known mitochondrial targets (MitoExome) to probands presenting with adult-onset mitochondrial myopathy and harbouring mitochondrial DNA multiple deletions in skeletal muscle. We identified autosomal recessive mutations in the DGUOK gene (encoding mitochondrial deoxyguanosine kinase), which has previously been associated with an infantile hepatocerebral form of mitochondrial DNA depletion. Mutations in DGUOK occurred in five independent subjects, representing 5.6% of our cohort of patients with mitochondrial DNA multiple deletions, and impaired both muscle DGUOK activity and protein stability. Clinical presentations were variable, including mitochondrial myopathy with or without progressive external ophthalmoplegia, recurrent rhabdomyolysis in a young female who had received a liver transplant at 9 months of age and adult-onset lower motor neuron syndrome with mild cognitive impairment. These findings reinforce the concept that mutations in genes involved in deoxyribonucleotide metabolism can cause diverse clinical phenotypes and suggest that DGUOK should be screened in patients harbouring mitochondrial DNA deletions in skeletal muscle.

Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing

Advances in next-generation sequencing (NGS) promise to facilitate diagnosis of inherited disorders. Although in research settings NGS has pinpointed causal alleles using segregation in large families, the key challenge for clinical diagnosis is application to single individuals. To explore its diagnostic use, we performed targeted NGS in 42 unrelated infants with clinical and biochemical evidence of mitochondrial oxidative phosphorylation disease. These devastating mitochondrial disorders are characterized by phenotypic and genetic heterogeneity, with more than 100 causal genes identified to date. We performed "MitoExome" sequencing of the mitochondrial DNA (mtDNA) and exons of ~1000 nuclear genes encoding mitochondrial proteins and prioritized rare mutations predicted to disrupt function. Because patients and healthy control individuals harbored a comparable number of such heterozygous alleles, we could not prioritize dominant-acting genes. However, patients showed a fivefold enrichment of genes with two such mutations that could underlie recessive disease. In total, 23 of 42 (55%) patients harbored such recessive genes or pathogenic mtDNA variants. Firm diagnoses were enabled in 10 patients (24%) who had mutations in genes previously linked to disease. Thirteen patients (31%) had mutations in nuclear genes not previously linked to disease. The pathogenicity of two such genes, NDUFB3 and AGK, was supported by complementation studies and evidence from multiple patients, respectively. The results underscore the potential and challenges of deploying NGS in clinical settings.

Evolutionary diversity of the mitochondrial calcium uniporter

Calcium uptake into mitochondria occurs via a recently identified ion channel called the uniporter. Here, we characterize the phylogenomic distribution of the uniporter's membrane-spanning pore subunit (MCU) and regulatory partner (MICU1). Homologs of both components tend to co-occur in all major branches of eukaryotic life, but both have been lost along certain protozoan and fungal lineages. Several bacterial genomes also contain putative MCU homologs that may represent prokaryotic calcium channels. The analyses indicate that the uniporter may have been an early feature of mitochondria.

Atypical case of Wolfram syndrome revealed through targeted exome sequencing in a patient with suspected mitochondrial disease

BACKGROUND

Mitochondrial diseases comprise a diverse set of clinical disorders that affect multiple organ systems with varying severity and age of onset. Due to their clinical and genetic heterogeneity, these diseases are difficult to diagnose. We have developed a targeted exome sequencing approach to improve our ability to properly diagnose mitochondrial diseases and apply it here to an individual patient. Our method targets mitochondrial DNA (mtDNA) and the exons of 1,600 nuclear genes involved in mitochondrial biology or Mendelian disorders with multi-system phenotypes, thereby allowing for simultaneous evaluation of multiple disease loci.

CASE PRESENTATION

Targeted exome sequencing was performed on a patient initially suspected to have a mitochondrial disorder. The patient presented with diabetes mellitus, diffuse brain atrophy, autonomic neuropathy, optic nerve atrophy, and a severe amnestic syndrome. Further work-up revealed multiple heteroplasmic mtDNA deletions as well as profound thiamine deficiency without a clear nutritional cause. Targeted exome sequencing revealed a homozygous c.1672C > T (p.R558C) missense mutation in exon 8 of WFS1 that has previously been reported in a patient with Wolfram syndrome.

CONCLUSION

This case demonstrates how clinical application of next-generation sequencing technology can enhance the diagnosis of patients suspected to have rare genetic disorders. Furthermore, the finding of unexplained thiamine deficiency in a patient with Wolfram syndrome suggests a potential link between WFS1 biology and thiamine metabolism that has implications for the clinical management of Wolfram syndrome patients.

Mutations in MTFMT underlie a human disorder of formylation causing impaired mitochondrial translation

The metazoan mitochondrial translation machinery is unusual in having a single tRNA(Met) that fulfills the dual role of the initiator and elongator tRNA(Met). A portion of the Met-tRNA(Met) pool is formylated by mitochondrial methionyl-tRNA formyltransferase (MTFMT) to generate N-formylmethionine-tRNA(Met) (fMet-tRNA(met)), which is used for translation initiation; however, the requirement of formylation for initiation in human mitochondria is still under debate. Using targeted sequencing of the mtDNA and nuclear exons encoding the mitochondrial proteome (MitoExome), we identified compound heterozygous mutations in MTFMT in two unrelated children presenting with Leigh syndrome and combined OXPHOS deficiency. Patient fibroblasts exhibit severe defects in mitochondrial translation that can be rescued by exogenous expression of MTFMT. Furthermore, patient fibroblasts have dramatically reduced fMet-tRNA(Met) levels and an abnormal formylation profile of mitochondrially translated COX1. Our findings demonstrate that MTFMT is critical for efficient human mitochondrial translation and reveal a human disorder of Met-tRNA(Met) formylation.

The molecular basis of human complex I deficiency

Disorders of oxidative phosphorylation (OXPHOS) have a birth prevalence of ∼1/5,000 and are the most common inborn errors of metabolism. The most common OXPHOS disorder is complex I deficiency. Patients with complex I deficiency present with variable symptoms, such as muscle weakness, cardiomyopathy, developmental delay or regression, blindness, seizures, failure to thrive, liver dysfunction or ataxia. Molecular diagnosis of patients with complex I deficiency is a challenging task due to the clinical heterogeneity of patients and the large number of candidate disease genes, both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded. In this review, we have thoroughly surveyed the literature to identify 149 patients described with both isolated complex I deficiency and pathogenic mutations within nuclear genes. In total, 115 different pathogenic mutations have been reported in 22 different nuclear genes encoding complex I subunits or assembly factors, highlighting the allelic and locus heterogeneity of this disorder. Missense mutations predominate in genes encoding core subunits and some assembly factors while null-type mutations are common in the genes encoding supernumerary subunits and other assembly factors. Despite developments in molecular technology, many patients do not receive molecular diagnosis and no gene has yet been identified that accounts for more than 5% of cases, suggesting that there are likely many disease genes that await discovery.

The mitochondrial proteome and human disease

For nearly three decades, the sequence of the human mitochondrial genome (mtDNA) has provided a molecular framework for understanding maternally inherited diseases. However, the vast majority of human mitochondrial disorders are caused by nuclear genome defects, which is not surprising since the mtDNA encodes only 13 proteins. Advances in genomics, mass spectrometry, and computation have only recently made it possible to systematically identify the complement of over 1,000 proteins that comprise the mammalian mitochondrial proteome. Here, we review recent progress in characterizing the mitochondrial proteome and highlight insights into its complexity, tissue heterogeneity, evolutionary origins, and biochemical versatility. We then discuss how this proteome is being used to discover the genetic basis of respiratory chain disorders as well as to expand our definition of mitochondrial disease. Finally, we explore future prospects and challenges for using the mitochondrial proteome as a foundation for systems analysis of the organelle.

FOXRED1, encoding an FAD-dependent oxidoreductase complex-I-specific molecular chaperone, is mutated in infantile-onset mitochondrial encephalopathy

Complex I is the first and largest enzyme in the respiratory chain and is located in the inner mitochondrial membrane. Complex I deficiency is the most commonly reported mitochondrial disorder presenting in childhood, but the molecular basis of most cases remains elusive. We describe a patient with complex I deficiency caused by mutation of the molecular chaperone FOXRED1. A combined homozygosity mapping and bioinformatics approach in a consanguineous Iranian-Jewish pedigree led to the identification of a homozygous mutation in FOXRED1 in a child who presented with infantile-onset encephalomyopathy. Silencing of FOXRED1 in human fibroblasts resulted in reduced complex I steady-state levels and activity, while lentiviral-mediated FOXRED1 transgene expression rescued complex I deficiency in the patient fibroblasts. This FAD-dependent oxidoreductase, which has never previously been associated with human disease, is now shown to be a complex I-specific molecular chaperone. The discovery of the c.1054C>T; p.R352W mutation in the FOXRED1 gene is a further contribution towards resolving the complex puzzle of the genetic basis of human mitochondrial disease.

High-throughput, pooled sequencing identifies mutations in NUBPL and FOXRED1 in human complex I deficiency

Discovering the molecular basis of mitochondrial respiratory chain disease is challenging given the large number of both mitochondrial and nuclear genes that are involved. We report a strategy of focused candidate gene prediction, high-throughput sequencing and experimental validation to uncover the molecular basis of mitochondrial complex I disorders. We created seven pools of DNA from a cohort of 103 cases and 42 healthy controls and then performed deep sequencing of 103 candidate genes to identify 151 rare variants that were predicted to affect protein function. We established genetic diagnoses in 13 of 60 previously unsolved cases using confirmatory experiments, including cDNA complementation to show that mutations in NUBPL and FOXRED1 can cause complex I deficiency. Our study illustrates how large-scale sequencing, coupled with functional prediction and experimental validation, can be used to identify causal mutations in individual cases.

A mitochondrial protein compendium elucidates complex I disease biology

Mitochondria are complex organelles whose dysfunction underlies a broad spectrum of human diseases. Identifying all of the proteins resident in this organelle and understanding how they integrate into pathways represent major challenges in cell biology. Toward this goal, we performed mass spectrometry, GFP tagging, and machine learning to create a mitochondrial compendium of 1098 genes and their protein expression across 14 mouse tissues. We link poorly characterized proteins in this inventory to known mitochondrial pathways by virtue of shared evolutionary history. Using this approach, we predict 19 proteins to be important for the function of complex I (CI) of the electron transport chain. We validate a subset of these predictions using RNAi, including C8orf38, which we further show harbors an inherited mutation in a lethal, infantile CI deficiency. Our results have important implications for understanding CI function and pathogenesis and, more generally, illustrate how our compendium can serve as a foundation for systematic investigations of mitochondria.

Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae

The aspergilli comprise a diverse group of filamentous fungi spanning over 200 million years of evolution. Here we report the genome sequence of the model organism Aspergillus nidulans, and a comparative study with Aspergillus fumigatus, a serious human pathogen, and Aspergillus oryzae, used in the production of sake, miso and soy sauce. Our analysis of genome structure provided a quantitative evaluation of forces driving long-term eukaryotic genome evolution. It also led to an experimentally validated model of mating-type locus evolution, suggesting the potential for sexual reproduction in A. fumigatus and A. oryzae. Our analysis of sequence conservation revealed over 5,000 non-coding regions actively conserved across all three species. Within these regions, we identified potential functional elements including a previously uncharacterized TPP riboswitch and motifs suggesting regulation in filamentous fungi by Puf family genes. We further obtained comparative and experimental evidence indicating widespread translational regulation by upstream open reading frames. These results enhance our understanding of these widely studied fungi as well as provide new insight into eukaryotic genome evolution and gene regulation.

The genome sequence of the filamentous fungus Neurospora crassa

Neurospora crassa is a central organism in the history of twentieth-century genetics, biochemistry and molecular biology. Here, we report a high-quality draft sequence of the N. crassa genome. The approximately 40-megabase genome encodes about 10,000 protein-coding genes--more than twice as many as in the fission yeast Schizosaccharomyces pombe and only about 25% fewer than in the fruitfly Drosophila melanogaster. Analysis of the gene set yields insights into unexpected aspects of Neurospora biology including the identification of genes potentially associated with red light photobiology, genes implicated in secondary metabolism, and important differences in Ca2+ signalling as compared with plants and animals. Neurospora possesses the widest array of genome defence mechanisms known for any eukaryotic organism, including a process unique to fungi called repeat-induced point mutation (RIP). Genome analysis suggests that RIP has had a profound impact on genome evolution, greatly slowing the creation of new genes through genomic duplication and resulting in a genome with an unusually low proportion of closely related genes.