Pathogenic variants in the transmembrane sulfate transporter protein SLC26A2 are associated with different phenotypes of inherited chondrodysplasias. As limited data is published from India, in this study we sought to elucidate the molecular basis of inherited chondrodysplasias in an Indian cohort. Molecular screening of 32 fetuses with antenatally diagnosed lethal skeletal dysplasia was performed by next generation sequencing and Sanger sequencing. The genotype-protein phenotype characterization was done using computational biology techniques like homology modelling, stability and pathogenicity predictions. We identified five rare autosomal recessive SLC26A2 [NM_000112.4] variants, including three homozygous c.796dupA(p.Thr266Asnfs*12), c.1724delA(p.Lys575Serfs*10), and c.1375_1377dup(p.Val459dup) and two heterozygous variants (c.532C > T(p.Arg178*)) and (c.1382C > T(p.Ala461Val)) in compound heterozygous form in a total of four foetuses. Genotype-protein phenotype annotations highlighted that the clinically severe achondrogenesis 1B causative c.796dupA(p.Thr266Asnfs*12) and c.1724delA(p.Lys575Serfs*10)variants impact SLC26A2 protein structure by deletion of the protein core and transmembrane STAS domains, respectively. In clinically moderate atelosteogenesis type 2 phenotype, the c.1382C > T(p.Ala461Val) variant is predicted to distort alpha helix conformation and alter the bonding properties and free energy dynamics of transmembrane domains and the c.532C > T(p.Arg178*) variant results in loss of both core transmembrane and STAS domains of the SLC26A2 protein. The c.1375_1377dup(p.Val459dup) variant identified in clinically milder atelosteogenesis type II-diastrophic dysplasia spectrum lethal phenotype is predicted to decrease the Qualitative Model Energy Analysis (QMean), which affects major geometrical aspects of the SLC26A2 protein structure. We expand the spectrum of SLC26A2 related lethal chondrodysplasia and report three novel variants correlating clinical severity and protein phenotype within the lethal spectrum of this rare dysplasia. We demonstrate the relevance of structural characterization to aid novel variant reclassification to provide better prenatal management and reproductive options to families with lethal antenatal skeletal disorder.

Mutations in the SLC26A2 gene cause a spectrum of currently incurable human chondrodysplasias. However, genotype-phenotype relationships of SLC26A2-deficient chondrodysplasias are still perplexing and thus stunt therapeutic development. To investigate the causative role of SLC26A2 deficiency in chondrodysplasias and confirm its skeleton-specific pathology, we generated and analyzed slc26a2-/- and Col2a1-Cre; slc26a2fl/fl mice. The therapeutic effect of NVP-BGJ398, an FGFR inhibitor, was tested with both explant cultures and timed pregnant females. Two lethal forms of human SLC26A2-related chondrodysplasias, achondrogenesis type IB (ACG1B) and atelosteogenesis type II (AO2), are phenocopied by slc26a2-/- mice. Unexpectedly, slc26a2-/- chondrocytes are defective for collagen secretion, exhibiting intracellular retention and compromised extracellular deposition of ColII and ColIX. As a consequence, the ATF6 arm of the unfolded protein response (UPR) is preferentially triggered to overactivate FGFR3 signaling by inducing excessive FGFR3 in slc26a2-/- chondrocytes. Consistently, suppressing FGFR3 signaling by blocking either FGFR3 or phosphorylation of the downstream effector favors the recovery of slc26a2-/- cartilage cultures from impaired growth and unbalanced cell proliferation and apoptosis. Moreover, administration of an FGFR inhibitor to pregnant females shows therapeutic effects on pathological features in slc26a2-/- newborns. Finally, we confirm the skeleton-specific lethality and pathology of global SLC26A2 deletion through analyzing the Col2a1-Cre; slc26a2fl/fl mouse line. Our study unveils a previously unrecognized pathogenic mechanism underlying ACG1B and AO2, and supports suppression of FGFR3 signaling as a promising therapeutic approach for SLC26A2-related chondrodysplasias. FUND: This work was supported by National Natural Science Foundation of China (81871743, 81730065 and 81772377).