Biallelic loss-of-function variants in RABGAP1 cause a novel neurodevelopmental syndrome

PurposeRABGAP1 is a GTPase-activating protein implicated in a variety of cellular and molecular processes, including mitosis, cell migration, vesicular trafficking, and mTOR signaling. There are no known Mendelian diseases caused by variants in RABGAP1.MethodsThrough GeneMatcher, we identified 5 patients from 3 unrelated families with homozygous variants in the RABGAP1 gene found on exome sequencing. We established lymphoblastoid cells lines derived from an affected individual and her parents and performed RNA sequencing and functional studies. Rabgap1 knockout mice were generated and phenotyped.ResultsWe report 5 patients presenting with a common constellation of features, including global developmental delay/intellectual disability, microcephaly, bilateral sensorineural hearing loss, and seizures, as well as overlapping dysmorphic features. Neuroimaging revealed common features, including delayed myelination, white matter volume loss, ventriculomegaly, and thinning of the corpus callosum. Functional analysis of patient cells revealed downregulated mTOR signaling and abnormal localization of early endosomes and lysosomes. Rabgap1 knockout mice exhibited several features in common with the patient cohort, including microcephaly, thinning of the corpus callosum, and ventriculomegaly.ConclusionCollectively, our results provide evidence of a novel neurodevelopmental syndrome caused by biallelic loss-of-function variants in RABGAP1.     

Biallelic variants in PPP1R13L cause paediatric dilated cardiomyopathy

Childhood dilated cardiomyopathy (DCM) is a leading cause of heart failure requiring cardiac transplantation and approximately 5% of cases result in sudden death. Knowledge of the underlying genetic cause can aid prognostication and clinical management and enables accurate recurrence risk counselling for the family. Here we used genomic sequencing to identify the causative genetic variant(s) in families with children affected by severe DCM. In an international collaborative effort facilitated by GeneMatcher, biallelic variants in PPP1R13L were identified in seven children with severe DCM from five unrelated families following exome or genome sequencing and inheritance-based variant filtering. PPP1R13L encodes inhibitor of apoptosis-stimulating protein of p53 protein (iASPP). In addition to roles in apoptosis, iASPP acts as a regulator of desmosomes and has been implicated in inflammatory pathways. DCM presented early (mean: 2 years 10 months; range: 3 months-9 years) and was progressive, resulting in death (n = 3) or transplant (n = 3), with one child currently awaiting transplant. Genomic sequencing technologies are valuable for the identification of novel and emerging candidate genes. Biallelic variants in PPP1R13L were previously reported in a single consanguineous family with paediatric DCM. The identification here of a further five families now provides sufficient evidence to support a robust gene-disease association between PPP1R13L and severe paediatric DCM. The PPP1R13L gene should be included in panel-based genetic testing for paediatric DCM.     

Biallelic loss-of-function variants of GFRA1 cause lethal bilateral renal agenesis

Bilateral renal agenesis belongs to a group of perinatal lethal renal diseases. To date, pathogenic variants in three genes (ITGA8, GREB1L, and FGF20) have been shown to cause renal agenesis in humans. Recently GFRA1 has been linked to a phenotype consistent with a nonsyndromic form of bilateral renal agenesis. GFRA1 encodes a member of the glial cell line-derived neurotrophic factor receptor family of proteins. The receptor on the Wolffian duct regulates ureteric bud outgrowth in developing a functional renal system. We report on four additional affected neonates from a consanguineous family who presented with a similar lethal phenotype whereby whole exome sequencing identified a homozygous deleterious sequence variant in GFRA1 (NM_005264.8:c.628G > T:p.[Gly210Ter]). The current study represents a second confirmation report on the causal association of GFRA1 pathogenic variants with lethal nonsyndromic bilateral renal agenesis in humans.     

Biallelic loss-of-function variants in KCNJ16 presenting with hypokalemic metabolic acidosis

KCNJ16 encodes K_ir5.1 and acts in combination with K_ir4.1, encoded by KCNJ10, to form an inwardly rectifying K⁺ channel expressed at the basolateral membrane of epithelial cells in the distal nephron. This K_ir4.1/K_ir5.1 channel is critical for controlling basolateral membrane potential and K⁺ recycling, the latter coupled to Na-K-ATPase activity, which determines renal Na⁺ handling. Previous work has shown that Kcnj16^−/− mice and SS^{Kcnj16−/−} rats demonstrate hypokalemic, hyperchloremic metabolic acidosis. Here, we present the first report of a patient identified to have biallelic loss-of-function variants in KCNJ16 by whole exome sequencing who presented with chronic metabolic acidosis with exacerbations triggered by minor infections.Subject terms: Paediatric kidney disease, Genetics research     

Biallelic variants in CRIPT cause a Rothmund-Thomson-like syndrome with increased cellular senescence

PurposeRothmund-Thomson syndrome (RTS) is characterized by poikiloderma, sparse hair, small stature, skeletal defects, cancer, and cataracts, resembling features of premature aging. RECQL4 and ANAPC1 are the 2 known disease genes associated with RTS in >70% of cases. We describe RTS-like features in 5 individuals with biallelic variants in CRIPT (OMIM 615789).MethodsTwo newly identified and 4 published individuals with CRIPT variants were systematically compared with those with RTS using clinical data, computational analysis of photographs, histologic analysis of skin, and cellular studies on fibroblasts.ResultsAll CRIPT individuals fulfilled the diagnostic criteria for RTS and additionally had neurodevelopmental delay and seizures. Using computational gestalt analysis, CRIPT individuals showed greatest facial similarity with individuals with RTS. Skin biopsies revealed a high expression of senescence markers (p53/p16/p21) and the senescence-associated ß-galactosidase activity was elevated in CRIPT-deficient fibroblasts. RECQL4- and CRIPT-deficient fibroblasts showed an unremarkable mitotic progression and unremarkable number of mitotic errors and no or only mild sensitivity to genotoxic stress by ionizing radiation, mitomycin C, hydroxyurea, etoposide, and potassium bromate.ConclusionCRIPT causes an RTS-like syndrome associated with neurodevelopmental delay and epilepsy. At the cellular level, RECQL4- and CRIPT-deficient cells display increased senescence, suggesting shared molecular mechanisms leading to the clinical phenotypes.     

Biallelic variants in four genes underlying recessive osteogenesis imperfecta

Osteogenesis imperfecta (OI) is an inherited heterogeneous rare skeletal disorder characterized by increased bone fragility and low bone mass. The disorder mostly segregates in an autosomal dominant manner. However, several rare autosomal recessive and X-linked forms, caused by mutations in 18 different genes, have also been described in the literature.Here, we present five consanguineous families segregating OI in an autosomal recessive pattern. Affected individuals in the five families presented severe forms of skeletal deformities. It included frequent bone fractures with abnormal healing, short stature, facial dysmorphism, osteopenia, joint laxity, and severe scoliosis. In order to search for the causative variants, DNA of at least one affected individual in three families (A-C) were subjected to whole exome sequencing (WES). In two other families (D-E), linkage analysis using highly polymorphic microsatellite markers was followed by Sanger sequencing. Sequence analysis revealed two novels and three previously reported disease-causing variants. The two novel homozygous variants including [c.824G > A; p.(Cys275Tyr)] in the SP7 gene and [c.397C > T, p.(Gln133*)] in the SERPINF1 gene were identified in families A and B, respectively. The three previously reported homozygous variants including [c.497G > A; p.(Arg166His)] in the SPARC gene, (c.359-3C > G; intron 2) and [c.677C > T; p.(Ser226Leu)] in the WNT1 gene were identified in family C, D, and E.In conclusion, our findings provided additional evidence of involvement of homozygous sequence variants in the SP7, SERPINF1, SPARC and WNT1 genes causing severe OI. It also highlights the importance of extensive genetic investigations to search for the culprit gene in each case of skeletal deformity.     

Biallelic variants in CEP164 cause a motile ciliopathy-like syndrome

Ciliopathies may be classed as primary or motile depending on the underlying ciliary defect and are usually considered distinct clinical entities. Primary ciliopathies are associated with multisystem syndromes typically affecting the brain, kidney, and eye, as well as other organ systems such as the liver, skeleton, auditory system, and metabolism. Motile ciliopathies are a heterogenous group of disorders with defects in specialised motile ciliated tissues found within the lung, brain, and reproductive system, and are associated with primary ciliary dyskinesia, bronchiectasis, infertility and rarely hydrocephalus. Primary and motile cilia share defined core ultra-structures with an overlapping proteome, and human disease phenotypes can reflect both primary and motile ciliopathies. CEP164 encodes a centrosomal distal appendage protein vital for primary ciliogenesis. Human CEP164 mutations are typically described in patients with nephronophthisis-related primary ciliopathies but have also been implicated in motile ciliary dysfunction. Here we describe a patient with an atypical motile ciliopathy phenotype and biallelic CEP164 variants. This work provides further evidence that CEP164 mutations can contribute to both primary and motile ciliopathy syndromes, supporting their functional and clinical overlap, and informs the investigation and management of CEP164 ciliopathy patients.     

Biallelic variants in the SLC13A1 sulfate transporter gene cause hyposulfatemia with a mild spondylo-epi-metaphyseal dysplasia

Sulfate is the fourth most abundant anion in human plasma but is not measured in clinical practice and little is known about the consequences of sulfate deficiency. Nevertheless, sulfation plays an essential role in the modulation of numerous compounds, including proteoglycans and steroids. We report the first patient with a homozygous loss-of-function variant in the SLC13A1 gene, encoding a renal and intestinal sulfate transporter, which is essential for maintaining plasma sulfate levels. The homozygous (Arg12Ter) variant in SLC13A1 was found by exome sequencing performed in a patient with unexplained skeletal dysplasia. The main clinical features were enlargement of joints and spondylo-epi-metaphyseal radiological abnormalities in early childhood, which improved with age. In addition, autistic features were noted. We found profound hyposulfatemia due to complete loss of renal sulfate reabsorption. Cholesterol sulfate was reduced. Intravenous N-acetylcysteine administration temporarily restored plasma sulfate levels. We conclude that loss of the SLC13A1 gene leads to profound hypersulfaturia and hyposulfatemia, which is mainly associated with abnormal skeletal development, possibly predisposing to degenerative bone and joint disease. The diagnosis might be easily missed and more frequent.     

Biallelic loss-of-function variants in WDR11 are associated with microcephaly and intellectual disability

Heterozygous missense variants in the WD repeat domain 11 (WDR11) gene are associated with hypogonadotropic hypogonadism in humans. In contrast, knockout of both alleles of Wdr11 in mice results in a more severe phenotype with growth and developmental delay, features of holoprosencephaly, heart defects and reproductive disorders. Similar developmental defects known to be associated with aberrant hedgehog signaling and ciliogenesis have been found in zebrafish after Wdr11 knockdown. We here report biallelic loss-of-function variants in the WDR11 gene in six patients from three independent families with intellectual disability, microcephaly and short stature. The findings suggest that biallelic WDR11 variants in humans result in an overlapping but milder phenotype compared to Wdr11-deficient animals. However, the observed human phenotype differs significantly from dominantly inherited variants leading to hypogonadotropic hypogonadism, suggesting that recessive WDR11 variants result in a clinically distinct entity.Subject terms: Disease genetics, Diseases     

Biallelic variants in KARS1 are associated with neurodevelopmental disorders and hearing loss recapitulated by the knockout zebrafish

PurposePathogenic variants in Lysyl-tRNA synthetase 1 (KARS1) have increasingly been recognized as a cause of early-onset complex neurological phenotypes. To advance the timely diagnosis of KARS1-related disorders, we sought to delineate its phenotype and generate a disease model to understand its function in vivo.MethodsThrough international collaboration, we identified 22 affected individuals from 16 unrelated families harboring biallelic likely pathogenic or pathogenic in KARS1 variants. Sequencing approaches ranged from disease-specific panels to genome sequencing. We generated loss-of-function alleles in zebrafish.ResultsWe identify ten new and four known biallelic missense variants in KARS1 presenting with a moderate-to-severe developmental delay, progressive neurological and neurosensory abnormalities, and variable white matter involvement. We describe novel KARS1-associated signs such as autism, hyperactive behavior, pontine hypoplasia, and cerebellar atrophy with prevalent vermian involvement. Loss of kars1 leads to upregulation of p53, tissue-specific apoptosis, and downregulation of neurodevelopmental related genes, recapitulating key tissue-specific disease phenotypes of patients. Inhibition of p53 rescued several defects of kars1^−/− knockouts.ConclusionOur work delineates the clinical spectrum associated with KARS1 defects and provides a novel animal model for KARS1-related human diseases revealing p53 signaling components as potential therapeutic targets.     

Biallelic loss of function variants in PPP1R21 cause a neurodevelopmental syndrome with impaired endocytic function

Next‐generation sequencing (NGS) has been instrumental in solving the genetic basis of rare inherited diseases, especially neurodevelopmental syndromes. However, functional workup is essential for precise phenotype definition and to understand the underlying disease mechanisms. Using whole exome (WES) and whole genome sequencing (WGS) in four independent families with hypotonia, neurodevelopmental delay, facial dysmorphism, loss of white matter, and thinning of the corpus callosum, we identified four previously unreported homozygous truncating PPP1R21 alleles: c.347delT p.(Ile116Lysfs*25), c.2170_2171insGGTA p.(Ile724Argfs*8), c.1607dupT p.(Leu536Phefs*7), c.2063delA p.(Lys688Serfs*26) and found that PPP1R21 was absent in fibroblasts of an affected individual, supporting the allele's loss of function effect. PPP1R21 function had not been studied except that a large scale affinity proteomics approach suggested an interaction with PIBF1 defective in Joubert syndrome. Our co‐immunoprecipitation studies did not confirm this but in contrast defined the localization of PPP1R21 to the early endosome. Consistent with the subcellular expression pattern and the clinical phenotype exhibiting features of storage diseases, we found patient fibroblasts exhibited a delay in clearance of transferrin‐488 while uptake was normal. In summary, we delineate a novel neurodevelopmental syndrome caused by biallelic PPP1R21 loss of function variants, and suggest a role of PPP1R21 within the endosomal sorting process or endosome maturation pathway.     

Biallelic and monoallelic variants in PLXNA1 are implicated in a novel neurodevelopmental disorder with variable cerebral and eye anomalies

PurposeTo investigate the effect of PLXNA1 variants on the phenotype of patients with autosomal dominant and recessive inheritance patterns and to functionally characterize the zebrafish homologs plxna1a and plxna1b during development.MethodsWe assembled ten patients from seven families with biallelic or de novo PLXNA1 variants. We describe genotype–phenotype correlations, investigated the variants by structural modeling, and used Morpholino knockdown experiments in zebrafish to characterize the embryonic role of plxna1a and plxna1b.ResultsShared phenotypic features among patients include global developmental delay (9/10), brain anomalies (6/10), and eye anomalies (7/10). Notably, seizures were predominantly reported in patients with monoallelic variants. Structural modeling of missense variants in PLXNA1 suggests distortion in the native protein. Our zebrafish studies enforce an embryonic role of plxna1a and plxna1b in the development of the central nervous system and the eye.ConclusionWe propose that different biallelic and monoallelic variants in PLXNA1 result in a novel neurodevelopmental syndrome mainly comprising developmental delay, brain, and eye anomalies. We hypothesize that biallelic variants in the extracellular Plexin-A1 domains lead to impaired dimerization or lack of receptor molecules, whereas monoallelic variants in the intracellular Plexin-A1 domains might impair downstream signaling through a dominant-negative effect.     

Biallelic variants in NSUN6 cause an autosomal recessive neurodevelopmental disorder

Purpose5-methylcytosine RNA modifications are driven by NSUN methyltransferases. Although variants in NSUN2 and NSUN3 were associated with neurodevelopmental diseases, the physiological role of NSUN6 modifications on transfer RNAs and messenger RNAs remained elusive.MethodsWe combined exome sequencing of consanguineous families with functional characterization to identify a new neurodevelopmental disorder gene.ResultsWe identified 3 unrelated consanguineous families with deleterious homozygous variants in NSUN6. Two of these variants are predicted to be loss-of-function. One maps to the first exon and is predicted to lead to the absence of NSUN6 via nonsense-mediated decay, whereas we showed that the other maps to the last exon and encodes a protein that does not fold correctly. Likewise, we demonstrated that the missense variant identified in the third family has lost its enzymatic activity and is unable to bind the methyl donor S-adenosyl-L-methionine. The affected individuals present with developmental delay, intellectual disability, motor delay, and behavioral anomalies. Homozygous ablation of the NSUN6 ortholog in Drosophila led to locomotion and learning impairment.ConclusionOur data provide evidence that biallelic pathogenic variants in NSUN6 cause one form of autosomal recessive intellectual disability, establishing another link between RNA modification and cognition.     

Biallelic variants in LINGO1 are associated with autosomal recessive intellectual disability,microcephaly, speech and motor delay

PurposeTo elucidate the novel molecular cause in two unrelated consanguineous families with autosomal recessive intellectual disability.MethodsA combination of homozygosity mapping and exome sequencing was used to locate the plausible genetic defect in family F162, while only exome sequencing was followed in the family PKMR65. The protein 3D structure was visualized with the University of California–San Francisco Chimera software.ResultsAll five patients from both families presented with severe intellectual disability, aggressive behavior, and speech and motor delay. Four of the five patients had microcephaly. We identified homozygous missense variants in LINGO1, p.(Arg290His) in family F162 and p.(Tyr288Cys) in family PKMR65. Both variants were predicted to be pathogenic, and segregated with the phenotype in the respective families. Molecular modeling of LINGO1 suggests that both variants interfere with the glycosylation of the protein.ConclusionLINGO1 is a transmembrane receptor, predominantly found in the central nervous system. Published loss-of-function studies in mouse and zebrafish have established a crucial role of LINGO1 in normal neuronal development and central nervous system myelination by negatively regulating oligodendrocyte differentiation and neuronal survival. Taken together, our results indicate that biallelic LINGO1 missense variants cause autosomal recessive intellectual disability in humans.     

Biallelic TMEM251 variants in patients with severe skeletal dysplasia and extreme short stature

Skeletal dysplasias are a heterogeneous group of disorders ranging from mild to lethal skeletal defects. We investigated two unrelated families with individuals presenting with a severe skeletal disorder. In family NMD02, affected individuals had a dysostosis multiplex‐like skeletal dysplasia and severe short stature (<?8.5 SD). They manifested increasingly coarse facial features, protruding abdomens, and progressive skeletal changes, reminiscent of mucopolysaccharidosis. The patients gradually lost mobility and the two oldest affected individuals died in their twenties. The affected child in family ID01 had coarse facial features and severe skeletal dysplasia with clinical features similar to mucopolysaccharidosis. She had short stature, craniosynostosis, kyphoscoliosis, and hip‐joint subluxation. She died at the age of 5 years. Whole‐exome sequencing identified two homozygous variants c.133C>T; p.(Arg45Trp) and c.215dupA; p.(Tyr72Ter), respectively, in the two families, affecting an evolutionary conserved gene TMEM251 (NM_001098621.1). Immunofluorescence and confocal studies using human osteosarcoma cells indicated that TMEM251 is localized to the Golgi complex. However, p.Arg45Trp mutant TMEM251 protein was targeted less efficiently and the localization was punctate. Tmem251 knockdown by small interfering RNA induced dedifferentiation of rat primary chondrocytes. Our work implicates TMEM251 in the pathogenesis of a novel disorder and suggests its potential function in chondrocyte differentiation.     

Biallelic TMEM260 variants cause truncus arteriosus,with or without renal defects

Only two families have been reported with biallelic TMEM260 variants segregating with structural heart defects and renal anomalies syndrome (SHDRA). With a combination of genome, exome sequencing and RNA studies, we identified eight individuals from five families with biallelic TMEM260 variants. Variants included one multi-exon deletion, four nonsense/frameshifts, two splicing changes and one missense change. Together with the published cases, analysis of clinical data revealed ventricular septal defects (12/12), mostly secondary to truncus arteriosus (10/12), elevated creatinine levels (6/12), horse-shoe kidneys (1/12) and renal cysts (1/12) in patients. Three pregnancies were terminated on detection of severe congenital anomalies. Six patients died between the ages of 6 weeks and 5 years. Using a range of stringencies, carrier frequency for SHDRA was estimated at 0.0007–0.007 across ancestries. In conclusion, this study confirms the genetic basis of SHDRA, expands its known mutational spectrum and clarifies its clinical features. We demonstrate that SHDRA is a severe condition associated with substantial mortality in early childhood and characterised by congenital cardiac malformations with a variable renal phenotype.     

Heterozygous variants in MYBPC1 are associated with an expanded neuromuscular phenotype beyond arthrogryposis

Encoding the slow skeletal muscle isoform of myosin binding protein‐C, MYBPC1 is associated with autosomal dominant and recessive forms of arthrogryposis. The authors describe a novel association for MYBPC1 in four patients from three independent families with skeletal muscle weakness, myogenic tremors, and hypotonia with gradual clinical improvement. The patients carried one of two de novo heterozygous variants in MYBPC1, with the p.Leu263Arg variant seen in three individuals and the p.Leu259Pro variant in one individual. Both variants are absent from controls, well conserved across vertebrate species, predicted to be damaging, and located in the M‐motif. Protein modeling studies suggested that the p.Leu263Arg variant affects the stability of the M‐motif, whereas the p.Leu259Pro variant alters its structure. In vitro biochemical and kinetic studies demonstrated that the p.Leu263Arg variant results in decreased binding of the M‐motif to myosin, which likely impairs the formation of actomyosin cross‐bridges during muscle contraction. Collectively, our data substantiate that damaging variants in MYBPC1 are associated with a new form of an early‐onset myopathy with tremor, which is a defining and consistent characteristic in all affected individuals, with no contractures. Recognition of this expanded myopathic phenotype can enable identification of individuals with MYBPC1 variants without arthrogryposis.     

Biallelic variants in NUDCD2 associated with a multiple malformation syndrome with cholestasis and renal failure

NudC-like protein 2 (NUDCD2) is a 4-exon protein-coding gene at 5q34. The protein appears to act in concert with other genes regulating cell migration and microtubule extension. Early studies in model organisms show associations with LIS1, HERC2, and cohesin subunits via a co-chaperone function with Heat shock protein 90 (Hsp90). It is a candidate gene for human pathology. We present two unrelated patients with biallelic variants in NUDCD2. Their phenotypes comprise similar dysmorphic facies, midline brain hypoplasia, hypothyroidism, pulmonary and aortic valve stenosis, severe dysfunction of the liver and kidneys, profound hypotonia, and early death. The cellular analysis demonstrates the absence of the NUDCD2 protein in fibroblasts of one patient with biallelic loss-of-function variants. The data suggest that NUDCD2 deficiency causes this recognizable syndrome that has features of a ciliopathy with additional complications.