A Panel of rSNPs Demonstrating Allelic Asymmetry in Both ChIP-seq and RNA-seq Data and the Search for Their Phenotypic Outcomes through Analysis of DEGs

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

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 [...].

Identification of Candidate Genes in Breast Cancer Induced by Estrogen Plus Progestogens Using Bioinformatic Analysis

Menopausal hormone therapy (MHT) was widely used to treat menopause-related symptoms in menopausal women. However, MHT therapies were controversial with the increased risk of breast cancer because of different estrogen and progestogen combinations, and the molecular basis behind this phenomenon is currently not understood. To address this issue, we identified differentially expressed genes (DEGs) between the estrogen plus progestogens treatment (EPT) and estrogen treatment (ET) using the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) data. As a result, a total of 96 upregulated DEGs were first identified. Seven DEGs related to the cell cycle (CCNE2, CDCA5, RAD51, TCF19, KNTC1, MCM10, and NEIL3) were validated by RT-qPCR. Specifically, these seven DEGs were increased in EPT compared to ET (p < 0.05) and had higher expression levels in breast cancer than adjacent normal tissues (p < 0.05). Next, we found that estrogen receptor (ER)-positive breast cancer patients with a higher CNNE2 expression have a shorter overall survival time (p < 0.05), while this effect was not observed in the other six DEGs (p > 0.05). Interestingly, the molecular docking results showed that CCNE2 might bind to 17β-estradiol (−6.791 kcal/mol), progesterone (−6.847 kcal/mol), and medroxyprogesterone acetate (−6.314 kcal/mol) with a relatively strong binding affinity, respectively. Importantly, CNNE2 protein level could be upregulated with EPT and attenuated by estrogen receptor antagonist, acolbifene and had interactions with cancer driver genes (AKT1 and KRAS) and high mutation frequency gene (TP53 and PTEN) in breast cancer patients. In conclusion, the current study showed that CCNE2, CDCA5, RAD51, TCF19, KNTC1, MCM10, and NEIL3 might contribute to EPT-related tumorigenesis in breast cancer, with CCNE2 might be a sensitive risk indicator of breast cancer risk in women using MHT.

ANANASTRA: annotation and enrichment analysis of allele-specific transcription factor binding at SNPs

We present ANANASTRA, https://ananastra.autosome.org, a web server for the identification and annotation of regulatory single-nucleotide polymorphisms (SNPs) with allele-specific binding events. ANANASTRA accepts a list of dbSNP IDs or a VCF file and reports allele-specific binding (ASB) sites of particular transcription factors or in specific cell types, highlighting those with ASBs significantly enriched at SNPs in the query list. ANANASTRA is built on top of a systematic analysis of allelic imbalance in ChIP-Seq experiments and performs the ASB enrichment test against background sets of SNPs found in the same source experiments as ASB sites but not displaying significant allelic imbalance. We illustrate ANANASTRA usage with selected case studies and expect that ANANASTRA will help to conduct the follow-up of GWAS in terms of establishing functional hypotheses and designing experimental verification.