Search

Discovering the genetic basis of a Mendelian phenotype establishes a causal link between genotype and phenotype, making possible carrier and population screening and direct diagnosis. Such discoveries also contribute to our knowledge of gene function, gene regulation, development, and biological mechanisms that can be used for developing new therapeutics. As of February 2015, 2,937 genes underlying 4,163 Mendelian phenotypes have been discovered, but the genes underlying ∼50% (i.e., 3,152) of all known Mendelian phenotypes are still unknown, and many more Mendelian conditions have yet to be recognized. This is a formidable gap in biomedical knowledge. Accordingly, in December 2011, the NIH established the Centers for Mendelian Genomics (CMGs) to provide the collaborative framework and infrastructure necessary for undertaking large-scale whole-exome sequencing and discovery of the genetic variants responsible for Mendelian phenotypes. In partnership with 529 investigators from 261 institutions in 36 countries, the CMGs assessed 18,863 samples from 8,838 families representing 579 known and 470 novel Mendelian phenotypes as of January 2015. This collaborative effort has identified 956 genes, including 375 not previously associated with human health, that underlie a Mendelian phenotype. These results provide insight into study design and analytical strategies, identify novel mechanisms of disease, and reveal the extensive clinical variability of Mendelian phenotypes. Discovering the gene underlying every Mendelian phenotype will require tackling challenges such as worldwide ascertainment and phenotypic characterization of families affected by Mendelian conditions, improvement in sequencing and analytical techniques, and pervasive sharing of phenotypic and genomic data among researchers, clinicians, and families. [ABSTRACT FROM AUTHOR]

Despite release of the GRCh38 human reference genome more than seven years ago, GRCh37 remains more widely used by most research and clinical laboratories. To date, no study has quantified the impact of utilizing different reference assemblies for the identification of variants associated with rare and common diseases from large-scale exome-sequencing data. By calling variants on both the GRCh37 and GRCh38 references, we identified single-nucleotide variants (SNVs) and insertion-deletions (indels) in 1,572 exomes from participants with Mendelian diseases and their family members. We found that a total of 1.5% of SNVs and 2.0% of indels were discordant when different references were used. Notably, 76.6% of the discordant variants were clustered within discrete discordant reference patches (DISCREPs) comprising only 0.9% of loci targeted by exome sequencing. These DISCREPs were enriched for genomic elements including segmental duplications, fix patch sequences, and loci known to contain alternate haplotypes. We identified 206 genes significantly enriched for discordant variants, most of which were in DISCREPs and caused by multi-mapped reads on the reference assembly that lacked the variant call. Among these 206 genes, eight are implicated in known Mendelian diseases and 53 are associated with common phenotypes from genome-wide association studies. In addition, variant interpretations could also be influenced by the reference after lifting-over variant loci to another assembly. Overall, we identified genes and genomic loci affected by reference assembly choice, including genes associated with Mendelian disorders and complex human diseases that require careful evaluation in both research and clinical applications. [ABSTRACT FROM AUTHOR]

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints. [ABSTRACT FROM AUTHOR]

Lissencephaly comprises a spectrum of malformations of cortical development. This spectrum includes agyria, pachygyria, and subcortical band heterotopia; each represents anatomical malformations of brain cortical development caused by neuronal migration defects. The molecular etiologies of neuronal migration anomalies are highly enriched for genes encoding microtubules and microtubule-associated proteins, and this enrichment highlights the critical role for these genes in cortical growth and gyrification. Using exome sequencing and family based rare variant analyses, we identified a homozygous variant (c.997C>T [p.Arg333Cys]) in TUBGCP2, encoding gamma-tubulin complex protein 2 (GCP2), in two individuals from a consanguineous family; both individuals presented with microcephaly and developmental delay. GCP2 forms the multiprotein γ-tubulin ring complex (γ-TuRC) together with γ-tubulin and other GCPs to regulate the assembly of microtubules. By querying clinical exome sequencing cases and through GeneMatcher-facilitated collaborations, we found three additional families with bi-allelic variation and similarly affected phenotypes including a homozygous variant (c.1843G>C [p.Ala615Pro]) in two families and compound heterozygous variants consisting of one missense variant (c.889C>T [p.Arg297Cys]) and one splice variant (c.2025-2A>G) in another family. Brain imaging from all five affected individuals revealed varying degrees of cortical malformations including pachygyria and subcortical band heterotopia, presumably caused by disruption of neuronal migration. Our data demonstrate that pathogenic variants in TUBGCP2 cause an autosomal recessive neurodevelopmental trait consisting of a neuronal migration disorder, and our data implicate GCP2 as a core component of γ-TuRC in neuronal migrating cells. [ABSTRACT FROM AUTHOR]

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders. [ABSTRACT FROM AUTHOR]

Premature termination codon (PTC)-bearing transcripts are often degraded by nonsense-mediated decay (NMD) resulting in loss-of-function (LoF) alleles. However, not all PTCs result in LoF mutations, i.e., some such transcripts escape NMD and are translated to truncated peptide products that result in disease due to gain-of-function (GoF) effects. Since the location of the PTC is a major factor determining transcript fate, we hypothesized that depletion of protein-truncating variants (PTVs) within the gene region predicted to escape NMD in control databases could provide a rank for genic susceptibility for disease through GoF versus LoF. We developed an NMD escape intolerance score to rank genes based on the depletion of PTVs that would render them able to escape NMD using the Atherosclerosis Risk in Communities Study (ARIC) and the Exome Aggregation Consortium (ExAC) control databases, which was further used to screen the Baylor-Center for Mendelian Genomics disease database. This analysis revealed 1,996 genes significantly depleted for PTVs that are predicted to escape from NMD, i.e., PTVesc; further studies provided evidence that revealed a subset as candidate genes underlying Mendelian phenotypes. Importantly, these genes have characteristically low pLI scores, which can cause them to be overlooked as candidates for dominant diseases. Collectively, we demonstrate that this NMD escape intolerance score is an effective and efficient tool for gene discovery in Mendelian diseases due to production of truncated or altered proteins. More importantly, we provide a complementary analytical tool to aid identification of genes associated with dominant traits through a mechanism distinct from LoF. [ABSTRACT FROM AUTHOR]

The proteasome processes proteins to facilitate immune recognition and host defense. When inherently defective, it can lead to aberrant immunity resulting in a dysregulated response that can cause autoimmunity and/or autoinflammation. Biallelic or digenic loss-of-function variants in some of the proteasome subunits have been described as causing a primary immunodeficiency disease that manifests as a severe dysregulatory syndrome: chronic atypical neutrophilic dermatosis with lipodystrophy and elevated temperature (CANDLE). Proteasome maturation protein (POMP) is a chaperone for proteasome assembly and is critical for the incorporation of catalytic subunits into the proteasome. Here, we characterize and describe POMP-related autoinflammation and immune dysregulation disease (PRAID) discovered in two unrelated individuals with a unique constellation of early-onset combined immunodeficiency, inflammatory neutrophilic dermatosis, and autoimmunity. We also begin to delineate a complex genetic mechanism whereby de novo heterozygous frameshift variants in the penultimate exon of POMP escape nonsense-mediated mRNA decay (NMD) and result in a truncated protein that perturbs proteasome assembly by a dominant-negative mechanism. To our knowledge, this mechanism has not been reported in any primary immunodeficiencies, autoinflammatory syndromes, or autoimmune diseases. Here, we define a unique hypo- and hyper-immune phenotype and report an immune dysregulation syndrome caused by frameshift mutations that escape NMD. [ABSTRACT FROM AUTHOR]

Hereditary gingival fibromatosis (HGF) is the most common genetic form of gingival fibromatosis that develops as a slowly progressive, benign, localized or generalized enlargement of keratinized gingiva. HGF is a genetically heterogeneous disorder and can be transmitted either as an autosomal-dominant or autosomal-recessive trait or appear sporadically. To date, four loci (2p22.1, 2p23.3–p22.3, 5q13–q22, and 11p15) have been mapped to autosomes and one gene ( SOS1 ) has been associated with the HGF trait observed to segregate in a dominant inheritance pattern. Here we report 11 individuals with HGF from three unrelated families. Whole-exome sequencing (WES) revealed three different truncating mutations including two frameshifts and one nonsense variant in RE1-silencing transcription factor ( REST ) in the probands from all families and further genetic and genomic analyses confirmed the WES-identified findings. REST is a transcriptional repressor that is expressed throughout the body; it has different roles in different cellular contexts, such as oncogenic and tumor-suppressor functions and hematopoietic and cardiac differentiation. Here we show the consequences of germline final-exon-truncating mutations in REST for organismal development and the association with the HGF phenotype. [ABSTRACT FROM AUTHOR]

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A , including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of bor R534W , the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of bor WT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2 , OPA1 , DNM1L , and STAT2 mutations. [ABSTRACT FROM AUTHOR]

Disruption of the establishment of left-right (L-R) asymmetry leads to situs anomalies ranging from situs inversus totalis (SIT) to situs ambiguus (heterotaxy). The genetic causes of laterality defects in humans are highly heterogeneous. Via whole-exome sequencing (WES), we identified homozygous mutations in PKD1L1 from three affected individuals in two unrelated families. PKD1L1 encodes a polycystin-1-like protein and its loss of function is known to cause laterality defects in mouse and medaka fish models. Family 1 had one fetus and one deceased child with heterotaxy and complex congenital heart malformations. WES identified a homozygous splicing mutation, c.6473+2_6473+3delTG, which disrupts the invariant splice donor site in intron 42, in both affected individuals. In the second family, a homozygous c.5072G>C (p.Cys1691Ser) missense mutation was detected in an individual with SIT and congenital heart disease. The p.Cys1691Ser substitution affects a highly conserved cysteine residue and is predicted by molecular modeling to disrupt a disulfide bridge essential for the proper folding of the G protein-coupled receptor proteolytic site (GPS) motif. Damaging effects associated with substitutions of this conserved cysteine residue in the GPS motif have also been reported in other genes, namely GPR56 , BAI3 , and PKD1 in human and lat-1 in C. elegans , further supporting the likely pathogenicity of p.Cys1691Ser in PKD1L1. The identification of bi-allelic PKD1L1 mutations recapitulates previous findings regarding phenotypic consequences of loss of function of the orthologous genes in mice and medaka fish and further expands our understanding of genetic contributions to laterality defects in humans. [ABSTRACT FROM AUTHOR]

Robinow syndrome is a rare congenital disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features. Recent reports have identified, in individuals with dominant Robinow syndrome, a specific type of variant characterized by being uniformly located in the penultimate exon of DVL1 and resulting in a −1 frameshift allele with a premature termination codon that escapes nonsense-mediated decay. Here, we studied a cohort of individuals who had been clinically diagnosed with Robinow syndrome but who had not received a molecular diagnosis from variant studies of DVL1, WNT5A , and ROR2 . Because of the uniform location of frameshift variants in DVL1 -mediated Robinow syndrome and the functional redundancy of DVL1, DVL2, and DVL3, we elected to pursue direct Sanger sequencing of the penultimate exon of DVL1 and its paralogs DVL2 and DVL3 to search for potential disease-associated variants. Remarkably, targeted sequencing identified five unrelated individuals harboring heterozygous, de novo frameshift variants in DVL3 , including two splice acceptor mutations and three 1 bp deletions. Similar to the variants observed in DVL1- mediated Robinow syndrome, all variants in DVL3 result in a −1 frameshift, indicating that these highly specific alterations might be a common cause of dominant Robinow syndrome. Here, we review the current knowledge of these peculiar variant alleles in DVL1 - and DVL3 -mediated Robinow syndrome and further elucidate the phenotypic features present in subjects with DVL1 and DVL3 frameshift mutations. [ABSTRACT FROM AUTHOR]

The paradigm of a single gene associated with one specific phenotype and mode of inheritance has been repeatedly challenged. Genotype-phenotype correlations can often be traced to different mutation types, localization of the variants in distinct protein domains, or the trigger of or escape from nonsense-mediated decay. Using whole-exome sequencing, we identified homozygous variants in EMC1 that segregated with a phenotype of developmental delay, hypotonia, scoliosis, and cerebellar atrophy in three families. In addition, a de novo heterozygous EMC1 variant was seen in an individual with a similar clinical and MRI imaging phenotype. EMC1 encodes a member of the endoplasmic reticulum (ER)-membrane protein complex (EMC), an evolutionarily conserved complex that has been proposed to have multiple roles in ER-associated degradation, ER-mitochondria tethering, and proper assembly of multi-pass transmembrane proteins. Perturbations of protein folding and organelle crosstalk have been implicated in neurodegenerative processes including cerebellar atrophy. We propose EMC1 as a gene in which either biallelic or monoallelic variants might lead to a syndrome including intellectual disability and preferential degeneration of the cerebellum. [ABSTRACT FROM AUTHOR]

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

Summary Development of the human nervous system involves complex interactions among fundamental cellular processes and requires a multitude of genes, many of which remain to be associated with human disease. We applied whole exome sequencing to 128 mostly consanguineous families with neurogenetic disorders that often included brain malformations. Rare variant analyses for both single nucleotide variant (SNV) and copy number variant (CNV) alleles allowed for identification of 45 novel variants in 43 known disease genes, 41 candidate genes, and CNVs in 10 families, with an overall potential molecular cause identified in >85% of families studied. Among the candidate genes identified, we found PRUNE , VARS , and DHX37 in multiple families and homozygous loss-of-function variants in AGBL2 , SLC18A2 , SMARCA1 , UBQLN1 , and CPLX1 . Neuroimaging and in silico analysis of functional and expression proximity between candidate and known disease genes allowed for further understanding of genetic networks underlying specific types of brain malformations. Video Abstract [ABSTRACT FROM AUTHOR]

Robinow syndrome is a genetically heterogeneous disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features and for which both autosomal-recessive and autosomal-dominant inheritance patterns have been described. Causative variants in the non-canonical signaling gene WNT5A underlie a subset of autosomal-dominant Robinow syndrome (DRS) cases, but most individuals with DRS remain without a molecular diagnosis. We performed whole-exome sequencing in four unrelated DRS-affected individuals without coding mutations in WNT5A and found heterozygous DVL1 exon 14 mutations in three of them. Targeted Sanger sequencing in additional subjects with DRS uncovered DVL1 exon 14 mutations in five individuals, including a pair of monozygotic twins. In total, six distinct frameshift mutations were found in eight subjects, and all were heterozygous truncating variants within the penultimate exon of DVL1 . In five families in which samples from unaffected parents were available, the variants were demonstrated to represent de novo mutations. All variant alleles are predicted to result in a premature termination codon within the last exon, escape nonsense-mediated decay (NMD), and most likely generate a C-terminally truncated protein with a distinct −1 reading-frame terminus. Study of the transcripts extracted from affected subjects’ leukocytes confirmed expression of both wild-type and variant alleles, supporting the hypothesis that mutant mRNA escapes NMD. Genomic variants identified in our study suggest that truncation of the C-terminal domain of DVL1, a protein hypothesized to have a downstream role in the Wnt-5a non-canonical pathway, is a common cause of DRS. [ABSTRACT FROM AUTHOR]

Human phosphoglucomutase 3 (PGM3) catalyzes the conversion of N-acetyl-glucosamine (GlcNAc)-6-phosphate into GlcNAc-1-phosphate during the synthesis of uridine diphosphate (UDP)-GlcNAc, a sugar nucleotide critical to multiple glycosylation pathways. We identified three unrelated children with recurrent infections, congenital leukopenia including neutropenia, B and T cell lymphopenia, and progression to bone marrow failure. Whole-exome sequencing demonstrated deleterious mutations in PGM3 in all three subjects, delineating their disease to be due to an unsuspected congenital disorder of glycosylation (CDG). Functional studies of the disease-associated PGM3 variants in E. coli cells demonstrated reduced PGM3 activity for all mutants tested. Two of the three children had skeletal anomalies resembling Desbuquois dysplasia: short stature, brachydactyly, dysmorphic facial features, and intellectual disability. However, these additional features were absent in the third child, showing the clinical variability of the disease. Two children received hematopoietic stem cell transplantation of cord blood and bone marrow from matched related donors; both had successful engraftment and correction of neutropenia and lymphopenia. We define PGM3-CDG as a treatable immunodeficiency, document the power of whole-exome sequencing in gene discoveries for rare disorders, and illustrate the utility of genomic analyses in studying combined and variable phenotypes. [ABSTRACT FROM AUTHOR]

Summary: CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis. [Copyright &y& Elsevier]

Summary DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes. Graphical Abstract Highlights • Orthogonal DNA sequencing approaches are required to observe all variant types • De novo SNVs and indels accompany non-recurrent structural variants (SVs) in cis • SV-associated SNVs primarily occur within genes and over megabase-sized distances • MMBIR involves extensive DNA replication resulting in regional hypermutation Newly occurring structural variants within in human genomes spawn extensive, local single-nucleotide changes leading to an enhanced mutational burden within proximal genes. [ABSTRACT FROM AUTHOR]

Summary De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the “CNV-mutator state,” and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology. [ABSTRACT FROM AUTHOR]

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9 -/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis., Competing Interests: Declaration of interests J.R.L. has stock ownership in 23andMe, is a paid consultant for Genome International, and is a co-inventor on multiple US and European patents related to molecular diagnostics for inherited neuropathies, eye diseases, genomic disorders, and bacterial genomic fingerprinting. The Department of Molecular and Human Genetics at the Baylor College of Medicine receives revenue from clinical genetic testing conducted at Baylor Genetics (BG) Laboratories. F.M. and T.S.-S. are employees of GeneDx., (Copyright © 2023 American Society of Human Genetics. All rights reserved.)

The leucine-rich glioma-inactivated (LGI) family consists of four highly conserved paralogous genes, LGI1-4, that are highly expressed in mammalian central and/or peripheral nervous systems. LGI1 antibodies are detected in subjects with autoimmune limbic encephalitis and peripheral nerve hyperexcitability syndromes (PNHSs) such as Isaacs and Morvan syndromes. Pathogenic variations of LGI1 and LGI4 are associated with neurological disorders as disease traits including familial temporal lobe epilepsy and neurogenic arthrogryposis multiplex congenita 1 with myelin defects, respectively. No human disease has been reported associated with either LGI2 or LGI3. We implemented exome sequencing and family-based genomics to identify individuals with deleterious variants in LGI3 and utilized GeneMatcher to connect practitioners and researchers worldwide to investigate the clinical and electrophysiological phenotype in affected subjects. We also generated Lgi3-null mice and performed peripheral nerve dissection and immunohistochemistry to examine the juxtaparanode LGI3 microarchitecture. As a result, we identified 16 individuals from eight unrelated families with loss-of-function (LoF) bi-allelic variants in LGI3. Deep phenotypic characterization showed LGI3 LoF causes a potentially clinically recognizable PNHS trait characterized by global developmental delay, intellectual disability, distal deformities with diminished reflexes, visible facial myokymia, and distinctive electromyographic features suggestive of motor nerve instability. Lgi3-null mice showed reduced and mis-localized Kv1 channel complexes in myelinated peripheral axons. Our data demonstrate bi-allelic LoF variants in LGI3 cause a clinically distinguishable disease trait of PNHS, most likely caused by disturbed Kv1 channel distribution in the absence of LGI3., Competing Interests: Declaration of interests J.R.L. has stock ownership in 23andMe; is a paid consultant for Regeneron Genetics Center; and is a co-inventor on multiple United States and European patents related to molecular diagnostics for inherited neuropathies, eye diseases, genomic disorders, and bacterial genomic fingerprinting. The Department of Molecular and Human Genetics at Baylor College of Medicine receives revenue from clinical genetic testing conducted at Baylor Genetics (BG); J.R.L. serves on the Scientific Advisory Board (SAB) of BG., (Copyright © 2022. Published by Elsevier Inc.)

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey., Competing Interests: Declaration of interests J.R.L. has stock ownership in 23andMe, is a paid consultant for the Regeneron Genetics Center, and is a co-inventor on multiple United States and European patents related to molecular diagnostics for inherited neuropathies, eye diseases, and bacterial genomic fingerprinting. The Department of Molecular and Human Genetics at Baylor College of Medicine receives revenue from clinical genetic testing conducted at Baylor Genetics (BG) Laboratories. J.R.L. serves on the Scientific Advisory Board of BG. Other authors have no potential conflicts to report., (Copyright © 2021 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Coatomer complexes function in the sorting and trafficking of proteins between subcellular organelles. Pathogenic variants in coatomer subunits or associated factors have been reported in multi-systemic disorders, i.e., coatopathies, that can affect the skeletal and central nervous systems. We have identified loss-of-function variants in COPB2, a component of the coatomer complex I (COPI), in individuals presenting with osteoporosis, fractures, and developmental delay of variable severity. Electron microscopy of COPB2-deficient subjects' fibroblasts showed dilated endoplasmic reticulum (ER) with granular material, prominent rough ER, and vacuoles, consistent with an intracellular trafficking defect. We studied the effect of COPB2 deficiency on collagen trafficking because of the critical role of collagen secretion in bone biology. COPB2 siRNA-treated fibroblasts showed delayed collagen secretion with retention of type I collagen in the ER and Golgi and altered distribution of Golgi markers. copb2-null zebrafish embryos showed retention of type II collagen, disorganization of the ER and Golgi, and early larval lethality. Copb2 +/- mice exhibited low bone mass, and consistent with the findings in human cells and zebrafish, studies in Copb2 +/- mouse fibroblasts suggest ER stress and a Golgi defect. Interestingly, ascorbic acid treatment partially rescued the zebrafish developmental phenotype and the cellular phenotype in Copb2 +/- mouse fibroblasts. This work identifies a form of coatopathy due to COPB2 haploinsufficiency, explores a potential therapeutic approach for this disorder, and highlights the role of the COPI complex as a regulator of skeletal homeostasis., Competing Interests: Declaration of interests The Department of Molecular and Human Genetics at Baylor College of Medicine receives financial support from Baylor Genetics. Dr. Brendan Lee serves on the Board of Directors of Baylor Genetics and chairs its Scientific Advisory Board but receives no personal income from these positions., (Copyright © 2021 American Society of Human Genetics. All rights reserved.)

NTNG2 encodes netrin-G2, a membrane-anchored protein implicated in the molecular organization of neuronal circuitry and synaptic organization and diversification in vertebrates. In this study, through a combination of exome sequencing and autozygosity mapping, we have identified 16 individuals (from seven unrelated families) with ultra-rare homozygous missense variants in NTNG2; these individuals present with shared features of a neurodevelopmental disorder consisting of global developmental delay, severe to profound intellectual disability, muscle weakness and abnormal tone, autistic features, behavioral abnormalities, and variable dysmorphisms. The variants disrupt highly conserved residues across the protein. Functional experiments, including in silico analysis of the protein structure, in vitro assessment of cell surface expression, and in vitro knockdown, revealed potential mechanisms of pathogenicity of the variants, including loss of protein function and decreased neurite outgrowth. Our data indicate that appropriate expression of NTNG2 plays an important role in neurotypical development., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Arthrogryposis is a clinical finding that is present either as a feature of a neuromuscular condition or as part of a systemic disease in over 400 Mendelian conditions. The underlying molecular etiology remains largely unknown because of genetic and phenotypic heterogeneity. We applied exome sequencing (ES) in a cohort of 89 families with the clinical sign of arthrogryposis. Additional molecular techniques including array comparative genomic hybridization (aCGH) and Droplet Digital PCR (ddPCR) were performed on individuals who were found to have pathogenic copy number variants (CNVs) and mosaicism, respectively. A molecular diagnosis was established in 65.2% (58/89) of families. Eleven out of 58 families (19.0%) showed evidence for potential involvement of pathogenic variation at more than one locus, probably driven by absence of heterozygosity (AOH) burden due to identity-by-descent (IBD). RYR3, MYOM2, ERGIC1, SPTBN4, and ABCA7 represent genes, identified in two or more families, for which mutations are probably causative for arthrogryposis. We also provide evidence for the involvement of CNVs in the etiology of arthrogryposis and for the idea that both mono-allelic and bi-allelic variants in the same gene cause either similar or distinct syndromes. We were able to identify the molecular etiology in nine out of 20 families who underwent reanalysis. In summary, our data from family-based ES further delineate the molecular etiology of arthrogryposis, yielded several candidate disease-associated genes, and provide evidence for mutational burden in a biological pathway or network. Our study also highlights the importance of reanalysis of individuals with unsolved diagnoses in conjunction with sequencing extended family members., (Copyright © 2019 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Acetylation of the lysine residues in histones and other DNA-binding proteins plays a major role in regulation of eukaryotic gene expression. This process is controlled by histone acetyltransferases (HATs/KATs) found in multiprotein complexes that are recruited to chromatin by the scaffolding subunit transformation/transcription domain-associated protein (TRRAP). TRRAP is evolutionarily conserved and is among the top five genes intolerant to missense variation. Through an international collaboration, 17 distinct de novo or apparently de novo variants were identified in TRRAP in 24 individuals. A strong genotype-phenotype correlation was observed with two distinct clinical spectra. The first is a complex, multi-systemic syndrome associated with various malformations of the brain, heart, kidneys, and genitourinary system and characterized by a wide range of intellectual functioning; a number of affected individuals have intellectual disability (ID) and markedly impaired basic life functions. Individuals with this phenotype had missense variants clustering around the c.3127G>A p.(Ala1043Thr) variant identified in five individuals. The second spectrum manifested with autism spectrum disorder (ASD) and/or ID and epilepsy. Facial dysmorphism was seen in both groups and included upslanted palpebral fissures, epicanthus, telecanthus, a wide nasal bridge and ridge, a broad and smooth philtrum, and a thin upper lip. RNA sequencing analysis of skin fibroblasts derived from affected individuals skin fibroblasts showed significant changes in the expression of several genes implicated in neuronal function and ion transport. Thus, we describe here the clinical spectrum associated with TRRAP pathogenic missense variants, and we suggest a genotype-phenotype correlation useful for clinical evaluation of the pathogenicity of the variants., (Copyright © 2019 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

SPONASTRIME dysplasia is an autosomal-recessive spondyloepimetaphyseal dysplasia characterized by spine (spondylar) abnormalities, midface hypoplasia with a depressed nasal bridge, metaphyseal striations, and disproportionate short stature. Scoliosis, coxa vara, childhood cataracts, short dental roots, and hypogammaglobulinemia have also been reported in this disorder. Although an autosomal-recessive inheritance pattern has been hypothesized, pathogenic variants in a specific gene have not been discovered in individuals with SPONASTRIME dysplasia. Here, we identified bi-allelic variants in TONSL, which encodes the Tonsoku-like DNA repair protein, in nine subjects (from eight families) with SPONASTRIME dysplasia, and four subjects (from three families) with short stature of varied severity and spondylometaphyseal dysplasia with or without immunologic and hematologic abnormalities, but no definitive metaphyseal striations at diagnosis. The finding of early embryonic lethality in a Tonsl -/- murine model and the discovery of reduced length, spinal abnormalities, reduced numbers of neutrophils, and early lethality in a tonsl -/- zebrafish model both support the hypomorphic nature of the identified TONSL variants. Moreover, functional studies revealed increased amounts of spontaneous replication fork stalling and chromosomal aberrations, as well as fewer camptothecin (CPT)-induced RAD51 foci in subject-derived cell lines. Importantly, these cellular defects were rescued upon re-expression of wild-type (WT) TONSL; this rescue is consistent with the hypothesis that hypomorphic TONSL variants are pathogenic. Overall, our studies in humans, mice, zebrafish, and subject-derived cell lines confirm that pathogenic variants in TONSL impair DNA replication and homologous recombination-dependent repair processes, and they lead to a spectrum of skeletal dysplasia phenotypes with numerous extra-skeletal manifestations., (Copyright © 2019 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Primary defects in lung branching morphogenesis, resulting in neonatal lethal pulmonary hypoplasias, are incompletely understood. To elucidate the pathogenetics of human lung development, we studied a unique collection of samples obtained from deceased individuals with clinically and histopathologically diagnosed interstitial neonatal lung disorders: acinar dysplasia (n = 14), congenital alveolar dysplasia (n = 2), and other lethal lung hypoplasias (n = 10). We identified rare heterozygous copy-number variant deletions or single-nucleotide variants (SNVs) involving TBX4 (n = 8 and n = 2, respectively) or FGF10 (n = 2 and n = 2, respectively) in 16/26 (61%) individuals. In addition to TBX4, the overlapping ∼2 Mb recurrent and nonrecurrent deletions at 17q23.1q23.2 identified in seven individuals with lung hypoplasia also remove a lung-specific enhancer region. Individuals with coding variants involving either TBX4 or FGF10 also harbored at least one non-coding SNV in the predicted lung-specific enhancer region, which was absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies. The occurrence of rare coding variants involving TBX4 or FGF10 with the putative hypomorphic non-coding SNVs implies a complex compound inheritance of these pulmonary hypoplasias. Moreover, they support the importance of TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling in human lung organogenesis and help to explain the histopathological continuum observed in these rare lethal developmental disorders of the lung., (Copyright © 2018 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Locus heterogeneity characterizes a variety of skeletal dysplasias often due to interacting or overlapping signaling pathways. Robinow syndrome is a skeletal disorder historically refractory to molecular diagnosis, potentially stemming from substantial genetic heterogeneity. All current known pathogenic variants reside in genes within the noncanonical Wnt signaling pathway including ROR2, WNT5A, and more recently, DVL1 and DVL3. However, ∼70% of autosomal-dominant Robinow syndrome cases remain molecularly unsolved. To investigate this missing heritability, we recruited 21 families with at least one family member clinically diagnosed with Robinow or Robinow-like phenotypes and performed genetic and genomic studies. In total, four families with variants in FZD2 were identified as well as three individuals from two families with biallelic variants in NXN that co-segregate with the phenotype. Importantly, both FZD2 and NXN are relevant protein partners in the WNT5A interactome, supporting their role in skeletal development. In addition to confirming that clustered -1 frameshifting variants in DVL1 and DVL3 are the main contributors to dominant Robinow syndrome, we also found likely pathogenic variants in candidate genes GPC4 and RAC3, both linked to the Wnt signaling pathway. These data support an initial hypothesis that Robinow syndrome results from perturbation of the Wnt/PCP pathway, suggest specific relevant domains of the proteins involved, and reveal key contributors in this signaling cascade during human embryonic development. Contrary to the view that non-allelic genetic heterogeneity hampers gene discovery, this study demonstrates the utility of rare disease genomic studies to parse gene function in human developmental pathways., (Copyright © 2017 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

The underlying genetic etiology of rhabdomyolysis remains elusive in a significant fraction of individuals presenting with recurrent metabolic crises and muscle weakness. Using exome sequencing, we identified bi-allelic mutations in TANGO2 encoding transport and Golgi organization 2 homolog (Drosophila) in 12 subjects with episodic rhabdomyolysis, hypoglycemia, hyperammonemia, and susceptibility to life-threatening cardiac tachyarrhythmias. A recurrent homozygous c.460G>A (p.Gly154Arg) mutation was found in four unrelated individuals of Hispanic/Latino origin, and a homozygous ∼34 kb deletion affecting exons 3-9 was observed in two families of European ancestry. One individual of mixed Hispanic/European descent was found to be compound heterozygous for c.460G>A (p.Gly154Arg) and the deletion of exons 3-9. Additionally, a homozygous exons 4-6 deletion was identified in a consanguineous Middle Eastern Arab family. No homozygotes have been reported for these changes in control databases. Fibroblasts derived from a subject with the recurrent c.460G>A (p.Gly154Arg) mutation showed evidence of increased endoplasmic reticulum stress and a reduction in Golgi volume density in comparison to control. Our results show that the c.460G>A (p.Gly154Arg) mutation and the exons 3-9 heterozygous deletion in TANGO2 are recurrent pathogenic alleles present in the Latino/Hispanic and European populations, respectively, causing considerable morbidity in the homozygotes in these populations., (Copyright © 2016 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Charcot-Marie-Tooth (CMT) disease is a clinically and genetically heterogeneous distal symmetric polyneuropathy. Whole-exome sequencing (WES) of 40 individuals from 37 unrelated families with CMT-like peripheral neuropathy refractory to molecular diagnosis identified apparent causal mutations in ∼ 45% (17/37) of families. Three candidate disease genes are proposed, supported by a combination of genetic and in vivo studies. Aggregate analysis of mutation data revealed a significantly increased number of rare variants across 58 neuropathy-associated genes in subjects versus controls, confirmed in a second ethnically discrete neuropathy cohort, suggesting that mutation burden potentially contributes to phenotypic variability. Neuropathy genes shown to have highly penetrant Mendelizing variants (HPMVs) and implicated by burden in families were shown to interact genetically in a zebrafish assay exacerbating the phenotype established by the suppression of single genes. Our findings suggest that the combinatorial effect of rare variants contributes to disease burden and variable expressivity., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

Gall-forming arthropods are highly specialized herbivores that, in combination with their hosts, produce extended phenotypes with unique morphologies [1]. Many are economically important, and others have improved our understanding of ecology and adaptive radiation [2]. However, the mechanisms that these arthropods use to induce plant galls are poorly understood. We sequenced the genome of the Hessian fly (Mayetiola destructor; Diptera: Cecidomyiidae), a plant parasitic gall midge and a pest of wheat (Triticum spp.), with the aim of identifying genic modifications that contribute to its plant-parasitic lifestyle. Among several adaptive modifications, we discovered an expansive reservoir of potential effector proteins. Nearly 5% of the 20,163 predicted gene models matched putative effector gene transcripts present in the M. destructor larval salivary gland. Another 466 putative effectors were discovered among the genes that have no sequence similarities in other organisms. The largest known arthropod gene family (family SSGP-71) was also discovered within the effector reservoir. SSGP-71 proteins lack sequence homologies to other proteins, but their structures resemble both ubiquitin E3 ligases in plants and E3-ligase-mimicking effectors in plant pathogenic bacteria. SSGP-71 proteins and wheat Skp proteins interact in vivo. Mutations in different SSGP-71 genes avoid the effector-triggered immunity that is directed by the wheat resistance genes H6 and H9. Results point to effectors as the agents responsible for arthropod-induced plant gall formation., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Invertebrate model systems are powerful tools for studying human disease owing to their genetic tractability and ease of screening. We conducted a mosaic genetic screen of lethal mutations on the Drosophila X chromosome to identify genes required for the development, function, and maintenance of the nervous system. We identified 165 genes, most of whose function has not been studied in vivo. In parallel, we investigated rare variant alleles in 1,929 human exomes from families with unsolved Mendelian disease. Genes that are essential in flies and have multiple human homologs were found to be likely to be associated with human diseases. Merging the human data sets with the fly genes allowed us to identify disease-associated mutations in six families and to provide insights into microcephaly associated with brain dysgenesis. This bidirectional synergism between fly genetics and human genomics facilitates the functional annotation of evolutionarily conserved genes involved in human health., (Copyright © 2014 Elsevier Inc. All rights reserved.)