Mutations in PDLIM5 are rare in dilated cardiomyopathy but are emerging as potential disease modifiers.

BACKGROUND: A causal genetic mutation is found in 40% of families with dilated cardiomyopathy (DCM), leaving a large percentage of families genetically unsolved. This prevents adequate counseling and clear recommendations in these families. We aim to identify novel genes or modifiers associated with DCM. METHODS: We performed computational ranking of human genes based on coexpression with a predefined set of genes known to be associated with DCM, which allowed us to prioritize gene candidates for their likelihood of being involved in DCM. Top candidates will be checked for variants in the available whole-exome sequencing data of 142 DCM patients. RNA was isolated from cardiac biopsies to investigate gene expression. RESULTS: PDLIM5 was classified as the top candidate. An interesting heterozygous variant (189_190delinsGG) was found in a DCM patient with a known pathogenic truncating TTN-variant. The PDLIM5 loss-of-function (LoF) variant affected all cardiac-specific isoforms of PDLIM5 and no LoF variants were detected in the same region in a control cohort of 26,000 individuals. RNA expression of PDLIM5 and its direct interactors (MYOT, LDB3, and MYOZ2) was increased in cardiac tissue of this patient, indicating a possible compensatory mechanism. The PDLIM5 variant cosegregated with the TTN-variant and the phenotype, leading to a high disease penetrance in this family. A second patient was an infant with a homozygous 10 kb-deletion of exon 2 in PDLIM5 resulting in early-onset cardiac disease, showing the importance of PDLIM5 in cardiac function. CONCLUSIONS: Heterozygous PDLIM5 variants are rare and therefore will not have a major contribution in DCM. Although they likely play a role in disease development as this gene plays a major role in contracting cardiomyocytes and homozygous variants lead to early-onset cardiac disease. Other environmental and/or genetic factors are probably necessary to unveil the cardiac phenotype in PDLIM5 mutation carriers.

CiliaCarta: An integrated and validated compendium of ciliary genes.

The cilium is an essential organelle at the surface of mammalian cells whose dysfunction causes a wide range of genetic diseases collectively called ciliopathies. The current rate at which new ciliopathy genes are identified suggests that many ciliary components remain undiscovered. We generated and rigorously analyzed genomic, proteomic, transcriptomic and evolutionary data and systematically integrated these using Bayesian statistics into a predictive score for ciliary function. This resulted in 285 candidate ciliary genes. We generated independent experimental evidence of ciliary associations for 24 out of 36 analyzed candidate proteins using multiple cell and animal model systems (mouse, zebrafish and nematode) and techniques. For example, we show that OSCP1, which has previously been implicated in two distinct non-ciliary processes, causes ciliogenic and ciliopathy-associated tissue phenotypes when depleted in zebrafish. The candidate list forms the basis of CiliaCarta, a comprehensive ciliary compendium covering 956 genes. The resource can be used to objectively prioritize candidate genes in whole exome or genome sequencing of ciliopathy patients and can be accessed at http://bioinformatics.bio.uu.nl/john/syscilia/ciliacarta/.

Whole Exome Sequencing Is the Preferred Strategy to Identify the Genetic Defect in Patients With a Probable or Possible Mitochondrial Cause.

Mitochondrial disorders, characterized by clinical symptoms and/or OXPHOS deficiencies, are caused by pathogenic variants in mitochondrial genes. However, pathogenic variants in some of these genes can lead to clinical manifestations which overlap with other neuromuscular diseases, which can be caused by pathogenic variants in non-mitochondrial genes as well. Mitochondrial pathogenic variants can be found in the mitochondrial DNA (mtDNA) or in any of the 1,500 nuclear genes with a mitochondrial function. We have performed a two-step next-generation sequencing approach in a cohort of 117 patients, mostly children, in whom a mitochondrial disease-cause could likely or possibly explain the phenotype. A total of 86 patients had a mitochondrial disorder, according to established clinical and biochemical criteria. The other 31 patients had neuromuscular symptoms, where in a minority a mitochondrial genetic cause is present, but a non-mitochondrial genetic cause is more likely. All patients were screened for pathogenic variants in the mtDNA and, if excluded, analyzed by whole exome sequencing (WES). Variants were filtered for being pathogenic and compatible with an autosomal or X-linked recessive mode of inheritance in families with multiple affected siblings and/or consanguineous parents. Non-consanguineous families with a single patient were additionally screened for autosomal and X-linked dominant mutations in a predefined gene-set. We identified causative pathogenic variants in the mtDNA in 20% of the patient-cohort, and in nuclear genes in 49%, implying an overall yield of 68%. We identified pathogenic variants in mitochondrial and non-mitochondrial genes in both groups with, obviously, a higher number of mitochondrial genes affected in mitochondrial disease patients. Furthermore, we show that 31% of the disease-causing genes in the mitochondrial patient group were not included in the MitoCarta database, and therefore would have been missed with MitoCarta based gene-panels. We conclude that WES is preferable to panel-based approaches for both groups of patients, as the mitochondrial gene-list is not complete and mitochondrial symptoms can be secondary. Also, clinically and genetically heterogeneous disorders would require sequential use of multiple different gene panels. We conclude that WES is a comprehensive and unbiased approach to establish a genetic diagnosis in these patients, able to resolve multi-genic disease-causes.

Genetic defects in mtDNA-encoded protein translation cause pediatric, mitochondrial cardiomyopathy with early-onset brain disease.

This study aims to identify gene defects in pediatric cardiomyopathy and early-onset brain disease with oxidative phosphorylation (OXPHOS) deficiencies. We applied whole-exome sequencing in three patients with pediatric cardiomyopathy and early-onset brain disease with OXPHOS deficiencies. The brain pathology was studied by MRI analysis. In consanguineous patient 1, we identified a homozygous intronic variant (c.850-3A > G) in the QRSL1 gene, which was predicted to cause abnormal splicing. The variant segregated with the disease and affected the protein function, which was confirmed by complementation studies, restoring OXPHOS function only with wild-type QRSL1. Patient 2 was compound heterozygous for two novel affected and disease-causing variants (c.[253G > A];[938G > A]) in the MTO1 gene. In patient 3, we detected one unknown affected and disease-causing variants (c.2872C > T) and one known disease-causing variant (c.1774C > T) in the AARS2 gene. The c.1774C > T variant was present in the paternal copy of the AARS2 gene, the c.2872C > T in the maternal copy. All genes were involved in translation of mtDNA-encoded proteins. Defects in mtDNA-encoded protein translation lead to severe pediatric cardiomyopathy and brain disease with OXPHOS abnormalities. This suggests that the heart and brain are particularly sensitive to defects in mitochondrial protein synthesis during late embryonic or early postnatal development, probably due to the massive mitochondrial biogenesis occurring at that stage. If both the heart and brain are involved, the prognosis is poor with a likely fatal outcome at young age.

Selection and Characterization of Palmitic Acid Responsive Patients with an OXPHOS Complex I Defect.

Mitochondrial disorders are genetically and clinically heterogeneous, mainly affecting high energy-demanding organs due to impaired oxidative phosphorylation (OXPHOS). Currently, effective treatments for OXPHOS defects, with complex I deficiency being the most prevalent, are not available. Yet, clinical practice has shown that some complex I deficient patients benefit from a high-fat or ketogenic diet, but it is unclear how these therapeutic diets influence mitochondrial function and more importantly, which complex I patients could benefit from such treatment. Dietary studies in a complex I deficient patient with exercise intolerance showed increased muscle endurance on a high-fat diet compared to a high-carbohydrate diet. We performed whole-exome sequencing to characterize the genetic defect. A pathogenic homozygous p.G212V missense mutation was identified in the TMEM126B gene, encoding an early assembly factor of complex I. A complementation study in fibroblasts confirmed that the p.G212V mutation caused the complex I deficiency. The mechanism turned out to be an incomplete assembly of the peripheral arm of complex I, leading to a decrease in the amount of mature complex I. The patient clinically improved on a high-fat diet, which was supported by the 25% increase in maximal OXPHOS capacity in TMEM126B defective fibroblast by the saturated fatty acid palmitic acid, whereas oleic acid did not have any effect in those fibroblasts. Fibroblasts of other patients with a characterized complex I gene defect were tested in the same way. Patient fibroblasts with complex I defects in NDUFS7 and NDUFAF5 responded to palmitic acid, whereas ACAD9, NDUFA12, and NDUFV2 defects were non-responding. Although the data are too limited to draw a definite conclusion on the mechanism, there is a tendency that protein defects involved in early assembly complexes, improve with palmitic acid, whereas proteins defects involved in late assembly, do not. Our data show at a clinical and biochemical level that a high fat diet can be beneficial for complex I patients and that our cell line assay will be an easy tool for the selection of patients, who might potentially benefit from this therapeutic diet.

A comprehensive strategy for exome-based preconception carrier screening.

PURPOSE: Whole-exome sequencing (WES) provides the possibility of genome-wide preconception carrier screening (PCS). Here, we propose a filter strategy to rapidly identify the majority of relevant pathogenic mutations. METHODS: Our strategy was developed using WES data from eight consanguineous and five fictive nonconsanguineous couples and was subsequently applied to 20 other fictive nonconsanguineous couples. Presumably pathogenic variants based on frequency and database annotations or generic characteristics and mutation type were selected in genes shared by the couple and in the female's X-chromosome. Unclassified variants were not included. RESULTS: This yielded an average of 29 (19-51) variants in genes shared by the consanguineous couples and 15 (6-30) shared by the nonconsanguineous couples. For X-linked variants, the numbers per female were 3 (1-5) and 1 (0-3), respectively. Remaining variants were verified manually. The majority were able to be quickly discarded, effectively leaving true pathogenic variants. CONCLUSION: We conclude that WES is applicable for PCS, both for consanguineous and nonconsanguineous couples, with the remaining number of variants being manageable in a clinical setting. The addition of gene panels for filtering was not favorable because it resulted in missing pathogenic variants. It is important to develop and continuously curate databases with pathogenic mutations to further increase the sensitivity of WES-based PCS.Genet Med advance online publication 27 October 2016.

Network topology of NaV1.7 mutations in sodium channel-related painful disorders.

BACKGROUND: Gain-of-function mutations in SCN9A gene that encodes the voltage-gated sodium channel NaV1.7 have been associated with a wide spectrum of painful syndromes in humans including inherited erythromelalgia, paroxysmal extreme pain disorder and small fibre neuropathy. These mutations change the biophysical properties of NaV1.7 channels leading to hyperexcitability of dorsal root ganglion nociceptors and pain symptoms. There is a need for better understanding of how gain-of-function mutations alter the atomic structure of Nav1.7. RESULTS: We used homology modeling to build an atomic model of NaV1.7 and a network-based theoretical approach, which can predict interatomic interactions and connectivity arrangements, to investigate how pain-related NaV1.7 mutations may alter specific interatomic bonds and cause connectivity rearrangement, compared to benign variants and polymorphisms. For each amino acid substitution, we calculated the topological parameters betweenness centrality (B ct ), degree (D), clustering coefficient (CC ct ), closeness (C ct ), and eccentricity (E ct ), and calculated their variation (Δ value = mutant value -WT value ). Pathogenic NaV1.7 mutations showed significantly higher variation of |ΔB ct | compared to benign variants and polymorphisms. Using the cut-off value ±0.26 calculated by receiver operating curve analysis, we found that ΔB ct correctly differentiated pathogenic NaV1.7 mutations from variants not causing biophysical abnormalities (nABN) and homologous SNPs (hSNPs) with 76% sensitivity and 83% specificity. CONCLUSIONS: Our in-silico analyses predict that pain-related pathogenic NaV1.7 mutations may affect the network topological properties of the protein and suggest |ΔB ct | value as a potential in-silico marker.

A homozygous FITM2 mutation causes a deafness-dystonia syndrome with motor regression and signs of ichthyosis and sensory neuropathy.

A consanguineous family from Pakistan was ascertained to have a novel deafness-dystonia syndrome with motor regression, ichthyosis-like features and signs of sensory neuropathy. By applying a combined strategy of linkage analysis and whole-exome sequencing in the presented family, a homozygous nonsense mutation, c.4G>T (p.Glu2*), in FITM2 was identified. FITM2 and its paralog FITM1 constitute an evolutionary conserved protein family involved in partitioning of triglycerides into cellular lipid droplets. Despite the role of FITM2 in neutral lipid storage and metabolism, no indications for lipodystrophy were observed in the affected individuals. In order to obtain independent evidence for the involvement of FITM2 in the human pathology, downregulation of the single Fitm ortholog, CG10671, in Drosophila melanogaster was pursued using RNA interference. Characteristics of the syndrome, including progressive locomotor impairment, hearing loss and disturbed sensory functions, were recapitulated in Drosophila, which supports the causative nature of the FITM2 mutation. Mutation-based genetic counseling can now be provided to the family and insight is obtained into the potential impact of genetic variation in FITM2.

Rapid Resolution of Blended or Composite Multigenic Disease in Infants by Whole-Exome Sequencing.

Whole-exome sequencing identified multiple genetic causes in 2 infants with heterogeneous disease. Three gene defects in the first patient explained all symptoms, but manifestations were overlapping (blended phenotype). Two gene defects in the second patient explained nonoverlapping symptoms (composite phenotype). Whole-exome sequencing rapidly and comprehensively resolves heterogeneous genetic disease.

Specific MRI Abnormalities Reveal Severe Perrault Syndrome due to CLPP Defects.

In establishing a genetic diagnosis in heterogeneous neurological disease, clinical characterization and whole exome sequencing (WES) go hand-in-hand. Clinical data are essential, not only to guide WES variant selection and define the clinical severity of a genetic defect but also to identify other patients with defects in the same gene. In an infant patient with sensorineural hearing loss, psychomotor retardation, and epilepsy, WES resulted in identification of a novel homozygous CLPP frameshift mutation (c.21delA). Based on the gene defect and clinical symptoms, the diagnosis Perrault syndrome type 3 (PRLTS3) was established. The patient's brain-MRI revealed specific abnormalities of the subcortical and deep cerebral white matter and the middle blade of the corpus callosum, which was used to identify similar patients in the Amsterdam brain-MRI database, containing over 3000 unclassified leukoencephalopathy cases. In three unrelated patients with similar MRI abnormalities the CLPP gene was sequenced, and in two of them novel missense mutations were identified together with a large deletion that covered part of the CLPP gene on the other allele. The severe neurological and MRI abnormalities in these young patients were due to the drastic impact of the CLPP mutations, correlating with the variation in clinical manifestations among previously reported patients. Our data show that similarity in brain-MRI patterns can be used to identify novel PRLTS3 patients, especially during early disease stages, when only part of the disease manifestations are present. This seems especially applicable to the severely affected cases in which CLPP function is drastically affected and MRI abnormalities are pronounced.

WeGET: predicting new genes for molecular systems by weighted co-expression.

We have developed the Weighted Gene Expression Tool and database (WeGET, http://weget.cmbi.umcn.nl) for the prediction of new genes of a molecular system by correlated gene expression. WeGET utilizes a compendium of 465 human and 560 murine gene expression datasets that have been collected from multiple tissues under a wide range of experimental conditions. It exploits this abundance of expression data by assigning a high weight to datasets in which the known genes of a molecular system are harmoniously up- and down-regulated. WeGET ranks new candidate genes by calculating their weighted co-expression with that system. A weighted rank is calculated for human genes and their mouse orthologs. Then, an integrated gene rank and p-value is computed using a rank-order statistic. We applied our method to predict novel genes that have a high degree of co-expression with Gene Ontology terms and pathways from KEGG and Reactome. For each query set we provide a list of predicted novel genes, computed weights for transcription datasets used and cell and tissue types that contributed to the final predictions. The performance for each query set is assessed by 10-fold cross-validation. Finally, users can use the WeGET to predict novel genes that co-express with a custom query set.

Integrative Genomics-Based Discovery of Novel Regulators of the Innate Antiviral Response.

The RIG-I-like receptor (RLR) pathway is essential for detecting cytosolic viral RNA to trigger the production of type I interferons (IFNα/ß) that initiate an innate antiviral response. Through systematic assessment of a wide variety of genomics data, we discovered 10 molecular signatures of known RLR pathway components that collectively predict novel members. We demonstrate that RLR pathway genes, among others, tend to evolve rapidly, interact with viral proteins, contain a limited set of protein domains, are regulated by specific transcription factors, and form a tightly connected interaction network. Using a Bayesian approach to integrate these signatures, we propose likely novel RLR regulators. RNAi knockdown experiments revealed a high prediction accuracy, identifying 94 genes among 187 candidates tested (~50%) that affected viral RNA-induced production of IFNß. The discovered antiviral regulators may participate in a wide range of processes that highlight the complexity of antiviral defense (e.g. MAP3K11, CDK11B, PSMA3, TRIM14, HSPA9B, CDC37, NUP98, G3BP1), and include uncharacterized factors (DDX17, C6orf58, C16orf57, PKN2, SNW1). Our validated RLR pathway list (http://rlr.cmbi.umcn.nl/), obtained using a combination of integrative genomics and experiments, is a new resource for innate antiviral immunity research.

Transcriptome analysis of complex I-deficient patients reveals distinct expression programs for subunits and assembly factors of the oxidative phosphorylation system.

BACKGROUND: Transcriptional control of mitochondrial metabolism is essential for cellular function. A better understanding of this process will aid the elucidation of mitochondrial disorders, in particular of the many genetically unsolved cases of oxidative phosphorylation (OXPHOS) deficiency. Yet, to date only few studies have investigated nuclear gene regulation in the context of OXPHOS deficiency. In this study we performed RNA sequencing of two control and two complex I-deficient patient cell lines cultured in the presence of compounds that perturb mitochondrial metabolism: chloramphenicol, AICAR, or resveratrol. We combined this with a comprehensive analysis of mitochondrial and nuclear gene expression patterns, co-expression calculations and transcription factor binding sites. RESULTS: Our analyses show that subsets of mitochondrial OXPHOS genes respond opposingly to chloramphenicol and AICAR, whereas the response of nuclear OXPHOS genes is less consistent between cell lines and treatments. Across all samples nuclear OXPHOS genes have a significantly higher co-expression with each other than with other genes, including those encoding mitochondrial proteins. We found no evidence for complex-specific mRNA expression regulation: subunits of different OXPHOS complexes are similarly (co-)expressed and regulated by a common set of transcription factors. However, we did observe significant differences between the expression of nuclear genes for OXPHOS subunits versus assembly factors, suggesting divergent transcription programs. Furthermore, complex I co-expression calculations identified 684 genes with a likely role in OXPHOS biogenesis and function. Analysis of evolutionarily conserved transcription factor binding sites in the promoters of these genes revealed almost all known OXPHOS regulators (including GABP, NRF1/2, SP1, YY1, E-box factors) and a set of novel candidates (ELK1, KLF7, SP4, EHF, ZNF143, and TEL2). CONCLUSIONS: OXPHOS genes share an expression program distinct from other genes encoding mitochondrial proteins, indicative of targeted nuclear regulation of a mitochondrial sub-process. Within the subset of OXPHOS genes we established a difference in expression between mitochondrial and nuclear genes, and between nuclear genes encoding subunits and assembly factors. Most transcription regulators of genes that co-express with complex I are well-established factors for OXPHOS biogenesis. For the remaining six factors we here suggest for the first time a link with transcription regulation in OXPHOS deficiency.

Olfactory receptor genes cooperate with protocadherin genes in human extreme obesity.

Worldwide, the incidence of obesity has increased dramatically over the past decades. More knowledge about the complex etiology of obesity is needed in order to find additional approaches for treatment and prevention. Investigating the exome sequencing data of 30 extremely obese subjects (BMI 45-65 kg/m(2)) shows that predicted damaging missense variants in olfactory receptor genes on chromosome 1q and rare predicted damaging variants in the protocadherin (PCDH) beta-cluster genes on chromosome 5q31, reported in our previous work, co-localize in subjects with extreme obesity. This implies a synergistic effect between genetic variation in these gene clusters in the predisposition to extreme obesity. Evidence for a general involvement of the olfactory transduction pathway on itself could not be found. Bioinformatic analysis indicates a specific involvement of the PCDH beta-cluster genes in controlling tissue development. Further mechanistic insight needs to await the identification of the ligands of the 1q olfactory receptors. Eventually, this may provide the possibility to manipulate food flavor in a way to reduce the risk of overeating and of extreme obesity in genetically predisposed subjects.

CLPB mutations cause 3-methylglutaconic aciduria, progressive brain atrophy, intellectual disability, congenital neutropenia, cataracts, movement disorder.

We studied a group of individuals with elevated urinary excretion of 3-methylglutaconic acid, neutropenia that can develop into leukemia, a neurological phenotype ranging from nonprogressive intellectual disability to a prenatal encephalopathy with progressive brain atrophy, movement disorder, cataracts, and early death. Exome sequencing of two unrelated individuals and subsequent Sanger sequencing of 16 individuals with an overlapping phenotype identified a total of 14 rare, predicted deleterious alleles in CLPB in 14 individuals from 9 unrelated families. CLPB encodes caseinolytic peptidase B homolog ClpB, a member of the AAA+ protein family. To evaluate the relevance of CLPB in the pathogenesis of this syndrome, we developed a zebrafish model and an in vitro assay to measure ATPase activity. Suppression of clpb in zebrafish embryos induced a central nervous system phenotype that was consistent with cerebellar and cerebral atrophy that could be rescued by wild-type, but not mutant, human CLPB mRNA. Consistent with these data, the loss-of-function effect of one of the identified variants (c.1222A>G [p.Arg408Gly]) was supported further by in vitro evidence with the mutant peptides abolishing ATPase function. Additionally, we show that CLPB interacts biochemically with ATP2A2, known to be involved in apoptotic processes in severe congenital neutropenia (SCN) 3 (Kostmann disease [caused by HAX1 mutations]). Taken together, mutations in CLPB define a syndrome with intellectual disability, congenital neutropenia, progressive brain atrophy, movement disorder, cataracts, and 3-methylglutaconic aciduria.

A mutation in the human CBP4 ortholog UQCC3 impairs complex III assembly, activity and cytochrome b stability.

Complex III (cytochrome bc1) is a protein complex of the mitochondrial inner membrane that transfers electrons from ubiquinol to cytochrome c. Its assembly requires the coordinated expression of mitochondrial-encoded cytochrome b and nuclear-encoded subunits and assembly factors. Complex III deficiency is a severe multisystem disorder caused by mutations in subunit genes or assembly factors. Sequence-profile-based orthology predicts C11orf83, hereafter named UQCC3, to be the ortholog of the fungal complex III assembly factor CBP4. We describe a homozygous c.59T>A missense mutation in UQCC3 from a consanguineous patient diagnosed with isolated complex III deficiency, displaying lactic acidosis, hypoglycemia, hypotonia and delayed development without dysmorphic features. Patient fibroblasts have reduced complex III activity and lower levels of the holocomplex and its subunits than controls. They have no detectable UQCC3 protein and have lower levels of cytochrome b protein. Furthermore, in patient cells, cytochrome b is absent from a high-molecular-weight complex III. UQCC3 is reduced in cells depleted for the complex III assembly factors UQCC1 and UQCC2. Conversely, absence of UQCC3 in patient cells does not affect UQCC1 and UQCC2. This suggests that UQCC3 functions in the complex III assembly pathway downstream of UQCC1 and UQCC2 and is consistent with what is known about the function of Cbp4 and of the fungal orthologs of UQCC1 and UQCC2, Cbp3 and Cbp6. We conclude that UQCC3 functions in complex III assembly and that the c.59T>A mutation has a causal role in complex III deficiency.

Control of mitochondrial integrity in ageing and disease.

Various molecular and cellular pathways are active in eukaryotes to control the quality and integrity of mitochondria. These pathways are involved in keeping a 'healthy' population of this essential organelle during the lifetime of the organism. Quality control (QC) systems counteract processes that lead to organellar dysfunction manifesting as degenerative diseases and ageing. We discuss disease- and ageing-related pathways involved in mitochondrial QC: mtDNA repair and reorganization, regeneration of oxidized amino acids, refolding and degradation of severely damaged proteins, degradation of whole mitochondria by mitophagy and finally programmed cell death. The control of the integrity of mtDNA and regulation of its expression is essential to remodel single proteins as well as mitochondrial complexes that determine mitochondrial functions. The redundancy of components, such as proteases, and the hierarchies of the QC raise questions about crosstalk between systems and their precise regulation. The understanding of the underlying mechanisms on the genomic, proteomic, organellar and cellular levels holds the key for the development of interventions for mitochondrial dysfunctions, degenerative processes, ageing and age-related diseases resulting from impairments of mitochondria.

Loss, replacement and gain of proteins at the origin of the mitochondria.

We review what has been inferred about the changes at the level of the proteome that accompanied the evolution of the mitochondrion from an alphaproteobacterium. We regard these changes from an alphaproteobacterial perspective: which proteins were lost during mitochondrial evolution? And, of the proteins that were lost, which ones have been replaced by other, non-orthologous proteins with a similar function? Combining literature-supported replacements with quantitative analyses of mitochondrial proteomics data we infer that most of the loss and replacements that separate current day mitochondria in mammals from alphaproteobacteria took place before the radiation of the eukaryotes. Recent analyses show that also the acquisition of new proteins to the large protein complexes of the oxidative phosphorylation and the mitochondrial ribosome occurred mainly before the divergence of the eukaryotes. These results indicate a significant number of pivotal evolutionary events between the acquisition of the endosymbiont and the radiation of the eukaryotes and therewith support an early acquisition of mitochondria in eukaryotic evolution. Technically, advancements in the reconstruction of the evolutionary trajectories of loss, replacement and gain of mitochondrial proteins depend on using profile-based homology detection methods for sequence analysis. We highlight the mitochondrial Holliday junction resolvase endonuclease, for which such methods have detected new "family members" and in which function differentiation is accompanied by the loss of catalytic residues for the original enzymatic function and the gain of a protein domain for the new function. This article is part of a Special Issue entitled: The evolutionary aspects of bioenergetic systems.

BOLA1 is an aerobic protein that prevents mitochondrial morphology changes induced by glutathione depletion.

AIMS: The BolA protein family is widespread among eukaryotes and bacteria. In Escherichia coli, BolA causes a spherical cell shape and is overexpressed during oxidative stress. Here we aim to elucidate the possible role of its human homolog BOLA1 in mitochondrial morphology and thiol redox potential regulation. RESULTS: We show that BOLA1 is a mitochondrial protein that counterbalances the effect of L-buthionine-(S,R)-sulfoximine (BSO)-induced glutathione (GSH) depletion on the mitochondrial thiol redox potential. Furthermore, overexpression of BOLA1 nullifies the effect of BSO and S-nitrosocysteine on mitochondrial morphology. Conversely, knockdown of the BOLA1 gene increases the oxidation of mitochondrial thiol groups. Supporting a role of BOLA1 in controlling the mitochondrial thiol redox potential is that BOLA1 orthologs only occur in aerobic eukaryotes. A measured interaction of BOLA1 with the mitochondrial monothiol glutaredoxin GLRX5 provides hints for potential mechanisms behind BOLA1's effect on mitochondrial redox potential. Nevertheless, we have no direct evidence for a role of GLRX5 in BOLA1's function. INNOVATION: We implicate a new protein, BOLA1, in the regulation of the mitochondrial thiol redox potential. CONCLUSION: BOLA1 is an aerobic, mitochondrial protein that prevents mitochondrial morphology aberrations induced by GSH depletion and reduces the associated oxidative shift of the mitochondrial thiol redox potential.

A mutation in the FAM36A gene, the human ortholog of COX20, impairs cytochrome c oxidase assembly and is associated with ataxia and muscle hypotonia.

The mitochondrial respiratory chain complex IV (cytochrome c oxidase) is a multi-subunit enzyme that transfers electrons from cytochrome c to molecular oxygen, yielding water. Its biogenesis requires concerted expression of mitochondria- and nuclear-encoded subunits and assembly factors. In this report, we describe a homozygous missense mutation in FAM36A from a patient who displays ataxia and muscle hypotonia. The FAM36A gene is a remote, putative ortholog of the fungal complex IV assembly factor COX20. Messenger RNA (mRNA) and protein co-expression analyses support the involvement of FAM36A in complex IV function in mammals. The c.154A>C mutation in the FAM36A gene, a mutation that is absent in sequenced exomes, leads to a reduced activity and lower levels of complex IV and its protein subunits. The FAM36A protein is nearly absent in patient's fibroblasts. Cells affected by the mutation accumulate subassemblies of complex IV that contain COX1 but are almost devoid of COX2 protein. We observe co-purification of FAM36A and COX2 proteins, supporting that the FAM36A defect hampers the early step of complex IV assembly at the incorporation of the COX2 subunit. Lentiviral complementation of patient's fibroblasts with wild-type FAM36A increases the complex IV activity as well as the amount of holocomplex IV and of individual subunits. These results establish the function of the human gene FAM36A/COX20 in complex IV assembly and support a causal role of the gene in complex IV deficiency.