Non-synonymous Single Nucleotide Polymorphisms Research Articles

Relationships between genomic dissipation and de novo SNP evolution.

Patterns of single nucleotide polymorphisms (SNPs) in eukaryotic DNA are traditionally attributed to selective pressure, drift, identity descent, or related factors-without accounting for ways in which bias during de novo SNP formation, itself, might contribute. A functional and phenotypic analysis based on evolutionary resilience of DNA points to decreased numbers of non-synonymous SNPs in human and other genomes, with a predominant component of SNP depletion in the human gene pool caused by robust preferences during de novo SNP formation (rather than selective constraint). Ramifications of these findings are broad, belie a number of concepts regarding human evolution, and point to a novel interpretation of evolving DNA across diverse species.

Identification of Breed-Specific SNPs of Danish Large White Pig in Comparison with Four Chinese Local Pig Breed Genomes.

Genetic variation facilitates the evolution, environmental adaptability, and biodiversity of organisms. Danish Large White (LW) pigs have more desirable phenotypes compared with local Chinese pigs, which have difficulty adapting to the modern swine industry. However, the genome-wide mutational differences between these pig breeds are yet to be evaluated. Therefore, this study aimed to evaluate genomic variation and identify breed-specific SNPs in Danish LW pigs. Here, 43 LW, 15 Diqing Tibetan (DQZ), and 15 Diannan small-ear (DN) pigs whose genomes were re-sequenced with 5× depth were selected. This was followed by a conjoined analysis of our previous resequencing data of 24 Anqing six-end white (AQ) and six Asian wild (SS) pigs. In total, 39,158,378 SNPs and 13,143,989 insertion-deletions were obtained in all breeds. The variation number of LW pigs was the lowest, with 287,194 breed-specific and 1289 non-synonymous SNPs compared with Chinese breeds. Functional analysis of the breed-specific non-synonymous SNPs indicated that these mutations were mainly associated with the reproductive performance, feed intake, and feed conversion ratio of LW pigs. These findings provide a theoretical basis for genetic improvements in the Chinese swine industry.

Novel polymorphisms and genetic studies of the shadow of prion protein gene (SPRN) in pheasants.

Prion diseases in mammals are caused by the structural conversion of the natural prion protein (PrPC) to a pathogenic isoform, the "scrapie form of prion protein (PrPSc)." Several studies reported that the shadow of prion protein (Sho), encoded by the shadow of prion protein gene (SPRN), is involved in prion disease development by accelerating the conformational conversion of PrPC to PrPSc. Until now, genetic polymorphisms of the SPRN gene and the protein structure of Sho related to fragility to prion disease have not been investigated in pheasants, which are a species of poultry. Here, we identified the SPRN gene sequence by polymerase chain reaction (PCR) and compared the SPRN gene and Sho protein sequences among various prion disease-susceptible and -resistant species to identify the distinctive genetic features of pheasant Sho using Clustal Omega. In addition, we investigated genetic polymorphisms of the SPRN gene in pheasants and analyzed genotype, allele, and haplotype frequencies, as well as linkage disequilibrium among the genetic polymorphisms. Furthermore, we used in silico programs, namely Mutpred2, MUpro and AMYCO, to investigate the effect of non-synonymous single nucleotide polymorphisms (SNPs). Finally, the predicted secondary and tertiary structures of Sho proteins from various species were analyzed by Alphafold2. In the present study, we reported pheasant SPRN gene sequences for the first time and identified a total of 14 novel SNPs, including 7 non-synonymous and 4 synonymous SNPs. In addition, the pheasant Sho protein sequence showed 100% identity with the chicken Sho protein sequence. Furthermore, amino acid substitutions were predicted to affect the hydrogen bond distribution in the 3D structure of the pheasant Sho protein. To the best of our knowledge, this is the first report of the genetic and structural features of the pheasant SPRN gene.

Plasmodium falciparum AMA1 and CSP antigen diversity in parasite isolates from southern Ghana.

Diversity in malarial antigens is an immune evasion mechanism that gives malaria parasites an edge over the host. Immune responses against one variant of a polymorphic antigen are usually not fully effective against other variants due to altered epitopes. This study aimed to evaluate diversity in the Plasmodium falciparum antigens apical membrane antigen 1 (PfAMA1) and circumsporozoite protein (PfCSP) from circulating parasites in a malaria-endemic community in southern Ghana and to determine the effects of polymorphisms on antibody response specificity. The study involved 300 subjects, whose P. falciparum infection status was determined by microscopy and PCR. Diversity within the two antigens was evaluated by msp2 gene typing and molecular gene sequencing, while the host plasma levels of antibodies against PfAMA1, PfCSP, and two synthetic 24mer peptides from the conserved central repeat region of PfCSP, were measured by ELISA. Of the 300 subjects, 171 (57%) had P. falciparum infection, with 165 of the 171 (96.5%) being positive for either or both of the msp2 allelic families. Gene sequencing of DNA from 55 clonally infected samples identified a total of 56 non-synonymous single nucleotide polymorphisms (SNPs) for the Pfama1 gene and these resulted in 44 polymorphic positions, including two novel positions (363 and 365). Sequencing of the Pfcsp gene from 69 clonal DNA samples identified 50 non-synonymous SNPs that resulted in 42 polymorphic positions, with half (21) of these polymorphic positions being novel. Of the measured antibodies, only anti-PfCSP antibodies varied considerably between PCR parasite-positive and parasite-negative persons. These data confirm the presence of a considerable amount of unique, previously unreported amino acid changes, especially within PfCSP. Drivers for this diversity in the Pfcsp gene do not immediately seem apparent, as immune pressure will be expected to drive a similar level of diversity in the Pfama1 gene.

Anti-bacterial compounds from Indian curry-leaf tree Murraya koenigii have potential to inhibit carbapenem-resistant Streptococcus pneumoniae

Problems consideredSpike in the mortality rates due to pneumococcal infections post-COVID-19 has raised clinical concerns globally. β-lactam antibiotic-susceptibility in the causative organism, Streptococcus pneumoniae has been compromised due to non-synonymous single nucleotide polymorphisms (nsSNPs) in different penicillin-binding protein (PBP) targets. Carbapenem-resistant S. pneumoniae strains have further constricted the pharmacologically safe therapeutic options. MethodsIn-house sequenced whole genomes of meropenem and cefotaxime resistant Indian and global isolates were used to retrieve parent and mutant proteins. Anti-bacterial phtytochemicals (ABPs) from Indian curry-leaf tree Murraya koenigii were identified using database-mining. ABPs were screeened based on pharmacokinetics, and toxicity profiles. Site-specific molecular docking studies revealed ABPs‘ affinity towards the targets as compared to control drug meropenem. Coarse-grained and normal-mode molecular dynamics simulation further asserted the interactive stability of the chosen compounds. ResultsFrom the whole genomes (n = 28) of meropenem and cefotaxime resistant (MIC ≥0.5 μg/ml) Indian and global isolates, ∼50 nsSNPs in the active-site regions of PBP1A, PBP2B and PBP2X were identified. Based on pharmacological safety and high drug-likeliness, 10 ABPs were funnelled down from 226 ABPs in M. koenigii. A pyrano-carbazole alkaloid koenine was further prioritized as the most efficient candidate based on consistently high binding affinity (binding energies −7.7 to −8.9 kcal/mol), low inhibition constant (0.69–3.5 μM) and low structural fluctuations (<1.1 Å), against mutant targets as compared to meropenem. ConclusionThe study encourages experimental validation of koenine against pneumococci and also promotes the use of similar phytochemical resources for evaluation against emerging drug-resistant pathogens harbouring mutant targets.

Polymorphism of spore germination receptor genes in Bacillus anthracis strains of major genetic lineages

Introduction. The genetic structure of the global population of Bacillus anthracis characterized by an unequal distribution of isolates of the main genetic lineages A, B and C, the reason for which is unknown. Determining the characteristics of genes encoding factors that determine the existence of this pathogen at the intra- and extra organismal stages of the life cycle, which can influence the prevalence of strains, is relevant. Aim — сharacterization of the genes and proteins of spore germination in strains of the anthrax pathogen of different genetic lineages. Materials and methods. Whole genome sequences of 46 B. anthracis strains and the CI strain of B. cereus biovar anthracis from the NCBI GenBank database studied. In silico analysis carried out using the programs “BLASTn”, “MEGA X”, “Tandem Repeat Finder”. Results. The number of SNPs, indels and pseudogenes in B. anthracis strains of line B was 2,7–25,6 times higher, in line C 2–23,5 times, and in the B. cereus biovar anthracis strain was 20–2841 times higher than in strains of line A. Significant substitutions in genes leading to changes in the amino acid composition of 10 germination receptor proteins were also significantly more common in B. anthracis strains of lines B, C and the B. cereus biovar anthracis strain. Undescribed VNTRs within the gerHA gene with repeat units of 78 and 117 bp identified, varying between and within isolates of different genetic lineages. Six germination receptor genes have been shown to have rare starting codons. Conclusion. A larger number of non-synonymous SNPs in the genes of spore germination receptors with changes in the amino acid composition of proteins in B. anthracis strains of the main genetic lines B, C and B. cereus biovar anthracis than in strains of line A suggests their limited adaptive capabilities and may be one of the explanations for the lower prevalence compared to line A. Differences in the gerHA and gerM genes make it possible to differentiate the major B and C genetic lineages from A.

The implications of APOBEC3-mediated C-to-U RNA editing for human disease

Intra-organism biodiversity is thought to arise from epigenetic modification of constituent genes and post-translational modifications of translated proteins. Here, we show that post-transcriptional modifications, like RNA editing, may also contribute. RNA editing enzymes APOBEC3A and APOBEC3G catalyze the deamination of cytosine to uracil. RNAsee (RNA site editing evaluation) is a computational tool developed to predict the cytosines edited by these enzymes. We find that 4.5% of non-synonymous DNA single nucleotide polymorphisms that result in cytosine to uracil changes in RNA are probable sites for APOBEC3A/G RNA editing; the variant proteins created by such polymorphisms may also result from transient RNA editing. These polymorphisms are associated with over 20% of Medical Subject Headings across ten categories of disease, including nutritional and metabolic, neoplastic, cardiovascular, and nervous system diseases. Because RNA editing is transient and not organism-wide, future work is necessary to confirm the extent and effects of such editing in humans.

Identification and analysis of oncogenic non-synonymous single nucleotide polymorphisms in the human NRAS gene: An exclusive in silico study

BackgroundN-ras protein is encoded by the NRAS gene and operates as GDP-GTP-controlled on/off switching. N-ras interacts with cellular signaling networks that regulate various cellular activities including cell proliferation and survival. The nonsynonymous single nucleotide polymorphism (nsSNPs)-mediated alteration can substantially disrupt the structure and activity of the corresponding protein. N-ras has been reported to be associated with numerous diseases including cancers due to the nsSNPs. A comprehensive study on the NRAS gene to unveil the potentially damaging and oncogenic nsSNPs is yet to be accomplished. Hence, this extensive in silico study is intended to identify the disease-associated, specifically oncogenic nsSNPs of the NRAS gene. ResultsOut of 140 missense variants, 7 nsSNPs (I55R, G60E, G60R, Y64D, L79F, D119G, and V152F) were identified to be damaging utilizing 10 computational tools that works based on different algorithms with high accuracy. Among those, G60E, G60R, and D119G variants were further filtered considering their location in the highly conserved region and later identified as oncogenic variants. Interestingly, G60E and G60R variants were revealed to be particularly associated with lung adenocarcinoma, rhabdomyosarcoma, and prostate adenocarcinoma. Therefore, D119G could be subjected to detailed investigation for identifying its association with specific cancer. ConclusionThis in silico study identified the deleterious and oncogenic missense variants of the human NRAS gene that could be utilized for designing further experimental investigation. The outcomes of this study would be worthwhile in future research for developing personalized medicine.

Effects of nonsynonymous single nucleotide polymorphisms of the KIAA1217, SNTA1 and LTBP1 genes on the growth traits of Ujumqin sheep.

Sheep body size can directly reflect the growth rates and fattening rates of sheep and is also an important index for measuring the growth performance of meat sheep. In this study, high-resolution resequencing data from four sheep breeds (Dorper sheep, Suffolk sheep, Ouessant sheep, and Shetland sheep) were analyzed. The nonsynonymous single nucleotide polymorphisms of three candidate genes (KIAA1217, SNTA1, and LTBP1) were also genotyped in 642 healthy Ujumqin sheep using MALDI-TOFMS and the genotyping results were associated with growth traits. The results showed that different genotypes of the KIAA1217 g.24429511T>C locus had significant effects on the chest circumferences of Ujumqin sheep. The SNTA1 g.62222626C>A locus had different effects on the chest depths, shoulder widths and rump widths of Ujumqin sheep. This study showed that these two sites can be used for marker-assisted selection, which will be beneficial for future precision molecular breeding.

Strong Homology Between Colonizing and Bloodstream Carbapenem-Resistant Acinetobacter Spp.: Implications for Empiric Antibiotic Therapy in Hematological Patients.

This study aimed to assess the impact of colonization status on the outcomes of Acinetobacter spp. bloodstream infection (BSI) and investigate the homology and within-host evolution between colonizing and bloodstream carbapenem-resistant Acinetobacter spp. (CRA) to inform antibiotic therapeutic decisions. We analyzed clinical outcomes of 46 hematological patients with Acinetobacter spp. BSI and performed whole-genome sequencing on the remaining CRA isolates. Among the patients, 39.1% (n=18) had prior Acinetobacter spp. colonization. Colonized patients had higher rates of polymicrobial BSI (50.0% vs 21.4%, P=0.044) and CRA BSI (72.2% vs 17.9%, P<0.001), resulting in elevated inflammatory markers and increased 30-day mortality. Each of the eight pairs of the remaining respiratory colonizing and bloodstream CRA strains belonged to the same genomospecies. Each pair exhibited definitive agreement in at least 21 of the 22 most representative antibiotic susceptibility tests. The minimum spanning tree based on multilocus sequence typing (MLST) and phylogenetic trees based on MLST and single nucleotide polymorphism (SNP) all indicated that each pair shared the same minimum branch. Very few non-synonymous SNPs in genic regions were identified during the transition from respiratory colonization to bloodstream infection, with minimal changes in virulence genes. Homology analysis suggested that CRA BSI originated from colonizing isolates in the respiratory tract. Strict infection control measures are needed to manage Acinetobacter spp. colonisation in hematological patients. Appropriate empirical therapy can be administered for suspected CRA BSI based on the antimicrobial minimum inhibitory concentration of CRA colonising the respiratory tract.

Characterization of missense nonsynonymous single-nucleotide polymorphism of runt-related transcription factor-2 gene - An in silico approach.

Single-nucleotide polymorphism (SNP) codes for multiple amino acids, impacting protein functions and disease prognosis. Runt-related transcription factor-2 (RUNX2), a transcription factor linked to osteoblast differentiation, regulates cell proliferation in endothelium and osteoblastic cells. Understanding Runx2's role in nonosseous tissues is rapidly advancing. This study aims to identify harmful SNPs of the RUNX2 gene that may alter disease susceptibility using computational techniques. The study uses various in silico methods to identify nonsynonymous SNPs (nsSNPs) of the RUNX2 gene, which could potentially alter protein structure and functions, with further analyses by I-Mutant, ConSurf, Netsurf 3.0, GeneMANIA, and Have (y)Our Protein Explained. Six missense nsSNPs were identified as potentially harmful, disease-causing, and damaging. Four were found to be unstable, while five were conserved. All six nsSNPs had a coiled secondary structure. Five nsSNPs were found to be destabilized. The RUNX2 gene's deleterious missense nsSNPs were identified by this study, and they may be exploited in future experimental studies. These high-risk nsSNPs might be considered target molecules in therapeutic and diagnostic therapies in teeth and bone development.

Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries.

We develop a method, SBayesRC, that integrates genome-wide association study (GWAS) summary statistics with functional genomic annotations to improve polygenic prediction of complex traits. Our method is scalable to whole-genome variant analysis and refines signals from functional annotations by allowing them to affect both causal variant probability and causal effect distribution. We analyze 50 complex traits and diseases using ∼7 million common single-nucleotide polymorphisms (SNPs) and 96 annotations. SBayesRC improves prediction accuracy by 14% in European ancestry and up to 34% in cross-ancestry prediction compared to the baseline method SBayesR, which does not use annotations, and outperforms other methods, including LDpred2, LDpred-funct, MegaPRS, PolyPred-S and PRS-CSx. Investigation of factors affecting prediction accuracy identifies a significant interaction between SNP density and annotation information, suggesting whole-genome sequence variants with annotations may further improve prediction. Functional partitioning analysis highlights a major contribution of evolutionary constrained regions to prediction accuracy and the largest per-SNP contribution from nonsynonymous SNPs.

Human P2X7 receptor variants Gly150Arg and Arg276His polymorphisms have differential effects on risk association and cellular functions in pancreatic cancer

BackgroundThe purinergic P2X7 receptor (P2X7R) plays an important role in the crosstalk between pancreatic stellate cells (PSCs) and cancer cells, thus promoting progression of pancreatic ductal adenocarcinoma (PDAC). Single nucleotide polymorphisms (SNPs) in the P2X7R have been reported for several cancers, but have not been explored in PDAC.Materials and methodsBlood samples from PDAC patients and controls were genotyped for 11 non-synonymous SNPs in P2X7R and a risk analysis was performed. Relevant P2X7R-SNP GFP variants were expressed in PSCs and cancer cells and their function was assayed in the following tests. Responses in Ca2+ were studied with Fura-2 and dye uptake with YO-PRO-1. Cell migration was monitored by fluorescence microscopy. Released cytokines were measured with MSD assay.ResultsRisk analysis showed that two SNPs 474G>A and 853G>A (rs28360447, rs7958316), that lead to the Gly150Arg and Arg276His variants, had a significant but opposite risk association with PDAC development, protecting against and predisposing to the disease, respectively. In vitro experiments performed on cancer cells and PSCs expressing the Gly150Arg variant showed reduced intracellular Ca2+ response, fluorescent dye uptake, and cell migration, while the Arg276His variant reduced dye uptake but displayed WT-like Ca2+ responses. As predicted, P2X7R was involved in cytokine release (IL-6, IL-1β, IL-8, TNF-α), but the P2X7R inhibitors displayed varied effects.ConclusionIn conclusion, we provide evidence for the P2X7R SNPs association with PDAC and propose that they could be considered as potential biomarkers.

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.

Computational Analysis of Deleterious nsSNPs in INS Gene Associated with Permanent Neonatal Diabetes Mellitus.

Insulin gene mutations affect the structure of insulin and are considered a leading cause of neonatal diabetes and permanent neonatal diabetes mellitus PNDM. These mutations can affect the production and secretion of insulin, resulting in inadequate insulin levels and subsequent hyperglycemia. Early discovery or prediction of PNDM can aid in better management and treatment. The current study identified potential deleterious non-synonymous single nucleotide polymorphisms nsSNPs in the INS gene. The analysis of the nsSNPs in the INS gene was conducted using bioinformatics tools by implementing computational algorithms including SIFT, PolyPhen2, SNAP2, SNPs & GO, PhD-SNP, MutPred2, I-Mutant, MuPro, and HOPE tools to investigate the prediction of the potential association between nsSNPs in the INS gene and PNDM. Three mutations, C96Y, P52R, and C96R, were shown to potentially reduce the stability and function of the INS protein. These mutants were subjected to MDSs for structural analysis. Results suggested that these three potential pathogenic mutations may affect the stability and functionality of the insulin protein encoded by the INS gene. Therefore, these changes may influence the development of PNDM. Further researches are required to fully understand the various effects of mutations in the INS gene on insulin synthesis and function. These data can aid in genetic testing for PNDM to evaluate its risk and create treatment and prevention strategies in personalized medicine.

Prediction of deleterious non-synonymous SNPs of human MDC1 gene: an in silico approach

MDC1 (Mediator of DNA damage Checkpoint protein 1) functions to facilitate the localization of numerous DNA damage response (DDR) components to DNA double-strand break sites. MDC1 is an integral component in preserving genomic stability and appropriate DDR regulation. There haven’t been systematic investigations of MDC1 mutations that induce cancer and genomic instability. Variations in nsSNPs have the potential to modify the protein chemistry and their function. Describing functional SNPs in disease-associated genes presents a significant conundrum for investigators, it is possible to assess potential functional SNPs before conducting larger population examinations. Multiple sequences and structure-based bioinformatics strategies were implemented in the current in-silico investigation to discern potential nsSNPs of the MDC1 genes. The nsSNPs were identified with SIFT, SNAP2, Align GVGD, PolyPhen-2, and PANTHER, and their stability was determined with MUpro. The conservation, solvent accessibility, and structural effects of the mutations were identified with ConSurf, NetSurfP-2.0, and SAAFEC-SEQ respectively. Cancer-related analysis of the nsSNPs was conducted using cBioPortal and TCGA web servers. The present study appraised five nsSNPs (P1426T, P69S, P194R, P203L, and H131Y) as probably mutilating due to their existence in highly conserved regions and propensity to deplete protein stability. The nsSNPs P194R, P203L, and H131Y were concluded as deleterious and possibly damaging from the 5 prediction tools. The functional nsSNP P194R mutation is associated with skin cutaneous melanoma while no significant records were found for other nsSNPs. The present study concludes that the highly deleterious P194R mutations can potentially induce genomic instability and contribute to various cancers’ pathogenesis. Developing drugs targeting these mutations can undoubtedly be advantageous in large population-based studies, particularly in the development of precision medicine.

Comprehensive in silico analysis of prolactin receptor (PRLR) gene nonsynonymous single nucleotide polymorphisms (nsSNPs) reveals multifaceted impact on protein structure, function, and interactions

This study delves into the functional and structural implications of non-synonymous single nucleotide polymorphisms (nsSNPs) within the Prolactin Receptor (PRLR) gene. Thirteen deleterious nsSNPs were identified through bioinformatics tools, with SIFT predicting 168 out of 395 nsSNPs as detrimental, exhibiting tolerance index (TI) scores ranging from 0 to 0.05. Polyphen2 assigned likelihood scores >0.99 to all 13 nsSNPs, indicating high probability of harm, while Panther scores classified most nsSNPs as ‘probably damaging’, with specific mutations like W218R scoring 0.74, suggesting a higher impact. Stability analysis using DDG I-Mutant and DDG Mupro consistently predicted decreased stability for all mutations, with CUPSAT indicating mutations like V125G and W218R significantly decreasing stability. Structural analysis through DynaMut predicted destabilization for all mutations except L196I and L292H. MutPred2 highlighted structural alterations for all nsSNPs except L196I, L293V, R315W, and S353N. Domain analysis revealed key mutations within essential functional domains, with five nsSNPs located within Fibronectin type-III domains. Bayesian analysis through ConSurf identified 9 critical residues, with 11 nsSNPs exhibiting notably high conservation. STRING analysis unveiled a complex interaction network, indicating involvement in vital biological processes like lactation. Molecular dynamics (MD) simulations, spanning 100 nanoseconds, elucidated structural dynamics induced by detrimental missense SNPs. Post-translational modification (PTM) analysis identified specific mutations, such as R351, involved in methylation, while S353 was implicated in phosphorylation and glycosylation. These findings offer comprehensive insights into the molecular and phenotypic effects of deleterious nsSNPs in the PRLR gene, crucial for selective breeding. Communicated by Ramaswamy H. Sarma

A Comprehensive &lt;i&gt;in Silico&lt;/i&gt; Analysis of the Functional and Structural Consequences of the Deleterious Missense Nonsynonymous SNPs in Human &lt;i&gt;GABRA6&lt;/i&gt; Gene

Epilepsy, a prevalent neurological disorder, affects more than 50 million individuals worldwide and is characterized by recurring seizures. Nonsynonymous single nucleotide polymorphisms (nsSNPs) found within coding regions of epilepsy-related genes are believed to have significant impacts on protein function. This is due to their tendency to cause mutations in the encoded amino acids, which can subsequently lead to pathogenic alterations in protein structure and function. Consequently, nsSNPs have the potential to serve as diagnostic markers for epilepsy and other neuropsychiatric conditions. The primary objective of this study is to evaluate the harmful effects of missense nsSNP mutations on the &lt;em&gt;GABRA6&lt;/em&gt; gene. The &lt;em&gt;GABRA6&lt;/em&gt; gene encodes the alpha-6 subunit of the GABA&lt;sub&gt;A&lt;/sub&gt; receptor, and previous research showed one case substitution mutation in the GABRA6 gene is associated with childhood absence epilepsy (CAE) and atonic seizures. To achieve this, we employed various computational tools, including SIFT, PolyPhen-2, PROVEAN, Condel, SNPs &amp; GO, PMut, SNAP2, MutPred2, and SNPeffect4.0, for predicting missense nsSNPs. Additionally, we used I-Mutant3.0 and MUpro to analyze protein stability, ConSurf to assess evolutionary conservation, FTSite and COACH to predict ligand binding sites, SOPMA and PSIPRED to analyze protein secondary structures, project HOPE to predict structural changes, and I-TASSER to model the 3D structure. Furthermore, structural validation was conducted using the PROCHECK and ERRAT servers. At the same time, molecular dynamics simulations were performed using GROMACS to gain a better understanding of the effects of mutations on protein structure. Among the 451 missense nsSNPs identified within the &lt;em&gt;GABRA6&lt;/em&gt; gene, three were found to have pathogenic effects on the structure and function of the protein, potentially, there may be a contribution to the development of seizures or other neuropsychiatric disorders. Notably, two of these missense nsSNPs (&lt;em&gt;W87S&lt;/em&gt; and &lt;em&gt;W112R&lt;/em&gt;) were located within the ligand-binding domain, while the third (&lt;em&gt;C310R&lt;/em&gt;) was situated in the transmembrane domain. It is crucial to acknowledge that despite their predicted pathogenicity, these variants are currently classified as of uncertain significance in clinical and genomic databases worldwide due to the lack of correlation with epilepsy in empirical studies. Without experimental data to validate these predictions, caution is warranted in interpreting the findings.

A gain-of-function mutation in zinc cluster transcription factor Rob1 drives Candida albicans adaptive growth in the cystic fibrosis lung environment.

Candida albicans chronically colonizes the respiratory tract of patients with Cystic Fibrosis (CF). It competes with CF-associated pathogens (e.g. Pseudomonas aeruginosa) and contributes to disease severity. We hypothesize that C. albicans undergoes specific adaptation mechanisms that explain its persistence in the CF lung environment. To identify the underlying genetic and phenotypic determinants, we serially recovered 146 C. albicans clinical isolates over a period of 30 months from the sputum of 25 antifungal-naive CF patients. Multilocus sequence typing analyses revealed that most patients were individually colonized with genetically close strains, facilitating comparative analyses between serial isolates. We strikingly observed differential ability to filament and form monospecies and dual-species biofilms with P. aeruginosa among 18 serial isolates sharing the same diploid sequence type, recovered within one year from a pediatric patient. Whole genome sequencing revealed that their genomes were highly heterozygous and similar to each other, displaying a highly clonal subpopulation structure. Data mining identified 34 non-synonymous heterozygous SNPs in 19 open reading frames differentiating the hyperfilamentous and strong biofilm-former strains from the remaining isolates. Among these, we detected a glycine-to-glutamate substitution at position 299 (G299E) in the deduced amino acid sequence of the zinc cluster transcription factor ROB1 (ROB1G299E), encoding a major regulator of filamentous growth and biofilm formation. Introduction of the G299E heterozygous mutation in a co-isolated weak biofilm-former CF strain was sufficient to confer hyperfilamentous growth, increased expression of hyphal-specific genes, increased monospecies biofilm formation and increased survival in dual-species biofilms formed with P. aeruginosa, indicating that ROB1G299E is a gain-of-function mutation. Disruption of ROB1 in a hyperfilamentous isolate carrying the ROB1G299E allele abolished hyperfilamentation and biofilm formation. Our study links a single heterozygous mutation to the ability of C. albicans to better survive during the interaction with other CF-associated microbes and illuminates how adaptive traits emerge in microbial pathogens to persistently colonize and/or infect the CF-patient airways.

Allele mining for Granule Bound Starch Synthase1 (GBSS1) gene governing amylose content in aromatic rice (Oryza sativa L.) germplasm

The current research investigated sequence polymorphism and allelic variation in the GBSS1 gene responsible for determining the apparent amylose content (AAC) in an aromatic rice germplasm panel. Analysis of 271 germplasm accessions from India revealed a wide range of variation in AAC, from 4.81% (Manipur Black Rice) to 35.88% (IGSR -3-1-5). Re-sequencing and sequence analysis of a representative subset of 28 accessions unveiled a nucleotide diversity (π) of 0.00374, identifying 72 sequence variants. Two novel non-synonymous SNPs were identified: a G/A transition in exon 2, resulting in the substitution of glutamic acid at the 91st position with lysine, and a C/T transition in exon 11, substituting alanine at the 468thposition with valine. Additionally, two previously reported SNPs were validated: A/C in exon 6 (Tyr224Ser) and C/T in exon 10 (Pro415Ser). Furthermore, a 23 bp exon 2 duplication was found in Manipur Black Rice (4.81%), Bhogali Bora (6.20%), and NJ-72 (7.36%), indicating the prevalence of this allele in the very low AAC category (3–9%) in addition to the glutinous or waxy rice category with 0–2% AAC. These findings are significant for manipulating starch content in aromatic rice germplasm.