Development, validation and application of single molecule molecular inversion probe based novel integrated genetic screening method for 29 common lysosomal storage disorders in India.

BACKGROUND: Current clinical diagnosis pathway for lysosomal storage disorders (LSDs) involves sequential biochemical enzymatic tests followed by DNA sequencing, which is iterative, has low diagnostic yield and is costly due to overlapping clinical presentations. Here, we describe a novel low-cost and high-throughput sequencing assay using single-molecule molecular inversion probes (smMIPs) to screen for causative single nucleotide variants (SNVs) and copy number variants (CNVs) in genes associated with 29 common LSDs in India. RESULTS: 903 smMIPs were designed to target exon and exon-intron boundaries of targeted genes (n = 23; 53.7 kb of the human genome) and were equimolarly pooled to create a sequencing library. After extensive validation in a cohort of 50 patients, we screened 300 patients with either biochemical diagnosis (n = 187) or clinical suspicion (n = 113) of LSDs. A diagnostic yield of 83.4% was observed in patients with prior biochemical diagnosis of LSD. Furthermore, diagnostic yield of 73.9% (n = 54/73) was observed in patients with high clinical suspicion of LSD in contrast with 2.4% (n = 1/40) in patients with low clinical suspicion of LSD. In addition to detecting SNVs, the assay could detect single and multi-exon copy number variants with high confidence. Critically, Niemann-Pick disease type C and neuronal ceroid lipofuscinosis-6 diseases for which biochemical testing is unavailable, could be diagnosed using our assay. Lastly, we observed a non-inferior performance of the assay in DNA extracted from dried blood spots in comparison with whole blood. CONCLUSION: We developed a flexible and scalable assay to reliably detect genetic causes of 29 common LSDs in India. The assay consolidates the detection of multiple variant types in multiple sample types while having improved diagnostic yield at same or lower cost compared to current clinical paradigm.

DISP1 deficiency: Monoallelic and biallelic variants cause a spectrum of midline craniofacial malformations.

PURPOSE: DISP1 encodes a transmembrane protein that regulates the secretion of the morphogen, Sonic hedgehog, a deficiency of which is a major cause of holoprosencephaly (HPE). This disorder covers a spectrum of brain and midline craniofacial malformations. The objective of the present study was to better delineate the clinical phenotypes associated with division transporter dispatched-1 (DISP1) variants. METHODS: This study was based on the identification of at least 1 pathogenic variant of the DISP1 gene in individuals for whom detailed clinical data were available. RESULTS: A total of 23 DISP1 variants were identified in heterozygous, compound heterozygous or homozygous states in 25 individuals with midline craniofacial defects. Most cases were minor forms of HPE, with craniofacial features such as orofacial cleft, solitary median maxillary central incisor, and congenital nasal pyriform aperture stenosis. These individuals had either monoallelic loss-of-function variants or biallelic missense variants in DISP1. In individuals with severe HPE, the DISP1 variants were commonly found associated with a variant in another HPE-linked gene (ie, oligogenic inheritance). CONCLUSION: The genetic findings we have acquired demonstrate a significant involvement of DISP1 variants in the phenotypic spectrum of midline defects. This underlines its importance as a crucial element in the efficient secretion of Sonic hedgehog. We also demonstrated that the very rare solitary median maxillary central incisor and congenital nasal pyriform aperture stenosis combination is part of the DISP1-related phenotype. The present study highlights the clinical risks to be flagged up during genetic counseling after the discovery of a pathogenic DISP1 variant.

Short stature and dysmorphic features in Asian Indian siblings with DAAM2-associated steroid-resistant nephrotic syndrome: Expansion of the phenotypic spectrum or a blended phenotype?

Variants in more than 60 different genes, most of which code for podocyte-related proteins, have been found to be associated with monogenic forms of nephrotic syndrome (NS). Biallelic variants in DAAM2, a member of the formin family, were recently identified to cause autosomal recessive (AR) NS type 24 in four unrelated families with steroid-resistant nephrotic syndrome (SRNS). This case report represents only the fifth reported family of DAAM2-associated NS and the first from India, with two sibs who presented with a complex phenotype characterized by steroid-resistant nephrotic syndrome, short stature, dysmorphic facial features, deep-set toenails, myopia, increased thickness of the calvarium of the skull, and sloping ribs. Both sibs were found to have a homozygous likely pathogenic nonsense variant c.196C>T (p.Arg66Ter; NM_001201427.2) in exon 3 of the DAAM2 gene through whole exome sequencing. The dysmorphic features could possibly be part of the DAAM2-related phenotype which has hitherto not been reported or could represent a blended phenotype, with the extrarenal manifestations resulting from a yet to be identified coexisting genetic condition.

SERPINA11 related novel serpinopathy - A perinatal lethal disorder.

SERPINA11 is a hitherto poorly characterised gene belonging to Clade A of the SERPIN superfamily, with unknown expression pattern and functional significance. We report a perinatal lethal phenotype in two foetuses from the same family associated with a biallelic loss of function variant in SERPINA11, and provide functional evidence to support its candidature as a Mendelian disorder. The SERPINA11 variant-associated foetal phenotype is characterised by gross and histopathological features of extracellular matrix disruption. Western blot and immunofluorescence analyses revealed SERPINA11 expression in multiple mouse tissues, with pronounced expression in the bronchiolar epithelium. We observed a significant decrease in SERPINA11 immunofluorescence in the affected foetal lung compared with a healthy gestation-matched foetus. Protein expression data from HEK293T cell lines following site-directed mutagenesis support the loss of function nature of the variant. Transcriptome analysis from the affected foetal liver indicated the possibility of reduced SERPINA11 transcript abundance. This novel serpinopathy appears to be a consequence of the loss of inhibition of serine proteases involved in extracellular matrix remodelling, revealing SERPINA11 as a protease inhibitor critical for embryonic development.

A unique case of hyperammonemia due to CA5A deficiency: Impact of coexisting gene mutations, pseudogene, and microdeletion.

Carbonic anhydrase 5A (CA5A) belongs to a family of carbonic anhydrases which are zinc metalloenzymes involved in the reversible hydration of CO2 to bicarbonate. Mutations in CA5A are very rare and known to cause Carbonic anhydrase 5A deficiency (CA5AD), an autosomal recessive inborn error of metabolism characterized clinically by acute onset of encephalopathy in infancy or early childhood. CA5A also has two very identical pseudogenes whose interference may result in compromised accuracy in targeted sequencing. We report a unique case of CA5AD caused by compound heterozygous variant (NM_001739.2: c.721G>A: p.Glu241Lys & NM_001739.2: c.619-3420_c.774 + 502del4078bp) in an infant in order to expand the phenotypic spectrum and underscore the impact of pseudogenes, which can introduce complexities in molecular genetic analysis.

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.

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.

A novel homozygous synonymous splicing variant in SELENOI gene causes spastic paraplegia 81.

BACKGROUND: Hereditary spastic paraplegia 81 is a recently identified, rare autosomal recessive disease, caused by biallelic pathogenic variants in the SELENOI gene, with only two families reported to date. The features documented in the two previous affected families include sensorineural deafness, blindness, cleft palate, delayed motor development, regression of motor skills, impaired intellectual development, poor speech and language acquisition, spasticity, hyperreflexia, white matter abnormalities and cerebral and cerebellar atrophy. METHODS: In the present study, we performed exome sequencing analysis in a single family with two affected siblings to identify the genetic cause of complicated hereditary spastic paraplegia. The results were further confirmed by Sanger sequencing, cDNA analysis and 3D protein modelling. RESULTS: Exome sequencing identified a homozygous, synonymous variant in the SELENOI gene (NM_033505.4:c.126G>A:p.(Lys42Lys)) in both of the siblings. Sanger sequencing confirmed the heterozygous status in both parents consistent with the autosomal recessive inheritance. This variant has been found to disrupt normal splicing and lead to skipping of exon 2, causing in-frame deletion of SELENOI N-terminal 23 amino acids [NM_033505.4:c.57_126del:p.(Tyr20_Lys42del)] and further leading to structural changes in the protein. CONCLUSIONS: We report a novel homozygous synonymous variant in the SELENOI gene causing abnormal splicing in two patients affected with hereditary spastic paraplegia 81. This report further expands the phenotypic and genotypic spectrum of hereditary spastic paraplegia 81.

Monoallelic loss-of-function BMP2 variants result in BMP2-related skeletal dysplasia spectrum.

PURPOSE: Bone morphogenic proteins (BMPs) regulate gene expression that is related to many critical developmental processes, including osteogenesis for which they are named. In addition, BMP2 is widely expressed in cells of mesenchymal origin, including bone, cartilage, skeletal and cardiac muscle, and adipose tissue. It also participates in neurodevelopment by inducing differentiation of neural stem cells. In humans, BMP2 variants result in a multiple congenital anomaly syndrome through a haploinsufficiency mechanism. We sought to expand the phenotypic spectrum and highlight phenotypes of patients harboring monoallelic missense variants in BMP2. METHODS: We used retrospective chart review to examine phenotypes from an international cohort of 18 individuals and compared these with published cases. Patient-derived missense variants were modeled in zebrafish to examine their effect on the ability of bmp2b to promote embryonic ventralization. RESULTS: The presented cases recapitulated existing descriptions of BMP2-related disorders, including craniofacial, cardiac, and skeletal anomalies and exhibit a wide phenotypic spectrum. We also identified patients with neural tube defects, structural brain anomalies, and endocrinopathies. Missense variants modeled in zebrafish resulted in loss of protein function. CONCLUSION: We use this expansion of reported phenotypes to suggest multidisciplinary medical monitoring and management of patients with BMP2-related skeletal dysplasia spectrum.

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.

Muscle spasms as presenting feature of Nivelon-Nivelon-Mabile syndrome.

Hedgehog acyltransferase gene (HHAT)-associated Nivelon-Nivelon-Mabile syndrome (NNMS) is a rare genetic disorder of multiple system involvement with microcephaly, central nervous system malformations, skeletal dysplasia, and 46,XY sex reversal. Other variable and inconsistent features reported in this disorder are muscle spasms, facial dysmorphism, prenatal onset growth restriction, microphthalmia, and holoprosencephaly. This is the sixth postnatal reported patient with biallelic variants in HHAT gene, who presented with microcephaly, short stature, muscle hypertrophy, muscle spasms, and facial dysmorphism. The most prominent and presenting finding in this patient were muscle hypertrophy and muscle spasms which had a clinical response to phenytoin and acetazolamide treatment. Our report emphasizes the phenotypic variability of NNMS and further reiterates muscle spasms as an important clinical manifestation of this extremely rare condition.

Long-range PCR amplification-based targeted enrichment & next generation sequencing: A cost-effective testing strategy for lysosomal storage disorders.

Background & objectives: Lysosomal storage disorders (LSDs) are genetic metabolic disorders which result from deficiency of lysosomal enzymes or defects in other lysosomal components. Molecular genetic testing of LSDs is required for diagnostic confirmation when lysosomal enzyme assays are not available or not feasible to perform, and for the identification of the disease causing genetic variants. The aim of this study was to develop a cost-effective, readily customizable and scalable molecular genetic testing strategy for LSDs. Methods: A testing method was designed based on the in-house creation of selective amplicons through long range PCR amplification for targeted capture and enrichment of different LSD genes of interest, followed by next generation sequencing of pooled samples. Results: In the first phase of the study, standardization and validation of the study protocol were done using 28 samples of affected probands and/or carrier parents (group A) with previously identified variants in seven genes, and in the second phase of the study, 30 samples of enzymatically confirmed or biopsy-proven patients with LSDs and/or their carrier parents who had not undergone any prior mutation analysis (group B) were tested and the sequence variants identified in them through the study method were validated by targeted Sanger sequencing. Interpretation & conclusions: This testing approach was found to be reliable, easily customizable and cost-effective for the molecular genetic evaluation of LSDs. The same strategy may be applicable, especially in resource poor settings, for developing cost-effective multigene panel tests for other conditions with genetic heterogeneity.

Next Generation Sequencing in a Case of Early Onset Hydrops: Closing the Loop on the Diagnostic Odyssey!

Non-immune fetal hydrops (NIFH) is an etiologically heterogeneous condition. Cardiac anomalies are one of the common causes of NIFH. Cardiac anomalies can be isolated, multifactorial malformations or have a genetic basis. PLD1 variants have been associated with developmental defects involving the right heart. We present a NIFH with a PLD1 associated right heart malformation.We describe a spontaneously aborted 14 weeks old NIFH fetus with a rudimentary right ventricle, pulmonary valve atresia and pulmonary artery stenosis found at fetopsy. After a normal microarray, whole exome sequencing revealed a homozygous missense variant c.2023 C > T (p. Arg675Trp) in the PLD1 gene. Conclusion: Detailed fetopsy and genetic evaluation in this NIFH allowed an etiological explanation, further corroborated the association of PLD1 gene variants and developmental right heart defects, and that this defect can be associated with NIHF.

Microcephalic primordial dwarfism with predominant Meier-Gorlin phenotype, ichthyosis, and multiple joint deformities-Further expansion of DONSON Cell Cycle-opathy phenotypic spectrum.

We report a patient with microcephalic primordial dwarfism with predominant Meier-Gorlin syndrome phenotype with ichthyosis and disabling multiple joint deformities in addition to classic features of the syndrome. The patient was a 10.5-year-old girl referred in view of short stature, joint deformities, and facial dysmorphism. There was history of intrauterine growth restriction and collodion like skin abnormality at birth. She had normal developmental milestones and intellect. On clinical evaluation, anthropometry was suggestive of proportionate short stature and microcephaly. There was abnormal posture due to spine and peripheral joint deformities, along with ichthyosis, facial, and digital dysmorphism. Skeletal radiographs showed radial subluxation, acetabular dysplasia and hip dislocation, bilateral knee joint dislocation, absent patellae, slender long bones with delayed bone age, and subluxation of small joints of hands and feet. Work up for metabolic bone disease and peripheral blood karyotype was normal. Whole exome sequencing revealed a pathogenic homozygous variant c.C1297T (p.Pro433Ser) in the exon 8 of DONSON gene. This report further expands the genotypic-phenotypic spectrum of the group of disorders known as Cell Cycle-opathies.

Fetal presentation of chondrodysplasia with joint dislocations, GPAPP type, caused by novel biallelic IMPAD1 variants.

Biallelic IMPAD1 pathogenic variants leads to deficiency of GPAPP (Golgi 3-prime phosphoadenosine 5-prime phosphate 3-prime phosphatase) protein and clinically causes chondrodysplasia, which is characterized by short stature with short limbs, craniofacial malformations, cleft palate, hand and foot anomalies, and various radiographic skeletal manifestations. Here we describe prenatal presentation of GPAPP deficiency caused by novel biallelic pathogenic variants, 2 base pair duplication in exon 2 of IMAPD1 gene in a patient of Asian-Indian origin. Further we report on diagnostic clues of prenatal presentation of GPAPP deficiency through ultrasonography, fetal MRI, and postmortem findings. We also provide evidence of pathophysiology of underlying GPAPP deficiency in the form of disorganization and dysplastic chondrocytes and reduced sulfation of glycoproteins through histopathology of cartilage similar to that described in mice IMPAD1 homozygous mutant model.

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.

Prenatal phenotype of FBXL4-associated encephalomyopathic mitochondrial DNA depletion syndrome-13.

FBXL4 -associated encephalomyopathic mitochondrial DNA depletion syndrome-13 (MTDPS13) is a rare genetic disorder characterized by early neonatal onset of encephalopathy, seizures, lactic acidosis, hypotonia, dysmorphism, and severe global developmental delay. Prenatal phenotype of molecularly confirmed MTDPS13 has not been well studied. This is the case report of a non-consanguineously conceived fetus ascertained first at 20 weeks of gestation with multiple soft markers. Follow-up fetal ultrasonogram at 26 weeks revealed periventricular cysts, periventricular echogenicity, ventriculomegaly, thin corpus callosum, mega cisterna magna, and large cavum. Fetal MRI confirmed these findings. Postnatally, the baby had clinical and biochemical findings indicative of a mitochondriopathy and died on neonatal day 3. Whole exome sequencing on stored amniotic fluid DNA confirmed the diagnosis of encephalomyopathic mitochondrial DNA depletion syndrome-13 (MTDPS13). This report presents the prenatal phenotype of this rare mitochondriopathy, which has been recognized primarily in postnatal patients. The brain imaging findings in the reported fetus indicate that MTDPS13 is associated with progressive neurological involvement and brain tissue destructive changes starting as early as the second trimester of pregnancy. The case also raises concerns regarding the association of so-called soft markers, which were the only initial finding in this case, with severe monogenic diseases.

Fetal phenotypes of Mendelian disorders: A descriptive study from India.

OBJECTIVE: Exome sequencing (ES)-based diagnosis of Mendelian diseases in the fetus is limited by paucity of phenotypic information. This study reports the comprehensive phenotypes of some fetuses with Mendelian disorders. METHODS: Next generation technology-based sequencing of all coding regions of the genome (Exome sequencing) or targeted gene sequencing using Sanger or next generation platforms was performed in a cohort of deeply phenotyped, cytogenetically normal fetuses with morphological defects. Prenatal ultrasonographic phenotypes and postmortem details including dysmorphology, histopathology, and radiography were ascertained. Novel candidate genes, novel/unusual findings, and unusual genotypes in cases with confirmed Mendelian disorders are described. RESULTS: Of the 102 fetuses sequenced, 45 (44%) achieved definitive diagnosis of a Mendelian disorder with 50 pathogenic/likely pathogenic variants. The majority (87%) were autosomal recessive, 69% families were consanguineous, and 54% variants were novel. Dysmorphic syndromes, skeletal dysplasias, and metabolic disorders were the commonest disease categories, ciliopathies and dystroglycanopathies, commonest molecular categories. We describe the first fetal description of six monogenic diseases, and nine cases with novel histological findings. Nineteen cases had novel/unusual findings. CONCLUSION: This cohort demonstrates how deep fetal phenotypes of some Mendelian disorders can show novel/unusual findings, which have important implications for prenatal diagnosis of these conditions.

A data set of variants derived from 1455 clinical and research exomes is efficient in variant prioritization for early-onset monogenic disorders in Indians.

Given the genomic uniqueness, a local data set is most desired for Indians, who are underrepresented in existing public databases. We hypothesize patients with rare monogenic disorders and their family members can provide a reliable source of common variants in the population. Exome sequencing (ES) data from families with rare Mendelian disorders was aggregated from five centers in India. The dataset was refined by excluding related individuals and removing the disease-causing variants (refined cohort). The efficiency of these data sets was assessed in a new set of 50 exomes against gnomAD and GenomeAsia. Our original cohort comprised 1455 individuals from 1203 families. The refined cohort had 836 unrelated individuals that retained 1,251,064 variants with 181,125 population-specific and 489,618 common variants. The allele frequencies from our cohort helped to define 97,609 rare variants in gnomAD and 44,520 rare variants in GenomeAsia as common variants in our population. Our variant dataset provided an additional 1.7% and 0.1% efficiency for prioritizing heterozygous and homozygous variants respectively for rare monogenic disorders. We observed additional 19 genes/human knockouts. We list carrier frequency for 142 recessive disorders. This is a large and useful resource of exonic variants for Indians. Despite limitations, datasets from patients are efficient tools for variant prioritization in a resource-limited setting.

Functional characterization of novel variants in SMPD1 in Indian patients with acid sphingomyelinase deficiency.

Pathogenic variations in SMPD1 lead to acid sphingomyelinase deficiency (ASMD), that is, Niemann-Pick disease (NPD) type A and B (NPA, NPB), which is a recessive lysosomal storage disease. The knowledge of variant spectrum in Indian patients is crucial for early and accurate NPD diagnosis and genetic counseling of families. In this study, we recruited 40 unrelated pediatric patients manifesting symptoms of ASMD and subnormal ASM enzyme activity. Variations in SMPD1 were studied using Sanger sequencing for all exons, followed by interpretation of variants based on American College of Medical Genetics and Genomics & Association for Molecular Pathology (ACMG/AMP) criteria. We identified 18 previously unreported variants and 21 known variants, including missense, nonsense, deletions, duplications, and splice site variations with disease-causing potential. Eight missense variants were functionally characterized using in silico molecular dynamic simulation and in vitro transient transfection in HEK293T cells, followed by ASM enzyme assay, immunoblot, and immunofluorescence studies. All the variants showed reduced ASM activity in transfected cells confirming their disease-causing potential. The study provides data for efficient prenatal diagnosis and genetic counseling of families with NPD type A and B.