Genome sequencing of C. elegans balancer strains reveals previously unappreciated complex genomic rearrangements.

Genetic balancers in Caenorhabditis elegans are complex variants that allow lethal or sterile mutations to be stably maintained in a heterozygous state by suppressing crossover events. Balancers constitute an invaluable tool in the C. elegans scientific community and have been widely used for decades. The first/traditional balancers were created by applying X-rays, UV, or gamma radiation on C. elegans strains, generating random genomic rearrangements. Their structures have been mostly explored with low-resolution genetic techniques (e.g., fluorescence in situ hybridization or PCR), before genomic mapping and molecular characterization through sequencing became feasible. As a result, the precise nature of most chromosomal rearrangements remains unknown, whereas, more recently, balancers have been engineered using the CRISPR-Cas9 technique for which the structure of the chromosomal rearrangement has been predesigned. Using short-read whole-genome sequencing (srWGS) and tailored bioinformatic analyses, we previously interpreted the structure of four chromosomal balancers randomly created by mutagenesis processes. Here, we have extended our analyses to five CRISPR-Cas9 balancers and 17 additional traditional balancing rearrangements. We detected and experimentally validated their breakpoints and have interpreted the balancer structures. Many of the balancers were found to be more intricate than previously described, being composed of complex genomic rearrangements (CGRs) such as chromoanagenesis-like events. Furthermore, srWGS revealed additional structural variants and CGRs not known to be part of the balancer genomes. Altogether, our study provides a comprehensive resource of complex genomic variations in C. elegans and highlights the power of srWGS to study the complexity of genomes by applying tailored analyses.

GPAD: a natural language processing-based application to extract the gene-disease association discovery information from OMIM.

BACKGROUND: Thousands of genes have been associated with different Mendelian conditions. One of the valuable sources to track these gene-disease associations (GDAs) is the Online Mendelian Inheritance in Man (OMIM) database. However, most of the information in OMIM is textual, and heterogeneous (e.g. summarized by different experts), which complicates automated reading and understanding of the data. Here, we used Natural Language Processing (NLP) to make a tool (Gene-Phenotype Association Discovery (GPAD)) that could syntactically process OMIM text and extract the data of interest. RESULTS: GPAD applies a series of language-based techniques to the text obtained from OMIM API to extract GDA discovery-related information. GPAD can inform when a particular gene was associated with a specific phenotype, as well as the type of validation-whether through model organisms or cohort-based patient-matching approaches-for such an association. GPAD extracted data was validated with published reports and was compared with large language model. Utilizing GPAD's extracted data, we analysed trends in GDA discoveries, noting a significant increase in their rate after the introduction of exome sequencing, rising from an average of about 150-250 discoveries each year. Contrary to hopes of resolving most GDAs for Mendelian disorders by now, our data indicate a substantial decline in discovery rates over the past five years (2017-2022). This decline appears to be linked to the increasing necessity for larger cohorts to substantiate GDAs. The rising use of zebrafish and Drosophila as model organisms in providing evidential support for GDAs is also observed. CONCLUSIONS: GPAD's real-time analyzing capacity offers an up-to-date view of GDA discovery and could help in planning and managing the research strategies. In future, this solution can be extended or modified to capture other information in OMIM and scientific literature.

CIAO1 and MMS19 deficiency: A lethal neurodegenerative phenotype caused by cytosolic Fe-S cluster protein assembly disorders.

PURPOSE: The functionality of many cellular proteins depends on cofactors; yet, they have only been implicated in a minority of Mendelian diseases. Here, we describe the first 2 inherited disorders of the cytosolic iron-sulfur protein assembly system. METHODS: Genetic testing via genome sequencing was applied to identify the underlying disease cause in 3 patients with microcephaly, congenital brain malformations, progressive developmental and neurologic impairments, recurrent infections, and a fatal outcome. Studies in patient-derived skin fibroblasts and zebrafish models were performed to investigate the biochemical and cellular consequences. RESULTS: Metabolic analysis showed elevated uracil and thymine levels in body fluids but no pathogenic variants in DPYD, encoding dihydropyrimidine dehydrogenase. Genome sequencing identified compound heterozygosity in 2 patients for missense variants in CIAO1, encoding cytosolic iron-sulfur assembly component 1, and homozygosity for an in-frame 3-nucleotide deletion in MMS19, encoding the MMS19 homolog, cytosolic iron-sulfur assembly component, in the third patient. Profound alterations in the proteome, metabolome, and lipidome were observed in patient-derived fibroblasts. We confirmed the detrimental effect of deficiencies in CIAO1 and MMS19 in zebrafish models. CONCLUSION: A general failure of cytosolic and nuclear iron-sulfur protein maturation caused pleiotropic effects. The critical function of the cytosolic iron-sulfur protein assembly machinery for antiviral host defense may well explain the recurrent severe infections occurring in our patients.

Identification of a mosaic MTOR variant in purified neuronal DNA in a patient with focal cortical dysplasia using a novel depth electrode harvesting technique.

OBJECTIVE: Recent studies have identified brain somatic variants as a cause of focal epilepsy. These studies relied on resected tissue from epilepsy surgery, which is not available in most patients. The use of trace tissue adherent to depth electrodes used for stereo electroencephalography (EEG) has been proposed as an alternative but is hampered by the low cell quality and contamination by nonbrain cells. Here, we use our improved depth electrode harvesting technique that purifies neuronal nuclei to achieve molecular diagnosis in a patient with focal cortical dysplasia (FCD). METHODS: Depth electrode tips were collected, pooled by brain region and seizure onset zone, and nuclei were isolated and sorted using fluorescence-activated nuclei sorting (FANS). Somatic DNA was amplified from neuronal and astrocyte nuclei using primary template amplification followed by exome sequencing of neuronal DNA from the affected pool, unaffected pool, and saliva. The identified variant was validated using droplet digital polymerase chain reaction (PCR). RESULTS: An 11-year-old male with drug-resistant genetic-structural epilepsy due to left anterior insula FCD had seizures from age 3 years. Stereo EEG confirmed seizure onset in the left anterior insula. The two anterior insula electrodes were combined as the affected pool and three frontal electrodes as the unaffected pool. FANS isolated 140 neuronal nuclei from the affected and 245 neuronal nuclei from the unaffected pool. A novel somatic missense MTOR variant (p.Leu489Met, CADD score 23.7) was identified in the affected neuronal sample. Droplet digital PCR confirmed a mosaic gradient (variant allele frequency = .78% in affected neuronal sample; variant was absent in all other samples). SIGNIFICANCE: Our findings confirm that harvesting neuronal DNA from depth electrodes followed by molecular analysis to identify brain somatic variants is feasible. Our novel method represents a significant improvement compared to the previous method by focusing the analysis on high-quality cells of the cell type of interest.

Bi-allelic GOT2 Mutations Cause a Treatable Malate-Aspartate Shuttle-Related Encephalopathy.

Early-infantile encephalopathies with epilepsy are devastating conditions mandating an accurate diagnosis to guide proper management. Whole-exome sequencing was used to investigate the disease etiology in four children from independent families with intellectual disability and epilepsy, revealing bi-allelic GOT2 mutations. In-depth metabolic studies in individual 1 showed low plasma serine, hypercitrullinemia, hyperlactatemia, and hyperammonemia. The epilepsy was serine and pyridoxine responsive. Functional consequences of observed mutations were tested by measuring enzyme activity and by cell and animal models. Zebrafish and mouse models were used to validate brain developmental and functional defects and to test therapeutic strategies. GOT2 encodes the mitochondrial glutamate oxaloacetate transaminase. GOT2 enzyme activity was deficient in fibroblasts with bi-allelic mutations. GOT2, a member of the malate-aspartate shuttle, plays an essential role in the intracellular NAD(H) redox balance. De novo serine biosynthesis was impaired in fibroblasts with GOT2 mutations and GOT2-knockout HEK293 cells. Correcting the highly oxidized cytosolic NAD-redox state by pyruvate supplementation restored serine biosynthesis in GOT2-deficient cells. Knockdown of got2a in zebrafish resulted in a brain developmental defect associated with seizure-like electroencephalography spikes, which could be rescued by supplying pyridoxine in embryo water. Both pyridoxine and serine synergistically rescued embryonic developmental defects in zebrafish got2a morphants. The two treated individuals reacted favorably to their treatment. Our data provide a mechanistic basis for the biochemical abnormalities in GOT2 deficiency that may also hold for other MAS defects.

Glutaminase Deficiency Caused by Short Tandem Repeat Expansion in GLS.

We report an inborn error of metabolism caused by an expansion of a GCA-repeat tract in the 5' untranslated region of the gene encoding glutaminase (GLS) that was identified through detailed clinical and biochemical phenotyping, combined with whole-genome sequencing. The expansion was observed in three unrelated patients who presented with an early-onset delay in overall development, progressive ataxia, and elevated levels of glutamine. In addition to ataxia, one patient also showed cerebellar atrophy. The expansion was associated with a relative deficiency of GLS messenger RNA transcribed from the expanded allele, which probably resulted from repeat-mediated chromatin changes upstream of the GLS repeat. Our discovery underscores the importance of careful examination of regions of the genome that are typically excluded from or poorly captured by exome sequencing.

Driving mosaicism: somatic variants in reference population databases and effect on variant interpretation in rare genetic disease.

BACKGROUND: Genetic variation databases provide invaluable information on the presence and frequency of genetic variants in the 'untargeted' human population, aggregated with the primary goal to facilitate the interpretation of clinically important variants. The presence of somatic variants in such databases can affect variant assessment in undiagnosed rare disease (RD) patients. Previously, the impact of somatic mosaicism was only considered in relation to two Mendelian disease-associated genes. Here, we expand the analyses to identify additional mosaicism-prone genes in blood-derived reference population databases. RESULTS: To identify additional mosaicism-prone genes relevant to RDs, we focused on known/previously established ClinVar pathogenic and likely pathogenic single-nucleotide variants, residing in genes associated with early onset, severe autosomal dominant diseases. We asked whether any of these variants are present in a higher-than-expected frequency in the reference population databases and whether there is evidence of somatic origin (i.e., allelic imbalance) rather than germline heterozygosity (~ half of the reads supporting alternative allele). The mosaicism-prone genes identified were further categorized according to the processes they are involved in. Beyond the previously reported ASXL1 and DNMT3A, we identified 7 additional autosomal dominant RD-associated genes with known pathogenic single-nucleotide variants present in the reference population databases and good evidence of allelic imbalance: BRAF, CBL, FGFR3, IDH2, KRAS, PTPN11 and SETBP1. From this group of 9 genes, the majority (n = 7) was important for hematopoiesis. In addition, 4 of these genes were involved in cell proliferation. Further assessment of the known 156 hematopoietic genes led to identification of 48 genes (21 not yet associated with RDs) with at least some evidence of mosaicism detectable in reference population databases. CONCLUSIONS: These results stress the importance of considering genes involved in hematopoiesis and cell proliferation when interpreting the presence and frequency of genetic variants in blood-derived reference population databases, both public and private. This is especially important when considering new variants of uncertain significance in known hematopoietic/cell proliferation RD genes and future novel gene-disease associations involving this class of genes.

De novo stop-loss variants in CLDN11 cause hypomyelinating leukodystrophy.

Claudin-11, a tight junction protein, is indispensable in the formation of the radial component of myelin. Here, we report de novo stop-loss variants in the gene encoding claudin-11, CLDN11, in three unrelated individuals presenting with an early-onset spastic movement disorder, expressive speech disorder and eye abnormalities including hypermetropia. Brain MRI showed a myelin deficit with a discrepancy between T1-weighted and T2-weighted images and some progress in myelination especially involving the central and peripheral white matter. Exome sequencing identified heterozygous stop-loss variants c.622T>C, p.(*208Glnext*39) in two individuals and c.622T>G, p.(*208Gluext*39) in one individual, all occurring de novo. At the RNA level, the variant c.622T>C did not lead to a loss of expression in fibroblasts, indicating this transcript is not subject to nonsense-mediated decay and most likely translated into an extended protein. Extended claudin-11 is predicted to form an alpha helix not incorporated into the cytoplasmic membrane, possibly perturbing its interaction with intracellular proteins. Our observations suggest that stop-loss variants in CLDN11 expand the genetically heterogeneous spectrum of hypomyelinating leukodystrophies.

Whole genome sequencing facilitates intragenic variant interpretation following modifier screening in C. elegans.

BACKGROUND: Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds. RESULTS: Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes. CONCLUSIONS: Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds.

Secondary biogenic amine deficiencies: genetic etiology, therapeutic interventions, and clinical effects.

Monoamine neurotransmitter disorders present predominantly with neurologic features, including dystonic or dyskinetic cerebral palsy and movement disorders. Genetic conditions that lead to secondary defects in the synthesis, catabolism, transport, and metabolism of biogenic amines can lead to neurotransmitter abnormalities, which can present with similar features. Eleven patients with secondary neurotransmitter abnormalities were enrolled between 2011 and 2015. All patients underwent research-based whole exome and/or whole genome sequencing (WES/WGS). A trial of treatment with levodopa/carbidopa and 5-hydroxytryptophan was initiated. In six families with abnormal neurotransmitter profiles and neurological phenotypes, variants in known disease-causing genes (KCNJ6, SCN2A, CSTB in 2 siblings, NRNX1, KIF1A and PAK3) were identified, while one patient had a variant of uncertain significance in a candidate gene (DLG4) that may explain her phenotype. In 3 patients, no compelling candidate genes were identified. A trial of neurotransmitter replacement therapy led to improvement in motor and behavioral symptoms in all but two patients. The patient with KCNJ6 variant did not respond to L-dopa therapy, but rather experienced increased dyskinetic movements even at low dose of medication. The patient's symptoms harboring the NRNX1 deletion remained unaltered. This study demonstrates the utility of genome-wide sequencing in further understanding the etiology and pathophysiology of neurometabolic conditions, and the potential of secondary neurotransmitter deficiencies to serve as novel therapeutic targets. As there was a largely favorable response to therapy in our case series, a careful trial of neurotransmitter replacement therapy should be considered in patients with cerebrospinal fluid (CSF) monoamines below reference range.

Germline De Novo Mutations in ATP1A1 Cause Renal Hypomagnesemia, Refractory Seizures, and Intellectual Disability.

Over the last decades, a growing spectrum of monogenic disorders of human magnesium homeostasis has been clinically characterized, and genetic studies in affected individuals have identified important molecular components of cellular and epithelial magnesium transport. Here, we describe three infants who are from non-consanguineous families and who presented with a disease phenotype consisting of generalized seizures in infancy, severe hypomagnesemia, and renal magnesium wasting. Seizures persisted despite magnesium supplementation and were associated with significant intellectual disability. Whole-exome sequencing and conventional Sanger sequencing identified heterozygous de novo mutations in the catalytic Na+, K+-ATPase α1 subunit (ATP1A1). Functional characterization of mutant Na+, K+-ATPase α1 subunits in heterologous expression systems revealed not only a loss of Na+, K+-ATPase function but also abnormal cation permeabilities, which led to membrane depolarization and possibly aggravated the effect of the loss of physiological pump activity. These findings underline the indispensable role of the α1 isoform of the Na+, K+-ATPase for renal-tubular magnesium handling and cellular ion homeostasis, as well as maintenance of physiologic neuronal activity.

PLPHP deficiency: clinical, genetic, biochemical, and mechanistic insights.

Biallelic pathogenic variants in PLPBP (formerly called PROSC) have recently been shown to cause a novel form of vitamin B6-dependent epilepsy, the pathophysiological basis of which is poorly understood. When left untreated, the disease can progress to status epilepticus and death in infancy. Here we present 12 previously undescribed patients and six novel pathogenic variants in PLPBP. Suspected clinical diagnoses prior to identification of PLPBP variants included mitochondrial encephalopathy (two patients), folinic acid-responsive epilepsy (one patient) and a movement disorder compatible with AADC deficiency (one patient). The encoded protein, PLPHP is believed to be crucial for B6 homeostasis. We modelled the pathogenicity of the variants and developed a clinical severity scoring system. The most severe phenotypes were associated with variants leading to loss of function of PLPBP or significantly affecting protein stability/PLP-binding. To explore the pathophysiology of this disease further, we developed the first zebrafish model of PLPHP deficiency using CRISPR/Cas9. Our model recapitulates the disease, with plpbp-/- larvae showing behavioural, biochemical, and electrophysiological signs of seizure activity by 10 days post-fertilization and early death by 16 days post-fertilization. Treatment with pyridoxine significantly improved the epileptic phenotype and extended lifespan in plpbp-/- animals. Larvae had disruptions in amino acid metabolism as well as GABA and catecholamine biosynthesis, indicating impairment of PLP-dependent enzymatic activities. Using mass spectrometry, we observed significant B6 vitamer level changes in plpbp-/- zebrafish, patient fibroblasts and PLPHP-deficient HEK293 cells. Additional studies in human cells and yeast provide the first empirical evidence that PLPHP is localized in mitochondria and may play a role in mitochondrial metabolism. These models provide new insights into disease mechanisms and can serve as a platform for drug discovery.

SimPEL: Simulation-based power estimation for sequencing studies of low-prevalence conditions.

Power estimations are important for optimizing genotype-phenotype association study designs. However, existing frameworks are designed for common disorders, and thus ill-suited for the inherent challenges of studies for low-prevalence conditions such as rare diseases and infrequent adverse drug reactions. These challenges include small sample sizes and the need to leverage genetic annotation resources in association analyses for the purpose of ranking potential causal genes. We present SimPEL, a simulation-based program providing power estimations for the design of low-prevalence condition studies. SimPEL integrates the usage of gene annotation resources for association analyses. Customizable parameters, including the penetrance of the putative causal allele and the employed pathogenic scoring system, allow SimPEL to realistically model a large range of study designs. To demonstrate the effects of various parameters on power, we estimated the power of several simulated designs using SimPEL and captured power trends in agreement with observations from current literature on low-frequency condition studies. SimPEL, as a tool, provides researchers studying low-frequency conditions with an intuitive and highly flexible avenue for statistical power estimation. The platform-independent "batteries included" executable and default input files are available at https://github.com/precisionomics/SimPEL.

Biallelic Mutations in UNC80 Cause Persistent Hypotonia, Encephalopathy, Growth Retardation, and Severe Intellectual Disability.

Ion channel proteins are required for both the establishment of resting membrane potentials and the generation of action potentials. Hundreds of mutations in genes encoding voltage-gated ion channels responsible for action potential generation have been found to cause severe neurological diseases. In contrast, the roles of voltage-independent "leak" channels, important for the establishment and maintenance of resting membrane potentials upon which action potentials are generated, are not well established in human disease. UNC80 is a large component of the NALCN sodium-leak channel complex that regulates the basal excitability of the nervous system. Loss-of-function mutations of NALCN cause infantile hypotonia with psychomotor retardation and characteristic facies (IHPRF). We report four individuals from three unrelated families who have homozygous missense or compound heterozygous truncating mutations in UNC80 and persistent hypotonia, encephalopathy, growth failure, and severe intellectual disability. Compared to control cells, HEK293T cells transfected with an expression plasmid containing the c.5098C>T (p.Pro1700Ser) UNC80 mutation found in one individual showed markedly decreased NALCN channel currents. Our findings demonstrate the fundamental significance of UNC80 and basal ionic conductance to human health.

Exome Sequencing and the Management of Neurometabolic Disorders.

BACKGROUND: Whole-exome sequencing has transformed gene discovery and diagnosis in rare diseases. Translation into disease-modifying treatments is challenging, particularly for intellectual developmental disorder. However, the exception is inborn errors of metabolism, since many of these disorders are responsive to therapy that targets pathophysiological features at the molecular or cellular level. METHODS: To uncover the genetic basis of potentially treatable inborn errors of metabolism, we combined deep clinical phenotyping (the comprehensive characterization of the discrete components of a patient's clinical and biochemical phenotype) with whole-exome sequencing analysis through a semiautomated bioinformatics pipeline in consecutively enrolled patients with intellectual developmental disorder and unexplained metabolic phenotypes. RESULTS: We performed whole-exome sequencing on samples obtained from 47 probands. Of these patients, 6 were excluded, including 1 who withdrew from the study. The remaining 41 probands had been born to predominantly nonconsanguineous parents of European descent. In 37 probands, we identified variants in 2 genes newly implicated in disease, 9 candidate genes, 22 known genes with newly identified phenotypes, and 9 genes with expected phenotypes; in most of the genes, the variants were classified as either pathogenic or probably pathogenic. Complex phenotypes of patients in five families were explained by coexisting monogenic conditions. We obtained a diagnosis in 28 of 41 probands (68%) who were evaluated. A test of a targeted intervention was performed in 18 patients (44%). CONCLUSIONS: Deep phenotyping and whole-exome sequencing in 41 probands with intellectual developmental disorder and unexplained metabolic abnormalities led to a diagnosis in 68%, the identification of 11 candidate genes newly implicated in neurometabolic disease, and a change in treatment beyond genetic counseling in 44%. (Funded by BC Children's Hospital Foundation and others.).

Atypical cerebral palsy: genomics analysis enables precision medicine.

PURPOSE: The presentation and etiology of cerebral palsy (CP) are heterogeneous. Diagnostic evaluation can be a prolonged and expensive process that might remain inconclusive. This study aimed to determine the diagnostic yield and impact on management of next-generation sequencing (NGS) in 50 individuals with atypical CP (ACP). METHODS: Patient eligibility criteria included impaired motor function with onset at birth or within the first year of life, and one or more of the following: severe intellectual disability, progressive neurological deterioration, other abnormalities on neurological examination, multiorgan disease, congenital anomalies outside of the central nervous system, an abnormal neurotransmitter profile, family history, brain imaging findings not typical for cerebral palsy. Previous assessment by a neurologist and/or clinical geneticist, including biochemical testing, neuroimaging, and chromosomal microarray, did not yield an etiologic diagnosis. RESULTS: A precise molecular diagnosis was established in 65% of the 50 patients. We also identified candidate disease genes without a current OMIM disease designation. Targeted intervention was enabled in eight families (~15%). CONCLUSION: NGS enabled a molecular diagnosis in ACP cases, ending the diagnostic odyssey, improving genetic counseling and personalized management, all in all enhancing precision medicine practices.

Genome sequencing reveals a novel genetic mechanism underlying dihydropyrimidine dehydrogenase deficiency: A novel missense variant c.1700G>A and a large intragenic inversion in DPYD spanning intron 8 to intron 12.

Dihydropyrimidine dehydrogenase (DPD) deficiency is associated with a variable clinical presentation. A family with three DPD-deficient patients presented with unusual clinical phenotypes including pregnancy-induced symptoms, transient visual impairment, severe developmental delay, cortical blindness, and delayed myelination in the brain. DPYD Sanger sequencing showed heterozygosity for the c.1905+1G>A mutation and a novel missense variant c.1700G>A (p.G567E). The recombinantly expressed p.G567E DPD variant showed increased temperature lability probably caused by structural rearrangements within the DPD protein. Genome sequencing of the affected son established compound heterozygosity for the c.1700G>A and an imperfect 115,731 bp inversion with breakpoints at chr1: 98,113,121 (intron 8) and chr1: 97,997,390 (intron 12) of the DPYD associated with a 4 bp deletion (chr1: 97,997,386_97,997,389del). Whole exome and mitochondrial DNA analyses for the mother and daughter did not reveal additional mutated genes of significance. Thus, an inversion in DPYD should be considered in patients with an inconclusive genotype or unusual clinical phenotype.

The genotypic and phenotypic spectrum of MTO1 deficiency.

BACKGROUND: Mitochondrial diseases, a group of multi-systemic disorders often characterized by tissue-specific phenotypes, are usually progressive and fatal disorders resulting from defects in oxidative phosphorylation. MTO1 (Mitochondrial tRNA Translation Optimization 1), an evolutionarily conserved protein expressed in high-energy demand tissues has been linked to human early-onset combined oxidative phosphorylation deficiency associated with hypertrophic cardiomyopathy, often referred to as combined oxidative phosphorylation deficiency-10 (COXPD10). MATERIAL AND METHODS: Thirty five cases of MTO1 deficiency were identified and reviewed through international collaboration. The cases of two female siblings, who presented at 1 and 2years of life with seizures, global developmental delay, hypotonia, elevated lactate and complex I and IV deficiency on muscle biopsy but without cardiomyopathy, are presented in detail. RESULTS: For the description of phenotypic features, the denominator varies as the literature was insufficient to allow for complete ascertainment of all data for the 35 cases. An extensive review of all known MTO1 deficiency cases revealed the most common features at presentation to be lactic acidosis (LA) (21/34; 62% cases) and hypertrophic cardiomyopathy (15/34; 44% cases). Eventually lactic acidosis and hypertrophic cardiomyopathy are described in 35/35 (100%) and 27/34 (79%) of patients with MTO1 deficiency, respectively; with global developmental delay/intellectual disability present in 28/29 (97%), feeding difficulties in 17/35 (49%), failure to thrive in 12/35 (34%), seizures in 12/35 (34%), optic atrophy in 11/21 (52%) and ataxia in 7/34 (21%). There are 19 different pathogenic MTO1 variants identified in these 35 cases: one splice-site, 3 frameshift and 15 missense variants. None have bi-allelic variants that completely inactivate MTO1; however, patients where one variant is truncating (i.e. frameshift) while the second one is a missense appear to have a more severe, even fatal, phenotype. These data suggest that complete loss of MTO1 is not viable. A ketogenic diet may have exerted a favourable effect on seizures in 2/5 patients. CONCLUSION: MTO1 deficiency is lethal in some but not all cases, and a genotype-phenotype relation is suggested. Aside from lactic acidosis and cardiomyopathy, developmental delay and other phenotypic features affecting multiple organ systems are often present in these patients, suggesting a broader spectrum than hitherto reported. The diagnosis should be suspected on clinical features and the presence of markers of mitochondrial dysfunction in body fluids, especially low residual complex I, III and IV activity in muscle. Molecular confirmation is required and targeted genomic testing may be the most efficient approach. Although subjective clinical improvement was observed in a small number of patients on therapies such as ketogenic diet and dichloroacetate, no evidence-based effective therapy exists.

The role of the clinician in the multi-omics era: are you ready?

Since Garrod's first description of alkaptonuria in 1902, and newborn screening for phenylketonuria introduced in the 1960s, P4 medicine (preventive, predictive, personalized, and participatory) has been a reality for the clinician serving patients with inherited metabolic diseases. The era of high-throughput technologies promises to accelerate its scale dramatically. Genomics, transcriptomics, epigenomics, proteomics, glycomics, metabolomics, and lipidomics offer an amazing opportunity for holistic investigation and contextual pathophysiologic understanding of inherited metabolic diseases for precise diagnosis and tailored treatment. While each of the -omics technologies is important to systems biology, some are more mature than others. Exome sequencing is emerging as a reimbursed test in clinics around the world, and untargeted metabolomics has the potential to serve as a single biochemical testing platform. The challenge lies in the integration and cautious interpretation of these big data, with translation into clinically meaningful information and/or action for our patients. A daunting but exciting task for the clinician; we provide clinical cases to illustrate the importance of his/her role as the connector between physicians, laboratory experts and researchers in the basic, computer, and clinical sciences. Open collaborations, data sharing, functional assays, and model organisms play a key role in the validation of -omics discoveries. Having all the right expertise at the table when discussing the diagnostic approach and individualized management plan according to the information yielded by -omics investigations (e.g., actionable mutations, novel therapeutic interventions), is the stepping stone of P4 medicine. Patient participation and the adjustment of the medical team's plan to his/her and the family's wishes most certainly is the capstone. Are you ready?

Integration of genomics and metabolomics for prioritization of rare disease variants: a 2018 literature review.

Many inborn errors of metabolism (IEMs) are amenable to treatment; therefore, early diagnosis and treatment is imperative. Despite recent advances, the genetic basis of many metabolic phenotypes remains unknown. For discovery purposes, whole exome sequencing (WES) variant prioritization coupled with clinical and bioinformatics expertise is the primary method used to identify novel disease-causing variants; however, causation is often difficult to establish due to the number of plausible variants. Integrated analysis of untargeted metabolomics (UM) and WES or whole genome sequencing (WGS) data is a promising systematic approach for identifying disease-causing variants. In this review, we provide a literature-based overview of UM methods utilizing liquid chromatography mass spectrometry (LC-MS), and assess approaches to integrating WES/WGS and LC-MS UM data for the discovery and prioritization of variants causing IEMs. To embed this integrated -omics approach in the clinic, expansion of gene-metabolite annotations and metabolomic feature-to-metabolite mapping methods are needed.