Multi-gene panel sequencing in highly consanguineous families and patients with congenital forms of skeletal dysplasias.

Skeletal dysplasias (SKDs) are a heterogeneous group of more than 750 genetic disorders characterized by abnormal development, growth, and maintenance of bones or cartilage in the human skeleton. SKDs are often caused by variants in early patterning genes and in many cases part of multiple malformation syndromes and occur in combination with non-skeletal phenotypes. The aim of this study was to investigate the underlying genetic cause of congenital SKDs in highly consanguineous Pakistani families, as well as in sporadic and familial SKD cases from India using multigene panel sequencing analysis. Therefore, we performed panel sequencing of 386 bone-related genes in 7 highly consanguineous families from Pakistan and 27 cases from India affected with SKDs. In the highly consanguineous families, we were able to identify the underlying genetic cause in five out of seven families, resulting in a diagnostic yield of 71%. Whereas, in the sporadic and familial SKD cases, we identified 12 causative variants, corresponding to a diagnostic yield of 44%. The genetic heterogeneity in our cohorts was very high and we were able to detect various types of variants, including missense, nonsense, and frameshift variants, across multiple genes known to cause different types of SKDs. In conclusion, panel sequencing proved to be a highly effective way to decipher the genetic basis of SKDs in highly consanguineous families as well as sporadic and or familial cases from South Asia. Furthermore, our findings expand the allelic spectrum of skeletal dysplasias.

Two sisters with RSPRY1-related spondyloepimetaphyseal dysplasia.

Biallelic variants in RSPRY1 have been found to result in spondyloepimetaphyseal dysplasia. Two siblings presenting with short stature, facial dysmorphism, progressive vertebral defects, small epiphysis, cupping and fraying of metaphyses, brachydactyly, and short metatarsals harbored a homozygous missense variant c.1652G>A;p.(Cys551Tyr) in the RSPRY1 gene. The phenotype in our patients resembles spondyloepimetaphyseal dysplasia, Faden-Alkuraya type. Thus, our study provides further evidence to support the association of RSPRY1 variants with spondyloepimetaphyseal dysplasia. We observed joint dislocation as a novel clinical feature of this condition.

Indian patients with CHST3-related chondrodysplasia with congenital joint dislocations.

CHST3-related chondrodysplasia with congenital joint dislocations (CDCJD, #MIM 143095), is a rare genetic skeletal disorder caused by biallelic loss of function variants in CHST3. CHST3 is critical for the sulfation of chondroitin sulfate. This study delineates the clinical presentation of nine individuals featuring the key symptoms of CDCJD; congenital joint (knee and elbow) dislocations, short trunk short stature progressive vertebral anomalies, and metacarpal shortening. Additional manifestations include irregular distal femoral epiphysis, supernumerary carpal ossification centers, bifid humerus, club foot, and cardiac abnormalities. Sanger sequencing was carried out to investigate molecular etiology in eight patients and exome sequencing in one. Genetic testing revealed five homozygous variants in CHST3 (four were novel and one was previously reported). All these variants are located on sulfotransferase domain of CHST3 protein and were classified as pathogenic/ likely pathogenic. We thus report on nine individuals with CHST3-related chondrodysplasia with congenital joint dislocations from India and suggest monitoring the health of cardiac valves in this condition.

PRKACA-related, atrial defects-polydactyly-multiple congenital malformation syndrome in an Indian patient.

PRKACA-related, atrial defects-polydactyly-multiple congenital malformation syndrome is a recently described skeletal ciliopathy, which is caused by disease-causing variants in PRKACA. The primary phenotypic description includes atrial septal defects, and limb anomalies including polydactyly and short limbs. To date, only four molecularly proven patients have been reported in the literature with a recurrent variant, c.409G>A p.Gly137Arg in PRKACA. In this study, we report the fifth affected individual with the same variant and review the clinical features and radiographic findings of this rare syndrome.

Clinical, radiological and molecular studies in 24 individuals with Dyggve-Melchior-Clausen dysplasia and Smith-McCort dysplasia from India.

BACKGROUND: Dyggve-Melchior-Clausen dysplasia (DMC) and Smith-McCort dysplasia (SMC types 1 and 2) are rare spondyloepimetaphyseal dysplasias with identical radiological findings. The presence of intellectual disability in DMC and normal intellect in SMC differentiates the two. DMC and SMC1 are allelic and caused by homozygous or compound heterozygous variants in DYM. SMC2 is caused by variations in RAB33B. Both DYM and RAB33B are important in intravesicular transport and function in the Golgi apparatus. METHODS: Detailed clinical phenotyping and skeletal radiography followed by molecular testing were performed in all affected individuals. Next-generation sequencing and Sanger sequencing were used to confirm DYM and RAB33B variants. Sanger sequencing of familial variants was done in all parents. RESULTS: 24 affected individuals from seven centres are described. 18 had DMC and 6 had SMC2. Parental consanguinity was present in 15 of 19 (79%). Height <3 SD and gait abnormalities were seen in 20 and 14 individuals, respectively. The characteristic radiological findings of lacy iliac crests and double-humped vertebral bodies were seen in 96% and 88% of the affected. Radiological findings became attenuated with age. 23 individuals harboured biallelic variants in either DYM or RAB33B. Fourteen different variants were identified, out of which 10 were novel. The most frequently occurring variants in this group were c.719 C>A (3), c.1488_1489del (2), c.1484dup (2) and c.1563+2T>C (2) in DYM and c.400C>T (2) and c.186del (2) in RAB33B. The majority of these have not been reported previously. CONCLUSION: This large cohort from India contributes to the increasing knowledge of clinical and molecular findings in these rare 'Golgipathies'.

Homozygous loss-of-function variants in FILIP1 cause autosomal recessive arthrogryposis multiplex congenita with microcephaly.

Arthrogryposis multiplex congenita forms a broad group of clinically and etiologically heterogeneous disorders characterized by congenital joint contractures that involve at least two different parts of the body. Neurological and muscular disorders are commonly underlying arthrogryposis. Here, we report five affected individuals from three independent families sharing an overlapping phenotype with congenital contractures affecting shoulder, elbow, hand, hip, knee and foot as well as scoliosis, reduced palmar and plantar skin folds, microcephaly and facial dysmorphism. Using exome sequencing, we identified homozygous truncating variants in FILIP1 in all patients. FILIP1 is a regulator of filamin homeostasis required for the initiation of cortical cell migration in the developing neocortex and essential for the differentiation process of cross-striated muscle cells during myogenesis. In summary, our data indicate that bi-allelic truncating variants in FILIP1 are causative of a novel autosomal recessive disorder and expand the spectrum of genetic factors causative of arthrogryposis multiplex congenita.

Genome sequencing identifies a large non-coding region deletion of SNX10 causing autosomal recessive osteopetrosis.

Autosomal recessive osteopetrosis (ARO) is a rare genetic disorder caused by impaired osteoclast activity. In this study, we describe a 4-year-old boy with increased bone density due to osteopetrosis, autosomal recessive 8. Using genome sequencing, we identified a large deletion in the 5'-untranslated region (UTR) of SNX10 (sorting nexin 10), where the regulatory region of this gene is located. This large deletion resulted in the absence of the SNX10 transcript and led to abnormal osteoclast activity. SNX10 is one of the nine genes known to cause ARO, shown to interact with V-ATPase (vacuolar type H( + )-ATPase), as it plays an important role in bone resorption. Our study highlights the importance of regulatory regions in the 5'-UTR of SNX10 for its expression while also demonstrating the importance of genome sequencing for detecting large deletion of the regulatory region of SNX10.

Novel subtype of mucopolysaccharidosis caused by arylsulfatase K (ARSK) deficiency.

BACKGROUND: Mucopolysaccharidoses (MPS) are monogenic metabolic disorders that significantly affect the skeleton. Eleven enzyme defects in the lysosomal degradation of glycosaminoglycans (GAGs) have been assigned to the known MPS subtypes (I-IX). Arylsulfatase K (ARSK) is a recently characterised lysosomal hydrolase involved in GAG degradation that removes the 2-O-sulfate group from 2-sulfoglucuronate. Knockout of Arsk in mice was consistent with mild storage pathology, but no human phenotype has yet been described. METHODS: In this study, we report four affected individuals of two unrelated consanguineous families with homozygous variants c.250C>T, p.(Arg84Cys) and c.560T>A, p.(Leu187Ter) in ARSK, respectively. Functional consequences of the two ARSK variants were assessed by mutation-specific ARSK constructs derived by site-directed mutagenesis, which were ectopically expressed in HT1080 cells. Urinary GAG excretion was analysed by dimethylene blue and electrophoresis, as well as liquid chromatography/mass spectrometry (LC-MS)/MS analysis. RESULTS: The phenotypes of the affected individuals include MPS features, such as short stature, coarse facial features and dysostosis multiplex. Reverse phenotyping in two of the four individuals revealed additional cardiac and ophthalmological abnormalities. Mild elevation of dermatan sulfate was detected in the two subjects investigated by LC-MS/MS. Human HT1080 cells expressing the ARSK-Leu187Ter construct exhibited absent protein levels by western blot, and cells with the ARSK-Arg84Cys construct showed markedly reduced enzyme activity in an ARSK-specific enzymatic assay against 2-O-sulfoglucuronate-containing disaccharides as analysed by C18-reversed-phase chromatography followed by MS. CONCLUSION: Our work provides a detailed clinical and molecular characterisation of a novel subtype of mucopolysaccharidosis, which we suggest to designate subtype X.

Ectodysplasin pathogenic variants affecting the furin-cleavage site and unusual clinical features define X-linked hypohidrotic ectodermal dysplasia in India.

Hypohidrotic ectodermal dysplasia (HED) is a rare genetic disorder caused by mutational inactivation of a developmental pathway responsible for generation of tissues of ectodermal origin. The X-linked form accounts for the majority of HED cases and is caused by Ectodysplasin (EDA) pathogenic variants. We performed a combined analysis of 29 X-linked hypohidrotic ectodermal dysplasia (XLHED) families (including 12 from our previous studies). In addition to the classical triad of symptoms including loss (or reduction) of ectodermal structures, such as hair, teeth, and sweat glands, we detected additional HED-related clinical features including facial dysmorphism and hyperpigmentation in several patients. Interestingly, global developmental delay was identified as an unusual clinical symptom in many patients. More importantly, we identified 22 causal pathogenic variants that included 15 missense, four small in-dels, and one nonsense, splice site, and large deletion each. Interestingly, we detected 12 unique (India-specific) pathogenic variants. Of the 29 XLHED families analyzed, 11 (38%) harbored pathogenic variant localized to the furin cleavage site. A comparison with HGMD revealed significant differences in the frequency of missense pathogenic variants; involvement of specific exons and/or protein domains and transition/transversion ratios. A significantly higher proportion of missense pathogenic variants (33%) localized to the EDA furin cleavage when compared to HGMD (7%), of which p.R155C, p.R156C, and p.R156H were detected in three families each. Therefore, the first comprehensive analysis of XLHED from India has revealed several unique features including unusual clinical symptoms and high frequency of furin cleavage site pathogenic variants.

Bosley-Salih-Alorainy syndrome in patients from India.

Bi-allelic HOXA1 pathogenic variants clinically manifest as two distinct syndromes, Bosley-Salih-Alorainy syndrome (BSAS) and Athabascan brainstem dysgenesis syndrome, mainly reported in two different populations from Saudi Arabia and southwest North America, respectively. Here we report two siblings of Indian origin with BSAS phenotype caused by a novel homozygous exon 2 HOXA1 pathogenic variants.

The homozygous variant c.797G>A/p.(Cys266Tyr) in PISD is associated with a Spondyloepimetaphyseal dysplasia with large epiphyses and disturbed mitochondrial function.

Spondyloepimetaphyseal dysplasias (SEMD) are a group of genetically heterogeneous skeletal disorders characterized by abnormal vertebral bodies and epimetaphyseal abnormalities. We investigated two families with a new SEMD type with one proband each. They showed mild facial dysmorphism, flat vertebral bodies (platyspondyly), large epiphyses, metaphyseal dysplasia, and hallux valgus as common clinical features. By trio-exome sequencing, the homozygous missense variant c.797G>A/p.(Cys266Tyr) in PISD was found in both affected individuals. Based on exome data analyses for homozygous regions, the two patients shared a single homozygous block on chromosome 22 including PISD, indicating their remote consanguinity. PISD encodes phosphatidylserine (PS) decarboxylase that is localized in the inner mitochondrial membrane and catalyzes the decarboxylation of PS to phosphatidylethanolamine (PE) in mammalian cells. PE occurs at high abundance in mitochondrial membranes. Patient-derived fibroblasts showed fragmented mitochondrial morphology. Treatment of patient cells with MG-132 or staurosporine to induce activation of the intrinsic apoptosis pathway revealed significantly decreased cell viability with increased caspase-3 and caspase-7 activation. Remarkably, ethanolamine (Etn) supplementation largely restored cell viability and enhanced apoptosis in MG-132-stressed patient cells. Our data demonstrate that the biallelic hypomorphic PISD variant p.(Cys266Tyr) is associated with a novel SEMD form, which may be treatable with Etn administration.

Homozygous variant, p.(Arg643Trp) in VAC14 causes striatonigral degeneration: report of a novel variant and review of VAC14-related disorders.

VAC14-related disorders include two distinct phenotypes, striatonigral degeneration [MIM# 617054] and Yunis-Varon syndrome. Striatonigral degeneration is a recently described childhood onset dystonia caused by pathogenic variants in VAC14. It is characterized by a period of apparent normalcy followed by abrupt onset neuroregression, dystonia, involuntary movements and degenerative brain lesions involving caudate nucleus, putamen and substantia nigra. Yunis-Varon syndrome is a well described severe condition characterised by skeletal findings and dysmorphism along with neuronal degeneration. Pathogenic variants in FIG4 have been previously reported to cause Yunis-Varon syndrome. Recently, loss of function variants in VAC14 were also reported in an individual affected with Yunis-Varon syndrome. Total seven individuals from four families are reported to have VAC14-related disorders till date. Here, we report another individual with clinical and radiological features suggestive of striatonigral degeneration with homozygous missense variant in VAC14. The patient fibroblasts showed extensive vacuolization, characteristic of VAC14-related disorders. We also review the phenotype and genotype associated with these disorders.

Additional three patients with Smith-McCort dysplasia due to novel RAB33B mutations.

Smith-McCort dysplasia (SMC OMIM 615222) and Dyggve-Melchior-Clausen dysplasia (DMC OMIM 223800) are allelic skeletal dysplasias caused by homozygous or compound heterozygous mutations in DYM (OMIM 607461). Both disorders share the same skeletal phenotypes characterized by spondylo-epi-metaphyseal dysplasia with distinctive lacy ilia. The difference rests on the presence or absence of intellectual disability, that is, intellectual disability in DMC and normal cognition in SMC. However, genetic heterogeneity was suspected in SMC. Recently, RAB33B (OMIM 605950) has been identified as the second gene for SMC. Nevertheless, only two affected families have been reported so far. Here we present three SMC patients with four novel pathogenic variants in RAB33B, including homozygosity for c.211C>T (p.R71*), homozygosity for c.365T>C (p.F122S), and compound heterozygosity for c.48delCGGGGCAG (p.G17Vfs*58) and c.490C>T (p.Q164*). We also summarize the clinical, radiological, and mutation profile of RAB33B after literature mining. This report ascertains the pathogenic relationship between RAB33B and SMC. © 2017 Wiley Periodicals, Inc.

Clinical and mutation profile of multicentric osteolysis nodulosis and arthropathy.

​Multicentric osteolysis nodulosis and arthropathy (MONA) is an infrequently described autosomal recessive skeletal dysplasia characterized by progressive osteolysis and arthropathy. Inactivating mutations in MMP2, encoding matrix metalloproteinase-2, are known to cause this disorder. Fifteen families with mutations in MMP2 have been reported in literature. In this study we screened thirteen individuals from eleven families for MMP2 mutations and identified eight mutations (five novel and three known variants). We characterize the clinical, radiographic and molecular findings in all individuals with molecularly proven MONA from the present cohort and previous reports, and provide a comprehensive review of the MMP2 related disorders.

Spectrum of SMPD1 mutations in Asian-Indian patients with acid sphingomyelinase (ASM)-deficient Niemann-Pick disease.

Acid sphingomyelinase (ASM)-deficient Niemann-Pick disease is an autosomal recessive lysosomal storage disorder caused by biallelic mutations in the SMPD1 gene. To date, around 185 mutations have been reported in patients with ASM-deficient NPD world-wide, but the mutation spectrum of this disease in India has not yet been reported. The aim of this study was to ascertain the mutation profile in Indian patients with ASM-deficient NPD. We sequenced SMPD1 in 60 unrelated families affected with ASM-deficient NPD. A total of 45 distinct pathogenic sequence variants were found, of which 14 were known and 31 were novel. The variants included 30 missense, 4 nonsense, and 9 frameshift (7 single base deletions and 2 single base insertions) mutations, 1 indel, and 1 intronic duplication. The pathogenicity of the novel mutations was inferred with the help of the mutation prediction software MutationTaster, SIFT, Polyphen-2, PROVEAN, and HANSA. The effects of the identified sequence variants on the protein structure were studied using the structure modeled with the help of the SWISS-MODEL workspace program. The p. (Arg542*) (c.1624C>T) mutation was the most commonly identified mutation, found in 22% (26 out of 120) of the alleles tested, but haplotype analysis for this mutation did not identify a founder effect for the Indian population. To the best of our knowledge, this is the largest study on mutation analysis of patients with ASM-deficient Niemann-Pick disease reported in literature and also the first study on the SMPD1 gene mutation spectrum in India. © 2016 Wiley Periodicals, Inc.

A syndrome of facial dysmorphism, cubital pterygium, short distal phalanges, swan neck deformity of fingers, and scoliosis.

We report on an adolescent girl with sparse scalp hair, wide columella extending below alae nasi, webbing at elbows, broad finger tips, short distal phalanx of fingers, swan neck deformity of fingers, scoliosis, tall vertebrae, short fibulae, short fourth metatarsal bone, abnormal distal humeri, and unilateral clubfoot at birth. The combination of these features represents a novel phenotype. We sequenced the protein-coding regions of the FLNA and FLNB genes and did not observe any pathogenic sequence variation. Chromosomal microarray revealed a de novo copy number variation of uncertain clinical significance on 7p22.3.

Nuclear Mitochondrial Disorder Due to a Variant in NAXE in Two Unrelated Indian Children.

Progressive encephalopathy with brain edema and/or leukoencephalopathy type 1 (PEBEL1) is a nuclear mitochondrial disorder involving the NAD(P)HX repair mechanism due to a NAXE variation. PEBEL1 is characterized by rapid neurologic deterioration culminating in death following high-grade fever during infancy. Currently, 23 patients from 14 families are described in the literature, with only three survivors. The authors report two living children from unrelated families with PEBEL1. Both children presented in infancy with ptosis, squint, and ataxia with no skin manifestations. Whole-exome sequencing revealed previously reported c.804_807delInsA (p.Lys270del) variation in exon 6 of NAXE. This is the first Indian report of PEBEL1.

Thirteen Indians with camptodactyly-arthropathy-coxa vara-pericarditis syndrome.

Camptodactyly-arthropathy-coxa vara-pericarditis (CACP) syndrome (MIM# 208250) is a rare monogenic disorder, characterized by early onset of camptodactyly, progressive coxa vara, bilateral arthropathy and constrictive pericarditis. The syndrome is caused by biallelic loss-of-function variants in PRG4. Deficiency of PRG4 results in progressive worsening of joint deformity with age. Thirteen individuals with CACP syndrome from eight consanguineous Indian families were evaluated. We used exome sequencing to elucidate disease-causing variants in all the probands. These variants were further validated and segregated by Sanger sequencing, confirming the diagnosis of CACP syndrome in them. Seven females and six males aged 2-23 years were studied. Camptodactyly (13/13), coxa vara (11/13), short femoral neck (11/13) and arthritis in large joints (12/13) [wrists (11/13), ankle (11/13), elbow (10/13) and knee (10/13)] were observed commonly. Five novel disease-causing variants (c.3636G>T, c.1935del, c.1134dup, c.1699del and c.962T>A) and two previously reported variants (c.1910_1911del and c.2816_2817del) were identified in homozygous state in PRG4. We describe the phenotype and mutations in one of the large cohorts of patients with CACP syndrome, from India.

Mutant MESD links cellular stress to type I collagen aggregation in osteogenesis imperfecta type XX.

Aberrant forms of endoplasmic reticulum (ER)-resident chaperones are implicated in loss of protein quality control in rare diseases. Here we report a novel mutation (p.Asp233Asn) in the ER retention signal of MESD by whole exome sequencing of an individual diagnosed with osteogenesis imperfecta (OI) type XX. While MESDD233N has similar stability and chaperone activity as wild-type MESD, its mislocalization to cytoplasm leads to imbalance of ER proteostasis, resulting in improper folding and aggregation of proteins, including LRP5 and type I collagen. Aggregated LRP5 loses its plasma membrane localization to disrupt the expression of WNT-responsive genes, such as BMP2, BMP4, in proband fibroblasts. We show that MESD is a direct chaperone of pro-α1(I) [COL1A1], and absence of MESDD233N in ER results in cytosolic type I collagen aggregates that remain mostly not secreted. While cytosolic type I collagen aggregates block the intercellular nanotubes, decreased extracellular type I collagen also results in loss of interaction of ITGB1 with type I collagen and weaker attachment of fibroblasts to matrix. Although proband fibroblasts show increased autophagy to degrade the aggregated type I collagen, an overall cellular stress overwhelms the proband fibroblasts. In summary, we present an essential chaperone function of MESD for LRP5 and type I collagen and demonstrating how the D233N mutation in MESD correlates with impaired WNT signaling and proteostasis in OI.

Deep intronic mutation in CRTAP results in unstable isoforms of the protein to induce type I collagen aggregation in a lethal type of osteogenesis imperfecta type VII.

Genetic mutations are involved in Mendelian disorders. Unbuffered intronic mutations in gene variants can generate aberrant splice sites in mutant transcripts, resulting in mutant isoforms of proteins with modulated expression, stability, and function in diseased cells. Here, we identify a deep intronic variant, c.794_1403A>G, in CRTAP by genome sequencing of a male fetus with osteogenesis imperfecta (OI) type VII. The mutation introduces cryptic splice sites in intron-3 of CRTAP, resulting in two mature mutant transcripts with cryptic exons. While transcript-1 translates to a truncated isoform (277 amino acids) with thirteen C-terminal non-wild-type amino acids, transcript-2 translates to a wild-type protein sequence, except that this isoform contains an in-frame fusion of non-wild-type twenty-five amino acids in a tetratricopeptide repeat sequence. Both mutant isoforms of CRTAP are unstable due to the presence of a unique 'GWxxI' degron, which finally leads to loss of proline hydroxylation and aggregation of type I collagen. Although type I collagen aggregates undergo autophagy, the overall proteotoxicity resulted in death of the proband cells by senescence. In summary, we present a genetic disease pathomechanism by linking a novel deep intronic mutation in CRTAP to unstable mutant isoforms of the protein in lethal OI type VII.