CRISPR/Cas9-based functional genomics strategy to decipher the pathogenicity of genetic variants in inherited metabolic disorders.

The determination of the functional impact of variants of uncertain significance (VUS) is one of the major bottlenecks in the diagnostic workflow of inherited genetic diseases. To face this problem, we set up a CRISPR/Cas9-based strategy for knock-in cellular model generation, focusing on inherited metabolic disorders (IMDs). We selected variants in seven IMD-associated genes, including seven reported disease-causing variants and four benign/likely benign variants. Overall, 11 knock-in cell models were generated via homology-directed repair in HAP1 haploid cells using CRISPR/Cas9. The functional impact of the variants was determined by analyzing the characteristic biochemical alterations of each disorder. Functional studies performed in knock-in cell models showed that our approach accurately distinguished the functional effect of pathogenic from non-pathogenic variants in a reliable manner in a wide range of IMDs. Our study provides a generic approach to assess the functional impact of genetic variants to improve IMD diagnosis and this tool could emerge as a promising alternative to invasive tests, such as muscular or skin biopsies. Although the study has been performed only in IMDs, this strategy is generic and could be applied to other genetic disorders.

Functional Evidence of CCDC186 as a New Disease-Associated Gene with Endocrine and Central Nervous System Alterations.

CCDC186 protein is involved in the maturation of dense-core vesicles (DCVs) in the trans-Golgi network in neurons and endocrine cells. Mutations in genes involved in DCV regulation, other than CCDC186, have been described in patients with neurodevelopmental disorders. To date, only one patient, within a large sequencing study of 1000 cases, and a single case report with variants in CCDC186, had previously been described. However, no functional studies in any of these two cases had been performed. We identified three patients from two gypsy families, unrelated to each other, with mutations in the CCDC186 gene. Clinically, all patients presented with seizures, frontotemporal atrophy, hypomyelination, recurrent infections, and endocrine disturbances such as severe non-ketotic hypoglycemia. Low levels of cortisol, insulin, or growth hormone could only be verified in one patient. All of them had a neonatal onset and died between 7 months and 4 years of age. Whole exome sequencing identified a homozygous variant in the CCDC186 gene (c.2215C>T, p.Arg739Ter) in the index patients of both families. Protein expression studies demonstrated that CCDC186 was almost undetectable in fibroblasts and muscle tissue. These observations correlated with the transcriptomic analysis performed in fibroblasts in one of the patients, which showed a significant reduction of CCDC186 mRNA levels. Our study provides functional evidence that mutations in this gene have a pathogenic effect on the protein and reinforces CCDC186 as a new disease-associated gene. In addition, mutations in CCDC186 could explain the combined endocrine and neurologic alterations detected in our patients.

Leigh syndrome is the main clinical characteristic of PTCD3 deficiency.

Mitochondrial translation defects are a continuously growing group of disorders showing a large variety of clinical symptoms including a wide range of neurological abnormalities. To date, mutations in PTCD3, encoding a component of the mitochondrial ribosome, have only been reported in a single individual with clinical evidence of Leigh syndrome. Here, we describe three additional PTCD3 individuals from two unrelated families, broadening the genetic and phenotypic spectrum of this disorder, and provide definitive evidence that PTCD3 deficiency is associated with Leigh syndrome. The patients presented in the first months of life with psychomotor delay, respiratory insufficiency and feeding difficulties. The neurologic phenotype included dystonia, optic atrophy, nystagmus and tonic-clonic seizures. Brain MRI showed optic nerve atrophy and thalamic changes, consistent with Leigh syndrome. WES and RNA-seq identified compound heterozygous variants in PTCD3 in both families: c.[1453-1G>C];[1918C>G] and c.[710del];[902C>T]. The functional consequences of the identified variants were determined by a comprehensive characterization of the mitochondrial function. PTCD3 protein levels were significantly reduced in patient fibroblasts and, consistent with a mitochondrial translation defect, a severe reduction in the steady state levels of complexes I and IV subunits was detected. Accordingly, the activity of these complexes was also low, and high-resolution respirometry showed a significant decrease in the mitochondrial respiratory capacity. Functional complementation studies demonstrated the pathogenic effect of the identified variants since the expression of wild-type PTCD3 in immortalized fibroblasts restored the steady-state levels of complexes I and IV subunits as well as the mitochondrial respiratory capacity. Additionally, minigene assays demonstrated that three of the identified variants were pathogenic by altering PTCD3 mRNA processing. The fourth variant was a frameshift leading to a truncated protein. In summary, we provide evidence of PTCD3 involvement in human disease confirming that PTCD3 deficiency is definitively associated with Leigh syndrome.

Diagnostic Odyssey in an Adult Patient with Ophthalmologic Abnormalities and Hearing Loss: Contribution of RNA-Seq to the Diagnosis of a PEX1 Deficiency.

Peroxisomal biogenesis disorders (PBDs) are a heterogeneous group of genetic diseases. Multiple peroxisomal pathways are impaired, and very long chain fatty acids (VLCFA) are the first line biomarkers for the diagnosis. The clinical presentation of PBDs may range from severe, lethal multisystemic disorders to milder, late-onset disease. The vast majority of PBDs belong to Zellweger Spectrum Disordes (ZSDs) and represents a continuum of overlapping clinical symptoms, with Zellweger syndrome being the most severe and Heimler syndrome the less severe disease. Mild clinical conditions frequently present normal or slight biochemical alterations, making the diagnosis of these patients challenging. In the present study we used a combined WES and RNA-seq strategy to diagnose a patient presenting with retinal dystrophy as the main clinical symptom. Results showed the patient was compound heterozygous for mutations in PEX1. VLCFA were normal, but retrospective analysis of lysosphosphatidylcholines (LPC) containing C22:0-C26:0 species was altered. This simple test could avoid the diagnostic odyssey of patients with mild phenotype, such as the individual described here, who was diagnosed very late in adult life. We provide functional data in cell line models that may explain the mild phenotype of the patient by demonstrating the hypomorphic nature of a deep intronic variant altering PEX1 mRNA processing.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics.

BACKGROUND: Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics. METHODS: We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage. RESULTS: We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions. CONCLUSION: Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.

Over-Mutated Mitochondrial, Lysosomal and TFEB-Regulated Genes in Parkinson's Disease.

The association between Parkinson's disease (PD) and mutations in genes involved in lysosomal and mitochondrial function has been previously reported. However, little is known about the involvement of other genes or cellular mechanisms. We aim to identify novel genetic associations to better understand the pathogenesis of PD. We performed WES in a cohort of 32 PD patients and 30 age-matched controls. We searched for rare variants in 1667 genes: PD-associated, related to lysosomal function and mitochondrial function and TFEB-regulated. When comparing the PD patient cohort with that of age matched controls, a statistically significant burden of rare variants in the previous group of genes were identified. In addition, the Z-score calculation, using the European population database (GnomAD), showed an over-representation of particular variants in 36 genes. Interestingly, 11 of these genes are implicated in mitochondrial function and 18 are TFEB-regulated genes. Our results suggest, for the first time, an involvement of TFEB-regulated genes in the genetic susceptibility to PD. This is remarkable as TFEB factor has been reported to be sequestered inside Lewy bodies, pointing to a role of TFEB in the pathogenesis of PD. Our data also reinforce the involvement of lysosomal and mitochondrial mechanisms in PD.

Clinico-radiological features, molecular spectrum, and identification of prognostic factors in developmental and epileptic encephalopathy due to inosine triphosphate pyrophosphatase (ITPase) deficiency.

Developmental and epileptic encephalopathy 35 (DEE 35) is a severe neurological condition caused by biallelic variants in ITPA, encoding inosine triphosphate pyrophosphatase, an essential enzyme in purine metabolism. We delineate the genotypic and phenotypic spectrum of DEE 35, analyzing possible predictors for adverse clinical outcomes. We investigated a cohort of 28 new patients and reviewed previously described cases, providing a comprehensive characterization of 40 subjects. Exome sequencing was performed to identify underlying ITPA pathogenic variants. Brain MRI (magnetic resonance imaging) scans were systematically analyzed to delineate the neuroradiological spectrum. Survival curves according to the Kaplan-Meier method and log-rank test were used to investigate outcome predictors in different subgroups of patients. We identified 18 distinct ITPA pathogenic variants, including 14 novel variants, and two deletions. All subjects showed profound developmental delay, microcephaly, and refractory epilepsy followed by neurodevelopmental regression. Brain MRI revision revealed a recurrent pattern of delayed myelination and restricted diffusion of early myelinating structures. Congenital microcephaly and cardiac involvement were statistically significant novel clinical predictors of adverse outcomes. We refined the molecular, clinical, and neuroradiological characterization of ITPase deficiency, and identified new clinical predictors which may have a potentially important impact on diagnosis, counseling, and follow-up of affected individuals.

Implementation of second-tier tests in newborn screening for the detection of vitamin B12 related acquired and genetic disorders: results on 258,637 newborns.

BACKGROUND: Alteration of vitamin B12 metabolism can be genetic or acquired, and can result in anemia, failure to thrive, developmental regression and even irreversible neurologic damage. Therefore, early diagnosis and intervention is critical. Most of the neonatal cases with acquired vitamin B12 deficiency have been detected by clinical symptoms and only few of them trough NBS programs. We aim to assess the usefulness of the second-tier test: methylmalonic acid (MMA), methylcitric acid (MCA) and homocysteine (Hcys) in our newborn screening program and explore the implications on the detection of cobalamin (vitamin B12) related disorders, both genetic and acquired conditions. METHODS: A screening strategy using the usual primary markers followed by the analysis of MMA, MCA and Hcys as second tier-test in the first dried blood spot (DBS) was developed and evaluated. RESULTS: During the period 2015-2018 a total of 258,637 newborns were screened resulting in 130 newborns with acquired vitamin B12 deficiency (incidence 1:1989), 19 with genetic disorders (incidence 1:13,613) and 13 were false positive. No false negatives were notified. Concerning the second-tier test, the percentage of cases with MMA above the cut-off levels, both for genetic and acquired conditions was very similar (58% and 60%, respectively). Interestingly, the percentage of cases with increased levels of Hcys was higher in acquired conditions than in genetic disorders (87% and 47%, respectively). In contrast, MCA was high only in 5% of the acquired conditions versus in 53% of the genetic disorders, and it was always very high in all patients with propionic acidemia. CONCLUSIONS: When screening for methylmalonic acidemia and homocystinuria, differential diagnosis with acquired vitamin B12 deficiency should be done. The results of our strategy support the inclusion of this acquired condition in the NBS programs, as it is easily detectable and allows the adoption of corrective measures to avoid the consequences of its deficiency.

Delineating the neurological phenotype in children with defects in the ECHS1 or HIBCH gene.

The neurological phenotype of 3-hydroxyisobutyryl-CoA hydrolase (HIBCH) and short-chain enoyl-CoA hydratase (SCEH) defects is expanding and natural history studies are necessary to improve clinical management. From 42 patients with Leigh syndrome studied by massive parallel sequencing, we identified five patients with SCEH and HIBCH deficiency. Fourteen additional patients were recruited through collaborations with other centres. In total, we analysed the neurological features and mutation spectrum in 19 new SCEH/HIBCH patients. For natural history studies and phenotype to genotype associations we also included 70 previously reported patients. The 19 newly identified cases presented with Leigh syndrome (SCEH, n = 11; HIBCH, n = 6) and paroxysmal dystonia (SCEH, n = 2). Basal ganglia lesions (18 patients) were associated with small cysts in the putamen/pallidum in half of the cases, a characteristic hallmark for diagnosis. Eighteen pathogenic variants were identified, 11 were novel. Among all 89 cases, we observed a longer survival in HIBCH compared to SCEH patients, and in HIBCH patients carrying homozygous mutations on the protein surface compared to those with variants inside/near the catalytic region. The SCEH p.(Ala173Val) change was associated with a milder form of paroxysmal dystonia triggered by increased energy demands. In a child harbouring SCEH p.(Ala173Val) and the novel p.(Leu123Phe) change, an 83.6% reduction of the protein was observed in fibroblasts. The SCEH and HIBCH defects in the catabolic valine pathway were a frequent cause of Leigh syndrome in our cohort. We identified phenotype and genotype associations that may help predict outcome and improve clinical management.

International consensus guidelines for phosphoglucomutase 1 deficiency (PGM1-CDG): Diagnosis, follow-up, and management.

Phosphoglucomutase 1 (PGM1) deficiency is a rare genetic disorder that affects glycogen metabolism, glycolysis, and protein glycosylation. Previously known as GSD XIV, it was recently reclassified as a congenital disorder of glycosylation, PGM1-CDG. PGM1-CDG usually manifests as a multisystem disease. Most patients present as infants with cleft palate, liver function abnormalities and hypoglycemia, but some patients present in adulthood with isolated muscle involvement. Some patients develop life-threatening cardiomyopathy. Unlike most other CDG, PGM1-CDG has an effective treatment option, d-galactose, which has been shown to improve many of the patients' symptoms. Therefore, early diagnosis and initiation of treatment for PGM1-CDG patients are crucial decisions. In this article, our group of international experts suggests diagnostic, follow-up, and management guidelines for PGM1-CDG. These guidelines are based on the best available evidence-based data and experts' opinions aiming to provide a practical resource for health care providers to facilitate successful diagnosis and optimal management of PGM1-CDG patients.

Leigh syndrome associated with TRMU gene mutations.

tRNA 5-methylaminomethyl-2-thiouridylate methyltransferase (TRMU) deficiency causes an early onset potentially reversible acute liver failure, so far reported in less than 30 patients. We describe two new unrelated patients with an acute liver failure and a neuroimaging compatible with Leigh syndrome (LS) due to TRMU deficiency, a combination not previously reported. Our report enlarges the phenotypical spectrum of TRMU disease.

Biallelic mutations in NDUFA8 cause complex I deficiency in two siblings with favorable clinical evolution.

Isolated complex I (CI) deficiency is the most common cause of oxidative phosphorylation (OXPHOS) dysfunction. Whole-exome sequencing identified biallelic mutations in NDUFA8 (c.[293G > T]; [293G > T], encoding for an accessory subunit of CI, in two siblings with a favorable clinical evolution. The individuals reported here are practically asymptomatic, with the exception of slight failure to thrive and some language difficulties at the age of 6 and 9 years, respectively. These observations are remarkable since the vast majority of patients with CI deficiency, including the only NDUFA8 patient reported so far, showed an extremely poor clinical outcome. Western blot studies demonstrated that NDUFA8 protein was strongly reduced in the patients' fibroblasts and muscle extracts. In addition, there was a marked and specific decrease in the steady-state levels of CI subunits. BN-PAGE demonstrated an isolated defect in the assembly and the activity of CI with impaired supercomplexes formation and abnormal accumulation of CI subassemblies. Confocal microscopy analysis in fibroblasts showed rounder mitochondria and diminished branching degree of the mitochondrial network. Functional complementation studies demonstrated disease-causality for the identified mutation as lentiviral transduction with wild-type NDUFA8 cDNA restored the steady-state levels of CI subunits and completely recovered the deficient enzymatic activity in immortalized mutant fibroblasts. In summary, we provide additional evidence of the involvement of NDUFA8 as a mitochondrial disease-causing gene associated with altered mitochondrial morphology, CI deficiency, impaired supercomplexes formation, and very mild progression of the disease.

Complex I deficiency, due to NDUFAF4 mutations, causes severe mitochondrial dysfunction and is associated to early death and dysmorphia.

Pathogenic mutations in NDUFAF4 have been reported in very few cases. Here we present new data to further delineate the phenotypic spectrum of NDUFAF4 deficiency. We describe two siblings presenting with facial dysmorphia and lactic acidosis in the neonatal period. Later on, they developed fatal early encephalopathy with apneic episodes, irritability, central hypoventilation, liver involvement and hyperammonemia. Abnormality of the cerebral white matter was demonstrated in one case, and cardiomyopathy in the other. Urine organic acid profile showed an increased excretion of lactate, Krebs cycle metabolites and 3-methylglutaconate. Whole-exome sequencing identified a novel homozygous nonsense mutation in NDUFAF4 (c.478G > T; p.Glu160Ter), encoding a mitochondrial complex I assembly factor. The disruptive effect of the mutation was corroborated by the absence of NDUFAF4 expression in patient fibroblasts. OXPHOS assembly studies demonstrated almost undetectable levels of fully assembled complex I and complex I-containing supercomplexes and an abnormal accumulation of SCIII2IV1 supercomplexes. Morphologically, fibroblasts showed rounder mitochondria and a diminished degree of branching of the mitochondrial network. Cellular respiratory capacity in fibroblasts was also markedly reduced. In sum, we provide insights into the physiopathological mechanisms underlying NDUFAF4 deficiency and expand the knowledge about the clinical and biochemical spectrum of this disorder.

Physiopathological Bases of the Disease Caused by HACE1 Mutations: Alterations in Autophagy, Mitophagy and Oxidative Stress Response.

Recessive HACE1 mutations are associated with a severe neurodevelopmental disorder (OMIM: 616756). However, the physiopathologycal bases of the disease are yet to be completely clarified. Whole-exome sequencing identified homozygous HACE1 mutations (c.240C>A, p.Cys80Ter) in a patient with brain atrophy, psychomotor retardation and 3-methylglutaconic aciduria, a biomarker of mitochondrial dysfunction. To elucidate the pathomechanisms underlying HACE1 deficiency, a comprehensive molecular analysis was performed in patient fibroblasts. Western Blot demonstrated the deleterious effect of the mutation, as the complete absence of HACE1 protein was observed. Immunofluorescence studies showed an increased number of LC3 puncta together with the normal initiation of the autophagic cascade, indicating a reduction in the autophagic flux. Oxidative stress response was also impaired in HACE1 fibroblasts, as shown by the reduced NQO1 and Hmox1 mRNA levels observed in H2O2-treated cells. High levels of lipid peroxidation, consistent with accumulated oxidative damage, were also detected. Although the patient phenotype could resemble a mitochondrial defect, the analysis of the mitochondrial function showed no major abnormalities. However, an important increase in mitochondrial oxidative stress markers and a strong reduction in the mitophagic flux were observed, suggesting that the recycling of damaged mitochondria might be targeted in HACE1 cells. In summary, we demonstrate for the first time that the impairment of autophagy, mitophagy and oxidative damage response might be involved in the pathogenesis of HACE1 deficiency.

Clinical presentation and proteomic signature of patients with TANGO2 mutations.

Transport And Golgi Organization protein 2 (TANGO2) deficiency has recently been identified as a rare metabolic disorder with a distinct clinical and biochemical phenotype of recurrent metabolic crises, hypoglycemia, lactic acidosis, rhabdomyolysis, arrhythmias, and encephalopathy with cognitive decline. We report nine subjects from seven independent families, and we studied muscle histology, respiratory chain enzyme activities in skeletal muscle and proteomic signature of fibroblasts. All nine subjects carried autosomal recessive TANGO2 mutations. Two carried the reported deletion of exons 3 to 9, one homozygous, one heterozygous with a 22q11.21 microdeletion inherited in trans. The other subjects carried three novel homozygous (c.262C>T/p.Arg88*; c.220A>C/p.Thr74Pro; c.380+1G>A), and two further novel heterozygous (c.6_9del/p.Phe6del); c.11-13delTCT/p.Phe5del mutations. Immunoblot analysis detected a significant decrease of TANGO2 protein. Muscle histology showed mild variation of fiber diameter, no ragged-red/cytochrome c oxidase-negative fibers and a defect of multiple respiratory chain enzymes and coenzyme Q10 (CoQ10 ) in two cases, suggesting a possible secondary defect of oxidative phosphorylation. Proteomic analysis in fibroblasts revealed significant changes in components of the mitochondrial fatty acid oxidation, plasma membrane, endoplasmic reticulum-Golgi network and secretory pathways. Clinical presentation of TANGO2 mutations is homogeneous and clinically recognizable. The hemizygous mutations in two patients suggest that some mutations leading to allele loss are difficult to detect. A combined defect of the respiratory chain enzymes and CoQ10 with altered levels of several membrane proteins provides molecular insights into the underlying pathophysiology and may guide rational new therapeutic interventions.

Mutations in TIMM50 cause severe mitochondrial dysfunction by targeting key aspects of mitochondrial physiology.

3-Methylglutaconic aciduria (3-MGA-uria) syndromes comprise a heterogeneous group of diseases associated with mitochondrial membrane defects. Whole-exome sequencing identified compound heterozygous mutations in TIMM50 (c.[341 G>A];[805 G>A]) in a boy with West syndrome, optic atrophy, neutropenia, cardiomyopathy, Leigh syndrome, and persistent 3-MGA-uria. A comprehensive analysis of the mitochondrial function was performed in fibroblasts of the patient to elucidate the molecular basis of the disease. TIMM50 protein was severely reduced in the patient fibroblasts, regardless of the normal mRNA levels, suggesting that the mutated residues might be important for TIMM50 protein stability. Severe morphological defects and ultrastructural abnormalities with aberrant mitochondrial cristae organization in muscle and fibroblasts were found. The levels of fully assembled OXPHOS complexes and supercomplexes were strongly reduced in fibroblasts from this patient. High-resolution respirometry demonstrated a significant reduction of the maximum respiratory capacity. A TIMM50-deficient HEK293T cell line that we generated using CRISPR/Cas9 mimicked the respiratory defect observed in the patient fibroblasts; notably, this defect was rescued by transfection with a plasmid encoding the TIMM50 wild-type protein. In summary, we demonstrated that TIMM50 deficiency causes a severe mitochondrial dysfunction by targeting key aspects of mitochondrial physiology, such as the maintenance of proper mitochondrial morphology, OXPHOS assembly, and mitochondrial respiratory capacity.

Muscle Involvement in a Large Cohort of Pediatric Patients with Genetic Diagnosis of Mitochondrial Disease.

Mitochondrial diseases (MD) are a group of genetic and acquired disorders which present significant diagnostic challenges. Here we report the disease characteristics of a large cohort of pediatric MD patients (n = 95) with a definitive genetic diagnosis, giving special emphasis on clinical muscle involvement, biochemical and histopathological features. Of the whole cohort, 51 patients harbored mutations in nuclear DNA (nDNA) genes and 44 patients had mutations in mitochondrial DNA (mtDNA) genes. The nDNA patients were more likely to have a reduction in muscle fiber succinate dehydrogenase (SDH) stains and in SDH-positive blood vessels, while a higher frequency of mtDNA patients had ragged red (RRF) and blue fibers. The presence of positive histopathological features was associated with ophthalmoplegia, myopathic facies, weakness and exercise intolerance. In 17 patients younger than two years of age, RRF and blue fibers were observed only in one case, six cases presented cytochrome c oxidase (COX) reduction/COX-fibers, SDH reduction was observed in five and all except one presented SDH-positive blood vessels. In conclusion, muscle involvement was a frequent finding in our series of MD patients, especially in those harboring mutations in mtDNA genes.