Search

Genetic balancers inCaenorhabditis elegansare 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 theC. elegansscientific community and have been widely used for decades. The first/traditional balancers were created by applying X-rays, UV, or gamma radiation onC. elegansstrains, 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 inC. elegansand highlights the power of srWGS to study the complexity of genomes by applying tailored analyses.

Genetic balancers in

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.

Complex genetic disease mechanisms, such as structural or non-coding variants, currently pose a substantial difficulty in frontline diagnostic tests. They thus may account for most unsolved rare disease patients regardless of the clinical phenotype. However, the clinical diagnosis can narrow the genetic focus to just a couple of genes for patients with well-established syndromes defined by prominent physical and/or unique biochemical phenotypes, allowing deeper analyses to consider complex genetic origin. Then, clinical-diagnosis-driven genome sequencing strategies may expedite the development of testing and analytical methods to account for complex disease mechanisms as well as to advance functional assays for the confirmation of complex variants, clinical management, and the development of new therapies.

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.

Genome-wide sequencing (GWS) is a standard of care for diagnosis of suspected genetic disorders, but the proportion of patients found to have pathogenic or likely pathogenic variants ranges from less than 30% to more than 60% in reported studies. It has been suggested that the diagnostic rate can be improved by interpreting genomic variants in the context of each affected individual's full clinical picture and by regular follow-up and reinterpretation of GWS laboratory results. Trio exome sequencing was performed in 415 families and trio genome sequencing in 85 families in the CAUSES study. The variants observed were interpreted by a multidisciplinary team including laboratory geneticists, bioinformaticians, clinical geneticists, genetic counselors, pediatric subspecialists, and the referring physician, and independently by a clinical laboratory using standard American College of Medical Genetics and Genomics (ACMG) criteria. Individuals were followed for an average of 5.1 years after testing, with clinical reassessment and reinterpretation of the GWS results as necessary. The multidisciplinary team established a diagnosis of genetic disease in 43.0% of the families at the time of initial GWS interpretation, and longitudinal follow-up and reinterpretation of GWS results produced new diagnoses in 17.2% of families whose initial GWS interpretation was uninformative or uncertain. Reinterpretation also resulted in rescinding a diagnosis in four families (1.9%). Of the families studied, 33.6% had ACMG pathogenic or likely pathogenic variants related to the clinical indication. Close collaboration among clinical geneticists, genetic counselors, laboratory geneticists, bioinformaticians, and individuals’ primary physicians, with ongoing follow-up, reanalysis, and reinterpretation over time, can improve the clinical value of GWS.

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Background The diagnostic journey for many rare disease patients remains challenging despite use of latest genetic technological advancements. We hypothesize that some patients remain undiagnosed due to more complex diagnostic scenarios that are currently not considered in genome analysis pipelines. To better understand this, we characterized the rare disorder (RD) spectrum using various bioinformatics resources (e.g., Orphanet/Orphadata, Human Phenotype Ontology, Reactome pathways) combined with custom-made R scripts. Results Our in silico characterization led to identification of 145 borderline-common, 412 rare and 2967 ultra-rare disorders. Based on these findings and point prevalence, we would expect that approximately 6.53%, 0.34%, and 0.30% of individuals in a randomly selected population have a borderline-common, rare, and ultra-rare disorder, respectively (equaling to 1 RD patient in 14 people). Importantly, our analyses revealed that (1) a higher proportion of borderline-common disorders were caused by multiple gene defects and/or other factors compared with the rare and ultra-rare disorders, (2) the phenotypic expressivity was more variable for the borderline-common disorders than for the rarer disorders, and (3) unique clinical characteristics were observed across the disorder categories forming the spectrum. Conclusions Recognizing that RD patients who remain unsolved even after genome sequencing might belong to the more common end of the RD spectrum support the usage of computational pipelines that account for more complex genetic and phenotypic scenarios.

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.

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.

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.

Many inborn errors of metabolism (IEMs) are amenable to treatment, therefore early diagnosis is imperative. Whole-exome sequencing (WES) variant prioritization coupled with phenotype-guided clinical and bioinformatics expertise is typically used to identify disease-causing variants; however, it can be challenging to identify the causal candidate gene when a large number of rare and potentially pathogenic variants are detected. Here, we present a network-based approach, metPropagate, that uses untargeted metabolomics (UM) data from a single patient and a group of controls to prioritize candidate genes in patients with suspected IEMs. We validate metPropagate on 107 patients with IEMs diagnosed in Miller et al. (2015) and 11 patients with both CNS and metabolic abnormalities. The metPropagate method ranks candidate genes by label propagation, a graph-smoothing algorithm that considers each gene’s metabolic perturbation in addition to the network of interactions between neighbors. metPropagate was able to prioritize at least one causative gene in the top 20th percentile of candidate genes for 92% of patients with known IEMs. Applied to patients with suspected neurometabolic disease, metPropagate placed at least one causative gene in the top 20th percentile in 9/11 patients, and ranked the causative gene more highly than Exomiser’s phenotype-based ranking in 6/11 patients. Interestingly, ranking by a weighted combination of metPropagate and Exomiser scores resulted in improved prioritization. The results of this study indicate that network-based analysis of UM data can provide an additional mode of evidence to prioritize causal genes in patients with suspected IEMs.

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.

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.

Here, we describe a fourth case of a human with a de novo KCNJ6 (GIRK2) mutation, who presented with clinical findings of severe hyperkinetic movement disorder and developmental delay, similar to the Keppen–Lubinsky syndrome but without lipodystrophy. Whole-exome sequencing of the patient’s DNA revealed a heterozygous de novo variant in the KCNJ6 (c.512T>G, p.Leu171Arg). We conducted in vitro functional studies to determine if this Leu-to-Arg mutation alters the function of GIRK2 channels. Heterologous expression of the mutant GIRK2 channel alone produced an aberrant basal inward current that lacked G protein activation, lost K+ selectivity and gained Ca2+ permeability. Notably, the inward current was inhibited by the Na+ channel blocker QX-314, similar to the previously reported weaver mutation in murine GIRK2. Expression of a tandem dimer containing GIRK1 and GIRK2(p.Leu171Arg) did not lead to any currents, suggesting heterotetramers are not functional. In neurons expressing p.Leu171Arg GIRK2 channels, these changes in channel properties would be expected to generate a sustained depolarization, instead of the normal G protein-gated inhibitory response, which could be mitigated by expression of other GIRK subunits. The identification of the p.Leu171Arg GIRK2 mutation potentially expands the Keppen–Lubinsky syndrome phenotype to include severe dystonia and ballismus. Our study suggests screening for dominant KCNJ6 mutations in the evaluation of patients with severe movement disorders, which could provide evidence to support a causal role of KCNJ6 in neurological channelopathies.

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.

Dystonic storm or status dystonicus is a life-threatening hyperkinetic movement disorder with biochemical alterations due to the excessive muscle contractions. The medical management can require pediatric intensive care unit admission and a combination of medications while the underlying trigger is managed. Severe cases may require general anesthesia and paralytic agents with intubation and may relapse when these drugs are weaned. Deep brain stimulation of the globus pallidum has been reported to terminate dystonic storm in several pediatric cases. We present a 10-year-old boy with a de novo GNAO1 mutation–induced dystonic storm who required a 2-month pediatric intensive care unit admission and remained refractory to all medical treatments. Deep brain stimulation was performed under general anesthetic without complication. His dyskinetic movements stopped with initiation of stimulation. He was discharged from the pediatric intensive care unit after 4 days. We present prospectively evaluated changes in dystonia symptoms and quality of life for a patient with GNAO1 mutation treated with deep brain stimulation.

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?

PAK3-related intellectual disability is caused by mutations in the gene encoding the p21-activated kinase (PAK) protein. It is characterized by mild to moderate cognitive impairment, micro/normocephaly, and a neurobehavioral phenotype characterized by short attention span, anxiety, restlessness, aggression, and self-abusive behaviors. The authors report a patient with a novel PAK3 mutation, who presented with intellectual disability, severe automutilation, and epilepsy. His magnetic resonance imaging changes were most likely secondary to lacerations from parenchymal contusions. His behavior was difficult to manage with behavior interventions or multiple medications. After finding low levels of dopamine and borderline low serotonin metabolites in the spinal fluid, treatment with low dose L-dopa/carbidopa and 5-hydroxytryptophan significantly improved his self-injurious behavior. This is the first case of PAK3-related intellectual disability presenting with severe self-injury with improvement following treatment. The patient’s response to neurotransmitter replacement therapy raises the question if this treatment intervention might help other individuals suffering genetic syndromes and self-injurious behaviors.

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.

Purpose: We analyzed the Exome Aggregation Consortium (ExAC) data set for the presence of individuals with pathogenic genotypes implicated in Mendelian pediatric disorders. Methods: ClinVar likely/pathogenic variants supported by at least one peer-reviewed publication were assessed within the ExAC database to identify individuals expected to exhibit a childhood disorder based on concordance with disease inheritance modes: heterozygous (for dominant), homozygous (for recessive) or hemizygous (for X-linked recessive conditions). Variants from 924 genes reported to cause Mendelian childhood disorders were considered. Results: We identified ExAC individuals with candidate pathogenic genotypes for 190 previously published likely/pathogenic variants in 128 genes. After curation, we determined that 113 of the variants have sufficient support for pathogenicity and identified 1,717 ExAC individuals (similar to 2.8% of the ExAC population) with corresponding possible/disease-associated genotypes implicated in rare Mendelian disorders, ranging from mild (e.g., due to SCN2A deficiency) to severe pediatric conditions (e.g., due to FGFR1 deficiency). Conclusion: Large-scale sequencing projects and data aggregation consortia provide unprecedented opportunities to determine the prevalence of pathogenic genotypes in unselected populations. This knowledge is crucial for understanding the penetrance of disease-associated variants, phenotypic variability, somatic mosaicism, as well as published literature curation for variant classification procedures and predicted clinical outcomes

Many inborn errors of metabolism (IEMs) are amenable to treatment, therefore early diagnosis before irreversible damage occurs 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 phenotype-guided clinical and bioinformatics expertise is currently the primary method used to identify novel disease-causing variants; however, it can be challenging to identify the causal candidate gene given the large number of plausible variants. Using untargeted metabolomics (UM) to prioritize metabolically relevant candidate genes is a promising approach to diagnosing known or novel IEMs in a single patient. Here, we present a network-based bioinformatics approach, metPropagate, that uses UM data from a single patient and a group of controls to prioritize candidate genes. We validate metProp on 107 patients with diagnosed IEMs and 11 patients with novel IEMs diagnosed through the TIDE gene discovery project at BC Children’s Hospital. The metPropagate method ranks candidate genes by considering the network of interactions between them. This is done by using a graph smoothing algorithm called label propagation to quantify the metabolic disruption in genes’ local neighbourhood. metPropagate was able to prioritize the causative gene in the top 20th percentile of candidate genes for 91% of patients with known IEM disorders. For novel IEMs, metPropagate placed the causative gene in the top 20th percentile in 9/11 patients. Using metPropagate, the causative gene was ranked higher than Exomiser’s phenotype-based ranking in 6/11 patients. The results of this study indicate that for diagnostic and gene discovery purposes, network-based analysis of metabolomics data can lend support to WES gene-discovery methods by providing an additional mode of evidence to help identify causal genes.

Background Profiling the entire genome at base pair resolution in a single test offers novel insights into disease by means of dissection of genetic contributors to phenotypic features. Methods We performed genome sequencing for a patient who presented with atypical hereditary sensory and autonomic neuropathy, severe epileptic encephalopathy, global developmental delay, and growth hormone deficiency. Results Assessment of the variants detected by mapped sequencing reads followed by Sanger confirmation revealed that the proband is a compound heterozygote for rare variants within RETREG1 (FAM134B), a gene associated with a recessive form of hereditary sensory and autonomic neuropathy, but not with epileptic encephalopathy or global developmental delay. Further analysis of the data also revealed a heterozygous missense variant in DNM1L, a gene previously implicated in an autosomal dominant encephalopathy, epilepsy, and global developmental delay and confirmed by Sanger sequencing to be a de novo variant not present in parental genomes. Conclusions Our findings emphasize the importance of genome‐wide sequencing in patients with a well‐characterized genetic disease with atypical presentation. This approach reduces the potential for misdiagnoses.

The problem of ‘missing heritability’ affects both common and rare diseases hindering: discovery, diagnosis, and patient care. The ‘missing heritability’ concept has been mainly associated with common and complex diseases where promising modern technological advances, like genome-wide association studies (GWAS), were unable to uncover the complete genetic mechanism of the disease/trait. Although rare diseases (RDs) have low prevalence individually, collectively they are common. Furthermore, multi-level genetic and phenotypic complexity when combined with the individual rarity of these conditions poses an important challenge in the quest to identify causative genetic changes in RD patients. In recent years, high throughput sequencing has accelerated discovery and diagnosis in RDs. However, despite the several-fold increase (from ~10% using traditional to ~40% using genome-wide genetic testing) in finding genetic causes of these diseases in RD patients, as is the case in common diseases—the majority of RDs are also facing the ‘missing heritability’ problem. This review outlines the key role of high throughput sequencing in uncovering genetics behind RDs, with a particular focus on genome sequencing. We review current advances and challenges of sequencing technologies, bioinformatics approaches, and resources.

Objectives: Clinicians need to be aware of the growing list of defined monogenic etiologies of autoimmune diseases. This is particularly relevant when evaluating children, as these rare monogenic forms of autoimmunity tend to present very early in life. Methods and results: By harnessing the transformative power of next generation sequencing, we made the unifying diagnosis of RAS-associated autoimmune leukoproliferative disease (RALD), caused by the somatic gain -of function p.G13C KRAS mutation, in a boy with the seemingly unrelated immune dysregulatory conditions of Rosai-Dorfman and systemic lupus erythematosus (SLE). Conclusions: This case expands our understanding of the clinical phenotypes associated with the extremely rare condition of RALD, and emphasizes the importance of always considering the possibility of a monogenic cause for autoimmunity, particularly when the disease manifestations begin early in life and do not follow a typical clinical course. Crown Copyright (C) 2016 Published by Elsevier Inc All rights reserved

Here we report a 12 year old male with an extreme presentation of spastic paraplegia along with autism and dysmorphisms. Whole exome sequencing identified a predicted pathogenic pair of missense variants in SPAST at the same chromosomal location, each with a different alternative allele, while a chromosome microarray identified a 1.73 Mb paternally inherited copy gain of 1q21.1q21.2 resulting in a blended phenotype of both Spastic paraplegia 4 and 1q21.1 microduplication syndrome. We believe that the extreme phenotype observed is likely caused by the presence of cells which contain only mutant SPAST, but that the viability of the patient is possible due mosaicism of mutant alleles observed in different proportions across tissues.

Introduction: Episodic ataxia (EA) is characterized by paroxysmal attacks of ataxia interspersed by asymptomatic periods. Dominant mutations or copy number variants in CACNA1A are a well-known cause of EA. Clinical presentation: This boy presented with clinical features of episodic ataxia, and also showed cerebellar atrophy, hypotonia, autism and global developmental delay at age 4 years. Acetazolamide prevented further episodes of ataxia, dystonia and encephalopathy. Extensive biochemical and genetic tests were unrevealing; whole exome sequencing found a previously unreported variant in SCN2A, proven to be de novo and predicted to be protein damaging. Conclusion: Considered alongside previous reports of episodic ataxia in SCN2A mutation positive patients, our case further illustrates the genetic heterogeneity of episodic ataxia. In addition, this case suggests that acetazolamide may be an effective treatment for some aspects of the phenotype in a broader range of channelopathy-related conditions. (C) 2016 European Paediatric Neurology Society. Published by Elsevier Ltd. All rights reserved

We report on 19 individuals with a recurrent de novo c.607C>T mutation in PACS1. This specific mutation gives rise to a recognizable intellectual disability syndrome. There is a distinctive facial appearance (19/19), characterized by full and arched eyebrows, hypertelorism with downslanting palpebral fissures, long eye lashes, ptosis, low set and simple ears, bulbous nasal tip, wide mouth with downturned corners and a thin upper lip with an unusual "wavy" profile, flat philtrum, and diastema of the teeth. Intellectual disability, ranging from mild to moderate, was present in all. Hypotonia is common in infancy (8/19). Seizures are frequent (12/19) and respond well to anticonvulsive medication. Structural malformations are common, including heart (10/19), brain (12/16), eye (10/19), kidney (3/19), and cryptorchidism (6/12 males). Feeding dysfunction is presenting in infancy with failure to thrive (5/19), gastroesophageal reflux (6/19), and gastrostomy tube placement (4/19). There is persistence of oral motor dysfunction. We provide suggestions for clinical work-up and management and hope that the present study will facilitate clinical recognition of further cases. Canadian Institutes of Health Research (301221) HICF-1009-003 WT098051 2.259 JCR (2016) Q3, 96/167 Genetics & Heredity UEM

Despite advances in high-throughput sequencing that have revolutionized the discovery of gene defects in rare Mendelian diseases, there are still gaps in translating individual genome variation to observed phenotypic outcomes. While we continue to improve genomics approaches to identify primary disease-causing variants, it is evident that no genetic variant acts alone. In other words, some other variants in the genome (genetic modifiers) may alleviate (suppress) or exacerbate (enhance) the severity of the disease, resulting in the variability of phenotypic outcomes. Thus, to truly understand the disease, we need to consider how the disease-causing variants interact with the rest of the genome in an individual. Here, we review the current state-of-the-field in the identification of genetic modifiers in rare Mendelian diseases and discuss the potential for future approaches that could bridge the existing gap.

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.

PEHO syndrome (OMIM no. 260565) is characterized by myoclonic jerking and infantile spasms, profound psychomotor retardation with the absence of motor milestones and speech, absence or early loss of visual fixation with atrophy of optic discs by 2 years of age and progressive brain atrophy on neuroimaging. We describe the results of a genomic study of a girl with PEHO syndrome and review the literature on cases with a disease-causing variant in the same gene. Exome sequencing of the index and unaffected parents followed by Sanger confirmation identified nine candidate genes harboring nonsynonymous rare variants identified by trio whole-exome sequencing. The de novo variant, a missense variant (c.296C>T, p.(T99M)), affecting the motor domain of KIF1A was considered the pathogenic mutation. The literature review revealed 24 cases with disease-causing variants in the motor domain of KIF1A, of which three met all the criteria for PEHO syndrome and an additional patient with incomplete clinical data met four of the five criteria. If the criteria were modified to include cases with any convulsive disorder and less profound intellectual disability, a total of six patients met all five of the criteria, three patients met four of the criteria and six met three of the criteria. Our results indicate that the molecular basis for PEHO syndrome, in at least a subset of patients, is a dominant KIF1A variant affecting the motor domain of the protein. Variable expressivity is seen with recurrent variants causing the full phenotype of PEHO syndrome in some patients and in other patients, a partial or milder PEHO phenotype.

Advances in next-generation sequencing (NGS) technologies have helped reveal causal variants for genetic diseases. In order to establish causality, it is often necessary to compare genomes of unrelated individuals with similar disease phenotypes to identify common disrupted genes. When working with cases of rare genetic disorders, finding similar individuals can be extremely difficult. We introduce a web tool, GeneYenta, which facilitates the matchmaking process, allowing clinicians to coordinate detailed comparisons for phenotypically similar cases. Importantly, the system is focused on phenotype annotation, with explicit limitations on highly confidential data that create barriers to participation. The procedure for matching of patient phenotypes, inspired by online dating services, uses an ontology­based semantic case matching algorithm with attribute weighting. We evaluate the capacity of the system using a curated reference data set and 19 clinician entered cases comparing four matching algorithms. We find that the inclusion of clinician weights can augment phenotype matching.

Sialic acids are important components of glycoproteins and glycolipids essential for cellular communication, infection, and metastasis. The importance of sialic acid biosynthesis in human physiology is well illustrated by the severe metabolic disorders in this pathway. However, the biological role of sialic acid catabolism in humans remains unclear. Here, we present evidence that sialic acid catabolism is important for heart and skeletal muscle function and development in humans and zebrafish. In two siblings, presenting with sialuria, exercise intolerance/muscle wasting, and cardiac symptoms in the brother, compound heterozygous mutations [chr1:182775324C>T (c.187C>T; p.Arg63Cys) and chr1:182772897A>G (c.133A>G; p.Asn45Asp)] were found in the N-acetylneuraminate pyruvate lyase gene (NPL). In vitro, NPL activity and sialic acid catabolism were affected, with a cell-type-specific reduction of N-acetyl mannosamine (ManNAc). A knockdown of NPL in zebrafish resulted in severe skeletal myopathy and cardiac edema, mimicking the human phenotype. The phenotype was rescued by expression of wild-type human NPL but not by the p.Arg63Cys or p.Asn45Asp mutants. Importantly, the myopathy phenotype in zebrafish embryos was rescued by treatment with the catabolic products of NPL: N-acetyl glucosamine (GlcNAc) and ManNAc; the latter also rescuing the cardiac phenotype. In conclusion, we provide the first report to our knowledge of a human defect in sialic acid catabolism, which implicates an important role of the sialic acid catabolic pathway in mammalian muscle physiology, and suggests opportunities for monosaccharide replacement therapy in human patients.

Objective The objective of this study was to investigate the genetic etiology of the X-linked disorder “Hypomyelination of Early Myelinating Structures” (HEMS). Methods We included 16 patients from 10 families diagnosed with HEMS by brain MRI criteria. Exome sequencing was used to search for causal mutations. In silico analysis of effects of the mutations on splicing and RNA folding was performed. In vitro gene splicing was examined in RNA from patients’ fibroblasts and an immortalized immature oligodendrocyte cell line after transfection with mutant minigene splicing constructs. Results All patients had unusual hemizygous mutations of PLP1 located in exon 3B (one deletion, one missense and two silent), which is spliced out in isoform DM20, or in intron 3 (five mutations). The deletion led to truncation of PLP1, but not DM20. Four mutations were predicted to affect PLP1/DM20 alternative splicing by creating exonic splicing silencer motifs or new splice donor sites or by affecting the local RNA structure of the PLP1 splice donor site. Four deep intronic mutations were predicted to destabilize a long-distance interaction structure in the secondary PLP1 RNA fragment involved in regulating PLP1/DM20 alternative splicing. Splicing studies in fibroblasts and transfected cells confirmed a decreased PLP1/DM20 ratio. Interpretation Brain structures that normally myelinate early are poorly myelinated in HEMS, while they are the best myelinated structures in Pelizaeus–Merzbacher disease, also caused by PLP1 alterations. Our data extend the phenotypic spectrum of PLP1-related disorders indicating that normal PLP1/DM20 alternative splicing is essential for early myelination and support the need to include intron 3 in diagnostic sequencing.

Background Dramatic improvements in DNA-sequencing technologies and computational analyses have led to wide use of whole exome sequencing (WES) to identify the genetic basis of Mendelian disorders. More than 180 novel rare-disease-causing genes with Mendelian inheritance patterns have been discovered through sequencing the exomes of just a few unrelated individuals or family members. As rare/novel genetic variants continue to be uncovered, there is a major challenge in distinguishing true pathogenic variants from rare benign mutations. Methods We used publicly available exome cohorts, together with the dbSNP database, to derive a list of genes (n = 100) that most frequently exhibit rare (

Corrigendum: NANS-mediated synthesis of sialic acid is required for brain and skeletal development

Contains fulltext : 202680.pdf (Publisher’s version ) (Closed access) 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.