METHODS: Detailed phenotyping and next-generation sequencing (panel and exome).

RESULTS: Our analysis revealed 224 pathogenic/likely pathogenic variants (54 (24%) of which are novel) in 123 genes with established or tentative links to skeletal dysplasia. In addition, we propose 5 genes as candidate disease genes with suggestive biological links (WNT3A, SUCO, RIN1, DIP2C, and PAN2). Phenotypically, we note that our cohort spans 36 established phenotypic categories by the International Skeletal Dysplasia Nosology, as well as 18 novel skeletal dysplasia phenotypes that could not be classified under these categories, e.g., the novel C3orf17-related skeletal dysplasia. We also describe novel phenotypic aspects of well-known disease genes, e.g., PGAP3-related Toriello-Carey syndrome-like phenotype. We note a strong founder effect for many genes in our cohort, which allowed us to calculate a minimum disease burden for the autosomal recessive forms of skeletal dysplasia in our population (7.16E-04), which is much higher than the global average.

METHODS: A total of 42 affected individuals with homozygous SLC18A2 variant alleles were identified. We evaluated genotype-phenotype correlations and the missense variants in the affected individuals based on the structural modeling of rat VMAT2 encoded by Slc18a2, with cytoplasm- and lumen-facing conformations. A Caenorhabditis elegans model was created for functional studies.

RESULTS: A total of 19 homozygous SLC18A2 variants, including 3 recurrent variants, were identified using exome sequencing. The affected individuals typically showed global developmental delay, hypotonia, dystonia, oculogyric crisis, and autonomic nervous system involvement (temperature dysregulation/sweating, hypersalivation, and gastrointestinal dysmotility). Among the 58 affected individuals described to date, 16 (28%) died before the age of 13 years. Of the 17 patients with p.Pro237His, 9 died, whereas all 14 patients with p.Pro387Leu survived. Although a dopamine agonist mildly improved the disease symptoms in 18 of 21 patients (86%), some affected individuals with p.Ile43Phe and p.Pro387Leu showed milder phenotypes and presented prolonged survival even without treatment. The C. elegans model showed behavioral abnormalities.