Methods: The genomes of 24 MTBC isolated from sputum and pus samples were sequenced. The phenotypic drug susceptibility testing (DST) of the isolates was determined for ten anti-TB drugs. Bioinformatic analysis comprising genome assembly and annotation and single-nucleotide polymorphism (SNP) analysis in genes associated with resistance to the ten anti-TB drugs were done on each sequenced genome.

Results: The draft assemblies covered an average of 97% of the expected genome size. Eleven isolates were aligned to the Indo-Oceanic lineage, eight were East-Asian lineage, three were East African-Indian lineage, and one was of Euro-American and Bovis lineages, respectively. Twelve of the 24 MTBC isolates were phenotypically MDR M. tuberculosis: one is polyresistance and another one is monoresistance. Twenty-six SNPs across nine genes associated with resistance toward ten anti-TB drugs were detected where some of the mutations were found in isolates that were previously reported as pan-susceptible using DST. A haplotype consisting of 65 variants was also found among the MTBC isolates with drug-resistance traits.

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.