Identification of novel homozygous variants in FOXE3 and AP4M1 underlying congenital syndromic anophthalmia and microphthalmia.

BACKGROUND: Anophthalmia and microphthalmia are severe developmental ocular disorders that affect the size of the ocular globe and can be unilateral or bilateral. The disease is found in syndromic as well as non-syndromic forms. It is genetically caused by chromosomal aberrations, copy number variations and single gene mutations, along with non-genetic factors such as viral infections, deficiency of vitamin A and an exposure to alcohol or drugs during pregnancy. To date, more than 30 genes having different modes of inheritance patterns are identified as causing anophthalmia and microphthalmia. METHODS: In the present study, a clinical and genetic analysis was performed of six patients with anophthalmia and microphthalmia and/or additional phenotypes of intellectual disability, developmental delay and cerebral palsy from a large consanguineous Pakistani family. Whole exome sequencing followed by data analysis for variants prioritization and validation through Sanger sequencing was performed to identify the disease causing variant(s). American College of Medical Genetics and Genomics (ACMG) guidelines were applied to classify clinical interpretation of the prioritized variants. RESULTS: Clinical investigations revealed that the affected individuals are afflicted with anophthalmia. Three of the patients showed additional phenotype of intellectual disability, developmental delays and other neurological symptoms. Whole exome sequencing of the DNA samples of the affected members in the family identified a novel homozygous stop gain mutation (NM_012186: c.106G>T: p.Glu36*) in Forkhead Box E3 (FOXE3) gene shared by all affected individuals. Moreover, patients segregating additional phenotypes of spastic paraplegia, intellectual disability, hearing loss and microcephaly showed an additional homozygous sequence variant (NM_004722: c.953G>A: p.Arg318Gln) in AP4M1. Sanger sequencing validated the correct segregation of the identified variants in the affected family. ACMG guidelines predicted the variants to be pathogenic. CONCLUSIONS: We have investigated first case of syndromic anophthalmia caused by variants in the FOXE3 and AP4M1. The present findings are helpful for understanding pathological role of the mutations of the genes in syndromic forms of anophthalmia. Furthermore, the study signifies searching for the identification of second variant in families with patients exhibiting variable phenotypes. In addition, the findings will help clinical geneticists, genetic counselors and the affected family with respect to prenatal testing, family planning and genetic counseling.

Mutational spectrum of CFTR in cystic fibrosis patients with gastrointestinal and hepatobiliary manifestations.

BACKGROUND: Cystic fibrosis (CF) is a rare and debilitating autosomal recessive disorder. It hampers the normal function of various organs and causes severe damage to the lungs, and digestive system leading to recurring pneumonia. Cf also affects reproductive health eventually may cause infertility. The disease manifests due to genetic aberrations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. This study aimed to screen for CFTR gene variants in Pakistani CF patients representing variable phenotypes. METHODS: Clinical exome and Sanger sequencing were performed after clinical characterization of 25 suspected cases of CF (CF1-CF25). ACMG guidelines were followed to interpret the clinical significance of the identified variants. RESULTS: Clinical investigations revealed common phenotypes such as pancreatic insufficiency, chest infections, chronic liver and lung diseases. Some patients also displayed symptoms like gastroesophageal reflux disease (GERD), neonatal cholestasis, acrodermatitis, diabetes mellitus, and abnormal malabsorptive stools. Genetic analysis of the 25 CF patients identified deleterious variants in the CFTR gene. Notably, 12% of patients showed compound heterozygous variants, while 88% had homozygous variants. The most prevalent variant was p. (Met1Thr or Met1?) at 24%, previously not reported in the Pakistani population. The second most common variant was p. (Phe508del) at 16%. Other variants, including p. (Leu218*), p. (Tyr569Asp), p. (Glu585Ter), and p. (Arg1162*) were also identified in the present study. Genetic analysis of one of the present patients showed a pathogenic variant in G6PD in addition to CFTR. CONCLUSION: The study reports novel and reported variants in the CFTR gene in CF patients in Pakistani population having distinct phenotypes. It also emphasizes screening suspected Pakistani CF patients for the p. (Met1Thr) variant because of its increased observance and prevalence in the study. Moreover, the findings also signify searching for additional pathogenic variants in the genome of CF patients, which may modify the phenotypes. The findings contribute valuable information for the diagnosis, genetic counseling, and potential therapeutic strategies for CF patients in Pakistan.

Comprehensive in-silico analysis of deleterious SNPs in APOC2 and APOA5 and their differential expression in cancer and cardiovascular diseases conditions.

Genetic variations in APOC2 and APOA5 genes involve activating lipoprotein lipase (LPL), responsible for the hydrolysis of triglycerides (TG) in blood and whose impaired functions affect the TG metabolism and are associated with metabolic diseases. In this study, we investigate the biological significance of genetic variations at the DNA sequence and structural level using various computational tools. Subsequently, 8 (APOC2) and 17 (APOA5) non-synonymous SNPs (nsSNPs) were identified as high-confidence deleterious SNPs based on the effects of the mutations on protein conservation, stability, and solvent accessibility. Furthermore, based on our docking results, the interaction of native and mutant forms of the corresponding proteins with LPL depicts differences in root mean square deviation (RMSD), and binding affinities suggest that these mutations may affect their function. Furthermore, in vivo, and in vitro studies have shown that differential expression of these genes in disease conditions due to the influence of nsSNPs abundance may be associated with promoting the development of cancer and cardiovascular diseases. Preliminary screening using computational methods can be a helpful start in understanding the effects of mutations in APOC2 and APOA5 on lipid metabolism; however, further wet-lab experiments would further strengthen the conclusions drawn from the computational study.

A novel biallelic variant in the Popeye domain-containing protein 1 (POPDC1) underlies limb girdle muscle dystrophy type 25.

POPDC1 also known as BVES, is a highly conserved transmembrane protein, important for striated muscle function and homeostasis. Pathogenic variants in the POPDC1 gene are associated with limb-girdle muscular dystrophy type 25 (LGMDR25). In the present study, we performed trio-whole exome sequencing (WES) followed by Sanger sequencing on a single family having LGMD clinical features. Protein modeling of all POPDC1 missense variants (POPDC1Pro134Leu , POPDC1Ile193Ser , and POPDC1Ser201Phe ) associated with LGMDR25 were performed using Molecular Dynamics (MD) simulation. We identified a homozygous missense variant (c.401C>T; p.Pro134Leu) in the POPDC1 gene. Altered 3D structure, disruptive fluctuation, less compactness, and instability were observed in all the three variants of POPDC1 protein models. In comparison, POPDC1Ser201Phe protein dynamics were more unstable than other variants. Functional study of newly identified variant would add key answers to underlying mechanisms of the disease.

Molecular Dynamic Simulation Analysis of a Novel Missense Variant in CYB5R3 Gene in Patients with Methemoglobinemia.

Background and Objective: Mutations in the CYB5R3 gene cause reduced NADH-dependent cytochrome b5 reductase enzyme function and consequently lead to recessive congenital methemoglobinemia (RCM). RCM exists as RCM type I (RCM1) and RCM type II (RCM2). RCM1 leads to higher methemoglobin levels causing only cyanosis, while in RCM2, neurological complications are also present along with cyanosis. Materials and Methods: In the current study, a consanguineous Pakistani family with three individuals showing clinical manifestations of cyanosis, chest pain radiating to the left arm, dyspnea, orthopnea, and hemoptysis was studied. Following clinical assessment, a search for the causative gene was performed using whole exome sequencing (WES) and Sanger sequencing. Various variant effect prediction tools and ACMG criteria were applied to interpret the pathogenicity of the prioritized variants. Molecular dynamic simulation studies of wild and mutant systems were performed to determine the stability of the mutant CYB5R3 protein. Results: Data analysis of WES revealed a novel homozygous missense variant NM_001171660.2: c.670A > T: NP_001165131.1: p.(Ile224Phe) in exon 8 of the CYB5R3 gene located on chromosome 22q13.2. Sanger sequencing validated the segregation of the identified variant with the disease phenotype within the family. Bioinformatics prediction tools and ACMG guidelines predicted the identified variant p.(Ile224Phe) as disease-causing and likely pathogenic, respectively. Molecular dynamics study revealed that the variant p.(Ile224Phe) in the CYB5R3 resides in the NADH domain of the protein, the aberrant function of which is detrimental. Conclusions: The present study expanded the variant spectrum of the CYB5R3 gene. This will facilitate genetic counselling of the same and other similar families carrying mutations in the CYB5R3 gene.

Biallelic loss-of-function variants in NEMF cause central nervous system impairment and axonal polyneuropathy.

We aimed to detect the causative gene in five unrelated families with recessive inheritance pattern neurological disorders involving the central nervous system, and the potential function of the NEMF gene in the central nervous system. Exome sequencing (ES) was applied to all families and linkage analysis was performed on family 1. A minigene assay was used to validate the splicing effect of the relevant discovered variants. Immunofluorescence (IF) experiment was performed to investigate the role of the causative gene in neuron development. The large consanguineous family confirms the phenotype-causative relationship with homozygous frameshift variant (NM_004713.6:c.2618del) as revealed by ES. Linkage analysis of the family showed a significant single-point LOD of 4.5 locus. Through collaboration in GeneMatcher, four additional unrelated families' likely pathogenic NEMF variants for a spectrum of central neurological disorders, two homozygous splice-site variants (NM_004713.6:c.574+1G>T and NM_004713.6:c.807-2A>C) and a homozygous frameshift variant (NM_004713.6: c.1234_1235insC) were subsequently identified and segregated with all affected individuals. We further revealed that knockdown (KD) of Nemf leads to impairment of axonal outgrowth and synapse development in cultured mouse primary cortical neurons. Our study demonstrates that disease-causing biallelic NEMF variants result in central nervous system impairment and other variable features. NEMF is an important player in mammalian neuron development.

Recalling the pathology of Parkinson's disease; lacking exact figure of prevalence and genetic evidence in Asia with an alarming outcome: A time to step-up.

Parkinson's disease (PD) is the second most common and progressive neurodegenerative disease globally, with major symptoms like bradykinesia, impaired posture, and tremor. Several genetic and environmental factors have been identified but elucidating the main factors have been challenging due to the disease's complex nature. Diagnosis, prognosis, and management of such diseases are challenging and require effective targeted attention in developing countries. Recently, PD is growing rapidly in many crowded Asian countries as an alarming threat with inadequate knowledge of its prevalence, genetic architecture, and geographic distribution. This study gave an in-depth overview of the prevalence, incidence and genomic/genetics studies published so far in the Asian population. To the best of our knowledge, PD has increased significantly in several Asian countries, including China, South Korea, Japan, Thailand, and Israel over the past few years, requiring a greater level of care and attention. Genetic screening of families with PD at national levels and establishing an official database of PD cases are essential to get a comprehensive and conclusive view of the exact prevalence and genetic diversity of PD in the Asian population to properly manage and treat the disease.

The deleterious variants of N-acetylgalactosamine-6-sulfatase (GalN6S) enzyme trigger Morquio a syndrome by disrupting protein foldings.

Lysosomal enzymes degrade cellular macromolecules, while their inactivation causes human hereditary metabolic disorders. Mucopolysaccharidosis IVA (MPS IVA; Moquio A syndrome) is one of the lysosomal storage disorders caused by a defective Galactosamine-6-sulfatase (GalN6S) enzyme. In several populations, disease incidence is elevated due to missense mutations brought on by non-synonymous allelic variation in the GalN6S enzyme. Here, we studied the effect of non-synonymous single nucleotide polymorphism (nsSNPs) on the structural dynamics of the GalN6S enzyme and its binding with N-acetylgalactosamine (GalNAc) using all-atom molecular dynamics simulation and an essential dynamics approach. Consequently, in this study, we have identified three functionally disruptive mutations in domain-I and domain-II, that is, S80L, R90W, and S162F, which presumably contribute to post-translational modifications. The study delineated that both domains work cooperatively, and alteration in domain II (S80L, R90W) leads to conformational changes in the catalytic site in domain-I, while mutation S162F mainly provokes higher residual flexibility of domain II. These results show that these mutations impair the hydrophobic core, implying that Morquio A syndrome is caused by misfolding of the GalN6S enzyme. The results also show the instability of the GalN6S-GalNAc complex upon substitution. Overall, the structural dynamics resulting from point mutations give the molecular rationale for Moquio A syndrome and, more importantly, the Mucopolysaccharidoses (MPS) family of diseases, re-establishing MPS IVA as a protein-folding disease.Communicated by Ramaswamy H. Sarma.

Structural and dynamics insights into the GBA variants associated with Parkinson's disease.

The GBA1 gene encodes for the lysosomal enzyme glucocerebrosidase (GCase), which maintains glycosphingolipid homeostasis and regulates the autophagy process. Genomic variants of GBA1 are associated with Goucher disease; however, several heterozygous variants of GBA (E326K, T369M, N370S, L444P) are frequent high-risk factors for Parkinson's disease (PD). The underlying mechanism of these variants has been revealed through functional and patient-centered research, but the structural and dynamical aspects of these variants have not yet been thoroughly investigated. In the current study, we used a thorough computational method to pinpoint the structural changes that GBA underwent because of genomic variants and drug binding mechanisms. According to our findings, PD-linked nsSNP variants of GBA showed structural variation and abnormal dynamics when compared to wild-typ. The docking analysis demonstrated that the mutants E326K, N370S, and L444P have higher binding affinities for Ambroxol. Root means square deviation (RMSD), Root mean square fluctuation analysis (RMSF), and MM-GBSA analysis confirmed that the Ambroxol are more stable in the binding site of N370S and L444P, and that their binding affinities are stronger as compared to the wild-type and T369M variants of GBA. The evaluation of hydrogen bonds and the calculation of the free binding energy provided additional evidence in favor of this conclusion. When docked with Ambroxol, GBA demonstrated an increase in binding affinity and catalytic activity. Understanding the therapeutic efficacy and potential against the aforementioned changes in the GBA will be beneficial in order to use more efficient methods for developing novel drugs.Communicated by Ramaswamy H. Sarma.

Truncated DNM1 variant underlines developmental delay and epileptic encephalopathy.

Background: Developmental and epileptic encephalopathies (DEEs) signify a group of heterogeneous neurodevelopmental disorder associated with early-onset seizures accompanied by developmental delay, hypotonia, mild to severe intellectual disability, and developmental regression. Variants in the DNM1 gene have been associated with autosomal dominant DEE type 31A and autosomal recessive DEE type 31B. Methods: In the current study, a consanguineous Pakistani family consisting of a proband (IV-2) was clinically evaluated and genetically analyzed manifesting in severe neurodevelopmental phenotypes. WES followed by Sanger sequencing was performed to identify the disease-causing variant. Furthermore, 3D protein modeling and dynamic simulation of wild-type and mutant proteins along with reverse transcriptase (RT)-based mRNA expression were checked using standard methods. Results: Data analysis of WES revealed a novel homozygous non-sense variant (c.1402G>T; p. Glu468*) in exon 11 of the DNM1 gene that was predicted as pathogenic class I. Variants in the DNM1 gene have been associated with DEE types 31A and B. Different bioinformatics prediction tools and American College of Medical Genetics guidelines were used to verify the identified variant. Sanger sequencing was used to validate the disease-causing variant. Our approach validated the pathogenesis of the variant as a cause of heterogeneous neurodevelopmental disorders. In addition, 3D protein modeling showed that the mutant protein would lose most of the amino acids and might not perform the proper function if the surveillance non-sense-mediated decay mechanism was skipped. Molecular dynamics analysis showed varied trajectories of wild-type and mutant DNM1 proteins in terms of root mean square deviation, root mean square fluctuation and radius of gyration. Similarly, RT-qPCR revealed a substantial reduction of the DNM1 gene in the index patient. Conclusion: Our finding further confirms the association of homozygous, loss-of-function variants in DNM1 associated with DEE type 31B. The study expands the genotypic and phenotypic spectrum of pathogenic DNM1 variants related to DNM1-associated pathogenesis.

Molecular Insights into the Role of Pathogenic nsSNPs in GRIN2B Gene Provoking Neurodevelopmental Disorders.

The GluN2B subunit of N-methyl-D-aspartate receptors plays an important role in the physiology of different neurodevelopmental diseases. Genetic variations in the GluN2B coding gene (GRIN2B) have consistently been linked to West syndrome, intellectual impairment with focal epilepsy, developmental delay, macrocephaly, corticogenesis, brain plasticity, as well as infantile spasms and Lennox-Gastaut syndrome. It is unknown, however, how GRIN2B genetic variation impacts protein function. We determined the cumulative pathogenic impact of GRIN2B variations on healthy participants using a computational approach. We looked at all of the known mutations and calculated the impact of single nucleotide polymorphisms on GRIN2B, which encodes the GluN2B protein. The pathogenic effect, functional impact, conservation analysis, post-translation alterations, their driving residues, and dynamic behaviors of deleterious nsSNPs on protein models were then examined. Four polymorphisms were identified as phylogenetically conserved PTM drivers and were related to structural and functional impact: rs869312669 (p.Thr685Pro), rs387906636 (p.Arg682Cys), rs672601377 (p.Asn615Ile), and rs1131691702 (p.Ser526Pro). The combined impact of protein function is accounted for by the calculated stability, compactness, and total globularity score. GluN2B hydrogen occupancy was positively associated with protein stability, and solvent-accessible surface area was positively related to globularity. Furthermore, there was a link between GluN2B protein folding, movement, and function, indicating that both putative high and low local movements were linked to protein function. Multiple GRIN2B genetic variations are linked to gene expression, phylogenetic conservation, PTMs, and protein instability behavior in neurodevelopmental diseases. These findings suggest the relevance of GRIN2B genetic variations in neurodevelopmental problems.

Clinical and genetic characterization of patients segregating variants in KPTN, MINPP1, NGLY1, AP4B1, and SON underlying neurodevelopmental disorders: Genetic and phenotypic expansion.

Neurodevelopmental disorders (NDDs) are heterogeneous genetic conditions of the central nervous system (CNS). Primary phenotypes of NDDs include epilepsy, loss of developmental skills, abnormal movements, muscle weakness, ocular anomalies, hearing problems, and macro- or microcephaly. NDDs occur due to variants in genes encoding proteins involved in the structure and function of CNS, thus interrupting its normal physiological role. In the study presented here, four consanguineous families (A-D), with members showing neurodevelopmental symptoms, were recruited for clinical and genetic characterization of the phenotypes. Clinical examinations, including Seguin Form Board Test (SFBT), Vineland Social Maturity Scale (VSMS), brain Magnetic Resonance Imaging (MRI), Electroencephalogram (EEG), Electromyography (EMG), Nerve Conduction Velocity (NCV), and Magnetic Resonance Spectroscopy, were employed to characterize the disease phenotypes. Whole exome sequencing (WES) followed by Sanger sequencing was employed to search for the genetic basis of the neurological symptoms observed in four families (A-D). Two of these families (A, B) were of Saudi Arabian origin, and two others (C, D) were of Pakistan origin. Two homozygous missense (KPTN: NM_007059.4:c.301T>G: NP_008990.2; p.(Phe101Val) and MINPP1:NM_001178118.2:c.1202G>A: NP_001171588.1; p.(Arg401Gln)) variants in families A and B, respectively, and two homozygous nonsense (NGLY1:NM_018297.3:c.1534_1541dup: NP_060767.2; p.(Ser515LysfsTer51) and AP4B1:NM_001253852:c.1668G>A: NP_001240781.1; p.(Trp556X)) variants in families C and D, respectively, were identified. Interestingly, additional heterozygous nonsense variant in SON: NM_138927.2: c.5753_5756del: NP_620305.3; p.(Val1918GlufsTer87) and a homozygous variant in FLG (FLG: NM_002016.2:c.7339C>T: NP_002007.1; p.(Arg2447X) were detected in families A and D, respectively. Further, we determined the deleteriousness of each variant through computational approaches. The present study expands the phenotypic and genetic spectrum of NDD-associated genes (KPTN, MINPP1, NGLY1, and AP4B1). Moreover, additional nonsense variants (SON: c.5753_5756del and FLG: c.7339C>T) identified in two families segregating with the phenotype might explain the phenotypic variability and severity in our patients.

A splice-site variant (c.3289-1G>T) in OTOF underlies profound hearing loss in a Pakistani kindred.

BACKGROUND: Hearing loss/deafness is a common otological disorder found in the Pakistani population due to the high prevalence of consanguineous unions, but the full range of genetic causes is still unknown. METHODS: A large consanguineous Pakistani kindred with hearing loss was studied. Whole-exome sequencing and Sanger sequencing were performed to search for the candidate gene underlying the disease phenotype. A minigene assay and reverse transcription polymerase chain reaction was used to assess the effect of splicing variants. RESULTS: The splicing variants of OTOF (NM_194248, c.3289-1G>T) cosegregated with the disease phenotype in this Pakistani family. The substitution of a single base pair causes the deletion of 10 bp (splicing variant 1) or 13 bp (splicing variant 2) from exon 27, which results in truncated proteins of 1141 and 1140 amino acids, respectively. CONCLUSION: Our findings reveal an OTOF splice-site variant as pathogenic for profound hearing loss in this family.

Excess of RALGAPB de novo variants in neurodevelopmental disorders.

Autism spectrum disorder is a neurodevelopmental disorder (NDD) with complex genetic architecture marked primarily by social and communication impairments along with deficits in restrictive and repetitive behaviors. Due to the complex nature and genetic heterogeneity of the disease, genotype and phenotype correlation remains challenging. Prior studies have implicated RALGAPB as a candidate gene for ASD, but stringent analysis is required to determine the pathogenicity. By targeted sequencing, we identified a new de novo RALGAPB missense variant (c.1238C> T; p.T413M) in an ASD family. By leveraging published large-scale genome sequencing studies, we curated five de novo likely gene-disruptive (LGD) variants and 5 de novo missense variants in ASD and related NDDs and revealed a genome-wide significant excess of RALGAPB de novo LGD variants (P_adjust = 0.0053). Quantitative reverse transcription PCR revealed that the frameshift variant c.1927dupA; p.N643fs*3 reduced mRNA expression levels confirming the loss-of-function effect. Co-expression analysis using human brain transcriptome data provide the potential functional link of RALGAPB and 38 ASD and/or NDD genes. Our study suggests RALGAPB as a new NDD risk gene which should be considered in clinical diagnosis of ASD and related NDDs.