Computational biology insights into genotype-clinical phenotype-protein phenotype relationships between novel SLC26A2 variants identified in inherited skeletal dysplasias

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSION

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.