Predicting the functional consequences of non-synonymous single nucleotide polymorphisms in IL8 gene

Rapid discrimination between deleterious and benign missense mutations in the CAGI 6 experiment

We describe the machine learning tool that we applied in the CAGI 6 experiment to predict whether single residue mutations in proteins are deleterious or benign. This tool was trained using only single sequences, i.e., without multiple sequence alignments or structural information. Instead, we used global characterizations of the protein sequence. Training and testing data for human gene mutations was obtained from ClinVar (ncbi.nlm.nih.gov/pub/ClinVar/), and for non-human gene mutations from Uniprot (www.uniprot.org). Testing was done on post-training data from ClinVar. This testing yielded high AUC and Matthews correlation coefficient (MCC) for well trained examples but low generalizability. For genes with either sparse or unbalanced training data, the prediction accuracy is poor. The resulting prediction server is available online at http://www.mamiris.com/Shoni.cagi6.

Analysis of potential genes, immunological and antioxidant profiles associated with trypanosomiasis susceptibility in dromedary camels

Trypanosomiasis is associated with tissue damage and may trigger an immunological response. These tissue lesions are linked to metabolic issues and oxidative stress. The current study aimed to investigate the immunological, antioxidant, and metabolic changes that may be connected to camel trypanosomiasis. Blood samples were collected from 54 camels and allocated into two groups: The control group (35 camels) and the infected group (19 camels). The genes TLR2, TLR5, IL-17, MARCHF3, RASGRP1, EPS15L1, PPIE, ASB16, CMPK2, LPCAT1, FPGT, GPHN, TNNI3K, DIO3, keap1, and OXSR1 were significantly up-regulated in trypanosomiasis camels. However, down-regulation was observed for the genes Nrf2, PRDX6, and NDUFS5. PCR-DNA sequencing was used to identify nucleotide sequence polymorphisms in the immune (TLR2, TLR5, IL-17, MARCHF3, RASGRP1, and EPS15L1), metabolic (PPIE, ASB16, CMPK2, LPCAT1, FPGT, GPHN, TNNI3K, and DIO3), and antioxidant (Nrf2, Keap1, PRDX6, NDUFS5, and OXSR1) genes between healthy and trypanosomiasis-affected camels. Exploring the serum profile also showed a significant (P ˂ 0.05) increase in Hp, SAA, Cp, IL-1β, IL-6, IL 10, TNF-α, and MDA, with significant (P ˂ 0.05) reduction in the serum levels of CAT, SOD, GSH, T3, and T4 in diseased camels compared with healthy ones. Our findings confirm the significance of nucleotide variations, gene expression patterns, and the biochemical profile of the investigated markers as indicators for the susceptibility of trypanosomiasis in dromedary camels and may be utilized to create management strategies.

In-silico identification of deleterious non-synonymous SNPs of TBX1 gene: Functional and structural impact towards 22q11.2DS

The TBX1 gene plays a critical role in the development of 22q11.2 deletion syndrome (22q11.2DS), a complex genetic disorder associated with various phenotypic manifestations. In this study, we performed in-silico analysis to identify potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) within the TBX1 gene and evaluate their functional and structural impact on 22q11.2DS. A comprehensive analysis pipeline involving multiple computational tools was employed to predict the pathogenicity of nsSNPs. This study assessed protein stability and explored potential alterations in protein-protein interactions. The results revealed the rs751339103(C>A), rs780800634(G>A), rs1936727304(T>C), rs1223320618(G>A), rs1248532217(T>C), rs1294927055 (C>T), rs1331240435 (A>G, rs1601289406 (A>C), rs1936726164 (G>A), and rs911796187(G>A) with a high-risk potential for affecting protein function and stability. These nsSNPs were further analyzed for their impact on post-translational modifications and structural characteristics, indicating their potential disruption of molecular pathways associated with TBX1 and its interacting partners. These findings provide a foundation for further experimental studies and elucidation of potential therapeutic targets and personalized treatment approaches for individuals affected by 22q11.2DS.

Prediction of the most deleterious non-synonymous SNPs in the human IL1B gene: evidence from bioinformatics analyses

BACKGROUND

Polymorphisms in IL1B play a significant role in depression, multiple inflammatory-associated disorders, and susceptibility to infection. Functional non-synonymous SNPs (nsSNPs) result in changes in the encoded amino acids, potentially leading to structural and functional alterations in the mutant proteins. So far, most genetic studies have concentrated on SNPs located in the IL1B promoter region, without addressing nsSNPs and their association with multifactorial diseases. Therefore, this study aimed to explore the impact of deleterious nsSNPs retrieved from the dbSNP database on the structure and functions of the IL1B protein.

RESULTS

Six web servers (SIFT, PolyPhen-2, PROVEAN, SNPs&GO, PHD-SNP, PANTHER) were used to analyze the impact of 222 missense SNPs on the function and structure of IL1B protein. Five novel nsSNPs (E100K, T240I, S53Y, D128Y, and F228S) were found to be deleterious and had a mutational impact on the structure and function of the IL1B protein. The I-mutant v2.0 and MUPro servers predicted that these mutations decreased the stability of the IL1B protein. Additionally, these five mutations were found to be conserved, underscoring their significance in protein structure and function. Three of them (T240I, D128Y, and F228S) were predicted to be cancer-causing nsSNPs. To analyze the behavior of the mutant structures under physiological conditions, we conducted a 50 ns molecular dynamics simulation using the WebGro online tool. Our findings indicate that the mutant values differ from those of the IL1B wild type in terms of RMSD, RMSF, Rg, SASA, and the number of hydrogen bonds.

CONCLUSIONS

This study provides valuable insights into nsSNPs located in the coding regions of IL1B, which lead to direct deleterious effects on the functional and structural aspects of the IL1B protein. Thus, these nsSNPs could be considered significant candidates in the pathogenesis of disorders caused by IL1B dysfunction, contributing to effective drug discovery and the development of precision medications. Thorough research and wet lab experiments are required to verify our findings. Moreover, bioinformatic tools were found valuable in the prediction of deleterious nsSNPs.

Resistance to African swine fever virus among African domestic pigs appears to be associated with a distinct polymorphic signature in the RelA gene and upregulation of RelA transcription

African swine fever virus (ASFV) is a highly contagious and fatal hemorrhagic disease of domestic pigs, which poses a major threat to the swine industry worldwide. Studies have shown that indigenous African pigs tolerate ASFV infection better than European pigs. The porcine v-rel avian reticuloendotheliosis viral oncogene homolog A (RelA) encoding a p65 kD protein, a major subunit of the NF-kB transcription factor, plays important roles in controlling both innate and adaptive immunity during infection with ASFV. In the present study, RelA genes from ASFV-surviving and symptomatic pigs were sequenced and found to contain polymorphisms revealing two discrete RelA amino acid sequences. One was found in the surviving pigs, and the other in symptomatic pigs. In total, 16 nonsynonymous SNPs (nsSNPs) resulting in codon changes were identified using bioinformatics software (SIFT and Polyphen v2) and web-based tools (MutPre and PredictSNP). Seven nsSNPs (P374-S, T448-S, P462-R, V464-P, Q478-H, L495-E, and P499-Q) were predicted to alter RelA protein function and stability, while 5 of these (P374-S, T448-S, P462-R, L495-E, and Q499-P) were predicted as disease-related SNPs.Additionally, the inflammatory cytokine levels of IFN-α, IL-10, and TNF-α at both the protein and the mRNA transcript levels were measured using ELISA and Real-Time PCR, respectively. The resulting data was used in correlation analysis to assess the association between cytokine levels and the RelA gene expression. Higher levels of IFN-α and detectable levels of IL-10 protein and RelA mRNA were observed in surviving pigs compared to healthy (non-infected). A positive correlation of IFN-α cytokine levels with RelA mRNA expression was also obtained. In conclusion, 7 polymorphic events in the coding region of the RelA gene may contribute to the tolerance of ASFV in pigs.

Predisposing deleterious variants in the cancer-associated human kinases in the global populations

Human kinases play essential and diverse roles in the cellular activities of maintaining homeostasis and growth. Genetic mutations in the genes encoding the kinases (or phosphotransferases) have been linked with various types of cancers. In this study, we cataloged mutations in 500 kinases genes in >65,000 individuals of global populations from the Human Genetic Diversity Project (HGDP) and ExAC databases, and assessed their potentially deleterious impact by using the in silico tools SIFT, Polyphen2, and CADD. The analysis highlighted 35 deleterious non-synonymous SNVs in the ExAC and 5 SNVs in the HGDP project. Notably, a higher number of deleterious mutations was observed in the Non-Finnish Europeans (26 SNVs), followed by the Africans (14 SNVs), East Asians (13 SNVs), and South Asians (12 SNVs). The gene set enrichment analysis highlighted NTRK1 and FGFR3 being most significantly enriched among the kinases. The gene expression analysis revealed over-expression of NTRK1 in liver cancer, whereas, FGFR3 was found over-expressed in lung, breast, and liver cancers compared to their expression in the respective normal tissues. Also, 13 potential drugs were identified that target the NTRK1 protein, whereas 6 potential drugs for the FGFR3 target were identified. Taken together, the study provides a framework for exploring the predisposing germline mutations in kinases to suggest the underlying pathogenic mechanisms in cancers. The potential drugs are also suggested for personalized cancer management.

Investigating the Influence of ANTXR2 Gene Mutations on Protective Antigen Binding for Heightened Anthrax Resistance

Bacillus anthracis is the bacterium responsible for causing the zoonotic disease called anthrax. The disease presents itself in different forms like gastrointestinal, inhalation, and cutaneous. Bacterial spores are tremendously adaptable, can persist for extended periods and occasionally endanger human health. The Anthrax Toxin Receptor-2 (ANTXR2) gene acts as membrane receptor and facilitates the entry of the anthrax toxin into host cells. Additionally, mutations in the ANTXR2 gene have been linked to various autoimmune diseases, including Hyaline Fibromatosis Syndrome (HFS), Ankylosing Spondylitis (AS), Juvenile Hyaline Fibromatosis (JHF), and Infantile Systemic Hyalinosis (ISH). This study delves into the genetic landscape of ANTXR2, aiming to comprehend its associations with diverse disorders, elucidate the impacts of its mutations, and pinpoint minimal non-pathogenic mutations capable of reducing the binding affinity of the ANTXR2 gene with the protective antigen. Recognizing the pivotal role of single-nucleotide polymorphisms (SNPs) in shaping genetic diversity, we conducted computational analyses to discern highly deleterious and tolerated non-synonymous SNPs (nsSNPs) in the ANTXR2 gene. The Mutpred2 server determined that the Arg465Trp alteration in the ANTXR2 gene leads to altered DNA binding (p = 0.22) with a probability of a deleterious mutation of 0.808; notably, among the identified deleterious SNPs, rs368288611 (Arg465Trp) stands out due to its significant impact on altering the DNA-binding ability of ANTXR2. We propose these SNPs as potential candidates for hypertension linked to the ANTXR2 gene, which is implicated in blood pressure regulation. Noteworthy among the tolerated substitutions is rs200536829 (Ala33Ser), recognized as less pathogenic; this highlights its potential as a valuable biomarker, potentially reducing side effects on the host while also reducing binding with the protective antigen protein. Investigating these SNPs holds the potential to correlate with several autoimmune disorders and mitigate the impact of anthrax disease in humans.

Analysis of nonsynonymous SNPs in candidate genes that influence bovine temperament and evaluation of their effect in Brahman cattle

BACKGROUND

Temperament is an important production trait in cattle and multiple strategies had been developed to generate molecular markers to assist animal selection. As nonsynonymous single nucleotide polymorphisms are markers with the potential to affect gene functions, they could be useful to predict phenotypic effects. Genetic selection of less stress-responsive, temperamental animals is desirable from an economic and welfare point of view.

METHODS AND RESULTS

Two nonsynonymous single nucleotide polymorphisms identified in HTR1B and SLC18A2 candidate genes for temperament were analyzed in silico to determine their effects on protein structure. Those nsSNPs allowing changes in proteins were selected for a temperament association analysis in a Brahman population. Transversion effects on protein structure were evaluated in silico for each amino acid change model, revealing structural changes in the proteins of the HTR1B and SLC18A2 genes. The selected nsSNPs were genotyped in a Brahman population (n = 138), and their genotypic effects on three temperament traits were analyzed: exit velocity, pen score, and temperament score. Only the SNP rs209984404-HTR1B (C/A) showed a significant association (P = 0.0144) with pen score. The heterozygous genotype showed a pen score value 1.17 points lower than that of the homozygous CC genotype.

CONCLUSION

The results showed that in silico analysis could direct the selection of nsSNPs with the potential to change the protein. Non-synonymous single nucleotide polymorphisms causing structural changes and reduced protein stability were identified. Only rs209984404-HTR1B shows that the allele affecting protein stability was associated with the genotype linked to docility in cattle.

Association of Four Interleukin-8 Polymorphisms (-251 A>T, +781 C>T, +1633 C>T, +2767 A>T) with Ovarian Cancer Risk: Focus on Menopausal Status and Endometriosis-Related Subtypes

Interleukin-8 (IL-8) is involved in the regulation of inflammatory processes and carcinogenesis. Single-nucleotide polymorphisms (SNPs) within the IL-8 gene have been shown to alter the risks of lung, gastric, or hepatocellular carcinomas. To date, only one study examined the role of IL-8 SNPs in ovarian cancer (OC), suggesting an association between two IL-8 SNPs and OC risk. In this study, we investigated four common IL-8 SNPs, rs4073 (-251 A>T), rs2227306 (+781 C>T), rs2227543 (+1633 C>T), and rs1126647 (+2767 A>T), using the restriction fragment length polymorphism (PCR-RFLP) technique. Our study included a cohort of 413 women of Central European descent, consisting of 200 OC patients and 213 healthy controls. The most common (73.5%) histological type was high-grade serous OC (HGSOC), whereas 28/200 (14%) patients had endometriosis-related (clear cell or endometrioid) OC subtypes (EROC). In postmenopausal women, three of the four investigated SNPs, rs4073 (-251 A>T), rs2227306 (+781 C>T), and rs2227543 (+1633 C>T), were associated with OC risk. Furthermore, we are the first to report a significant relationship between the T allele or TT genotype of SNP rs1126647 (+2767 A>T) and the EROC subtype (p = 0.02 in the co-dominant model). The TT homozygotes were found more than twice as often in EROC compared to other OC subtypes (39% vs. 19%, p = 0.015). None of the examined SNPs appeared to influence OC risk in premenopausal women, nor were they associated with the aggressive HGSOC subtype or the stage of disease at the initial diagnosis.

In silico screening of non-synonymous SNPs in human TUFT1 gene

BACKGROUND

Tuftelin 1 (TUFT1) gene is important in the development and mineralization of dental enamel. The study aimed to identify potential functionally deleterious non-synonymous SNPs (nsSNPs) in the TUFT1 gene by using different in silico tools. The deleterious missense SNPs were identified from SIFT, PolyPhen-2, PROVEAN, SNPs & GO, PANTHER, and SNAP2. The stabilization, conservation, and three-dimensional modeling of mutant proteins were analyzed by I-Mutant 3.0, Consurf, and Project HOPE, respectively. The protein-protein interaction using STRING, GeneMANIA for gene-gene interaction, and DynaMut for evaluating the impact of the mutation on protein stability, conformation, and flexibility.

RESULTS

Eight deleterious nsSNPs (E242A, R303W, K182N, K123N, R117W, H289Q, R203W, and Q107R) out of 304 were found to have high-risk damaging effects using six in silico tools. Among them, K182N and K123N alone had increased stability, whereas E242A, R303W, R117W, H289Q, Q107R, and R203W exhibited a decrease in protein stability, based on DDG values. Meanwhile, all the eight deleterious nsSNPs altered the size, charge, hydrophobicity, and spatial organization of the amino acids and predominantly had alpha helix domains. These deleterious variants were located in highly conserved regions except R203W. Protein-protein interaction predicted that TUFT1 interacted with ten proteins that are involved in enamel mineralization and odontogenesis. Gene-gene interaction network showed that TUFT1 is involved in physical interactions, gene co-localization, and pathway interactions. DynaMut ΔΔG values predicted that five nsSNPs were destabilizing the protein, ΔΔG ENCoM values showed a destabilizing effect for all mutants, and seven nsSNPs increased the molecular flexibility of TUFT1.

CONCLUSION

Our study predicted eight functional SNPs that had detrimental effects on the structure and function of the TUFT1 gene. This will aid in the development of candidate deleterious markers as a potential target for disease diagnosis and therapeutic interventions.

DNA Polymorphisms and mRNA Levels of Immune Biomarkers as Candidates for Inflammatory Postpartum Disorders Susceptibility in Italian Buffaloes

The immunological genes that may interact with inflammatory postpartum diseases in Italian buffaloes were examined in this study. A total number of 120 female Italian buffaloes (60 normal and 60 with inflammatory reproductive diseases) were employed. Each buffalo's jugular vein was pierced to get five milliliters of blood. To obtain whole blood and extract DNA and RNA, the blood was placed within tubes containing sodium fluoride or EDTA anticoagulants. The immunological (IKBKG, LGALS, IL1B, CCL2, RANTES, MASP2, HMGB1, and S-LZ) genes' nucleotide sequence differences between healthy buffaloes and buffaloes affected by inflammatory reproductive diseases were found by employing PCR-DNA sequencing. According to Fisher's exact test (p ˂ 0.01), there were noticeably different probabilities of all major nucleotide changes spreading among buffalo groups with and without reproductive problems. Buffaloes were significantly more likely to express the examined genes when they had inflammatory reproductive diseases. The outcomes might support the significance of these markers' nucleotide variations and gene expression patterns as indicators of the prevalence of inflammatory reproductive disorders and provide a workable buffalo management policy.

In Silico Evaluation of the Potential Association of the Pathogenic Mutations of Alpha Synuclein Protein with Induction of Synucleinopathies

Alpha synuclein (α-Syn) is a neuronal protein encoded by the SNCA gene and is involved in the development of Parkinson's disease (PD). The objective of this study was to examine in silico the functional implications of non-synonymous single nucleotide polymorphisms (nsSNPs) in the SNCA gene. We used a range of computational algorithms such as sequence conservation, structural analysis, physicochemical properties, and machine learning. The sequence of the SNCA gene was analyzed, resulting in the mapping of 42,272 SNPs that are classified into different functional categories. A total of 177 nsSNPs were identified within the coding region; there were 20 variants that may influence the α-Syn protein structure and function. This identification was made by employing different analytical tools including SIFT, PolyPhen2, Mut-pred, SNAP2, PANTHER, PhD-SNP, SNP&Go, MUpro, Cosurf, I-Mut, and HOPE. Three mutations, V82A, K80E, and E46K, were selected for further examinations due to their spatial positioning within the α-Syn as determined by PyMol. Results indicated that these mutations may affect the stability and function of α-Syn. Then, a molecular dynamics simulation was conducted for the SNCA wildtype and the four mutant variants (p.A18G, p.V82A, p.K80E, and p.E46K). The simulation examined temperature, pressure, density, root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), solvent-accessible surface area (SASA), and radius of gyration (Rg). The data indicate that the mutations p.V82A, p.K80E, and p.E46K reduce the stability and functionality of α-Syn. These findings highlight the importance of understanding the impact of nsSNPs on α-syn structure and function. Our results required verifications in further protein functional and case-control studies. After being verified these findings can be used in genetic testing for the early diagnosis of PD, the evaluation of the risk factors, and therapeutic approaches.

Genetic association of Nuclear factor of activated T cells' 3'UTR and structural polymorphisms with susceptibility to generalized vitiligo in Gujarat population

Generalized vitiligo(GV) is an autoimmune skin depigmenting disease characterized by loss of functional melanocytes. Nuclear factor of activated T cells(NFATs) play a key role in regulatory T cells' (Tregs) activation and function. Our previous studies have highlighted the role of reduced NFATs expression and activity in impaired Tregs suppressive capacity, leading to GV pathogenesis. 3'UTR region and structural single nucleotide polymorphisms(SNPs) could lead to reduced NFAT expression and activity. Therefore, we studied the association of NFATs 3'UTR [NFATC2 rs4811198(T>G) &NFATC4 rs11848279(A>G)] and structural [NFATC1 rs754093(T>G) &NFATC2 rs12479626(T>C)] SNPs in 427 GV patients and 415 controls from Gujarat population by Polymerase chain reaction-restriction fragment length polymorphism(PCR-RFLP). Additionally, we carried out genotype-phenotype correlation and in silico analysis to assess the effect of NFATs SNPs on NFATs expression and structure. NFATC2 rs4811198(T>G) 3' UTR &NFATC2 rs12479626(T>C) structural SNPs were significantly associated with GV(p<0.0001). Interestingly, for NFATC2 rs4811198(T>G) SNP, there was a significant difference in the TT vs GG genotypes' frequencies (p=0.0034; Table 2), and for NFATC2 rs12479626(T>C) SNP there was a significant difference between TT vs TC and CC genotypes' frequencies(p<0.0001 & p=0.0002) between GV patients and controls. Furthermore, Odds ratio suggested that the susceptible alleles for NFATC2 rs4811198(T>G) &NFATC2 rs12479626(T>C) SNPs increased the risk of GV by 1.38 & 3.043 fold. However, the NFAT 3' UTR [NFATC2 rs4811198 (T>G)] and structural [NFATC1 rs754093(T>G)] SNPs were not significantly associated with GV. Interestingly, the genotype-phenotype correlation suggested that the susceptible 'G' allele of NFATC2 rs4811198(T>G) &NFATC4 rs11848279(A>G) 3' UTR SNPs lead to reduced NFATC2 and NFATC4 expression(p<0.0001). Furthermore, in silico analysis suggested that hsa-miR-3183 & hsa-miR-6720-3p miRNAs specifically bound to 'G' allele of NFATC2 rs4811198 SNP and has-miR-4652-3p miRNA specifically bound to 'G' allele of NFATC4 rs11848279 SNP. Overall, our study suggests that NFATC2 rs4811198(T>G) 3' UTR &NFATC2 rs12479626(T>C) structural SNPs may be associated with GV susceptibility in Gujarat population. Moreover, the susceptible alleles for the 3' UTR SNPs could lead to reduced NFATs levels, which may further possibly, affect the Treg suppressive function leading to GV.

Non-synonymous SNPs variants of PRKCG and its association with oncogenes predispose to hepatocellular carcinoma

BACKGROUND

PRKCG encodes PKC γ, which is categorized under the classical protein kinase C family. No studies have specifically established the relationship between PRKCG nsSNPs with structural and functional variations in PKC γ in the context of hepatocellular carcinoma (HCC). The present study aims to uncover this link through in-silico and experimental studies.

METHODS

The 3D structure of PKC γ was predicted. Molecular Dynamic (MD) Simulations were run and estimates were made for interactions, stability, conservation and post-translational alterations between wild and mutant structures. The association of PRKCG levels with HCC survival rate was determined. Genotyping analyses were conducted to investigate the deleterious PRKCG nsSNP association with HCC. mRNA expression of PKC γ, HIF-1 alpha, AKT, SOCS3 and VEGF in the blood of controls and HCC patients was analyzed and a genetic cascade was constructed depicting these interactions.

RESULTS

The expression level of studied oncogenes was compared to tumour suppressor genes. Through Alphafold, the 3D structure of PKC γ was explored. Fifteen SNPs were narrowed down for in-silico analyses that were identified in exons 5, 10 and 18 and the regulatory and kinase domain of PKC γ. Root mean square deviation and fluctuation along with the radius of gyration unveiled potential changes between the wild and mutated variant structures. Mutant genotype AA (homozygous) corresponding to nsSNP, rs386134171 had more frequency in patients with OR (2.446), RR (1.564) and P-values (< 0.0029) that highlights its significant association with HCC compared to controls in which the wild genotype GG was found more prevalent.

CONCLUSION

nsSNP rs386134171 can be a genetic marker for HCC diagnosis and therapeutic studies. This study has laid down a road map for future studies to be conducted on HCC.

Impact of highly deleterious non-synonymous polymorphisms on GRIN2A protein's structure and function

GRIN2A is a gene that encodes NMDA receptors found in the central nervous system and plays a pivotal role in excitatory synaptic transmission, plasticity and excitotoxicity in the mammalian central nervous system. Changes in this gene have been associated with a spectrum of neurodevelopmental disorders such as epilepsy. Previous studies on GRIN2A suggest that non-synonymous single nucleotide polymorphisms (nsSNPs) can alter the protein's structure and function. To gain a better understanding of the impact of potentially deleterious variants of GRIN2A, a range of bioinformatics tools were employed in this study. Out of 1320 nsSNPs retrieved from the NCBI database, initially 16 were predicted as deleterious by 9 tools. Further assessment of their domain association, conservation profile, homology models, interatomic interaction, and Molecular Dynamic Simulation revealed that the variant I463S is likely to be the most deleterious for the structure and function of the protein. Despite the limitations of computational algorithms, our analyses have provided insights that can be a valuable resource for further in vitro and in vivo research on GRIN2A-associated diseases.

Evaluation of SNPs from human IGFBP6 associated with gene expression: an in-silico study

In the bloodstream and in local tissues, most IGF molecules are associated with the one of the members of the IGF-binding protein (IGFBP) family, which are divided into six distinct types. IGF-binding proteins have been demonstrated to either decrease or increase the growth-promoting effects of IGFs on cell culture, by extending their half-life. They alter how IGFs interact with the receptors on their cell surfaces. IGFBP6 gene is associated with disease in-situ carcinoma. Upregulation or downregulation of IGFBP6 gene has been implicated in different types of cancer in humans. Nonsynonymous SNPs changes have the potential to affect the protein's structure and function. Potential functional SNPs can be assessed before undertaking studies in larger populations because validation of these functional SNPs can be a crucial problem. So, in this in-silico investigation, we used a variety of sequence- and structure-based bioinformatics methods to separate the potential nsSNPs of the IGFBP6 gene from the neutral ones. In total of 216 nsSNPs, 5 were found to have potential effects using 5 prediction tools. From which, 2 nsSNPs (R128G and R164H) were selected as potentially damaging due to their presence in highly conserved region and ability to decrease protein stability. Among these 2 nsSNPs, only R164H was found to be associated with Uterine corpus endometrial carcinoma. It was also found that both, upregulation or downregulation of IGFBP6 gene can lead to the different types of cancers. The findings of the present study will certainly be valuable in the future large population-based investigations as well as drug discovery, especially developing personalized medicine.Communicated by Ramaswamy H. Sarma.

High risk genetic variants of human insulin receptor substrate 1(IRS1) infer structural instability and functional interference

Insulin receptor substrate 1(IRS1) is a signaling adapter protein encoded by the IRS1 gene. This protein delivers signals from insulin and insulin-like growth factor-1(IGF-1) receptors to the phosphatidylinositol 3-kinases (P13K)/protein kinase B (Akt) and Extracellular signal-regulated kinases (Erk) - Mitogen-activated protein (MAP) kinase pathways, which regulate particular cellular processes. Mutations in this gene have been linked to type 2 diabetes mellitus, a heightened risk of insulin resistance, and an increased likelihood of developing a number of different malignancies. The structure and function of IRS1 could be severely compromised as a result of single nucleotide polymorphism (SNP) type genetic variants. In this study, we focused on identification of the most harmful non-synonymous SNPs (nsSNPs) of the IRS1 gene as well as prediction of their structural and functional consequences. Six different algorithms made the initial prediction that 59 of the 1142 IRS1 nsSNPs would have a negative impact on the protein structure. In-depth evaluations detected 26 nsSNPs located inside the functional domains of IRS1. Following that, 16 nsSNPs were identified as more harmful based on conservation profile, hydrophobic interaction, surface accessibility, homology modelling, and inter-atomic interactions. Following an in-depth analysis of protein stability, M249T (rs373826433), I223T (rs1939785175) and V204G (rs1574667052) were identified as three most deleterious SNPs and were subjected to molecular dynamics simulation for further insights. These findings will help us understand the implications for disease susceptibility, cancer progression, and the efficacy of therapeutic development against IRS1 gene mutants.Communicated by Ramaswamy H. Sarma.

In silico scanning of structural and functional deleterious nsSNPs in Arabidopsis thaliana's SOG1 protein, using molecular dynamic simulation approaches

Suppressor of gamma response 1 (SOG1) is a member of the NAC domain family transcription factors of the DNA damage response (DDR) signaling in the plant's genome. SOG1 is directly involved in transcriptional response to DNA damage, cell cycle checkpoints and ATR or ATM-mediated activation of the DNA damage responses and repair functioning in programmed cell death and regulation of end reduplication. Different mutations in the SOG1 protein lead to severe diseases and, ultimately, cell death. Single nucleotide polymorphisms (SNPs) are an important type of genetic alteration that cause different diseases or programmed cell death. The current study applied different computational approaches to Arabidopsis thaliana L. SOG1 protein to identify the potential deleterious nsSNPs and monitor their impact on the structure, function and protein stability. Various bioinformatics tools were applied to analyze the retrieved 34 nsSNPs and interestingly extracted four deleterious nsSNPs, that is, ensvath13968004 (Q166L), tmp18998388 (P159L), ensvath01103049 (K199N) and tmp18998295 (Y190F). For example, homology modeling, conservation and conformational analysis of the mutant's models were considered to scrutinize the deviations of these variants from the native SOG1 structure. All atoms molecular dynamic simulation confirmed the significance of these mutations on the protein stability, residual and structural conformation, compactness, surface conformation, dominant motion, Gibbs free energy distribution and dynamic effects. Similarly, protein-protein interaction revealed that SOG1 operates as a hub-linking cluster of various proteins, and any changes in the SOG1 might result in the disassociation of several signal transduction cascades.Communicated by Ramaswamy H. Sarma.

Structural Consequences of IRS-2 nsSNPs and Implication for Insulin Receptor Substrate-2 Protein Stability

Single-Nucleotide Polymorphisms (SNPs) are common genetic variations implicated in human diseases. The non-synonymous SNPs (nsSNPs) affect the proteins' structures and their molecular interactions with other interacting proteins during the accomplishment of biochemical processes. This ultimately causes proteins functional perturbation and disease phenotypes. The Insulin receptor substrate-2 (IRS-2) protein promotes glucose absorption and participates in the biological regulation of glucose metabolism and energy production. Several IRS-2 SNPs are reported in association with type 2 diabetes and obesity in human populations. However, there are no comprehensive reports about the protein structural consequences of these nsSNPs. Keeping in view the pathophysiological consequences of the IRS-2 nsSNPs, we designed the current study to understand their possible structural impact on coding protein. The prioritized list of the deleterious IRS-2 nsSNPs was acquired from multiple bioinformatics resources, including VEP (SIFT, PolyPhen, and Condel), PROVEAN, SNPs&GO, PMut, and SNAP2. The protein structure stability assessment of these nsSNPs was performed by MuPro and I-Mutant-3.0 servers via structural modeling approaches. The atomic-level structural and molecular dynamics (MD) impact of these nsSNPs were examined using GROMACS 2019.2 software package. The analyses initially predicted 8 high-risk nsSNPs located in the highly conserved regions of IRS-2. The MD simulation analysis eventually prioritized the N232Y, R218C, and R104H nsSNPs that predicted to significantly compromise the structure stability and may affect the biological function of IRS-2. These nsSNPs are predicted as high-risk candidates for diabetes and obesity. The validation of protein structural impact of these shortlisted nsSNPs may provide biochemical insight into the IRS-2-mediated type-2 diabetes.

Predicting the effects of rare genetic variants on oncogenic signaling pathways: A computational analysis of HRAS protein function

The HRAS gene plays a crucial role in regulating essential cellular processes for life, and this gene's misregulation is linked to the development of various types of cancers. Nonsynonymous single nucleotide polymorphisms (nsSNPs) within the coding region of HRAS can cause detrimental mutations that disrupt wild-type protein function. In the current investigation, we have employed in-silico methodologies to anticipate the consequences of infrequent genetic variations on the functional properties of the HRAS protein. We have discovered a total of 50 nsSNPs, of which 23 were located in the exon region of the HRAS gene and denoting that they were expected to cause harm or be deleterious. Out of these 23, 10 nsSNPs ([G60V], [G60D], [R123P], [D38H], [I46T], [G115R], [R123G], [P11OL], [A59L], and [G13R]) were identified as having the most delterious effect based on results of SIFT analysis and PolyPhen2 scores ranging from 0.53 to 69. The DDG values -3.21 kcal/mol to 0.87 kcal/mol represent the free energy change associated with protein stability upon mutation. Interestingly, we identified that the three mutations (Y4C, T58I, and Y12E) were found to improve the structural stability of the protein. We performed molecular dynamics (MD) simulations to investigate the structural and dynamic effects of HRAS mutations. Our results showed that the stable model of HRAS had a significantly lower energy value of -18756 kj/mol compared to the initial model of -108915 kj/mol. The RMSD value for the wild-type complex was 4.40 Å, and the binding energies for the G60V, G60D, and D38H mutants were -107.09 kcal/mol, -109.42 kcal/mol, and -107.18 kcal/mol, respectively as compared to wild-type HRAS protein had -105.85 kcal/mol. The result of our investigation presents convincing corroboration for the potential functional significance of nsSNPs in augmenting HRAS expression and adding to the activation of malignant oncogenic signalling pathways.

Association of IL-8-251T/A (rs4073) gene polymorphism with Systemic lupus erythematosus in a cohort of Egyptian patients

OBJECTIVE

Systemic lupus erythematosus (SLE) susceptibility was found to be correlated with genetic polymorphisms. Therefore, the current study aimed to investigate the association of the IL-8-251 T/A polymorphism with the risk of SLE.

METHODS

A total of 135 SLE patients and 75 controls were enrolled in our study. The IL-8-251 T/A polymorphism was analysed by PCR-RFLP. Also, the serum concentration of IL-8 was measured using ELISA. Finally, investigate possible IL-8 pathways in SLE pathogenesis by using the STRING database.

RESULTS

Our results revealed that the risk of having SLE in AA genotype carriers was significantly increased (OR = 1.92, 95 % CI = 1.23-3.10, p = 0.006) when compared with TT genotype carriers. Patients with SLE had a significantly higher frequency of the A allele (OR = 1.37, 95 % CI = 1.09-1.7; P = 0.01) than controls. Serum IL-8 levels were significantly increased in SLE patients (77.81 ± 21.27; p < 0.001) when compared to healthy controls (48.85 ± 7.89). Also, it was found that the serum IL-8 level had significant positive correlations with proteinuria, ESR, ANA, urea and SELADI, and significant negative correlations with RBCs count, C3 and hemoglobin. According to ROC curve analysis, serum IL-8 levels are a good biomarker for the detection of SLE disease, with 87.5 % sensitivity and 85 % specificity. STRING analysis revealed that IL-8 interacts with key SLE signaling pathway members such as TNF-α, IL-1β, IL-6, and IL-10.

CONCLUSIONS

There was a correlation between the IL-8-251 T/A polymorphism and the risk of SLE. Our findings also suggest that the IL8-251 A allele may be an important risk factor for the development of SLE.

Analysis of amplification and association polymorphisms in the bovine beta-defensin 129 (BBD129) gene revealed its function in bull fertility

β-defensins are adsorbable on the sperm surface in the male reproductive tract (MRT) and enhance sperm functional characteristics. The beta-defensin 129 (DEFB129) antimicrobial peptide is involved in sperm maturation, motility, and fertilization. However, its role in bovine fertility has not been well investigated. This study examines the relationship between the bovine BBD129 gene and Bos indicus x Bos taurus bull fertility. The complete coding sequence of BBD129 mRNA was identified by RNA Ligase Mediated-Rapid Amplification of cDNA End (RLM-RACE) and Sanger sequencing methodologies. It consisted of 582 nucleotides (nts) including 5' untranslated region (UTR) (46nts) and 3'UTR (23nts). It conserves all beta-defensin-like features. The expression level of BBD129 was checked by RT-qPCR and maximal expression was detected in the corpus-epididymis region compared to other parts of MRT. Polymorphism in BBD129 was also confirmed by Sanger sequencing of 254 clones from 5 high fertile (HF) and 6 low fertile (LF) bulls at two positions, 169 T > G and 329A > G, which change the S57A and N110S in the protein sequence respectively. These two mutations give rise to four types of BBD129 haplotypes. The non-mutated TA-BBD129 (169 T/329A) haplotype was substantially more prevalent among high-fertile bulls (P < 0.005), while the double-site mutated GG-BBD129 (169 T > G/329A > G) haplotype was significantly more prevalent among low-fertile bulls (P < 0.005). The in silico analysis confirmed that the polymorphism in BBD129 results in changes in mRNA secondary structure, protein conformations, protein stability, extracellular-surface availability, post-translational modifications (O-glycosylation and phosphorylation), and affects antibacterial and immunomodulatory capabilities. In conclusion, the mRNA expression of BBD129 in the MRT indicates its region-specific dynamics in sperm maturation. BBD129 polymorphisms were identified as the deciding elements accountable for the changed proteins with impaired functionality, contributing to cross-bred bulls' poor fertility.

Genetic Factors for Coronary Heart Disease and Their Mechanisms: A Meta-Analysis and Comprehensive Review of Common Variants from Genome-Wide Association Studies

Genome-wide association studies (GWAS) have discovered 163 loci related to coronary heart disease (CHD). Most GWAS have emphasized pathways related to single-nucleotide polymorphisms (SNPs) that reached genome-wide significance in their reports, while identification of CHD pathways based on the combination of all published GWAS involving various ethnicities has yet to be performed. We conducted a systematic search for articles with comprehensive GWAS data in the GWAS Catalog and PubMed, followed by a meta-analysis of the top recurring SNPs from ≥2 different articles using random or fixed-effect models according to Cochran Q and I² statistics, and pathway enrichment analysis. Meta-analyses showed significance for 265 of 309 recurring SNPs. Enrichment analysis returned 107 significant pathways, including lipoprotein and lipid metabolisms (rs7412, rs6511720, rs11591147, rs1412444, rs11172113, rs11057830, rs4299376), atherogenesis (rs7500448, rs6504218, rs3918226, rs7623687), shared cardiovascular pathways (rs72689147, rs1800449, rs7568458), diabetes-related pathways (rs200787930, rs12146487, rs6129767), hepatitis C virus infection/hepatocellular carcinoma (rs73045269/rs8108632, rs56062135, rs188378669, rs4845625, rs11838776), and miR-29b-3p pathways (rs116843064, rs11617955, rs146092501, rs11838776, rs73045269/rs8108632). In this meta-analysis, the identification of various genetic factors and their associated pathways associated with CHD denotes the complexity of the disease. This provides an opportunity for the future development of novel CHD genetic risk scores relevant to personalized and precision medicine.

A novel germline HAVCR2 (TIM-3) compound heterozygous mutation is related to hemophagocytic lymphohistiocytic syndrome in EBV-positive peripheral T-cell lymphoma (NOS) with down-regulated TIM-3 signaling

Recently, it have been reported that Hepatitis A Virus-Cellular Receptor 2(HAVCR2,encoding T-cell immunoglobulin and Mucin-Containing Protein 3[TIM3]) mutations are associated with severe hemophagocytic syndrome(HLH) in subcutaneous panniculitis-like T-cell lymphoma(SPTCL),and there are also frequent mutations in sporadic SPTCL, suggesting the individuals harboring HAVCR2(TIM-3) germline mutations are highly susceptible to familial or sporadic SPTCL. Here, we identify a novel germline compound heterozygous mutation of TIM-3 gene,c.245A>G (p.Tyr82Cys) and c.265C>T(p.Arg89Cys) variations in a single familial case with EBV-positive peripheral T-cell lymphoma(NOS),accompanied HLH;we also detected Tyr82Cys germline mutation in TIM-3 gene in one sporadic patient with cutaneous T cell lymphoma. We screened the distributive frequencies for TIM-3 mutations in healthy controls(n=87), B-(n=79) or T-cell lymphoma(n=25) not SPTCL, and the results showed that the mutation was found in two out of 25 patients with T-cell lymphoma but was not detected in 79 patients with B-cell lymphoma nor in a group of 87 controls. The mRNA expression of TIM-3 on primary cells and transfected HEK293 cells reduced significantly, indicating Tyr82Cys and Arg89Cys mutations is a loss-of function mutations on TIM-3,resulting in a weakened TIM-3 signaling. Our results suggest Tyr82Cys TIM-3 germline mutations are not only limited in SPTCL, and also occurred in other types of T-cell lymphoma, especially complicated HLH. TIM-3 mutations may be an predisposing factor for T-cell lymphoma and molecular marker for auxiliary diagnosis in T cell lymphoma,especially complicated with HLH.

In silico analyses of Wnt1 nsSNPs reveal structurally destabilizing variants, altered interactions with Frizzled receptors and its deregulation in tumorigenesis

Wnt1 is the first mammalian Wnt gene, which is discovered as proto-oncogene and in human the gene is located on the chromosome 12q13. Mutations in Wnt1 are reported to be associated with various cancers and other human diseases. The structural and functional consequences of most of the non-synonymous SNPs (nsSNPs), present in the human Wnt1 gene, are not known. In the present work, extensive bioinformatics analyses are used to screen 292 nsSNPs of Wnt1 for predicting pathogenic and harmless polymorphisms. We have identified 10 highly deleterious nsSNPs among which 7 are located within the highly conserved areas. These 10 nsSNPs are also predicted to affect the post-translational modifications of Wnt1. Further, structure based stability analyses of these 10 highly deleterious nsSNPs revealed 8 variants as highly destabilizing. These 8 highly destabilizing variants were shown to have high BC score and high RMSIP score from normal mode analyses. Based on the deformation energies, obtained from the normal mode analyses, variants like G169A, G169S, G331R and G331S were found to be unstable. Molecular Dynamics (MD) simulations revealed structural stability and fluctuation of WT Wnt1 and its prioritized variants. RMSD remained fluctuating mostly between 4 and 5 Å and occasionally between 3.5 and 5.5 Å ranges. RMSF in the CTD region (residues 330–360) of the binding pocket were lower compared to that of WT. Studying the impacts of nsSNPs on the binding interface of Wnt1 and seven Frizzled receptors have predicted substitutions which can stabilize or destabilize the binding interface. We have found that Wnt1 and FZD8-CRD is the best docked complex in our study. MD simulation based analyses of wild type Wnt1-FZD8-CRD complex and the 8 prioritized variants revealed that RMSF was higher in the unstructured regions and RMSD remained fluctuating in the region of 5 Å ± 1 Å. We have also observed differential Wnt1 gene expression pattern in normal, tumor and metastatic conditions across different tissues. Wnt1 gene expression was significantly higher in metastatic tissues of lungs, colon and skin; and was significantly lower in metastatic tissues of breast, esophagus and kidney. We have also found that Wnt1 deregulation is associated with survival outcome in patients with gastric and breast cancer. Furthermore, these computationally screened highly deleterious nsSNPs of Wnt1 can be analyzed in population based genetic studies and may help understand the Wnt1 associated diseases.

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.

Structural and Dynamic Investigation of non-synonymous variations in Renin-AGT complex revealed altered binding via hydrogen bonding network reprogramming to accelerate the hypertension pathway

Hypertension is one of the major issues worldwide and one of the main factors involved in heart and kidney failure. Angiotensinogen and renin are key components of the renin-angiotensin-aldosterone system (RAAS), which plays an indispensable role in hypertension. The aimed of this study to find out the non-synonymous mutations and structure-based mutation-function correlation in the Renin-AGT complex and reveal the most deleterious mutations to accelerated hypertension. In the current study, we employed computational modelling and molecular simulation approaches to demonstrate the impact of specific mutations in the REN-AGT interface in hypertension. Computational algorithms i.e. PhD-SNP, PolyPhen-1, MAPP, SIFT, SNAP, PredictSNP, PolyPhen-2, and PANTHER predicted 20 mutations as deleterious in AGT while only five mutations were conformed as deleterious in the Renin protein. Investigation of the bonding analysis revealed that two mutations S107L and V193F in Renin altered the hydrogen-bonding paradigm at the interface site. Furthermore, exploration of structural-dynamic behaviors demonstrated by that these mutations also increases the structural stability to regulate the expression of disease pathway. The flexibility index of each residues and structural compactness analysis further validated the findings by portraying the difference in the dynamic behavior in contrast to the wild type. Binding energy calculations based on molecular mechanics/generalized Born surface area (MM/GBSA) methods were used which further established the binding differences between the wild type, S107L, and V193F mutant variants. The total binding energy for wild type, S107L, and V193F were reported to be -27.79 kcal/mol, -47.72 kcal/mol, and -38.25 kcal/mol respectively. In conclusion, these two mutations increase the binding free energy alongside the docking score to enhance the binding between Renin and AGT to overexpress this pathway in a hypertension disease condition. Patients with these mutations may be screened for potential therapeutic intervention.

Genetic determinants of lung cancer: Understanding the oncogenic potential of somatic missense mutations

BACKGROUND

Treatment of lung cancer is getting more personalized nowadays and medical practitioners are moving away from conventional histology-driven empirical treatments, platinum-based chemotherapy, and other invasive surgical resections and have started adopting alternate therapies in which therapeutic targets are patient's molecular oncogenic drivers.

AIM

The aim of the current study is to extract meaningful information from the online somatic mutation data (retrieved from cBioPortal) of 16 most significantly mutated oncogenes in non-small-cell lung cancer (NSCLC), namely EGFR, NRAS, KRAS, HER2 (ERBB2), RET, MET, ROS1, FGFR1, BRAF, AKT1, MEK1 (MAP2K1), PIK3CA, PTEN, DDR2, LKB1 (STK11) and ALK, for improving our understanding of the pathobiology of the lung cancer that can aid decision-making on critical clinical and therapeutic considerations.

METHODS

Using an integrated approach comprising 4 steps, the oncogenic potential of 661 missense non-synonymous single nucleotide polymorphisms (nsSNPs) in 16 genes was ascertained using 2059 NSCLC (1575 lung adenocarcinomas, 484 lung squamous cell carcinomas) patients' online mutation data. The steps used comprise sequence/structure homology-based prediction, scoring of conservation of mutated residues and positions, prediction of resulting molecular and functional consequences using machine-learning and structure-guided approach.

RESULTS

Out of a total of 661 nsSNPs analyzed, a set of 29 nsSNPs has been identified as conserved high confidence mutations in 10 of 16 genes relevant to the under study. Out of 29 conserved high confidence nsSNPs, 4 nsSNPs (EGFR N1094Y, BRAF M620I, DDR2 R307L, ALK P1350T) have been found to be putative novel rare genetic markers for NSCLC.

CONCLUSIONS

The current study, the first of its kind, has provided a list of deleterious non-synonymous somatic mutations in a selected pool of oncogenes that can be considered as a promising target for future drug design and therapy for patients with lung adenocarcinomas and squamous cell carcinomas.

PSnpBind: a database of mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow

A key concept in drug design is how natural variants, especially the ones occurring in the binding site of drug targets, affect the inter-individual drug response and efficacy by altering binding affinity. These effects have been studied on very limited and small datasets while, ideally, a large dataset of binding affinity changes due to binding site single-nucleotide polymorphisms (SNPs) is needed for evaluation. However, to the best of our knowledge, such a dataset does not exist. Thus, a reference dataset of ligands binding affinities to proteins with all their reported binding sites' variants was constructed using a molecular docking approach. Having a large database of protein-ligand complexes covering a wide range of binding pocket mutations and a large small molecules' landscape is of great importance for several types of studies. For example, developing machine learning algorithms to predict protein-ligand affinity or a SNP effect on it requires an extensive amount of data. In this work, we present PSnpBind: A large database of 0.6 million mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow. It provides a web interface to explore and visualize the protein-ligand complexes and a REST API to programmatically access the different aspects of the database contents. PSnpBind is open source and freely available at https://psnpbind.org .

Discovery of DNA polymorphisms via whole genome resequencing and their functional relevance in salinity stress response in chickpea

Salinity stress is one of the major constraints for plant growth and yield. The salinity stress response of different genotypes of crop plants may largely be governed by DNA polymorphisms. To determine the molecular genetic factors involved in salinity stress tolerance in chickpea, we performed a whole genome resequencing data analysis of three each of salinity-sensitive and salinity-tolerant genotypes. A total of 6173 single nucleotide polymorphisms and 920 insertions and deletions differentiating the chickpea genotypes with contrasting salinity stress responses were identified. Gene ontology analysis revealed the enrichment of functional terms related to stress response and development among the genes harboring DNA polymorphisms in their promoter and/or coding regions. DNA polymorphisms located within the cis-regulatory motifs of the quantitative trait loci (QTL)-associated and abiotic stress related genes were identified, which may influence salinity stress response via modulating binding affinity of the transcription factors. Several genes including QTL-associated and abiotic stress response related genes harboring DNA polymorphisms exhibited differential expression in response to salinity stress especially at the reproductive stage of development in the salinity-tolerant genotype. Furthermore, effects of non-synonymous DNA polymorphisms on mutational sensitivity and structural integrity of the encoded proteins by the candidate QTL-associated and abiotic stress response related genes were revealed. The results suggest that DNA polymorphisms may determine salinity stress response via influencing differential gene expression in genotype and/or stage-dependent manner. Altogether, we provide a high-quality set of DNA polymorphisms and candidate genes that may govern salinity stress tolerance in chickpea.

Computational and Mass Spectrometry-Based Approach Identify Deleterious Non-Synonymous Single Nucleotide Polymorphisms (nsSNPs) in JMJD6

The jumonji domain-containing protein 6 (JMJD6) gene catalyzes the arginine demethylation and lysine hydroxylation of histone and a growing list of its known substrate molecules, including p53 and U2AF65, suggesting a possible role in mRNA splicing and transcription in cancer progression. Mass spectrometry-based technology offers the opportunity to detect SNP variants accurately and effectively. In our study, we conducted a combined computational and filtration workflow to predict the nonsynonymous single nucleotide polymorphisms (nsSNPs) present in JMJD6, followed by a liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and validation. The computational approaches SIFT, PolyPhen-2, SNAP, I-Mutant 2.0, PhD-SNP, PANTHER, and SNPS&GO were integrated to screen out the predicted damaging/deleterious nsSNPs. Through the three-dimensional structure of JMJD6, H187R (rs1159480887) was selected as a candidate for validation. The validation experiments showed that the mutation of this nsSNP in JMJD6 obviously affected mRNA splicing or the transcription of downstream genes through the reduced lysyl-hydroxylase activity of its substrates, U2AF65 and p53, further indicating the accuracy of this prediction method. This research provides an effective computational workflow for researchers with an opportunity to select prominent deleterious nsSNPs and, thus, remains promising for examining the dysfunction of proteins.

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein-protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein-protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.

In silico analysis of deleterious SNPs of human MTUS1 gene and their impacts on subsequent protein structure and function

The mitochondrial tumor suppressor 1 (MTUS1) gene acts as a crucial tumor suppressor by inhibiting growth and proliferation of eukaryotic cells including tumor cell lines. Down regulation of MTUS1 gene has been implicated in a wide range of cancers as well as various human diseases. Alteration through nsSNPs can potentially damage the structure and/or function of the protein. As characterization of functional SNPs in such disease linked genes is a major challenge, it is feasible to analyze putative functional SNPs prior to performing larger population studies. Hence, in this in silico study we differentiated the potentially harmful nsSNPs of the MTUS1 gene from the neutral ones by using various sequence and structure based bioinformatic tools. In a total of 215 nsSNPs, 9 were found to be most likely to exert deleterious effect using 7 prediction tools. From which, 5nsSNPs (S1259L, E960K, P503T, L1084V and L1143Q) were selected as potentially damaging due to their presence in the highly conserved region and ability to decrease protein stability. In fact, 2 nsSNPs (S1259L and E960K) among these 5 were found to be individually associated with two distinctive cancers named Stomach adenocarcinoma and Uterine corpus endometrial carcinoma. As this is the first comprehensive study analyzing the functional nsSNPs of MTUS1, the results of the current study would certainly be helpful in future prospects concerning large population-based studies as well as drug discovery, especially developing individualized medicine.

Common polymorphic OTC variants can act as genetic modifiers of enzymatic activity

Understanding the role of common polymorphisms in modulating the clinical phenotype when they co‐occur with a disease‐causing lesion is of critical importance in medical genetics. We explored the impact of apparently neutral common polymorphisms, using the gene encoding the urea cycle enzyme, ornithine transcarbamylase (OTC), as a model system. Distinct combinations of genetic backgrounds embracing two missense polymorphisms were created in cis with the pathogenic p.Arg40His replacement. In vitro enzymatic assays revealed that the polymorphic variants were able to modulate OTC activity both in the presence or absence of the pathogenic lesion. First, we found that the combination of the minor alleles of polymorphisms p.Lys46Arg and p.Gln270Arg significantly enhanced enzymatic activity in the wild‐type protein. Second, enzymatic assays revealed that the minor allele of the p.Gln270Arg polymorphism was capable of ameliorating OTC activity when combined in cis with the pathogenic p.Arg40His replacement. Structural analysis predicted that the minor allele of the p.Gln270Arg polymorphism would serve to stabilize the OTC wild‐type protein, thereby corroborating the results of the experimental assays. Our findings demonstrate the potential importance of cis‐interactions between common polymorphic variants and pathogenic missense mutations and illustrate how standing genetic variation can modulate protein function.

HLA-DQ Gene Polymorphisms Associated with Hepatitis B Virus in Patients with Hepatocellular Carcinoma Progression

Background:

The development and progression of hepatitis B (HBV)-related disease can lead to liver cirrhosis and hepatocellular carcinoma (LC and HCC, respectively). Human leukocyte antigen (HLA) DQ polymorphism has been reported in other recent studies to deal with the association between HBV and liver disease. Our study on the Egyptian population was introduced to assess the strong association between HLA-DQ polymorphism and HBV infection in addition to the progression of HCC.

Aim:

The aim of this work was to estimate HLA-DQ gene polymorphisms in HBV and HCC.

Methods:

HLA-DQ genotype polymorphism was assayed by using the ABI Taq Man allelic discrimination assay in different groups in this study. According to the relevant HLA Class II single nucleotide polymorphism (SNP) literature, one single nucleotide polymorphism (SNPs) was selected as the candidate site; it was an HLA-DQ, which showed minor allele frequencies AA, GA, and GG.

Results:

Haplotype analysis was performed on all the subjects in the study; AA haplotype was the most frequent haplotype in HCC cases (18%) in comparison with HBV and healthy individuals (3%). The haplotype GA was more frequent in the HCC group and slightly more frequent in LC when compared to HBV only cases and also when compared to the control group. In contrast, the GG haplotype was recorded less frequently in HCC individuals, but the HBV and LC groups showed more frequency of this haplotype compared with the HCC group. There was a correlation between AFP serum levels and the frequency of GA and AA polymorphism in HCC cases.

Conclusion:

We found that AA and GA haplotype was significantly most frequent in HCC. Our findings suggest that HLA-DQ AA and GG polymorphism might serve as a novel potential predictive marker for HCC and may function in tumorigenesis of HBV.

The Unique Biology behind the Early Onset of Breast Cancer

Breast cancer commonly affects women of older age; however, in developing countries, up to 20% of breast cancer cases present in young women (younger than 40 years as defined by oncology literature). Breast cancer in young women is often defined to be aggressive in nature, usually of high histological grade at the time of diagnosis and negative for endocrine receptors with poor overall survival rate. Several researchers have attributed this aggressive nature to a hidden unique biology. However, findings in this aspect remain controversial. Thus, in this article, we aimed to review published work addressing somatic mutations, chromosome copy number variants, single nucleotide polymorphisms, differential gene expression, microRNAs and gene methylation profile of early-onset breast cancer, as well as its altered pathways resulting from those aberrations. Distinct biology behind early-onset of breast cancer was clear among estrogen receptor-positive and sporadic cases. However, further research is needed to determine and validate specific novel markers, which may help in customizing therapy for this group of patients.

Epigenetic Regulatory Enzymes: mutation Prevalence and Coexistence in Cancers

Epigenetic regulation is an important layer of transcriptional control with the particularity to affect the broad spectrum of genome. Over the years, largely due to the substantial number of recurrent mutations, there have been hundreds of novel driver genes characterized in various cancers. Additionally, the relative contribution of two dysregulated epigenomic entities (DNA methylation and histone modifications) that gradually drive the cancer phenotype remains in the research focus. However, a complex scenario arises when the disease phenotype does not harbor any relevant mutation or an abnormal transcription level. Although the cancer landscape involves the contribution of multiple genetic and non-genetic factors, herein, we discuss specifically the mutation spectrum of epigenetically-related enzymes in cancer. In addition, we address the coexistence of these two epigenetic entities in malignant human diseases, especially cancer. We suggest that the study of epigenetically-related somatic mutations in the early cellular differentiation stage of embryonic development might help to understand their later-staged footprints in the cancer genome. Furthermore, understanding the co-occurrence and/or inverse association of different disease types and redefining the general definition of "healthy" controls could provide insights into the genome reorganization.

Transcriptome analysis and association mapping reveal the genetic regulatory network response to cadmium stress in Populus tomentosa

Long non-coding RNAs (lncRNAs) play essential roles in plant abiotic stress responses, but the response of lncRNA-mediated genetic networks to cadmium (Cd) treatment remain elusive in trees, the promising candidates for phytoremediation of Cd contamination. We identified 172 Cd-responsive lncRNAs and 295 differentially expressed target genes in the leaves of Cd-treated Populus tomentosa. Functional annotation revealed that these lncRNAs were involved in various processes, including photosynthesis, hormone regulation, and phenylalanine metabolism. Association studies identified 78 significant associations, representing 14 Cd-responsive lncRNAs and 28 target genes for photosynthetic and leaf physiological traits. Epistasis uncovered 83 pairwise interactions among these traits, revealing Cd-responsive lncRNA-mediated genetic networks for photosynthesis and leaf physiology in P. tomentosa. We focused on the roles of two Cd-responsive lncRNA-gene pairs, MSTRG.22608.1-PtoMYB73 and MSTRG.5634.1-PtoMYB27, in Cd tolerance of Populus, and detected insertions/deletions within lncRNAs as polymorphisms driving target gene expression. Genotype analysis of lncRNAs and heterologous overexpression of PtoMYB73 and PtoMYB27 in Arabidopsis indicated their effects on enhancing Cd tolerance, photosynthetic rate, and leaf growth, and the potential interaction mechanisms of PtoMYB73 with abiotic stresses. Our study identifies the genetic basis for the response of Populus to Cd treatment, facilitating genetic improvement of Cd tolerance in trees.

Functional Analysis of Single Nucleotide Polymorphism in ZUFSP Protein and Implication in Pathogenesis

Researches have revealed that functional non-synonymous Single Nucleotide Polymorphism (nsSNPs) present in the Zinc-finger with UFM1-Specific Peptidase domain protein (ZUFSP) may be involved in genetic instability and carcinogenesis. For the first time, we employed in-silico approach using predictive tools to identify and validate potential nsSNPs that could be pathogenic. Our result revealed that 8 nsSNPs (rs 112738382, rs 140094037, rs 201652589, rs 201847265, rs 202076827, rs 373634906, rs 375114528, rs 772591104) are pathogenic after being subjected to rigorous filtering process. The structural impact of the nsSNPs on ZUFSP structure indicated that the nsSNPs affect the stability of the protein by lowering ZUFSP protein stability. Furthermore, conservation analysis showed that rs 201652589, rs 140094037, rs 201847265, and rs 772591104 were highly conserved. Interestingly, the protein-protein affinity between ZUFSP and Ubiquitin was altered rs 201652589, rs 140094037, rs 201847265, and rs 772591104 had a binding affinity of - 0.46, - 0.83, - 1.62, and - 1.12 kcal/mol respectively. Our study has been able to identify potential nsSNPs that could be used as genetic biomarkers for some diseases arising as a result of aberration in the ZUFSP structure, however, being a predictive study, the identified nsSNPs need to be experimentally investigated.

Polymorphisms of the genes CTLA4, PTPN22, CD40, and PPARG and their roles in Graves' disease: susceptibility and clinical features

PURPOSE

CTLA4, PTPN22, and CD40 are immune-regulatory genes strongly associated with GD, as well as PPARG, but their clinical significance in different populations is still uncertain.

METHODS

We genotyped 282 Brazilian GD patients (234 women and 48 men, 39.80 ± 11.69 years old), including 144 patients with GO, and 308 healthy control individuals (246 women and 62 men, 36.86 ± 12.95 years old).

RESULTS

A multivariate analysis demonstrated that the inheritance of the GG genotype rs3087243 of CTLA4 (OR = 2.593; 95% CI = 1.630-4.123; p < 0.0001) and the CC genotype of rs3789607 of PTPN22 (OR = 2.668; 95% CI = 1.399-5.086; p = 0.0029) consisted in factors independent of the susceptibility to GD. The inheritance of polymorphic genotypes of rs5742909 of CTLA4 was associated with older age at the time of diagnosis (42.90 ± 10.83 versus 38.84 ± 11.81 years old; p = 0.0105), with higher TRAb levels (148.17 ± 188.90 U/L versus 112.14 ± 208.54 U/L; p = 0.0229) and the need for higher therapeutic doses of radioiodine (64.23 ± 17.16 versus 50.22 ± 16.86; p = 0.0237). The inheritance of the CC genotype of rs1883832 CD40 gene was more frequent among women (69.65%) than men (52.00%; p = 0.0186). The polymorphic genotype of PPARG gene (rs1801282) was associated with TPOAb positivity (p = 0.0391), and the GG genotype of rs2476601 of PTPN22 gene was associated with positivity for both TgAb (p = 0.0360) and TPOAb (p < 0.0001). Both polymorphic genotypes rs2476601 and rs3789607 of the PTPN22 gene were more frequent among nonsmoking patients (p = 0.0102 and p = 0.0124, respectively).

CONCLUSIONS

Our data confirm the important role of CTLA4 polymorphisms in GD susceptibility; demonstrate the role of PTPN22 polymorphisms in patients' clinical features; and suggest these genes may influence the severity of the disease.

Predicting the functional consequences of genetic variants in co-stimulatory ligand B7-1 using in-silico approaches

The purpose of this research is to identify and characterize deleterious genetic variants in the co-stimulatory ligand B7-1, also known as the human cluster of differentiation CD80 marker. The B7-1 ligand and the major histocompatibility complex class II (MHC II) molecules are the main determinants that provide B-cells the required competency to act as antigen presenting cells. For this, participation of both MHC class II molecules and CD80 is required. The interaction of the CD80 ligand with CD28 on the surface 7 of T_H cells plays a key role in the activation of TH cells and progression of B cells through the S phase, hence, leading to their proliferation in mitosis. A set of 2313 genetic variants in the B7-1 ligand have been mapped and retrieved from dbSNP database. Subsequently, 150 non-synonymous single nucleotide polymorphisms (nsSNPs) were mapped and subjected to the sequence and structural homology based predictions, which were further analyzed for protein stability and the disease phenotypes. Finally, we identified 7 potentially damaging nsSNPs in the B7-1 ligand that may affect its interaction with the cognitive receptor CD28, hence, may also interfere with TH cell activation and B cell proliferation. We propose that subsequent experimental analyses (stability, expression and interactions) on these proteins can provide a deep understanding about the effect of these variants on the structure and function of CD80.

Transcription-translation error: In-silico investigation of the structural and functional impact of deleterious single nucleotide polymorphisms in GULP1 gene

Nonsynonymous single nucleotide polymorphisms (nsSNPs) are one of the most common forms of mutations known to disrupt the product of translation thereby altering the protein structure-function relationship. GULP1 (PTB domain-containing engulfment adaptor protein 1) is an evolutionarily conserved adaptor protein that has been associated with glycated hemoglobin (HbA1c) in Genome-Wide Association Studies (GWAS). In order to understand the role of GULP1 in the etiology of diabetes, it is important to study some functional nsSNPs present within the GULP1 protein. We, therefore, used a SNPinformatics approach to retrieve, classify, and determine the stability effect of some nsSNPs. Y27C, G142D, A144T, and Y149C were jointly predicted by the pathogenic-classifying tools to be disease-causing, however, only G142D, A144T, and Y149C had their structural architecture perturbed as predicted by I-MUTANT and MuPro. Interestingly, G142D and Y149C occur at positions 142 and 149 of GULP1 which coincidentally are found within the binding site of GULP1. Protein-Protein interaction analysis also revealed that GULP1 interacted with 10 proteins such as Cell division cycle 5-like protein (CDC5L), ADP-ribosylation factor 6 (ARF6), Arf-GAP with coiled-coil (ACAP1), and Multiple epidermal growth factor-like domains protein 10 (MEGF10), etc. Taken together, rs1357922096, rs1264999716, and rs128246649 could be used as genetic biomarkers for the diagnosis of diabetes. However, being a computational study, these nsSNPs require experimental validation to explore their metabolic involvement in the pathogenesis of diseases.

Exploring the effect of nsSNPs in human YPEL3 gene in cellular senescence

YPEL3 that induces cellular senescence in both normal and tumour cells of humans may show altered expression under the influence of incidental mutations. In this study, we proposed the first structure of Native YPEL3 protein and its five possible deleterious mutants—V40M, C61Y, G98R, G108S, and A131T and predicted their deleterious effects to alter stability, flexibility and conformational changes in the protein. The MD simulation (RMSD, RMSF, Rg, h-bond and SASA) analysis revealed that the variants V40M, G98R and G108S increased the flexibility in protein, and variant V40M imparted more compactness to the protein.. In general, variants attributed changes in the native conformation and structure of the YPEL3 protein which might affect the native function of cellular senescence. The study provides opportunities for health professionals and practitioners in formulating précised medicines to effectively cure various cancers. We propose in-vitro or in-vivo studies should consider these reported nsSNPs while examining any malfunction in the YPEL3 protein.

Predicting the most deleterious missense nsSNPs of the protein isoforms of the human HLA-G gene and in silico evaluation of their structural and functional consequences

BACKGROUND

The Human Leukocyte Antigen G (HLA-G) protein is an immune tolerogenic molecule with 7 isoforms. The change of expression level and some polymorphisms of the HLA-G gene are involved in various pathologies. Therefore, this study aimed to predict the most deleterious missense non-synonymous single nucleotide polymorphisms (nsSNPs) in HLA-G isoforms via in silico analyses and to examine structural and functional effects of the predicted nsSNPs on HLA-G isoforms.

RESULTS

Out of 301 reported SNPs in dbSNP, 35 missense SNPs in isoform 1, 35 missense SNPs in isoform 5, 8 missense SNPs in all membrane-bound HLA-G isoforms and 8 missense SNPs in all soluble HLA-G isoforms were predicted as deleterious by all eight servers (SIFT, PROVEAN, PolyPhen-2, I-Mutant 3.0, SNPs&GO, PhD-SNP, SNAP2, and MUpro). The Structural and functional effects of the predicted nsSNPs on HLA-G isoforms were determined by MutPred2 and HOPE servers, respectively. Consurf analyses showed that the majority of the predicted nsSNPs occur in conserved sites. I-TASSER and Chimera were used for modeling of the predicted nsSNPs. rs182801644 and rs771111444 were related to creating functional patterns in 5'UTR. 5 SNPs in 3'UTR of the HLA-G gene were predicted to affect the miRNA target sites. Kaplan-Meier analysis showed the HLA-G deregulation can serve as a prognostic marker for some cancers.

CONCLUSIONS

The implementation of in silico SNP prioritization methods provides a great framework for the recognition of functional SNPs. The results obtained from the current study would be called laboratory investigations.

In Silico Evaluation of Acetylation Mimics in the 27 Lysine Residues of Human Tau Protein

BACKGROUND

Various neurodegenerative diseases, including Alzheimer's disease (AD), are related to abnormal hyperphosphorylated microtubule-associated protein tau accumulation in brain lesions. Recent studies have focused on toxicity caused by another post-translational modification (PTM), acetylation of the lysine (K) residues of tau protein. Because there are numerous acetylation sites, several studies have introduced mimics of tau acetylation using amino acid substitutions from lysine to glutamine (Q). However, human tau protein contains over 20 acetylation sites; thus, investigation of the effects of an acetylated tau is difficult.

OBJECTIVE

Here, the authors in silico evaluated acetylation effects using SIFT, PolyPhen-2 and PROVEAN which can estimate the effects of amino acid substitutions based on the sequence homology or protein structure in tau isoforms. In addition, they also investigated 27 acetylation effects on the amyloid formation of tau proteins using Waltz.

RESULTS

15 acetylation mimics were estimated to be the most detrimental, which indicates that there may be novel pathogenic acetylation sites in the human tau protein. Interestingly, the deleterious effect of acetylation mimics was different according to the type of isoforms. Furthermore, all acetylation mimics were predicted to be a region of amyloid formation at the codons 274-279 of human tau protein. Notably, acetylation mimic of codon 311 (K311Q) induced the formation of an additional amyloid region located on codons 306-311 of the human tau protein.

CONCLUSION

To the best of our knowledge, this is the first simultaneous in-silico evaluation of the acetylation state of 27 human tau protein residues.

Prediction and evaluation of deleterious and disease causing non-synonymous SNPs (nsSNPs) in human NF2 gene responsible for neurofibromatosis type 2 (NF2)

The majority of genetic variations in the human genome that lead to variety of different diseases are caused by non-synonymous single nucleotide polymorphisms (nsSNPs). Neurofibromatosis type 2 (NF2) is a deadly disease caused by nsSNPs in the NF2 gene that encodes for a protein called merlin. This study used various in silico methods, SIFT, Polyphen-2, PhD-SNP and MutPred, to investigate the pathogenic effect of 14 nsSNPs in the merlin FERM domain. The G197C and L234R mutations were found to be two deleterious and disease mutations associated with the mild and severe forms of NF2, respectively. Molecular dynamics (MD) simulations were conducted to understand the stability, structure and dynamics of these mutations. Both mutant structures experienced larger flexibility compared to the wildtype. The L234R mutant suffered from more prominent structural instability, which may help to explain why it is associated with the more severe form of NF2. The intramolecular hydrogen bonding in L234R mutation decreased from the wildtype, while intermolecular hydrogen bonding of L234R mutation with solvent greatly increased. The native contacts were also found to be important. Protein-protein docking revealed that L234R mutation decreased the binding complementarity and binding affinity of LATS2 to merlin, which may have an impact on merlin's ability to regulate the Hippo signaling pathway. The calculated binding affinity of the LATS2 to L234R mutant and wildtype merlin protein is found to be 21.73 and -11 kcal/mol, respectively. The binding affinity of the wildtype merlin agreed very well with the experimental value, -8 kcal/mol.Communicated by Ramaswamy H. Sarma.