Non-synonymous Single Nucleotide Polymorphisms Research Articles

Evaluating the Genetic Landscape of (Plasmodium falciparum) PfKelch13 Gene Polymorphisms in Côte d'Ivoire Following a Decade of Artemisinin-based Combination Therapy

Introduction: Artemisinin-based combination therapies (ACTs) are the mainstay of malaria treatment globally. However, their effectiveness is threatened by the emergence of resistance in Plasmodium falciparum (P.f), particularly in Southeast Asia (SEA). Specific mutations within the pfKelch13 gene, such as Cys-580-Tyr, Arg-539-Thr, Tyr-493-His, and Ile-543-Thr, have been strongly linked to delayed parasite clearance following ACT treatment. This study aimed to investigate polymorphisms within the pfKelch13 gene (also known as the K13-propeller or K13 gene) in four regions of Côte d'Ivoire. Côte d'Ivoire experiences year-round malaria transmission and has utilized ACTs for over a decade. Methods: From September 2013 to July 2014, samples were collected from patients residing in Abidjan (south), Ayamé (southeast), Man (west), and Korhogo (north) who presented with microscopically confirmed uncomplicated malaria. Following Plasmodium falciparum DNA extraction, nested PCR was employed to amplify an 849 bp fragment of the pfKelch13 gene. Amplicons were subsequently sequenced and analyzed using BioEdit software. Results: 531 DNA sequences were analyzed including 301 (58.7%) from Abidjan, 61(11.5%) from Ayamé, 93 (17.5%) from Man and 76 (13.4%) from Korhogo. Only 20 isolates carrying 22 mutations were observed including 6 non-synonymous single-nucleotide polymorphisms (nsSNP): 4 in Abidjan (Asp-535-Met; Ala-578-Ser; Phe-583-Ser and Ile-601-Thr) and 2 in Korhogo (Asp-559-Asn and Val-510-Met). Only synonymous SNP (sSNP) were identified in the two other Towns. The proportion of mutated pfK13 sequences is 3.8% (20/531). Conclusion: The identification of non-synonymous mutations in this study underscores the importance of heightened surveillance for potential ACT resistance in Plasmodium falciparum within Côte d'Ivoire. Combining in vitro assays, such as the Ring-stage Survival Assay, with molecular testing will be crucial for definitively determining the phenotypic impact of these mutations on parasite susceptibility to ACTs.

Computational Analysis of S1PR1 SNPs Reveals Drug Binding Modes Relevant to Multiple Sclerosis Treatment

Background/Objectives: Multiple sclerosis (MS) is an autoimmune disorder of the central nervous system (CNS) characterized by myelin and axonal damage with a globally rising incidence. While there is no known cure for MS, various disease-modifying treatments (DMTs) exist, including those targeting Sphingosine-1-Phosphate Receptors (S1PRs), which play important roles in immune response, CNS function, and cardiovascular regulation. This study focuses on understanding how nonsynonymous single nucleotide polymorphisms (rs1299231517, rs1323297044, rs1223284736, rs1202284551, rs1209378712, rs201200746, and rs1461490142) in the S1PR1’s active site affect the binding of endogenous ligands, as well as different drugs used in MS management. Methods: Extensive molecular dynamics simulations and linear interaction energy (LIE) calculations were employed to predict binding affinities and potentially guide future personalized medicinal therapies. The empirical parameters of the LIE method were optimized using the binding free energies calculated from experimentally determined IC50 values. These optimized parameters were then applied to calculate the binding free energies of S1P to mutated S1PR1, which correlated well with experimental values, confirming their validity for assessing the impact of SNPs on S1PR1 binding affinities. Results: The binding free energies varied from the least favorable −8.2 kcal/mol for the wild type with ozanimod to the most favorable −16.7 kcal/mol for the combination of siponimod with the receptor carrying the F2055.42L mutation. Conclusions: We successfully demonstrated the differences in the binding modes, interactions, and affinities of investigated MS drugs in connection with SNPs in the S1PR1 binding site, resulting in several viable options for personalized therapies depending on the present mutations.

The landscape of sequence variations between resistant and susceptible hot peppers to predict functional candidate genes against bacterial wilt disease

BackgroundBacterial wilt (BW), caused by Ralstonia solanacearum (Ral), results in substantial yield losses in pepper crops. Developing resistant pepper varieties through breeding is the most effective strategy for managing BW. To achieve this, a thorough understanding of the genetic information connected with resistance traits is essential. Despite identifying three major QTLs for bacterial wilt resistance in pepper, Bw1 on chromosome 8, qRRs-10.1 on chromosome 10, and pBWR-1 on chromosome 1, the genetic information of related BW pepper varieties has not been sufficiently studied. Here, we resequenced two pepper inbred lines, C. annuum ‘MC4’ (resistant) and C. annuum ‘Subicho’ (susceptible), and analyzed genomic variations through SNPs and Indels to identify candidate genes for BW resistance.ResultsAn average of 139.5 Gb was generated among the two cultivars, with coverage ranging from 44.94X to 46.13X. A total of 8,815,889 SNPs was obtained between ‘MC4’ and ‘Subicho’. Among them, 31,190 (0.35%) were non-synonymous SNPs (nsSNPs) corresponding to 10,926 genes, and these genes were assigned to 142 Gene Ontology (GO) terms across the two cultivars. We focused on three known BW QTL regions by identifying genes with sequence variants through gene set enrichment analysis and securing those belonging to high significant GO terms. Additionally, we found 310 NLR genes with nsSNP variants between ‘MC4’ (R) and ‘Subicho’ (S) within these regions. Also, we performed an Indel analysis on these genes. By integrating all this data, we identified eight candidate BW resistance genes, including two NLR genes with nsSNPs variations in qRRs-10.1 on chromosome 10.ConclusionWe identified genomic variations in the form of SNPs and Indels by re-sequencing two pepper cultivars with contrasting traits for bacterial wilt. Specifically, the four genes associated with pBWR-1 and Bw1 that exhibit both nsSNP and Indel variations simultaneously in ‘Subicho’, along with the two NLR genes linked to qRRs-10.1, which are known for their direct involvement in immune responses, are identified as most likely BW resistance genes. These variants in leading candidate genes associated with BW resistance can be used as important markers for breeding pepper varieties.

Impact of missense mutations on the structure–function relationship of human succinyl-CoA synthetase using in silico analysis

The encephalomyopathic mtDNA depletion syndrome with methylmalonic aciduria is associated with succinyl-CoA synthetase (SCS) deficiency caused by pathogenic variants in genes encoding its two subunits. SCS is a mitochondrial enzyme involved in several metabolic pathways and acts as a heterodimer composed of α and β subunits encoded by SUCLG1 and SUCLA2 genes, respectively. The purpose of this study was to analyze the effects of the most pathogenic non-synonymous single nucleotide polymorphisms (nsSNPs) by applying, using different prediction tools, a filtering strategy, on the 343 and 365 nsSNPs found in SUCLG1 and SUCLA2 genes, respectively, retrieved from the databases, then to evaluate their structural and functional effects using homology modeling and molecular docking. Results showed that most deleterious mutations selected for structural analysis were located in loop regions critical for protein stability and function, especially, variants altering glycine and proline residues in these regions supporting their importance. We also showed that variants leading to hydrophobic and hydrophilic residues can destabilize the folding and binding of the protein. Molecular docking has also been used to identify the most important regions of ligand binding site (CoA binding site, ADP-Mg2+ binding site and phosphate ion binding site) and between the two subunits themselves, which mainly involving the ligase CoA domain. Our structural analysis, performed on selected nsSNP, are in accordance with experimental studies reported in the literature and predicted that they would responsible to either nonfunctional protein, subunit instability resulting in reduced amounts of misassembled protein, or in a protein unable to phosphorylate ADP.

Transcriptional variation and RNA polymorphism among different Lentinula edodes (Berk.) Pegler strains.

Lentinula edodes is a commercially important mushroom known for its nutritional and therapeutic values. However, the molecular mechanisms underlying the distinct nutritional and physiological attributes of various L. edodes strains are not well understood. This study focused on three Lentinula strains (DMRO-356, DMRO-623, and DMRO-388s) with different nutritional and productivity profiles. Illumina sequencing was used to perform a whole-transcriptome analysis, conducting 100-base pair paired-end sequencing of total messenger RNA (mRNA) in duplicate, resulting in 28-48 million sequencing reads per strain. After rigorous data filtering, over 99% of high-quality reads were retained, and more than 95% were aligned to the Lentinula genome. Differential gene expression analyses identified 2210 differentially expressed genes between DMRO-356 and DMRO-623, 862 between DMRO-356 and DMRO-388s, and 2212 between DMRO-623 and DMRO-388s. Significant genetic variations were found among the strains, including 7753 single nucleotide polymorphisms (SNPs) in DMRO-356 versus DMRO-623 and 4080 SNPs in DMRO-356 versus DMRO-388s. Additionally, 349 insertions/deletions (InDels) were found in DMRO-356/DMRO-623 and 218 in DMRO-356/DMRO-388 s. Non-synonymous SNPs, which alter amino acid compositions, were analyzed, showing a preference for polar over charged amino acids. These differentially expressed genes were associated with various nutritional and developmental processes, highlighting the importance of genetic variations in shaping amino acid composition and potentially affecting protein function. This study is the first comprehensive exploration of transcriptional differences among Lentinula strains available for its cultivation, providing valuable insights to enhance mushroom quality and productivity. © 2024 Society of Chemical Industry.

Elucidating the role of MLL1 nsSNPs: Structural and functional alterations and their contribution to leukemia development.

The Mixed lineage leukemia 1 (MLL1) gene, located on chromosome 11q23, plays a pivotal role in histone lysine-specific methylation and is consistently associated with various types of leukemia. Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) have been tied to numerous diseases, including cancers, and have become valuable cancer biomarkers. There's a notable gap in studies probing the influence of SNPs on MLL1 protein structure, function, and subsequent modifications. We utilized an array of bioinformatics tools, including PredictSNP, InterPro, ConSurf, I-Mutant2.0, MUpro, Musitedeep, Project HOPE, RegulomeDB, Mutpred2, and both CScape and CScape Somatic, to meticulously analyze the consequences of nsSNPs in the MLL1 gene. Out of 2,097 nsSNPs analyzed, 62 were determined to be significantly pathogenic by the PredictSNP tool, with ten crucial MLL1 functional domains identified using InterPro. Additionally, 50 of these nsSNPs had high conservation scores, hinting at potential effects on protein structure and function, while 32 were found to undermine MLL1 protein stability. Notably, four nsSNPs were deemed oncogenic, with two identified as cancer drivers. The nsSNP, D2724G, between the MLL1 protein's FY-rich domains, could disrupt proteolytic cleavage, altering gene expression patterns and potentially promoting cancer. Our research provides a comprehensive assessment of nsSNPs' impact in the MLL1 protein structure and function and consequently on leukemia development, suggesting potential avenues for personalized treatment, early detection, improved prognosis, and a deeper understanding of hematological malignancy genesis.

Bioinformatic Characterization of the Functional and Structural Effect of Single Nucleotide Mutations in Patients with High-Grade Glioma.

Background: Gliomas are neoplasms of the central nervous system that originate in glial cells. The genetic characteristics of this type of neoplasm are the loss of function of tumor suppressor genes such as TP53 and somatic mutations in genes such as IDH1/2. Additionally, in clinical cases, de novo single nucleotide polymorphisms (SNP) are reported, of which their pathogenicity and their effects on the function and stability of the protein are known. Methodology: Non-synonymous SNPs were analyzed for their structural and functional effect on proteins using a set of bioinformatics tools such as SIFT, PolyPhen-2, PhD-SNP, I-Mutant 3.0, MUpro, and mutation3D. A structural comparison between normal and mutated residues for disease-associated coding SNPs was performed using TM-aling and the SWISS MODEL. Results: A total of 13 SNPs were obtained for the TP53 gene, 1 SNP for IDH1, and 1 for IDH2, which would be functionally detrimental and associated with disease. Additionally, these changes compromise the structure and function of the protein; the A161S SNP for TP53 that has not been reported in any databases was classified as detrimental. Conclusions: All non-synonymous SNPs reported for TP53 were in the region of the deoxyribonucleic acid (DNA) binding domain and had a great impact on the function and stability of the protein. In addition, the two polymorphisms detected in IDH1 and IDH2 genes compromise the structure and activity of the protein. Both genes are related to the development of high-grade gliomas. All the data obtained in this study must be validated through experimental approaches.

Whole-Genome Resequencing Reveals Significant Genetic Differentiation Between Exserohilum turcicum Populations from Maize and Sorghum and Candidate Effector Genes Related to Host Specificity.

Exserohilum turcicum is a devastating fungal pathogen that infects both maize and sorghum, leading to severe leaf diseases of the two crops. According to host specificity, pathogenic isolates of E. turcicum are divided into two formae speciales, namely E. turcicum f. sp. zeae and E. turcicum f. sp. sorghi. To date, the molecular mechanism underlying the host specificity of E. turcicum is marginally known. In this study, the whole genomes of 60 E. turcicum isolates collected from both maize and sorghum were resequenced, which enabled identification of 233,022 single-nucleotide polymorphisms (SNPs) in total. Phylogenetic analysis indicated that all isolates are clustered into four genetic groups that have a close relationship with host source. This observation is validated by the result of principal component analysis. Analysis of population structure revealed that there is obvious genetic differentiation between two populations from maize and sorghum. Further analysis showed that 5,431 SNPs, including 612 nonsynonymous SNPs, are completely co-segregated with the host source. These nonsynonymous SNPs are located in 539 genes, among which 18 genes are predicted to encode secretory proteins, including six putative effector genes named SIX13-like, Ecp6, GH12, GH28-1, GH28-2, and CHP1. Sequence polymorphism analysis revealed various numbers of SNPs in the coding regions of these genes. These findings provide new insights into the molecular basis of host specificity in E. turcicum.

Genetic polymorphism of Duffy binding protein in Pakistan Plasmodium vivax isolates

Plasmodium vivax Duffy binding protein (PvDBP) is crucial for erythrocyte invasion, interacting with the Duffy Antigen Receptor for Chemokines (DARC) on the erythrocyte surface. The amino-terminal cysteine-rich region II of PvDBP (PvDBPII) is a promising blood stage vaccine candidate, yet the genetic polymorphisms of this protein in global P. vivax isolates complicate the design of effective vaccines against vivax malaria. This study analyzed the genetic polymorphism of PvDBPII in Pakistan P. vivax isolates. A total of 29 single nucleotide polymorphisms (SNPs), including 22 nonsynonymous SNPs, were identified in 118 Pakistan PvDBPII. Most amino acid substitutions occurred in subdomains II and III, with six commonly observed in the global PvDBPII population. The amino acid change patterns in Pakistan PvDBPII generally mirrored those in global PvDBPII, although the frequencies of amino acid changes varied by country. Nucleotide diversity in Pakistan PvDBPII was comparable to that found in global PvDBPII. Evidence of natural selection and recombination in Pakistan PvDBPII aligned with observations in global PvDBPII. Analysis of the haplotype network of global PvDBPII revealed a complexed network of 167 haplotypes, but no geographical clustering was observed. The findings are crucial for understanding the genetic characteristics of Pakistan PvDBPII. A comprehensive analysis of nucleotide diversity and evolutionary trends in the global PvDBPII population offers valuable insights for the development of vivax malaria vaccines based on this antigen.

Molecular characteristics of drug-susceptible Mycobacterium tuberculosis clinical isolates based on treatment duration

Objectives: In this study, we performed comparative genomic and transcriptomic analysis of clinical isolates of Mycobacterium tuberculosis collected from patients with drug-susceptible tuberculosis (DS-TB). The clinical isolates were categorized based on treatment duration: standard 6 months or >6 months.Methods: Study participants were recruited from a 2016 to 2018 tuberculosis cohort, and clinical M. tuberculosis isolates were collected from the sputum of patients with tuberculosis. We analyzed the genome and transcriptome of the isolated M. tuberculosis.Results: Genomic analysis revealed a specific non-synonymous single-nucleotide polymorphism in pe_pgrs9 and ppe34, exclusive to the group treated for >6 months. Transcriptomic analysis revealed increased expression of various virulence-associated protein family genes and decreased expression of ribosomal protein genes and ppe38 genes in the group treated for >6 months.Conclusion: The identified genetic variation and gene expression patterns may influence treatment outcomes by modulating host immune responses, increasing virulence, and potentially contributing to persister cell formation in M. tuberculosis. This study provides insights into the genetic and transcriptomic factors associated with prolonged DS-TB treatment. However, our study identified molecular characteristics using a small sample size, and further detailed studies are warranted.

Identification and Differentiation of the Fusarium graminearum NX-2 Chemotype Using High-Resolution Melting (HRM).

Fusarium head blight causes significant yield losses in wheat and other cereals and contaminates grain products with trichothecene mycotoxins. F. graminearum isolates are classified into different chemotypes depending on the type of mycotoxin produced, including the type B trichothecenes 3-acetyl deoxynivalenol (3-ADON), 15-acetyl deoxynivalenol (15-ADON), nivalenol (NIV), and the recently identified type A trichothecene NX-2. Molecular tools to differentiate NX-2 producers from other chemotypes have remained relatively laborious and time consuming. In this study, we developed and validated a high-resolution melting (HRM) assay that can identify NX-2 producers quickly and cost-effectively. By analyzing TRI1 coding sequences from 183 geographically diverse isolates representing all four F. graminearum chemotypes, we selected a 75-base pair region containing four non-synonymous single nucleotide polymorphisms (SNPs) that are specific to the NX-2 genotypes. The amplicon generated two HRM profiles, one of which was specific for only NX-2. We confirmed that the assay is robust across qPCR platforms and unambiguously differentiates NX-2 from other chemotypes using a panel of 72 diverse isolates previously collected from North America. The HRM assay was also successful in identifying NX-2 producers directly from DNA extracted from infected wheat spikes with varying levels of disease severity and fungal DNA. The assay can detect as little as 0.01 ng of fungal DNA in a background of 50 ng of plant DNA. This new diagnostic assay can be used for high-throughput molecular detection of the NX-2 chemotype of F. graminearum from infected plant samples and culture collections, thus making it a valuable tool for surveys of contemporary and historical FHB pathogen populations.

Missense variant rs75603675 within TMPRSS2 gene is associated with the increased risk of severe form of COVID-19

Host genetics among other factors play an important role in COVID-19 disease severity. TMPRSS2, a serine protease, facilitates the priming of the SARS-CoV-2 spike protein, which is essential for the virus to enter the host cell. Several studies had targeted the TMPRSS2 polymorphisms with respect to SARS-CoV-2 infection and COVID-19 disease severity. Initially, the Whole Exome Sequencing (WES) of 7 healthy individuals and 22 COVID-19 patients with different degrees of disease severity was conducted to find the mutational landscape in different genes. A total of 26 single nucleotide polymorphisms (SNPs) were identified of which six were within the exons of TMPRSS2 (four synonymous and two nonsynonymous variants) while the rest of the 20 SNPs were recognized within the flanking intronic regions of TMPRSS2 gene. The nonsynonymous SNPs, rs75603675 and rs12329760, were further evaluated for association with disease severity in a larger sample size of 120 individuals by PCR followed by Sanger sequencing. Neither allelic nor genotypic distributions of rs12329760 were significantly associated with COVID-19 disease severity. However, individuals harboring the A allele of rs75603675 was found to have higher risk of severe COVID-19 compared to the C allele [OR (95%CI): 1.95 (1.11–3.39), p = 0.019]. Also, the genotype A/A [OR (95%CI): 5.13 (1.00–26.38), p = 0.033] of rs75603675 was associated with increased risk of severe COVID-19 under the recessive genetic inheritance model. Although the impact of COVID-19 pandemic has waned due to vaccination and public health measures, continued research on the association of COVID-19 disease severity and infection susceptibility with host genetics is required to shed valuable insights on the future long-term health effects of COVID-19 and impact of new variants on different populations and enable implication of proper informed healthcare strategies during future public health crises.

Unlocking the Door for Precision Medicine in Rare Conditions: Structural and Functional Consequences of Missense ACVR1 Variants.

Rare diseases and conditions have thus far received relatively less attention in the field of precision/personalized medicine than common chronic diseases. There is a dire need for orphan drug discovery and therapeutics in ways that are informed by the precision/personalized medicine scholarship. Moreover, people with rare conditions, when considered collectively across diseases worldwide, impact many communities. In this overarching context, Activin A Receptor Type 1 (ACVR1) is a transmembrane kinase from the transforming growth factor-β superfamily and plays a critical role in modulating the bone morphogenetic protein signaling. Missense variants of the ACVR1 gene result in modifications in structure and function and, by extension, abnormalities and have been predominantly linked with two rare conditions: fibrodysplasia ossificans progressiva and diffuse intrinsic pontine glioma. We report here an extensive bioinformatic analyses assessing the pool of 50,951 variants and forecast seven highly destabilizing mutations (R206H, G356D, R258S, G328W, G328E, R375P, and R202I) that can significantly alter the structure and function of the native protein. Protein-protein interaction and ConSurf analyses revealed the crucial interactions and localization of highly deleterious mutations in highly conserved domains that may impact the binding and functioning of the protein. cBioPortal, CanSAR Black, and existing literature affirmed the association of these destabilizing mutations with posterior fossa ependymoma, uterine corpus carcinoma, and pediatric brain cancer. The current findings suggest these deleterious nonsynonymous single nucleotide polymorphisms as potential candidates for future functional annotations and validations associated with rare conditions, further aiding the development of precision medicine in rare diseases.

In Silico Screening, Molecular Dynamics Simulation and Binding Free Energy Identify Single-Point Mutations That Destabilize p53 and Reduce Binding to DNA.

Considering p53's pivotal role as a tumor suppressor protein, proactive identification and characterization of potentially harmful p53 mutations are crucial before they appear in the population. To address this, four computational prediction tools-SIFT, Polyphen-2, PhD-SNP, and MutPred2-utilizing sequence-based and machine-learning algorithms, were employed to identify potentially deleterious p53 nsSNPs (nonsynonymous single nucleotide polymorphisms) that may impact p53 structure, dynamics, and binding with DNA. These computational methods identified three variants, namely, C141Y, C238S, and L265P, as detrimental to p53 stability. Furthermore, molecular dynamics (MD) simulations revealed that all three variants exhibited heightened structural flexibility compared to the native protein, especially the C141Y and L265P mutations. Consequently, due to the altered structure of mutant p53, the DNA-binding affinity of all three variants decreased by approximately 1.8 to 9.7 times compared to wild-type p53 binding with DNA (14 μM). Notably, the L265P mutation exhibited an approximately ten-fold greater reduction in binding affinity. Consequently, the presence of the L265P mutation in p53 could pose a substantial risk to humans. Given that p53 regulates abnormal tumor growth, this research carries significant implications for surveillance efforts and the development of anticancer therapies.

Fine genetic mapping and transcriptomic analysis revealed major gene modulating the clear stripe margin pattern of watermelon peel.

The peel stripe margin pattern is one of the most important quality traits of watermelon. In this study, two contrasted watermelon lines [slb line (P1) with a clear peel stripe margin pattern and GWAS-38 line (P2) with a blurred peel stripe margin pattern] were crossed, and biparental F2 mapping populations were developed. Genetic segregation analysis revealed that a single recessive gene is modulating the main-effect genetic locus (Clcsm) of the clear stripe margin pattern of peel. Bulked segregant analysis-based sequencing (BSA-Seq) and fine genetic mapping exposed the delimited Clcsm locus to a 19.686-kb interval on chromosome 6, and the Cla97C06G126680 gene encoding the MYB transcription factor family was identified. The gene mutation analysis showed that two non-synonymous single-nucleotide polymorphism (nsSNP) sites [Chr6:28438793 (A-T) and Chr6:28438845 (A-C)] contribute to the clear peel stripe margin pattern, and quantitative real-time polymerase chain reaction (qRT-PCR) also showed a higher expression trend in the slb line than in the GWAS-38 line. Further, comparative transcriptomic analysis identified major differentially expressed genes (DEGs) in three developmental periods [4, 12, and 20 days after pollination (DAP)] of both parental lines. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses indicated highly enriched DEGs involved in metabolic processes and catalytic activity. A total of 44 transcription factor families and candidate genes belonging to the ARR-B transcription factor family are believed to regulate the clear stripe margin trait of watermelon peel. The gene structure, sequence polymorphism, and expression trends depicted significant differences in the peel stripe margin pattern of both parental lines. The ClMYB36 gene showed a higher expression trend for regulating the clear peel stripe margin of the slb line, and the ClAPRR5 gene depicted a higher expression for modulating the blurred peel stripe margin in the GWAS-38 line. Overall, our fine genetic mapping and transcriptomic analysis revealed candidate genes differentiating the clear and blurred peel stripe patterns of watermelon fruit.

Surviving high temperatures: The crucial role of vesicular inhibitory amino acid transporter in Asian honeybee, Apis cerana

Asian honeybees (Apis cerana) play a crucial role as pollinators to service for the ecological stability. However, their proliferation and growth are significantly impacted by environmental temperature stress. This study delves into the function of the Apis cerana vesicular inhibitory amino acid transporter gene (AcVIAAT) in safeguarding Asian honeybees against high-temperature stress. The AcVIAAT gene exhibits positive responsiveness in honeybees subjected to varying thermal conditions by triggering the genes associated with oxidative stress. Molecular docking, co-immunoprecipitation, and ELISA verify the capacity of the AcVIAAT protein to interact with γ-aminobutyric acid (GABA), a key inhibitory neurotransmitter. Administering GABA to honeybees significantly improves their survival rate under high-temperature stress and also simultaneously upregulating oxidative stress-related genes. Therefore, these findings reveal that the AcVIAAT gene enhances the thermoregulatory capacity of honeybees by modulating oxidative stress-related genes through facilitating GABA transport. The characterization of six non-synonymous SNPs in the AcVIAAT gene among A.cerana populations distributed across both the northern and southern regions indicates a potential association between gene variation and environmental adaptation. Our results contribute to elucidating the molecular mechanisms underlying high-temperature tolerance in Asian honeybees and provide a promising genetic marker for enhancing heat tolerance through genetic improvement.

Unraveling the potential effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on the Protein structure and function of the human SLC30A8 gene on type 2 diabetes and colorectal cancer: An In silico approach

Background and aimsThe single nucleotide polymorphisms (SNPs) in SLC30A8 gene have been recognized as contributing to type 2 diabetes (T2D) susceptibility and colorectal cancer. This study aims to predict the structural stability, and functional impacts on variations in non-synonymous SNPs (nsSNPs) in the human SLC30A8 gene using various computational techniques. Materials and methodsSeveral in silico tools, including SIFT, Predict-SNP, SNPs&GO, MAPP, SNAP2, PhD-SNP, PANTHER, PolyPhen-1,PolyPhen-2, I-Mutant 2.0, and MUpro, have been used in our study. ResultsAfter data analysis, out of 336 missenses, the eight nsSNPs, namely R138Q, I141N, W136G, I349N, L303R, E140A, W306C, and L308Q, were discovered by ConSurf to be in highly conserved regions, which could affect the stability of their proteins. Project HOPE determines any significant molecular effects on the structure and function of eight mutated proteins and the three-dimensional (3D) structures of these proteins. The two pharmacologically significant compounds, Luzonoid B and Roseoside demonstrate strong binding affinity to the mutant proteins, and they are more efficient in inhibiting them than the typical SLC30A8 protein using Autodock Vina and Chimera. Increased binding affinity to mutant SLC30A8 proteins has been determined not to influence drug resistance. Ultimately, the Kaplan-Meier plotter study revealed that alterations in SLC30A8 gene expression notably affect the survival rates of patients with various cancer types. ConclusionFinally, the study found eight highly deleterious missense nsSNPs in the SLC30A8 gene that can be helpful for further proteomic and genomic studies for T2D and colorectal cancer diagnosis. These findings also pave the way for personalized treatments using biomarkers and more effective healthcare strategies.

Altered O-glycosylation of β1-adrenergic receptor N-terminal single-nucleotide variants modulates receptor processing and functional activity.

N-terminal nonsynonymous single-nucleotide polymorphisms (SNPs) of G protein-coupled receptors (GPCRs) are common and often affect receptor post-translational modifications. Their functional implications are, however, largely unknown. We have previously shown that the human β1-adrenergic receptor (β1AR) is O-glycosylated in the N-terminal extracellular domain by polypeptide GalNAc transferase-2 that co-regulates receptor proteolytic cleavage. Here, we demonstrate that the common S49G and the rare A29T and R31Q SNPs alter these modifications, leading to distinct effects on receptor processing. This was achieved by in vitro O-glycosylation assays, analysis of native receptor N-terminal O-glycopeptides, and expression of receptor variants in cell lines and neonatal rat ventricular cardiomyocytes deficient in O-glycosylation. The SNPs eliminated (S49G) or introduced (A29T) regulatory O-glycosites that enhanced or inhibited cleavage at the adjacent sites (P52↓L53 and R31↓L32), respectively, or abolished the major site at R31↓L32 (R31Q). The inhibition of proteolysis of the T29 and Q31 variants correlated with increased full-length receptor levels at the cell surface. Furthermore, the S49 variant showed increased isoproterenol-mediated signaling in an enhanced bystander bioluminescence energy transfer β-arrestin2 recruitment assay in a coordinated manner with the common C-terminal R389G polymorphism. As Gly at position 49 is ancestral in placental mammals, the results suggest that its exchange to Ser has created a β1AR gain-of-function phenotype in humans. This study provides evidence for regulatory mechanisms by which GPCR SNPs outside canonical domains that govern ligand binding and activation can alter receptor processing and function. Further studies on other GPCR SNPs with clinical importance as drug targets are thus warranted.

Comparative Analysis of HPV16 Variants in the Untranslated Regulatory Region, L1, and E6 Genes among Vaccinated and Unvaccinated Young Women: Assessing Vaccine Efficacy and Viral Diversity.

HPV16 is occasionally detected in vaccinated women who received the bivalent HPV16/18 vaccine, usually at low viral loads. This study explored potential differences in HPV16 variants between vaccinated and unvaccinated women. HPV16-postive viral loads were detected in 1.9% (17/875) and 13% (162/760) of vaccinated and unvaccinated women, respectively, showcasing the vaccine's high efficacy. The L1, E6, and URR regions of HPV16 were sequenced from genital swabs from 16 vaccinated and 25 unvaccinated women in the HAVANA (HPV Among Vaccinated And Non-vaccinated Adolescents) study. The majority of HPV16 variants from vaccinated and unvaccinated women clustered similarly with sub-lineages A1 and A2. Additionally, a separate cluster within lineage A was found, with the variants sharing the L1-located SNP A753G (synonymous) and the URR-located SNP T340C, which did not occur in the other variants. Furthermore, four variants from vaccinated women had relatively long branches, but were not characterized by specific SNPs. The frequency of G712A in the URR was the only SNP observed to be marginally higher among vaccinated women than unvaccinated women. Non-synonymous SNPs T266A in the FG-loop of L1 and L83V in E6 were common among variants from vaccinated and unvaccinated women, but present in similar frequencies. In conclusion, the detection of HPV16 in vaccinated (and unvaccinated) women seemed to be the result of random circulation within this study population.

In silico analysis of the effects of non-synonymous SNPs associated with human GSK3B gene on its structure and function

Single Nucleotide Polymorphisms (SNPs) are abundantly identified by next generation sequencing (NGS) technology. Glycogen synthase kinase-3 beta (GSK3B), a widely expressed protein kinase, plays pivotal roles in cellular pathways. However, study on SNPs associated with GSK3B and their functional consequences is lacking. In this study, we analysed non-synonymous SNPs of GSK3B gene and their implications using computational tools. From NCBI dbSNP, 103,087 SNPs of GSK3B were initially gathered, later narrowed down to 255 unique nsSNPs. Around one-third of the nsSNPs resulted in charge and polarity change of the amino acids of the protein. 41 nsSNPs were found to significantly alter the stability of GSK3B protein (ΔΔG ≤ -1 or ≥ 1 kcal/mol) and few of them also affected the disorderness at the mutation site. Evolutionary conservation of the nsSNPs in the protein revealed 25 nsSNP may be deleterious to GSK3B protein function. Finally, 4 critical nsSNPs (Y161C, R167G, P225L and Y234D) were identified that can significantly alter both the stability and function of GSK3B. Furthermore, this study predicted 60 post-translational modification sites in GSK3B among which 26 sites contained nsSNPs. Interestingly, 7 upstream ORFs (uORFs) with high ribosomal occupancy were also detected in GSK3B mRNA that can reduce the expression of GSK3B protein. Altogether, this study has employed various in silico methods to characterize GSK3B nsSNPs, but they have limitations. These tools often overlook the cellular context, interacting partners, PTMs and the dynamic nature of proteins, which can affect protein behaviour and function. Despite these limitations, in silico tools are valuable for initial screening and prioritizing SNPs. The prioritized SNPs obtained in this study (Y161C, R167G, P225L and Y234D) should be experimentally validated using techniques like genome editing, biochemical assays, interactome analysis in cell lines and animal models to confirm their biological relevance.Clinical trial registration: Not Applicable.