- Open Access
Deletion mutation in BSCL2 gene underlies congenital generalized lipodystrophy in a Pakistani family
Diagnostic Pathology volume 8, Article number: 78 (2013)
Congenital generalized lipodystrophy (CGL) also known as Berardinelli-Seip Congenital Lipodystrophy (BSCL) is a genetically heterogeneous disorder characterized by loss of adipose tissues, Acanthosis nigricans, diabetes mellitus, muscular hypertrophy, hepatomegaly and hypertriglyceridemia. There are four subclinical phenotypes of CGL (CGL1-4) and mutations in four genes AGPAT2, BSCL2, CAV1 and PTRF have been assigned to each type.
The study included clinical and molecular investigations of CGL disease in a consanguineous Pakistani family. For mutation screening all the coding exons including splice junctions of AGPAT2, BSCL2, CAV1 and PTRF genes were PCR amplified and sequenced directly using an automated DNA sequencer ABI3730.
Sequence analysis revealed a single base pair deletion mutation (c.636delC; p.Tyr213ThrfsX20) in exon 5 of BSCL2 gene causing a frame shift and premature termination codon.
Mutation identified here in BSCL2 gene causing congenital generalized lipodystrophy is the first report in Pakistani population. The patients exhibited characteristic features of generalized lipodystrophy, Acanthosis nigricans, diabetes mellitus and hypertrophic cardiomyopathy.
The virtual slide(s) for this article can be found here: http://www.diagnosticpathology.diagnomx.eu/vs/1913913076864247.
Congenital generalized lipodystrophy (CGL), or Berardinelli-Seip congenital lipodystrophy (BSCL), is a genetically heterogeneous disorder characterized by loss of adipose tissues and sub-cutaneous fats, enlarged fatty liver, hypertrophic muscles, Acanthosis nigricans, increased serum triglyceride level, insulin intolerance or diabetes mellitus. Pathogenicity of four genes has been reported in four clinically overlapping phenotypes so far.
The affected individuals of congenital generalized lipodystrophy type 1 (CGL1, MIM 608594) are characterized by typical poor fat accumulation in the metabolically active and mechanical adipose tissues. In addition these patients have special features of lytic bone lesions that are absent in other forms of CGL [1, 2]. Loss-of-function mutations in 1-acylglycerol-3-phosphate-O-Acyltransferase (AGPAT2, MIM 603100) gene at chromosome 9q34.3 have been assigned to cause CGL1 phenotype [3, 4].
Pathogenic variants of homologous to mouse gamma-3-linked (BSCL2, MIM 606158) gene, located at chromosome 11q13, cause congenital generalized lipodystrophy type 2 (CGL2, MIM 269700) or a seipin-deficient phenotype [3, 5]. The subcutaneous biopsies examinations of patients have revealed scattered groups of small adipocytes with low but detectable lipid content . The seipin-deficient individuals have generalized congenital lipodystrophy, hypertrophic cardiomyopathy, higher rates of mild mental retardation and an earlier onset of diabetes [7–9].
Mutations in caveolin 1 (CAV1, MIM 601047) gene, located at chromosome 7q31.3, cause congenital generalized lipodystrophy type 3 (CGL3, MIM 612526). Caveolin-1-deficient patients exhibit unique manifestation of short stature, hypocalcemia and vitamin D resistance . Recently another form congenital generalized lipodystrophy type 4 (CGL4, MIM 613327) has also been described caused by mutations in RNA polymerase 1 and transcript release factor (PTRF, MIM 603198) gene at chromosome 17q21 .
In the present study clinical and molecular analysis of a four generations consanguineous Pakistani family demonstrating autosomal recessive CGL phenotype was performed. Direct sequencing of four candidate genes in a Pakistani family revealed a single base pair deletion mutation in BSCL2 gene.
The present four-generation family, exhibiting features of CGL or BSCL (Figure 1), belonged to Peshawar city of Khyber Pakhtunkhwa province, Pakistan. The affected individuals were born to first degree cousins, suggesting that affected individuals are homozygous for a mutant allele. All affected and unaffected individuals underwent examination at the local government hospitals. Prior to start of study, approval was obtained from Institutional Review Board (IRB), Khyber Medical University Peshawar, Pakistan. Informed consent for the study including presentation of photographs was obtained from affected individuals and their parents for publication. Genomic DNA extraction and polymerase chain reaction was performed as described earlier .
Entire coding region and splice junction sites of AGPAT2, BSCL2, CAV1 and PTRF genes were amplified by PCR and screened by DNA sequencing for potential sequence variants. Primer sequences (available on request) were designed for each exon using Primer3 software  and checked for specificity using Basic Local Alignment Search Tool (BLAST; http://www.ncbi.nlm.nih.gov/blast). Purification of the PCR-amplified DNA was performed with commercially available kits (Marligen Biosciences, Rockville MD, USA). Sequencing of the potential candidate genes was performed using Big Dye Terminator v3.1 Cycle Sequencing Kit, together with an ABI Prism 3730 Genetic Analyzer (Applera, Foster City, CA, USA). Sequence variants were identified via Bioedit sequence alignment editor version 6.0.7 (http://www.mbio.ncsu.edu/BioEdit/bioedit.html).
Clinical features of the affected individuals
The affected members (IV-1 and IV-3) of the family showed characteristic features of congenital generalized lipodystrophy. Both the affected individuals had acanthosis nigricans (velvety thickening and hyperpigmentation of the skin) which was more prominent around neck and in body folds including axillae, anticubital fossae and popliteal fossae. Both the affected individuals had prominent veins, rough dry skin and umbilical protrusion. Curly and dry hair was present with mild hypertrichosis over the scalp (Figure 2a).
Muscular hypertrophy was observed in skeletal muscles more prominently at arms and shin areas in both the affected individuals (Figure 2b). Radiological examinations of bones in the affected individuals showed relatively hyper-density of metacarpals, phalanges, sclerotic hip bones and adjacent portion of both femurs. Slight hyper-density was observed in skull bone, however, vault thickness and para-nasal sinuses were normal in both the affected individuals. Chest X-rays revealed normal lung fields, normal cardiac shadow with slight hyper-density of scapulae, humerus and clavicles. Magnetic Resonance Imaging examination (MRI) of the brain does not reveal generalized brain atrophy (Figure 2c).
Abdominal and cardiac examinations
Ultrasonography of the abdomen revealed moderate hepatomegaly with mild splenomegaly in both the affected individuals. There was no symptom of polycystic ovarian syndrome in 7 years affected female (VI-3), however her 11 years affected brother (VI-1) had moderate genital enlargement. Both the patients exhibited hypertrophic cardiomyopathy. Electrocardiography (ECG) of both the affected individuals showed sinus rhythm and right bundle branch pattern. Echocardiography of IV-1 revealed patent foramen ovale and mild left ventricular hypertrophy with intact inter-ventricular septum, while IV-3 showed small secundum atrial septal defect with left to right shunt, minimal right ventricular dilatation, and an intact inter-ventricular septum. There was no significant atrio-ventricular valve re-gurge. The left ventricular function was found normal (data not shown).
The total serum bilirubin and serum electrolytes levels of both the affected individuals were within the normal range. Serum glutamate pyruvate, blood sugar, alkaline phosphatase and triglyceride levels were raised in both the affected individuals. High density lipoprotein levels were low in both the affected individuals than normal in both of them. The glycosylated hemoglobin (HbA1C) level of the affected individual IV-1 was slight higher than controlled diabetic range while that of IV-3 was in the non-diabetic level (Table 1).
The patients did not cooperate for tissue biopsy therefore histo-pathological examinations of the skin and sural nerves were not performed.
To search for potential sequence variants, all coding exons and splice junction sites of AGPAT2, BSCL2, CAV1 and PTRF genes were sequenced initially in two affected and one unaffected individuals of the family. A homozygous deletion mutation of a single base cytosine at complementary DNA position 636 (c.636delC) was detected in exon 5 of BSCL2 gene in both the affected individuals (Figure 3). This deletion probably shifted the reading frame leading to a premature stop codon and adding 20 non-specific amino acid residues BSCL2 protein (p.Tyr213ThrfsX20). Mutation analysis of BSCL2 exon 5 revealed c.636del in heterozygous state in obligate carriers and phenotypically unaffected individuals of the family. To ensure that c.636delC does not represent a neutral polymorphism in this population, a panel of 100 unrelated and ethnically matched control individuals was screened for this mutation, thus confirming that mutation was not present outside the family.
Loss-of-function mutations in the BSCL2 gene cause a severe form of lipodystrophy, whilst characteristic gain-of-function mutations are believed to be associated with aggregation of unfolded protein in endoplasmic reticulum resulting in neurodegeneration leading to a heterogeneous group of neuropathies . The study presented here was a four generations consanguineous family, originated from Peshawar city Pakistan. Clinical features of the affected individuals resembled CGL2 phenotype. Moreover the first degree consanguinity of the parents and autosomal recessive mode of inheritance of the disease phenotype was a definite clue for a homozygous mutation. Mutations in BSCL2 leading to CGL2 phenotype have been identified worldwide [5, 14–16].
The generalized loss of adipose tissue, increased triglyceride levels and steatosis of the liver were comparable to a homozygous mutation p.Tyr213ThrfsX20 identified in an Indian family . Early onset diabetes mellitus present in our patients were comparable with homozygous nonsense mutations identified in Chinese and severe insulin resistance in Japanese and Brazilian patients [9, 14, 15], however we did not perform the insulin resistance test in our patients. The hypertrophic cardiomyopathy identified in our cases was less severe as compared to Chinese patient caused by homozygous nonsense mutation . There was mild mental retardation (IQ score 65–75) in both of our patients but the brain MRI did not reveal brain atrophy.
The variable features of neuropathies are mostly associated with heterozygous mutations in BSCL2 dHMN, Charcot-Marie-Tooth (CMT) and Silver syndromes [16, 18, 19]. In 2005, based on disease allele penetrance and severity of phenotype caused by heterozygous mutations in BSCL2 gene p.N88S in Austrian and German families, these patients were classified into six subclinical groups with overlapping features of HMN, CMT and Silver syndrome . However, in 2009, Brusse et al. discovered a digenic inheritance of hereditary motor neuropathy in a large Dutch family with heterozygous mutation p.N88S in the BSCL2 gene and segregating autosomal dominant disease haplotype at chromosome 16p . The phenotypic variability ranged from strictly neuropathic weakness to a spastic paraplegia with hereditary motor neuropathy presenting clinical phenotype similar to Silver syndrome . Very recently another heterozygous mutation p.S90W identified in two Korean CMT type 2 disease patients  was associated with increased density of myelinated fibers. However, this feature was slightly different from the previously reported mutation p.S90L in three Italian patients representing CMT type 2 phenotype with pyramidal signs and subclinical sensory involvement on sural nerve biopsy .
Human BSCL2 gene encodes 398 or 462 amino acids seipin protein from either of the three transcripts 1.6 kb, 1.8 kb and 2.2 kb . The 1.8 kb transcript is exclusively expressed in brain and testis while the other two are ubiquitously expressed . Seipin protein is located in endoplasmic reticulum acting as a regulator of lipid catabolism and is essential for the differentiation of fat cells [18, 24]. Studies have revealed role of seipin in proper lipid storage and regulation of cAMP/PKA-mediated lipolysis in adipose differentiation [25, 26]. In adult mice brain bscl2 expression studies have suggested the possible involvement of seipin in central regulation of energy balance.
According to UniprotKB database (http://www.uniprot.org/uniprot/Q96G97), 398 amino acids seipin protein has five domains including 2 cytoplasmic [1-26aa and 264-398aa], 2 transmembrane [27-47aa and 243-263aa] and a luminal [48-242aa]. So far 31 mutations are reported in BSCL2 gene including 12 missense/nonsense, 4 small insertions, 4 small deletions, 7 splice-site mutations and 4 complex rearrangements causing related phenotypes as summarized in Table 2.
More than three hundred cases of BSCL have been reported in the medical literature with an estimated prevalence of 1 in 10 million people in USA . However this condition is more common in other populations around the world like Lebanon, Brazil, Portugal and Sultanate of Oman with an estimated prevalence of 1:25000 to 1:1000000 . The incidence data of BSCL and many other rare disorders are not available from Pakistan partly due to lack of disease registry database systems and non-availability of neonatal screening programs in the hospitals.
Generally, in families segregating autosomal recessive disorders, high frequency of consanguineous marriages may increase the frequency of homozygotes in the population leading to increased incidence of certain lesions, their founder effect and also the appearance of new mutations . Although most causative mutations in Mendelian diseases are reported in single families, certain mutations may occur more frequently in some populations than others. Few examples of such prevalent mutations may be observed in Wilson’s disease , hemophilia , Hurler disease  and thalassemia . Knowledge of the differences in the worldwide distribution of particular mutations may help to design shortcuts for genetic diagnosis and screening of relevant inherited diseases.
The identified mutation c.636delC (p.Tyr213ThrfsX20) in exon 5 of BSCL2 gene is the first report of CGL2 from Pakistan. The current study extends the body of evidence that describes the role of BSCL2 in congenital generalized lipodystrophy and association of loss-of-function mutations with severe lipodystrophy phenotype. The clinical features in our patients are similar to CGL2 phenotype but lack of motor neuropathies that are usually associated with Silver syndrome and Charcot-Marie-Tooth disease were not observed in our case. Understanding the molecular defects of BSCL2 gene in the patients can be helpful in the genetic counseling and prenatal diagnosis of affected families and may help to improve specific therapeutic interventions.
Van Maldergem L, Magre J, Khallouf TE, Gedde-Dahl T, Delepine M, Trygstad O, Seemanova E, Stephenson T, Albott CS, Bonnici F: Genotype-phenotype relationships in Berardinelli-Seip congenital lipodystrophy. J Med Genet. 2002, 39 (10): 722-733.
Simha V, Garg A: Phenotypic heterogeneity in body fat distribution in patients with congenital generalized lipodystrophy caused by mutations in the AGPAT2 or seipin genes. J Clin Endocrinol Metab. 2003, 88 (11): 5433-5437.
Agarwal AK, Arioglu E, De Almeida S, Akkoc N, Taylor SI, Bowcock AM, Barnes RI, Garg A: AGPAT2 is mutated in congenital generalized lipodystrophy linked to chromosome 9q34. Nat Genet. 2002, 31 (1): 21-23.
Fu M, Kazlauskaite R, Baracho Mde F, Santos MG, Brandao-Neto J, Villares S, Celi FS, Wajchenberg BL, Shuldiner AR: Mutations in Gng3lg and AGPAT2 in Berardinelli-Seip congenital lipodystrophy and Brunzell syndrome: phenotype variability suggests important modifier effects. J Clin Endocrinol Metab. 2004, 89 (6): 2916-2922.
Magre J, Delepine M, Khallouf E, Gedde-Dahl T, Van Maldergem L, Sobel E, Papp J, Meier M, Megarbane A, Bachy A: Identification of the gene altered in Berardinelli-Seip congenital lipodystrophy on chromosome 11q13. Nat Genet. 2001, 28 (4): 365-370.
Rognum TO, Bjerve KS, Seip M, Trygstad O, Oseid S: Fat cell size and lipid content of subcutaneous tissue in congenital generalized lipodystrophy. Acta Endocrinol. 1978, 88 (1): 182-189.
Agarwal AK, Simha V, Oral EA, Moran SA, Gorden P, O'Rahilly S, Zaidi Z, Gurakan F, Arslanian SA, Klar A: Phenotypic and genetic heterogeneity in congenital generalized lipodystrophy. J Clin Endocrinol Metab. 2003, 88 (10): 4840-4847.
Agarwal AK, Barnes RI, Garg A: Genetic basis of congenital generalized lipodystrophy. Int J Obes Relat Metab Disord. 2004, 28 (2): 336-339.
Gomes KB, Pardini VC, Ferreira AC, Fernandes AP: Phenotypic heterogeneity in biochemical parameters correlates with mutations in AGPAT2 or Seipin genes among Berardinelli-Seip congenital lipodystrophy patients. J Inherit Metab Dis. 2005, 28 (6): 1123-1131.
Kim CA, Delepine M, Boutet E, El Mourabit H, Le Lay S, Meier M, Nemani M, Bridel E, Leite CC, Bertola DR: Association of a homozygous nonsense caveolin-1 mutation with Berardinelli-Seip congenital lipodystrophy. J Clin Endocrinol Metab. 2008, 93 (4): 1129-1134.
Hayashi YK, Matsuda C, Ogawa M, Goto K, Tominaga K, Mitsuhashi S, Park YE, Nonaka I, Hino-Fukuyo N, Haginoya K: Human PTRF mutations cause secondary deficiency of caveolins resulting in muscular dystrophy with generalized lipodystrophy. J Clin Invest. 2009, 119 (9): 2623-2633.
Jelani M, Chishti MS, Ahmad W: Mutation in PVRL4 gene encoding nectin-4 underlies ectodermal-dysplasia-syndactyly syndrome (EDSS1). J Hum Genet. 2011, 56 (5): 352-357.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.
Jin J, Cao L, Zhao Z, Shen S, Kiess W, Zhi D, Ye R, Cheng R, Chen L, Yang Y: Novel BSCL2 gene mutation E189X in Chinese congenital generalized lipodystrophy child with early onset diabetes mellitus. Eur J Endocrinol/Eur Fed Endocrine Soc. 2007, 157 (6): 783-787.
Nishiyama A, Yagi M, Awano H, Okizuka Y, Maeda T, Yoshida S, Takeshima Y, Matsuo M: Two Japanese infants with congenital generalized lipodystrophy due to BSCL2 mutations. Pediatr Int. 2009, 51 (6): 775-779.
van de Warrenburg BP, Scheffer H, van Eijk JJ, Versteeg MH, Kremer H, Zwarts MJ, Schelhaas HJ, van Engelen BG: BSCL2 mutations in two Dutch families with overlapping Silver syndrome-distal hereditary motor neuropathy. Neuromuscul Disord. 2006, 16 (2): 122-125.
Friguls B, Coroleu W, del Alcazar R, Hilbert P, Van Maldergem L, Pintos-Morell G: Severe cardiac phenotype of Berardinelli-Seip congenital lipodystrophy in an infant with homozygous E189X BSCL2 mutation. Eur J Med Genet. 2009, 52 (1): 14-16.
Windpassinger C, Auer-Grumbach M, Irobi J, Patel H, Petek E, Horl G, Malli R, Reed JA, Dierick I, Verpoorten N: Heterozygous missense mutations in BSCL2 are associated with distal hereditary motor neuropathy and Silver syndrome. Nat Genet. 2004, 36 (3): 271-276.
Cafforio G, Calabrese R, Morelli N, Mancuso M, Piazza S, Martinuzzi A, Bassi MT, Crippa F, Siciliano G: The first Italian family with evidence of pyramidal impairment as phenotypic manifestation of Silver syndrome BSCL2 gene mutation. Neurol Sci. 2008, 29 (3): 189-191.
Auer-Grumbach M, Schlotter-Weigel B, Lochmuller H, Strobl-Wildemann G, Auer-Grumbach P, Fischer R, Offenbacher H, Zwick EB, Robl T, Hartl G: Phenotypes of the N88S Berardinelli-Seip congenital lipodystrophy 2 mutation. Ann Neurol. 2005, 57 (3): 415-424.
Brusse E, Majoor-Krakauer D, de Graaf BM, Visser GH, Swagemakers S, Boon AJ, Oostra BA, Bertoli-Avella AM: A novel 16p locus associated with BSCL2 hereditary motor neuronopathy: a genetic modifier?. Neurogenetics. 2009, 10 (4): 289-297.
Choi BO, Park MH, Chung KW, Woo HM, Koo H, Chung HK, Choi KG, Park KD, Lee HJ, Hyun YS: Clinical and histopathological study of Charcot-Marie-Tooth neuropathy with a novel S90W mutation in BSCL2. Neurogenetics. 2013, 14 (1): 35-42.
Luigetti M, Fabrizi GM, Madia F, Ferrarini M, Conte A, Delgrande A, Tonali PA, Sabatelli M: Seipin S90L mutation in an Italian family with CMT2/dHMN and pyramidal signs. Muscle Nerve. 2010, 42 (3): 448-451.
Lundin C, Nordstrom R, Wagner K, Windpassinger C, Andersson H, von Heijne G, Nilsson I: Membrane topology of the human seipin protein. FEBS Lett. 2006, 580 (9): 2281-2284.
Boutet E, El Mourabit H, Prot M, Nemani M, Khallouf E, Colard O, Maurice M, Durand-Schneider AM, Chretien Y, Gres S: Seipin deficiency alters fatty acid Delta9 desaturation and lipid droplet formation in Berardinelli-Seip congenital lipodystrophy. Biochimie. 2009, 91 (6): 796-803.
Chen W, Chang B, Saha P, Hartig SM, Li L, Reddy VT, Yang Y, Yechoor V, Mancini MA, Chan L: Berardinelli-seip congenital lipodystrophy 2/seipin is a cell-autonomous regulator of lipolysis essential for adipocyte differentiation. Mol Cell Biol. 2012, 32 (6): 1099-1111.
Miranda DM, Wajchenberg BL, Calsolari MR, Aguiar MJ, Silva JM, Ribeiro MG, Fonseca C, Amaral D, Boson WL, Resende BA: Novel mutations of the BSCL2 and AGPAT2 genes in 10 families with Berardinelli-Seip congenital generalized lipodystrophy syndrome. Clin Endocrinol. 2009, 71 (4): 512-517.
Yagi T, Ito D, Nihei Y, Ishihara T, Suzuki N: N88S seipin mutant transgenic mice develop features of seipinopathy/BSCL2-related motor neuron disease via endoplasmic reticulum stress. Hum Mol Genet. 2011, 20 (19): 3831-3840.
Cho HJ, Sung DH, Ki CS: Identification of de novo BSCL2 Ser90Leu mutation in a Korean family with Silver syndrome and distal hereditary motor neuropathy. Muscle Nerve. 2007, 36 (3): 384-386.
Ito D, Suzuki N: Seipin/BSCL2-related motor neuron disease: Seipinopathy is a novel conformational disease associated with endoplasmic reticulum stress. Clinical neurology. 2007, 47 (6): 329-335.
Ebihara K, Kusakabe T, Masuzaki H, Kobayashi N, Tanaka T, Chusho H, Miyanaga F, Miyazawa T, Hayashi T, Hosoda K: Gene and phenotype analysis of congenital generalized lipodystrophy in Japanese: a novel homozygous nonsense mutation in seipin gene. J Clin Endocrinol Metab. 2004, 89 (5): 2360-2364.
Shirwalkar HU, Patel ZM, Magre J, Hilbert P, Van Maldergem L, Mukhopadhyay RR, Maitra A: Congenital generalized lipodystrophy in an Indian patient with a novel mutation in BSCL2 gene. J Inherit Metab Dis. 2008, 31 (Suppl 2): S317-S322.
Mandal K, Aneja S, Seth A, Khan A: Berardinelli-Seip congenital lipodystrophy. Indian Pediatr. 2006, 43 (5): 440-445.
Jeninga EH, de Vroede M, Hamers N, Breur JM, Verhoeven-Duif NM, Berger R, Kalkhoven E: A Patient with Congenital Generalized Lipodystrophy Due To a Novel Mutation in BSCL2: Indications for Secondary Mitochondrial Dysfunction. JIMD reports. 2012, 4: 47-54.
Agarwal AK, Garg A: Genetic disorders of adipose tissue development, differentiation, and death. Annu Rev Genomics Hum Genet. 2006, 7: 175-199.
Rajab A, Bappal B, Al-Shaikh H, Al-Khusaibi S, Mohammed AJ: Common autosomal recessive diseases in Oman derived from a hospital-based registry. Community Genet. 2005, 8 (1): 27-30.
Chkioua L, Khedhiri S, Ben Turkia H, Chahed H, Ferchichi S, Ben Dridi MF, Laradi S, Miled A: Hurler disease (mucopolysaccharidosis type IH): clinical features and consanguinity in Tunisian population. Diagn Pathol. 2011, 6: 113-
Ferenci P: Regional distribution of mutations of the ATP7B gene in patients with Wilson disease: impact on genetic testing. Hum Genet. 2006, 120 (2): 151-159.
Elmahmoudi H, Khodjet-el-khil H, Wigren E, Jlizi A, Zahra K, Pellechia D, Vinciguerra C, Meddeb B, Elggaaied AB, Gouider E: First report of molecular diagnosis of Tunisian hemophiliacs A: identification of 8 novel causative mutations. Diagn Pathol. 2012, 7: 93-
Khateeb B, Moatter T, Shaghil AM, Haroon S, Kakepoto GN: Genetic diversity of beta-thalassemia mutations in Pakistani population. J Pak Med Assoc. 2000, 50 (9): 293-296.
The authors declare that they have no competing interests.
OR performed the molecular studies, NK and KK performed clinical diagnosis, JA performed radiological examination and data interpretation, MAK performed sequencing alignment and manuscript writing, JYA and MN critically reviewed the manuscript, MJ planned the study and finalized the manuscript. All authors agreed upon the final manuscript.