Bi-allelic loss-of-function variants in WBP4, encoding a spliceosome protein, result in a variable neurodevelopmental syndrome.

Over two dozen spliceosome proteins are involved in human diseases, also referred to as spliceosomopathies. WW domain-binding protein 4 (WBP4) is part of the early spliceosomal complex and has not been previously associated with human pathologies in the Online Mendelian Inheritance in Man (OMIM) database. Through GeneMatcher, we identified ten individuals from eight families with a severe neurodevelopmental syndrome featuring variable manifestations. Clinical manifestations included hypotonia, global developmental delay, severe intellectual disability, brain abnormalities, musculoskeletal, and gastrointestinal abnormalities. Genetic analysis revealed five different homozygous loss-of-function variants in WBP4. Immunoblotting on fibroblasts from two affected individuals with different genetic variants demonstrated a complete loss of protein, and RNA sequencing analysis uncovered shared abnormal splicing patterns, including in genes associated with abnormalities of the nervous system, potentially underlying the phenotypes of the probands. We conclude that bi-allelic variants in WBP4 cause a developmental disorder with variable presentations, adding to the growing list of human spliceosomopathies.

A neurodevelopmental disorder associated with a loss-of-function missense mutation in RAB35.

Rab35 (Ras-associated binding protein) is a small GTPase that regulates endosomal membrane trafficking and functions in cell polarity, cytokinesis, and growth factor signaling. Altered Rab35 function contributes to progression of glioblastoma, defects in primary cilia formation, and altered cytokinesis. Here, we report a pediatric patient with global developmental delay, hydrocephalus, a Dandy-Walker malformation, axial hypotonia with peripheral hypertonia, visual problems, and conductive hearing impairment. Exome sequencing identified a homozygous missense variant in the GTPase fold of RAB35 (c.80G>A; p.R27H) as the most likely candidate. Functional analysis of the R27H-Rab35 variant protein revealed enhanced interaction with its guanine-nucleotide exchange factor, DENND1A and decreased interaction with a known effector, MICAL1, indicating that the protein is in an inactive conformation. Cellular expression of the variant drives the activation of Arf6, a small GTPase under negative regulatory control of Rab35. Importantly, variant expression leads to delayed cytokinesis and altered length, number, and Arl13b composition of primary cilia, known factors in neurodevelopmental disease. Our findings provide evidence of altered Rab35 function as a causative factor of a neurodevelopmental disorder.

Rare germline heterozygous missense variants in BRCA1-associated protein 1, BAP1, cause a syndromic neurodevelopmental disorder.

Nuclear deubiquitinase BAP1 (BRCA1-associated protein 1) is a core component of multiprotein complexes that promote transcription by reversing the ubiquitination of histone 2A (H2A). BAP1 is a tumor suppressor whose germline loss-of-function variants predispose to cancer. To our knowledge, there are very rare examples of different germline variants in the same gene causing either a neurodevelopmental disorder (NDD) or a tumor predisposition syndrome. Here, we report a series of 11 de novo germline heterozygous missense BAP1 variants associated with a rare syndromic NDD. Functional analysis showed that most of the variants cannot rescue the consequences of BAP1 inactivation, suggesting a loss-of-function mechanism. In T cells isolated from two affected children, H2A deubiquitination was impaired. In matching peripheral blood mononuclear cells, histone H3 K27 acetylation ChIP-seq indicated that these BAP1 variants induced genome-wide chromatin state alterations, with enrichment for regulatory regions surrounding genes of the ubiquitin-proteasome system (UPS). Altogether, these results define a clinical syndrome caused by rare germline missense BAP1 variants that alter chromatin remodeling through abnormal histone ubiquitination and lead to transcriptional dysregulation of developmental genes.

Neurodevelopmental disorder mutations in the purine biosynthetic enzyme IMPDH2 disrupt its allosteric regulation.

Inosine 5' monophosphate dehydrogenase (IMPDH) is a critical regulatory enzyme in purine nucleotide biosynthesis that is inhibited by the downstream product GTP. Multiple point mutations in the human isoform IMPDH2 have recently been associated with dystonia and other neurodevelopmental disorders, but the effect of the mutations on enzyme function has not been described. Here, we report the identification of two additional missense variants in IMPDH2 from affected individuals and show that all of the disease-associated mutations disrupt GTP regulation. Cryo-EM structures of one IMPDH2 mutant suggest this regulatory defect arises from a shift in the conformational equilibrium toward a more active state. This structural and functional analysis provides insight into IMPDH2-associated disease mechanisms that point to potential therapeutic approaches and raises new questions about fundamental aspects of IMPDH regulation.

Unbiased phenotype and genotype matching maximizes gene discovery and diagnostic yield.

PURPOSE: Widespread application of next-generation sequencing, combined with data exchange platforms, has provided molecular diagnoses for countless families. To maximize diagnostic yield, we implemented an unbiased semi-automated genematching algorithm based on genotype and phenotype matching. METHODS: Rare homozygous variants identified in 2 or more affected individuals, but not in healthy individuals, were extracted from our local database of ∼12,000 exomes. Phenotype similarity scores (PSS), based on human phenotype ontology terms, were assigned to each pair of individuals matched at the genotype level using HPOsim. RESULTS: 33,792 genotype-matched pairs were discovered, representing variants in 7567 unique genes. There was an enrichment of PSS ≥0.1 among pathogenic/likely pathogenic variant-level pairs (94.3% in pathogenic/likely pathogenic variant-level matches vs 34.75% in all matches). We highlighted founder or region-specific variants as an internal positive control and proceeded to identify candidate disease genes. Variant-level matches were particularly helpful in cases involving inframe indels and splice region variants beyond the canonical splice sites, which may otherwise have been disregarded, allowing for detection of candidate disease genes, such as KAT2A, RPAIN, and LAMP3. CONCLUSION: Semi-automated genotype matching combined with PSS is a powerful tool to resolve variants of uncertain significance and to identify candidate disease genes.

Mutations in MYLPF Cause a Novel Segmental Amyoplasia that Manifests as Distal Arthrogryposis.

We identified ten persons in six consanguineous families with distal arthrogryposis (DA) who had congenital contractures, scoliosis, and short stature. Exome sequencing revealed that each affected person was homozygous for one of two different rare variants (c.470G>T [p.Cys157Phe] or c.469T>C [p.Cys157Arg]) affecting the same residue of myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF). In a seventh family, a c.487G>A (p.Gly163Ser) variant in MYLPF arose de novo in a father, who transmitted it to his son. In an eighth family comprised of seven individuals with dominantly inherited DA, a c.98C>T (p.Ala33Val) variant segregated in all four persons tested. Variants in MYLPF underlie both dominant and recessively inherited DA. Mylpf protein models suggest that the residues associated with dominant DA interact with myosin whereas the residues altered in families with recessive DA only indirectly impair this interaction. Pathological and histological exam of a foot amputated from an affected child revealed complete absence of skeletal muscle (i.e., segmental amyoplasia). To investigate the mechanism for this finding, we generated an animal model for partial MYLPF impairment by knocking out zebrafish mylpfa. The mylpfa mutant had reduced trunk contractile force and complete pectoral fin paralysis, demonstrating that mylpf impairment most severely affects limb movement. mylpfa mutant muscle weakness was most pronounced in an appendicular muscle and was explained by reduced myosin activity and fiber degeneration. Collectively, our findings demonstrate that partial loss of MYLPF function can lead to congenital contractures, likely as a result of degeneration of skeletal muscle in the distal limb.

A recurrent de novo variant in NUSAP1 escapes nonsense-mediated decay and leads to microcephaly, epilepsy, and developmental delay.

NUSAP1 encodes a cell cycle-dependent protein with key roles in mitotic progression, spindle formation, and microtubule stability. Both over- and under-expression of NUSAP1 lead to dysregulation of mitosis and impaired cell proliferation. Through exome sequencing and Matchmaker Exchange, we identified two unrelated individuals with the same recurrent, de novo heterozygous variant (NM_016359.5 c.1209C > A; p.(Tyr403Ter)) in NUSAP1. Both individuals had microcephaly, severe developmental delay, brain abnormalities, and seizures. The gene is predicted to be tolerant of heterozygous loss-of-function mutations, and we show that the mutant transcript escapes nonsense mediated decay, suggesting that the mechanism is likely dominant-negative or toxic gain of function. Single-cell RNA-sequencing of an affected individual's post-mortem brain tissue indicated that the NUSAP1 mutant brain contains all main cell lineages, and that the microcephaly could not be attributed to loss of a specific cell type. We hypothesize that pathogenic variants in NUSAP1 lead to microcephaly possibly by an underlying defect in neural progenitor cells.

Recurrent de novo missense variants in GNB2 can cause syndromic intellectual disability.

PURPOSE: Binding proteins (G-proteins) mediate signalling pathways involved in diverse cellular functions and comprise Gα and Gßγ units. Human diseases have been reported for all five Gß proteins. A de novo missense variant in GNB2 was recently reported in one individual with developmental delay/intellectual disability (DD/ID) and dysmorphism. We aim to confirm GNB2 as a neurodevelopmental disease gene, and elucidate the GNB2-associated neurodevelopmental phenotype in a patient cohort. METHODS: We discovered a GNB2 variant in the index case via exome sequencing and sought individuals with GNB2 variants via international data-sharing initiatives. In silico modelling of the variants was assessed, along with multiple lines of evidence in keeping with American College of Medical Genetics and Genomics guidelines for interpretation of sequence variants. RESULTS: We identified 12 unrelated individuals with five de novo missense variants in GNB2, four of which are recurrent: p.(Ala73Thr), p.(Gly77Arg), p.(Lys89Glu) and p.(Lys89Thr). All individuals have DD/ID with variable dysmorphism and extraneurologic features. The variants are located at the universally conserved shared interface with the Gα subunit, which modelling suggests weaken this interaction. CONCLUSION: Missense variants in GNB2 cause a congenital neurodevelopmental disorder with variable syndromic features, broadening the spectrum of multisystem phenotypes associated with variants in genes encoding G-proteins.

[NEXT-GENERATION SEQUENCING PERFORMED IN PATIENTS RAISING THE SUSPICION OF AN INBORN ERROR OF METABOLISM UNCOVERED A HOMOZYGOUS VARIANT IN YARS1 ALLOWING A NOVEL THERAPEUTIC TRIAL].

INTRODUCTION: Inborn-Errors of Metabolism (IEM) are genetic disorders resulting from mutations in genes encoding proteins involved in biochemical-metabolic pathways. However, some IEMs lack specific biochemical markers. Early incorporation of next-generation-sequencing (NGS) including whole exome sequencing (WES) into the diagnostic algorithm of IEMs herein provided, increases diagnostic accuracy, permits genetic counseling and improves therapeutic options. This is exemplified by diseases affecting aminoacyl-tRNA synthetases (ARSs), enzymes involved in protein translation. Recent studies showed that supplementing amino-acids to cell-culture and patients with ARSs deficiencies resulted in improvement of biochemical and clinical parameters, respectively.

Combining cytogenetic and genomic technologies for deciphering challenging complex chromosomal rearrangements.

Complex chromosomal rearrangements (CCRs), a class of structural variants (SVs) involving more than two chromosome breaks, were classically thought to be extremely rare. As advanced technologies become more available, it has become apparent that CCRs are more common than formerly thought, and are a substantial cause of genetic disorders. We attempted a novel approach for solving the mechanism of challenging CCRs, which involve repetitive sequences, by precisely identifying sequence-level changes and their order. Chromosomal microarray (CMA) and FISH analyses were used for interpretation of SVs detected by whole exome sequencing (WES). Breakpoint junctions were analyzed by Nanopore sequencing, a novel long-read whole genome sequencing tool. A large deletion identified by WES, encompassing the FOXF1 enhancer, was the cause of alveolar capillary dysplasia and respiratory insufficiency, resulting in perinatal death. CMA analysis of the newborn's mother revealed two duplications encompassing the deleted region in the proband, raising our hypothesis that the deletion resulted from the mother's CCR. Breakpoint junctions of complex SVs were determined at the nucleotide level using Nanopore long-read sequencing. According to sequencing results of breakpoint junctions, the CCR in the newborn was considered the consequence of at least one double-strand break during meiosis, and reassembly of DNA fragments by intra-chromosomal homologous recombination. Our comprehensive approach, combining cytogenetics and long-read sequencing, enabled delineation of the exact breakpoints in a challenging CCR, and proposal of a mechanism in which it arises. We suggest applying our integrative approach combining technologies for deciphering future challenging CCRs, enabling risk assessment in families.

Bi-allelic PAGR1 variants are associated with microcephaly and a severe neurodevelopmental disorder: Genetic evidence from two families.

Exome and genome sequencing were used to identify the genetic etiology of a severe neurodevelopmental disorder in two unrelated Ashkenazi Jewish families with three affected individuals. The clinical findings included a prenatal presentation of microcephaly, polyhydramnios and clenched hands while postnatal findings included microcephaly, severe developmental delay, dysmorphism, neurologic deficits, and death in infancy. A shared rare homozygous, missense variant (c.274A > G; p.Ser92Gly, NM_024516.4) was identified in PAGR1, a gene currently not associated with a Mendelian disease. PAGR1 encodes a component of the histone methyltransferase MLL2/MLL3 complex and may function in the DNA damage response pathway. Complete knockout of the murine Pagr1a is embryonic-lethal. Given the available evidence, PAGR1 is a strong candidate gene for a novel autosomal recessive severe syndromic neurodevelopmental disorder.

Expanding the phenotypic spectrum of COLEC10-Related 3MC syndrome: A glimpse into COLEC10-Related 3MC syndrome in the Ashkenazi Jewish population.

Bi-allelic variants in COLEC11 and MASP1 have been associated with 3MC syndrome, a clinical entity made of up four rare autosomal recessive disorders: Carnevale, Mingarelli, Malpuech, and Michels syndromes, characterized by variable expression of facial dysmorphia, cleft lip/palate, postnatal growth deficiency, hearing loss, cognitive impairment, craniosynostosis, radioulnar synostosis, and genital and vesicorenal anomalies. More recently, bi-allelic variants in COLEC10 have been described to be associated with 3MC syndrome. Syndromic features seen in 3MC syndrome are thought to be due to disruption of the chemoattractant properties that influence neural crest cell migration. We identified nine individuals from five families of Ashkenazi Jewish descent with homozygosity of the c.311G > T (p.Gly104Val) variant in COLEC10 and phenotype consistent with 3MC syndrome. Carrier frequency was calculated among 52,278 individuals of Jewish descent. Testing revealed 400 carriers out of 39,750 individuals of Ashkenazi Jewish descent, giving a carrier frequency of 1 in 99 or 1.01%. Molecular protein modeling suggested that the p.Gly104Val substitution alters local conformation. The c.311G > T (p.Gly104Val) variant likely represents a founder variant, and homozygosity is associated with features of 3MC syndrome. 3MC syndrome should be in the differential diagnosis for individuals with short stature, radioulnar synostosis, cleft lip and cleft palate.

The many etiologies of nonimmune hydrops fetalis diagnosed by exome sequencing.

OBJECTIVE: To explain the importance of identifying an etiology for the pathological finding of nonimmune hydrops fetalis (NIHF) and to explore the impact of exome sequencing in recurrent NIHF. In addition, we present two cases of pregnancies affected with recurrent NIHF, in which genetic investigation was advantageous. METHODS: Our study aimed to investigate the genetic background, if available, of all fetuses with NIHF referred to our tertiary medical center from January 2013 to August 2020. We summarized the etiology of NIHF if known, sonographic findings, genetic investigation and the pregnancies' outcomes. RESULTS: We encountered 144 families with NIHF. Genetic investigation was performed by chromosomal microarray analysis (CMA) in 63 (63/144. 44%) fetuses. Seventeen of 63 (27%) had a positive CMA result. In the negative CMA group, 15 (15/46, 33%) opted for exome sequencing, of which seven exomes were positive (47%). Among these, there were four couples with recurrent pregnancies affected by hydrops. Among the remaining 11 exome investigations for non-recurrent hydrops, another three were diagnostic. CONCLUSION: As identifying the etiology of the NIHF is an invaluable tool for the prognosis of the pregnancy, exome sequencing can provide further elucidation of the underlying pathogenesis of NIHF. Thus, genetic investigation should be recommended for cases of NIHF.

A Zebrafish Model for a Rare Genetic Disease Reveals a Conserved Role for FBXL3 in the Circadian Clock System.

The circadian clock, which drives a wide range of bodily rhythms in synchrony with the day-night cycle, is based on a molecular oscillator that ticks with a period of approximately 24 h. Timed proteasomal degradation of clock components is central to the fine-tuning of the oscillator's period. FBXL3 is a protein that functions as a substrate-recognition factor in the E3 ubiquitin ligase complex, and was originally shown in mice to mediate degradation of CRY proteins and thus contribute to the mammalian circadian clock mechanism. By exome sequencing, we have identified a FBXL3 mutation in patients with syndromic developmental delay accompanied by morphological abnormalities and intellectual disability, albeit with a normal sleep pattern. We have investigated the function of FBXL3 in the zebrafish, an excellent model to study both vertebrate development and circadian clock function and, like humans, a diurnal species. Loss of fbxl3a function in zebrafish led to disruption of circadian rhythms of promoter activity and mRNA expression as well as locomotor activity and sleep-wake cycles. However, unlike humans, no morphological effects were evident. These findings point to an evolutionary conserved role for FBXL3 in the circadian clock system across vertebrates and to the acquisition of developmental roles in humans.

Clinical and Functional Study of a De Novo Variant in the PVP Motif of Kv1.1 Channel Associated with Epilepsy, Developmental Delay and Ataxia.

Mutations in the KCNA1 gene, encoding the voltage-gated potassium channel Kv1.1, have been associated with a spectrum of neurological phenotypes, including episodic ataxia type 1 and developmental and epileptic encephalopathy. We have recently identified a de novo variant in KCNA1 in the highly conserved Pro-Val-Pro motif within the pore of the Kv1.1 channel in a girl affected by early onset epilepsy, ataxia and developmental delay. Other mutations causing severe epilepsy are located in Kv1.1 pore domain. The patient was initially treated with a combination of antiepileptic drugs with limited benefit. Finally, seizures and ataxia control were achieved with lacosamide and acetazolamide. The aim of this study was to functionally characterize Kv1.1 mutant channel to provide a genotype-phenotype correlation and discuss therapeutic options for KCNA1-related epilepsy. To this aim, we transfected HEK 293 cells with Kv1.1 or P403A cDNAs and recorded potassium currents through whole-cell patch-clamp. P403A channels showed smaller potassium currents, voltage-dependent activation shifted by +30 mV towards positive potentials and slower kinetics of activation compared with Kv1.1 wild-type. Heteromeric Kv1.1+P403A channels, resembling the condition of the heterozygous patient, confirmed a loss-of-function biophysical phenotype. Overall, the functional characterization of P403A channels correlates with the clinical symptoms of the patient and supports the observation that mutations associated with severe epileptic phenotype cluster in a highly conserved stretch of residues in Kv1.1 pore domain. This study also strengthens the beneficial effect of acetazolamide and sodium channel blockers in KCNA1 channelopathies.

Biallelic deletion in a minimal CAPN15 intron in siblings with a recognizable syndrome of congenital malformations and developmental delay.

Calpainopathies constitute a heterogeneous group of disorders resulting from deficiencies in calpains, calcium-specific proteases that modulate substrates by limited proteolysis. Clinical manifestations depend on tissue-specific expression of the defective calpain and substrate specificity. CAPN15, encoding the Drosophila small optic lobes (sol) homolog, was recently found to cause various eye defects in individuals carrying bi-allelic missense variants. Here we report on two siblings with manifestations reminiscent of Johanson-Blizzard syndrome including failure to thrive, microcephaly, global developmental delay, dysmorphic features, endocrine abnormalities and congenital malformations, in addition to eye abnormalities. Exome sequencing identified a homozygous 47 base-pair deletion in a minimal intron of CAPN15, including the splice donor site. Sequencing of cDNA revealed single exon skipping, resulting in an out-of-frame deletion with a predicted premature termination codon. These findings expand the phenotypic spectrum associated with CAPN15 variants, and suggest that complete loss-of-function is associated with a recognizable syndrome of congenital malformations and developmental delay, overlapping Johanson-Blizzard syndrome and the recently observed brain defects in Capn15 knockout (KO) mice. Moreover, the data highlight the unique opportunity for indel detection in minimal introns.

De novo variant in AMOTL1 in infant with cleft lip and palate, imperforate anus and dysmorphic features.

AMOTL1 belongs to the Motin family of proteins that are involved in organogenesis and tumorigenesis through regulation of cellular migration, tube formation, and angiogenesis. While involvement of all AMOTs in development or suppression of cancers is relatively well described, little is known about the congenital phenotype of pathogenic variants in these genes in humans. Recently, a heterozygous variant in AMOTL1 was published in association with orofacial clefts and cardiac abnormalities in an affected father and his daughter. However, studies in mice did not recapitulate the human phenotype and the case was summarized as inconclusive. We present a female infant with cleft lip and palate, imperforate anus and dysmorphic features, in whom trio exome sequencing revealed a de novo variant in AMOTL1 affecting a highly conserved amino acid (c.479C>T; p.[Pro160Leu]). Bioinformatic predictions and in silico modeling supported pathogenicity. This case reinforces the conjecture regarding the disruptive effect of pathogenic variants in AMOTL1 on organ formation in humans. Studies of additional families will reveal the full phenotypic spectrum associated with this multiple malformation syndrome.

Grandparental genotyping enhances exome variant interpretation.

Trio exome sequencing is a powerful tool in the molecular investigation of monogenic disorders and provides an incremental diagnostic yield over proband-only sequencing, mainly due to the rapid identification of de novo disease-causing variants. However, heterozygous variants inherited from unaffected parents may be inadvertently dismissed, although multiple explanations are available for such scenarios including mosaicism in the parent, incomplete penetrance, imprinting, or skewed X-inactivation. We report three probands, in which a pathogenic or likely pathogenic variant was identified upon exome sequencing, yet was inherited from an unaffected parent. Segregation of the variants (in NOTCH1, PHF6, and SOX10) in the grandparent generation revealed that the variant was de novo in each case. Additionally, one proband had skewed X-inactivation. We discuss the possible genetic mechanism in each case, and urge caution in data interpretation of exome sequencing data. We illustrate the utility of expanding segregation studies to the grandparent generation and demonstrate the impact on exome interpretation strategies, by showing that objective genotype data can overcome subjective parental report of lack of symptoms.

Homozygous variants in MAPRE2 and CDON in individual with skin folds, growth delay, retinal coloboma, and pyloric stenosis.

Cases with multiple molecular diagnoses are challenging to diagnose clinically, yet may be resolved by unbiased exome sequencing analysis. We report an infant with developmental delay, severe growth delay, dysmorphic features, and multiple congenital anomalies including retinal coloboma, congenital pyloric stenosis, and circumferential skin creases. Exome sequencing identified a homozygous missense variant in MAPRE2 and a homozygous stopgain (nonsense) variant in CDON. Variants in MAPRE2, encoding a regulator of microtubule dynamics, lead to congenital symmetric circumferential skin creases type 2, with associated dysmorphism, small growth parameters, and congenital cardiac and genital anomalies. Monoallelic variants in CDON, encoding a coreceptor for sonic hedgehog, have been associated with autosomal dominant pituitary stalk interruption syndrome and holoprosencephaly. Cdon-/- mice have multiple eye defects including coloboma, consistent with the observed human phenotype. Thus, the complex phenotypic presentation of the infant may potentially be attributed to a dual molecular diagnosis. Furthermore, we present CDON as a candidate gene for coloboma formation in addition to the known holoprosencephaly phenotype, and propose to expand the allelic spectrum of CDON to variants associated with autosomal recessive inheritance in addition to dominant inheritance.