The Pex1-G844D mouse: A model for mild human Zellweger spectrum disorder

Zellweger spectrum disorder (ZSD) is a disease continuum that results from inherited defects in PEX genes essential for normal peroxisome assembly. These autosomal recessive disorders impact brain development and also cause postnatal liver, adrenal, and kidney dysfunction, as well as loss of vision and hearing. The hypomorphic PEX1-G843D missense allele, observed in approximately 30% of ZSD patients, is associated with milder clinical and biochemical phenotypes, with some homozygous individuals surviving into early adulthood. Nonetheless, affected children with the PEX1-G843D allele have intellectual disability, failure to thrive, and significant sensory deficits. To enhance our ability to test candidate therapies that improve human PEX1-G843D function, we created the novel Pex1-G844D knock-in mouse model that represents the murine equivalent of the common human mutation. We show that Pex1-G844D homozygous mice recapitulate many classic features of mild ZSD cases, including growth retardation and fatty livers with cholestasis. In addition, electrophysiology, histology, and gene expression studies provide evidence that these animals develop a retinopathy similar to that observed in human patients, with evidence of cone photoreceptor cell death. Similar to skin fibroblasts obtained from ZSD patients with a PEX1-G843D allele, we demonstrate that murine cells homozygous for the Pex1-G844D allele respond to chaperone-like compounds, which normalizes peroxisomal β-oxidation. Thus, the Pex1-G844D mouse provides a powerful model system for testing candidate therapies that address the most common genetic cause of ZSD. In addition, this murine model will enhance studies focused on mechanisms of pathogenesis.     

A homozygous Gly470Ala variant in PEX6 causes severe Zellweger spectrum disorder

Zellweger spectrum disorder (ZSD) is a group of autosomal recessive disorders caused by biallelic pathogenic variants in any one of the 13 PEX genes essential for peroxisomal biogenesis. We report a cohort of nine infants who presented at birth with severe neonatal features suggestive of ZSD and found to be homozygous for a variant in PEX6 (NM_000287.4:c.1409G > C[p.Gly470Ala]). All were of Mixtec ancestry and identified by the California Newborn Screening (NBS) Program to have elevated C26:0-lysophosphatidylcholine but no reportable variants in ABCD1. The clinical and biochemical features of this cohort are described within. Gly470Ala may represent a founder variant in the Mixtec population of Central California. ZSD should be considered in patients who present at birth with severe hypotonia and enlarged fontanelles, especially in the setting of an abnormal NBS, Mixtec ancestry, or family history of infant death. There is a need to further characterize the natural history of ZSD, the Gly470Ala variant, and expand upon possible genotype–phenotype correlations.     

Spectrum of PEX6 mutations in Zellweger syndrome spectrum patients

The autosomal recessive Zellweger syndrome spectrum (ZSS) disorders comprise a main subgroup of the peroxisome biogenesis disorders. The ZSS disorders can be caused by mutations in any of 12 different currently identified PEX genes resulting in severe, often lethal, multi‐systemic disorders. Defects in the PEX6 gene are the second most common cause for ZSS disorders. The encoded protein PEX6 belongs to the AAA ATPase family and contains two AAA cassettes and an AAA protein family signature. The PEX6 gene consists of 17 exons and previously mutations in the PEX6 gene were found to be scattered over all exons. We developed a post‐PCR high‐resolution melting (HRM) curve assay to scan the PEX6 gene for potential sequence variations followed by selective sequencing to identify these. We analyzed the PEX6 genes of 75 patients assigned to the PEX6 complementation group. We identified a total of 77 different mutations of which 47 mutations have not been reported previously, and 14 polymorphic variants. © 2009 Wiley‐Liss, Inc.     

Carrier frequency estimation of Zellweger spectrum disorder using ExAC database and bioinformatics tools

PurposeWe aimed to estimate the carrier frequency of Zellweger spectrum disorder (ZSD), a rare autosomal recessive disease, and the associated disease incidence based on data from the Exome Aggregation Consortium (ExAC) of approximately 60,000 individuals.MethodsWe obtained variants from ExAC in 13 PEX genes associated with ZSD. Potentially pathogenic missense variants were identified with computational variant analysis tools according to three stringency levels. Using variants classified as potentially pathogenic, we estimated the carrier frequency and the associated incidence for the entire ExAC population and its subpopulations. We also evaluated variants based on pathogenicity criteria for sequence variant interpretation outlined by the American College of Medical Genetics and Genomics (ACMG) and calculated the carrier frequency and incidence based on those variants.ResultsThe bioinformatically estimated incidence rate of ZSD in the ExAC population is 1 in 83,841 using our least stringent pathogenicity cutoff. Under clinical guidelines outlined by ACMG, the estimated incidence is 1 in 3,275,751 births.ConclusionWe outlined a process for estimating the ZSD disease carrier frequency and incidence in a large consortium using bioinformatics tools. Our results are close to current newborn screening estimates in New York of 1 in 90,000 births, estimated from 1.08 million screenings.     

Heterozygous UCHL1 loss-of-function variants cause a neurodegenerative disorder with spasticity,ataxia, neuropathy,and optic atrophy

PurposeBiallelic variants in UCHL1 have been associated with a progressive early-onset neurodegenerative disorder, autosomal recessive spastic paraplegia type 79. In this study, we investigated heterozygous UCHL1 variants on the basis of results from cohort-based burden analyses.MethodsGene-burden analyses were performed on exome and genome data of independent cohorts of patients with hereditary ataxia and spastic paraplegia from Germany and the United Kingdom in a total of 3169 patients and 33,141 controls. Clinical data of affected individuals and additional independent families were collected and evaluated. Patients’ fibroblasts were used to perform mass spectrometry-based proteomics.ResultsUCHL1 was prioritized in both independent cohorts as a candidate gene for an autosomal dominant disorder. We identified a total of 34 cases from 18 unrelated families, carrying 13 heterozygous loss-of-function variants (15 families) and an inframe insertion (3 families). Affected individuals mainly presented with spasticity (24/31), ataxia (28/31), neuropathy (11/21), and optic atrophy (9/17). The mass spectrometry-based proteomics showed approximately 50% reduction of UCHL1 expression in patients’ fibroblasts.ConclusionOur bioinformatic analysis, in-depth clinical and genetic workup, and functional studies established haploinsufficiency of UCHL1 as a novel disease mechanism in spastic ataxia.     

Mutations in PEX10 is the cause of Zellweger peroxisome deficiency syndrome of complementation group B

Peroxisome biogenesis disorders (PBD), such as Zellweger syndrome, are autosomal recessive diseases caused by a deficiency in peroxisome assembly as well as a malfunction of the peroxisomes, where at least 10 genotypes have been reported. We have isolated a human PEX10 cDNA (HsPEX10) by an expressed sequence tag homology search on a human DNA database using yeast PEX10 from Hansenula polymorpha, followed by screening of a human liver cDNA library. This cDNA encodes a peroxisomal protein (a peroxin Pex10p) comprising 326 amino acids, with two putative transmembrane segments and a C3HC4zinc finger RING motif. Both the N- and C-terminal regions of Pex10p are exposed to the cytosol, as assessed by an expression study of epitope-tagged Pex10p. HsPEX10 expression morphologically and biochemically restored peroxisome biogenesis in fibroblasts from Zellweger patients of complementation group B in Japan (complementation group VII in the USA). One patient (PBDB-01) possessed a homozygous, inactivating mutation, a 2 bp deletion immediately upstream of the RING motif, which resulted in a frameshift, altering 65 amino acids from the normal. This implies that the C-terminal part, including the RING finger, is required for biological function of Pex10p. PEX10 cDNA derived from patient PBDB-01 was defective in peroxisome-restoring activity when expressed in patient fibroblasts. These results demonstrate that mutation in PEX10 is the genetic cause of complementation group B PBD.      

High incidence of hyperoxaluria in generalized peroxisomal disorders

The Zellweger spectrum disorders (ZSDs) are characterized by a generalized loss of peroxisomal functions caused by deficient peroxisomal assembly. Clinical presentation and survival are heterogeneous. Although most peroxisomal enzymes are unstable in the cytosol of peroxisome-deficient cells of ZSD patients, a few enzymes remain stable among which alanine:glyoxylate aminotransferase (AGT). Its deficiency causes primary hyperoxaluria type 1 (PH1, MIM 259900), an inborn error of glyoxylate metabolism characterized by hyperoxaluria, nephrocalcinosis, and renal insufficiency. Despite the normal level of AGT activity in ZSD patients, hyperoxaluria has been reported in several ZSD patients. We observed the unexpected occurrence of renal stones in a cohort of ZSD patients. This led us to perform a study in this cohort to determine the prevalence of hyperoxaluria in ZSDs and to find clinically relevant clues that correlate with the urinary oxalate load. We reviewed medical charts of 31 Dutch ZSD patients with prolonged survival (>1 year). Urinary oxalate excretion was assessed in 23 and glycolate in 22 patients. Hyperoxaluria was present in 19 (83%), and hyperglycolic aciduria in 14 (64%). Pyridoxine treatment in six patients did not reduce the oxalate excretion as in some PH1 patients. Renal involvement with urolithiasis and nephrocalcinosis was present in five of which one developed end-stage renal disease. The presence of hyperoxaluria, potentially leading to severe renal involvement, was statistically significant correlated with the severity of neurological dysfunction. ZSD patients should be screened by urinalysis for hyperoxaluria and renal ultrasound for nephrocalcinosis in order to take timely measures to prevent renal insufficiency.     

A novel PEX12 mutation identified as the cause of a peroxisomal biogenesis disorder with mild clinical phenotype, mild biochemical abnormalities in fibroblasts and a mosaic catalase immunofluorescence pattern, even at 40 °C

Mutations in 12 different PEX genes can cause a generalized peroxisomal biogenesis disorder with clinical phenotypes ranging from Zellweger syndrome to infantile Refsum disease. To identify the specific PEX gene to be sequenced, complementation analysis is first performed in fibroblasts using catalase immunofluorescence. A patient with a relatively mild phenotype of infantile cholestasis, hypotonia and motor delay had elevated plasma very long-chain fatty acids and bile acid precursors, but fibroblast studies revealed normal or only mildly abnormal peroxisomal parameters and mosaic catalase immunofluorescence. This mosaicism persisted even when the incubation temperature was increased from 37 °C to 40 °C, a maneuver previously shown to abolish mosaicism by exacerbating peroxisomal dysfunction. As mosaicism precludes complementation analysis, a candidate gene approach was employed. After PEX1 sequencing was unrewarding, PEX12 sequencing revealed homozygosity for a novel c.102A>T (p.R34S) missense mutation affecting a partially conserved residue in the N-terminal region important for localization to peroxisomes. Transfection of patient fibroblasts with wild-type PEX12 cDNA confirmed that a PEX12 defect was the basis for the PBD. Homozygosity for c.102A>T was identified in a second patient of similar ethnic origin also presenting with a mild phenotype. PEX12 is a highly probable candidate gene for direct sequencing in the context of a mild clinical phenotype with mosaicism and minimally abnormal peroxisomal parameters in fibroblasts.     

Genomic structure and identification of 11 novel mutations of the PEX6 (peroxisome assembly factor‐2) gene in patients with peroxisome biogenesis disorders

The PEX6 (peroxisome assembly factor‐2, PAF‐2) gene encodes a member of the AAA protein (ATPases associated with diverse cellular activities) family and restores peroxisome assembly in fibroblasts from peroxisome biogenesis disorder patients belonging to complementation group C (group 4 in the United States). We have now clarified the genomic DNA structure of human PEX6 and identified mutations in patients from various ethnic groups. The human PEX6 gene consists of 17 exons and 16 introns, spanning about 14kb. The largest exon, exon 1, has at least 952 bp nucleotides. Eleven novel mutations (18 alleles) were identified by direct sequencing of the PEX6 cDNA from 10 patients. All these mutations have been confirmed in the corresponding genomic DNA. There was no common mutation, but an exon skip was identified in two unrelated Japanese patients. Most of the mutations led to premature termination or large deletions of the PEX6 protein and resulted in the most severe peroxisome biogenesis disorder phenotype of Zellweger syndrome. A patient with an atypical Zellweger syndrome had a missense mutation that was shown to disrupt the cell's ability to form peroxisomes. This mutation analysis will aid in understanding the functions of the PEX6 protein in peroxisomal biogenesis. Hum Mutat 13:487–496, 1999. © 1999 Wiley‐Liss, Inc.     

De novo retinoic acid receptor beta (RARB) variant associated with microphthalmia and dystonia

BackgroundDefinition of the individual genotypes that cause a Mendelian phenotype is of great importance both to clinical diagnostics and disease characterization. Heterozygous de novo gain-of-function missense variants in RARB are associated with syndromic microphthalmia 12 (MCOPS12), a developmental disorder characterized by eye malformations and variable involvement of other organs. A subset of patients were described with poorly delineated movement disorders. Additionally, RARB bi-allelic loss-of-function variants, inherited from asymptomatic heterozygous carrier parents, have been found in a recessive family with four MCOPS12-affected members.Patient/methodsWe used trio whole-exome sequencing to explore the molecular basis of disease in an individual with congenital eye abnormality and movement disorder. All patients with reported RARB variants were reviewed.ResultsWe report on identification of a heterozygous de novo RARB nonsense variant in a girl with microphthalmia and progressive generalized dystonia. Public database entries indicate that the de novo variant is recurrently present in clinically affected subjects but a literature report has not yet been available.ConclusionsWe provide the first detailed evidence for a role of dominant RARB truncating alterations in congenital eye-brain disease, expanding the spectrum of MCOPS12-associated mutations. Considered together with the published family with bi-allelic variants, the data suggest manifestation and non-manifestation of disease in relation to almost identical RARB loss-of-function variations, an apparent paradox that is seen in a growing number of human genetic conditions associated with both recessive and dominant inheritance patterns.     

Mutations in the mitochondrial complex I assembly factor NDUFAF6 cause isolated bilateral striatal necrosis and progressive dystonia in childhood

AimTo perform a deep phenotype characterisation in a pedigree of 3 siblings with Leigh syndrome and compound heterozygous NDUFAF6 mutations.MethodA multi-gene panel of childhood-onset basal ganglia neurodegeneration inherited conditions was analysed followed by functional studies in fibroblasts.ResultsThree siblings developed gait dystonia in infancy followed by rapid progression to generalised dystonia and psychomotor regression. Brain magnetic resonance showed symmetric and bilateral cytotoxic lesions in the putamen and proliferation of the lenticular-striate arteries, latter spreading to the caudate and progressing to cavitation and volume loss. We identified a frameshift novel change (c.554_558delTTCTT; p.Tyr187AsnfsTer65) and a pathogenic missense change (c.371T>C; p.Ile124Thr) in the NDUFAF6 gene, which segregated with an autosomal recessive inheritance within the family. Patient mutations were associated with the absence of the NDUFAF6 protein and reduced activity and assembly of mature complex I in fibroblasts. By functional complementation assay, the mutant phenotype was rescued by the canonical version of the NDUFAF6. A literature review of 14 NDUFAF6 patients showed a consistent phenotype of an early childhood insidious onset neurological regression with prominent dystonia associated with basal ganglia degeneration and long survival.InterpretationNDUFAF6-related Leigh syndrome is a relevant cause of childhood onset dystonia and isolated bilateral striatal necrosis. By genetic complementation, we could demonstrate the pathogenicity of novel genetic variants in NDUFAF6.     

Zellweger syndrome resulting from maternal isodisomy of chromosome 1

Zellweger syndrome (ZS) is an autosomal recessive peroxisomal disorder that results from mutations in one of the peroxisome biogenesis (PEX) genes. This is the first patient reported with uniparental disomy (UPD) resulting in ZS, in this case maternal isodisomy of chromosome 1 involving reduction to homoallelism of a frameshift mutation within PEX 10. Other reported cases of UPD1, and evidence for the imprinting of genes on chromosome 1, are reviewed. The molecular findings in this patient have important implications for molecular testing and genetic counseling in ZS.     

Molecular autopsy by trio exome sequencing (ES) and postmortem examination in fetuses and neonates with prenatally identified structural anomalies

PurposeTo determine the diagnostic yield of combined exome sequencing (ES) and autopsy in fetuses/neonates with prenatally identified structural anomalies resulting in termination of pregnancy, intrauterine, neonatal, or early infant death.MethodsES was undertaken in 27 proband/parent trios following full autopsy. Candidate pathogenic variants were classified by a multidisciplinary clinical review panel using American College of Medical Genetics and Genomics (ACMG) guidelines.ResultsA genetic diagnosis was established in ten cases (37%). Pathogenic/likely pathogenic variants were identified in nine different genes including four de novo autosomal dominant, three homozygous autosomal recessive, two compound heterozygous autosomal recessive, and one X-linked. KMT2D variants (associated with Kabuki syndrome postnatally) occurred in two cases. Pathogenic variants were identified in 5/13 (38%) cases with multisystem anomalies, in 2/4 (50%) cases with fetal akinesia deformation sequence, and in 1/4 (25%) cases each with cardiac and brain anomalies and hydrops fetalis. No pathogenic variants were detected in fetuses with genitourinary (1), skeletal (1), or abdominal (1) abnormalities.ConclusionThis cohort demonstrates the clinical utility of molecular autopsy with ES to identify an underlying genetic cause in structurally abnormal fetuses/neonates. These molecular findings provided parents with an explanation of the developmental abnormality, delineated the recurrence risks, and assisted the management of subsequent pregnancies.     

Heterozygous variants in MYH10 associated with neurodevelopmental disorders and congenital anomalies with evidence for primary cilia-dependent defects in Hedgehog signaling

PurposeNonmuscle myosin II complexes are master regulators of actin dynamics that play essential roles during embryogenesis with vertebrates possessing 3 nonmuscle myosin II heavy chain genes, MYH9, MYH10, and MYH14. As opposed to MYH9 and MYH14, no recognizable disorder has been associated with MYH10. We sought to define the clinical characteristics and molecular mechanism of a novel autosomal dominant disorder related to MYH10.MethodsAn international collaboration identified the patient cohort. CAS9-mediated knockout cell models were used to explore the mechanism of disease pathogenesis.ResultsWe identified a cohort of 16 individuals with heterozygous MYH10 variants presenting with a broad spectrum of neurodevelopmental disorders and variable congenital anomalies that affect most organ systems and were recapitulated in animal models of altered MYH10 activity. Variants were typically de novo missense changes with clustering observed in the motor domain. MYH10 knockout cells showed defects in primary ciliogenesis and reduced ciliary length with impaired Hedgehog signaling. MYH10 variant overexpression produced a dominant-negative effect on ciliary length.ConclusionThese data presented a novel genetic cause of isolated and syndromic neurodevelopmental disorders related to heterozygous variants in the MYH10 gene with implications for disrupted primary cilia length control and altered Hedgehog signaling in disease pathogenesis.     

Genotype-phenotype correlations in RHOBTB2-associated neurodevelopmental disorders

PurposeMissense variants clustering in the BTB domain region of RHOBTB2 cause a developmental and epileptic encephalopathy with early-onset seizures and severe intellectual disability.MethodsBy international collaboration, we assembled individuals with pathogenic RHOBTB2 variants and a variable spectrum of neurodevelopmental disorders. By western blotting, we investigated the consequences of missense variants in vitro.ResultsIn accordance with previous observations, de novo heterozygous missense variants in the BTB domain region led to a severe developmental and epileptic encephalopathy in 16 individuals. Now, we also identified de novo missense variants in the GTPase domain in 6 individuals with apparently more variable neurodevelopmental phenotypes with or without epilepsy. In contrast to variants in the BTB domain region, variants in the GTPase domain do not impair proteasomal degradation of RHOBTB2 in vitro, indicating different functional consequences. Furthermore, we observed biallelic splice-site and truncating variants in 9 families with variable neurodevelopmental phenotypes, indicating that complete loss of RHOBTB2 is pathogenic as well.ConclusionBy identifying genotype-phenotype correlations regarding location and consequences of de novo missense variants in RHOBTB2 and by identifying biallelic truncating variants, we further delineate and expand the molecular and clinical spectrum of RHOBTB2-related phenotypes, including both autosomal dominant and recessive neurodevelopmental disorders.     

Genome and RNA sequencing in patients with methylmalonic aciduria of unknown cause

PurposeOur laboratory has classified patients with methylmalonic aciduria using somatic cell studies for over four decades. We have accumulated 127 fibroblast lines from patients with persistent elevated methylmalonic acid (MMA) levels in which no genetic cause could be identified. Cultured fibroblasts from 26 of these patients had low [¹⁴C]propionate incorporation into macromolecules, possibly reflecting decreased methylmalonyl-CoA mutase function.MethodsGenome sequencing (GS), copy-number variation (CNV) analysis, and RNA sequencing were performed on genomic DNA and complementary DNA (cDNA) from these 26 patients.ResultsNo patient had two pathogenic variants in any gene associated with cobalamin metabolism. Nine patients had heterozygous variants of unknown significance previously identified by a next-generation sequencing (NGS) panel targeting cobalamin metabolic genes. Three patients had pathogenic changes in genes not associated with cobalamin metabolism (PCCA, EPCAM, and a 17q12 duplication) that explain parts of their phenotypes other than elevated MMA.ConclusionGenome and RNA sequencing did not detect any additional putative causal genetic defects in known cobalamin genes following somatic cell studies and the use of a targeted NGS panel. They did detect pathogenic variants in other genes in three patients that explained some aspects of their clinical presentation.     

Genotype-phenotype correlation in PEX5-deficient peroxisome biogenesis defective cell lines

Proteins destined for the peroxisomal matrix are targeted by virtue of a peroxisomal targeting sequence type 1 (PTS1) or type 2 (PTS2). In humans, targeting of either class of proteins relies on a cytosolic receptor protein encoded by the PEX5 gene. Alternative splicing of PEX5 results in two protein variants, PEX5S and PEX5L. PEX5S is exclusively involved in PTS1 protein import, whereas PEX5L mediates the import of both PTS1 and PTS2 proteins. Genetic complementation testing with over 500 different fibroblast cell lines from patients diagnosed with a peroxisome biogenesis disorder (PBD) identified 11 cell lines with a defect in PEX5. The aim of this study was to characterize these cell lines at a biochemical and genetic level. To this end, the cultured fibroblasts were analyzed for very long chain fatty acid (VLCFA) concentrations, peroxisomal beta-and alpha-oxidation, dihydroxyacetone-phosphate acyltransferase (DHAPAT) activity, peroxisomal thiolase, and catalase immunofluorescence. Mutation analysis of the PEX5 gene revealed 11 different mutations, eight of which are novel. PTS1- and PTS2-protein import capacity was assessed by transfection of the cells with green fluorescent protein (GFP) tagged with either PTS1 or PTS2. Six cell lines showed a defect in both PTS1 and PTS2 protein import, whereas four cell lines only showed a defect in PTS1 protein import. The location of the different mutations within the PEX5 amino acid sequence correlates rather well with the peroxisomal protein import defect observed in the cell lines.     

Histologic and ultrastructural features in early and advanced phases of Zellweger spectrum disorder (infantile Refsum disease)

Zellweger spectrum disorders (ZSD) are rare autosomal recessive inherited metabolic disorders and include severe (Zellweger syndrome) and milder phenotypes [neonatal adrenoleukodystrophy and infantile Refsum disease (IRD)]. ZSD are characterized by impaired peroxisomal functions and lack of peroxisomes detected by electron microscopy (EM). ZSD are caused by mutations in any of the 14 PEX genes. Patients with ZSD commonly demonstrate nonspecific hepatic symptoms within the first year, often without clinical suspicion of ZSD. Thus, recognition of pathologic findings in the liver is critical for the early diagnosis.

We herein demonstrate the histologic and ultrastructural features in liver biopsies in the early and advanced phases from a 16-year-old male with IRD. The initial biopsy at 5 months of age showed a lack of peroxisomes by EM, and this finding played a critical role in the early diagnosis. In contrast, the second biopsy at 14 years of age, after long-term diet therapy, demonstrated significant disease progression with near-cirrhotic liver. In addition to lack of peroxisomes, EM revealed abundant trilamellar inclusions within large angulated lysosomes in many of the hepatocytes and Kupffer cells. Mitochondrial abnormalities were identified only in the second biopsy and were mainly identified in damaged cells; thus they were likely nonspecific secondary changes.

This is the first report demonstrating histological and ultrastructural features of liver biopsies in the early and advanced phases from a child with ZSD. Trilamellar inclusions are considered to be an ultrastructural hallmark of ZSD, but they may not be apparent in the early phases.     

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.