New tools for Mendelian disease gene identification: PhenoDB variant analysis module; and GeneMatcher, a web-based tool for linking investigators with an interest in the same gene

Exome sequencing (ES) of a pediatric cohort with chronic endocrine diseases: a single-center study (within the framework of the TRANSLATE-NAMSE project)

BACKGROUND

Endocrine disorders are heterogeneous and include a significant number of rare monogenic diseases.

METHODS

We performed exome sequencing (ES) in 106 children recruited from a single center within the TRANSLATE‑NAMSE project. They were categorized into subgroups: proportionate short stature (PSS), disproportionate short stature (DSS), hypopituitarism (H), differences in sexual development (DSD), syndromic diseases (SD) and others.

RESULTS

The overall diagnostic yield was 34.9% (n = 37/106), including 5 patients with variants in candidate genes, which have contributed to collaborations to identify gene-disease associations. The diagnostic yield varied significantly between subgroups: PSS: 16.6% (1/6); DSS: 18.8% (3/16); H: 17.1% (6/35); DSD: 37.5% (3/8); SD: 66.6% (22/33); others: 25% (2/8). Confirmed diagnoses included 75% ultrarare diseases. Three patients harbored more than one disease-causing variant, resulting in dual diagnoses.

CONCLUSIONS

ES is an effective tool for genetic diagnosis in pediatric patients with complex endocrine diseases. An accurate phenotypic description, including comprehensive endocrine diagnostics, as well as the evaluation of variants in multidisciplinary case conferences involving geneticists, are necessary for personalized diagnostic care. Here, we illustrate the broad spectrum of genetic endocrinopathies that have led to the initiation of specific treatment, surveillance, and family counseling.

Rare exonic CELSR3 variants identified in Bladder Exstrophy Epispadias Complex

Introduction/background

Bladder exstrophy epispadias complex (BEEC) is a rare congenital anomaly of unknown etiology, although, genetic and environmental factors have been associated with its development. Variants in several genes expressed in the urogenital pathway have been reported as causative for bladder exstrophy in human and murine models. The expansion of next-generation sequencing and molecular genomics has improved our ability to identify the underlying genetic causes of similarly complex diseases and could thus assist with the investigation of the molecular basis of BEEC.

Objective

The objective was to identify the presence of rare heterozygous variants in genes previously implicated in bladder exstrophy and correlate them with the presence or absence of bladder regeneration in our study population.

Patients and Methods

We present a case series of 12 patients with BEEC who had bladder biopsies performed by pediatric urology during bladder neck reconstruction or bladder augmentation. Cases were classified as "sufficient" or "insufficient" (n = 5 and 7, respectively) based on a bladder volume of greater than or less than 40% of expected bladder size. Control bladder tissue specimens were obtained from patients (n = 6) undergoing biopsies for conditions other than bladder exstrophy. Whole exome sequencing was performed on DNA isolated from the bladder specimens. Based on the hypothesis of de novo mutations, as well as the potential implications of autosomal dominant conditions with incomplete penetrance, each case was evaluated for autosomal dominant variants in a set of genes previously implicated in BEEC.

Results

Our review of the literature identified 44 genes that have been implicated in human models of bladder exstrophy. Our whole exome sequencing data analysis identified rare variants in two of these genes among the cases classified as sufficient, and seven variants in five of these genes among the cases classified as insufficient.

Conclusion

We identified rare variants in seven previously implicated genes in our BEEC specimens. Additional research is needed to further understand the cellular signaling underlying this potentially genetically heterogeneous embryological condition.

Whole genome sequencing as a first-tier diagnostic test for infants in neonatal intensive care units: A pilot study in Brazil

In this pilot study, we aimed to evaluate the feasibility of whole genome sequencing (WGS) as a first-tier diagnostic test for infants hospitalized in neonatal intensive care units in the Brazilian healthcare system. The cohort presented here results from a joint collaboration between private and public hospitals in Brazil considering the initiative of a clinical laboratory to provide timely diagnosis for critically ill infants. We performed trio (proband and parents) WGS in 21 infants suspected of a genetic disease with an urgent need for diagnosis to guide medical care. Overall, the primary indication for genetic testing was dysmorphic syndromes (n = 14, 67%) followed by inborn errors of metabolism (n = 6, 29%) and skeletal dysplasias (n = 1, 5%). The diagnostic yield in our cohort was 57% (12/21) based on cases that received a definitive or likely definitive diagnostic result from WGS analysis. A total of 16 pathogenic/likely pathogenic variants and 10 variants of unknown significance were detected, and in most cases inherited from an unaffected parent. In addition, the reported variants were of different types, but mainly missense (58%) and associated with autosomal diseases (19/26); only three were associated with X-linked diseases, detected in hemizygosity in the proband an inherited from an unaffected mother. Notably, we identified 10 novel variants, absent from public genomic databases, in our cohort. Considering the entire diagnostic process, the average turnaround time from enrollment to medical report in our study was 53 days. Our findings demonstrate the remarkable utility of WGS as a diagnostic tool, elevating the potential of transformative impact since it outperforms conventional genetic tests. Here, we address the main challenges associated with implementing WGS in the medical care system in Brazil, as well as discuss the potential benefits and limitations of WGS as a diagnostic tool in the neonatal care setting.

Loss of TBC1D2B causes a progressive neurological disorder with gingival overgrowth

Biallelic loss-of-function variants in TBC1D2B have been reported in five subjects with cognitive impairment and seizures with or without gingival overgrowth. TBC1D2B belongs to the family of Tre2-Bub2-Cdc16 (TBC)-domain containing RAB-specific GTPase activating proteins (TBC/RABGAPs). Here, we report five new subjects with biallelic TBC1D2B variants, including two siblings, and delineate the molecular and clinical features in the ten subjects known to date. One of the newly reported subjects was compound heterozygous for the TBC1D2B variants c.2584C>T; p.(Arg862Cys) and c.2758C>T; p.(Arg920*). In subject-derived fibroblasts, TBC1D2B mRNA level was similar to control cells, while the TBC1D2B protein amount was reduced by about half. In one of two siblings with a novel c.360+1G>T splice site variant, TBC1D2B transcript analysis revealed aberrantly spliced mRNAs and a drastically reduced TBC1D2B mRNA level in leukocytes. The molecular spectrum included 12 different TBC1D2B variants: seven nonsense, three frameshifts, one splice site, and one missense variant. Out of ten subjects, three had fibrous dysplasia of the mandible, two of which were diagnosed as cherubism. Most subjects developed gingival overgrowth. Half of the subjects had developmental delay. Seizures occurred in 80% of the subjects. Six subjects showed a progressive disease with mental deterioration. Brain imaging revealed cerebral and/or cerebellar atrophy with or without lateral ventricle dilatation. The TBC1D2B disorder is a progressive neurological disease with gingival overgrowth and abnormal mandible morphology. As TBC1D2B has been shown to positively regulate autophagy, defects in autophagy and the endolysosomal system could be associated with neuronal dysfunction and the neurodegenerative disease in the affected individuals.

Coordination among frequent genetic variants imparts substance use susceptibility and pathogenesis

Determining the key genetic variants is a crucial step to comprehensively understand substance use disorders (SUDs). In this study, utilizing whole exome sequences of five multi-generational pedigrees with SUDs, we used an integrative omics-based approach to uncover candidate genetic variants that impart susceptibility to SUDs and influence addition traits. We identified several SNPs and rare, protein-function altering variants in genes, GRIA3, NCOR1, and SHANK1; compound heterozygous variants in LNPEP, LRP1, and TBX2, that play a significant role in the neurotransmitter-neuropeptide axis, specifically in the dopaminergic circuits. We also noted a greater frequency of heterozygous and recessive variants in genes involved in the structural and functional integrity of synapse receptors, CHRNA4, CNR2, GABBR1, DRD4, NPAS4, ADH1B, ADH1C, OPRM1, and GABBR2. Variant analysis in upstream promoter regions revealed regulatory variants in NEK9, PRRX1, PRPF4B, CELA2A, RABGEF1, and CRBN, crucial for dopamine regulation. Using family-and pedigree-based data, we identified heterozygous recessive alleles in LNPEP, LRP1 (4 frameshift deletions), and TBX2 (2 frameshift deletions) linked to SUDs. GWAS overlap identified several SNPs associated with SUD susceptibility, including rs324420 and rs1229984. Furthermore, miRNA variant analysis revealed notable variants in mir-548 U and mir-532. Pathway studies identified the presence of extensive coordination among these genetic variants to impart substance use susceptibility and pathogenesis. This study identified variants that were found to be overrepresented among genes of dopaminergic circuits participating in the neurotransmitter-neuropeptide axis, suggesting pleiotropic influences in the development and sustenance of chronic substance use. The presence of a diverse set of haploinsufficient variants in varying frequencies demonstrates the existence of extraordinary coordination among them in attributing risk and modulating severity to SUDs.

De novo variants identified by trio whole exome sequencing of bladder exstrophy epispadias complex

Bladder exstrophy epispadias complex (BEEC) encompasses a spectrum of conditions ranging from mild epispadias to the most severe form: omphalocele-bladder exstrophy-imperforate anus-spinal defects (OEIS). BEEC involves abnormalities related to anatomical structures that are proposed to have a similar underlying etiology and pathogenesis. In general, BEEC, is considered to arise from a sequence of events in embryonic development and is believed to be a multi-etiological disease with contributions from genetic and environmental factors. Several genes have been implicated and mouse models have been generated, including a knockout model of p63, which is involved in the synthesis of stratified epithelium. Mice lacking p63 have undifferentiated ventral urothelium. MNX1 has also been implicated. In addition, cigarette smoking, diazepam and clomid have been implied as environmental factors due to their relative association. By in large, the etiology and pathogenesis of human BEEC is unknown. We performed de novo analysis of whole exome sequencing (WES) of germline samples from 31 unrelated trios where the probands have a diagnosis of BEEC syndrome. We also evaluated the DECIPHER database to identify copy number variants (CNVs) in genes in individuals with the search terms "bladder exstrophy" in an attempt to identify additional candidate genes within these regions. Several de novo variants were identified; however, a candidate gene is still unclear. This data further supports the multi-etiological nature of BEEC.

Anxiety and dysautonomia symptoms in patients with a NaV1.7 mutation and the potential benefits of low-dose short-acting guanfacine

PURPOSE

Guanfacine is an α2A-adrenergic receptor agonist, FDA-approved to treat attention-deficit hyperactivity disorder and high blood pressure, typically as an extended-release formulation up to 7 mg/day. In our dysautonomia clinic, we observed that off-label use of short-acting guanfacine at 1 mg/day facilitated symptom relief in two families with multiple members presenting with severe generalized anxiety. We also noted anecdotal improvements in associated dysautonomia symptoms such as hyperhidrosis, cognitive impairment, and palpitations. We postulated that a genetic deficit existed in these patients that might augment guanfacine susceptibility.

METHODS

We used whole-exome sequencing to identify mutations in patients with shared generalized anxiety and dysautonomia symptoms. Guanfacine-induced changes in the function of voltage-gated Na+ channels were investigated using voltage-clamp electrophysiology.

RESULTS

Whole-exome sequencing uncovered the p.I739V mutation in SCN9A in the proband of two nonrelated families. Moreover, guanfacine inhibited ionic currents evoked by wild-type and mutant NaV1.7 encoded by SCN9A, as well as other NaV channel subtypes to a varying degree.

CONCLUSION

Our study provides further evidence for a possible pathophysiological role of NaV1.7 in anxiety and dysautonomia. Combined with off-target effects on NaV channel function, daily administration of 1 mg short-acting guanfacine may be sufficient to normalize NaV channel mutation-induced changes in sympathetic activity, perhaps aided by partial inhibition of NaV1.7 or other channel subtypes. In a broader context, expanding genetic and functional data about ion channel aberrations may enable the prospect of stratifying patients in which mutation-induced increased sympathetic tone normalization by guanfacine can support treatment strategies for anxiety and dysautonomia symptoms.

Molecular and Phenotypic Characterization of the RORB-Related Disorder

BACKGROUND AND OBJECTIVES

Heterozygous variants in RAR-related orphan receptor B (RORB) have recently been associated with susceptibility to idiopathic generalized epilepsy. However, few reports have been published so far describing pathogenic variants of this gene in patients with epilepsy and intellectual disability (ID). In this study, we aimed to delineate the epilepsy phenotype associated with RORB pathogenic variants and to provide arguments in favor of the pathogenicity of variants.

METHODS

Through an international collaboration, we analyzed seizure characteristics, EEG data, and genotypes of a cohort of patients with heterozygous variants in RORB. To gain insight into disease mechanisms, we performed ex vivo cortical electroporation in mouse embryos of 5 selected variants, 2 truncating and 3 missense, and evaluated on expression and quantified changes in axonal morphology.

RESULTS

We identified 35 patients (17 male, median age 10 years, range 2.5-23 years) carrying 32 different heterozygous variants in RORB, including 28 single-nucleotide variants or small insertions/deletions (12 missense, 12 frameshift or nonsense, 2 splice-site variants, and 2 in-frame deletions), and 4 microdeletions; de novo in 18 patients and inherited in 10. Seizures were reported in 31/35 (89%) patients, with a median age at onset of 3 years (range 4 months-12 years). Absence seizures occurred in 25 patients with epilepsy (81%). Nineteen patients experienced a single seizure type: absences, myoclonic absences, or absences with eyelid myoclonia and focal seizures. Nine patients had absence seizures combined with other generalized seizure types. One patient had presented with absences associated with photosensitive occipital seizures. Three other patients had generalized tonic-clonic seizures without absences. ID of variable degree was observed in 85% of the patients. Expression studies in cultured neurons showed shorter axons for the 5 tested variants, both truncating and missense variants, supporting an impaired protein function.

DISCUSSION

In most patients, the phenotype of the RORB-related disorder associates absence seizures with mild-to-moderate ID. In silico and in vitro evaluation of the variants in our cohort, including axonal morphogenetic experiments in cultured neurons, supports their pathogenicity, showing a hypomorphic effect.

The best of both worlds: Blending cutting-edge research with clinical processes for a productive exome clinic

Genomic medicine has been transformed by next-generation sequencing (NGS), inclusive of exome sequencing (ES) and genome sequencing (GS). Currently, ES is offered widely in clinical settings, with a less prevalent alternative model consisting of hybrid programs that incorporate research ES along with clinical patient workflows. We were among the earliest to implement a hybrid ES clinic, have provided diagnoses to 45% of probands, and have identified several novel candidate genes. Our program is enabled by a cost-effective investment by the health system and is unique in encompassing all the processes that have been variably included in other hybrid/clinical programs. These include careful patient selection, utilization of a phenotype-agnostic bioinformatics pipeline followed by manual curation of variants and phenotype integration by clinicians, close collaborations between the clinicians and the bioinformatician, pursuit of interesting variants, communication of results to patients in categories that are predicated upon the certainty of a diagnosis, and tracking changes in results over time and the underlying mechanisms for such changes. Due to its effectiveness, scalability to GS and its resource efficiency, specific elements of our paradigm can be incorporated into existing clinical settings, or the entire hybrid model can be implemented within health systems that have genomic medicine programs, to provide NGS in a scientifically rigorous, yet pragmatic setting.

De novo PHF5A variants are associated with craniofacial abnormalities, developmental delay, and hypospadias

PURPOSE

The SF3B splicing complex is composed of SF3B1-6 and PHF5A. We report a developmental disorder caused by de novo variants in PHF5A.

METHODS

Clinical, genomic, and functional studies using subject-derived fibroblasts and a heterologous cellular system were performed.

RESULTS

We studied 9 subjects with congenital malformations, including preauricular tags and hypospadias, growth abnormalities, and developmental delay who had de novo heterozygous PHF5A variants, including 4 loss-of-function (LOF), 3 missense, 1 splice, and 1 start-loss variant. In subject-derived fibroblasts with PHF5A LOF variants, wild-type and variant PHF5A mRNAs had a 1:1 ratio, and PHF5A mRNA levels were normal. Transcriptome sequencing revealed alternative promoter use and downregulated genes involved in cell-cycle regulation. Subject and control fibroblasts had similar amounts of PHF5A with the predicted wild-type molecular weight and of SF3B1-3 and SF3B6. SF3B complex formation was unaffected in 2 subject cell lines.

CONCLUSION

Our data suggest the existence of feedback mechanisms in fibroblasts with PHF5A LOF variants to maintain normal levels of SF3B components. These compensatory mechanisms in subject fibroblasts with PHF5A or SF3B4 LOF variants suggest disturbed autoregulation of mutated splicing factor genes in specific cell types, that is, neural crest cells, during embryonic development rather than haploinsufficiency as pathomechanism.

COL11A2 as a candidate gene for vertebral malformations and congenital scoliosis

Human vertebral malformations (VMs) have an estimated incidence of 1/2000 and are associated with significant health problems including congenital scoliosis (CS) and recurrent organ system malformation syndromes such as VACTERL (vertebral anomalies; anal abnormalities; cardiac abnormalities; tracheo-esophageal fistula; renal anomalies; limb anomalies). The genetic cause for the vast majority of VMs are unknown. In a CS/VM patient cohort, three COL11A2 variants (R130W, R1407L and R1413H) were identified in two patients with cervical VM. A third patient with a T9 hemivertebra and the R130W variant was identified from a separate study. These substitutions are predicted to be damaging to protein function, and R130 and R1407 residues are conserved in zebrafish Col11a2. To determine the role for COL11A2 in vertebral development, CRISPR/Cas9 was used to create a nonsense mutation (col11a2L642*) as well as a full gene locus deletion (col11a2del) in zebrafish. Both col11a2L642*/L642* and col11a2del/del mutant zebrafish exhibit vertebral fusions in the caudal spine, which form due to mineralization across intervertebral segments. To determine the functional consequence of VM-associated variants, we assayed their ability to suppress col11a2del VM phenotypes following transgenic expression within the developing spine. While wildtype col11a2 expression suppresses fusions in col11a2del/+ and col11a2del/del backgrounds, patient missense variant-bearing col11a2 failed to rescue the loss-of-function phenotype in these animals. These results highlight an essential role for COL11A2 in vertebral development and support a pathogenic role for two missense variants in CS.

A holistic approach to maximise diagnostic output in trio exome sequencing

Introduction

Rare genetic diseases are a major cause for severe illness in children. Whole exome sequencing (WES) is a powerful tool for identifying genetic causes of rare diseases. For a better and faster assessment of the vast number of variants that are identified in the index patient in WES, parental sequencing can be applied ("trio WES").

Methods

We assessed the diagnostic rate of routine trio WES including analysis of copy number variants in 224 pediatric patients during an evaluation period of three years.

Results

Trio WES provided a diagnosis in 67 (30%) of all 224 analysed children. The turnaround time of trio WES analysis has been reduced significantly from 41 days in 2019 to 23 days in 2021. Copy number variants could be identified to be causative in 10 cases (4.5%), underlying the importance of copy number variant analysis. Variants in three genes which were previously not associated with a clinical condition (GAD1, TMEM222 and ZNFX1) were identified using the matching tool GeneMatcher and were part of the first description of a new syndrome.

Discussion

Trio WES has proven to have a high diagnostic yield and to shorten the process of identifying the correct diagnosis in paediatric patients. Re-evaluation of all 224 trio WES 1-3 years after initial analysis did not establish new diagnoses. Initiating (trio) WES as a first-tier diagnostics including copy number variant detection should be considered as early as possible, especially for children treated in ICU, if a monogenetic disease is suspected.

Homozygous loss-of-function variants in FILIP1 cause autosomal recessive arthrogryposis multiplex congenita with microcephaly

Arthrogryposis multiplex congenita forms a broad group of clinically and etiologically heterogeneous disorders characterized by congenital joint contractures that involve at least two different parts of the body. Neurological and muscular disorders are commonly underlying arthrogryposis. Here, we report five affected individuals from three independent families sharing an overlapping phenotype with congenital contractures affecting shoulder, elbow, hand, hip, knee and foot as well as scoliosis, reduced palmar and plantar skin folds, microcephaly and facial dysmorphism. Using exome sequencing, we identified homozygous truncating variants in FILIP1 in all patients. FILIP1 is a regulator of filamin homeostasis required for the initiation of cortical cell migration in the developing neocortex and essential for the differentiation process of cross-striated muscle cells during myogenesis. In summary, our data indicate that bi-allelic truncating variants in FILIP1 are causative of a novel autosomal recessive disorder and expand the spectrum of genetic factors causative of arthrogryposis multiplex congenita.

Genetic testing in the evaluation of individuals with clinical diagnosis of atypical Sturge-Weber syndrome

Sturge-Weber Syndrome (SWS) is a rare vascular malformation disorder characterized by abnormal blood vessels in the brain, skin, and eye. SWS is most commonly caused by a somatic mosaic GNAQ-p.Arg183Gln variant. In this series, 12 patients presented for clinical evaluation of SWS but were noted to have atypical features, and therefore germline and/or somatic genetic testing was performed. Atypical features included extensive capillary malformation on the body as well as the face, frontal bossing, macrocephaly, telangiectasia, overgrowth of extremities, absence of neurologic signs and symptoms, and family history of vascular malformations. Five patients had a somatic GNAQ or GNA11 pathogenic variant, one patient had a somatic mosaic likely-pathogenic variant in PIK3CA, and another one had a somatic mosaic deletion that disrupted PTPRD. The other five patients had germline variants in RASA1, EPHB4, or KIT. Our findings suggest that patients presenting for SWS evaluation who have atypical clinical characteristics may have pathogenic germline or somatic variants in genes other than GNAQ or GNA11. Broad germline and somatic genetic testing in these patients with atypical findings may have implications for medical care, prognosis, and trial eligibility.

Mechanism of KMT5B haploinsufficiency in neurodevelopment in humans and mice

Pathogenic variants in KMT5B, a lysine methyltransferase, are associated with global developmental delay, macrocephaly, autism, and congenital anomalies (OMIM# 617788). Given the relatively recent discovery of this disorder, it has not been fully characterized. Deep phenotyping of the largest (n = 43) patient cohort to date identified that hypotonia and congenital heart defects are prominent features that were previously not associated with this syndrome. Both missense variants and putative loss-of-function variants resulted in slow growth in patient-derived cell lines. KMT5B homozygous knockout mice were smaller in size than their wild-type littermates but did not have significantly smaller brains, suggesting relative macrocephaly, also noted as a prominent clinical feature. RNA sequencing of patient lymphoblasts and Kmt5b haploinsufficient mouse brains identified differentially expressed pathways associated with nervous system development and function including axon guidance signaling. Overall, we identified additional pathogenic variants and clinical features in KMT5B-related neurodevelopmental disorder and provide insights into the molecular mechanisms of the disorder using multiple model systems.

Venous malformation may be a feature of EXT1-related hereditary multiple exostoses: A report of two unrelated probands

Hereditary multiple exostoses (HME), also known as hereditary multiple osteochondroma (HMO), is an autosomal dominant disorder caused by pathogenic variants in exostosin-1 or -2 (EXT1 or EXT2). It is characterized by the formation of multiple benign growing osteochondromas (exostoses) that most commonly affect the long bones; however, it may also occur throughout the body. Although many of these lesions are clinically asymptomatic, some can lead to chronic pain and skeletal deformities and interfere with adjacent neurovascular structures. Here, we report two unrelated probands that presented with a clinical and molecular diagnosis of HME with venous malformation, a clinical feature not previously reported in individuals with HME.

Systemic artery to pulmonary artery aneurysm malformations associated with variants at MCF2L

Arteriovenous malformations (AVM) are characterized by abnormal vessels connecting arteries and veins resulting in a disruption of normal blood flow. Hereditary hemorrhagic telangiectasia (HHT) is the most common cause of pulmonary AVM characterized by a right to left shunt. Here we describe a distinct malformation where the flow of blood was from a systemic artery to the pulmonary artery (PA) resulting in a left to right shunt instead of the right to left shunt seen in individuals with HHT. This distinct malformation was identified in seven probands, one from a multiplex family containing 10 affected individuals from five generations. To identify the molecular basis of this distinct malformation, we performed exome sequencing (ES) on the seven probands and the affected paternal female cousin from the multiplex family. PhenoDB was used to prioritize candidate causative variants along with burden analysis. We describe the clinical and radiological details of the new systemic artery to PA malformation with or without pulmonary artery aneurysm (SA-PA(A)) and recommend distinct treatment techniques. Moreover, ES analysis revealed possible causative variants identified in three families with variants in a novel candidate disease gene, MCF2L. Further functional studies will be necessary to better understand the molecular mechanisms involved on SA-PA(A) malformation, however our findings suggest that MCF2L is a novel disease gene associated with SA-PA(A).

A Missense Variant in CASKIN1's Proline-Rich Region Segregates with Psychosis in a Three-Generation Family

The polygenic nature of schizophrenia (SCZ) implicates many variants in disease development. Rare variants of high penetrance have been shown to contribute to the disease prevalence. Whole-exome sequencing of a large three-generation family with SCZ and bipolar disorder identified a single segregating novel, rare, non-synonymous variant in the gene CASKIN1. The variant D1204N is absent from all databases, and CASKIN1 has a gnomAD missense score Z = 1.79 and pLI = 1, indicating its strong intolerance to variation. We find that introducing variants in the proline-rich region where the D1204N resides results in significant cellular changes in iPSC-derived neurons, consistent with CASKIN1’s known functions. We observe significant transcriptomic changes in 368 genes (padj < 0.05) involved in neuronal differentiation and nervous system development. We also observed nominally significant changes in the frequency of action potentials during differentiation, where the speed at which the edited and unedited cells reach the same level of activity differs. Our results suggest that CASKIN1 is an excellent gene candidate for psychosis development with high penetrance in this family.

Collaborative Platforms and Matchmaking Algorithms for Research and Education, Establishment, and Optimization of Consortia

Matchmaking has a great position in the rational allocation of resources in several fields, ranging from market operation to people's daily lives. Matchmakers have evolved through artificial intelligence technologies and are being introduced in numerous aspects of industry, research, and academia in solving decision issues, research innovation design, and building robust and efficient networks. The goal of this report is to describe the collaborative platforms and matchmaking algorithms for research and education, as well as the establishment and optimization of consortia.

X-linked variations inSHROOM4are implicated in congenital anomalies of the urinary tract and the anorectal, cardiovascular and central nervous systems

Background

SHROOM4is thought to play an important role in cytoskeletal modification and development of the early nervous system. Previously, single-nucleotide variants (SNVs) or copy number variations (CNVs) inSHROOM4have been associated with the neurodevelopmental disorder Stocco dos Santos syndrome, but not with congenital anomalies of the urinary tract and the visceral or the cardiovascular system.

Methods

Here, exome sequencing and CNV analyses besides expression studies in zebrafish and mouse andknockdown(KD) experiments using a splice blocking morpholino in zebrafish were performed to study the role ofSHROOM4during embryonic development.

Results

In this study, we identified putative disease-causing SNVs and CNVs inSHROOM4in six individuals from four families with congenital anomalies of the urinary tract and the anorectal, cardiovascular and central nervous systems (CNS). Embryonic mouse and zebrafish expression studies showedShroom4expression in the upper and lower urinary tract, the developing cloaca, the heart and the cerebral CNS. KD studies in zebrafish larvae revealed pronephric cysts, anomalies of the cloaca and the heart, decreased eye-to-head ratio and higher mortality compared with controls. These phenotypes could be rescued by co-injection of human wild-typeSHROOM4mRNA and morpholino.

Conclusion

The identified SNVs and CNVs in affected individuals with congenital anomalies of the urinary tract, the anorectal, the cardiovascular and the central nervous systems, and subsequent embryonic mouse and zebrafish studies suggestSHROOM4as a developmental gene for different organ systems.

Heterozygous pathogenic variants involving CBFB cause a new skeletal disorder resembling cleidocranial dysplasia

BACKGROUND

Cleidocranial dysplasia (CCD) is a rare skeletal dysplasia with significant clinical variability. Patients with CCD typically present with delayed closure of fontanels and cranial sutures, dental anomalies, clavicular hypoplasia or aplasia and short stature. Runt-related transcription factor 2 (RUNX2) is currently the only known disease-causing gene for CCD, but several studies have suggested locus heterogeneity.

METHODS

The cohort consists of eight subjects from five unrelated families partially identified through GeneMatcher. Exome or genome sequencing was applied and in two subjects the effect of the variant was investigated at RNA level.

RESULTS

In each subject a heterozygous pathogenic variant in CBFB was detected, whereas no genomic alteration involving RUNX2 was found. Three CBFB variants (one splice site alteration, one nonsense variant, one 2 bp duplication) were shown to result in a premature stop codon. A large intragenic deletion was found to delete exon 4, without affecting CBFB expression. The effect of a second splice site variant could not be determined but most likely results in a shortened or absent protein. Affected individuals showed similarities with RUNX2-related CCD, including dental and clavicular abnormalities. Normal stature and neurocognitive problems were however distinguishing features. CBFB encodes the core-binding factor β subunit, which can interact with all RUNX proteins (RUNX1, RUNX2, RUNX3) to form heterodimeric transcription factors. This may explain the phenotypic differences between CBFB-related and RUNX2-related CCD.

CONCLUSION

We confirm the previously suggested locus heterogeneity for CCD by identifying five pathogenic variants in CBFB in a cohort of eight individuals with clinical and radiographic features reminiscent of CCD.

ModelMatcher: A scientist-centric online platform to facilitate collaborations between stakeholders of rare and undiagnosed disease research

Next-generation sequencing is a prevalent diagnostic tool for undiagnosed diseases and has played a significant role in rare disease gene discovery. Although this technology resolves some cases, others are given a list of possibly damaging genetic variants necessitating functional studies. Productive collaborations between scientists, clinicians, and patients (affected individuals) can help resolve such medical mysteries and provide insights into in vivo function of human genes. Furthermore, facilitating interactions between scientists and research funders, including nonprofit organizations or commercial entities, can dramatically reduce the time to translate discoveries from bench to bedside. Several systems designed to connect clinicians and researchers with a shared gene of interest have been successful. However, these platforms exclude some stakeholders based on their role or geography. Here we describe ModelMatcher, a global online matchmaking tool designed to facilitate cross-disciplinary collaborations, especially between scientists and other stakeholders of rare and undiagnosed disease research. ModelMatcher is integrated into the Rare Diseases Models and Mechanisms Network and Matchmaker Exchange, allowing users to identify potential collaborators in other registries. This living database decreases the time from when a scientist or clinician is making discoveries regarding their genes of interest, to when they identify collaborators and sponsors to facilitate translational and therapeutic research.

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.

Hnrnpul1 controls transcription, splicing, and modulates skeletal and limb development in vivo

Mutations in RNA-binding proteins can lead to pleiotropic phenotypes including craniofacial, skeletal, limb, and neurological symptoms. Heterogeneous nuclear ribonucleoproteins (hnRNPs) are involved in nucleic acid binding, transcription, and splicing through direct binding to DNA and RNA, or through interaction with other proteins in the spliceosome. We show a developmental role for Hnrnpul1 in zebrafish, resulting in reduced body and fin growth and missing bones. Defects in craniofacial tendon growth and adult-onset caudal scoliosis are also seen. We demonstrate a role for Hnrnpul1 in alternative splicing and transcriptional regulation using RNA-sequencing, particularly of genes involved in translation, ubiquitination, and DNA damage. Given its cross-species conservation and role in splicing, it would not be surprising if it had a role in human development. Whole-exome sequencing detected a homozygous frameshift variant in HNRNPUL1 in 2 siblings with congenital limb malformations, which is a candidate gene for their limb malformations. Zebrafish Hnrnpul1 mutants suggest an important developmental role of hnRNPUL1 and provide motivation for exploring the potential conservation of ancient regulatory circuits involving hnRNPUL1 in human development.

Centers for Mendelian Genomics: A decade of facilitating gene discovery

PURPOSE

Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration.

METHODS

Over the past 10 years, the National Institutes of Health-supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution.

RESULTS

We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher.

CONCLUSION

The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

Novel biallelic variants expand the SLC5A6-related phenotypic spectrum

The sodium (Na+):multivitamin transporter (SMVT), encoded by SLC5A6, belongs to the sodium:solute symporter family and is required for the Na+-dependent uptake of biotin (vitamin B7), pantothenic acid (vitamin B5), the vitamin-like substance α-lipoic acid, and iodide. Compound heterozygous SLC5A6 variants have been reported in individuals with variable multisystemic disorder, including failure to thrive, developmental delay, seizures, cerebral palsy, brain atrophy, gastrointestinal problems, immunodeficiency, and/or osteopenia. We expand the phenotypic spectrum associated with biallelic SLC5A6 variants affecting function by reporting five individuals from three families with motor neuropathies. We identified the homozygous variant c.1285 A > G [p.(Ser429Gly)] in three affected siblings and a simplex patient and the maternally inherited c.280 C > T [p.(Arg94*)] variant and the paternally inherited c.485 A > G [p.(Tyr162Cys)] variant in the simplex patient of the third family. Both missense variants were predicted to affect function by in silico tools. 3D homology modeling of the human SMVT revealed 13 transmembrane helices (TMs) and Tyr162 and Ser429 to be located at the cytoplasmic facing region of TM4 and within TM11, respectively. The SLC5A6 missense variants p.(Tyr162Cys) and p.(Ser429Gly) did not affect plasma membrane localization of the ectopically expressed multivitamin transporter suggesting reduced but not abolished function, such as lower catalytic activity. Targeted therapeutic intervention yielded clinical improvement in four of the five patients. Early molecular diagnosis by exome sequencing is essential for timely replacement therapy in affected individuals.

The impact of GeneMatcher on international data sharing and collaboration

GeneMatcher (genematcher.org) is a tool designed to connect individuals with an interest in the same gene. Now used around the world to create collaborations and generate the evidence needed to support novel disease gene identification, GeneMatcher is a founding member of the Matchmaker Exchange (MME; matchmakerexchange.org) and strongest possible advocate for global data sharing including those in resource‐limited environments. As of October 1, 2021, there are 12,531 submitters from 94 countries who have submitted 58,134 submissions with 13,498 unique genes in the database. Among these genes, 8970 (64%) have matched at least once and the total number of matches is 378,806, growing by about 10,000 per month. GeneMatcher submitters increase by 80–120 each month and submissions grow by >800 per month, while unique genes and gene matches continue to grow steadily at rate of about 80 per month. The number of genes without a match peaked at 4371 in February of 2019 and despite the increase in the number of new submissions, the number of unique genes without a match continues to slowly decline, currently standing at 4,016. All submissions in GeneMatcher are available for matching across the MME.

Diagnostic testing laboratories are valuable partners for disease gene discovery: 5-year experience with GeneMatcher

Although the rates of disease gene discovery have steadily increased with the expanding use of genome and exome sequencing by clinical and research laboratories, only ~16% of genes in the genome have confirmed disease associations. Here we describe our clinical laboratory's experience utilizing GeneMatcher, an online portal designed to promote disease gene discovery and data sharing. Since 2016, we submitted 246 candidates from 243 unique genes to GeneMatcher, of which 111 (45%) are now clinically characterized. Submissions meeting our candidate gene‐reporting criteria based on a scoring system using patient and molecular‐weighted evidence were significantly more likely to be characterized as of October 2021 versus genes that did not meet our clinical‐reporting criteria (p = 0.025). We reported relevant findings related to these newly characterized gene–disease associations in 477 probands. In 218 (46%) instances, we issued reclassifications after an initial negative or candidate gene (uncertain) report. We coauthored 104 publications delineating gene–disease relationships, including descriptions of new associations (60%), additional supportive evidence (13%), subsequent descriptive cohorts (23%), and phenotypic expansions (4%). Clinical laboratories are pivotal for disease gene discovery efforts and can screen phenotypes based on genotype matches, contact clinicians of relevant cases, and issue proactive reclassification reports.

A clinician's guide to omics resources in dermatology

With recent advances in high‐throughput technologies, we are now in an era where the use of large‐scale datasets of biological samples and individual diseases can be analysed using omics methodologies. These include genomics, transcriptomics, proteomics, metabolomics, lipidomics and epigenomics. Omics approaches have been developed to deliver a holistic understanding of systems biology, to identify key biomarkers, and to aid in the interpretation of molecular, biochemical and environmental interactions. Navigating through the plethora of online datasets to find useful and concise information for comparison of data can be complex and overwhelming. The purpose of this article is to review the current repositories and databases, and to evaluate their application in dermatological research and their relevance to clinical practice. For this study, an extensive review of online platforms used in dermatology research was undertaken. Online resources for genetic disease information, genetic disease connection platforms for patients and researchers, clinical interpretation of variants, genome and DNA databases, and omics data repositories and resources were collected. This study provides a comprehensive overview of relevant databases that will aid clinicians and scientists using omics data in dermatology.

Bi-allelic variants in OGDHL cause a neurodevelopmental spectrum disease featuring epilepsy, hearing loss, visual impairment, and ataxia

The 2-oxoglutarate dehydrogenase-like (OGDHL) protein is a rate-limiting enzyme in the Krebs cycle that plays a pivotal role in mitochondrial metabolism. OGDHL expression is restricted mainly to the brain in humans. Here, we report nine individuals from eight unrelated families carrying bi-allelic variants in OGDHL with a range of neurological and neurodevelopmental phenotypes including epilepsy, hearing loss, visual impairment, gait ataxia, microcephaly, and hypoplastic corpus callosum. The variants include three homozygous missense variants (p.Pro852Ala, p.Arg244Trp, and p.Arg299Gly), three compound heterozygous single-nucleotide variants (p.Arg673Gln/p.Val488Val, p.Phe734Ser/p.Ala327Val, and p.Trp220Cys/p.Asp491Val), one homozygous frameshift variant (p.Cys553Leufs^∗16), and one homozygous stop-gain variant (p.Arg440Ter). To support the pathogenicity of the variants, we developed a novel CRISPR-Cas9-mediated tissue-specific knockout with cDNA rescue system for dOgdh, the Drosophila ortholog of human OGDHL. Pan-neuronal knockout of dOgdh led to developmental lethality as well as defects in Krebs cycle metabolism, which was fully rescued by expression of wild-type dOgdh. Studies using the Drosophila system indicate that p.Arg673Gln, p.Phe734Ser, and p.Arg299Gly are severe loss-of-function alleles, leading to developmental lethality, whereas p.Pro852Ala, p.Ala327Val, p.Trp220Cys, p.Asp491Val, and p.Arg244Trp are hypomorphic alleles, causing behavioral defects. Transcript analysis from fibroblasts obtained from the individual carrying the synonymous variant (c.1464T>C [p.Val488Val]) in family 2 showed that the synonymous variant affects splicing of exon 11 in OGDHL. Human neuronal cells with OGDHL knockout exhibited defects in mitochondrial respiration, indicating the essential role of OGDHL in mitochondrial metabolism in humans. Together, our data establish that the bi-allelic variants in OGDHL are pathogenic, leading to a Mendelian neurodevelopmental disease in humans.

Inhibition of G-protein signalling in cardiac dysfunction of intellectual developmental disorder with cardiac arrhythmia (IDDCA) syndrome

BACKGROUND

Pathogenic variants of GNB5 encoding the β5 subunit of the guanine nucleotide-binding protein cause IDDCA syndrome, an autosomal recessive neurodevelopmental disorder associated with cognitive disability and cardiac arrhythmia, particularly severe bradycardia.

METHODS

We used echocardiography and telemetric ECG recordings to investigate consequences of Gnb5 loss in mouse.

RESULTS

We delineated a key role of Gnb5 in heart sinus conduction and showed that Gnb5-inhibitory signalling is essential for parasympathetic control of heart rate (HR) and maintenance of the sympathovagal balance. Gnb5-/- mice were smaller and had a smaller heart than Gnb5+/+ and Gnb5+/- , but exhibited better cardiac function. Lower autonomic nervous system modulation through diminished parasympathetic control and greater sympathetic regulation resulted in a higher baseline HR in Gnb5-/- mice. In contrast, Gnb5-/- mice exhibited profound bradycardia on treatment with carbachol, while sympathetic modulation of the cardiac stimulation was not altered. Concordantly, transcriptome study pinpointed altered expression of genes involved in cardiac muscle contractility in atria and ventricles of knocked-out mice. Homozygous Gnb5 loss resulted in significantly higher frequencies of sinus arrhythmias. Moreover, we described 13 affected individuals, increasing the IDDCA cohort to 44 patients.

CONCLUSIONS

Our data demonstrate that loss of negative regulation of the inhibitory G-protein signalling causes HR perturbations in Gnb5-/- mice, an effect mainly driven by impaired parasympathetic activity. We anticipate that unravelling the mechanism of Gnb5 signalling in the autonomic control of the heart will pave the way for future drug screening.

Biallelic FRA10AC1 variants cause a neurodevelopmental disorder with growth retardation

The major spliceosome mediates pre-mRNA splicing by recognizing the highly conserved sequences at the 5' and 3' splice sites and the branch point. More than 150 proteins participate in the splicing process and are organized in the spliceosomal A, B, and C complexes. FRA10AC1 is a peripheral protein of the spliceosomal C complex and its ortholog in the green alga facilitates recognition or interaction with splice sites. We identified biallelic pathogenic variants in FRA10AC1 in five individuals from three consanguineous families. The two unrelated patients 1 and 2 with loss-of-function variants showed developmental delay, intellectual disability, and no speech, while three siblings with the c.494_496delAAG (p.Glu165del) variant had borderline to mild intellectual disability. All patients had microcephaly, hypoplasia or agenesis of the corpus callosum, growth retardation, and craniofacial dysmorphism. FRA10AC1 transcripts and proteins were drastically reduced or absent in fibroblasts of patients 1 and 2. In a heterologous expression system, the p. Glu165del variant impacts intrinsic stability of FRA10AC1 but does not affect its nuclear localization. By co-immunoprecipitation, we found ectopically expressed HA-FRA10AC1 in complex with endogenous DGCR14, another component of the spliceosomal C complex, while the splice factors CHERP, NKAP, RED, and SF3B2 could not be co-immunoprecipitated. Using an in vitro splicing reporter assay, we did not obtain evidence for FRA10AC1 deficiency to suppress missplicing events caused by mutations in the highly conserved dinucleotides of 5' and 3' splice sites in an in vitro splicing assay in patient-derived fibroblasts. Our data highlight the importance of specific peripheral spliceosomal C complex proteins for neurodevelopment. It remains possible that FRA10AC1 may have other and/or additional cellular functions, such as coupling of transcription and splicing reactions.

Phenotypic Trade-Offs: Deciphering the Impact of Neurodiversity on Drug Development in Fragile X Syndrome

Fragile X syndrome (FXS) is the most common single-gene cause of intellectual disability and autism spectrum disorder. Individuals with FXS present with a wide range of severity in multiple phenotypes including cognitive delay, behavioral challenges, sleep issues, epilepsy, and anxiety. These symptoms are also shared by many individuals with other neurodevelopmental disorders (NDDs). Since the discovery of the FXS gene, FMR1, FXS has been the focus of intense preclinical investigation and is placed at the forefront of clinical trials in the field of NDDs. So far, most studies have aimed to translate the rescue of specific phenotypes in animal models, for example, learning, or improving general cognitive or behavioral functioning in individuals with FXS. Trial design, selection of outcome measures, and interpretation of results of recent trials have shown limitations in this type of approach. We propose a new paradigm in which all phenotypes involved in individuals with FXS would be considered and, more importantly, the possible interactions between these phenotypes. This approach would be implemented both at the baseline, meaning when entering a trial or when studying a patient population, and also after the intervention when the study subjects have been exposed to the investigational product. This approach would allow us to further understand potential trade-offs underlying the varying effects of the treatment on different individuals in clinical trials, and to connect the results to individual genetic differences. To better understand the interplay between different phenotypes, we emphasize the need for preclinical studies to investigate various interrelated biological and behavioral outcomes when assessing a specific treatment. In this paper, we present how such a conceptual shift in preclinical design could shed new light on clinical trial results. Future clinical studies should take into account the rich neurodiversity of individuals with FXS specifically and NDDs in general, and incorporate the idea of trade-offs in their designs.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population

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.

Genome sequencing data analysis for rare disease gene discovery

Rare diseases occur in a smaller proportion of the general population, which is variedly defined as less than 200 000 individuals (US) or in less than 1 in 2000 individuals (Europe). Although rare, they collectively make up to approximately 7000 different disorders, with majority having a genetic origin, and affect roughly 300 million people globally. Most of the patients and their families undergo a long and frustrating diagnostic odyssey. However, advances in the field of genomics have started to facilitate the process of diagnosis, though it is hindered by the difficulty in genome data analysis and interpretation. A major impediment in diagnosis is in the understanding of the diverse approaches, tools and datasets available for variant prioritization, the most important step in the analysis of millions of variants to select a few potential variants. Here we present a review of the latest methodological developments and spectrum of tools available for rare disease genetic variant discovery and recommend appropriate data interpretation methods for variant prioritization. We have categorized the resources based on various steps of the variant interpretation workflow, starting from data processing, variant calling, annotation, filtration and finally prioritization, with a special emphasis on the last two steps. The methods discussed here pertain to elucidating the genetic basis of disease in individual patient cases via trio- or family-based analysis of the genome data. We advocate the use of a combination of tools and datasets and to follow multiple iterative approaches to elucidate the potential causative variant.

IFIH1 loss-of-function variants contribute to very early-onset inflammatory bowel disease

Genetic defects of innate immunity impairing intestinal bacterial sensing are linked to the development of Inflammatory Bowel Disease (IBD). Although much evidence supports a role of the intestinal virome in gut homeostasis, most studies focus on intestinal viral composition rather than on host intestinal viral sensitivity. To demonstrate the association between the development of Very Early Onset IBD (VEOIBD) and variants in the IFIH1 gene which encodes MDA5, a key cytosolic sensor for viral nucleic acids. Whole exome sequencing (WES) was performed in two independent cohorts of children with VEOIBD enrolled in Italy (n = 18) and USA (n = 24). Luciferase reporter assays were employed to assess MDA5 activity. An enrichment analysis was performed on IFIH1 comparing 42 VEOIBD probands with 1527 unrelated individuals without gastrointestinal or immunological issues. We identified rare, likely loss-of-function (LoF), IFIH1 variants in eight patients with VEOIBD from a combined cohort of 42 children. One subject, carrying a homozygous truncating variant resulting in complete LoF, experienced neonatal-onset, pan-gastrointestinal, IBD-like enteropathy plus multiple infectious episodes. The remaining seven subjects, affected by VEOIBD without immunodeficiency, were carriers of one LoF variant in IFIH1. Among these, two patients also carried a second hypomorphic variant, with partial function apparent when MDA5 was weakly stimulated. Furthermore, IFIH1 variants were significantly enriched in children with VEOIBD as compared to controls (p = 0.007). Complete and partial MDA5 deficiency is associated with VEOIBD with variable penetrance and expressivity, suggesting a role for impaired intestinal viral sensing in IBD pathogenesis.

PhenoDB, GeneMatcher and VariantMatcher, tools for analysis and sharing of sequence data

Background

With the advent of whole exome (ES) and genome sequencing (GS) as tools for disease gene discovery, rare variant filtering, prioritization and data sharing have become essential components of the search for disease genes and variants potentially contributing to disease phenotypes. The computational storage, data manipulation, and bioinformatic interpretation of thousands to millions of variants identified in ES and GS, respectively, is a challenging task. To aid in that endeavor, we constructed PhenoDB, GeneMatcher and VariantMatcher.

Results

PhenoDB is an accessible, freely available, web-based platform that allows users to store, share, analyze and interpret their patients’ phenotypes and variants from ES/GS data. GeneMatcher is accessible to all stakeholders as a web-based tool developed to connect individuals (researchers, clinicians, health care providers and patients) around the globe with interest in the same gene(s), variant(s) or phenotype(s). Finally, VariantMatcher was developed to enable public sharing of variant-level data and phenotypic information from individuals sequenced as part of multiple disease gene discovery projects. Here we provide updates on PhenoDB and GeneMatcher applications and implementation and introduce VariantMatcher.

Conclusion

Each of these tools has facilitated worldwide data sharing and data analysis and improved our ability to connect genes to phenotypic traits. Further development of these platforms will expand variant analysis, interpretation, novel disease-gene discovery and facilitate functional annotation of the human genome for clinical genomics implementation and the precision medicine initiative.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-021-01916-z.

Loss of function mutations in CCDC32 cause a congenital syndrome characterized by craniofacial, cardiac and neurodevelopmental anomalies

Despite the wide use of genomics to investigate the molecular basis of rare congenital malformations, a significant fraction of patients remains bereft of diagnosis. As part of our continuous effort to recruit and perform genomic and functional studies on such cohorts, we investigated the genetic and mechanistic cause of disease in two independent consanguineous families affected by overlapping craniofacial, cardiac, laterality and neurodevelopmental anomalies. Using whole exome sequencing, we identified homozygous frameshift CCDC32 variants in three affected individuals. Functional analysis in a zebrafish model revealed that ccdc32 depletion recapitulates the human phenotypes. Because some of the patient phenotypes overlap defects common to ciliopathies, we asked if loss of CCDC32 might contribute to the dysfunction of this organelle. Consistent with this hypothesis, we show that ccdc32 is required for normal cilia formation in zebrafish embryos and mammalian cell culture, arguing that ciliary defects are at least partially involved in the pathomechanism of this disorder.

Exome survey of individuals affected by VATER/VACTERL with renal phenotypes identifies phenocopies and novel candidate genes

The acronym VATER/VACTERL refers to the rare nonrandom association of the following component features (CFs): vertebral defects (V), anorectal malformations (ARM) (A), cardiac anomalies (C), tracheoesophageal fistula with or without esophageal atresia (TE), renal malformations (R), and limb anomalies (L). For the clinical diagnosis, the presence of at least three CFs is required, individuals presenting with only two CFs have been categorized as VATER/VACTERL-like. The majority of VATER/VACTERL individuals displays a renal phenotype. Hitherto, variants in FGF8, FOXF1, HOXD13, LPP, TRAP1, PTEN, and ZIC3 have been associated with the VATER/VACTERL association; however, large-scale re-sequencing could only confirm TRAP1 and ZIC3 as VATER/VACTERL disease genes, both associated with a renal phenotype. In this study, we performed exome sequencing in 21 individuals and their families with a renal VATER/VACTERL or VATER/VACTERL-like phenotype to identify potentially novel genetic causes. Exome analysis identified biallelic and X-chromosomal hemizygous potentially pathogenic variants in six individuals (29%) in B9D1, FREM1, ZNF157, SP8, ACOT9, and TTLL11, respectively. The online tool GeneMatcher revealed another individual with a variant in ZNF157. Our study suggests six biallelic and X-chromosomal hemizygous VATER/VACTERL disease genes implicating all six genes in the expression of human renal malformations.

Variants in LSM7 impair LSM complexes assembly, neurodevelopment in zebrafish and may be associated with an ultra-rare neurological disease

Leukodystrophies, genetic neurodevelopmental and/or neurodegenerative disorders of cerebral white matter, result from impaired myelin homeostasis and metabolism. Numerous genes have been implicated in these heterogeneous disorders; however, many individuals remain without a molecular diagnosis. Using whole-exome sequencing, biallelic variants in LSM7 were uncovered in two unrelated individuals, one with a leukodystrophy and the other who died in utero. LSM7 is part of the two principle LSM protein complexes in eukaryotes, namely LSM1-7 and LSM2-8. Here, we investigate the molecular and functional outcomes of these LSM7 biallelic variants in vitro and in vivo. Affinity purification-mass spectrometry of the LSM7 variants showed defects in the assembly of both LSM complexes. Lsm7 knockdown in zebrafish led to central nervous system defects, including impaired oligodendrocyte development and motor behavior. Our findings demonstrate that variants in LSM7 cause misassembly of the LSM complexes, impair neurodevelopment of the zebrafish, and may be implicated in human disease. The identification of more affected individuals is needed before the molecular mechanisms of mRNA decay and splicing regulation are added to the categories of biological dysfunctions implicated in leukodystrophies, neurodevelopmental and/or neurodegenerative diseases.

Germline ERBB2/HER2 Coding Variants Are Associated with Increased Risk of Myeloproliferative Neoplasms

Familial cases of myeloproliferative neoplasms (MPN) are relatively common, yet few inherited risk factors have been identified. Exome sequencing of a kindred with a familial cancer syndrome characterized by both MPN and melanoma produced a germline variant in the ERBB2/HER2 gene that co-segregates with disease. To further investigate whether germline ERBB2 variants contribute to MPN predisposition, the frequency of ERBB2 variants was analyzed in 1604 cases that underwent evaluation for hematologic malignancy, including 236 cases of MPN. MPN cases had a higher frequency of rare germline ERBB2 coding variants compared to non-MPN hematologic malignancies (8.9% vs. 4.1%, OR 2.4, 95% CI: 1.4 to 4.0, p = 0.0028) as well as cases without a blood cancer diagnosis that served as an internal control (8.9% vs. 2.7%, OR 3.5, 95% CI: 1.4 to 8.3, p = 0.0053). This finding was validated via comparison to an independent control cohort of 1587 cases without selection for hematologic malignancy (8.9% in MPN cases vs. 5.2% in controls, p = 0.040). The most frequent variant identified, ERBB2 c.1960A > G; p.I654V, was present in MPN cases at more than twice its expected frequency. These data indicate that rare germline coding variants in ERBB2 are associated with an increased risk for development of MPN. The ERBB2 gene is a novel susceptibility locus which likely contributes to cancer risk in combination with additional risk alleles.

Biallelic variants in VPS50 cause a neurodevelopmental disorder with neonatal cholestasis

Golgi-associated retrograde protein (GARP) and endosome-associated recycling protein (EARP) complexes are membrane-tethering heterotetramers located at the trans-Golgi network and recycling endosomes, respectively. GARP and EARP share the three subunits VPS51, VPS52, and VPS53, while VPS50 is unique to EARP and VPS54 to GARP. Retrograde transport of endosomal cargos to the TGN is mediated by GARP and endocytic recycling by EARP. Here we report two unrelated individuals with homozygous variants in VPS50, a splice variant (c.1978-1G>T) and an in-frame deletion (p.Thr608del). Both patients had severe developmental delay, postnatal microcephaly, corpus callosum hypoplasia, seizures and irritability, transient neonatal cholestasis, and failure to thrive. Light and transmission electron microscopy of liver from one revealed absence of gamma-glutamyltransferase at bile canaliculi, with mislocalization to basolateral membranes, and abnormal tight junctions. Using patient-derived fibroblasts, we identified reduced VPS50 protein accompanied by reduced levels of VPS52 and VPS53. While transferrin-receptor internalization rate was normal in cells of both patients, recycling of the receptor to the plasma membrane was significantly delayed. These data underscore the importance of VPS50 and/or the EARP complex in endocytic recycling and suggest an additional function in establishing cell polarity and trafficking between basolateral and apical membranes in hepatocytes. Individuals with biallelic hypomorphic variants in VPS50, VPS51 or VPS53 show an overarching neurodegenerative disorder with severe developmental delay, intellectual disability, microcephaly, early-onset epilepsy, and variable atrophy of the cerebellum, cerebrum, and/or brainstem. The term "GARP/EARP deficiency" designates disorders in such individuals.

Genotype-phenotype correlations and novel molecular insights into the DHX30-associated neurodevelopmental disorders

Background

We aimed to define the clinical and variant spectrum and to provide novel molecular insights into the DHX30-associated neurodevelopmental disorder.

Methods

Clinical and genetic data from affected individuals were collected through Facebook-based family support group, GeneMatcher, and our network of collaborators. We investigated the impact of novel missense variants with respect to ATPase and helicase activity, stress granule (SG) formation, global translation, and their effect on embryonic development in zebrafish. SG formation was additionally analyzed in CRISPR/Cas9-mediated DHX30-deficient HEK293T and zebrafish models, along with in vivo behavioral assays.

Results

We identified 25 previously unreported individuals, ten of whom carry novel variants, two of which are recurrent, and provide evidence of gonadal mosaicism in one family. All 19 individuals harboring heterozygous missense variants within helicase core motifs (HCMs) have global developmental delay, intellectual disability, severe speech impairment, and gait abnormalities. These variants impair the ATPase and helicase activity of DHX30, trigger SG formation, interfere with global translation, and cause developmental defects in a zebrafish model. Notably, 4 individuals harboring heterozygous variants resulting either in haploinsufficiency or truncated proteins presented with a milder clinical course, similar to an individual harboring a de novo mosaic HCM missense variant. Functionally, we established DHX30 as an ATP-dependent RNA helicase and as an evolutionary conserved factor in SG assembly. Based on the clinical course, the variant location, and type we establish two distinct clinical subtypes. DHX30 loss-of-function variants cause a milder phenotype whereas a severe phenotype is caused by HCM missense variants that, in addition to the loss of ATPase and helicase activity, lead to a detrimental gain-of-function with respect to SG formation. Behavioral characterization of dhx30-deficient zebrafish revealed altered sleep-wake activity and social interaction, partially resembling the human phenotype.

Conclusions

Our study highlights the usefulness of social media to define novel Mendelian disorders and exemplifies how functional analyses accompanied by clinical and genetic findings can define clinically distinct subtypes for ultra-rare disorders. Such approaches require close interdisciplinary collaboration between families/legal representatives of the affected individuals, clinicians, molecular genetics diagnostic laboratories, and research laboratories.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00900-3.

Loss of function mutations in GEMIN5 cause a neurodevelopmental disorder

GEMIN5, an RNA-binding protein is essential for assembly of the survival motor neuron (SMN) protein complex and facilitates the formation of small nuclear ribonucleoproteins (snRNPs), the building blocks of spliceosomes. Here, we have identified 30 affected individuals from 22 unrelated families presenting with developmental delay, hypotonia, and cerebellar ataxia harboring biallelic variants in the GEMIN5 gene. Mutations in GEMIN5 perturb the subcellular distribution, stability, and expression of GEMIN5 protein and its interacting partners in patient iPSC-derived neurons, suggesting a potential loss-of-function mechanism. GEMIN5 mutations result in disruption of snRNP complex assembly formation in patient iPSC neurons. Furthermore, knock down of rigor mortis, the fly homolog of human GEMIN5, leads to developmental defects, motor dysfunction, and a reduced lifespan. Interestingly, we observed that GEMIN5 variants disrupt a distinct set of transcripts and pathways as compared to SMA patient neurons, suggesting different molecular pathomechanisms. These findings collectively provide evidence that pathogenic variants in GEMIN5 perturb physiological functions and result in a neurodevelopmental delay and ataxia syndrome.

Next-generation sequencing and the evolution of data sharing

Disease gene identification often relies on identifying multiple affected individuals with similar phenotypes and candidate variants in the same gene. Phenotypic and genomic data sharing tools have facilitated connections that led to novel disease gene discoveries and better characterization and recognition of rare diseases. Additionally, data sharing has evolved. From gene-based matches to variant-level information with increasing use of phenotypic information. We expect that these initiatives will continue to expand in the future affording clinicians, researchers, and most importantly, patients and their families faster and more comprehensive answers.

Pathogenic alleles in microtubule, secretory granule and extracellular matrix-related genes in familial keratoconus

Keratoconus is a common corneal defect with a complex genetic basis. By whole exome sequencing of affected members from 11 multiplex families of European ancestry, we identified 23 rare, heterozygous, potentially pathogenic variants in 8 genes. These include nonsynonymous single amino acid substitutions in HSPG2, EML6 and CENPF in two families each, and in NBEAL2, LRP1B, PIK3CG and MRGPRD in three families each; ITGAX had nonsynonymous single amino acid substitutions in two families and an indel with a base substitution producing a nonsense allele in the third family. Only HSPG2, EML6 and CENPF have been associated with ocular phenotypes previously. With the exception of MRGPRD and ITGAX, we detected the transcript and encoded protein of the remaining genes in the cornea and corneal cell cultures. Cultured stromal cells showed cytoplasmic punctate staining of NBEAL2, staining of the fibrillar cytoskeletal network by EML6, while CENPF localized to the basal body of primary cilia. We inhibited the expression of HSPG2, EML6, NBEAL2 and CENPF in stromal cell cultures and assayed for the expression of COL1A1 as a readout of corneal matrix production. An upregulation in COL1A1 after siRNA inhibition indicated their functional link to stromal cell biology. For ITGAX, encoding a leukocyte integrin, we assayed its level in the sera of 3 affected families compared with 10 unrelated controls to detect an increase in all affecteds. Our study identified genes that regulate the cytoskeleton, protein trafficking and secretion, barrier tissue function and response to injury and inflammation, as being relevant to keratoconus.

De novo variants in neurodevelopmental disorders-experiences from a tertiary care center

Up to 40% of neurodevelopmental disorders (NDDs) such as intellectual disability, developmental delay, autism spectrum disorder, and developmental motor abnormalities have a documented underlying monogenic defect, primarily due to de novo variants. Still, the overall burden of de novo variants as well as novel disease genes in NDDs await discovery. We performed parent-offspring trio exome sequencing in 231 individuals with NDDs. Phenotypes were compiled using human phenotype ontology terms. The overall diagnostic yield was 49.8% (n = 115/231) with de novo variants contributing to more than 80% (n = 93/115) of all solved cases. De novo variants affected 72 different-mostly constrained-genes. In addition, we identified putative pathogenic variants in 16 genes not linked to NDDs to date. Reanalysis performed in 80 initially unsolved cases revealed a definitive diagnosis in two additional cases. Our study consolidates the contribution and genetic heterogeneity of de novo variants in NDDs highlighting trio exome sequencing as effective diagnostic tool for NDDs. Besides, we illustrate the potential of a trio-approach for candidate gene discovery and the power of systematic reanalysis of unsolved cases.

Early-infantile onset epilepsy and developmental delay caused by bi-allelic GAD1 variants

Gamma-aminobutyric acid (GABA) and glutamate are the most abundant amino acid neurotransmitters in the brain. GABA, an inhibitory neurotransmitter, is synthesized by glutamic acid decarboxylase (GAD). Its predominant isoform GAD67, contributes up to ∼90% of base-level GABA in the CNS, and is encoded by the GAD1 gene. Disruption of GAD1 results in an imbalance of inhibitory and excitatory neurotransmitters, and as Gad1-/- mice die neonatally of severe cleft palate, it has not been possible to determine any potential neurological dysfunction. Furthermore, little is known about the consequence of GAD1 disruption in humans. Here we present six affected individuals from six unrelated families, carrying bi-allelic GAD1 variants, presenting with developmental and epileptic encephalopathy, characterized by early-infantile onset epilepsy and hypotonia with additional variable non-CNS manifestations such as skeletal abnormalities, dysmorphic features and cleft palate. Our findings highlight an important role for GAD1 in seizure induction, neuronal and extraneuronal development, and introduce GAD1 as a new gene associated with developmental and epileptic encephalopathy.

Bi-allelic GAD1 variants cause a neonatal onset syndromic developmental and epileptic encephalopathy

Developmental and epileptic encephalopathies are a heterogeneous group of early-onset epilepsy syndromes dramatically impairing neurodevelopment. Modern genomic technologies have revealed a number of monogenic origins and opened the door to therapeutic hopes. Here we describe a new syndromic developmental and epileptic encephalopathy caused by bi-allelic loss-of-function variants in GAD1, as presented by 11 patients from six independent consanguineous families. Seizure onset occurred in the first 2 months of life in all patients. All 10 patients, from whom early disease history was available, presented with seizure onset in the first month of life, mainly consisting of epileptic spasms or myoclonic seizures. Early EEG showed suppression-burst or pattern of burst attenuation or hypsarrhythmia if only recorded in the post-neonatal period. Eight patients had joint contractures and/or pes equinovarus. Seven patients presented a cleft palate and two also had an omphalocele, reproducing the phenotype of the knockout Gad1-/- mouse model. Four patients died before 4 years of age. GAD1 encodes the glutamate decarboxylase enzyme GAD67, a critical actor of the γ-aminobutyric acid (GABA) metabolism as it catalyses the decarboxylation of glutamic acid to form GABA. Our findings evoke a novel syndrome related to GAD67 deficiency, characterized by the unique association of developmental and epileptic encephalopathies, cleft palate, joint contractures and/or omphalocele.