Neurofibromatosis type 1 (NF1) is one of the most common genetic tumor predisposition syndrome, caused by mutations in the NF1. To date, few genotype-phenotype correlations have been discerned in NF1, due to a highly variable clinical presentation. We aimed to study the molecular spectrum of NF1 and genotype-phenotype correlations in a monocentric study cohort of 85 NF1 patients (20 relatives, 65 sporadic cases). Clinical data were collected at the time of the mutation analysis and reviewed for accuracy in this investigation. An internal phenotypic categorization was applied. The 94% of the patients enrolled showed a severe phenotype with at least one systemic complication and a wide range of associated malignancies. Spine deformities were the most common complications in this cohort. We also reported 66 different NF1 mutations, of which 7 are novel mutations. Correlation analysis identified a slight significant inverse correlation between age at diagnosis and delayed acquisition of psychomotor skills with residual multi-domain cognitive impairment. Odds ratio with 95% confidence interval showed a higher prevalence of learning disabilities in patients carrying frameshift mutations. Overall, our results aim to offer an interesting contribution to studies on the genotype–phenotype of NF1 and in genetic management and counselling.