Prediction of the Impact of Deleterious Nonsynonymous Single Nucleotide Polymorphisms on the Human RRM2B Gene: A Molecular Modeling Study

Investigation of pathogenic germline variants in gastric cancer and development of "GasCanBase" database

BACKGROUND

Gastric cancer, which is also known as stomach cancer, can be influenced by both germline and somatic mutations. Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in germline have long been reported to play a pivotal role in cancer progression.

AIM

The aim of this study is to examine the nsSNP in GC-associated genes. The study also aims to develop a database with extensive information regarding the nsSNPs in the GC-associated genes and their impacts.

METHODS AND RESULTS

A total of 34,588 nsSNPs from 1,493,460 SNPs of the 40 genes were extracted from the available SNP database. Drug binding and energy minimization were examined by molecular docking and YASARA. To validate the existence of the germline CDH1 gene mutation (rs34466743) in the isolated blood DNA of gastric cancer (GC) patients, polymerase chain reaction (PCR) and DNA sequencing were performed. According to the results of the gene network analysis, 17 genes may interact with other types of cancer. A total of 11,363 nsSNPs were detected within the 40 GC genes. Among these, 474 nsSNPs were predicted to be damaging and 40 to be the most damaging. The SNPs in domain regions were thought to be strong candidates that alter protein functions. Our findings proposed that most of the selected nsSNPs were within the domains or motif regions. Free Energy Deviation calculation of protein structure pointed toward noteworthy changes in the structure of each protein that can demolish its natural function. Subsequently, drug binding confirmed the structural variation and the ineffectiveness of the drug against the mutant model in individuals with these germline variants. Furthermore, in vitro analysis of the rs34466743 germline variant from the CDH1 gene confirmed the strength and robustness of the pipeline that could expand the somatic alteration for causing cancer. In addition, a comprehensive gastric cancer polymorphism database named "GasCanBase" was developed to make data available to researchers.

CONCLUSION

The findings of this study and the "GasCanBase" database may greatly contribute to our understanding of molecular epidemiology and the development of precise therapeutics for gastric cancer. GasCanBase is available at: https://www.gascanbase.com/.

In silico exploration and molecular dynamics of deleterious SNPs on the human TERF1 protein triggering male infertility

By limiting chromosome erosion and end-to-end fusions, telomere integrity is critical for chromosome stability and cell survival. During mitotic cycles or due to environmental stresses, telomeres become progressively shorter and dysfunctional, thus triggering cellular senescence, genomic instability and cell death. To avoid such consequences, the telomerase action, as well as the Shelterin and CST complexes, assure the telomere's protection. Telomeric repeat binding factor 1 (TERF1), which is one of the primary components of the Shelterin complex, binds directly to the telomere and controls its length and function by regulating the telomerase activity. Several reports about TERF1 gene variations have been associated with different diseases, and some of them have linked these variations to male infertility. Hence, this paper can be advantageous to investigate the association between the missense variants of the TERF1 gene and the susceptibility to male infertility. The stepwise prediction of SNPs pathogenicity followed in this study was based on stability and conservation analysis, post-translational modification, secondary structure, functional interaction prediction, binding energy evaluation and finally molecular dynamic simulation. Prediction matching among the tools revealed that out of 18 SNPs, only four (rs1486407144, rs1259659354, rs1257022048 and rs1320180267) were predicted as the most damaging and highly deleterious SNPs affecting the TERF1 protein and its molecular dynamics when interacting with the TERB1 protein by influencing the function, structural stability, flexibility and compaction of the overall complex. Interestingly, these polymorphisms should be considered during genetic screening so they can be used effectively as genetic biomarkers for male infertility diagnosis.Communicated by Ramaswamy H. Sarma.

Mitochondrial DNA Depletion Syndrome and Its Associated Cardiac Disease

Mitochondria is a ubiquitous, energy-supplying (ATP-based) organelle found in nearly all eukaryotes. It acts as a “power plant” by producing ATP through oxidative phosphorylation, providing energy for the cell. The bioenergetic functions of mitochondria are regulated by nuclear genes (nDNA). Mitochondrial DNA (mtDNA) and respiratory enzymes lose normal structure and function when nuclear genes encoding the related mitochondrial factors are impaired, resulting in deficiency in energy production. Massive generation of reactive oxygen species and calcium overload are common causes of mitochondrial diseases. The mitochondrial depletion syndrome (MDS) is associated with the mutations of mitochondrial genes in the nucleus. It is a heterogeneous group of progressive disorders characterized by the low mtDNA copy number. TK2, FBXL4, TYPM, and AGK are genes known to be related to MDS. More recent studies identified new mutation loci associated with this disease. Herein, we first summarize the structure and function of mitochondria, and then discuss the characteristics of various types of MDS and its association with cardiac diseases.

Genome-wide discovery of DNA polymorphisms via resequencing of chickpea cultivars with contrasting response to drought stress

Drought stress limits plant growth, resulting in a significant yield loss in chickpea. The diversification in genome sequence and selective sweep of allele(s) in different genotypes of a crop plant may play an important role in the determination of agronomic traits, including drought stress response. We investigated, via whole genome resequencing, the DNA polymorphisms between two sets of chickpea genotypes with contrasting drought stress responses (3 drought-sensitive vs. 6 drought-tolerant). In total, 36,406 single nucleotide polymorphisms (SNPs) and 3407 insertions or deletions (InDels) differentiating drought-sensitive and drought-tolerant chickpea genotypes were identified. Interestingly, most (91%) of these DNA polymorphisms were located in chromosomes 1 and 4. The genes harboring DNA polymorphisms in their promoter and/or coding regions and exhibiting differential expression under control and/or drought stress conditions between/within the drought-sensitive and tolerant genotypes were found implicated in the stress response. Furthermore, we identified DNA polymorphisms within the cis-regulatory motifs in the promoter region of abiotic stress-related and QTL-associated genes, which might contribute to the differential expression of the candidate drought-responsive genes. In addition, we revealed the effect of nonsynonymous SNPs on mutational sensitivity and stability of the encoded proteins. Taken together, we identified DNA polymorphisms having relevance in drought stress response and revealed candidate genes to engineer drought tolerance in chickpea.

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.