Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

Single nucleotide polymorphism of Notch1 gene rs3124599 allele is associated with the severity of CVA6-related HFMD in the Chinese Han population

BACKGROUND

There is evidence suggesting that Notch1 signaling pathway contributes to the development of hand, foot, and mouth disease (HFMD); however, the role of Notch1 gene polymorphisms in the severity of coxsackievirus A6 (CVA6)-related HFMD remains unclear. This study aimed to investigate the correlation between Notch1 gene polymorphisms and the severity of CVA6-related HFMD.

METHODS

A total of 196 patients (Chinese Han population) diagnosed with CVA6-related HFMD through nucleic acid testing were included in this study. Among them, 97 patients were classified as severe cases, while 99 cases were categorized as mild. The mRNA levels of Notch1 in the peripheral blood leukocytes of HFMD patients were detected by quantitative real-time polymerase chain reaction (qRT-PCR), and the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique was utilized for genotyping of rs3124599, rs3124603, and rs3124591.

RESULTS

The frequencies of rs3124599 alleles were G (39.0%) and A (61.0%), while the frequencies of rs3124599 genotypes were GG (12.2%), GA (53.6%), and AA (34.2%), respectively. In the recessive model, the frequency of rs3124599 AA genotypes significantly increased in severe patients, compared to mild patients (P < 0.05). Due to the low frequency of alleles for rs3124591 and rs3124603 in patients, as well as the absence of any difference in their distribution between the two groups (P > 0.05), no additional statistical analysis was performed. After adjusting for age and sex, patients with rs3124599 AA genotype had a significantly higher risk of severe HFMD in comparison to G allele carriers (GA/GG), with an odds ratio (95% confidence interval) of 2.010 (1.094, 3.691). Meanwhile, the mRNA levels of Notch1 were found to be significantly higher in severe patients compared to mild patients (P < 0.05), and a positive correlation was observed between Notch1 mRNA levels and the peripheral blood monocyte count (r = 0.42, P < 0.001). Additionally, there were significant differences observed in Notch1 mRNA levels and peripheral blood monocyte counts between patients with the AA genotype of rs3124599 and those with the GA genotype or G allele carriers (P < 0.05).

CONCLUSION

In the Chinese Han population, there is a strong correlation between the Notch1 rs3124599 allele and the severity of CVA6-related HFMD. This correlation may be attributed to genetic polymorphism of rs3124599 regulating Notch1 transcription levels. These findings reveal the important role of Notch1 gene polymorphism in CVA6 infection, establishing a scientific foundation for the precise control of severe HFMD.

MotifbreakR v2: extended capability and database integration

MotifbreakR is a software tool that scans genetic variants against position weight matrices of transcription factors (TF) to determine the potential for the disruption of TF binding at the site of the variant. It leverages the Bioconductor suite of software packages and annotations to operate across a diverse array of genomes and motif databases. Initially developed to interrogate the effect of single nucleotide variants (common and rare SNVs) on potential TF binding sites, in motifbreakR v2, we have updated the functionality. New features include the ability to query other types of more complex genetic variants, such as short insertions and deletions (indels). This function allows modeling a more extensive array of variants that may have more significant effects on TF binding. Additionally, while TF binding is based partly on sequence preference, predictions of TF binding based on sequence preference alone can indicate many more potential binding events than observed. Adding information from DNA-binding sequencing datasets lends confidence to motif disruption prediction by demonstrating TF binding in cell lines and tissue types. Therefore, motifbreakR implements querying the ReMap2022 database for evidence that a TF matching the disrupted motif binds over the disrupting variant. Finally, in motifbreakR, in addition to the existing interface, we have implemented an R/Shiny graphical user interface to simplify and enhance access to researchers with different skill sets.

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

Methods for Functional Characterization of Genetic Polymorphisms of Non-Coding Regulatory Regions of the Human Genome

Currently, numerous associations between genetic polymorphisms and various diseases have been characterized through the Genome-Wide Association Studies. Majority of the clinically significant polymorphisms are localized in non-coding regions of the genome. While modern bioinformatic resources make it possible to predict molecular mechanisms that explain influence of the non-coding polymorphisms on gene expression, such hypotheses require experimental verification. This review discusses the methods for elucidating molecular mechanisms underlying dependence of the disease pathogenesis on specific genetic variants within the non-coding sequences. A particular focus is on the methods for identification of transcription factors with binding efficiency dependent on polymorphic variations. Despite remarkable progress in bioinformatic resources enabling prediction of the impact of polymorphisms on the disease pathogenesis, there is still the need for experimental approaches to investigate this issue.

Genetic analysis of PALB2 gene WD40 domain in canine mammary tumour patients

BACKGROUND

DNA repair mechanisms are essential for tumorigenesis and disruption of HR mechanism is an important predisposing factor of human breast cancers (BC). PALB2 is an important part of the HR. There are similarities between canine mammary tumours (CMT) and BCs. As its human counterpart, PALB2 mutations could be a predisposing factor of CMT.

OBJECTIVES

In this study, we aimed to investigate the impacts of PALB2 variants on tumorigenesis and canine mammary tumor (CMT) malignancy.

METHODS

We performed Sanger sequencing to detect germline mutations in the WD40 domain of the canine PALB2 gene in CMT patients. We conducted in silico analysis to investigate the variants, and compared the germline PALB2 mutations in humans that cause breast cancer (BC) with the variants detected in dogs with CMT.

RESULTS

We identified an intronic (c.3096+8C>G) variant, two exonic (p.A1050V and p.R1354R) variants, and a 3' UTR variant (c.4071T>C). Of these, p.R1354R and c.4071T>C novel variants were identified for the first time in this study. We found that the p.A1050V mutation had a significant effect. However, we could not determine sufficient similarity due to the differences in nucleotide/amino acid sequences between two species. Nonetheless, possible variants of human sequences in the exact location as their dog counterparts are associated with several cancer types, implying that the variants could be crucial for tumorigenesis in dogs. Our results did not show any effect of the variants on tumor malignancy.

CONCLUSIONS

The current project is the first study investigating the relationship between the PALB2 gene WD40 domain and CMTs. Our findings will contribute to a better understanding of the pathogenic mechanism of the PALB2 gene in CMTs. In humans, variant positions in canines have been linked to cancer-related phenotypes such as familial BC, endometrial tumor, and hereditary cancer predisposition syndrome. The results of bioinformatics analyses should be investigated through functional tests or case-control studies.

Complex regulatory networks influence pluripotent cell state transitions in human iPSCs

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

Parkinson's disease risk enhancers in microglia

Genome-wide association studies have identified thousands of single nucleotide polymorphisms that associate with increased risk for Parkinson's disease (PD), but the functions of most of them are unknown. Using assay for transposase-accessible chromatin (ATAC) and H3K27ac chromatin immunoprecipitation (ChIP) sequencing data, we identified 73 regulatory elements in microglia that overlap PD risk SNPs. To determine the target genes of a "risk enhancer" within intron two of SNCA, we used CRISPR-Cas9 to delete the open chromatin region where two PD risk SNPs reside. The loss of the enhancer led to reduced expression of multiple genes including SNCA and the adjacent gene MMRN1. It also led to expression changes of genes involved in glucose metabolism, a process that is known to be altered in PD patients. Our work expands the role of SNCA in PD and provides a connection between PD-associated genetic variants and underlying biology that points to a risk mechanism in microglia.

Identification of the SIRT1 gene's most harmful non-synonymous SNPs and their effects on functional and structural features-an in silico analysis

Introduction

The sirtuin (Silent mating type information regulation 2 homolog)1(SIRT1) protein plays a vital role in many disorders such as diabetes, cancer, obesity, inflammation, and neurodegenerative and cardiovascular diseases. The objective of this in silico analysis of SIRT1's functional single nucleotide polymorphisms (SNPs) was to gain valuable insight into the harmful effects of non-synonymous SNPs (nsSNPs) on the protein. The objective of the study was to use bioinformatics methods to investigate the genetic variations and modifications that may have an impact on the SIRT1 gene's expression and function.

Methods

nsSNPs of SIRT1 protein were collected from the dbSNP site, from its three (3) different protein accession IDs. These were then fed to various bioinformatic tools such as SIFT, Provean, and I- Mutant to find the most deleterious ones. Functional and structural effects were examined using the HOPE server and I-Tasser. Gene interactions were predicted by STRING software. The SIFT, Provean, and I-Mutant tools detected the most deleterious three nsSNPs (rs769519031, rs778184510, and rs199983221).

Results

Out of 252 nsSNPs, SIFT analysis showed that 94 were deleterious, Provean listed 67 dangerous, and I-Mutant found 58 nsSNPs resulting in lowered stability of proteins. HOPE modelling of rs199983221 and rs769519031 suggested reduced hydrophobicity due to Ile 4Thr and Ile223Ser resulting in decreased hydrophobic interactions. In contrast, on modelling rs778184510, the mutant protein had a higher hydrophobicity than the wild type.

Conclusions

Our study reports that three nsSNPs (D357A, I223S, I4T) are the most damaging mutations of the SIRT1 gene. Mutations may result in altered protein structure and functions. Such altered protein may be the basis for various disorders. Our findings may be a crucial guide in establishing the pathogenesis of various disorders.

NR3C1 and NR3C2 Genes Increase the Risk of Suicide Attempt in Psychiatric Disorder Patients with History of Childhood Trauma

Background

Hypothalamic-pituitary-adrenal axis gene variants and childhood trauma (CT) are considered risk factors for suicide attempt (SA). The aim of the present study was analyzed gene x environment (GxE) interaction of NR3C1, NR3C2, and CT, and NR3C1 and NR3C2 gene expression in the development of SA with CT.

Participants and Methods

A total of 516 psychiatric Mexican patients from Instituto Nacional de Psiquiatría Ramón de la Fuente Muñiz. Among them, 274 had SA at least once and 242 had not SA. Genetic variants of NR3C1 and NR3C2 were genotyped in all the patients, of which were obtained the CT information from medical records. Additionally, the gene expression of NR3C1 and NR3C2 was also analyzed for a subsample of 96 patients, obtaining the TC information from Childhood Trauma Questionnaire (CTQ).

Results

The analysis showed a GxE interaction of NR3C1, NR3C2, and CT (OR=2.8, 95% CI [1.9-3.9], p<0.0001). Interactions were also observed with neglect (OR=2.1, 95% CI [1.4-3.1], p<0.0001), emotional abuse (OR=2.1, 95% CI [1.5-3], p<0.0001), and sexual abuse (OR=2.4, 95% CI [1.4-2.9], p<0.0001) in the prediction of SA. The analysis of gene expression identified an overexpression of NR3C1 in SA patients with high scores for physical and sexual abuse (p<0.0001; p<0.0006, respectively) and emotional neglect (p=0.014). An underexpression was observed of NR3C2, associated with high scores of trauma subtypes (p<0.0001) except physical neglect. Additionally, we observed an overexpression of NR3C1 gene in patients with SA carriers of A allele of rs6191 (p=0.0015). Also, overexpression of NR3C1 gene in carriers of G allele of rs6198 and underexpression of NR3C2 gene in carriers of G allele of rs5522 (p<0.0001).

Conclusion

Our findings suggest that genetic variants of NR3C1 and NR3C2 differentially affect expression levels, increasing the susceptibility to SA in psychiatric patients with a history of CT.

Analysis of regulatory network modules in hundreds of human stem cell lines reveals complex epigenetic and genetic factors contribute to pluripotency state differences between subpopulations

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discovered 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which were highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlied the coordinated expression of genes in the GNMs. Epigenetic analyses revealed that regulatory networks underlying self-renewal and pluripotency have a surprising level of complexity. Genetic analyses identified thousands of regulatory variants that overlapped predicted transcription factor binding sites and were associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network were significantly enriched for regulatory variants with large effects, suggesting that they may play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work captures the coordinated activity of tens of thousands of regulatory elements in hiPSCs and bins these elements into discrete functionally characterized regulatory networks, shows that regulatory elements in pluripotency networks harbor variants with large effects, and provides a rich resource for future pluripotent stem cell research.

N -acetyltransferase 2 haplotype modifies risks for both dyslipidemia and urinary bladder cancer

A novel haplotype in N -acetyltransferase 2 ( NAT2 ) composed of seven non-coding variants (rs1495741, rs4921913, rs4921914, rs4921915, rs146812806, rs35246381, and rs35570672) has been linked to dyslipidemia by multiple, independent genome-wide association studies. The haplotype is located approximately 14 kb downstream of NAT2-coding region (ch8:18,272,377-18,272,881; GRCh38/hg38) and represents a non-coding, intergenic haplotype. Interestingly, the same dyslipidemia NAT2 haplotype is also linked to urinary bladder cancer risk. Dyslipidemia risk alleles are associated with rapid acetylator phenotype, whereas bladder cancer risk alleles are associated with slow acetylator, suggesting that the level of systemic NAT2 activity modifies the risk of these pathologies. We speculate that rs1495741 (and its associated haplotype) belongs to a distal regulatory element of human NAT2 gene (e.g., enhancer or silencer), and the genetic variation at the newly discovered haplotype results in a differential level of NAT2 gene expression. Understanding how this NAT2 haplotype contributes to not only urinary bladder cancer but also to dyslipidemia will ultimately help devise strategies to identify and protect susceptible individuals.

Arabidopsis Transcription Regulatory Factor Domain/Domain Interaction Analysis Tool-Liquid/Liquid Phase Separation, Oligomerization, GO Analysis: A Toolkit for Interaction Data-Based Domain Analysis

Although a large number of databases are available for regulatory elements, a bottleneck has been created by the lack of bioinformatics tools to predict the interaction modes of regulatory elements. To reduce this gap, we developed the Arabidopsis Transcription Regulatory Factor Domain/Domain Interaction Analysis Tool-liquid/liquid phase separation (LLPS), oligomerization, GO analysis (ART FOUNDATION-LOG), a useful toolkit for protein-nucleic acid interaction (PNI) and protein-protein interaction (PPI) analysis based on domain-domain interactions (DDIs). LLPS, protein oligomerization, the structural properties of protein domains, and protein modifications are major components in the orchestration of the spatiotemporal dynamics of PPIs and PNIs. Our goal is to integrate PPI/PNI information into the development of a prediction model for identifying important genetic variants in peaches. Our program unified interdatabase relational keys based on protein domains to facilitate inference from the model species. A key advantage of this program lies in the integrated information of related features, such as protein oligomerization, LOG analysis, structural characterizations of domains (e.g., domain linkers, intrinsically disordered regions, DDIs, domain-motif (peptide) interactions, beta sheets, and transmembrane helices), and post-translational modification. We provided simple tests to demonstrate how to use this program, which can be applied to other eukaryotic organisms.

Exploration of the association between the single-nucleotide polymorphism of co-stimulatory system and rheumatoid arthritis

Introduction

The human leukocyte antigen (HLA) has been linked to the majority of autoimmune diseases (ADs). However, non-HLA genes may be risk factors for ADs. A number of genes encoding proteins involved in regulating T-cell and B-cell function have been identified as rheumatoid arthritis (RA) susceptibility genes.

Methods

In this study, we investigated the association between RA and single-nucleotide polymorphisms (SNPs) of co-stimulatory or co-inhibitory molecules in 124 RA cases and 100 healthy controls without immune-related diseases [including tumor necrosis factor superfamily member 4 (TNFSF4), CD28, cytotoxic T-lymphocyte-associated protein 4 (CTLA4), and programmed cell death protein 1 (PDCD1)].

Results

The results showed that there were 13 SNPs associated with RA, including rs181758110 of TNFSF4 (CC vs. CT, p = 0.038); rs3181096 of CD28 (TT vs. CC + CT, p = 0.035; CC vs. TT, p = 0.047); rs11571315 (TT vs. CT, p = 0.045), rs733618 (CC vs. TT + CT, p = 0.043), rs4553808 (AA vs. AG vs. GG, p = 0.035), rs11571316 (GG vs. AG vs. AA, p = 0.048; GG vs. AG + AA, p = 0.026; GG vs. AG, p = 0.014), rs16840252 (CC vs. CT vs. TT, p = 0.007; CC vs. CT, p = 0.011), rs5742909 (CC vs. CT vs. TT, p = 0.040), and rs11571319 of CTLA4 (GG vs. AG vs. AA, p < 0.001; GG vs. AG + AA, p = 0.048; AA vs. GG + AG, p = 0.001; GG vs. AA, p = 0.008; GG vs. AG, p ≤ 0.001); and rs10204525 (TT vs. CT + CC, p = 0.024; TT vs. CT, p = 0.021), rs2227982 (AA vs. GG, p = 0.047), rs36084323 (TT vs. CT vs. CC, p = 0.022; TT vs. CT + CC, p = 0.013; CC vs. TT + CT, p = 0.048; TT vs. CC, p = 0.008), and rs5839828 of PDCD1 (DEL vs. DEL/G vs. GG, p = 0.014; DEL vs. DEL/G + GG, p = 0.014; GG vs. DEL + DEL/G, p = 0.025; DEL vs. GG, p = 0.007).

Discussion

Consequently, these SNPs may play an important role in immune regulation, and further research into the role of these SNPs of immune regulatory genes in the pathogenesis of RA is required.

Exploring the Genetic Predisposition to Epigenetic Changes in Alzheimer's Disease

Alzheimer's disease (AD) is a prevalent type of dementia in elderly populations with a significant genetic component. The accumulating evidence suggests that AD involves a reconfiguration of the epigenetic landscape, including DNA methylation, post-translational modification of histone proteins, and chromatin remodeling. Along with environmental factors, individual specific genetic features play a considerable role in the formation of epigenetic architecture. In this study, we attempt to identify the non-coding regulatory SNPs (rSNPs) able to affect the epigenetic mechanisms in AD. To this end, the multi-omics approach is used. The GEO (Gene Expression Omnibus) available data (GSE153875) for AD patients and controls are integrated to reveal the rSNPs that display allele-specific features in both ChIP-seq profiles of four histone modifications and RNA-seq. Furthermore, we analyze the presence of rSNPs in the promoters of genes reported to be differentially expressed between AD and the normal brain (AD-related genes) and involved in epigenetic regulation according to the EpiFactors database. We also searched for the rSNPs in the promoters of the genes coding for transcription regulators of the identified AD-related genes. These regulators were selected based on the corresponding ChIP-seq peaks (ENCODE) in the promoter regions of these genes. Finally, we formed a panel of rSNPs localized to the promoters of genes that contribute to the epigenetic landscape in AD and, thus, to the genetic predisposition for this disease.

Genetic Variation in Transcription Factor Binding Sites

The interaction between transcription factors (TFs) and DNA is the core process that determines the state of a cell's transcriptome [...].

An atrial fibrillation-associated regulatory region modulates cardiac Tbx5 levels and arrhythmia susceptibility

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation (AF). We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)) and one downstream (RE(down)) of Tbx5. In both lines, we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. Increased atrial Tbx5 expression was associated with induction of genes involved in development, ion transport and conduction, with increased susceptibility to atrial arrhythmias, and increased action potential duration of atrial cardiomyocytes. We identified an AF-associated variant in the human RE(int) that increases its transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in Tbx5^RE(int)KO cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in Tbx5^RE(int)KO mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.

SNPs in 3'UTR miRNA Target Sequences Associated with Individual Drug Susceptibility

The complementary interaction of microRNAs (miRNAs) with their binding sites in the 3'untranslated regions (3'UTRs) of target gene mRNAs represses translation, playing a leading role in gene expression control. MiRNA recognition elements (MREs) in the 3'UTRs of genes often contain single nucleotide polymorphisms (SNPs), which can change the binding affinity for target miRNAs leading to dysregulated gene expression. Accumulated data suggest that these SNPs can be associated with various human pathologies (cancer, diabetes, neuropsychiatric disorders, and cardiovascular diseases) by disturbing the interaction of miRNAs with their MREs located in mRNA 3'UTRs. Numerous data show the role of SNPs in 3'UTR MREs in individual drug susceptibility and drug resistance mechanisms. In this review, we brief the data on such SNPs focusing on the most rigorously proven cases. Some SNPs belong to conventional genes from the drug-metabolizing system (in particular, the genes coding for cytochromes P450 (CYP 450), phase II enzymes (SULT1A1 and UGT1A), and ABCB3 transporter and their expression regulators (PXR and GATA4)). Other examples of SNPs are related to the genes involved in DNA repair, RNA editing, and specific drug metabolisms. We discuss the gene-by-gene studies and genome-wide approaches utilized or potentially utilizable to detect the MRE SNPs associated with individual response to drugs.

PLK2 Single Nucleotide Variant in Gastric Cancer Patients Affects miR-23b-5p Binding

PURPOSE

Chromosomal instability is a hallmark of gastric cancer (GC). It can be driven by single nucleotide variants (SNVs) in cell cycle genes. We investigated the associations between SNVs in candidate genes, PLK2, PLK3, and ATM, and GC risk and clinicopathological features.

MATERIALS AND METHODS

The genotyping study included 542 patients with GC and healthy controls. Generalized linear models were used for the risk and clinicopathological association analyses. Survival analysis was performed using the Kaplan-Meier method. The binding of candidate miRs was analyzed using a luciferase reporter assay.

RESULTS

The PLK2 Crs15009-Crs963615 haplotype was under-represented in the GC group compared to that in the control group (Pcorr=0.050). Male patients with the PLK2 rs963615 CT genotype had a lower risk of GC, whereas female patients had a higher risk (P=0.023; P=0.026). The PLK2 rs963615 CT genotype was associated with the absence of vascular invasion (P=0.012). The PLK3 rs12404160 AA genotype was associated with a higher risk of GC in the male population (P=0.015). The ATM Trs228589-Ars189037-Grs4585 haplotype was associated with a higher risk of GC (P<0.001). The ATM rs228589, rs189037, and rs4585 genotypes TA+AA, AG+GG, and TG+GG were associated with the absence of perineural invasion (P=0.034). In vitro analysis showed that the cancer-associated miR-23b-5p mimic specifically bound to the PLK2 rs15009 G allele (P=0.0097). Moreover, low miR-23b expression predicted longer 10-year survival (P=0.0066) in patients with GC.

CONCLUSIONS

PLK2, PLK3, and ATM SNVs could potentially be helpful for the prediction of GC risk and clinicopathological features. PLK2 rs15009 affects the binding of miR-23b-5p. MiR-23b-5p expression status could serve as a prognostic marker for survival in patients with GC.

Regulation of cardiac ion channels by transcription factors: Looking for new opportunities of druggable targets for the treatment of arrhythmias

Cardiac electrical activity is governed by different ion channels that generate action potentials. Acquired or inherited abnormalities in the expression and/or function of ion channels usually result in electrophysiological changes that can cause cardiac arrhythmias. Transcription factors (TFs) control gene transcription by binding to specific DNA sequences adjacent to target genes. Linkage analysis, candidate-gene screening within families, and genome-wide association studies have linked rare and common genetic variants in the genes encoding TFs with genetically-determined cardiac arrhythmias. Besides its critical role in cardiac development, recent data demonstrated that they control cardiac electrical activity through the direct regulation of the expression and function of cardiac ion channels in adult hearts. This narrative review summarizes some studies showing functional data on regulation of the main human atrial and ventricular Na⁺, Ca²⁺, and K⁺ channels by cardiac TFs such as Pitx2c, Tbx20, Tbx5, Zfhx3, among others. The results have improved our understanding of the mechanisms regulating cardiac electrical activity and may open new avenues for therapeutic interventions in cardiac acquired or inherited arrhythmias through the identification of TFs as potential drug targets. Even though TFs have for a long time been considered as 'undruggable' targets, advances in structural biology have led to the identification of unique pockets in TFs amenable to be targeted with small-molecule drugs or peptides that are emerging as novel therapeutic drugs.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Methylation and expression quantitative trait locus rs6296 in the HTR1B gene is associated with susceptibility to opioid use disorder

Serotonin (5-HT) is implicated in the reward processes underlying substance use disorder. Epigenetic and transcriptional mechanisms contribute to the development of addictive states. To examine the potential mechanisms of 5-HT receptor genes in opioid use disorder, we first determined the associations between several single-nucleotide polymorphism (SNPs) in three representative 5-HT receptor genes (HTR1B, HTR2A, and HTR3B) and susceptibility to heroin use disorder in 1731 participants. Gene-gene interactions among these genes were analyzed. After identifying the susceptibility genes and SNPs for heroin use disorder, DNA methylation in the promoter region of these susceptibility genes was compared between 111 healthy controls and 120 patients with heroin use disorder. In addition, associations between the susceptibility SNPs and methylation of the CpG sites and gene promoters with differential methylation between groups were examined. Finally, the function of the susceptibility SNPs in the expression of the corresponding genes was screened. Our results demonstrated that rs6296 in the HTR1B gene was correlated with susceptibility to heroin use disorder. Gene-gene interactions between the HTR1B and HTR2A genes were identified. The CpG sites HTR1B_07 and HTR1B_26 and the promoter region of the HTR1B gene were hypermethylated in patients with heroin use disorder compared with healthy controls. Notably, rs6296 correlated in an allele-specific manner with methylation in the HTR1B gene promoter in the blood and gene expression of the HTR1B gene in the frontal cortex and hypothalamus. SNP rs6296 was associated with opioid use disorder by involving mechanisms of DNA methylation and expression of the HTR1B gene.

Positional weight matrices have sufficient prediction power for analysis of noncoding variants

The position weight matrix, also called the position-specific scoring matrix, is the commonly accepted model to quantify the specificity of transcription factor binding to DNA. Position weight matrices are used in thousands of projects and software tools in regulatory genomics, including computational prediction of the regulatory impact of single-nucleotide variants. Yet, recently Yan et al. reported that "the position weight matrices of most transcription factors lack sufficient predictive power" if applied to the analysis of regulatory variants studied with a newly developed experimental method, SNP-SELEX. Here, we re-analyze the rich experimental dataset obtained by Yan et al. and show that appropriately selected position weight matrices in fact can adequately quantify transcription factor binding to alternative alleles.

agReg-SNPdb-Plants: A Database of Regulatory SNPs for Agricultural Plant Species

Simple Summary

In breeding research, the investigation of regulatory SNPs (rSNPs) is becoming increasingly important due to their potential causal role for specific functional traits. Especially for crop species, there is still a lack of systematic analyses to detect rSNPs and their predicted effects on the binding of transcription factors. In this study, we present agReg-SNPdb-Plants, a database storing genome-wide collections of regulatory SNPs for agricultural plant species which can be queried via a web interface.

Abstract

Single nucleotide polymorphisms (SNPs) that are located in the promoter regions of genes and affect the binding of transcription factors (TFs) are called regulatory SNPs (rSNPs). Their identification can be highly valuable for the interpretation of genome-wide association studies (GWAS), since rSNPs can reveal the biologically causative variant and decipher the regulatory mechanisms behind a phenotype. In our previous work, we presented agReg-SNPdb, a database of regulatory SNPs for agriculturally important animal species. To complement this previous work, in this study we present the extension agReg-SNPdb-Plants storing rSNPs and their predicted effects on TF-binding for 13 agriculturally important plant species and subspecies (Brassica napus, Helianthus annuus, Hordeum vulgare, Oryza glaberrima, Oryza glumipatula, Oryza sativa Indica, Oryza sativa Japonica, Solanum lycopersicum, Sorghum bicolor, Triticum aestivum, Triticum turgidum, Vitis vinifera, and Zea mays). agReg-SNPdb-Plants can be queried via a web interface that allows users to search for SNP IDs, chromosomal regions, or genes. For a comprehensive interpretation of GWAS results or larger SNP-sets, it is possible to download the whole list of SNPs and their impact on transcription factor binding sites (TFBSs) from the website chromosome-wise.

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.

agReg-SNPdb: A Database of Regulatory SNPs for Agricultural Animal Species

Transcription factors (TFs) govern transcriptional gene regulation by specifically binding to short DNA motifs, known as transcription factor binding sites (TFBSs), in regulatory regions, such as promoters. Today, it is well known that single nucleotide polymorphisms (SNPs) in TFBSs can dramatically affect the level of gene expression, since they can cause a change in the binding affinity of TFs. Such SNPs, referred to as regulatory SNPs (rSNPs), have gained attention in the life sciences due to their causality for specific traits or diseases. In this study, we present agReg-SNPdb, a database comprising rSNP data of seven agricultural and domestic animal species: cattle, pig, chicken, sheep, horse, goat, and dog. To identify the rSNPs, we constructed a bioinformatics pipeline and identified a total of 10,623,512 rSNPs, which are located within TFBSs and affect the binding affinity of putative TFs. Altogether, we implemented the first systematic analysis of SNPs in promoter regions and their impact on the binding affinity of TFs for livestock and made it usable via a web interface.