Variants in FREM1 and trisomy 18 identified in a neonatal progeria patient.

BACKGROUND: Neonatal progeroid disorders are rare disorders with clinical features including low body mass index, proptosis, aged and dysmorphic facial features at the time of birth, prominent veins, sparse scalp hairs, and severe growth retardation. Very few cases have been identified with an unknown genetic cause. Here, we report clinical and genetic findings of a proband with hallmark features of neonatal progeria. METHODS: Microarray comparative genomic hybridization, whole exome sequencing (WES) and Sanger sequencing were performed using standard methods. RESULTS: Array combined genome hybridization data revealed trisomy 18 in the proband (II-1), and WES data identified novel compound heterozygous variants (c.247 C > T; p.H83Y and c.14769868InsA) in the FREM1 gene. CONCLUSION: We report a novel complex case of neonatal progeria with atrial septal defects, trisomy 18 without typical features of Edward syndrome. The phenotype of the patient was more consistent with neonatal progeria, thus we speculate it to be caused by the FREM1 variants.

Antecedent infections, recent developments and future directions in Guillain-Barré syndrome.

The Guillain-Barré syndrome is an autoimmune polyradiculoneuropathy causing symmetrical weakness of limbs. After poliomyelitis, it is the second most common cause of paralysis, with an annual incidence of 0.84-1.91 per 100,000 individuals. The syndrome affects both men and women, showing a male preponderance. Campylobacter jejuni, epstein-barr virus, cytomegalovirus, mycoplasma pneumoniae and haemophilus influenzae are amongst the most common causative agents of Guillain-Barré syndrome. Several immunological and genetic factors have been recognised as the risk factors. Human leukocyte antigen, cluster of differentiation 1, and tumour necrosis factor-alpha alleles are among the frequently investigated loci in Guillain-Barré syndrome. Genome-wide association studies have found no significant association of Guillain-Barré syndrome with common variants. Many vaccines against Campylobacter jejuni infection have been proposed, but there are concerns about the efficacy and safety of these vaccines. So far, there is no approved vaccine against Campylobacter jejuni.

Biallelic variants in CHST3 cause Spondyloepiphyseal dysplasia with joint dislocations in three Pakistani kindreds.

BACKGROUND: Skeletal dysplasia is a heterogeneous group of disorders. Spondyloepiphyseal dysplasias comprise one subgroup. Deficiency of carbohydrate sulfotransferase 3 has been reported in a small number of patients with recessively inherited spondyloepiphyseal dysplasia with joint dislocation, short stature and scoliosis. We report here molecular and clinical findings of affected individuals in three consanguineous Pakistani families. Affected individuals in all three families had a uniform phenotype including severe short stature, multiple dislocated joints, progressive scoliosis and facial dysmorphism. METHODS: Clinical evaluation was done for three unrelated families. Radiological survey of bones was completed for patients from two of the families. Whole exome sequencing index patients from each family was performed followed by Sanger sequencing for validation of segregation of identified variants in respective families. In-silico analysis for determining pathogenicity of identified variants and conservation was done. RESULTS: Whole-exome sequencing revealed biallelic variants c.590 T > C;p.(Leu197Pro), c.603C > A;p.(Tyr201Ter) and c.661C > T;p.(Arg221Cys) in CHST3 (NM_004273.5) in the three families with eight, five and two affected individuals, respectively. Contrary to previous reports, affected individuals in none of the families exhibited a hearing loss. CONCLUSION: We describe genotypic and phenotypic findings of three unrelated families with spondyloepiphyseal dysplasia. Our study confirms phenotypic variability and adds to the genotypic spectrum of spondyloepiphyseal dysplasia.

Study of Guillain-Barre syndrome etiology in Pakistani patients.

OBJECTIVE: To examine clinical features, biochemical markers, demographic features, antecedent infections, frequency and treatment strategies related to Guillain-Barré syndrome. METHODS: The case-control study was conducted at the Pakistan Institute of Medical Sciences, Islamabad, Pakistan, and the District Headquarters Hospital, Rawalpindi, Pakistan, from 2018 to 2020, and comprised Guillain-Barré syndrome patients in group A and healthy controls in group B. The patients were diagnosed on the basis of clinical presentation, nerve conduction study, electromyography, cerebrospinal fluid analysis and biochemical profile. Data was analysed using SPSS 23. RESULTS: Of the 167 subjects, 90(54%) were in group A and 77(46%) were in group B. The mean age of group A was 40.20±14.90 years, while there were 61(67.7%) males and 29(32.2%) females compared to 50 (64.93%) males and 27 (35.06%) females with mean age 38.40±12.34 years in group B. Acute inflammatory demyelinating polyneuropathy was the most common electrophysiological variant of Guillain-Barré syndrome 41(46%). There was significant difference in mean interleukin-17 levels between group A 23.12±3.41 pg/ml and group B 8.82±2.49 (p<0.05). Gastrointestinal infection was the most common preceding infection 51(56.66%). The mean cerebrospinal fluid protein was 100.83±51.32g/dl and albumiocytologic dissociation was found in all the four variants (p= 0.005). CONCLUSION: Guillain-Barré syndrome affected patients regardless of age, while males were more affected than females. Majority of the patients had an antecedent infection before disease onset. Increased levels of interleukin-17 showed involvement of autoimmunity. Albuminocytologic dissociation differentiated it from poliomyelitis.

Clinical and biochemical profile of Guillain Barré Syndrome in Pakistan.

OBJECTIVES: To evaluate the frequency, distribution, characteristics, and biochemical features of Guillain-Barré syndrome (GBS). METHODS: This retrospective, case-control study was conducted between November 2018 and February 2020 at PIMS Hospital, Islamabad and Arid Agriculture University, Rawalpindi, Pakistan. The GBS patients were diagnosed through physical examination and nerve conduction study (NCS). Hemoglobin (Hb), liver function tests, renal function tests, lipid profile, and interleukin-17 levels were investigated through the blood. Cerebrospinal fluid analysis (CSF) was performed to measure albuminocytological dissociation (ACD). RESULTS: One hundred and ten patients and 130 controls participated in this study. Fifty-nine (53.63 %) patients had gastrointestinal infection while 52 (47.27 %) had AIDP. IL-17 serum levels were associated with GBS showing elevated values 18.49 pg/mL (SD=4.10) in cases as compared to controls 10.66 (SD=2.09), p<0.001. The CSF proteins were statistically significant in GBS patients, 127.93 mg/dL (SD=51.28), p=0.002. Ascending weakness was observed in 82 (74.5 %) patients. The results showed that the mean age value of GBS patients was 41.27 years, showing males preponderance. The mean Hb value for males and females was 14.83 g/dL and 10.88 g/dL respectively. Seasonal trends in the disease showed that 40% of GBS patients had infections during the spring. CONCLUSION: The results of the current study suggest that IL-17 levels trigger autoimmunity in GBS patients. The ACD could be used as a diagnostic marker of GBS along with NCS. Antecedent infections were common in a majority of GBS patients.

Direct-acting antiviral agents in the treatment of chronic hepatitis C-Real-life experience from clinical practices in Pakistan.

This study aims to evaluate the clinical effectiveness in terms of sustained virological response (SVR), predictors of SVR and safety of available second-generation generic direct-acting antivirals in Pakistani chronic hepatitis C patients. This is a retrospective study conducted in multiple centers of Pakistan from January 2015 to January 2019. The samples include patients infected with chronic hepatitis C virus, regardless of virus genotype, cirrhosis, or prior treatment. A total of 993 patients were included in the present study, with the majority receiving sofosbuvir with daclatasvir (95%), sofosbuvir with daclatasvir and ribavirin (4%), and sofosbuvir with ribavirin (1%). There were 96% cases of chronic hepatitis, 3% cases compensated cirrhosis, and 1% cases of decompensated cirrhosis. Genotype 3 (99.6%) was the most common genotype. Overall SVR after 12 weeks was 98% for all treatment regimens. High SVR12 was observed with sofosbuvir in combination with daclatasvir (98.5%), then sofosbuvir in combination with daclatasvir and ribavirin (90.2%) and sofosbuvir in combination with ribavirin (75%). SVR rates were high in chronic hepatitis C patients (98.2%) as compared with cirrhotic patients (92.1%) and it was high in treatment-naive (98.8%) then interferon experienced patients (90.1%). In multivariate binary logistic regression analysis, patients' education status, treatment strategy, viral load, and alanine aminotransferase had a statistically significant association with SVR at 12 weeks. No major adverse events occurred which required treatment discontinuation. Generic oral direct acting antiviralss (sofosbuvir with daclatasvir) achieved higher SVR12 rates and were well tolerated in this large real-world cohort of genotype 3 infected patients.

REPORT - CYP2D6*4 null allele frequency in sixteen Pakistani ethnic groups.

CYP2D6 belongs to a family of Cytochrome P450 and is involved in metabolism of a number of commonly prescribed drugs. This study was designed to identify *4 allelic frequency of CYP2D6 in Pakistani population. The ethno-geographic variations in the CYP2D6 alleles are responsible for varied expression of this enzyme and thus influence the metabolic rate and efficacy of prescribed drugs. In total, 976 volunteers belonging to 16 different ethnic groups of Pakistan were screened for CYP2D6*4 polymorphism. The *4 allele was detected in all the ethnic groups with varied frequency ranging from 3.73%-13.64% and an overall average of 7.22% in different ethnic groups of the population. Maximum frequency was detected in northern Pakistani population including Meo (13.64%), Punjabi (11.96%) and Pathan (10.42%). Low frequency (<4%) of*4 polymorphism was observed in Kalash and Makrani groups, whereas an intermediate frequency (5-9%) was observed in all the other ethnic groups. The data indicates that despite ethnic diversity poor metabolizers in Pakistani population are expected to carry CYP2D6*4 allele at a relatively higher frequency than most other Asian populations. (Word count = 186).

Correction to: Novel mutation G324C in WNT1 mapped in a large Pakistani family with severe recessively inherited Osteogenesis Imperfecta.

In the original publication of this article [1], there are two errors in the article which the cDNA position of the pathogenic variant WNT1 p.Gly324Cys should be c.970G>T instead of c.1168G>T.

Novel mutation G324C in WNT1 mapped in a large Pakistani family with severe recessively inherited Osteogenesis Imperfecta.

INTRODUCTION: Osteogenesis imperfecta (OI) is a clinically and genetically heterogeneous disease with skeletal fragility and variable extra-skeletal manifestations. To date several point mutations in 18 different genes causing different types of OI have been identified. Mutations in WNT1 compromise activity of the osteoblasts leading to disturbed bone mass accrual, fragility fractures and progressive skeletal abnormalities. The present study was conducted to determine the underlying genetic cause of an autosomal recessive skeletal dysplasia in a large consanguineous family from Chinute, Pakistan. MATERIALS AND METHODS: Blood was collected from 24 individuals of affected family along with clinical data. Homozygosity mapping was performed to confirm consanguinity. SNPs were identified, followed by whole exome and Sanger sequencing. In silico characterization of WNT1 mutation was performed using multiple platforms. RESULTS: Nine affected family members exhibited severe bone deformities, recurrent fractures, short stature and low bone mineral density. SNP array data revealed homozygous segments > 1 Mb in length accounting for 2.1-12.7% of the genome in affected individuals and their siblings and a single 6,344,821 bp homozygous region in all affected individuals on chromosome 12q12-q13. This region includes two potential OI candidate genes WNT1 and VDR. We did whole-exome sequencing for both genes in two patients and identified a novel damaging missense mutation in exon 4 of WNT1: c.1168G > T (NM_005430) resulting in p.G324C. Sanger sequencing confirmed segregation of mutation with the disease in family. CONCLUSION: We report a novel mutation responsible for OI and our investigation expands the spectrum of disease-causing WNT1 mutations and the resulting OI phenotypes.

Genetics of vascular dementia - review from the ICVD working group.

BACKGROUND: Vascular dementia is a common disorder resulting in considerable morbidity and mortality. Determining the extent to which genes play a role in disease susceptibility and their pathophysiological mechanisms could improve our understanding of vascular dementia, leading to a potential translation of this knowledge to clinical practice. DISCUSSION: In this review, we discuss what is currently known about the genetics of vascular dementia. The identification of causal genes remains limited to monogenic forms of the disease, with findings for sporadic vascular dementia being less robust. However, progress in genetic research on associated phenotypes, such as cerebral small vessel disease, Alzheimer's disease, and stroke, have the potential to inform on the genetics of vascular dementia. We conclude by providing an overview of future developments in the field and how such work could impact patients and clinicians. CONCLUSION: The genetic background of vascular dementia is well established for monogenic disorders, but remains relatively obscure for the sporadic form. More work is needed for providing robust findings that might eventually lead to clinical translation.

Identification and in silico characterization of a novel p.P208PfsX1 mutation in V-ATPase a3 subunit associated with autosomal recessive osteopetrosis in a Pakistani family.

BACKGROUND: Osteopetrosis is a rare inherited bone disorder mainly described as an increased bone density caused by defective osteoclastic bone resorption. To date, genetic variants of eleven genes have been reported so far to be associated with different types of osteopetrosis. However, malignant infantile osteopetrosis, a lethal form of the disease, is mostly (50%) caused by mutation(s) in TCIRG1 gene. In this study, we investigated a consanguineous Pakistani family clinically and genetically to elucidate underlying molecular basis of the infantile osteopetrosis. METHODS: DNA samples from five family members were subjected to SNP-array based whole genome homozygosity mapping. Data was analyzed and potentially pathogenic mutation was identified by Sanger sequencing of two affected as well as three phenotypically healthy individuals in the family. The significance of identified pathogenic variation and its impact on protein structure and function was studied using various bioinformatics tools. RESULTS: DNA samples from five family members were subjected to genome-wide SNP array genotyping and homozygosity mapping which identified ~4 Mb region on chr11 harboring the TCIRG1 gene. Sanger sequencing unveiled a novel homozygous deletion c. 624delC in exon 6 of the TCIRG1 gene encodes a3 subunit of V-ATPase complex. The identified deletion resulted in a frame shift producing a truncated protein of 208 aa. In silico analysis of premature termination of the a3 subunit of V-ATPase complex revealed deleterious effects on the protein structure, predicting impaired or complete loss of V-ATPase function causing infantile osteopetrosis. CONCLUSIONS: Since a3 subunit of V-ATPase complex plays a crucial role in bone resorption process, structurally abnormal a3 subunit might have adversely affected bone resorption process, leading to infantile osteopetrosis in Pakistani family. Therefore, the present study not only expands the genotypic spectrum of osteopetrosis but also improve understandings of the role of V-ATPase a3 subunit in bone resorption process. Moreover, our findings should help in genetic counseling and provide further insight into the disease pathogenesis and potential targeted therapy.

A canonical splice site mutation in GIPC3 causes sensorineural hearing loss in a large Pakistani family.

With homozygosity mapping we have identified two large homozygous regions on chromosome 3q13.11-q13.31 and chromosome 19p13.3-q31.32 in a large Pakistani family suffering from autosomal recessive nonsyndromic hearing impairment (arNSHI). The region on chromosome 19 overlaps with the previously described deafness loci DFNB15, DFNB72 and DFNB95. Mutations in GIPC3 have been shown to underlie the nonsyndromic hearing impairment linked to these loci. Sequence analysis of all exons and exon-intron boundaries of GIPC3 revealed a homozygous canonical splice site mutation, c.226-1G>T, in GIPC3. This is the first mutation described in GIPC3 that affects splicing. The c.226-1G>T mutation is located in the acceptor splice site of intron 1 and is predicted to affect the normal splicing of exon 2. With a minigene assay it was shown to result in the use of an alternative acceptor site in exon 2, resulting in a frameshift and a premature stop codon. This study expands the mutational spectrum of GIPC3 in arNSHI.

Novel mutation in AAA domain of BCS1L causing Bjornstad syndrome.

Bjørnstad syndrome is an extremely rare condition characterized by pilitorti and nerve deafness. Only few large families have been reported worldwide. Here we describe a large Pakistani family with five affected individuals. The hair fibers of all the patients were twisted around their axis and devoid of any pigment. In addition the patients had a moderate-to-severe degree of hearing impairment. Genotyping with high-density single-nucleotide polymorphism arrays showed homozygosity in two intervals on chromosome 2. Linkage with one of these regions (genomic position 218745685-221025443, hg19) was confirmed. This region encompasses the BCS1L gene. Mutations in this gene have previously been associated with Bjørnstad's syndrome. We sequenced the BCS1L gene for identification of the causative mutation in the family. A novel homozygous missense mutation c.901T>A was identified, which segregated with the disease in the family. This mutation results in the amino acid change p.Tyr301Asn and was predicted to be pathogenic by bioinformatics tools.

Study of PKRBD in HCV genotype 3a infected patients in response to interferon therapy in Pakistani population.

BACKGROUND: Hepatitis C virus (HCV) is a major cause of liver cirrhosis and hepatocellular carcinoma and infects about 3% world population. Response to interferon therapy depends upon the genotype of the virus and factors associated with the host. Despite a good response to interferon therapy, a considerable number of genotype 3a infected patients remains unalleviated. RESULTS: In total forty-nine patients including twenty-five non-responders (non-SVR) and twenty-four responders (SVR) were recruited. Patients were tested for viral status at different intervals and the isolated RNA was sequenced for the NS5A region in both groups. The comparison of PKRBD of HCV between the SVR and non-SVR patients did not confirm any significant difference in the number of mutations. However, when the sequence downstream to the PKRBD of NS5A was compared, two important statistically significant mutations were observed; at positions 2309 (Ala to Ser) and 2326 (Gly to Ala). These mutations were then analysed for tertiary protein structure and important structural changes were observed. Statistically significant difference was also observed when age groups of patients were compared; younger patients showed better response than the older ones. CONCLUSIONS: The region between PKRBD and IRRDR may be important for prediction of response to IFN therapy for genotype 3a. ISDR and PKRBD have not shown any involvement in treatment response. Further functional analyses of these findings can help in understanding the involvement of the NS5A region in interferon treatment of HCV-3a infected patients.

Structural and functional insights into a novel homozygous missense pathogenic variant in CUL7 identified in consanguineous Pakistani family.

3M syndrome is a rare genetic familial disorder characterized by short stature, growth retardation, facial dysmorphism, skeletal abnormalities, fleshy protruding heels, and normal intelligence, caused by mutations in the CUL7, OBSL1 and CCDC8 genes. In the present study, a novel homozygous missense variant of CUL7 (NP_001161842.1, c.4493T > C, p.L1498P) has been identified in a consanguineous Pakistani family by whole exome sequencing. In silico structural evaluation, molecular docking and simulation studies of mutant CUL7 provides substantial evidence about its crucial role in the progression of discussed ailment. The newly discovered variant significantly altered the protein's three dimensional structure, leading to abnormal interaction with binding proteins. This computational and experimental investigation provides useful information to drug developers for the synthesis of novel therapeutics against the discussed ailment.Communicated by Ramaswamy H. Sarma.

Case report: Expanding the phenotype of ARHGEF17 mutations from increased intracranial aneurysm risk to a neurodevelopmental disease.

RhoGTPase regulators play a key role in the development of the nervous system, and their dysfunction can result in brain malformation and associated disorders. Several guanine nucleotide exchange factors (GEF) have been linked to neurodevelopmental disorders. In line with this, ARHGEF17 has been recently linked as a risk gene to intracranial aneurysms. Here we report siblings of a consanguineous Pakistani family with biallelic variants in the ARHGEF17 gene associated with a neurodevelopmental disorder with intellectual disability, speech delay and motor dysfunction but not aneurysms. Cranial MRI performed in one patient revealed generalized brain atrophy with an enlarged ventricular system, thin corpus callosum and microcephaly. Whole exome sequencing followed by Sanger sequencing in two of the affected individuals revealed a homozygous missense variant (g.11:73021307, c.1624C>T (NM_014786.4), p.R542W) in the ARHGEF17 gene. This variant is in a highly conserved DCLK1 phosphorylation consensus site (I/L/V/F/M]RRXX[pS/pT][I/L/M/V/F) of the protein. Our report expands the phenotypic spectrum of ARHGEF17 variants from increased intracranial aneurysm risk to neurodevelopmental disease and thereby add ARHGEF17 to the list of GEF genes involved in neurodevelopmental disorders.

Genome sequencing and analysis of genomic diversity in the locally transmitted SARS-CoV-2 in Pakistan.

Surveillance of genetic diversity of the SARS-CoV-2 is extremely important to detect the emergence of more infectious and deadly strains of the virus. In this study, we evaluated mutational events in the SARS-CoV-2 genomes through whole genome sequencing. The samples were collected from COVID-19 patients in different major cities of Pakistan during the four waves of the pandemic (May 2020 to July 2021) and subjected to whole genome sequencing. Using in silico and machine learning tools, the viral mutational events were analyzed, and variants of concern and of interest were identified during each of the four waves. The overall mutation frequency (mutations per genome) increased during the course of the pandemic from 12.19 to 23.63, 31.03, and 41.22 in the first, second, third, and fourth waves, respectively. We determined that the viral strains rose to higher frequencies in local transmission. The first wave had three most common strains B.1.36, B.1.160, and B.1.255, the second wave comprised B.1.36 and B.1.247 strains, the third wave had B.1.1.7 (Alpha variant) and B.1.36 strains, and the fourth waves comprised B.1.617.2 (Delta). Intriguingly, the B.1.36 variants were found in all the waves of the infection indicating their survival fitness. Through phylogenetic analysis, the probable routes of transmission of various strains in the country were determined. Collectively, our study provided an insight into the evolution of SARS-CoV-2 lineages in the spatiotemporal local transmission during different waves of the pandemic, which aided the state institutions in implementing adequate preventive measures.

Occult HCV or delayed viral clearance from lymphocytes of Chronic HCV genotype 3 patients after interferon therapy.

BACKGROUND: A recently discovered occult HCV entity reported by various investigators seems to be highly controversial. Especially, the clinical significance of these findings remains uncertain. For optimal outcome of antiviral therapy, investigation of occult HCV needs a broad-based probe in order to investigate the results of viral therapy and its host/viral interaction. The current study was aimed at determining the prevalence of occult HCV in peripheral blood lymphocytes of predominantly genotype 3 HCV-infected patients after completion of antiviral therapy and to investigate long term outcomes in the presence or absence of PBMC positivity. METHOD: A total of 151 chronic, antiHCV and serum RNA-positive patients were enrolled in the study. Patients with a complete virological response at the end of treatment were screened for the presence of viral RNA in their PBMCs and were followed for up to one year for the presence of serum and PBMC viral genomic RNA. RESULTS: Out of 151 patients, 104 (70%) responded to the prescribed interferon treatment and showed viral-clearance from serum. These were screened for the presence of genomic RNA in their PBMCs. Sixteen samples were PBMC-positive for viral RNA at the end of treatment (EOT). All these patients had also cleared the virus from peripheral blood cells after the 6-12 month follow-up study. CONCLUSION: True occult hepatitis C virus does not exist in our cohort. Residual viremia at the EOT stage merely reflects a difference in viral kinetics in various compartments that remains a target of immune response even after the end of antiviral therapy and is eventually cleared out at the sustained viral response (SVR).

Patient HLA-DRB1* and -DQB1* allele and haplotype association with hepatitis C virus persistence and clearance.

Hepatitis C virus (HCV) infection is prevalent throughout the world and interferon (IFN)-based treatments are currently the only therapeutic option. However, depending upon variations in their human leukocyte antigen (HLA), some patients do not respond well to IFN therapy. The current study evaluated the HLA allele and haplotype distribution of 204 HCV-seropositive individuals from Islamabad, Pakistan, who were receiving standard IFN therapy. In this cohort, 150 patients (74%) showed a sustained virological response to IFN therapy, whereas 54 (26%) did not. In addition to the HCV patients, 102 unrelated healthy volunteers were used as controls. DNA was isolated from the blood of the patients and controls for HLA-DRB1 and HLA-DQB1 allele typing, whilst plasma was used for HCV detection and genotyping. HLA-DRB1*04 was found to impart a significant protective advantage [Bonferroni-corrected P value (pc)=0.047] against HCV infection. In patients on IFN therapy, HLA-DRB1*11 and -DQB1*0301 (pc=0.044) were found to be associated with viral clearance. In contrast, HLA-DRB1*07 (pc=0.008) individually or in combination with HLA-DQB1*02 was found to be associated with viral persistence. These associations of HLA with HCV persistence or clearance will be beneficial in deciding the therapeutic regimen for Pakistani patients infected with HCV genotype 3a.

Biallelic variants in four genes underlying recessive osteogenesis imperfecta.

Osteogenesis imperfecta (OI) is an inherited heterogeneous rare skeletal disorder characterized by increased bone fragility and low bone mass. The disorder mostly segregates in an autosomal dominant manner. However, several rare autosomal recessive and X-linked forms, caused by mutations in 18 different genes, have also been described in the literature. Here, we present five consanguineous families segregating OI in an autosomal recessive pattern. Affected individuals in the five families presented severe forms of skeletal deformities. It included frequent bone fractures with abnormal healing, short stature, facial dysmorphism, osteopenia, joint laxity, and severe scoliosis. In order to search for the causative variants, DNA of at least one affected individual in three families (A-C) were subjected to whole exome sequencing (WES). In two other families (D-E), linkage analysis using highly polymorphic microsatellite markers was followed by Sanger sequencing. Sequence analysis revealed two novels and three previously reported disease-causing variants. The two novel homozygous variants including [c.824G > A; p.(Cys275Tyr)] in the SP7 gene and [c.397C > T, p.(Gln133*)] in the SERPINF1 gene were identified in families A and B, respectively. The three previously reported homozygous variants including [c.497G > A; p.(Arg166His)] in the SPARC gene, (c.359-3C > G; intron 2) and [c.677C > T; p.(Ser226Leu)] in the WNT1 gene were identified in family C, D, and E. In conclusion, our findings provided additional evidence of involvement of homozygous sequence variants in the SP7, SERPINF1, SPARC and WNT1 genes causing severe OI. It also highlights the importance of extensive genetic investigations to search for the culprit gene in each case of skeletal deformity.