ABSTRACT
The human PLAA gene encodes Phospholipase-A2-Activating-Protein (PLAA) involved in trafficking of membrane proteins. Through its PUL domain (PLAP, Ufd3p, and Lub1p), PLAA interacts with p97/VCP modulating synaptic vesicles recycling. Although few families carrying biallelic PLAA variants were reported with progressive neurodegeneration, consequences of monoallelic PLAA variants have not been elucidated. Using exome or genome sequencing we identified PLAA de-novo missense variants, affecting conserved residues within the PUL domain, in children affected with neurodevelopmental disorders (NDDs), including psychomotor regression, intellectual disability (ID) and autism spectrum disorders (ASDs). Computational and in-vitro studies of the identified variants revealed abnormal chain arrangements at C-terminal and reduced PLAA-p97/VCP interaction, respectively. These findings expand both allelic and phenotypic heterogeneity associated to PLAA-related neurological disorders, highlighting perturbed vesicle recycling as a potential disease mechanism in NDDs due to genetic defects of PLAA.
ABSTRACT
Sotos syndrome (SoS) is a neurodevelopmental disorder that results from NSD1 mutations that cause haploinsufficiency of NSD1. Here, we generated an induced pluripotent stem cell (iPSC) line from fibroblasts of a SoS patient carrying the pathogenic variant (c.1633delA). The cell line shows typical iPSC morphology, high expression of pluripotent markers, normal karyotype, and it differentiates into three germ layers in vitro. This line is a valuable resource for studying pathological pathways involved in SoS.
Subject(s)
Craniosynostoses , Induced Pluripotent Stem Cells , Intellectual Disability , Sotos Syndrome , Humans , Sotos Syndrome/genetics , Sotos Syndrome/metabolism , Sotos Syndrome/pathology , Induced Pluripotent Stem Cells/metabolism , Mutation , Exons , Histone-Lysine N-Methyltransferase/geneticsABSTRACT
ZEB2 is a protein-coding gene belonging to a very restricted family of transcription factors. ZEB2 acts mainly as a transcription repressor, is expressed in various tissues and its role is fundamental for the correct development of the nervous system. The best-known clinical picture associated with ZEB2 mutations is Mowat-Wilson syndrome, caused mostly by haploinsufficiency and characterized by possible multi-organ malformations, dysmorphic features, intellectual disability, and epilepsy. In this study we report the generation of IGGi004-A and IGGi005-A, iPSC clones from two patients carrying different heterozygous mutations in ZEB2, which can be used for disease modelling, pathophysiological studies and therapeutics testing.
Subject(s)
Facies , Hirschsprung Disease , Induced Pluripotent Stem Cells , Intellectual Disability , Microcephaly , Humans , Intellectual Disability/complications , Zinc Finger E-box Binding Homeobox 2/genetics , Mutation/genetics , Transcription Factors/genetics , Homeodomain Proteins/geneticsABSTRACT
The Rab family of guanosine triphosphatases (GTPases) includes key regulators of intracellular transport and membrane trafficking targeting specific steps in exocytic, endocytic, and recycling pathways. DENND5B (Rab6-interacting Protein 1B-like protein, R6IP1B) is the longest isoform of DENND5, an evolutionarily conserved DENN domain-containing guanine nucleotide exchange factor (GEF) that is highly expressed in the brain. Through exome sequencing and international matchmaking platforms, we identified five de novo variants in DENND5B in a cohort of five unrelated individuals with neurodevelopmental phenotypes featuring cognitive impairment, dysmorphism, abnormal behavior, variable epilepsy, white matter abnormalities, and cortical gyration defects. We used biochemical assays and confocal microscopy to assess the impact of DENND5B variants on protein accumulation and distribution. Then, exploiting fluorescent lipid cargoes coupled to high-content imaging and analysis in living cells, we investigated whether DENND5B variants affected the dynamics of vesicle-mediated intracellular transport of specific cargoes. We further generated an in silico model to investigate the consequences of DENND5B variants on the DENND5B-RAB39A interaction. Biochemical analysis showed decreased protein levels of DENND5B mutants in various cell types. Functional investigation of DENND5B variants revealed defective intracellular vesicle trafficking, with significant impairment of lipid uptake and distribution. Although none of the variants affected the DENND5B-RAB39A interface, all were predicted to disrupt protein folding. Overall, our findings indicate that DENND5B variants perturb intracellular membrane trafficking pathways and cause a complex neurodevelopmental syndrome with variable epilepsy and white matter involvement.
Subject(s)
Epilepsy , Intellectual Disability , Neurodevelopmental Disorders , Humans , Neurodevelopmental Disorders/genetics , Neurodevelopmental Disorders/metabolism , Brain/metabolism , Epilepsy/genetics , Epilepsy/metabolism , Guanine Nucleotide Exchange Factors/genetics , Guanine Nucleotide Exchange Factors/metabolism , Lipids , Intellectual Disability/genetics , Intellectual Disability/metabolism , rab GTP-Binding Proteins/metabolismABSTRACT
WDR44 prevents ciliogenesis initiation by regulating RAB11-dependent vesicle trafficking. Here, we describe male patients with missense and nonsense variants within the WD40 repeats (WDR) of WDR44, an X-linked gene product, who display ciliopathy-related developmental phenotypes that we can model in zebrafish. The patient phenotypic spectrum includes developmental delay/intellectual disability, hypotonia, distinct craniofacial features and variable presence of brain, renal, cardiac and musculoskeletal abnormalities. We demonstrate that WDR44 variants associated with more severe disease impair ciliogenesis initiation and ciliary signaling. Because WDR44 negatively regulates ciliogenesis, it was surprising that pathogenic missense variants showed reduced abundance, which we link to misfolding of WDR autonomous repeats and degradation by the proteasome. We discover that disease severity correlates with increased RAB11 binding, which we propose drives ciliogenesis initiation dysregulation. Finally, we discover interdomain interactions between the WDR and NH2-terminal region that contains the RAB11 binding domain (RBD) and show patient variants disrupt this association. This study provides new insights into WDR44 WDR structure and characterizes a new syndrome that could result from impaired ciliogenesis.
Subject(s)
Ciliopathies , Genes, X-Linked , WD40 Repeats , Animals , Humans , Male , Brain , Ciliopathies/genetics , Cognition , Zebrafish/geneticsABSTRACT
Cystic Fibrosis Transmembrane conductance Regulator (CFTR) is a chloride channel defective in cystic fibrosis (CF). Several CFTR mutations are causative of CF, among which G542X is a nonsense mutation introducing a premature stop codon which prevents CFTR protein synthesis. We generated a new iPSC line from nasal cells carrying G542X homozygous mutation for CFTR: IGGi002A. This cell line has normal female karyotype, express pluripotency markers and could differentiate into three germ layers in vitro. This iPSC line may be used for disease modeling (cell differentiation and organoid formation) and development of personalized treatments by genome editing or pharmacological screening.
Subject(s)
Cystic Fibrosis , Induced Pluripotent Stem Cells , Humans , Female , Cystic Fibrosis/genetics , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Induced Pluripotent Stem Cells/metabolism , Mutation/genetics , Codon, NonsenseABSTRACT
Neurofibromatosis type 1 (NF1) is a neurocutaneous disorder caused by mutations in NF1 gene, coding for neurofibromin 1. NF1 can be associated with Moyamoya disease (MMD), and this association, typical of paediatric patients, is referred to as Moyamoya syndrome (MMS). MMD is a cerebral arteriopathy characterized by the occlusion of intracranial arteries and collateral vessel formation, which increase the risk of ischemic and hemorrhagic events. RNF213 gene mutations have been associated with MMD, so we investigated whether rare variants of RNF213 could act as genetic modifiers of MMS phenotype in a pediatric cohort of 20 MMS children, 25 children affected by isolated MMD and 47 affected only by isolated NF1. By next-generation re-sequencing (NGS) of patients' DNA and gene burden tests, we found that RNF213 seems to play a role only for MMD occurrence, while it does not appear to be involved in the increased risk of Moyamoya for MMS patients. We postulated that the loss of neurofibromin 1 can be enough for the excessive proliferation of vascular smooth muscle cells, causing Moyamoya arteriopathy associated with NF1. Further studies will be crucial to support these findings and to elucidate the possible role of other genes, enhancing our knowledge about pathogenesis and treatment of MMS.
ABSTRACT
Sotos syndrome (SoS) is a neurodevelopmental disorder caused by haploinsufficiency of the NSD1 gene located on chromosome 5 region q35.3. In order to understand the pathogenesis of Sotos syndrome and in view of future therapeutic approaches for its efficient treatment, we generated two human induced pluripotent stem cells (iPSCs) lines from one SoS patient carrying a 5q35 microdeletion. The established iPSCs expressed pluripotency markers, showing the capacity to differentiate into the three germ layers.
Subject(s)
Induced Pluripotent Stem Cells , Sotos Syndrome , Humans , Sotos Syndrome/genetics , Sotos Syndrome/pathology , Induced Pluripotent Stem Cells/pathology , Histone Methyltransferases/genetics , Histone-Lysine N-Methyltransferase/genetics , HaploinsufficiencyABSTRACT
Cerebellar hypoplasia and dysplasia encompass a group of clinically and genetically heterogeneous disorders frequently associated with neurodevelopmental impairment. The Neuron Navigator 2 (NAV2) gene (MIM: 607,026) encodes a member of the Neuron Navigator protein family, widely expressed within the central nervous system (CNS), and particularly abundant in the developing cerebellum. Evidence across different species supports a pivotal function of NAV2 in cytoskeletal dynamics and neurite outgrowth. Specifically, deficiency of Nav2 in mice leads to cerebellar hypoplasia with abnormal foliation due to impaired axonal outgrowth. However, little is known about the involvement of the NAV2 gene in human disease phenotypes. In this study, we identified a female affected with neurodevelopmental impairment and a complex brain and cardiac malformations in which clinical exome sequencing led to the identification of NAV2 biallelic truncating variants. Through protein expression analysis and cell migration assay in patient-derived fibroblasts, we provide evidence linking NAV2 deficiency to cellular migration deficits. In model organisms, the overall CNS histopathology of the Nav2 hypomorphic mouse revealed developmental anomalies including cerebellar hypoplasia and dysplasia, corpus callosum hypo-dysgenesis, and agenesis of the olfactory bulbs. Lastly, we show that the NAV2 ortholog in Drosophila, sickie (sick) is widely expressed in the fly brain, and sick mutants are mostly lethal with surviving escapers showing neurobehavioral phenotypes. In summary, our results unveil a novel human neurodevelopmental disorder due to genetic loss of NAV2, highlighting a critical conserved role of the NAV2 gene in brain and cerebellar development across species.
Subject(s)
Brain , Nervous System Malformations , Animals , Female , Humans , Mice , Cerebellum/abnormalities , NeuronsABSTRACT
Segmental overgrowth syndromes include a group of clinical entities, all characterized by the abundant proliferation of tissues or organs in association with vascular abnormalities. These syndromes show a wide spectrum of severity ranging from limited involvement of only small areas of the body to complex cases with impressive distortions of multiple tissues and organs. It is now clear that somatic mutations in genes of the phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT)/mammalian target of rapamycin (mTOR) pathway (in brief "mTOR pathway") are responsible for such entities. Not all the cells of the body carry the same causative mutation, which is mosaic, appearing from two (or more) distinct cell lineages after fertilization. In this article, we reconsider the clinical spectrum and surveillance programs of patients with segmental overgrowth syndromes, based on the features of six patients with diverse clinical forms of overgrowth and pathogenic variants in genes of the mTOR pathway.
Subject(s)
Phosphatidylinositol 3-Kinases , TOR Serine-Threonine Kinases , Humans , Phosphatidylinositol 3-Kinases/genetics , Phosphatidylinositol 3-Kinases/metabolism , Class I Phosphatidylinositol 3-Kinases/genetics , Class I Phosphatidylinositol 3-Kinases/metabolism , TOR Serine-Threonine Kinases/genetics , TOR Serine-Threonine Kinases/metabolism , Genetic Testing , Mutation , SyndromeABSTRACT
Human-induced pluripotent stem cells (hiPSCs) represent one of the main and powerful tools for the in vitro modeling of neurological diseases. Standard hiPSC-based protocols make use of animal-derived feeder systems to better support the neuronal differentiation process. Despite their efficiency, such protocols may not be appropriate to dissect neuronal specific properties or to avoid interspecies contaminations, hindering their future translation into clinical and drug discovery approaches. In this work, we focused on the optimization of a reproducible protocol in feeder-free conditions able to generate functional glutamatergic neurons. This protocol is based on a generation of neuroprecursor cells differentiated into human neurons with the administration in the culture medium of specific neurotrophins in a Geltrex-coated substrate. We confirmed the efficiency of this protocol through molecular analysis (upregulation of neuronal markers and neurotransmitter receptors assessed by gene expression profiling and expression of the neuronal markers at the protein level), morphological analysis, and immunfluorescence detection of pre-synaptic and post-synaptic markers at synaptic boutons. The hiPSC-derived neurons acquired Ca2+-dependent glutamate release properties as a hallmark of neuronal maturation. In conclusion, our study describes a new methodological approach to achieve feeder-free neuronal differentiation from hiPSC and adds a new tool for functional characterization of hiPSC-derived neurons.
Subject(s)
Glutamic Acid , Induced Pluripotent Stem Cells , Animals , Cell Differentiation/genetics , Glutamic Acid/metabolism , Humans , Nerve Growth Factors/metabolism , Neurons/metabolism , Receptors, Neurotransmitter/metabolismABSTRACT
BACKGROUND: Pathogenic variants in PEX-genes can affect peroxisome assembly and function and cause Zellweger spectrum disorders (ZSDs), characterized by variable phenotypes in terms of disease severity, age of onset and clinical presentations. So far, defects in at least 15 PEX-genes have been implicated in Mendelian diseases, but in some of the ultra-rare ZSD subtypes genotype-phenotype correlations and disease mechanisms remain elusive. METHODS: We report five families carrying biallelic variants in PEX13. The identified variants were initially evaluated by using a combination of computational approaches. Immunofluorescence and complementation studies on patient-derived fibroblasts were performed in two patients to investigate the cellular impact of the identified mutations. RESULTS: Three out of five families carried a recurrent p.Arg294Trp non-synonymous variant. Individuals affected with PEX13-related ZSD presented heterogeneous clinical features, including hypotonia, developmental regression, hearing/vision impairment, progressive spasticity and brain leukodystrophy. Computational predictions highlighted the involvement of the Arg294 residue in PEX13 homodimerization, and the analysis of blind docking predicted that the p.Arg294Trp variant alters the formation of dimers, impairing the stability of the PEX13/PEX14 translocation module. Studies on muscle tissues and patient-derived fibroblasts revealed biochemical alterations of mitochondrial function and identified mislocalized mitochondria and a reduced number of peroxisomes with abnormal PEX13 concentration. CONCLUSIONS: This study expands the phenotypic and mutational spectrum of PEX13-related ZSDs and also highlight a variety of disease mechanisms contributing to PEX13-related clinical phenotypes, including the emerging contribution of secondary mitochondrial dysfunction to the pathophysiology of ZSDs.
Subject(s)
Zellweger Syndrome , Genetic Association Studies , Humans , Membrane Proteins/genetics , Mutation/genetics , Peroxisomes/genetics , Peroxisomes/pathology , Zellweger Syndrome/genetics , Zellweger Syndrome/pathologyABSTRACT
BACKGROUND: Cystic Fibrosis (CF) is a genetic disorder affecting around 1 in every 3000 newborns. In the most common mutation, F508del, the defective anion channel, CFTR, is prevented from reaching the plasma membrane (PM) by the quality check control of the cell. Little is known about how CFTR pharmacological rescue impacts the cell proteome. METHODS: We used high-resolution mass spectrometry, differential ultracentrifugation, machine learning and bioinformatics to investigate both changes in the expression and localization of the human bronchial epithelium CF model (F508del-CFTR CFBE41o-) proteome following treatment with VX-809 (Lumacaftor), a drug able to improve the trafficking of CFTR. RESULTS: The data suggested no stark changes in protein expression, yet subtle localization changes of proteins of the mitochondria and peroxisomes were detected. We then used high-content confocal microscopy to further investigate the morphological and compositional changes of peroxisomes and mitochondria under these conditions, as well as in patient-derived primary cells. We profiled several thousand proteins and we determined the subcellular localization data for around 5000 of them using the LOPIT-DC spatial proteomics protocol. CONCLUSIONS: We observed that treatment with VX-809 induces extensive structural and functional remodelling of mitochondria and peroxisomes that resemble the phenotype of healthy cells. Our data suggest additional rescue mechanisms of VX-809 beyond the correction of aberrant folding of F508del-CFTR and subsequent trafficking to the PM.
Subject(s)
Cystic Fibrosis , Aminopyridines , Benzodioxoles , Cystic Fibrosis/metabolism , Cystic Fibrosis Transmembrane Conductance Regulator/metabolism , Epithelium/metabolism , Humans , Infant, Newborn , Mitochondria/metabolism , Proteome/metabolismABSTRACT
Background and Objectives: Clinical manifestations in STXBP1 developmental and epileptic encephalopathy (DEE) vary in severity and outcome, and the genotypic spectrum is diverse. We aim to trace the neurodevelopmental trajectories in individuals with STXBP1-DEE and dissect the relationship between neurodevelopment and epilepsy. Methods: Retrospective standardized clinical data were collected through international collaboration. A composite neurodevelopmental score system compared the developmental trajectories in STXBP1-DEE. Results: Forty-eight patients with de novo STXBP1 variants and a history of epilepsy were included (age range at the time of the study: 10 months to 35 years, mean 8.5 years). At the time of inclusion, 65% of individuals (31/48) had active epilepsy, whereas 35% (17/48) were seizure free, and 76% of those (13/17) achieved remission within the first year of life. Twenty-two individuals (46%) showed signs of developmental impairment and/or neurologic abnormalities before epilepsy onset. Age at seizure onset correlated with severity of developmental outcome and the developmental milestones achieved, with a later seizure onset associated with better developmental outcome. In contrast, age at seizure remission and epilepsy duration did not affect neurodevelopmental outcomes. Overall, we did not observe a clear genotype-phenotype correlation, but monozygotic twins with de novo STXBP1 variant showed similar phenotype and parallel disease course. Discussion: The disease course in STXBP1-DEE presents with 2 main trajectories, with either early seizure remission or drug-resistant epilepsy, and a range of neurodevelopmental outcomes from mild to profound intellectual disability. Age at seizure onset is the only epilepsy-related feature associated with neurodevelopment outcome. These findings can inform future dedicated natural history studies and trial design.
ABSTRACT
Background: WOREE syndrome is a rare neurodevelopmental disorder featuring drug-resistant epilepsy and global developmental delay. The disease, caused by biallelic pathogenic variants in the WWOX gene, usually leads to severe disability or death within the first years of life. Clinicians have become more confident with the phenotypic picture of WOREE syndrome, allowing earlier clinical diagnosis. We report a boy with a peculiar clinic-radiological pattern supporting the diagnosis of WOREE syndrome. Methods: DNA was extracted from blood samples of the proband and his parents and subjected to Exome Sequencing (ES). Agarose gel electrophoresis, real-time quantitative PCR (Q-PCR), and array-CGH 180K were also performed. Results: ES detected a pathogenic stop variant (c.790C > T, p.Arg264*) in one allele of WWOX in the proband and his unaffected mother. A 180K array-CGH analysis revealed a 84,828-bp (g.chr16:78,360,803-78,445,630) deletion encompassing exon 6. The Q-PCR product showed that the proband and his father harbored the same deleted fragment, fusing exons 5 and 7 of WWOX. Conclusions: Genetic testing remains crucial in establishing the definitive diagnosis of WOREE syndrome and allows prenatal interventions/parental counseling. However, our findings suggest that targeted Next Generation Sequencing-based testing may occasionally show technical pitfalls, prompting further genetic investigation in selected cases with high clinical suspicion.
ABSTRACT
Wieacker-Wolff syndrome (WWS) is an X-linked Arthrogryposis Multiplex Congenita (AMC) disorder associated with broad neurodevelopmental impairment. The genetic basis of WWS lies in hemizygous pathogenic variants in ZC4H2, encoding a C4H2 type zinc-finger nuclear factor abundantly expressed in the developing human brain. The main clinical features described in WWS families carrying ZC4H2 pathogenic variants encompass having a short stature, microcephaly, birth respiratory distress, arthrogryposis, hypotonia, distal muscle weakness, and broad neurodevelopmental delay. We hereby report a Sicilian family with a boy clinically diagnosed with WWS and genetically investigated with exome sequencing (ES), leading to the identification of a c.593G>A (p. R198Q) hemizygous pathogenic variant in the ZC4H2 gene. During the first year of life, the onset of central hypoadrenalism led to recurrent hypoglycemic events, which likely contributed to seizure susceptibility. Also, muscle biopsy studies confirmed a pathology of the muscle tissue and revealed peculiar abnormalities of the neuromuscular junction. In conclusion, we expand the phenotypic spectrum of the WWS-related neurodevelopmental disorders and discuss the role of ZC4H2 in the context of the potential neuroendocrinological and neuromuscular features associated with this condition.
ABSTRACT
BACKGROUND: Currarino syndrome (CS) is a rare genetic condition characterized by the association of three major clinical signs: anorectal malformation (ARM), sacro-coccygeal bone defects, and presacral mass. Different kinds of ARM can be present such as anteriorly placed anus, imperforate anus, anorectal stenosis, rectal duplication, and fistulae. The presacral mass can be a benign teratoma, a dermoid or neurenteric cyst, anterior meningocele or hamartoma. Females are more frequently affected and usually present with associated gynecologic and urinary tract problems. CS is considered an autosomal dominant trait, with reduced penetrance and variable expressivity. CS is associated with mutations in the MNX1 gene (motor neuron and pancreas homeobox-1, previously known as HLXB9) mapped to chromosome 7q36. Heterozygous loss-of-function mutations in the coding sequence of MNX1 gene have been reported in nearly all familial CS cases and in approximately 30% of CS sporadic patients. CASE: Here, we present the case of a woman with features of CS carrying a mosaic mutation in the coding region of MNX1 gene. This is the only reported case of a CS diagnosis in which the mutation is present in less than 50% of cells. CONCLUSION: The lower detection rate of MNX1 mutations in sporadic cases could similarly be explained by somatic mosaicism, mutations occurring outside the coding regions, or genetic heterogeneity.
Subject(s)
Anal Canal , Genes, Homeobox , Adult , Anal Canal/abnormalities , Digestive System Abnormalities , Female , Homeodomain Proteins/genetics , Humans , Mutation/genetics , Rectum/abnormalities , Sacrum/abnormalities , Syringomyelia , Transcription Factors/geneticsABSTRACT
OBJECTIVE: To describe the clinical and genetic findings in a cohort of individuals with bathing epilepsy, a rare form of reflex epilepsy. METHODS: We investigated by Sanger and targeted resequencing the SYN1 gene in 12 individuals from 10 different families presenting with seizures triggered primarily by bathing or showering. An additional 12 individuals with hot-water epilepsy were also screened. RESULTS: In all families with bathing epilepsy, we identified 8 distinct pathogenic or likely pathogenic variants and 2 variants of unknown significance in SYN1, 9 of which are novel. Conversely, none of the individuals with hot-water epilepsy displayed SYN1 variants. In mutated individuals, seizures were typically triggered by showering or bathing regardless of the water temperature. Additional triggers included fingernail clipping, haircutting, or watching someone take a shower. Unprovoked seizures and a variable degree of developmental delay were also common. CONCLUSION: Bathing epilepsy is genetically distinct reflex epilepsy caused mainly by SYN1 mutations.
Subject(s)
Baths , Epilepsy, Reflex/genetics , Epilepsy, Reflex/physiopathology , Hygiene , Synapsins/genetics , Adolescent , Child , Child, Preschool , Female , Hot Temperature , Humans , Male , Middle Aged , Pedigree , WaterABSTRACT
Neurofibromatosis type 1 (NF1) is a proteiform genetic condition caused by pathogenic variants in NF1 and characterized by a heterogeneous phenotypic presentation. Relevant genotype-phenotype correlations have recently emerged, but only few pertinent studies are available. We retrospectively reviewed clinical, instrumental, and genetic data from a cohort of 583 individuals meeting at least 1 diagnostic National Institutes of Health (NIH) criterion for NF1. Of these, 365 subjects fulfilled ≥2 NIH criteria, including 235 pediatric patients. Genetic testing was performed through cDNA-based sequencing, Next Generation Sequencing (NGS), and Multiplex Ligation-dependent Probe Amplification (MLPA). Uni- and multivariate statistical analysis was used to investigate genotype-phenotype correlations. Among patients fulfilling ≥ 2 NIH criteria, causative single nucleotide variants (SNVs) and copy number variations (CNVs) were detected in 267/365 (73.2%) and 20/365 (5.5%) cases. Missense variants negatively correlated with neurofibromas (p = 0.005). Skeletal abnormalities were associated with whole gene deletions (p = 0.05) and frameshift variants (p = 0.006). The c.3721C>T; p.(R1241*) variant positively correlated with structural brain alterations (p = 0.031), whereas Lisch nodules (p = 0.05) and endocrinological disorders (p = 0.043) were associated with the c.6855C>A; p.(Y2285*) variant. We identified novel NF1 genotype-phenotype correlations and provided an overview of known associations, supporting their potential relevance in the implementation of patient management.