Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus.

The 2013-2015 West African epidemic of Ebola virus disease (EVD) reminds us of how little is known about biosafety level 4 viruses. Like Ebola virus, Lassa virus (LASV) can cause hemorrhagic fever with high case fatality rates. We generated a genomic catalog of almost 200 LASV sequences from clinical and rodent reservoir samples. We show that whereas the 2013-2015 EVD epidemic is fueled by human-to-human transmissions, LASV infections mainly result from reservoir-to-human infections. We elucidated the spread of LASV across West Africa and show that this migration was accompanied by changes in LASV genome abundance, fatality rates, codon adaptation, and translational efficiency. By investigating intrahost evolution, we found that mutations accumulate in epitopes of viral surface proteins, suggesting selection for immune escape. This catalog will serve as a foundation for the development of vaccines and diagnostics. VIDEO ABSTRACT.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Advanced variant classification framework reduces the false positive rate of predicted loss-of-function variants in population sequencing data.

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Missense variants in DPYSL5 cause a neurodevelopmental disorder with corpus callosum agenesis and cerebellar abnormalities.

The collapsin response mediator protein (CRMP) family proteins are intracellular mediators of neurotrophic factors regulating neurite structure/spine formation and are essential for dendrite patterning and directional axonal pathfinding during brain developmental processes. Among this family, CRMP5/DPYSL5 plays a significant role in neuronal migration, axonal guidance, dendrite outgrowth, and synapse formation by interacting with microtubules. Here, we report the identification of missense mutations in DPYSL5 in nine individuals with brain malformations, including corpus callosum agenesis and/or posterior fossa abnormalities, associated with variable degrees of intellectual disability. A recurrent de novo p.Glu41Lys variant was found in eight unrelated patients, and a p.Gly47Arg variant was identified in one individual from the first family reported with Ritscher-Schinzel syndrome. Functional analyses of the two missense mutations revealed impaired dendritic outgrowth processes in young developing hippocampal primary neuronal cultures. We further demonstrated that these mutations, both located in the same loop on the surface of DPYSL5 monomers and oligomers, reduced the interaction of DPYSL5 with neuronal cytoskeleton-associated proteins MAP2 and ßIII-tubulin. Our findings collectively indicate that the p.Glu41Lys and p.Gly47Arg variants impair DPYSL5 function on dendritic outgrowth regulation by preventing the formation of the ternary complex with MAP2 and ßIII-tubulin, ultimately leading to abnormal brain development. This study adds DPYSL5 to the list of genes implicated in brain malformation and in neurodevelopmental disorders.

A form of muscular dystrophy associated with pathogenic variants in JAG2.

JAG2 encodes the Notch ligand Jagged2. The conserved Notch signaling pathway contributes to the development and homeostasis of multiple tissues, including skeletal muscle. We studied an international cohort of 23 individuals with genetically unsolved muscular dystrophy from 13 unrelated families. Whole-exome sequencing identified rare homozygous or compound heterozygous JAG2 variants in all 13 families. The identified bi-allelic variants include 10 missense variants that disrupt highly conserved amino acids, a nonsense variant, two frameshift variants, an in-frame deletion, and a microdeletion encompassing JAG2. Onset of muscle weakness occurred from infancy to young adulthood. Serum creatine kinase (CK) levels were normal or mildly elevated. Muscle histology was primarily dystrophic. MRI of the lower extremities revealed a distinct, slightly asymmetric pattern of muscle involvement with cores of preserved and affected muscles in quadriceps and tibialis anterior, in some cases resembling patterns seen in POGLUT1-associated muscular dystrophy. Transcriptome analysis of muscle tissue from two participants suggested misregulation of genes involved in myogenesis, including PAX7. In complementary studies, Jag2 downregulation in murine myoblasts led to downregulation of multiple components of the Notch pathway, including Megf10. Investigations in Drosophila suggested an interaction between Serrate and Drpr, the fly orthologs of JAG1/JAG2 and MEGF10, respectively. In silico analysis predicted that many Jagged2 missense variants are associated with structural changes and protein misfolding. In summary, we describe a muscular dystrophy associated with pathogenic variants in JAG2 and evidence suggests a disease mechanism related to Notch pathway dysfunction.

A recurrent missense variant in the E3 ubiquitin ligase substrate recognition subunit FEM1B causes a rare syndromic neurodevelopmental disorder.

PURPOSE: Fem1 homolog B (FEM1B) acts as a substrate recognition subunit for ubiquitin ligase complexes belonging to the CULLIN 2-based E3 family. Several biological functions have been proposed for FEM1B, including a structurally resolved function as a sensor for redox cell status by controlling mitochondrial activity, but its implication in human disease remains elusive. METHODS: To understand the involvement of FEM1B in human disease, we made use of Matchmaker exchange platforms to identify individuals with de novo variants in FEM1B and performed their clinical evaluation. We performed functional validation using primary neuronal cultures and in utero electroporation assays, as well as experiments on patient's cells. RESULTS: Five individuals with a recurrent de novo missense variant in FEM1B were identified: NM_015322.5:c.377G>A NP_056137.1:p.(Arg126Gln) (FEM1BR126Q). Affected individuals shared a severe neurodevelopmental disorder with behavioral phenotypes and a variable set of malformations, including brain anomalies, clubfeet, skeletal abnormalities, and facial dysmorphism. Overexpression of the FEM1BR126Q variant but not FEM1B wild-type protein, during mouse brain development, resulted in delayed neuronal migration of the target cells. In addition, the individuals' cells exhibited signs of oxidative stress and induction of type I interferon signaling. CONCLUSION: Overall, our data indicate that p.(Arg126Gln) induces aberrant FEM1B activation, resulting in a gain-of-function mechanism associated with a severe syndromic developmental disorder in humans.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

PURPOSE: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). METHODS: We coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. RESULTS: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. CONCLUSION: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

O'Donnell-Luria-Rodan syndrome: description of a second multinational cohort and refinement of the phenotypic spectrum.

BACKGROUND: O'Donnell-Luria-Rodan syndrome (ODLURO) is an autosomal-dominant neurodevelopmental disorder caused by pathogenic, mostly truncating variants in KMT2E. It was first described by O'Donnell-Luria et al in 2019 in a cohort of 38 patients. Clinical features encompass macrocephaly, mild intellectual disability (ID), autism spectrum disorder (ASD) susceptibility and seizure susceptibility. METHODS: Affected individuals were ascertained at paediatric and genetic centres in various countries by diagnostic chromosome microarray or exome/genome sequencing. Patients were collected into a case cohort and were systematically phenotyped where possible. RESULTS: We report 18 additional patients from 17 families with genetically confirmed ODLURO. We identified 15 different heterozygous likely pathogenic or pathogenic sequence variants (14 novel) and two partial microdeletions of KMT2E. We confirm and refine the phenotypic spectrum of the KMT2E-related neurodevelopmental disorder, especially concerning cognitive development, with rather mild ID and macrocephaly with subtle facial features in most patients. We observe a high prevalence of ASD in our cohort (41%), while seizures are present in only two patients. We extend the phenotypic spectrum by sleep disturbances. CONCLUSION: Our study, bringing the total of known patients with ODLURO to more than 60 within 2 years of the first publication, suggests an unexpectedly high relative frequency of this syndrome worldwide. It seems likely that ODLURO, although just recently described, is among the more common single-gene aetiologies of neurodevelopmental delay and ASD. We present the second systematic case series of patients with ODLURO, further refining the mutational and phenotypic spectrum of this not-so-rare syndrome.

seqr: A web-based analysis and collaboration tool for rare disease genomics.

Exome and genome sequencing have become the tools of choice for rare disease diagnosis, leading to large amounts of data available for analyses. To identify causal variants in these datasets, powerful filtering and decision support tools that can be efficiently used by clinicians and researchers are required. To address this need, we developed seqr - an open-source, web-based tool for family-based monogenic disease analysis that allows researchers to work collaboratively to search and annotate genomic callsets. To date, seqr is being used in several research pipelines and one clinical diagnostic lab. In our own experience through the Broad Institute Center for Mendelian Genomics, seqr has enabled analyses of over 10,000 families, supporting the diagnosis of more than 3,800 individuals with rare disease and discovery of over 300 novel disease genes. Here, we describe a framework for genomic analysis in rare disease that leverages seqr's capabilities for variant filtration, annotation, and causal variant identification, as well as support for research collaboration and data sharing. The seqr platform is available as open source software, allowing low-cost participation in rare disease research, and a community effort to support diagnosis and gene discovery in rare disease.

Recessive variants in COL25A1 gene as novel cause of arthrogryposis multiplex congenita with ocular congenital cranial dysinnervation disorder.

A proper interaction between muscle-derived collagen XXV and its motor neuron-derived receptors protein tyrosine phosphatases σ and δ (PTP σ/δ) is indispensable for intramuscular motor innervation. Despite this, thus far, pathogenic recessive variants in the COL25A1 gene had only been detected in a few patients with isolated ocular congenital cranial dysinnervation disorders. Here we describe five patients from three unrelated families with recessive missense and splice site COL25A1 variants presenting with a recognizable phenotype characterized by arthrogryposis multiplex congenita with or without an ocular congenital cranial dysinnervation disorder phenotype. The clinical features of the older patients remained stable over time, without central nervous system involvement. This study extends the phenotypic and genotypic spectrum of COL25A1 related conditions, and further adds to our knowledge of the complex process of intramuscular motor innervation. Our observations indicate a role for collagen XXV in regulating the appropriate innervation not only of extraocular muscles, but also of bulbar, axial, and limb muscles in the human.

The importance of automation in genetic diagnosis: Lessons from analyzing an inherited retinal degeneration cohort with the Mendelian Analysis Toolkit (MATK).

PURPOSE: In Mendelian disease diagnosis, variant analysis is a repetitive, error-prone, and time consuming process. To address this, we have developed the Mendelian Analysis Toolkit (MATK), a configurable, automated variant ranking program. METHODS: MATK aggregates variant information from multiple annotation sources and uses expert-designed rules with parameterized weights to produce a ranked list of potentially causal solutions. MATK performance was measured by a comparison between MATK-aided and human-domain expert analyses of 1060 families with inherited retinal degeneration (IRD), analyzed using an IRD-specific gene panel (589 individuals) and exome sequencing (471 families). RESULTS: When comparing MATK-assisted analysis with expert curation in both the IRD-specific gene panel and exome sequencing (1060 subjects), 97.3% of potential solutions found by experts were also identified by the MATK-assisted analysis (541 solutions identified with MATK of 556 solutions found by conventional analysis). Furthermore, MATK-assisted analysis identified 114 additional potential solutions from the 504 cases unsolved by conventional analysis. CONCLUSION: MATK expedites the process of identification of likely solving variants in Mendelian traits, and reduces variability stemming from human error and researcher bias. MATK facilitates data reanalysis to keep up with the constantly improving annotation sources and next-generation sequencing processing pipelines. The software is open source and available at https://gitlab.com/matthew_maher/mendelanalysis.

Major brain malformations: corpus callosum dysgenesis, agenesis of septum pellucidum and polymicrogyria in patients with BCORL1-related disorders.

OBJECTIVE: BCORL1, a transcriptional co-repressor, has a role in cortical migration, neuronal differentiation, maturation, and cerebellar development. We describe BCORL1 as a new genetic cause for major brain malformations. METHODS AND RESULTS: We report three patients from two unrelated families with neonatal onset intractable epilepsy and profound global developmental delay. Brain MRI of two siblings from the first family depicted hypoplastic corpus callosum and septal agenesis (ASP) in the older brother and unilateral perisylvian polymicrogyria (PMG) in the younger one. MRI of the patient from the second family demonstrated complete agenesis of corpus callosum (CC). Whole Exome Sequencing revealed a novel hemizygous variant in NM_021946.5 (BCORL1):c.796C>T (p.Pro266Ser) in the two siblings from the first family and the NM_021946.5 (BCORL1): c.3376G>A; p.Asp1126Asn variant in the patient from the second family, both variants inherited from healthy mothers. We reviewed the patients' charts and MRIs and compared the phenotype to the other published BCORL1-related cases. Brain malformations have not been previously described in association with the BCORL1 phenotype. We discuss the potential influence of BCORL1 on brain development. CONCLUSIONS: We suggest that BCORL1 variants present with a spectrum of neurodevelopmental disorders and can lead to major brain malformations originating at different stages of fetal development. We suggest adding BCORL1 to the genetic causes of PMG, ASP, and CC dysgenesis.

Cochlear nerve deficiency in SOX11-related Coffin-Siris syndrome.

The phenotypic spectrum of SOX11-related Coffin-Siris syndrome (CSS) is expanding with reports of new associations. SOX11 is implicated in neurogenesis and inner ear development. Cochlear nerve deficiency, absence or hypoplasia, is commonly associated with cochlear canal stenosis or with CHARGE syndrome, a monogenic condition that affects inner ear development. SOX11 is a transcription factor essential for neuronal identity, highly correlated with the expression of CHD7, which regulates SOX11. We present two unrelated probands, each with novel de novo SOX11 likely pathogenic variants and phenotypic manifestations of CSS including global developmental delay, growth deficiency, and hypoplastic nails. They have unilateral sensorineural hearing loss due to cochlear nerve deficiency confirmed on MRI. SOX11 is implicated in sensory neuron survival and maturation. It is highly expressed in the developing inner ear. Homozygous ablation of SOX11 in a mouse model resulted in a reduction in sensory neuron survival and decreased axonal growth. A heterozygous knockout mice model had hearing impairment with grossly normal inner ear structures like the two probands reported. We propose cochlear nerve deficiency as a new phenotypic feature of SOX11-related CSS. Magnetic resonance imaging is useful in delineating the cochlear nerve deficiency and other CSS-related brain malformations.

PIGN encephalopathy: Characterizing the epileptology.

OBJECTIVE: Epilepsy is common in patients with PIGN diseases due to biallelic variants; however, limited epilepsy phenotyping data have been reported. We describe the epileptology of PIGN encephalopathy. METHODS: We recruited patients with epilepsy due to biallelic PIGN variants and obtained clinical data regarding age at seizure onset/offset and semiology, development, medical history, examination, electroencephalogram, neuroimaging, and treatment. Seizure and epilepsy types were classified. RESULTS: Twenty six patients (13 female) from 26 families were identified, with mean age 7 years (range = 1 month to 21 years; three deceased). Abnormal development at seizure onset was present in 25 of 26. Developmental outcome was most frequently profound (14/26) or severe (11/26). Patients presented with focal motor (12/26), unknown onset motor (5/26), focal impaired awareness (1/26), absence (2/26), myoclonic (2/26), myoclonic-atonic (1/26), and generalized tonic-clonic (2/26) seizures. Twenty of 26 were classified as developmental and epileptic encephalopathy (DEE): 55% (11/20) focal DEE, 30% (6/20) generalized DEE, and 15% (3/20) combined DEE. Six had intellectual disability and epilepsy (ID+E): two generalized and four focal epilepsy. Mean age at seizure onset was 13 months (birth to 10 years), with a lower mean onset in DEE (7 months) compared with ID+E (33 months). Patients with DEE had drug-resistant epilepsy, compared to 4/6 ID+E patients, who were seizure-free. Hyperkinetic movement disorder occurred in 13 of 26 patients. Twenty-seven of 34 variants were novel. Variants were truncating (n = 7), intronic and predicted to affect splicing (n = 7), and missense or inframe indels (n = 20, of which 11 were predicted to affect splicing). Seven variants were recurrent, including p.Leu311Trp in 10 unrelated patients, nine with generalized seizures, accounting for nine of the 11 patients in this cohort with generalized seizures. SIGNIFICANCE: PIGN encephalopathy is a complex autosomal recessive disorder associated with a wide spectrum of epilepsy phenotypes, typically with substantial profound to severe developmental impairment.

Haploinsufficiency of ARFGEF1 is associated with developmental delay, intellectual disability, and epilepsy with variable expressivity.

PURPOSE: ADP ribosylation factor guanine nucleotide exchange factors (ARFGEFs) are a family of proteins implicated in cellular trafficking between the Golgi apparatus and the plasma membrane through vesicle formation. Among them is ARFGEF1/BIG1, a protein involved in axon elongation, neurite development, and polarization processes. ARFGEF1 has been previously suggested as a candidate gene for different types of epilepsies, although its implication in human disease has not been well characterized. METHODS: International data sharing, in silico predictions, and in vitro assays with minigene study, western blot analyses, and RNA sequencing. RESULTS: We identified 13 individuals with heterozygous likely pathogenic variants in ARFGEF1. These individuals displayed congruent clinical features of developmental delay, behavioral problems, abnormal findings on brain magnetic resonance image (MRI), and epilepsy for almost half of them. While nearly half of the cohort carried de novo variants, at least 40% of variants were inherited from mildly affected parents who were clinically re-evaluated by reverse phenotyping. Our in silico predictions and in vitro assays support the contention that ARFGEF1-related conditions are caused by haploinsufficiency, and are transmitted in an autosomal dominant fashion with variable expressivity. CONCLUSION: We provide evidence that loss-of-function variants in ARFGEF1 are implicated in sporadic and familial cases of developmental delay with or without epilepsy.

Exome sequencing identifies novel missense and deletion variants in RTN4IP1 associated with optic atrophy, global developmental delay, epilepsy, ataxia, and choreoathetosis.

Inherited optic neuropathies (IONs) are neurodegenerative disorders characterized by optic atrophy with or without extraocular manifestations. Optic atrophy-10 (OPA10) is an autosomal recessive ION recently reported to be caused by mutations in RTN4IP1, which encodes reticulon 4 interacting protein 1 (RTN4IP1), a mitochondrial ubiquinol oxydo-reductase. Here we report novel compound heterozygous mutations in RTN4IP1 in a male proband with developmental delay, epilepsy, optic atrophy, ataxia, and choreoathetosis. Workup was notable for transiently elevated lactate and lactate-to-pyruvate ratio, brain magnetic resonance imaging with optic atrophy and T2 signal abnormalities, and a nondiagnostic initial genetic workup, including chromosomal microarray and mitochondrial panel testing. Exome sequencing identified a paternally inherited missense variant (c.263T>G, p.Val88Gly) predicted to be deleterious and a maternally inherited deletion encompassing RTN4IP1. To our knowledge, this is the first report of a non-single nucleotide pathogenic variant associated with OPA10. This case highlights the expanding phenotypic spectrum of OPA10, the association between "syndromic" cases and severe RTN4IP1 mutations, and the importance of nonbiased genetic testing, such as ES, to analyze multiple genes and variants types, in patients suspected of having genetic disease.

Alternative genomic diagnoses for individuals with a clinical diagnosis of Dubowitz syndrome.

Dubowitz syndrome (DubS) is considered a recognizable syndrome characterized by a distinctive facial appearance and deficits in growth and development. There have been over 200 individuals reported with Dubowitz or a "Dubowitz-like" condition, although no single gene has been implicated as responsible for its cause. We have performed exome (ES) or genome sequencing (GS) for 31 individuals clinically diagnosed with DubS. After genome-wide sequencing, rare variant filtering and computational and Mendelian genomic analyses, a presumptive molecular diagnosis was made in 13/27 (48%) families. The molecular diagnoses included biallelic variants in SKIV2L, SLC35C1, BRCA1, NSUN2; de novo variants in ARID1B, ARID1A, CREBBP, POGZ, TAF1, HDAC8, and copy-number variation at1p36.11(ARID1A), 8q22.2(VPS13B), Xp22, and Xq13(HDAC8). Variants of unknown significance in known disease genes, and also in genes of uncertain significance, were observed in 7/27 (26%) additional families. Only one gene, HDAC8, could explain the phenotype in more than one family (N = 2). All but two of the genomic diagnoses were for genes discovered, or for conditions recognized, since the introduction of next-generation sequencing. Overall, the DubS-like clinical phenotype is associated with extensive locus heterogeneity and the molecular diagnoses made are for emerging clinical conditions sharing characteristic features that overlap the DubS phenotype.