A complete pedigree-based graph workflow for rare candidate variant analysis.

Methods that use a linear genome reference for genome sequencing data analysis are reference-biased. In the field of clinical genetics for rare diseases, a resulting reduction in genotyping accuracy in some regions has likely prevented the resolution of some cases. Pangenome graphs embed population variation into a reference structure. Although pangenome graphs have helped to reduce reference mapping bias, further performance improvements are possible. We introduce VG-Pedigree, a pedigree-aware workflow based on the pangenome-mapping tool of Giraffe and the variant calling tool DeepTrio using a specially trained model for Giraffe-based alignments. We demonstrate mapping and variant calling improvements in both single-nucleotide variants (SNVs) and insertion and deletion (indel) variants over those produced by alignments created using BWA-MEM to a linear-reference and Giraffe mapping to a pangenome graph containing data from the 1000 Genomes Project. We have also adapted and upgraded deleterious-variant (DV) detecting methods and programs into a streamlined workflow. We used these workflows in combination to detect small lists of candidate DVs among 15 family quartets and quintets of the Undiagnosed Diseases Program (UDP). All candidate DVs that were previously diagnosed using the Mendelian models covered by the previously published methods were recapitulated by these workflows. The results of these experiments indicate that a slightly greater absolute count of DVs are detected in the proband population than in their matched unaffected siblings.

FOXR1 regulates stress response pathways and is necessary for proper brain development.

The forkhead box (Fox) family of transcription factors are highly conserved and play essential roles in a wide range of cellular and developmental processes. We report an individual with severe neurological symptoms including postnatal microcephaly, progressive brain atrophy and global developmental delay associated with a de novo missense variant (M280L) in the FOXR1 gene. At the protein level, M280L impaired FOXR1 expression and induced a nuclear aggregate phenotype due to protein misfolding and proteolysis. RNAseq and pathway analysis showed that FOXR1 acts as a transcriptional activator and repressor with central roles in heat shock response, chaperone cofactor-dependent protein refolding and cellular response to stress pathways. Indeed, FOXR1 expression is increased in response to cellular stress, a process in which it directly controls HSPA6, HSPA1A and DHRS2 transcripts. The M280L mutant compromises FOXR1's ability to respond to stress, in part due to impaired regulation of downstream target genes that are involved in the stress response pathway. Quantitative PCR of mouse embryo tissues show Foxr1 expression in the embryonic brain. Using CRISPR/Cas9 gene editing, we found that deletion of mouse Foxr1 leads to a severe survival deficit while surviving newborn Foxr1 knockout mice have reduced body weight. Further examination of newborn Foxr1 knockout brains revealed a decrease in cortical thickness and enlarged ventricles compared to littermate wild-type mice, suggesting that loss of Foxr1 leads to atypical brain development. Combined, these results suggest FOXR1 plays a role in cellular stress response pathways and is necessary for normal brain development.

Lysosomal Storage and Albinism Due to Effects of a De Novo CLCN7 Variant on Lysosomal Acidification.

Optimal lysosome function requires maintenance of an acidic pH maintained by proton pumps in combination with a counterion transporter such as the Cl-/H+ exchanger, CLCN7 (ClC-7), encoded by CLCN7. The role of ClC-7 in maintaining lysosomal pH has been controversial. In this paper, we performed clinical and genetic evaluations of two children of different ethnicities. Both children had delayed myelination and development, organomegaly, and hypopigmentation, but neither had osteopetrosis. Whole-exome and -genome sequencing revealed a de novo c.2144A>G variant in CLCN7 in both affected children. This p.Tyr715Cys variant, located in the C-terminal domain of ClC-7, resulted in increased outward currents when it was heterologously expressed in Xenopus oocytes. Fibroblasts from probands displayed a lysosomal pH approximately 0.2 units lower than that of control cells, and treatment with chloroquine normalized the pH. Primary fibroblasts from both probands also exhibited markedly enlarged intracellular vacuoles; this finding was recapitulated by the overexpression of human p.Tyr715Cys CLCN7 in control fibroblasts, reflecting the dominant, gain-of-function nature of the variant. A mouse harboring the knock-in Clcn7 variant exhibited hypopigmentation, hepatomegaly resulting from abnormal storage, and enlarged vacuoles in cultured fibroblasts. Our results show that p.Tyr715Cys is a gain-of-function CLCN7 variant associated with developmental delay, organomegaly, and hypopigmentation resulting from lysosomal hyperacidity, abnormal storage, and enlarged intracellular vacuoles. Our data supports the hypothesis that the ClC-7 antiporter plays a critical role in maintaining lysosomal pH.

De Novo Variants in WDR37 Are Associated with Epilepsy, Colobomas, Dysmorphism, Developmental Delay, Intellectual Disability, and Cerebellar Hypoplasia.

WD40 repeat-containing proteins form a large family of proteins present in all eukaryotes. Here, we identified five pediatric probands with de novo variants in WDR37, which encodes a member of the WD40 repeat protein family. Two probands shared one variant and the others have variants in nearby amino acids outside the WD40 repeats. The probands exhibited shared phenotypes of epilepsy, colobomas, facial dysmorphology reminiscent of CHARGE syndrome, developmental delay and intellectual disability, and cerebellar hypoplasia. The WDR37 protein is highly conserved in vertebrate and invertebrate model organisms and is currently not associated with a human disease. We generated a null allele of the single Drosophila ortholog to gain functional insights and replaced the coding region of the fly gene CG12333/wdr37 with GAL4. These flies are homozygous viable but display severe bang sensitivity, a phenotype associated with seizures in flies. Additionally, the mutant flies fall when climbing the walls of the vials, suggesting a defect in grip strength, and repeat the cycle of climbing and falling. Similar to wall clinging defect, mutant males often lose grip of the female abdomen during copulation. These phenotypes are rescued by using the GAL4 in the CG12333/wdr37 locus to drive the UAS-human reference WDR37 cDNA. The two variants found in three human subjects failed to rescue these phenotypes, suggesting that these alleles severely affect the function of this protein. Taken together, our data suggest that variants in WDR37 underlie a novel syndromic neurological disorder.

Biallelic HEPHL1 variants impair ferroxidase activity and cause an abnormal hair phenotype.

Maintenance of the correct redox status of iron is functionally important for critical biological processes. Multicopper ferroxidases play an important role in oxidizing ferrous iron, released from the cells, into ferric iron, which is subsequently distributed by transferrin. Two well-characterized ferroxidases, ceruloplasmin (CP) and hephaestin (HEPH) facilitate this reaction in different tissues. Recently, a novel ferroxidase, Hephaestin like 1 (HEPHL1), also known as zyklopen, was identified. Here we report a child with compound heterozygous mutations in HEPHL1 (NM_001098672) who presented with abnormal hair (pili torti and trichorrhexis nodosa) and cognitive dysfunction. The maternal missense mutation affected mRNA splicing, leading to skipping of exon 5 and causing an in-frame deletion of 85 amino acids (c.809_1063del; p.Leu271_ala355del). The paternal mutation (c.3176T>C; p.Met1059Thr) changed a highly conserved methionine that is part of a typical type I copper binding site in HEPHL1. We demonstrated that HEPHL1 has ferroxidase activity and that the patient's two mutations exhibited loss of this ferroxidase activity. Consistent with these findings, the patient's fibroblasts accumulated intracellular iron and exhibited reduced activity of the copper-dependent enzyme, lysyl oxidase. These results suggest that the patient's biallelic variants are loss-of-function mutations. Hence, we generated a Hephl1 knockout mouse model that was viable and had curly whiskers, consistent with the hair phenotype in our patient. These results enhance our understanding of the function of HEPHL1 and implicate altered ferroxidase activity in hair growth and hair disorders.

Bi-allelic CCDC47 Variants Cause a Disorder Characterized by Woolly Hair, Liver Dysfunction, Dysmorphic Features, and Global Developmental Delay.

Ca2+ signaling is vital for various cellular processes including synaptic vesicle exocytosis, muscle contraction, regulation of secretion, gene transcription, and cellular proliferation. The endoplasmic reticulum (ER) is the largest intracellular Ca2+ store, and dysregulation of ER Ca2+ signaling and homeostasis contributes to the pathogenesis of various complex disorders and Mendelian disease traits. We describe four unrelated individuals with a complex multisystem disorder characterized by woolly hair, liver dysfunction, pruritus, dysmorphic features, hypotonia, and global developmental delay. Through whole-exome sequencing and family-based genomics, we identified bi-allelic variants in CCDC47 that encodes the Ca2+-binding ER transmembrane protein CCDC47. CCDC47, also known as calumin, has been shown to bind Ca2+ with low affinity and high capacity. In mice, loss of Ccdc47 leads to embryonic lethality, suggesting that Ccdc47 is essential for early development. Characterization of cells from individuals with predicted likely damaging alleles showed decreased CCDC47 mRNA expression and protein levels. In vitro cellular experiments showed decreased total ER Ca2+ storage, impaired Ca2+ signaling mediated by the IP3R Ca2+ release channel, and reduced ER Ca2+ refilling via store-operated Ca2+ entry. These results, together with the previously described role of CCDC47 in Ca2+ signaling and development, suggest that bi-allelic loss-of-function variants in CCDC47 underlie the pathogenesis of this multisystem disorder.

GGPS1 Mutations Cause Muscular Dystrophy/Hearing Loss/Ovarian Insufficiency Syndrome.

OBJECTIVE: A hitherto undescribed phenotype of early onset muscular dystrophy associated with sensorineural hearing loss and primary ovarian insufficiency was initially identified in 2 siblings and in subsequent patients with a similar constellation of findings. The goal of this study was to understand the genetic and molecular etiology of this condition. METHODS: We applied whole exome sequencing (WES) superimposed on shared haplotype regions to identify the initial biallelic variants in GGPS1 followed by GGPS1 Sanger sequencing or WES in 5 additional families with the same phenotype. Molecular modeling, biochemical analysis, laser membrane injury assay, and the generation of a Y259C knock-in mouse were done. RESULTS: A total of 11 patients in 6 families carrying 5 different biallelic pathogenic variants in specific domains of GGPS1 were identified. GGPS1 encodes geranylgeranyl diphosphate synthase in the mevalonate/isoprenoid pathway, which catalyzes the synthesis of geranylgeranyl pyrophosphate, the lipid precursor of geranylgeranylated proteins including small guanosine triphosphatases. In addition to proximal weakness, all but one patient presented with congenital sensorineural hearing loss, and all postpubertal females had primary ovarian insufficiency. Muscle histology was dystrophic, with ultrastructural evidence of autophagic material and large mitochondria in the most severe cases. There was delayed membrane healing after laser injury in patient-derived myogenic cells, and a knock-in mouse of one of the mutations (Y259C) resulted in prenatal lethality. INTERPRETATION: The identification of specific GGPS1 mutations defines the cause of a unique form of muscular dystrophy with hearing loss and ovarian insufficiency and points to a novel pathway for this clinical constellation. ANN NEUROL 2020;88:332-347.

Compound heterozygous KCTD7 variants in progressive myoclonus epilepsy.

KCTD7 is a member of the potassium channel tetramerization domain-containing protein family and has been associated with progressive myoclonic epilepsy (PME), characterized by myoclonus, epilepsy, and neurological deterioration. Here we report four affected individuals from two unrelated families in which we identified KCTD7 compound heterozygous single nucleotide variants through exome sequencing. RNAseq was used to detect a non-annotated splicing junction created by a synonymous variant in the second family. Whole-cell patch-clamp analysis of neuroblastoma cells overexpressing the patients' variant alleles demonstrated aberrant potassium regulation. While all four patients experienced many of the common clinical features of PME, they also showed variable phenotypes not previously reported, including dysautonomia, brain pathology findings including a significantly reduced thalamus, and the lack of myoclonic seizures. To gain further insight into the pathogenesis of the disorder, zinc finger nucleases were used to generate kctd7 knockout zebrafish. Kctd7 homozygous mutants showed global dysregulation of gene expression and increased transcription of c-fos, which has previously been correlated with seizure activity in animal models. Together these findings expand the known phenotypic spectrum of KCTD7-associated PME, report a new animal model for future studies, and contribute valuable insights into the disease.

A Syndromic Neurodevelopmental Disorder Caused by De Novo Variants in EBF3.

Early B cell factor 3 (EBF3) is a member of the highly evolutionarily conserved Collier/Olf/EBF (COE) family of transcription factors. Prior studies on invertebrate and vertebrate animals have shown that EBF3 homologs are essential for survival and that loss-of-function mutations are associated with a range of nervous system developmental defects, including perturbation of neuronal development and migration. Interestingly, aristaless-related homeobox (ARX), a homeobox-containing transcription factor critical for the regulation of nervous system development, transcriptionally represses EBF3 expression. However, human neurodevelopmental disorders related to EBF3 have not been reported. Here, we describe three individuals who are affected by global developmental delay, intellectual disability, and expressive speech disorder and carry de novo variants in EBF3. Associated features seen in these individuals include congenital hypotonia, structural CNS malformations, ataxia, and genitourinary abnormalities. The de novo variants affect a single conserved residue in a zinc finger motif crucial for DNA binding and are deleterious in a fly model. Our findings indicate that mutations in EBF3 cause a genetic neurodevelopmental syndrome and suggest that loss of EBF3 function might mediate a subset of neurologic phenotypes shared by ARX-related disorders, including intellectual disability, abnormal genitalia, and structural CNS malformations.

Expanding the phenotype of COPA syndrome: a kindred with typical and atypical features.

BACKGROUND: Copa syndrome is a rare autosomal dominant disorder with abnormal intracellular vesicle trafficking. The objective of this work is to expand the knowledge about this disorder by delineating phenotypic features of an unreported COPA family. METHODS AND RESULTS: A heterozygous missense variant (c.698 G>A, p.Arg233His) in COPA was identified in four members of a three-generation kindred with lung, autoimmune and malignant disease of unknown aetiology. Ages of onset were 56, 26, 16 and 1 year, with earlier age of onset in successive generations. Presenting symptoms were cough and dyspnoea. Findings included small lung cysts, follicular bronchiolitis, interstitial lung disease, neuroendocrine cell hyperplasia, rheumatoid arthritis, avascular necrosis and select abnormal autoimmune serologies. Neither alveolar haemorrhage nor glomerular disease were present. Features not previously associated with Copa syndrome included neuromyelitis optica, pulmonary carcinoid tumour, clear cell renal carcinoma, renal cysts, hepatic cysts, nephrolithiasis, pyelonephritis and meningitis. Longitudinal evaluations demonstrated slow progression of lung disease and extrapulmonary cysts. CONCLUSIONS: Worsening severity with successive generations may be observed in Copa syndrome. Extrapulmonary cysts, malignancies, autoimmune neurological disorders and infections are clinical features that may be associated with Copa syndrome. Further studies are indicated to fully define the phenotypic spectrum of this disorder.

Early infantile-onset epileptic encephalopathy 28 due to a homozygous microdeletion involving the WWOX gene in a region of uniparental disomy.

The genetic etiologies of many rare disorders, including early infantile epileptic encephalopathies, are largely undiagnosed. A 6-year-old girl was admitted to the National Institutes of Health Undiagnosed Diseases Program with profound intellectual disability, infantile-onset seizures, chronic respiratory failure, facial dysmorphisms, skeletal abnormalities, and atrial septum defect. A large region of homozygosity was discovered on chromosome 16, spanning 16q22.1-16q24.3' caused by uniparental disomy (UPD) that included a maternally inherited homozygous microdeletion covering exon 6 of WWOX (NM_016373.3). mRNA expression analysis revealed that the deletion led to nonsense-mediated decay of the NM_016373.3 transcript; the exon 6 of an alternative transcript (NM_130791.3), lacking the short-chain dehydrogenase, was utilized. The microdeletion in WWOX explains the seizures and intellectual disability, while pathogenic variants in another gene, HSPG2, are likely responsible for the patient's skeletal abnormalities. This report describes a rare autosomal recessive disorder with multiple genetic etiologies, one of which involves UPD.

A suite of automated sequence analyses reduces the number of candidate deleterious variants and reveals a difference between probands and unaffected siblings.

PURPOSE: Develop an automated exome analysis workflow that can produce a very small number of candidate variants yet still detect different numbers of deleterious variants between probands and unaffected siblings. METHODS: Ninety-seven outbred nuclear families from the Undiagnosed Diseases Program/Network included single probands and the corresponding unaffected sibling(s). Single-nucleotide polymorphism (SNP) chip and exome analyses were performed on all, with proband and unaffected sibling considered independently as the target. The total burden of candidate genetic variants was summed for probands and siblings over all considered disease models. RESULTS: Exome analysis workflow include automated programs for ethnicity-matched genotype calling, salvage pathway for Mendelian inconsistency, compound heterozygous recessive detection, BAM file regional curation, population frequency filtering, pedigree-aware BAM file noise evaluation, and exon deletion filtration. This workflow relied heavily on BAM file analysis. A greater average pathogenic variant number was found compared with unaffected siblings. This was significant (p < 0.05) when using published recommended thresholds, and implies that causal variants are retained in many probands' lists. CONCLUSION: Using Mendelian and non-Mendelian models, this agnostic exome analysis shows a difference between a small group of probands and their unaffected siblings. This workflow produces candidate lists small enough to pursue with laboratory validation.

KCTD7 deficiency defines a distinct neurodegenerative disorder with a conserved autophagy-lysosome defect.

OBJECTIVE: Several small case series identified KCTD7 mutations in patients with a rare autosomal recessive disorder designated progressive myoclonic epilepsy (EPM3) and neuronal ceroid lipofuscinosis (CLN14). Despite the name KCTD (potassium channel tetramerization domain), KCTD protein family members lack predicted channel domains. We sought to translate insight gained from yeast studies to uncover disease mechanisms associated with deficiencies in KCTD7 of unknown function. METHODS: Novel KCTD7 variants in new and published patients were assessed for disease causality using genetic analyses, cell-based functional assays of patient fibroblasts and knockout yeast, and electron microscopy of patient samples. RESULTS: Patients with KCTD7 mutations can exhibit movement disorders or developmental regression before seizure onset, and are distinguished from similar disorders by an earlier age of onset. Although most published KCTD7 patient variants were excluded from a genome sequence database of normal human variations, most newly identified patient variants are present in this database, potentially challenging disease causality. However, genetic analysis and impaired biochemical interactions with cullin 3 support a causal role for patient KCTD7 variants, suggesting deleterious alleles of KCTD7 and other rare disease variants may be underestimated. Both patient-derived fibroblasts and yeast lacking Whi2 with sequence similarity to KCTD7 have impaired autophagy consistent with brain pathology. INTERPRETATION: Biallelic KCTD7 mutations define a neurodegenerative disorder with lipofuscin and lipid droplet accumulation but without defining features of neuronal ceroid lipofuscinosis or lysosomal storage disorders. KCTD7 deficiency appears to cause an underlying autophagy-lysosome defect conserved in yeast, thereby assigning a biological role for KCTD7. Ann Neurol 2018;84:774-788.

Recurrent Mutations in the Basic Domain of TWIST2 Cause Ablepharon Macrostomia and Barber-Say Syndromes.

Ablepharon macrostomia syndrome (AMS) and Barber-Say syndrome (BSS) are rare congenital ectodermal dysplasias characterized by similar clinical features. To establish the genetic basis of AMS and BSS, we performed extensive clinical phenotyping, whole exome and candidate gene sequencing, and functional validations. We identified a recurrent de novo mutation in TWIST2 in seven independent AMS-affected families, as well as another recurrent de novo mutation affecting the same amino acid in ten independent BSS-affected families. Moreover, a genotype-phenotype correlation was observed, because the two syndromes differed based solely upon the nature of the substituting amino acid: a lysine at TWIST2 residue 75 resulted in AMS, whereas a glutamine or alanine yielded BSS. TWIST2 encodes a basic helix-loop-helix transcription factor that regulates the development of mesenchymal tissues. All identified mutations fell in the basic domain of TWIST2 and altered the DNA-binding pattern of Flag-TWIST2 in HeLa cells. Comparison of wild-type and mutant TWIST2 expressed in zebrafish identified abnormal developmental phenotypes and widespread transcriptome changes. Our results suggest that autosomal-dominant TWIST2 mutations cause AMS or BSS by inducing protean effects on the transcription factor's DNA binding.

Novel variants in SPTAN1 without epilepsy: An expansion of the phenotype.

We describe two unrelated children with de novo variants in the non-erythrocytic alpha-II-spectrin (SPTAN1) gene who have hypoplastic brain structures, intellectual disability, and both fine and gross motor impairments. Using agnostic exome sequencing, we identified a nonsense variant creating a premature stop codon in exon 21 of SPTAN1, and in a second patient we identified an intronic substitution in SPTAN1 prior to exon 50 creating a new donor acceptor site. Neither of these variants has been described previously. Although some of these patients' features are consistent with the known SPTAN1 encephalopathy phenotype, these two children do not have epilepsy, in contrast to reports about nearly every other patient with heterozygous SPTAN1 variants and in all patients with a variant near the C-terminal coding region. Moreover, both children have abnormal thyroid function, which has not been previously reported in association with SPTAN1 variant. We present a detailed discussion of the clinical manifestations of these two unique SPTAN1 variants and provide evidence that both variants result in reduced mRNA expression despite different locations within the gene and clinical phenotypes. These findings expand the motor, cognitive, and behavioral spectrum of the SPTAN1-associated phenotype and invite speculation about underlying pathophysiologies.

Compound heterozygosity for loss-of-function GARS variants results in a multisystem developmental syndrome that includes severe growth retardation.

Aminoacyl-tRNA synthetases (ARSs) are ubiquitously expressed enzymes that ligate amino acids onto tRNA molecules. Genes encoding ARSs have been implicated in myriad dominant and recessive disease phenotypes. Glycyl-tRNA synthetase (GARS) is a bifunctional ARS that charges tRNAGly in the cytoplasm and mitochondria. GARS variants have been associated with dominant Charcot-Marie-Tooth disease but have not been convincingly implicated in recessive phenotypes. Here, we describe a patient from the NIH Undiagnosed Diseases Program with a multisystem, developmental phenotype. Whole-exome sequence analysis revealed that the patient is compound heterozygous for one frameshift (p.Glu83Ilefs*6) and one missense (p.Arg310Gln) GARS variant. Using in vitro and in vivo functional studies, we show that both GARS variants cause a loss-of-function effect: the frameshift variant results in depleted protein levels and the missense variant reduces GARS tRNA charging activity. In support of GARS variant pathogenicity, our patient shows striking phenotypic overlap with other patients having ARS-related recessive diseases, including features associated with variants in both cytoplasmic and mitochondrial ARSs; this observation is consistent with the essential function of GARS in both cellular locations. In summary, our clinical, genetic, and functional analyses expand the phenotypic spectrum associated with GARS variants.

Biallelic mutations in CAD, impair de novo pyrimidine biosynthesis and decrease glycosylation precursors.

In mitochondria, carbamoyl-phosphate synthetase 1 activity produces carbamoyl phosphate for urea synthesis, and deficiency results in hyperammonemia. Cytoplasmic carbamoyl-phosphate synthetase 2, however, is part of a tri-functional enzyme encoded by CAD; no human disease has been attributed to this gene. The tri-functional enzyme contains carbamoyl-phosphate synthetase 2 (CPS2), aspartate transcarbamylase (ATCase) and dihydroorotase (DHOase) activities, which comprise the first three of six reactions required for de novo pyrimidine biosynthesis. Here we characterize an individual who is compound heterozygous for mutations in different domains of CAD. One mutation, c.1843-1G>A, results in an in-frame deletion of exon 13. The other, c.6071G>A, causes a missense mutation (p.Arg2024Gln) in a highly conserved residue that is essential for carbamoyl-phosphate binding. Metabolic flux studies showed impaired aspartate incorporation into RNA and DNA through the de novo synthesis pathway. In addition, CTP, UTP and nearly all UDP-activated sugars that serve as donors for glycosylation were decreased. Uridine supplementation rescued these abnormalities, suggesting a potential therapy for this new glycosylation disorder.

Combined alpha-delta platelet storage pool deficiency is associated with mutations in GFI1B.

Combined alpha-delta platelet storage pool deficiency is characterized by the absence or reduction in the number of both alpha granules and dense bodies. This disorder can have variable severity as well as a variable inheritance pattern. We describe two patients from unrelated families with combined alpha-delta storage pool deficiency due to mutations in GFI1B, a zinc finger protein known to act as a transcriptional repressor of various genes. We demonstrate that this disease is associated with either a heterozygous mutation (de novo or familial) abrogating the binding of the zinc fingers with the promoter of its target genes, or by hypomorphic biallelic mutations in GFI1B leading to autosomal recessive inheritance.

Computational evaluation of exome sequence data using human and model organism phenotypes improves diagnostic efficiency.

PURPOSE: Medical diagnosis and molecular or biochemical confirmation typically rely on the knowledge of the clinician. Although this is very difficult in extremely rare diseases, we hypothesized that the recording of patient phenotypes in Human Phenotype Ontology (HPO) terms and computationally ranking putative disease-associated sequence variants improves diagnosis, particularly for patients with atypical clinical profiles. METHODS: Using simulated exomes and the National Institutes of Health Undiagnosed Diseases Program (UDP) patient cohort and associated exome sequence, we tested our hypothesis using Exomiser. Exomiser ranks candidate variants based on patient phenotype similarity to (i) known disease-gene phenotypes, (ii) model organism phenotypes of candidate orthologs, and (iii) phenotypes of protein-protein association neighbors. RESULTS: Benchmarking showed Exomiser ranked the causal variant as the top hit in 97% of known disease-gene associations and ranked the correct seeded variant in up to 87% when detectable disease-gene associations were unavailable. Using UDP data, Exomiser ranked the causative variant(s) within the top 10 variants for 11 previously diagnosed variants and achieved a diagnosis for 4 of 23 cases undiagnosed by clinical evaluation. CONCLUSION: Structured phenotyping of patients and computational analysis are effective adjuncts for diagnosing patients with genetic disorders.Genet Med 18 6, 608-617.