rs2072580T>A Polymorphism in the Overlapping Promoter Regions of the SART3 and ISCU Genes Associated with the Risk of Breast Cancer.

We analyzed association of potentially regulatory polymorphisms (rs590352, rs11542583, rs3829202, rs207258, and rs4796672) with breast cancer. A significant association was found between this disease and rs2072580T>A (p=0.001) located in the overlapping promoter regions of the SART3 and ISCU genes. In women with AA and AT genotypes, the risk of breast cancer is higher by 6.7 times (p=0.001) and 12 times (p=0.001), respectively, in comparison with TT genotype. Under a codominant model of inheritance (AT vs AA+TT), the risk of breast cancer was increased by 4.2 times (р=0.001) for the AT genotype. Under a recessive model of inheritance (TT vs AA+TT), the risk of disease was 10-fold higher (р=0.001) for the TT genotype. It has been demonstrated that the T>A substitution affects the binding properties of transcription factors CREB1 and REST.

Mechanisms of Brain Glucocorticoid Resistance in Stress-Induced Psychopathologies.

Exposure to stress activates the hypothalamic-pituitary-adrenal axis and leads to increased levels of glucocorticoid (GC) hormones. Prolonged elevation of GC levels causes neuronal dysfunction, decreases the density of synapses, and impairs neuronal plasticity. Decreased sensitivity to glucocorticoids (glucocorticoid resistance) that develops as a result of chronic stress is one of the characteristic features of stress-induced psychopathologies. In this article, we reviewed the published data on proposed molecular mechanisms that contribute to the development of glucocorticoid resistance in brain, including changes in the expression of the glucocorticoid receptor (GR) gene, biosynthesis of GR isoforms, and GR posttranslational modifications. We also present data on alterations in the expression of the FKBP5 gene encoding the main component of cell ultra-short negative feedback loop of GC signaling regulation. Recent discoveries on stress- and GR-induced changes in epigenetic modification patterns as well as normalizing action of antidepressants are discussed. GR and FKBP5 gene polymorphisms associated with stress-induced psychopathologies are described, and their role in glucocorticoid resistance is discussed.

The reasons of Trithorax-like expression disturbance in Trl 3609 allele of Drosophila melanogaster.

The regulatory region of the Trl gene was analyzed using the mutation Trl 3609 , resulting from the insertion of the P-element into the promoter region of the gene as well as mutations obtained on its basis. It is shown that two last transcription start sites, which are most often used in vitro in S2 cells, are almost not used in vivo. Experimental data indicate that transcription terminators in transposons play an important role in the decrease in the transcription level of the recipient gene.

[Regulatory transcription codes in eukaryotic genomes].

Key aspects of gene transcription regulation in multicellular organisms, including the characteristics of their promoters, transcription-factor binding sites, and composition elements are reviewed. The functional role of transcription regulatory proteins (basal factors and regulatory transcription factors), and the mechanisms responsible for regulation of their activity are also discussed. Furthermore, we describe the importance of DNA-encoded nucleosome organization and chromatin modifications in the course of transcription regulation, as well as some mechanisms that regulate the activity of transcription factors associated with genetic networks. The current outlook on regulatory gene expression codes in eukaryotes is presented.

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.

The functional insight into the genetics of cardiovascular disease: results from the post-GWAS study.

Cardiovascular diseases (CVDs), the leading cause of death worldwide, generally refer to a range of pathological conditions with the involvement of the heart and the blood vessels. A sizable fraction of the susceptibility loci is known, but the underlying mechanisms have been established only for a small proportion. Therefore, there is an increasing need to explore the functional relevance of trait-associated variants and, moreover, to search for novel risk genetic variation. We have reported the bioinformatic approach allowing effective identif ication of functional non-coding variants by integrated analysis of genome-wide data. Here, the analysis of 1361 previously identif ied regulatory SNPs (rSNPs) was performed to provide new insights into cardiovascular risk. We found 773,471 coding co-segregating markers for input rSNPs using the 1000 Genomes Project. The intersection of GWAS-derived SNPs with a relevance to cardiovascular traits with these markers was analyzed within a window of 10 Kbp. The effects on the transcription factor (TF) binding sites were explored by DeFine models. Functional pathway enrichment and protein-protein interaction (PPI) network analyses were performed on the targets and the extended genes by STRING and DAVID. Eighteen rSNPs were functionally linked to cardiovascular risk. A signif icant impact on binding sites of thirteen TFs including those involved in blood cells formation, hematopoiesis, macrophage function, inf lammation, and vasoconstriction was found in K562 cells. 21 rSNP gene targets and 5 partners predicted by PPI were enriched for spliceosome and endocytosis KEGG pathways, endosome sorting complex and mRNA splicing REACTOME pathways. Related Gene Ontology terms included mRNA splicing and processing, endosome transport and protein catabolic processes. Together, the f indings provide further insight into the biological basis of CVDs and highlight the importance of the precise regulation of splicing and alternative splicing.

OAT and 3'MeDAB azo compounds similarly cause liver tumors in GR mice, but differently modify activities of FoxA transcription factors.

Transcription factors of the FoxA family (forkhead box A) regulate cell metabolism and differentiation and maintain specificity of liver cell proteome and phenotype of mature hepatocytes. The relationship between hepatocarcinogenicity of azo compounds o-aminoazotoluene (OAT) and 3'-methyl-4-dimethylaminobenzene (3'MeDAB) for GR mice and one of the early events, modulation of the DNA-binding activity of FoxA transcription factor, was studied. Single injection of 3'MeDAB to 12-day-old mice caused liver tumors in 100% males and females similarly as OAT, a well-known mouse hepatocarcinogene. The DNA-binding activity of FoxA in the liver decreased 2.5-3 times by OAT, this resulting in a 40% reduction of glucocorticoid induction of tyrosine aminotransferase (liver-specific gene). In contrast to these, 3'MeDAB did not modify FoxA protein activities or the degree of glucocorticoid induction of tyrosine aminotransferase.

Application of alternative de novo motif recognition models for analysis of structural heterogeneity of transcription factor binding sites: a case study of FOXA2 binding sites.

The most popular model for the search of ChIP-seq data for transcription factor binding sites (TFBS) is the positional weight matrix (PWM). However, this model does not take into account dependencies between nucleotide occurrences in different site positions. Currently, two recently proposed models, BaMM and InMoDe, can do as much. However, application of these models was usually limited only to comparing their recognition accuracies with that of PWMs, while none of the analyses of the co-prediction and relative positioning of hits of different models in peaks has yet been performed. To close this gap, we propose the pipeline called MultiDeNA. This pipeline includes stages of model training, assessing their recognition accuracy, scanning ChIP-seq peaks and their classification based on scan results. We applied our pipeline to 22 ChIP-seq datasets of TF FOXA2 and considered PWM, dinucleotide PWM (diPWM), BaMM and InMoDe models. The combination of these four models allowed a significant increase in the fraction of recognized peaks compared to that for the sole PWM model: the increase was 26.3 %. The BaMM model provided the main contribution to the recognition of sites. Although the major fraction of predicted peaks contained TFBS of different models with coincided positions, the medians of the fraction of peaks containing the predictions of sole models were 1.08, 0.49, 4.15 and 1.73 % for PWM, diPWM, BaMM and InMoDe, respectively. Thus, FOXA2 BSs were not fully described by only a sole model, which indicates theirs heterogeneity. We assume that the BaMM model is the most successful in describing the structure of the FOXA2 BS in ChIP-seq datasets under study.

Social defeat stress in adult mice causes alterations in gene expression, alternative splicing, and the epigenetic landscape of H3K4me3 in the prefrontal cortex: An impact of early-life stress.

Chronic stress is the leading risk factor of a broad range of severe psychopathologies. Nonetheless, the molecular mechanisms triggering these pathological processes are not well understood. In our study, we investigated the effects of 15-day social defeat stress (SDS) on the genome-wide landscape of trimethylation at the 4th lysine residue of histone H3 (H3K4me3) and on the transcriptome in the prefrontal cortex of mice that were reared normally (group SDS) or subjected to maternal separation early in life (group MS+SDS). The mice with the history of stress early in life showed increased susceptibility to SDS in adulthood and demonstrated long-lasting genome-wide alterations in gene expression and splicing as well as in the H3K4me3 epigenetic landscape in the prefrontal cortex. Thus, the high-throughput techniques applied here allowed us to simultaneously detect, for the first time, genome-wide epigenetic and transcriptional changes in the murine prefrontal cortex that are associated with both chronic SDS and increased susceptibility to this stressor.

[Effect of hepatocarcinogenicity of estragole on the glucocorticoid-mediated induction of liver-specific enzymes and the activity of the transcription factors FOXA and HNF4 in the liver of mouse and rat].

The carcinogenic effects of estragole in mice of the earlier unexplored strain ICR has been studied. It has been shown that there is a distinct correlation between the extent of inhibition of glucocorticoid-mediated induction of tyrosine aminotransferase and trypthophan oxygenase after acute administration of estragole and the frequency of liver tumors after estragole exposure. Estragole inhibits the induction of these enzymes only in female mice, but not in male mice and rats. DNA-binding activities of liver-enriched transcription factors were investigated on carcinogen-susceptible and -resistant animals. Estragole decreases the HNF4 (hepatic nuclear factor 4) and FOXA DNA-binding activities only in susceptible female mice, but not in nonsusceptible male mice and rats and does not influence the C/EBP and HNF1 activities. Pentachlorophenol, which prevents the hepatocarcinogenic effect of estragole, abolishes its inhibitory effect on tyrosine aminotransferase and trypthophan oxygenase glucocorticoid induction and restores the FOXA and HNF4 DNA-binding activities. The parallelism between the hepatocarcinogenic effects of estragole and the inhibition of FOXA and HNF4 DNA-binding activities serves as an additional argument for the involvement of these factors in the mechanisms of tumor suppression in the liver.

Data of correlation analysis between the density of H3K4me3 in promoters of genes and gene expression: Data from RNA-seq and ChIP-seq analyses of the murine prefrontal cortex.

H3K4me3 is typically found in the promoter region of genes and is a mark associated with an open chromatin state and active gene transcription. Nonetheless, the role of H3K4me3 in the regulation of transcription is still debated. To improve the understanding of the connection between H3K4me3 density in promoters and gene expression, we assessed the correlation between these two parameters. We utilized genome-wide high-throughput RNA sequencing (RNA-seq) data and H3K4me3-based chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq), carried out on the same samples of the prefrontal cortex from 10 male C57Bl6 mice with different stress experience [Social defeat stress in adult mice causes alterations in gene expression, alternative splicing, and the epigenetic landscape of H3K4me3 in the prefrontal cortex: an impact of early-life stress, 1]. In addition, we assessed the correlation between H3K4me3 density and gene expression in datasets of cell-specific genes. Altogether, the results are useful for the elucidation of H3K4me3 involvement in the regulation of transcription in the murine prefrontal cortex.

Correlation between hepatocarcinogenic effect of estragole and its influence on glucocorticoid induction of liver-specific enzymes and activities of FOXA and HNF4 transcription factors in mouse and rat liver.

It is known that the carcinogenic effect of estragole, a component of essential oils of many spicy plants, is characterized by species, tissue, and sex specificity. It causes mainly liver tumors in female mice but is not carcinogenic for male mice and for rats. In this work, the estragole hepatocarcinogenicity was shown for female mice of previously not studied ICR line. The strict correlation between estragole hepatocarcinogenicity and its ability to decrease the level of glucocorticoid induction of liver-specific enzymes tyrosine aminotransferase (TAT) and tryptophan oxygenase (TO) was found. Inhibition of TAT and TO inducibility by estragole takes place only in female mice but not in male mice and in rats. Studying the estragole effect on DNA-binding activity of transcription factors, present mainly in liver and regulating expression of genes encoding liver-specific proteins, has shown that estragole decreases FOXA and HNF4 activities but not activities of C/EBP and HNF1, and this happens only in female mice, for which this substance is hepatocarcinogen, but not in male mice and in rats. Pentachlorophenol, preventing hepatocarcinogenic effect of estragole, abolishes inhibitory influence of the latter on the TAT and TO glucocorticoid induction and restores DNA-binding activity of FOXA and HNF4. Thus, a correlation was revealed between the estragole hepatocarcinogenic effect and decrease in DNA-binding activity of transcription factors FOXA and HNF4, which might be indicative of the role of these factors in tumor suppression mechanisms in liver.

rSNP_Guide, a database system for analysis of transcription factor binding to target sequences: application to SNPs and site-directed mutations.

rSNP_Guide is a novel curated database system for analysis of transcription factor (TF) binding to target sequences in regulatory gene regions altered by mutations. It accumulates experimental data on naturally occurring site variants in regulatory gene regions and site-directed mutations. This database system also contains the web tools for SNP analysis, i.e., active applet applying weight matrices to predict the regulatory site candidates altered by a mutation. The current version of the rSNP_Guide is supplemented by six sub-databases: (i) rSNP_DB, on DNA-protein interaction caused by mutation; (ii) SYSTEM, on experimental systems; (iii) rSNP_BIB, on citations to original publications; (iv) SAMPLES, on experimentally identified sequences of known regulatory sites; (v) MATRIX, on weight matrices of known TF sites; (vi) rSNP_Report, on characteristic examples of successful rSNP_Tools implementation. These databases are useful for the analysis of natural SNPs and site-directed mutations. The databases are available through the Web, http://wwwmgs.bionet.nsc.ru/mgs/systems/rsnp/.

Transcription Regulatory Regions Database (TRRD): its status in 2002.

Transcription Regulatory Regions Database (TRRD) is an informational resource containing an integrated description of the gene transcription regulation. An entry of the database corresponds to a gene and contains the data on localization and functions of the transcription regulatory regions as well as gene expression patterns. TRRD contains only experimental data that are inputted into the database through annotating scientific publication. TRRD release 6.0 comprises the information on 1167 genes, 5537 transcription factor binding sites, 1714 regulatory regions, 14 locus control regions and 5335 expression patterns obtained through annotating 3898 scientific papers. This information is arranged in seven databases: TRRDGENES (general gene description), TRRDLCR (locus control regions); TRRDUNITS (regulatory regions: promoters, enhancers, silencers, etc.), TRRDSITES (transcription factor binding sites), TRRDFACTORS (transcription factors), TRRDEXP (expression patterns) and TRRDBIB (experimental publications). Sequence Retrieval System (SRS) is used as a basic tool for navigating and searching TRRD and integrating it with external informational and software resources. The visualization tool, TRRD Viewer, provides the information representation in a form of maps of gene regulatory regions. The option allowing nucleotide sequences to be searched for according to their homology using BLAST is also included. TRRD is available at http://www.bionet.nsc.ru/trrd/.

[Recognition of the potential SF-1 binding sites by SiteGA method, their experimental verification and search for new SF-1 target genes].

The SF-1 (Steroidogenic Factor-1) is a transcription factor known as a key regulator of the steroidogenic gene expression. SF-1 is required for the development and functioning at all levels of the hypothalamic-pituitary-gonadal and adrenal axis. Also it plays an essential role in sex determination. SF-1 is a member of the nuclear receptor superfamily and it activates gene expression by binding to DNA in a monomeric form. Here, we report the results of potential SF-1 binding sites identification by using the SiteGA recognition method. The SiteGA method was implemented using a genetic algorithm (GA) involving a iterative discriminant analyses of local dinucleotide context characteristics. These characteristics were compiled not only over the core binding sites region but over its flanks as well. Developed SiteGA method is characterized by considerably better recognition accuracy when compared to that for the weight matrix method. The experimental tests demonstrated that 83% of the sites recognized by the SiteGA method in the regulatory regions of steroidogenic genes, indeed, interact with the SF-1 factor. We also estimated the density of predicted sites in regulatory region of genes, the members of different functional groups and developed the criterion to search for new SF-1 target genes in genome sequences.

[Method SiteGA for the recognition of transcription factor binding sites].

A new approach, SiteGA, for the prediction of functional transcription factor binding sites has been developed. The approach is based on the detection of locally positioned dinucleotides by the genetic algorithm and discriminant analysis. The approach has been applied to recognize transcription factor binding sites involved in the regulation of immune responses and cell growth (AP-1, IRF1, ISGF3, NFkappaB, STAT1), obesity and lipid metabolism (HNF4, PPAR, SREBP), and the expression of steroidogenesis genes (SF-1). SiteGA is far superior in accuracy to the traditionally used method of position weight matrices. The approach was implemented in the web tool, SiteGA http://wwwmgs2. bionet.nsc.ru/mgs/programs/sitega.

Hidden heterogeneity of transcription factor binding sites: A case study of SF-1.

Steroidogenic factor 1 (SF-1) belongs to a small group of the transcription factors that bind DNA only as a monomer. Three different approaches-Sitecon, SiteGA, and oPWM-constructed using the same training sample of experimentally confirmed SF-1 binding sites have been used to recognize these sites. The appropriate prediction thresholds for recognition models have been selected. Namely, the thresholds concordant by false positive or negative rates for various methods were used to optimize the discrimination of steroidogenic gene promoters from the datasets of non-specific promoters. After experimental verification, the models were used to analyze the ChIP-seq data for SF-1. It has been shown that the sets of sites recognized by different models overlap only partially and that an integration of these models allows for identification of SF-1 sites in up to 80% of the ChIP-seq loci. The structures of the sites detected using the three recognition models in the ChIP-seq peaks falling within the [-5000, +5000] region relative to the transcription start sites (TSS) extracted from the FANTOM5 project have been analyzed. The MATLIGN classified the frequency matrices for the sites predicted by oPWM, Sitecon, and SiteGA into two groups. The first group is described by oPWM/Sitecon and the second, by SiteGA. Gene ontology (GO) analysis has been used to clarify the differences between the sets of genes carrying different variants of SF-1 binding sites. Although this analysis in general revealed a considerable overlap in GO terms for the genes carrying the binding sites predicted by oPWM, Sitecon, or SiteGA, only the last method elicited notable trend to terms related to negative regulation and apoptosis. The results suggest that the SF-1 binding sites are different in both their structure and the functional annotation of the set of target genes correspond to the predictions by oPWM+Sitecon and SiteGA. Further application of Homer software for de novo identification of enriched motifs in ChIP-Seq data for SF-1ChIP-seq dataset gave the data similar to oPWM+Sitecon.

Nucleotide sequence of a fragment of the rat tryptophan oxygenase gene showing high affinity to glucocorticoid receptor in vitro.

A nucleotide sequence (3.2 kb) of a DNA region located approximately between introns 4 and 7 of the rat tryptophan oxygenase gene was determined. Using filter binding studies and monoclonal antibodies against the glucocorticoid receptor a high affinity binding of this region to the glucocorticoid receptor was demonstrated.

New elements of glucocorticoid-receptor binding sites of hormone-regulated genes.

The structure of the DNA regions recognized by glucocorticoid-receptor complexes (GIRC) was analyzed using frequency matrices and a modified perceptron method. Some complementary conservative elements which may modulate the efficiency of GIRC binding were found at both sides of the previously established conserved nucleotide sequence (core) (Beato, M. et al. (1987) J. Steroid Biochem. 27, 9-14). A criterion based on the concurrent use of several perceptron matrices to search for the potential GIRC binding site sequences has been worked out. By applying this criterion 73 sites were identified in 28 sequences of glucocorticoid regulated genes and 7 sites were identified in 26 sequences independent from glucocorticoid regulation.

Glucocorticoid-receptor binding sites at the 5' flanking region of a rat cytochrome CYP2B2 gene predicted with a novel computer method.

Using an original computer method to search for potential DNA binding sites for glucocorticoid-receptor complexes (GRC) (Seledtsov, I.A., Solovjev, V.V. and Merkulova, T.I. (1991) Biochim. Biophys. Acta 1089, 367-376), the presence of two such sites in the 5' flanking region of a rat cytochrome CYP2B2 gene has been predicted. This prediction has been confirmed by gel retardation experiments.