Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces. BMC Bioinformatics. 2010;11:548. [PMID: 21059217 PMCID: PMC2992548 DOI: 10.1186/1471-2105-11-548]

Expanding the genetic spectrum of the pyruvate carboxylase deficiency with novel missense, deep intronic and structural variants

Introduction

Pathogenic variants in the pyruvate carboxylase (PC) gene cause a wide spectrum of recessive phenotypes, ranging from the early-onset fatal encephalopathy to the adult-onset benign form.

Results

Patient 1 is a 6 y.o. boy with ataxia, hypoglycemia and episodes of lactic acidosis. WGS revealed the novel heterozygous missense variant c.1372A > G (p.Asn458Asp) in the PC gene. Additional analysis revealed discordant reads mapped to chromosomes 11 and 1, so a reciprocal translocation disrupted the PC gene was suspected. The translocation was validated via FISH-analysis and Sanger sequencing of its boundaries.

Patient 2 is a 13 y.o. girl with psychomotor delay, episodes of lactic acidosis and ketonuria. WES revealed the novel homozygous intronic variant c.1983-116C > T. The PC's mRNA analysis demonstrated the exonization of several intron 16 sequences and some residual amount of WT mRNA isoform.

Two other patients had more severe course of the disease. Their genotype represents missense variants in compound heterozygous and homozygous state (c.1876C > T (p.Arg626Trp), c.2606G > C (p.Gly869Ala), c.2435C > A (p.Ala812Asp).

Conclusion

In patients with metabolic crises, lactic acidosis and hypoglycemia analysis of PC gene is recommended. WGS with deep bioinformatic analysis should be taken into consideration when none or the only one pathogenic variant in the PC gene is found.

•

In patients with metabolic crises with lactic acidosis and hypoglycemia WGS with deep bioinformatic analysis should be taken when none or the only one heterozygous pathogenic variant in the PC gene is found.

•

PC gene may have more deep intronic pathogenic variants than included in known mutation databases.

Targeted Next Generation Sequencing Revealed Novel Variants in the PKD1 and PKD2 Genes of Iranian Patients with Autosomal Dominant Polycystic Kidney Disease

BACKGROUND

Autosomal dominant polycystic kidney disease (ADPKD), one of the common inherited disorders in humans, is characterized by the development and enlargement of renal cysts, often leading to end-stage renal disease (ESRD). In this study, Iranian ADPKD families were subjected to high-throughput DNA sequencing to find potential causative variants facilitating the way toward risk assessment and targeted therapy.

METHODS

Our protocol was based on the targeted next generation sequencing (NGS) panel previously developed in our center comprising 12 genes involved in PKD. This panel has been applied to investigate the genetic causes of 32 patients with a clinical suspicion of ADPKD.

RESULTS

We identified a total of 31 variants for 32 individuals, two of which were each detected in two individuals. Twenty-seven out of 31 detected variants were interpreted as pathogenic/likely pathogenic and the remaining 4 of uncertain significance with a molecular diagnostic success rate of 87.5%. Among these variants, 25 PKD1/2 pathogenic/likely pathogenic variants were detected in 32 index patients (78.1%), and variants of uncertain significance in four individuals (12.5% in PKD1/2). The majority of variants was identified in PKD1 (74.2%). Autosomal recessive PKD was identified in one patient, indicating the similarities between recessive and dominant PKD. In concordance with earlier studies, this biallelic PKD1 variant, p.Arg3277Cys, leads to rapidly progressive and severe disease with very early-onset ADPKD.

CONCLUSION

Our findings suggest that targeted gene panel sequencing is expected to be the method of choice to improve diagnostic and prognostic accuracy in PKD patients with heterogeneity in genetic background.

The TALE never ends: A comprehensive overview of the role of PBX1, a TALE transcription factor, in human developmental defects

PBX1 is a highly conserved atypical homeodomain transcription factor (TF) belonging to the TALE (three amino acid loop extension) family. Dimerized with other TALE proteins, it can interact with numerous partners and reach dozens of regulating sequences, suggesting its role as a pioneer factor. PBX1 is expressed throughout the embryonic stages (as early as the blastula stage) in vertebrates. In human, PBX1 germline variations are linked to syndromic renal anomalies (CAKUTHED). In this review, we summarized available data on PBX1 functions, PBX1-deficient animal models, and PBX1 germline variations in humans. Two types of genetic alterations were identified in PBX1 gene. PBX1 missense variations generate a severe phenotype including lung hypoplasia, cardiac malformations, and sexual development defects (DSDs). Conversely, truncating variants generate milder phenotypes (mainly cryptorchidism and deafness). We suggest that defects in PBX1 interactions with various partners, including proteins from the HOX (HOXA7, HOXA10, etc.), WNT (WNT9B, WNT3), and Polycomb (BMI1, EED) families are responsible for abnormal proliferation and differentiation of the embryonic mesenchyme. These alterations could explain most of the defects observed in humans. However, some phenotype variability (especially DSDs) remains poorly understood. Further studies are needed to explore the TALE family in greater depth.

Novel GRHL3 Variants in a South African Cohort With Cleft Lip and Palate

OBJECTIVE

The etiology of cleft palate (CP) is poorly understood compared with that of cleft lip with or without palate (CL ± P). Recently, variants in Grainyhead like transcription factor 3 (GRHL3) were reported to be associated with a risk for CP in European and some African populations including Nigeria, Ghana, and Ethiopia. In order to identify genetic variants that may further explain the etiology of CP, we sequenced GRHL3 in a South African population to determine if rare variants in GRHL3 are associated with the presence of syndromic or nonsyndromic CP.

DESIGN

We sequenced the exons of GRHL3 in 100 cases and where possible, we sequenced the parents of the individuals to determine the segregation pattern and presence of de novo variants.

SETTING

The cleft clinics from 2 public, tertiary hospitals in Durban, South Africa (SA), namely Inkosi Albert Luthuli Central Hospital and KwaZulu-Natal Children's Hospital.

PATIENTS, PARTICIPANTS

One hundred patients with CL ± P and their parents.

INTERVENTIONS

Saliva samples were collected.

MAIN OUTCOME MEASURES

To ascertain the genetic variants in the GRHL3 gene in patients with CL ± P in SA.

RESULTS

Five variants in GRHL3 were observed; 3 were novel and 2 were known variants. The novel variants were intronic variants (c.1062 + 77A>G and c.627 + 1G>A) and missense variant (p.Asp169Gly).

CONCLUSIONS

This study provides further evidence that variants in GRHL3 contribute to the risk of nonsyndromic CP in African populations, specifically, in the South African population.

RASopathies due to de novo pathogenic variants: clinical features, genetic findings and outcomes in nine neonates born with congenital heart defects

Background

There are limited information available related to neonatal characteristics of RASopathies, a group of autosomal dominant syndromes with considerable phenotypic overlap.

Methods

The retrospective review revealed 9 neonates born with congenital heart defects (CHDs) and diagnosed as RASopathies due to de novo mutations (DNMs) by trio-based exome sequencing (ES) between January 2017 and December 2020. We report in details of the neonatal course, molecular analysis and 180-days of age follow-up in affected individuals.

Results

The early clinical spectrum included various types of CHDs, less noticeable multiple extracardiac anomalies and unspecific symptoms like poor feeding. Of the 8 variants identified from 6 genes, 2 in RASA1 were novel: (NM_002890.2: c.2828 T > C (p.Leu943Pro)) and (NM_002890.2: c.2001del (p.Pro668Leufs*10)), which functionally impaired the protein structure. There was a relatively high mortality rate of 33.33% (3/9) for all the defects combined. A RAF1-deficient male and a RASA1-deficient male survived from severe heart failure by surgical interventions in early life.

Conclusions

Our results revealed that family-based ES was useful in identifying DNMs and causal genes for sporadic diseases and screening Rasopathies shortly after birth. We recommended a family-based ES and a full phenotypic evaluation including echocardiogram, magnetic resonance imaging, ultrasonography and coagulation screening in neonates with CHDs and a suspected genetic etiology.

Computational and structural investigation of Palmitoyl-Protein Thioesterase 1 (PPT1) protein causing Neuronal Ceroid Lipofuscinoses (NCL)

The Neuronal Ceroid Lipofuscinoses (NCL) are a group of progressive neurodegenerative disorders, associated with 14 Ceroid Lipofuscinosis Neuronal genes (CLN1-14). The mutations in the Palmitoyl-Protein Thioesterase 1 (PPT1) protein serve as one of the major reasons for the causative of NCL. The PPT1 involves degrading and modifying cysteine residues in proteins or peptides by removing thioester-linked fatty acyl groups like palmitate prefers acyl chains of 14-18 carbons in length. In this study, we have analyzed the impact of PPT1 mutations on the deleteriousness, stability, conservative nature of amino acid, and impact of mutations on the protein structure. We have also used molecular dynamics simulations using GROMACS to perceive the alteration in the dynamic behavior of the PPT1 at the residual level. In this study, we have retrieved 23 PPT1 mutations from the UniProt database, and these were subjected to a series of analyses using varied computer algorithms. From these analyses, out of 23 mutations, 16 mutations were identified as deleterious. Among 16, eight mutations were identified to destabilize the protein structure, and finally, two mutations (W38C and L222P) were found to be positioned in the highly conserved region. The structural impact study observed that the mutant proline could disrupt the alpha helix formed by the leucine at position 222. Finally, from the molecular dynamics simulations, we observed that due to the mutations (W38C and L222P), the protein had experienced higher deviation, fluctuation, and lower compactness. These structural changes elucidate that these mutations can impact the structure and function of the PPT1 protein.

Computational Analysis of Deleterious SNPs in NRAS to Assess Their Potential Correlation With Carcinogenesis

The NRAS gene is a well-known oncogene that acts as a major player in carcinogenesis. Mutations in the NRAS gene have been linked to multiple types of human tumors. Therefore, the identification of the most deleterious single nucleotide polymorphisms (SNPs) in the NRAS gene is necessary to understand the key factors of tumor pathogenesis and therapy. We aimed to retrieve NRAS missense SNPs and analyze them comprehensively using sequence and structure approaches to determine the most deleterious SNPs that could increase the risk of carcinogenesis. We also adopted structural biology methods and docking tools to investigate the behavior of the filtered SNPs. After retrieving missense SNPs and analyzing them using six in silico tools, 17 mutations were found to be the most deleterious mutations in NRAS. All SNPs except S145L were found to decrease NRAS stability, and all SNPs were found on highly conserved residues and important functional domains, except R164C. In addition, all mutations except G60E and S145L showed a higher binding affinity to GTP, implicating an increase in malignancy tendency. As a consequence, all other 14 mutations were expected to increase the risk of carcinogenesis, with 5 mutations (G13R, G13C, G13V, P34R, and V152F) expected to have the highest risk. Thermodynamic stability was ensured for these SNP models through molecular dynamics simulation based on trajectory analysis. Free binding affinity toward the natural substrate, GTP, was higher for these models as compared to the native NRAS protein. The Gly13 SNP proteins depict a differential conformational state that could favor nucleotide exchange and catalytic potentiality. A further application of experimental methods with all these 14 mutations could reveal new insights into the pathogenesis and management of different types of tumors.

Unravelling the genetic and phenotypic heterogeneity of SPTA1 gene variants in Hereditary Elliptocytosis and Hereditary Pyropoikilocytosis patients using next-generation sequencing

Hereditary Elliptocytosis (HE) and Hereditary Pyropoikilocytosis (HPP) are clinically and genetically heterogeneous red cell membranopathies that result from the defects in the horizontal linkage between RBC (red blood cell) membrane and cytoskeletal proteins affecting its mechanical stability and deformability thereby reducing its lifespan. The principal defect in HE and HPP is due to dysfunction or deficiency of RBC cytoskeletal proteins namely, α-spectrin (SPTA1), β-spectrin (SPTB) and protein 4.1R (EPB41R). This study reports the genetic and phenotypic heterogeneity of 10 Indian patients (5 with HE and 5 with HPP)harboringSPTA1 gene variants. We used targeted next-generation sequencing (t-NGS) to characterize the causative genetic variants in 10 HE/HPP suspected patients and studied the correlation between the identified variants with their corresponding phenotypic features.t-NGS detected 12 SPTA1 variants, out of which 8 are novel. Nearly all of the detected variants have a damaging effect on the protein stability and function, as shown by the insilico analysis. The possible effect of the detected variants on the protein structure was studied using the HOPE software and DynaMut tools wherever possible. To the best of our knowledge, this is the first report on HE/HPP cases confirmed by a genetic study from India. To conclude, HE is caused by monoallelic mutations while HPP, the more severe form, is typically caused by biallelic (homozygous or compound heterozygous) mutations justifying the phenotypic heterogeneity associated with patients. Moreover, analysis at the molecular level by NGS permits diagnosis in these disorders with highly variable heterogeneity requiring regular transfusions and may facilitate prognostic contemplations.

The T350G Variation of Human Papillomavirus 16 E6 Gene Prevails in Oropharyngeal Cancer from a Small Cohort of Greek Patients

Recent trends have shown a dramatic rise in the incidence of oropharyngeal squamous cell carcinoma strongly associated with high-risk human papillomavirus (HPV) of type 16. The genetic variability of HPV16 has been extensively studied in cervical cancer but there are very limited published data concerning the genetic variations of this HPV type in oropharyngeal cancer. In the present study, the genetic variations of HPV16 E6 gene sequences originated from a small cohort of Greek patients diagnosed with oropharyngeal cancer were assessed. The vast majority of the sequences clustered within the European variant branch. The T350G variation was found to be the predominant one. This finding may indicate the need for further studies that could explain the possible impact of this variant in the pathomechanisms of oropharyngeal cancer.

Pathogenic genetic variants identified in Australian families with paediatric cataract

Objective

Paediatric (childhood or congenital) cataract is an opacification of the normally clear lens of the eye and has a genetic basis in at least 18% of cases in Australia. This study aimed to replicate clinical gene screening to identify variants likely to be causative of disease in an Australian patient cohort.

Methods and analysis

Sixty-three reported isolated cataract genes were screened for rare coding variants in 37 Australian families using genome sequencing.

Results

Disease-causing variants were confirmed in eight families with variant classification as ‘likely pathogenic’. This included novel variants PITX3 p.(Ter303LeuextTer100), BFSP1 p.(Glu375GlyfsTer2), and GJA8 p.(Pro189Ser), as well as, previously described variants identified in genes GJA3, GJA8, CRYAA, BFSP1, PITX3, COL4A1 and HSF4. Additionally, eight variants of uncertain significance with evidence towards pathogenicity were identified in genes: GJA3, GJA8, LEMD2, PRX, CRYBB1, BFSP2, and MIP.

Conclusion

These findings expand the genotype–phenotype correlations of both pathogenic and benign variation in cataract-associated genes. They further emphasise the need to develop additional evidence such as functional assays and variant classification criteria specific to paediatric cataract genes to improve interpretation of variants and molecular diagnosis in patients.

De Novo ZMYND8 variants result in an autosomal dominant neurodevelopmental disorder with cardiac malformations

PURPOSE

ZMYND8 encodes a multidomain protein that serves as a central interactive hub for coordinating critical roles in transcription regulation, chromatin remodeling, regulation of super-enhancers, DNA damage response and tumor suppression. We delineate a novel neurocognitive disorder caused by variants in the ZMYND8 gene.

METHODS

An international collaboration, exome sequencing, molecular modeling, yeast two-hybrid assays, analysis of available transcriptomic data and a knockdown Drosophila model were used to characterize the ZMYND8 variants.

RESULTS

ZMYND8 variants were identified in 11 unrelated individuals; 10 occurred de novo and one suspected de novo; 2 were truncating, 9 were missense, of which one was recurrent. The disorder is characterized by intellectual disability with variable cardiovascular, ophthalmologic and minor skeletal anomalies. Missense variants in the PWWP domain of ZMYND8 abolish the interaction with Drebrin and missense variants in the MYND domain disrupt the interaction with GATAD2A. ZMYND8 is broadly expressed across cell types in all brain regions and shows highest expression in the early stages of brain development. Neuronal knockdown of the DrosophilaZMYND8 ortholog results in decreased habituation learning, consistent with a role in cognitive function.

CONCLUSION

We present genomic and functional evidence for disruption of ZMYND8 as a novel etiology of syndromic intellectual disability.

In silico Prediction of Deleterious Single Nucleotide Polymorphism in S100A4 Metastatic Gene: Potential Early Diagnostic Marker

S100A4 protein overexpression has been reported in different types of cancer and plays a key role by interacting with the tumor suppressor protein Tp53. Single nucleotide polymorphisms (SNP) in S100A4 could directly influence the biomolecular interaction with the tumor suppressor protein Tp53 due to their aberrant conformations. Hence, the study was designed to predict the deleterious SNP and its effect on the S100A4 protein structure and function. Twenty-one SNP data sets were screened for nonsynonymous mutations and subsequently subjected to deleterious mutation prediction using different computational tools. The screened deleterious mutations were analyzed for their changes in functionality and their interaction with the tumor suppressor protein Tp53 by protein-protein docking analysis. The structural effects were studied using the 3DMissense mutation tool to estimate the solvation energy and torsion angle of the screened mutations on the predicted structures. In our study, 21 deleterious nonsynonymous mutations were screened, including F72V, E74G, L5P, D25E, N65S, A28V, A8D, S20L, L58P, and K26N were found to be remarkably conserved by exhibiting the interaction either with the EF-hand 1 or EF-hand 2 domain. The solvation and torsion values significantly deviated for the mutant-type structures with S20L, N65S, and F72L mutations and showed a marked reduction in their binding affinity with the Tp53 protein. Hence, these deleterious mutations might serve as prospective targets for diagnosing and developing personalized treatments for cancer and other related diseases.

Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

In silico comprehensive analysis of coding and non-coding SNPs in human mTOR protein

The mammalian/mechanistic target of rapamycin (mTOR) protein is an important growth regulator and has been linked with multiple diseases including cancer and diabetes. Non-synonymous mutations of this gene have already been found in patients with renal clear cell carcinoma, melanoma, and acute lymphoid leukemia among many others. Such mutations can potentially affect a protein’s structure and hence its functions. In this study, therefore, the most deleterious SNPs of mTOR protein have been determined to identify potential biomarkers for various disease treatments. The aim is to generate a structured dataset of the mTOR gene’s SNPs that may prove to be an asset for the identification and treatment of multiple diseases associated with the target gene. Both sequence and structure-based approaches were adopted and a wide variety of bioinformatics tools were applied to analyze the SNPs of mTOR protein. In total 11 nsSNPs have been filtered out of 2178 nsSNPs along with two non-coding variations. All of the nsSNPs were found to destabilize the protein structure and disrupt its function. While R619C, A1513D, and T1977R mutations were shown to alter C alpha distances and bond angles of the mTOR protein, L509Q, R619C and N2043S were predicted to disrupt the mTOR protein’s interaction with NBS1 protein and FKBP1A/rapamycin complex. In addition, one of the non-coding SNPs was shown to alter miRNA binding sites. Characterizing nsSNPs and non-coding SNPs and their harmful effects on a protein’s structure and functions will enable researchers to understand the critical impact of mutations on the molecular mechanisms of various diseases. This will ultimately lead to the identification of potential targets for disease diagnosis and therapeutic interventions.

Computational Insights into the Structural and Functional Impacts of nsSNPs of Bone Morphogenetic Proteins

BMPs (bone morphogenetic proteins) are multipurpose (transforming growth factor)TGF-superfamily released cytokines. These glycoproteins, acting as disulfide-linked homo- or heterodimers, are highly potent regulators of bone and cartilage production and repair, cell proliferation throughout embryonic development, and bone homeostasis in the adults. Due to the fact that genetic variation might influence structural functions, this study is aimed to determine the pathogenic effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in BMP genes. The implications of these variations, investigated using computational analysis and molecular models of the mature TGF-β domain, revealed the impact of modifications on the function of BMP protein. The three-dimensional (3D) structure analysis was performed on the nsSNP Y316S, V386G, E387G, C389G, and C391G nsSNP in the TGF-β domain of chicken BMP2 and H344P, S347P, V357A nsSNP in the TGF-β domain of chicken BMP4 protein that was anticipated to be harmful and of high risk. The ability of the proteins to perform variety of tasks interact with other molecules depends on their tertiary structural composition. The current analysis revealed the four most damaging variants (Y316S, V386G, E387G, C389G, and C391G), highly conserved and functional and are located in the TGF-beta domain of BMP2 and BMP4. The amino acid substitutions E387G, C389G, and C391G are discovered in the binding region. It was observed that the mutations in the TGF-beta domain caused significant changes in its structural organization including the substrate binding sites. Current findings will assist future research focused on the role of these variants in BMP function loss and their role in skeletal disorders, and this will possibly help to develop practical strategies for treating bone-related conditions.

A novel missense variant in CEACAM16 gene causes autosomal dominant nonsyndromic hearing loss

Aim

Autosomal dominant non‐syndromic hearing loss is a common sensorineural disorder with extremely high genetic heterogeneity. CEA antigen‐related cell adhesion molecule 16(CEACAM16)is a secreted glycoprotein encoded by the CEACAM16 gene. Mutations in CEACAM16 lead to autosomal dominant non‐syndromic hearing loss in humans, due defects in the tectorial membrane of the inner ear. Here we reported a novel missense variant in CEACAM16 gene causes autosomal dominant non‐syndromic hearing loss.

Material and methods

A four‐generation Chinese family affected by late‐onset and progressive hearing loss was enrolled in this study. The proband was analyzed by targeted next‐generation sequencing and bioinformatic analysis. And in vitro experiments were performed in overexpressed transfected HEK293T cells to investigate the pathogenesis of the mutant protein.

Results

We identified a novel missense variant in the CEACAM16 gene c.763A>G; (p.Arg255Gly) as causing autosomal dominant non‐syndromic hearing loss in the Chinese family. Using Western blot analysis, ELISA, and immunofluorescence we found increased expression level of the secreted mutant CEACAM16 protein, both intracellularly and extracellularly, compared with wild type CEACAM16 protein.

Conclusion

Our study showed that the p.Arg255Gly variant leads to increased secretion of mutant CEACAM16 protein, with potential deleterious effect to the function of the protein. Our findings expand the mutation spectrum of CEACAM16, and further the understanding CEACAM16 function and implications in disease.

Determination of deleterious single-nucleotide polymorphisms of human LYZ C gene: an in silico study

Background

Single-nucleotide polymorphisms (SNPs) have a crucial function in affecting the susceptibility of individuals to diseases and also determine how an individual responds to different treatment options. The present study aimed to predict and characterize deleterious missense nonsynonymous SNPs (nsSNPs) of lysozyme C (LYZ C) gene using different computational methods. Lyz C is an important antimicrobial peptide capable of damaging the peptidoglycan layer of bacteria leading to osmotic shock and cell death. The nsSNPs were first analyzed by SIFT and PolyPhen v2 tools. The nsSNPs predicted as deleterious were then assessed by other in silico tools — SNAP, PROVEAN, PhD-SNP, and SNPs & GO. These SNPs were further examined by I-Mutant 3.0 and ConSurf. GeneMANIA and STRING tools were used to study the interaction network of the LYZ C gene. NetSurfP 2.0 was used to predict the secondary structure of Lyz C protein. The impact of variations on the structural characteristics of the protein was studied by HOPE analysis. The structures of wild type and variants were predicted by SWISS-MODEL web server, and energy minimization was carried out using XenoPlot software. TM-align tool was used to predict root-mean-square deviation (RMSD) and template modeling (TM) scores.

Results

Eight missense nsSNPs (T88N, I74T, F75I, D67H, W82R, D85H, R80C, and R116S) were found to be potentially deleterious. I-Mutant 3.0 determined that the variants decreased the stability of the protein. ConSurf predicted rs121913547, rs121913549, and rs387906536 nsSNPs to be conserved. Interaction network tools showed that LYZ C protein interacted with lactoferrin (LTF). HOPE tool analyzed differences in physicochemical properties between wild type and variants. TM-align tool predicted the alignment score, and the protein folding was found to be identical. PyMOL was used to visualize the superimposition of variants over wild type.

Conclusion

This study ascertained the deleterious missense nsSNPs of the LYZ C gene and could be used in further experimental analysis. These high-risk nsSNPs could be used as molecular targets for diagnostic and therapeutic interventions.

Identification of novel pathogenic variants in the GCDH gene and assessment of neurodevelopmental outcomes in 24 children with glutaric aciduria type 1

AIM

To evaluate the pathogenic variants in GCDH gene and to assess the neurodevelopmental outcomes in children with Glutaric aciduria type 1 (GA-1).

METHOD

Cross-sectional observational study between January 2019 and June 2020 in consecutive North Indian children with a clinical and biochemical suspicion of GA-1. Variants in the coding regions of GCDH gene were identified through Sanger sequencing. Neurodevelopmental and quality of life assessment was done using standardized scales.

RESULTS

24 children with GA-1 were identified. The median age at diagnosis was 12 months and the median delay in diagnosis was 3 months. Genetic analysis was done in 14 cases. It revealed 12 variants (11 missense and one nonsense) from 13 patients. Most of the pathogenic variants were in exon 9 and exon 5. Three novel variants were identified in three patients: two missense variants c.169G > A (p.Glu57Lys), c.1048T > C (p.Cys350Arg) and one nonsense variant c.331C > T (p.Lys111Ter). On neurodevelopmental assessment, majority of children with GA-1 were non ambulatory (62.5%), had limited hand skills (58.3%) and impaired communication (58.3%). Overall, poor global development was noted in 43.7%. A pre-existing developmental delay was significantly associated with impaired communication skills (p = 0.03), and the number of episodes of encephalopathy were significantly associated with impaired gross motor skill (p = 0.02). Presence of encephalopathy was significantly associated with poor performance in social emotional (p = 0.01) and cognitive (p = 0.03) domains of Developmental Profile-III scale and development of severe dystonia (p = 0.01).

CONCLUSION

Our findings highlight the clinical, biochemical, radiological and genetic spectrum of GA-1 in children in North India and report the presence of novel pathogenic variations.

Prediction and expression analysis of deleterious nonsynonymous SNPs of Arabidopsis ACD11 gene by combining computational algorithms and molecular docking approach

Accelerated cell death 11 (ACD11) is an autoimmune gene that suppresses pathogen infection in plants by preventing plant cells from becoming infected by any pathogen. This gene is widely known for growth inhibition, premature leaf chlorosis, and defense-related programmed cell death (PCD) in seedlings before flowering in Arabidopsis plant. Specific amino acid changes in the ACD11 protein’s highly conserved domains are linked to autoimmune symptoms including constitutive defensive responses and necrosis without pathogen awareness. The molecular aspect of the aberrant activity of the ACD11 protein is difficult to ascertain. The purpose of our study was to find the most deleterious mutation position in the ACD11 protein and correlate them with their abnormal expression pattern. Using several computational methods, we discovered PCD vulnerable single nucleotide polymorphisms (SNPs) in ACD11. We analysed the RNA-Seq data, identified the detrimental nonsynonymous SNPs (nsSNP), built genetically mutated protein structures and used molecular docking to assess the impact of mutation. Our results demonstrated that the A15T and A39D mutations in the GLTP domain were likely to be extremely detrimental mutations that inhibit the expression of the ACD11 protein domain by destabilizing its composition, as well as disrupt its catalytic effectiveness. When compared to the A15T mutant, the A39D mutant was more likely to destabilize the protein structure. In conclusion, these mutants can aid in the better understanding of the vast pool of PCD susceptibilities connected to ACD11 gene GLTP domain activation.

Non synonymous single nucleotide polymorphism (nsSNP) is a process in which amino acid sequence of a protein is altered as a result of single nucleotide alteration in the coding region (mRNA) of any living organism. Therefore, the entire protein structure, interactions and stability are altered, which may have a negative impact on living organisms. Hence, to completely comprehend this biological process, we must first solve the unresolved mutational protein structure and mutated protein interactions. The major goal of our research is to identify the most harmful mutation in our target protein structure and how it interacts within cells. However, it was discovered that only a few alterations in residues had the largest negative impact on the protein’s internal structure and also on the protein-ligand interactions. We show that based on the amino acid sequence of a protein computationally, it is feasible to discover mutational positions in the sequence, generate mutation protein structure and interactions with related ligands. Our findings show that the essential mechanisms underlying protein mutations generated by this process are identical. The capacity to correctly detect mutations from sequence allows the annotation and study of protein-ligand interactions throughout a whole organism, which might aid function prediction and gene expression.

Scrutinizing Deleterious Nonsynonymous SNPs and Their Effect on Human POLD1 Gene

POLD1 (DNA polymerase delta 1, catalytic subunit) is a protein-coding gene that encodes the large catalytic subunit (POLD1/p125) of the DNA polymerase delta (Polδ) complex. The consequence of missense or nonsynonymous SNPs (nsSNPs), which occur in the coding region of a specific gene, is the replacement of single amino acid. It may also change the structure, stability, and/or functions of the protein. Mutation in the POLD1 gene is associated with autosomal dominant predisposition to colonic adenomatous polyps, colon cancer, endometrial cancer (EDMC), breast cancer, and brain tumors. These de novo mutations in the POLD1 gene also result in autosomal dominant MDPL syndrome (mandibular hypoplasia, deafness, progeroid features, and lipodystrophy). In this study, genetic variations of POLD1 which may affect the structure and/or function were analyzed using different types of bioinformatics tools. A total of 17038 nsSNPs for POLD1 were collected from the NCBI database, among which 1317 were missense variants. Out of all missense nsSNPs, 28 were found to be deleterious functionally and structurally. Among these deleterious nsSNPs, 23 showed a conservation scale of >5, 2 were predicted to be associated with binding site formation, and one acted as a posttranslational modification site. All of them were involved in coil, extracellular structures, or helix formation, and some cause the change in size, charge, and hydrophobicity.

Dilated-Left Ventricular Non-Compaction Cardiomyopathy in a Pediatric Case with SPEG Compound Heterozygous Variants

Left Ventricular Non-Compaction (LVNC) is defined by the triad prominent myocardial trabecular meshwork, thin compacted layer, and deep intertrabecular recesses. LVNC associated with dilation is characterized by the coexistence of left ventricular dilation and systolic dysfunction. Pediatric cases with dilated-LVNC have worse outcomes than those with isolated dilated cardiomyopathy and adult patients. Herein, we report a clinical and genetic investigation using trio-based whole-exome sequencing of a pediatric case with early-onset dilated-LVNC. Compound heterozygous mutations were identified in the Striated Muscle Enriched Protein Kinase (SPEG) gene, a key regulator of cardiac calcium homeostasis. A paternally inherited mutation: SPEG; p.(Arg2470Ser) and the second variant, SPEG; p.(Pro2687Thr), is common and occurred de novo. Subsequently, Sanger sequencing was performed for the family in order to segregate the variants. Thus, the index case, his father, and both sisters carried the SPEG: p.(Arg2470Ser) variant. Only the index patient carried both SPEG variants. Both sisters, as well as the patient’s father, showed LVNC without cardiac dysfunction. The unaffected mother did not harbor any of the variants. The in silico analysis of the identified variants (rare and common) showed a decrease in protein stability with alterations of the physical properties as well as high conservation scores for the mutated residues. Interestingly, using the Project HOPE tool, the SPEG; p.(Pro2687Thr) variant is predicted to disturb the second fibronectin type III domain of the protein and may abolish its function. To our knowledge, the present case is the first description of compound heterozygous SPEG mutations involving a de novo variant and causing dilated-LVNC without neuropathy or centronuclear myopathy.

In silico analysis of missense variants of the C1qA gene related to infection and autoimmune diseases

Objectives

C1q is a key activator of the classical pathway of the complement system and exerts consequences relating to opsonization and phagocytosis. The C1qA gene is one of three genes encoding the C1q molecule. Defects in C1q, and especially in C1qA, have been linked to an increased susceptibility to infection, sepsis, and systemic lupus erythematosus. These defects could arise from missense single nucleotide polymorphisms (SNPs) and their deleterious impacts on protein structure and function. Thus, identifying high-risk missense SNPs in C1qA has become a necessity if we are to identify appropriate measures for prevention and management of affected patients.

Methods

A comprehensive in silico study was conducted to screen the 184 missense SNPs in the C1qA gene using different tools with different algorithms and approaches. We investigated the impact of SNPs on protein function, stability, and structure. In addition, we identified the location of the SNPs on protein domains, secondary structure alignment, and the phylogenetic conservation of their positions.

Results

Of the 184 missense SNPs, 10 SNPs were predicted to be the most damaging to protein function and structure.

Conclusion

Ten missense SNPs were predicted to have the highest risk of damaging protein function and structure, thus leading to infection, sepsis, and systemic lupus erythematosus. These 10 SNPs constitute the best candidates for further experimental investigations.

Investigation of GALNS variants and genotype-phenotype correlations in a large cohort of patients with mucopolysaccharidosis type IVA

Mucopolysaccharidosis type IVA (MPS IVA) is a rare autosomal recessive disorder resulting from the deficiency of N-acetylgalactosamine-6-sulfate sulfatase (GALNS) caused by pathogenic variants in the GALNS gene. A systematic analysis for genotype-phenotype correlation is essential due to hundreds of variants generating different levels of residual GALNS activity and causing a wide degree of clinical manifestation effects. Here, we retrospectively analyzed clinical and genetic data of 108 unrelated patients with MPS IVA to investigate the variants spectrum of GALNS and assess their clinical effects. In this cohort, 82 patients were classified as severe, 14 as intermediate, and 12 as mild. One hundred and one GALNS variants were identified, of which 47 were novel. Most patients with at least one GALNS null variant were classified as severe phenotype (92%, 33/36). Missense variants mapped to different residues of GALNS protein resulted in different phenotypes in patients with MPS IVA. Ninety-two percent of patients with two missense variants mapped to buried residues were classified as severe (92%, 24/26), while at least one missense variant mapped to surface residues was identified in patients with biallelic missense variants presenting intermediate MPS IVA (78%, 7/9) and presenting mild MPS IVA (86%, 6/7). Our study contributes to a better understanding of the molecular spectrum of GALNS variants and their clinical implications. Based on the data herein reported, we generated a systematic flowchart correlating the GALNS variants to assist in phenotype prediction and classification of patients with MPS IVA.

In silico analysis of a novel pathogenic variant c.7G > A in C14orf39 gene identified by WES in a Pakistani family with azoospermia

Infertility is a multifactorial disorder that affects approximately 12% of couples of childbearing ages worldwide. Few studies have been conducted to understand the genetic causes of infertility in depth. The synaptonemal complex (SC), which is essential for the progression of meiosis, is a conserved tripartite structure that binds homologous chromosomes together and is thus required for fertility. This study investigated genetic causes of infertility in a Pakistani consanguineous family containing two patients suffering from non-obstructive azoospermia (NOA). We performed whole-exome sequencing, followed by Sanger sequencing, and identified a novel pathogenic variant (c.7G > A [p.D3N]) in the SC coding gene C14orf39, which was recessively co-segregated with NOA. In silico analysis revealed that charges on wild-type residues were lost, which may result in loss of interactions with other molecules and residues, and a reduction in protein stability occurred, which was caused by the p.D3N mutation. The novel variant generated the mutant protein C14ORF39^D3N, and homozygous mutations in C14orf39 resulted in NOA. The transcriptome profile of C14ORF39 shows that it is specifically expressed in early brain development, which suggests that research in this area is required to study other functions of C14ORF39 in addition to its role in the germline. This research highlights the conserved role of C14orf39/SIX6OS1 in assembly of the SC and its indispensable role in facilitating genetic diagnosis in patients with infertility, which may enable the development of future treatments.

Novel Disease-Associated Missense Single-Nucleotide Polymorphisms Variants Predication by Algorithms Tools and Molecular Dynamics Simulation of Human TCIRG1 Gene Causing Congenital Neutropenia and Osteopetrosis

T Cell Immune Regulator 1, ATPase H + Transporting V0 Subunit A3 (TCIRG1 gene provides instructions for making one part, the a3 subunit, of a large protein complex known as a vacuolar H + -ATPase (V-ATPase). V-ATPases are a group of similar complexes that act as pumps to move positively charged hydrogen atoms (protons) across membranes. Single amino acid changes in highly conserved areas of the TCIRG1 protein have been linked to autosomal recessive osteopetrosis and severe congenital neutropenia. We used multiple computational approaches to classify disease-prone single nucleotide polymorphisms (SNPs) in TCIRG1. We used molecular dynamics analysis to identify the deleterious nsSNPs, build mutant protein structures, and assess the impact of mutation. Our results show that fifteen nsSNPs (rs199902030, rs200149541, rs372499913, rs267605221, rs374941368, rs375717418, rs80008675, rs149792489, rs116675104, rs121908250, rs121908251, rs121908251, rs149792489 and rs116675104) variants are likely to be highly deleterious mutations as by incorporating them into wild protein they destabilize the wild protein structure and function. They are also located in the V-ATPase I domain, which may destabilize the structure and impair TCIRG1 protein activation, as well as reduce its ATPase effectiveness. These mutants have not yet been identified in patients suffering from CN and osteopetrosis while (G405R, R444L, and D517N) reported in our study are already associated with osteopetrosis. Mutation V52L reported in our study was identified in a patient suspected for CN. Finally, these mutants can help to further understand the broad pool of illness susceptibilities associated with TCIRG1 catalytic kinase domain activation and aid in the development of an effective treatment for associated diseases.

In silico screening and analysis of single-nucleotide polymorphic variants of the ABCC2 gene affecting Dubin-Johnson syndrome

BACKGROUND AND STUDY AIMS

Dubin-Johnson syndrome (DJS) is a benevolent genetic disorder of the liver with autosomal inheritance. It is a rare disorder characterized by an increase in conjugated bilirubin and anomaly in coproporphyrin clearance. DJS is caused by deleterious mutations in the ABCC2 gene. A polymorphism in the ABCC2 gene causes malfunctions in its ability to regulate the efflux of different organic anions, such as bilirubin, from hepatocytes to the canaliculi. Multidrug resistance protein 2 (MRP2) encoded by the ABCC2 gene is one of the main regulators of the export of bilirubin to respective sites. ABCC2 gene mutations have widely drawn attention in the pathology of DJS in various populations.

PATIENTS AND METHODS

The ABCC2 gene was subjected to the National Center for Biotechnology Information (NCBI) database in 2020, and non-synonymous single-nucleotide polymorphisms (nsSNPs) and variants in untranslated regions were studied using different computational servers. SIFT, Protein variation effect analyzer, and PolyPhen-2 were used to retrieve the damaging Single-nucleotide polymorphisms (SNPs); PhD-SNP, SNPs&GO, and Protein Analysis Through Evolutionary Relationships were used to predict the association of nsSNPs with DJS; Mutation3D illustrated the location of variants in the protein; SNAP2, MutPred2, ELASPIC, and HOPE were used to predict the structural and functional effects of these mutations on MRP2; and I-mutant 3.0 and MuPro were used to determine the effects of polymorphism on the function of MRP2.

RESULTS

In this study, 18,947 SNPs were screened from the NCBI database, followed by a series of refinement of variants using online available servers. We concluded that 41 ABCC2 gene variants are vital etiological candidates for DJS in humans. These 41 variants had highly damaging effects on the MRP2 protein, which may lead to deficient transportation capacity, thereby affecting the efflux of bilirubin across the canalicular membrane.

CONCLUSION

In silico tools are an alternative approach for predicting the target SNPs. Hence, previously unreported variants can be considered strong etiological candidates for diseases related to MRP2.

Identification of the most damaging nsSNPs in the human CFL1 gene and their functional and structural impacts on cofilin-1 protein

The cofilin-1 protein, encoded by CFL1, is an actin-binding protein that regulates F-actin depolymerization and nucleation activity through phosphorylation and dephosphorylation. CFL1 has been implicated in the development of neurodegenerative diseases (Alzheimer's disease and Huntington's disease), neuronal migration disorders (lissencephaly, epilepsy, and schizophrenia), and neural tube closure defects. Mutations in CFL1 have been associated with impaired neural crest cell migration and neural tube closure defects. In our study, various computational approaches were utilized to explore single-nucleotide polymorphisms (SNPs) in CFL1. The Variation Viewer and gnomAD databases were used to retrieve CFL1 SNPs, including 46 nonsynonymous SNPs (nsSNPs). The functional and structural annotation of SNPs was performed using 12 sequence-based web applications, which identified 20 nsSNPs as being the most likely to be deleterious or disease-causing. The conservation of cofilin-1 protein structures was illustrated using the ConSurf and PROSITE web servers, which projected the 12 most deleterious nsSNPs onto conserved domains, with the potential to disrupt the protein's functionality. These 12 nsSNPs were selected for protein structure construction, and the DynaMut/DUET servers predicted that the protein variants V7G, L84P, and L99A were the most likely to be damaging to the cofilin-1 protein structure or function. The evaluation of molecular docking studies demonstrated that the L99A and L84P cofilin-1 variants reduce the binding affinity for actin compared with the native cofilin-1 structure, and molecular dynamic simulation studies confirmed that these variants might destabilize the protein structure. The consequences of putative mutations on protein-protein interactions and post-translational modification sites in the cofilin-1 protein structure were analyzed. This study represents the first complete approach to understanding the effects of nsSNPs within the actin-depolymerizing factor/cofilin family, which suggested that SNPs resulting in L84P (rs199716082) and L99A (rs267603119) variants represent significant CFL1 mutations associated with disease development.

Mutation profile of neurodegenerative mitochondriopathy - LHON in Southern India

BACKGROUND

Leber's Hereditary Optic Neuropathy (LHON) is a rare mitochondriopathy causing retinal ganglion cell degeneration resulting in central vision loss. It is caused by mitochondrial DNA (mtDNA) mutations and thus follows maternal inheritance pattern.

METHODS

We analysed the whole mitochondrial genome in 100 South Indian LHON patients by utilizing Sanger and Next Generation Sequencing approaches. Haplogroup analysis was performed using HaploGrep2 to predict the risk group. Methylation changes in the mtDNA D-loop region were investigated by performing methylation-specific polymerase chain reaction (MSP).

RESULTS

LHON associated mutations were detected in 55% of the patients of which 42% harboured the primary mutations and 13% harboured potentially pathogenic variants that were previously reported to cause LHON. The candidate mutations identified with confirmed pathogenicity are: m.11778G > A (38%), m.14484 T > C (3%), m.4171C > A (1%) and m.11696G > A (1%). MSP results demonstrated that the D-loop region was unmethylated in all the study subjects including mutation-positive patients, mutation-negative patients, asymptomatic carriers, and controls. Haplogroup-M was prevalent (69%) in the study cohort followed by R (14%), U (9%), N (3%), HV (2%), G (2%), and W (1%). The frequency of the predominant mutation m.11778G > A was found lower (̴ 11%) in haplogroup-U.

CONCLUSIONS

South Indian LHON cohort shows a unique profile of mtDNA mutations and haplogroup association presumably with no role of D-loop methylation. MT-ND4, MT-ND5, and MT-ND1 serve as the hotspot genes in this cohort. The presence of LHON associated mutations in patients lacking the common primary mutations insists on the necessity of mitochondrial genome sequencing in individuals suspected with LHON.

Co-segregation of variant NSUN2 Lue198Arg among Iranian family with intellectual disability: a case report

Background

Intellectual disability is characterized by impairments in adaptive behavior and cognitive functioning manifested during the developmental period. Since disabilities are heterogeneous, variant analysis can help us confirm and accurately diagnose children with intellectual disabilities. Some papers reported that bi-allelic variants of the NSUN2 gene caused a group of neurological disorders, including non-syndromic autosomal recessive intellectual disability (NS-ARID), Dubowitz syndrome, and familial restrictive cardiomyopathy 1 (RCM1). We report on a consanguineous family with three siblings diagnosed with intellectual disability.

Case presentation

The 7-year-old female was referred to Ali-Asghar hospital, Zahedan, Iran, with clinical manifestations comprising moderate intellectual disability, ptosis, long face, and short stature. Chromosome banding, metabolic testing, and magnetic resonance imaging examinations revealed no abnormalities. Accordingly, other affected siblings born of the same parents were considered. Whole-exome sequencing (WES) was conducted on the sufferer to consider NS-ARID variants. Findings identified a variant with uncertain significance (NM_017755.6: c.593 T > G) in the NSUN2 gene in the proband. This variant was confirmed through Sanger sequencing of the affected and unaffected family members. Besides, the computational results showed that the L198R exchange could change the interaction between wild-type and other residues in the protein. The affected patients with NS-ARID had similar clinical characteristics and genetic abnormalities.

Conclusion

Taken together, we described the variant in three Iranian siblings; further expanding of the other variants involved in the disease will be evident by using high-throughput sequencing technologies.

Novel biallelic variants expand the SLC5A6-related phenotypic spectrum

The sodium (Na+):multivitamin transporter (SMVT), encoded by SLC5A6, belongs to the sodium:solute symporter family and is required for the Na+-dependent uptake of biotin (vitamin B7), pantothenic acid (vitamin B5), the vitamin-like substance α-lipoic acid, and iodide. Compound heterozygous SLC5A6 variants have been reported in individuals with variable multisystemic disorder, including failure to thrive, developmental delay, seizures, cerebral palsy, brain atrophy, gastrointestinal problems, immunodeficiency, and/or osteopenia. We expand the phenotypic spectrum associated with biallelic SLC5A6 variants affecting function by reporting five individuals from three families with motor neuropathies. We identified the homozygous variant c.1285 A > G [p.(Ser429Gly)] in three affected siblings and a simplex patient and the maternally inherited c.280 C > T [p.(Arg94*)] variant and the paternally inherited c.485 A > G [p.(Tyr162Cys)] variant in the simplex patient of the third family. Both missense variants were predicted to affect function by in silico tools. 3D homology modeling of the human SMVT revealed 13 transmembrane helices (TMs) and Tyr162 and Ser429 to be located at the cytoplasmic facing region of TM4 and within TM11, respectively. The SLC5A6 missense variants p.(Tyr162Cys) and p.(Ser429Gly) did not affect plasma membrane localization of the ectopically expressed multivitamin transporter suggesting reduced but not abolished function, such as lower catalytic activity. Targeted therapeutic intervention yielded clinical improvement in four of the five patients. Early molecular diagnosis by exome sequencing is essential for timely replacement therapy in affected individuals.

Identification, analysis of deleterious SNPs of the human GSR gene and their effects on the structure and functions of associated proteins and other diseases

Hereditary glutathione reductase deficiency, caused by mutations of the GSR gene, is an autosomal recessive disorder characterized by decreased glutathione disulfide (GSSG) reduction activity and increased thermal instability. This study implemented computational analysis to screen the most likely mutation that might be associated with hereditary glutathione reductase deficiency and other diseases. Using ten online computational tools, the study revealed four nsSNPs among the 17 nsSNPs identified as most deleterious and disease associated. Structural analyses and evolutionary confirmation study of native and mutant GSR proteins using the HOPE project and ConSruf. HOPE revealed more flexibility in the native GSR structure than in the mutant structure. The mutation in GSR might be responsible for changes in the structural conformation and function of the GSR protein and might also play a significant role in inducing hereditary glutathione reductase deficiency. LD and haplotype studies of the gene revealed that the identified variations rs2978663 and rs8190955 may be responsible for obstructive heart defects (OHDs) and hereditary anemia, respectively. These interethnic differences in the frequencies of SNPs and haplotypes might help explain the unpredictability that has been reported in association studies and can contribute to predicting the pharmacokinetics and pharmacodynamics of drugs that make use of GSR.

Analysis of the APOB Gene and Apolipoprotein B Serum Levels in a Mexican Population with Acute Coronary Syndrome: Association with the Single Nucleotide Variants rs1469513, rs673548, rs676210, and rs1042034

Apolipoprotein B (APOB) is associated with the development of atherosclerosis and consequently in the acute coronary syndrome (ACS) physiopathology. Single number variants (SNVs) in apolipoprotein B gene (APOB) influence over the susceptibility for this syndrome. The aim of this study was to determine the impact of the rs1469513, rs673548, rs676210, and rs1042034 SNVs and serum levels of APOB in the risk of ACS in a population from western Mexico. We included 300 patients in the group of cases (ACSG) and 300 individuals in the control group (CG). APOB levels were evaluated by immunonephelometry, and SNVs were genotyped with TaqMan probes. We found significant allelic and genotypic differences between groups for rs673548 and rs676210 (OR = 1.33, p=0.030, OR = 2.69, p < 0.001) and rs1042034 (OR = 0.50, p=0.037) SNVs. We found a risk haplotype TAGT (OR: 2.14, IC 1.50–3.04, p < 0.001). Our findings support a significant risk association between rs673548 and rs676210 variants for ACS; meanwhile, rs1042034 could be considered protective factor in a western Mexican population. Also, in this population, haplotype TAGT may confer 2.14 times a higher risk. APOB serum levels were compared by genotype variants in both groups without any significant statistical difference.

Identification and In Silico Characterization of a Novel COLGALT2 Gene Variant in a Child with Mucosal Rectal Prolapse

Rectal prolapse is influenced by many factors including connective tissue dysfunction. Currently, there is no data about genetic contribution in the etiology of this disorder. In this study, we performed trio whole-exome sequencing in an 11-year-old girl with mucosal rectal prolapse and her parents and sibling. Genetic testing revealed a novel heterozygous missense variant c.1406G>T; p.G469V in exon 11 of the COLGALT2 gene encoding the GLT25 D2 enzyme. Sanger sequencing confirmed this variant only in the patient while the mother, father and sister showed a wild-type sequence. The pathogenicity of the novel variant was predicted using 10 different in silico tools that classified it as pathogenic. Further, in silico prediction, according to Phyre2, Project HOPE, I-Mutant3.0 and MutPred2 showed that the missense variant can decrease protein stability and cause alterations in the physical properties of amino acids resulting in structural and functional changes of the GLT25D2 protein. In conclusion, the present study identifies a previously unknown missense mutation in the COLGALT2 gene that encodes the enzyme involved in collagen glycosylation. The clinical features observed in the patient and the results of in silico analysis suggest that the new genetic variant can be pathogenic.

The Genomic Architecture of Hidradenitis Suppurativa-A Systematic Review

Hidradenitis suppurativa is a chronic, suppurative condition of the pilosebaceous unit manifesting as painful nodules, abscesses, and sinus tracts mostly in, but not limited to, intertriginous skin. Great strides have been made at elucidating the pathophysiology of hidradenitis suppurativa, which appears to be the product of hyperkeratinization and inflammation brought about by environmental factors and a genetic predisposition. The identification of familial hidradenitis suppurativa has sparked research aimed at identifying underlying pathogenic variants in patients who harbor them. The objective of this review is to provide a broad overview of the role of genetics in various aspects of hidradenitis suppurativa, specifically the pathophysiology, diagnosis, and clinical application.

Unraveling structural and conformational dynamics of DGAT1 missense nsSNPs in dairy cattle

Cattle are domestic animals that have been nourishing humans for thousands of years. Milk from cattle represents a key source of high-quality protein, fat, and other nutrients. The nutritional value of milk and dairy products is closely associated with the fat content, providing up to 30% of the total fat consumed in the human diet. The fat content in cattle milk represents a major concern for the scientific community due to its association with human health. The relationship between milk fat content and diacylglycerol o-acyltransferase 1 gene (DGAT1) is well described in literature. Several studies demonstrated the difference in fat contents and other milk production traits in a wide range of cattle breeds, to be associated with missense non-synonymous single nucleotide polymorphisms (nsSNPs) of the DGAT1 gene. As a result, an nsSNPs analysis is crucial for unraveling the DGAT1 structural and conformational dynamics linked to milk fat content. DGAT1-nsSNPs are yet to be studied in terms of their structural and functional impact. Therefore, state-of-the-art computational and structural genomic methods were used to analyze five selected variants (W128R, W214R, C215G, P245R, and W459G), along with the wild type DGAT1. Significant structural and conformational changes in the variants were observed. We illustrate how single amino acid substitutions affect DGAT1 function, how this contributes to our understanding of the molecular basis of variations in DGAT1, and ultimately its impact in improving fat quality in milk.

Effect of Cis-Compound Variants in MYH7 on Hypertrophic Cardiomyopathy With a Mild Phenotype

Patients with hypertrophic cardiomyopathy (HC) caused by compound variants have severe clinical manifestations, but significant clinical heterogeneity remains. Clinical diversity in these patients may result from different combinations of variants. We analyzed the role of cis-compound variants in a Chinese HC pedigree. Exome sequencing was performed in the proband. Identified variants were detected with bi-directional Sanger sequencing in a pedigree that comprised 3 generations and 28 family members. Follow-up was performed for 16 years. Two missense variants (c.2465T>C, p.Met822Thr; c.4258C>T, p.Arg1420Trp) were identified in the MYH7 gene. These variants were absent in our 761 in-house people without HC and predicted to be pathogenic.Both variants were detected in 11 family members, thus they were believed to inherit cis. In the 11 members, only 5 developed HC, the other 6 were asymptomatic variant carriers with an abnormal electrocardiogram. The HC members had mild hypertrophy with a maximum left ventricular wall thickness of 13 to 21 mm and showed a low incidence of cardiovascular events. In conclusion, the cis-compound variants of Met822Thr and Arg1420Trp in MYH7 are causal but relatively benign, variants associated with familial HC. This finding suggests that different types of compound variants might need to be analyzed for a genotype-phenotype study.

Structural modeling of human AKAP3 protein and in silico analysis of single nucleotide polymorphisms associated with sperm motility

AKAP3 is a member of the A-kinase anchoring proteins and it is a constituent of the sperm fibrous sheath. AKAP3 is needed for the formation of sperm flagellum structure, sperm motility, and male fertility. This study aims to model the AKAP3 tertiary structure and identify the probable impact of four mutations characterized in infertile men on the AKAP3 structure. The T464S, I500T, E525K, and I661T substitutions were analyzed using in silico methods. The secondary structure and three-dimensional model of AKAP3 were determined using PSI-BLAST based secondary structure prediction and Robetta servers. The TM-score was used to quantitatively measure the structural similarities between native and mutated models. All of the desired substitutions were classified as benign. I-Mutant results showed all of the substitutions decreased AKAP3 stability; however, the I500T and I661T were more effective. Superposition and secondary structure comparisons between native and mutants showed no dramatic deviations. Our study provided an appropriate model for AKAP3. Destabilization of AKAP3 caused by these substitutions did not appear to induce structural disturbances. As AKAP3 is involved in male infertility, providing more structural insights and the impact of mutations that cause protein functional diversity could elucidate the etiology of male fertility problems at molecular level.

Tryptophan depletion results in tryptophan-to-phenylalanine substitutants

Activated T cells secrete interferon-γ, which triggers intracellular tryptophan shortage by upregulating the indoleamine 2,3-dioxygenase 1 (IDO1) enzyme^1-4. Here we show that despite tryptophan depletion, in-frame protein synthesis continues across tryptophan codons. We identified tryptophan-to-phenylalanine codon reassignment (W>F) as the major event facilitating this process, and pinpointed tryptophanyl-tRNA synthetase (WARS1) as its source. We call these W>F peptides 'substitutants' to distinguish them from genetically encoded mutants. Using large-scale proteomics analyses, we demonstrate W>F substitutants to be highly abundant in multiple cancer types. W>F substitutants were enriched in tumours relative to matching adjacent normal tissues, and were associated with increased IDO1 expression, oncogenic signalling and the tumour-immune microenvironment. Functionally, W>F substitutants can impair protein activity, but also expand the landscape of antigens presented at the cell surface to activate T cell responses. Thus, substitutants are generated by an alternative decoding mechanism with potential effects on gene function and tumour immunoreactivity.

CNTNAP1-encephalopathy: Six novel patients surviving the neonatal period

CNTNAP1 encodes CASPR1, involved in the paranodal junction. Thirty-three patients, with CNTNAP1 biallelic mutations have been described previously. Most of them had a very severe neurological impairment and passed away in the first months of life. We identified four patients, from two unrelated families, who survived over the neonatal period. Exome sequencing showed compound heterozygous or homozygous variants. Severe hypotonia was a constant feature. When compared to previous reports, the most important clinical differences observed in our patients were the absence of antenatal problems and, in two of them, the lack of respiratory distress. Less commonly reported characteristics such as epileptic seizures, dystonia, and impaired communication skills were also observed. MRIs revealed hypomyelination or abnormal white matter signal, cerebral or cerebellar atrophy. The present observations support a wider than initially reported clinical spectrum, including survival after the neonatal period and additional neurological features. They contribute to better delineate the phenotype-genotype correlations for CNTNAP1. In addition, we report one more family with two sibs who carry a missense variant of uncertain significance which we propose could be associated with a milder phenotype.

Endoglin Wild Type and Variants Associated With Hereditary Hemorrhagic Telangiectasia Type 1 Undergo Distinct Cellular Degradation Pathways

Endoglin, also known as cluster of differentiation 105 (CD105), is an auxiliary receptor in the TGFβ signaling pathway. It is predominantly expressed in endothelial cells as a component of the heterotetrameric receptor dimers comprising type I, type II receptors and the binding ligands. Mutations in the gene encoding Endoglin (ENG) have been associated with hereditary hemorrhagic telangiectasia type 1 (HHT1), an autosomal dominant inherited disease that is generally characterized by vascular malformation. Secretory and many endomembrane proteins synthesized in the Endoplasmic reticulum (ER) are subjected to stringent quality control mechanisms to ensure that only properly folded and assembled proteins are trafficked forward through the secretory pathway to their sites of action. We have previously demonstrated that some Endoglin variants causing HHT1 are trapped in the ER and fail to traffic to their normal localization in plasma membrane, which suggested the possible involvement of ER associated protein degradation (ERAD) in their molecular pathology. In this study, we have investigated, for the first time, the degradation routes of Endoglin wild type and two mutant variants, P165L and V105D, and previously shown to be retained in the ER. Stably transfected HEK293 cells were treated with proteasomal and lysosomal inhibitors in order to elucidate the exact molecular mechanisms underlying the loss of function phenotype associated with these variants. Our results have shown that wild type Endoglin has a relatively short half-life of less than 2 hours and degrades through both the lysosomal and proteasomal pathways, whereas the two mutant disease-causing variants show high stability and predominantly degrades through the proteasomal pathway. Furthermore, we have demonstrated that Endoglin variants P165L and V105D are significantly accumulated in HEK293 cells deficient in HRD1 E3 ubiquitin ligase; a major ERAD component. These results implicate the ERAD mechanism in the pathology of HHT1 caused by the two variants. It is expected that these results will pave the way for more in-depth research studies that could provide new windows for future therapeutic interventions.

A novel homozygous missense variant p.D339N in the PKLR gene correlates with pyruvate kinase deficiency in a Pakistani family: a case report

Background

Pyruvate kinase deficiency is an exceptionally rare autosomal recessive Mendelian disorder caused by bi-allelic pathogenic variants in the PKLR gene. It is mainly characterized by chronic nonspherocytic hemolytic anemia though other symptoms such as splenomegaly, hepatomegaly, pallor, fatigue, iron overload, shortness of breath, hyperbilirubinemia, and gallstones might also prevail.

Case presentation

We present here a novel genetic defect in the PKLR gene that correlates with pyruvate kinase deficiency phenotype in a consanguineous family from North-Western Pakistan. The family included three affected individuals who were all born to consanguineous parents. The proband, a 13-year-old female of Pashtun ethnicity, showed chronic nonautoimmune hemolytic anemia since birth, extremely low hemoglobin (7.6 g/dL) and pyruvate kinase (12.4 U/g Hb) levels, splenomegaly, and hepatomegaly. Bone marrow aspirate showed a markedly decreased myeloid to erythroid ratio and hypercellular marrow particles due to hyperplasia of the erythroid elements. Molecular characterization of the proband’s genomic DNA uncovered a likely pathogenic homozygous missense variant p.[D339N] in exon 7 of the PKLR gene. In-depth in silico analysis and familial cosegregation implies p.[D339N] as the likely cause of pyruvate kinase deficiency in this family. Further in vitro or in vivo studies are required to validate the impact of p.[D339N] on protein structure and/or stability, and to determine its role in the disease pathophysiology.

Conclusions

In summary, these findings suggest a novel genetic defect in the PKLR gene as a likely cause of pyruvate kinase deficiency, thus further expanding the mutational landscape of this rare Mendelian disorder.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13256-022-03292-z.

Role of MBL2 Polymorphisms in Sepsis and Survival: A Pilot Study and In Silico Analysis

Sepsis is a serious infection-induced syndrome with serious ramifications, especially in intensive care units. Global concern motivated the investigation of the role of related genes’ polymorphism in predicting the liability to infection, sepsis, septic shock and survival. Among these genes is the gene encoding mannose-binding lectin (MBL), with its remarkable importance in the immune system. However, the previous studies showed conflicting results and ambiguity that urged us to engage with this issue in the Egyptian population. Prediction of functional and structural impacts of single nucleotide polymorphisms (SNPs) was done using in silico methods. A prospective observational study was conducted in intensive care units; one hundred and thirty patients were followed up. Genotyping was performed using real-time polymerase chain reaction (RT-PCR) technology. MBL SNPs showed a remarkable high frequency in our population, as well. No significant association was found between MBL2 genotypes and any of our analyses (sepsis, septic shock and survival). Only septic shock and age were independently associated with time of survival by Cox regression analysis. Our study may confirm the redundancy of MBL and the absence of significant impact on sepsis liability and mortality in adult patients.

Genotype-phenotype spectrum of 130 unrelated Indian families with Mucopolysaccharidosis type II

MPS II is an X linked recessive lysosomal storage disorder with multi-system involvement and marked molecular heterogeneity. In this study, we explored the clinical and molecular spectrum of 144 Indian patients with MPS II from 130 unrelated families. Clinical information was collected on a predesigned clinical proforma. Sanger method was employed to sequence all the exons and exon/intron boundaries of the IDS gene. In cases where causative variation was not detected by Sanger sequencing, MLPA and RFLP were performed to identify large deletions/duplications and complex rearrangements. Cytogenetic microarray was done in one patient to see the breakpoints and extent of deletion. In one patient with no detectable likely pathogenic or pathogenic variation, whole-genome sequencing was also performed. Novel variants were systematically assessed by in silico prediction software and protein modelling. The pathogenicity of variants was established based on ACMG criteria. An attempt was also made to establish a genotype-phenotype correlation. Positive family history was present in 31% (41/130) of patients. Developmental delay and intellectual disability were the main reasons for referral. Macrocephaly, coarse facies and dysostosis were present in almost all patients. Hepatosplenomegaly, joint contractures and short stature were the characteristic features, seen in 87% (101/116), 67.8% (74/109) and 41.4% (41/99) patients respectively. Attenuated phenotype was seen in 32.6% (47/144) patients, while severe phenotype was seen in 63% (91/144) patients. The detection rate for likely pathogenic or pathogenic variants in our cohort is 95.5% (107/112) by Sanger sequencing, MLPA and RFLP. We also found two variants of unknown significance, one each by Sanger sequencing and WGS. Total of 71 variants were identified by Sanger sequencing and 29 of these variants were found to be novel. Amongst the novel variants, there was a considerable proportion (51%) of frameshift variants (15/29). Almost half of the causative variants were located in exon 3,8 and 9. A significant genotype-phenotype correlation was also noted for both known and novel variants. This information about the genotype spectrum and phenotype will be helpful for diagnostic and prognostic purposes.

Evolution of Fertilization-Related Genes Provides Insights Into Reproductive Health in Natural Ascrotal Mammals

Cryptorchidism is the failure of one or both testes to descend into the bottom of the scrotum. This common congenital malformation in humans and domestic animals is the best characterized risk factor for abnormal sperm functions and infertility. However, current treatment approaches for cryptorchidism do not ensure paternity in all cases. Some lineages of mammals (such as elephants and cetaceans) have natural ascrotal testes (i.e., undescended or incompletely descended testes) and normal sperm motility and fertility, providing an opportunity to understand the genetic basis of cryptorchidism. In this study, we showed that genes associated with sperm motility and competition/fertility in ascrotal mammals experienced frequent, strong selective pressure. The fixation of specific amino acids and positive selection in ascrotal mammals could affect the physicochemical properties and functions of fertilization-related proteins. In a comparison between mammals with undescended testes and incompletely descended testes, discrepancies in genes showing evidence for adaptive evolution and in functional enrichment suggested that multiple molecular mechanisms contribute to the maintenance of fertility in the challenging testicular environment. Our findings revealed substantial heterogeneity in the divergence of fertilization-related genes between natural scrotal and ascrotal mammals and provide insight into molecular mechanisms underlying normal sperm motility and competition in natural ascrotal mammals. We provide a detailed theoretical basis for understanding the pathology of cryptorchidism from a molecular evolutionary perspective. This study may contribute to the establishment of diagnostic and therapeutic targets for sperm motility and fertility disorders due to congenital cryptorchidism in humans and domestic animals.

Molecular insights into the interaction of HPV-16 E6 variants against MAGI-1 PDZ1 domain

Oncogenic protein E6 from Human Papilloma Virus 16 (HPV-16) mediates the degradation of Membrane-associated guanylate kinase with inverted domain structure-1 (MAGI-1), throughout the interaction of its protein binding motif (PBM) with the Discs-large homologous regions 1 (PDZ1) domain of MAG1-1. Generic variation in the E6 gene that translates to changes in the protein’s amino acidic sequence modifies the interaction of E6 with the cellular protein MAGI-1. MAGI-1 is a scaffolding protein found at tight junctions of epithelial cells, where it interacts with a variety of proteins regulating signaling pathways. MAGI-1 is a multidomain protein containing two WW (rsp-domain-9), one guanylate kinase-like, and six PDZ domains. PDZ domains played an important role in the function of MAGI-1 and served as targets for several viral proteins including the HPV-16 E6. The aim of this work was to evaluate, with an in silico approach, employing molecular dynamics simulation and protein–protein docking, the interaction of the intragenic variants E-G350 (L83V), E-C188/G350 (E29Q/L83V), E-A176/G350 (D25N/L83V), E6-AAa (Q14H/H78Y/83V) y E6-AAc (Q14H/I27RH78Y/L83V) and E6-reference of HPV-16 with MAGI-1. We found that variants E-G350, E-C188/G350, E-A176/G350, AAa and AAc increase their affinity to our two models of MAGI-1 compared to E6-reference.