Sequence features that drive human promoter function and tissue specificity

Dynamic interplay of Sp1, YY1, and DUX4 in regulating FRG1 transcription with intricate balance

Maintaining precise levels of FRG1 is vital. It's over-expression is tied to muscular dystrophy, while reduced levels are linked to tumorigenesis. Despite extensive efforts to characterize FRG1 expression and downstream molecular signaling, a comprehensive understanding of its regulation has remained elusive. This study focused on unravelling the cis -regulatory elements within the FRG1 gene and their interplay. Employing a dual luciferase reporter assay on fragments of the FRG1 promoter upstream of the transcription start site, we observed variations in FRG1 transcription induction. Our in-silico analysis unveiled binding sequences for Sp1 and DUX4 within FRG1 promoter region showing an enhanced luciferase signal. Conversely, we identified a YY1 binding sequence in the FRG1 promoter fragment showing decreased luciferase signal. Confirming these binding sites through site-directed mutagenesis, chromatin immunoprecipitation, and EMSA provided concrete evidence of Sp1, YY1, and DUX4's interaction within the FRG1 promoter. Additionally, interaction between Sp1, YY1, and DUX4 was elucidated using sequential chromatin immunoprecipitation (ChIP re-ChIP) and co-immunoprecipitation assays. Furthermore, alterations in the expression levels of Sp1, YY1, and DUX4 resulted in parallel changes in FRG1 gene expression. Notably, YY1 exhibited the ability to suppress SP1 or DUX4-mediated FRG1 transcription activation, while Sp1 and DUX4 together could counteract YY1-mediated transcription suppression. Our cell proliferation and colony formation assay underscored the tumorigenic properties of these three transcription factors through the modulation of FRG1 expression levels. The in vitro results were verified in vivo using mouse xenograft model. Leveraging RNA sequencing data from various tissues in the GTEx portal, we established a correlation between FRG1, Sp1, and YY1. In essence, this study revealed the vital cis-regulatory components residing in the FRG1 promoter. The combined influence of Sp1, YY1, and DUX4 plays a central role in controlling FRG1 expression.

Construction of transcript regulation mechanism prediction models based on binding motif environment of transcription factor AoXlnR in Aspergillus oryzae

DNA-binding transcription factors (TFs) play a central role in transcriptional regulation mechanisms, mainly through their specific binding to target sites on the genome and regulation of the expression of downstream genes. Therefore, a comprehensive analysis of the function of these TFs will lead to the understanding of various biological mechanisms. However, the functions of TFs in vivo are diverse and complicated, and the identified binding sites on the genome are not necessarily involved in the regulation of downstream gene expression. In this study, we investigated whether DNA structural information around the binding site of TFs can be used to predict the involvement of the binding site in the regulation of the expression of genes located downstream of the binding site. Specifically, we calculated the structural parameters based on the DNA shape around the DNA binding motif located upstream of the gene whose expression is directly regulated by one TF AoXlnR from Aspergillus oryzae, and showed that the presence or absence of expression regulation can be predicted from the sequence information with high accuracy ([Formula: see text]-1.0) by machine learning incorporating these parameters.

Identification of teleost tnnc1a enhancers for specific pan-cardiac transcription

Troponin C regulates muscle contraction by forming the troponin complex with troponin I and troponin T. Different muscle types express different troponin C genes. The mechanisms of such differential transcription are not fully understood. The Zebrafish tnnc1a gene is restrictively expressed in cardiac muscles. We here identify the enhancers and promoters of the zebrafish and medaka tnnc1a genes, including intronic enhancers in zebrafish and medaka and an upstream enhancer in the medaka. The intronic and upstream enhancers are likely functionally redundant. The GFP transgenic reporter driven by these enhancers is expressed more strongly in the ventricle than in the atrium, recapitulating the expression pattern of the endogenous zebrafish tnnc1a gene. Our study identifies a new set of enhancers for cardiac-specific transgenic expression in zebrafish. These enhancers can serve as tools for future identification of transcription factor networks that drive cardiac-specific gene transcription.

Tissue-specific and tissue-agnostic effects of genome sequence variation modulating blood pressure

Genome-wide association studies (GWASs) have identified numerous variants associated with polygenic traits and diseases. However, with few exceptions, a mechanistic understanding of which variants affect which genes in which tissues to modulate trait variation is lacking. Here, we present genomic analyses to explain trait heritability of blood pressure (BP) through the genetics of transcriptional regulation using GWASs, multiomics data from different tissues, and machine learning approaches. Approximately 500,000 predicted regulatory variants across four tissues explain 33.4% of variant heritability: 2.5%, 5.3%, 7.7%, and 11.8% for kidney-, adrenal-, heart-, and artery-specific variants, respectively. Variation in the enhancers involved shows greater tissue specificity than in the genes they regulate, suggesting that gene regulatory networks perturbed by enhancer variants in a tissue relevant to a phenotype are the major source of interindividual variation in BP. Thus, our study provides an approach to scan human tissue and cell types for their physiological contribution to any trait.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Compatibility rules of human enhancer and promoter sequences

Gene regulation in the human genome is controlled by distal enhancers that activate specific nearby promoters1. A proposed model for this specificity is that promoters have sequence-encoded preferences for certain enhancers, for example, mediated by interacting sets of transcription factors or cofactors2. This 'biochemical compatibility' model has been supported by observations at individual human promoters and by genome-wide measurements in Drosophila3-9. However, the degree to which human enhancers and promoters are intrinsically compatible has not yet been systematically measured, and how their activities combine to control RNA expression remains unclear. Here we design a high-throughput reporter assay called enhancer × promoter self-transcribing active regulatory region sequencing (ExP STARR-seq) and applied it to examine the combinatorial compatibilities of 1,000 enhancer and 1,000 promoter sequences in human K562 cells. We identify simple rules for enhancer-promoter compatibility, whereby most enhancers activate all promoters by similar amounts, and intrinsic enhancer and promoter activities multiplicatively combine to determine RNA output (R2 = 0.82). In addition, two classes of enhancers and promoters show subtle preferential effects. Promoters of housekeeping genes contain built-in activating motifs for factors such as GABPA and YY1, which decrease the responsiveness of promoters to distal enhancers. Promoters of variably expressed genes lack these motifs and show stronger responsiveness to enhancers. Together, this systematic assessment of enhancer-promoter compatibility suggests a multiplicative model tuned by enhancer and promoter class to control gene transcription in the human genome.

Massively parallel reporter perturbation assays uncover temporal regulatory architecture during neural differentiation

Gene regulatory elements play a key role in orchestrating gene expression during cellular differentiation, but what determines their function over time remains largely unknown. Here, we perform perturbation-based massively parallel reporter assays at seven early time points of neural differentiation to systematically characterize how regulatory elements and motifs within them guide cellular differentiation. By perturbing over 2,000 putative DNA binding motifs in active regulatory regions, we delineate four categories of functional elements, and observe that activity direction is mostly determined by the sequence itself, while the magnitude of effect depends on the cellular environment. We also find that fine-tuning transcription rates is often achieved by a combined activity of adjacent activating and repressing elements. Our work provides a blueprint for the sequence components needed to induce different transcriptional patterns in general and specifically during neural differentiation.

How gene regulatory elements regulate gene expression during cellular differentiation remains largely unknown. Here the authors use perturbation-based massively parallel reporter assays at early time points of neural differentiation to systematically characterize how regulatory elements and motifs within them guide different transcriptional patterns.

Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.

Nuclear FGFR1 Regulates Gene Transcription and Promotes Antiestrogen Resistance in ER+ Breast Cancer

PURPOSE

FGFR1 overexpression has been associated with endocrine resistance in ER⁺ breast cancer. We found FGFR1 localized in the nucleus of breast cancer cells in primary tumors resistant to estrogen suppression. We investigated a role of nuclear FGFR1 on gene transcription and antiestrogen resistance.

EXPERIMENTAL DESIGN

Tumors from patients treated with letrozole were subjected to Ki67 and FGFR1 IHC. MCF7 cells were transduced with FGFR1(SP-)(NLS) to promote nuclear FGFR1 overexpression. FGFR1 genomic activity in ER⁺/FGFR1-amplified breast cancer cells ± FOXA1 siRNA or ± the FGFR tyrosine kinase inhibitor (TKI) erdafitinib was examined by chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq). The nuclear and chromatin-bound FGFR1 interactome was investigated by mass spectrometry (MS).

RESULTS

High nuclear FGFR1 expression in ER⁺ primary tumors positively correlated with post-letrozole Ki67 values. Nuclear FGFR1 overexpression influenced gene transcription and promoted resistance to estrogen suppression and to fulvestrant in vivo. A gene expression signature induced by nuclear FGFR1 correlated with shorter survival in the METABRIC cohort of patients treated with antiestrogens. ChIP-Seq revealed FGFR1 occupancy at transcription start sites, overlapping with active transcription histone marks. MS analysis of the nuclear FGFR1 interactome identified phosphorylated RNA-Polymerase II and FOXA1, with FOXA1 RNAi impairing FGFR1 recruitment to chromatin. Treatment with erdafitinib did not impair nuclear FGFR1 translocation and genomic activity.

CONCLUSIONS

These data suggest nuclear FGFR1 contributes to endocrine resistance by modulating gene transcription in ER⁺ breast cancer. Nuclear FGFR1 activity was unaffected by FGFR TKIs, thus supporting the development of treatment strategies to inhibit nuclear FGFR1 in ER+/FGFR1 overexpressing breast cancer.

The Long Road to Understanding RNAPII Transcription Initiation and Related Syndromes

In eukaryotes, transcription of protein-coding genes requires the assembly at core promoters of a large preinitiation machinery containing RNA polymerase II (RNAPII) and general transcription factors (GTFs). Transcription is potentiated by regulatory elements called enhancers, which are recognized by specific DNA-binding transcription factors that recruit cofactors and convey, following chromatin remodeling, the activating cues to the preinitiation complex. This review summarizes nearly five decades of work on transcription initiation by describing the sequential recruitment of diverse molecular players including the GTFs, the Mediator complex, and DNA repair factors that support RNAPII to enable RNA synthesis. The elucidation of the transcription initiation mechanism has greatly benefited from the study of altered transcription components associated with human diseases that could be considered transcription syndromes.

Genetic and epigenetic features of promoters with ubiquitous chromatin accessibility support ubiquitous transcription of cell-essential genes

Gene expression is controlled by regulatory elements within accessible chromatin. Although most regulatory elements are cell type-specific, a subset is accessible in nearly all the 517 human and 94 mouse cell and tissue types assayed by the ENCODE consortium. We systematically analyzed 9000 human and 8000 mouse ubiquitously-accessible candidate cis-regulatory elements (cCREs) with promoter-like signatures (PLSs) from ENCODE, which we denote ubi-PLSs. These are more CpG-rich than non-ubi-PLSs and correspond to genes with ubiquitously high transcription, including a majority of cell-essential genes. ubi-PLSs are enriched with motifs of ubiquitously-expressed transcription factors and preferentially bound by transcriptional cofactors regulating ubiquitously-expressed genes. They are highly conserved between human and mouse at the synteny level but exhibit frequent turnover of motif sites; accordingly, ubi-PLSs show increased variation at their centers compared with flanking regions among the ∼186 thousand human genomes sequenced by the TOPMed project. Finally, ubi-PLSs are enriched in genes implicated in Mendelian diseases, especially diseases broadly impacting most cell types, such as deficiencies in mitochondrial functions. Thus, a set of roughly 9000 mammalian promoters are actively maintained in an accessible state across cell types by a distinct set of transcription factors and cofactors to ensure the transcriptional programs of cell-essential genes.

Predicting mRNA Abundance Directly from Genomic Sequence Using Deep Convolutional Neural Networks

Algorithms that accurately predict gene structure from primary sequence alone were transformative for annotating the human genome. Can we also predict the expression levels of genes based solely on genome sequence? Here, we sought to apply deep convolutional neural networks toward that goal. Surprisingly, a model that includes only promoter sequences and features associated with mRNA stability explains 59% and 71% of variation in steady-state mRNA levels in human and mouse, respectively. This model, termed Xpresso, more than doubles the accuracy of alternative sequence-based models and isolates rules as predictive as models relying on chromatic immunoprecipitation sequencing (ChIP-seq) data. Xpresso recapitulates genome-wide patterns of transcriptional activity, and its residuals can be used to quantify the influence of enhancers, heterochromatic domains, and microRNAs. Model interpretation reveals that promoter-proximal CpG dinucleotides strongly predict transcriptional activity. Looking forward, we propose cell-type-specific gene-expression predictions based solely on primary sequences as a grand challenge for the field.

Transcription initiation of distant core promoters in a large-sized genome of an insect

Background

Core promoters have a substantial influence on various steps of transcription, including initiation, elongation, termination, polyadenylation, and finally, translation. The characterization of core promoters is crucial for exploring the regulatory code of transcription initiation. However, the current understanding of insect core promoters is focused on those of Diptera (especially Drosophila) species with small genome sizes.

Results

Here, we present an analysis of the transcription start sites (TSSs) in the migratory locust, Locusta migratoria, which has a genome size of 6.5 Gb. The genomic differences, including lower precision of transcription initiation and fewer constraints on the distance from transcription factor binding sites or regulatory elements to TSSs, were revealed in locusts compared with Drosophila insects. Furthermore, we found a distinct bimodal log distribution of the distances from the start codons to the core promoters of locust genes. We found stricter constraints on the exon length of mRNA leaders and widespread expression activity of the distant core promoters in locusts compared with fruit flies. We further compared core promoters in seven arthropod species across a broad range of genome sizes to reinforce our results on the emergence of distant core promoters in large-sized genomes.

Conclusions

In summary, our results provide novel insights into the effects of genome size expansion on distant transcription initiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01004-5.

Epigenetic potential affects immune gene expression in house sparrows

Epigenetic mechanisms may play a central role in mediating phenotypic plasticity, especially during range expansions, when populations face a suite of novel environmental conditions. Individuals may differ in their epigenetic potential (EP; their capacity for epigenetic modifications of gene expression), which may affect their ability to colonize new areas. One form of EP, the number of CpG sites, is higher in introduced house sparrows (Passer domesticus) than in native birds in the promoter region of a microbial surveillance gene, Toll-like Receptor 4 (TLR4), which may allow invading birds to fine-tune their immune responses to unfamiliar parasites. Here, we compared TLR4 gene expression from whole blood, liver and spleen in house sparrows with different EP, first challenging some birds with lipopolysaccharide (LPS), to increase gene expression by simulating a natural infection. We expected that high EP would predict high inducibility and reversibility of TLR4 expression in the blood of birds treated with LPS, but we did not make directional predictions regarding organs, as we could not repeatedly sample these tissues. We found that EP was predictive of TLR4 expression in all tissues. Birds with high EP expressed more TLR4 in the blood than individuals with low EP, regardless of treatment with LPS. Only females with high EP exhibited reversibility in gene expression. Further, the effect of EP varied between sexes and among tissues. Together, these data support EP as one regulator of TLR4 expression.

CD28 Gene Polymorphisms in the Promoter Region Are Associated with Transfusion Reactions: A Functional Study

Transfusion reactions are mainly induced by the interaction of an antigen and antibody. However, transfusion reactions still occur with the implementing of crossmatching and usage of pre-storage leukoreduced blood products. The roles of CD28 and CTLA4 gene polymorphisms in transfusion reaction have been shown, and subjects with certain single nucleotide polymorphisms (SNPs) of the CD28 or CTLA4 gene had a significantly higher risk of transfusion reactions. In total, 40 patients with transfusion reactions after receiving pre-storage leukoreduced blood products were enrolled in this study. We focused on the SNPs located in the CD28 promoter region (rs1879877, rs3181096, rs3181097, and rs3181098) to find out the significant SNP. A luciferase reporter assay was used to investigate the expression level of protein affected by promoter SNP variation. We found that the polymorphism of rs3181097 was associated with transfusion reactions (p = 0.003 in additive model and p = 0.015 in dominant model). Consequently, we investigated the biological function in the CD28 promoter polymorphisms (rs1879877 G > T, rs3181096 C > T, rs3181097 G > A, and rs3181098 G > A) by using dual-spectral luciferase reporter assay. The results showed that the ex-pression level of CD28 was decreased under the effect of rs3181097 with A-allele. This suggested that rs3181097 may regulate immune response through decreasing CD28 protein expression and then lead to development of transfusion reactions.

Long first exons and epigenetic marks distinguish conserved pachytene piRNA clusters from other mammalian genes

In the male germ cells of placental mammals, 26-30-nt-long PIWI-interacting RNAs (piRNAs) emerge when spermatocytes enter the pachytene phase of meiosis. In mice, pachytene piRNAs derive from ~100 discrete autosomal loci that produce canonical RNA polymerase II transcripts. These piRNA clusters bear 5' caps and 3' poly(A) tails, and often contain introns that are removed before nuclear export and processing into piRNAs. What marks pachytene piRNA clusters to produce piRNAs, and what confines their expression to the germline? We report that an unusually long first exon (≥ 10 kb) or a long, unspliced transcript correlates with germline-specific transcription and piRNA production. Our integrative analysis of transcriptome, piRNA, and epigenome datasets across multiple species reveals that a long first exon is an evolutionarily conserved feature of pachytene piRNA clusters. Furthermore, a highly methylated promoter, often containing a low or intermediate level of CG dinucleotides, correlates with germline expression and somatic silencing of pachytene piRNA clusters. Pachytene piRNA precursor transcripts bind THOC1 and THOC2, THO complex subunits known to promote transcriptional elongation and mRNA nuclear export. Together, these features may explain why the major sources of pachytene piRNA clusters specifically generate these unique small RNAs in the male germline of placental mammals.

DNMTs and Impact of CpG Content, Transcription Factors, Consensus Motifs, lncRNAs, and Histone Marks on DNA Methylation

DNA methyltransferases (DNMTs) play an essential role in DNA methylation and transcriptional regulation in the genome. DNMTs, along with other poorly studied elements, modulate the dynamic DNA methylation patterns of embryonic and adult cells. We summarize the current knowledge on the molecular mechanism of DNMTs’ functional targeting to maintain genome-wide DNA methylation patterns. We focus on DNMTs’ intrinsic characteristics, transcriptional regulation, and post-transcriptional modifications. Furthermore, we focus special attention on the DNMTs’ specificity for target sites, including key cis-regulatory factors such as CpG content, common motifs, transcription factors (TF) binding sites, lncRNAs, and histone marks to regulate DNA methylation. We also review how complexes of DNMTs/TFs or DNMTs/lncRNAs are involved in DNA methylation in specific genome regions. Understanding these processes is essential because the spatiotemporal regulation of DNA methylation modulates gene expression in health and disease.

Analysis of endogenous and exogenous tumor upregulated promoter expression in canine tumors

Gene therapy is a promising treatment option for cancer. However, its utility may be limited due to expression in off-target cells. Cancer-specific promoters such as telomerase reverse transcriptase (TERT), survivin, and chemokine receptor 4 (CXCR4) have enhanced activity in a variety of human and murine cancers, however, little has been published regarding these promoters in dogs. Given the utility of canine cancer models, the activity of these promoters along with adenoviral E2F enhanced E1a promoter (EEE) was evaluated in a variety of canine tumors, both from the endogenous gene and from exogenously administered constructs. Endogenous expression levels were measured for cTERT, cSurvivin, and cCXCR4 and were low for all three, with some non-malignant and some tumor cell lines and tissues expressing the gene. Expression levels from exogenously supplied promoters were measured by both the number of cells expressing the construct and the intensity of expression in individual cells. Exogenously supplied promoters were active in more cells in all tumor lines than in normal cells, with the EEE promoter being most active, followed by cTERT. The intensity of expression varied more with cell type than with specific promoters. Ultimately, no single promoter was identified that would result in reliable expression, regardless of the tumor type. Thus, these findings imply that identification of a pan-cancer promoter may be difficult. In addition, this data raises the concern that endogenous expression analysis may not accurately predict exogenous promoter activity.

Pervasive and CpG-dependent promoter-like characteristics of transcribed enhancers

The temporal and spatial expression of genes is controlled by promoters and enhancers. Findings obtained over the last decade that not only promoters but also enhancers are characterized by bidirectional, divergent transcription have challenged the traditional notion that promoters and enhancers represent distinct classes of regulatory elements. Over half of human promoters are associated with CpG islands (CGIs), relatively CpG-rich stretches of generally several hundred nucleotides that are often associated with housekeeping genes. Only about 6% of transcribed enhancers defined by CAGE-tag analysis are associated with CGIs. Here, we present an analysis of enhancer and promoter characteristics and relate them to the presence or absence of CGIs. We show that transcribed enhancers share a number of CGI-dependent characteristics with promoters, including statistically significant local overrepresentation of core promoter elements. CGI-associated enhancers are longer, display higher directionality of transcription, greater expression, a lesser degree of tissue specificity, and a higher frequency of transcription-factor binding events than non-CGI-associated enhancers. Genes putatively regulated by CGI-associated enhancers are enriched for transcription regulator activity. Our findings show that CGI-associated transcribed enhancers display a series of characteristics related to sequence, expression and function that distinguish them from enhancers not associated with CGIs.

A report on DNA sequence determinants in gene expression

The biased usage of nucleotides in coding sequence and its correlation with gene expression has been observed in several studies. A complex set of interactions between genes and other components of the expression system determine the amount of proteins produced from coding sequences. It is known that the elongation rate of polypeptide chain is affected by both codon usage bias and specific amino acid compositional constraints. Therefore, it is of interest to review local DNA-sequence elements and other positional as well as combinatorial constraints that play significant role in gene expression.

The Gossypium longicalyx Genome as a Resource for Cotton Breeding and Evolution

Cotton is an important crop that has made significant gains in production over the last century. Emerging pests such as the reniform nematode have threatened cotton production. The rare African diploid species Gossypium longicalyx is a wild species that has been used as an important source of reniform nematode immunity. While mapping and breeding efforts have made some strides in transferring this immunity to the cultivated polyploid species, the complexities of interploidal transfer combined with substantial linkage drag have inhibited progress in this area. Moreover, this species shares its most recent common ancestor with the cultivated A-genome diploid cottons, thereby providing insight into the evolution of long, spinnable fiber. Here we report a newly generated de novo genome assembly of G. longicalyx. This high-quality genome leveraged a combination of PacBio long-read technology, Hi-C chromatin conformation capture, and BioNano optical mapping to achieve a chromosome level assembly. The utility of the G. longicalyx genome for understanding reniform immunity and fiber evolution is discussed.

A comparison of curated gene sets versus transcriptomics-derived gene signatures for detecting pathway activation in immune cells

BACKGROUND

Despite the significant contribution of transcriptomics to the fields of biological and biomedical research, interpreting long lists of significantly differentially expressed genes remains a challenging step in the analysis process. Gene set enrichment analysis is a standard approach for summarizing differentially expressed genes into pathways or other gene groupings. Here, we explore an alternative approach to utilizing gene sets from curated databases. We examine the method of deriving custom gene sets which may be relevant to a given experiment using reference data sets from previous transcriptomics studies. We call these data-derived gene sets, "gene signatures" for the biological process tested in the previous study. We focus on the feasibility of this approach in analyzing immune-related processes, which are complicated in their nature but play an important role in the medical research.

RESULTS

We evaluate several statistical approaches to detecting the activity of a gene signature in a target data set. We compare the performance of the data-derived gene signature approach with comparable GO term gene sets across all of the statistical tests. A total of 61 differential expression comparisons generated from 26 transcriptome experiments were included in the analysis. These experiments covered eight immunological processes in eight types of leukocytes. The data-derived signatures were used to detect the presence of immunological processes in the test data with modest accuracy (AUC = 0.67). The performance for GO and literature based gene sets was worse (AUC = 0.59). Both approaches were plagued by poor specificity.

CONCLUSIONS

When investigators seek to test specific hypotheses, the data-derived signature approach can perform as well, if not better than standard gene-set based approaches for immunological signatures. Furthermore, the data-derived signatures can be generated in the cases that well-defined gene sets are lacking from pathway databases and also offer the opportunity for defining signatures in a cell-type specific manner. However, neither the data-derived signatures nor standard gene-sets can be demonstrated to reliably provide negative predictions for negative cases. We conclude that the data-derived signature approach is a useful and sometimes necessary tool, but analysts should be weary of false positives.

A Method for the Structure-Based, Genome-Wide Analysis of Bacterial Intergenic Sequences Identifies Shared Compositional and Functional Features

In this paper, we propose a computational strategy for performing genome-wide analyses of intergenic sequences in bacterial genomes. Following similar directions of a previous paper, where a method for genome-wide analysis of eucaryotic Intergenic sequences was proposed, here we developed a tool for implementing similar concepts in bacteria genomes. This allows us to (i) classify intergenic sequences into clusters, characterized by specific global structural features and (ii) draw possible relations with their functional features.

Bioinformatic and functional analysis of promoter region of human SLC25A13 gene

The human SLC25A13 gene encodes the liver type aspartate/glutamate carrier isoform 2 (AGC2, commonly named as citrin), which plays a key role in the main NADH-shuttle of human hepatocyte. Biallelic SLC25A13 mutations result in Citrin deficiency (CD). In order to identify the important regulatory region of SLC25A13 gene and elucidate the way how potential promoter mutations affect the citrin expression, we performed promoter deletion analysis and established the reporter constructs of luciferase gene-carrying SLC25A13 promoter containing several mutations located in putative transcription factor-binding sites. The luciferase activities of all promoter constructs were measured using a Dual-Luciferase Reporter Assay System. Bioinformatic analysis showed that the promoter of SLC25A13 gene lacks TATA box and obviously typical initiator element, but contains a CCAAT box and two GC box. Promoter deletion analysis confirmed the region from -221 to -1 upstream ATG was essential for SLC25A13 to maintain the promoter activity. We utilized dual-luciferase reporter system as function analytical model to tentatively assess the effect of artificially constructed promoter mutations on citrin expression, and our analysis revealed that mutated putative CCAAT box and GC box could significantly affect the citrin expression. Our study confirmed the important SLC25A13 promoter regions that influenced citrin expression in HL7702 cells, and constructed a function analytical model. This work may be useful to further identify the pathogenic mutations leading to CD in the promoter region.

Functional conservation of sequence determinants at rapidly evolving regulatory regions across mammals

Recent advances in epigenomics have made it possible to map genome-wide regulatory regions using empirical methods. Subsequent comparative epigenomic studies have revealed that regulatory regions diverge rapidly between genome of different species, and that the divergence is more pronounced in enhancers than in promoters. To understand genomic changes underlying these patterns, we investigated if we can identify specific sequence fragments that are over-enriched in regulatory regions, thus potentially contributing to regulatory functions of such regions. Here we report numerous sequence fragments that are statistically over-enriched in enhancers and promoters of different mammals (which we refer to as ‘sequence determinants’). Interestingly, the degree of statistical enrichment, which presumably is associated with the degree of regulatory impacts of the specific sequence determinant, was significantly higher for promoter sequence determinants than enhancer sequence determinants. We further used a machine learning method to construct prediction models using sequence determinants. Remarkably, prediction models constructed from one species could be used to predict regulatory regions of other species with high accuracy. This observation indicates that even though the precise locations of regulatory regions diverge rapidly during evolution, the functional potential of sequence determinants underlying regulatory sequences may be conserved between species.

Regions of the genome that do not encode genes but affect expression of other genes, such as enhancers and promoters, are referred to as regulatory regions. Because of their regulatory functions, it was thought that enhancers and promoters should be evolutionarily conserved. Regulatory regions can be now epigenomically identified because they are marked by specific modifications of histone tails at the chromatin level. Interestingly, when we compare epigenomically identified regulatory regions from different mammals, the specific positions of regulatory regions are often divergent between species. Enhancers in particular are highly divergent between species. In this study, we show that we can find sequence fragments that are statistically enriched in enhancers and promoters of different species, and that the degree of statistical enrichment can explain different levels of evolutionary sequence conservation between enhancers and promoters. We further constructed predictive models of enhancers and promoters using the enriched sequence fragments, and show that these models can not only accurately predict enhancers and promoters of the same species, but works comparably well when applied to other species. These results indicate that even though the specific positions of regulatory regions have diverged between species, the functions of sequence fragments that comprise those regions may be conserved.

CRX directs photoreceptor differentiation by accelerating chromatin remodeling at specific target sites

Background

Recent technological advances have delivered the genome-wide targets of many important transcription factors (TFs). However, increasing evidence suggests that not all target sites mediate regulatory function, raising the questions of how to determine which sites are active, what are the epigenetic consequences of TF binding at these sites, and how the specificity is coded. To address these questions, we focused on CRX, a disease-associated homeodomain TF required for photoreceptor gene expression and development. Since CRX binds more than 6000 sites across the genome in the retina, we profiled chromatin landscape changes at each binding site during normal development and in the absence of CRX and interpreted the results by thorough investigation of other epigenomic datasets and sequence features.

Results

CRX is required for chromatin remodeling at only a subset of its binding sites, which undergo retina or neuronal specific activation during photoreceptor differentiation. Genes near these “CRX Dependent” sites code for proteins important for photoreceptor physiology and function, and their transcription is significantly reduced in Crx deficient retinas. In addition, the nucleotide and motif content distinguish these CRX Dependent sites from other CRX-bound sites.

Conclusions

Together, our results suggest that CRX acts only at select, uniquely-coded binding sites to accelerate chromatin remodeling during photoreceptor differentiation. This study emphasizes the importance of connecting TF binding with its functional consequences and provides a framework for making such a connection using comparative analyses of available genomic datasets. Finally, this study prioritizes sets of non-coding DNA sites for future functional interrogation and identification of mutations associated with retinal disease.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0212-2) contains supplementary material, which is available to authorized users.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

Relationship between histone modifications and transcription factor binding is protein family specific

The very small fraction of putative binding sites (BSs) that are occupied by transcription factors (TFs) in vivo can be highly variable across different cell types. This observation has been partly attributed to changes in chromatin accessibility and histone modification (HM) patterns surrounding BSs. Previous studies focusing on BSs within DNA regulatory regions found correlations between HM patterns and TF binding specificities. However, a mechanistic understanding of TF-DNA binding specificity determinants is still not available. The ability to predict in vivo TF binding on a genome-wide scale requires the identification of features that determine TF binding based on evolutionary relationships of DNA binding proteins. To reveal protein family-dependent mechanisms of TF binding, we conducted comprehensive comparisons of HM patterns surrounding BSs and non-BSs with exactly matched core motifs for TFs in three cell lines: 33 TFs in GM12878, 37 TFs in K562, and 18 TFs in H1-hESC. These TFs displayed protein family-specific preferences for HM patterns surrounding BSs, with high agreement among cell lines. Moreover, compared to models based on DNA sequence and shape at flanking regions of BSs, HM-augmented quantitative machine-learning methods resulted in increased performance in a TF family-specific manner. Analysis of the relative importance of features in these models indicated that TFs, displaying larger HM pattern differences between BSs and non-BSs, bound DNA in an HM-specific manner on a protein family-specific basis. We propose that TF family-specific HM preferences reveal distinct mechanisms that assist in guiding TFs to their cognate BSs by altering chromatin structure and accessibility.

ZFX acts as a transcriptional activator in multiple types of human tumors by binding downstream from transcription start sites at the majority of CpG island promoters

High expression of the transcription factor ZFX is correlated with proliferation, tumorigenesis, and patient survival in multiple types of human cancers. However, the mechanism by which ZFX influences transcriptional regulation has not been determined. We performed ChIP-seq in four cancer cell lines (representing kidney, colon, prostate, and breast cancers) to identify ZFX binding sites throughout the human genome. We identified roughly 9000 ZFX binding sites and found that most of the sites are in CpG island promoters. Moreover, genes with promoters bound by ZFX are expressed at higher levels than genes with promoters not bound by ZFX. To determine if ZFX contributes to regulation of the promoters to which it is bound, we performed RNA-seq analysis after knockdown of ZFX by siRNA in prostate and breast cancer cells. Many genes with promoters bound by ZFX were down-regulated upon ZFX knockdown, supporting the hypothesis that ZFX acts as a transcriptional activator. Surprisingly, ZFX binds at +240 bp downstream from the TSS of the responsive promoters. Using Nucleosome Occupancy and Methylome Sequencing (NOMe-seq), we show that ZFX binds between the open chromatin region at the TSS and the first downstream nucleosome, suggesting that ZFX may play a critical role in promoter architecture. We have also shown that a closely related zinc finger protein ZNF711 has a similar binding pattern at CpG island promoters, but ZNF711 may play a subordinate role to ZFX. This functional characterization of ZFX provides important new insights into transcription, chromatin structure, and the regulation of the cancer transcriptome.

Cancer-specific promoters for expression-targeted gene therapy: ran, brms1 and mcm5

BACKGROUND

To expand the library of promoters that can be used for expression-targeted gene delivery to cancer cells, the specificity and strength of expression of three cancer-related gene promoters was evaluated: RAS-related nuclear protein ((P) ran), breast cancer metastasis suppressor 1 ((P) brms1) and minichromosome maintenance complex component 5 ((P) mcm5).

METHODS

The expression of reporter genes under the control of these promoters demonstrated selectivity in cancer cell lines of breast, prostate and ovarian origins versus a panel of normal cell types. The (P) ran was next used to regulate the expression of a bioactive exon (a constitutively active form of human caspase 3) to induce apoptosis in cancer cells. Further evaluation was performed in an orthotopic model of murine bladder cancer.

RESULTS

The average strengths of reporter expression had relative intensities of 99.8% ((P) ran), 87.7% ((P) brms1) and 55.8% ((P) mcm5) versus the strong (P) cmv-driven positive control. Comparisons of expression-targeted reporter gene expression for these three promoters versus the clinically interesting promoter for the human telomerase reverse transcriptase gene ((P) hTERT) yielded an improvement of two- to 15-fold. Following transfection, cell death was evident from morphologic observations and viability assays performed on the cancer cells lines, with little (if any) effects seen when the same genes were delivered to normal cells. Cell viability was reduced by up to 60% after one treatment, with cell death via apoptosis implied by caspase 3 detection. During the in vivo preclinical study, reduced tumor burden, lack of mineralization and decreased inflammation were demonstrated after only three treatments.

CONCLUSIONS

The ran, brms1, and mcm5 promoters have the specificity and strength needed for cancer-specific expression-targeted gene therapy. (p) ran in particular produced exciting results when coupled with a version of the caspase 3 exon to treat bladder cancer. Copyright © 2016 John Wiley & Sons, Ltd.

Local sequence features that influence AP-1 cis-regulatory activity

In the genome, most occurrences of transcription factor binding sites (TFBS) have no cis-regulatory activity, which suggests that flanking sequences contain information that distinguishes functional from nonfunctional TFBS. We interrogated the role of flanking sequences near Activator Protein 1 (AP-1) binding sites that reside in DNase I Hypersensitive Sites (DHS) and regions annotated as Enhancers. In these regions, we found that sequence features directly adjacent to the core motif distinguish high from low activity AP-1 sites. Some nearby features are motifs for other TFs that genetically interact with the AP-1 site. Other features are extensions of the AP-1 core motif, which cause the extended sites to match motifs of multiple AP-1 binding proteins. Computational models trained on these data distinguish between sequences with high and low activity AP-1 sites and also predict changes in cis-regulatory activity due to mutations in AP-1 core sites and their flanking sequences. Our results suggest that extended AP-1 binding sites, together with adjacent binding sites for additional TFs, encode part of the information that governs TFBS activity in the genome.

A new computational method to predict transcriptional activity of a DNA sequence from diverse datasets of massively parallel reporter assays

In recent years, the dramatic increase in the number of applications for massively parallel reporter assay (MPRA) technology has produced a large body of data for various purposes. However, a computational model that can be applied to decipher regulatory codes for diverse MPRAs does not exist yet. Here, we propose a new computational method to predict the transcriptional activity of MPRAs, as well as luciferase reporter assays, based on the TRANScription FACtor database. We employed regression trees and multivariate adaptive regression splines to obtain these predictions and considered a feature redundancy-dependent formula for conventional regression trees to enable adaptation to diverse data. The developed method was applicable to various MPRAs despite the use of different types of transfected cells, sequence lengths, construct numbers and sequence types. We demonstrate that this method can predict the transcriptional activity of promoters in HEK293 cells through predictive functions that were estimated by independent assays in eight tumor cell lines. The prediction was generally good (Pearson's r = 0.68) which suggested that common active transcription factor binding sites across different cell types make greater contributions to transcriptional activity and that known promoter activity could confer transcriptional activity of unknown promoters in some instances, regardless of cell type.

Baf60b-mediated ATM-p53 activation blocks cell identity conversion by sensing chromatin opening

Lineage conversion by expression of lineage-specific transcription factors is a process of epigenetic remodeling that has low efficiency. The mechanism by which a cell resists lineage conversion is largely unknown. Using hepatic-specific transcription factors Foxa3, Hnf1α and Gata4 (3TF) to induce hepatic conversion in mouse fibroblasts, we showed that 3TF induced strong activation of the ATM-p53 pathway, which led to proliferation arrest and cell death, and it further prevented hepatic conversion. Notably, ATM activation, independent of DNA damage, responded to chromatin opening during hepatic conversion. By characterizing the early molecular events during hepatic conversion, we found that Baf60b, a member of the SWI/SNF chromatin remodeling complex, links chromatin opening to ATM activation by facilitating ATM recruitment to the open chromatin regions of a panel of hepatic gene loci. These findings shed light on cellular responses to lineage conversion by revealing a function of the ATM-p53 pathway in sensing chromatin opening.

Systematic dissection of genomic features determining transcription factor binding and enhancer function

Enhancers regulate gene expression through the binding of sequence-specific transcription factors (TFs) to cognate motifs. Various features influence TF binding and enhancer function-including the chromatin state of the genomic locus, the affinities of the binding site, the activity of the bound TFs, and interactions among TFs. However, the precise nature and relative contributions of these features remain unclear. Here, we used massively parallel reporter assays (MPRAs) involving 32,115 natural and synthetic enhancers, together with high-throughput in vivo binding assays, to systematically dissect the contribution of each of these features to the binding and activity of genomic regulatory elements that contain motifs for PPARγ, a TF that serves as a key regulator of adipogenesis. We show that distinct sets of features govern PPARγ binding vs. enhancer activity. PPARγ binding is largely governed by the affinity of the specific motif site and higher-order features of the larger genomic locus, such as chromatin accessibility. In contrast, the enhancer activity of PPARγ binding sites depends on varying contributions from dozens of TFs in the immediate vicinity, including interactions between combinations of these TFs. Different pairs of motifs follow different interaction rules, including subadditive, additive, and superadditive interactions among specific classes of TFs, with both spatially constrained and flexible grammars. Our results provide a paradigm for the systematic characterization of the genomic features underlying regulatory elements, applicable to the design of synthetic regulatory elements or the interpretation of human genetic variation.

Expression of metastasis suppressor gene AES driven by a Yin Yang (YY) element in a CpG island promoter and transcription factor YY2

We recently found that the product of the AES gene functions as a metastasis suppressor of colorectal cancer (CRC) in both humans and mice. Expression of amino‐terminal enhancer of split (AES) protein is significantly decreased in liver metastatic lesions compared with primary colon tumors. To investigate its downregulation mechanism in metastases, we searched for transcriptional regulators of AES in human CRC and found that its expression is reduced mainly by transcriptional dysregulation and, in some cases, by additional haploidization of its coding gene. The AES promoter‐enhancer is in a typical CpG island, and contains a Yin‐Yang transcription factor recognition sequence (YY element). In human epithelial cells of normal colon and primary tumors, transcription factor YY2, a member of the YY family, binds directly to the YY element, and stimulates expression of AES. In a transplantation mouse model of liver metastases, however, expression of Yy2 (and therefore of Aes) is downregulated. In human CRC metastases to the liver, the levels of AES protein are correlated with those of YY2. In addition, we noticed copy‐number reduction for the AES coding gene in chromosome 19p13.3 in 12% (5/42) of human CRC cell lines. We excluded other mechanisms such as point or indel mutations in the coding or regulatory regions of the AES gene, CpG methylation in the AES promoter enhancer, expression of microRNAs, and chromatin histone modifications. These results indicate that Aes may belong to a novel family of metastasis suppressors with a CpG‐island promoter enhancer, and it is regulated transcriptionally.

Genome-wide mapping of autonomous promoter activity in human cells

Previous methods to systematically characterize sequence-intrinsic activity of promoters have been limited by relatively low throughput and the length of the sequences that could be tested. Here we present 'survey of regulatory elements' (SuRE), a method that assays more than 10⁸ DNA fragments, each 0.2-2 kb in size, for their ability to drive transcription autonomously. In SuRE, a plasmid library of random genomic fragments upstream of a 20-bp barcode is constructed, and decoded by paired-end sequencing. This library is used to transfect cells, and barcodes in transcribed RNA are quantified by high-throughput sequencing. When applied to the human genome, we achieve 55-fold genome coverage, allowing us to map autonomous promoter activity genome-wide in K562 cells. By computational modeling we delineate subregions within promoters that are relevant for their activity. We show that antisense promoter transcription is generally dependent on the sense core promoter sequences, and that most enhancers and several families of repetitive elements act as autonomous transcription initiation sites.

Genomic approaches for understanding the genetics of complex disease

There are thousands of known associations between genetic variants and complex human phenotypes, and the rate of novel discoveries is rapidly increasing. Translating those associations into knowledge of disease mechanisms remains a fundamental challenge because the associated variants are overwhelmingly in noncoding regions of the genome where we have few guiding principles to predict their function. Intersecting the compendium of identified genetic associations with maps of regulatory activity across the human genome has revealed that phenotype-associated variants are highly enriched in candidate regulatory elements. Allele-specific analyses of gene regulation can further prioritize variants that likely have a functional effect on disease mechanisms; and emerging high-throughput assays to quantify the activity of candidate regulatory elements are a promising next step in that direction. Together, these technologies have created the ability to systematically and empirically test hypotheses about the function of noncoding variants and haplotypes at the scale needed for comprehensive and systematic follow-up of genetic association studies. Major coordinated efforts to quantify regulatory mechanisms across genetically diverse populations in increasingly realistic cell models would be highly beneficial to realize that potential.

High-throughput functional comparison of promoter and enhancer activities

Promoters initiate RNA synthesis, and enhancers stimulate promoter activity. Whether promoter and enhancer activities are encoded distinctly in DNA sequences is unknown. We measured the enhancer and promoter activities of thousands of DNA fragments transduced into mouse neurons. We focused on genomic loci bound by the neuronal activity-regulated coactivator CREBBP, and we measured enhancer and promoter activities both before and after neuronal activation. We find that the same sequences typically encode both enhancer and promoter activities. However, gene promoters generate more promoter activity than distal enhancers, despite generating similar enhancer activity. Surprisingly, the greater promoter activity of gene promoters is not due to conventional core promoter elements or splicing signals. Instead, we find that particular transcription factor binding motifs are intrinsically biased toward the generation of promoter activity, whereas others are not. Although the specific biases we observe may be dependent on experimental or cellular context, our results suggest that gene promoters are distinguished from distal enhancers by specific complements of transcriptional activators.

Ovine HSP90AA1 gene promoter: functional study and epigenetic modifications

When environmental temperatures exceed a certain threshold, the upregulation of the ovine HSP90AA1 gene is produced to cope with cellular injuries caused by heat stress. It has been previously pointed out that several polymorphisms located at the promoter region of this gene seem to be the main responsible for the differences in the heat stress response observed among alternative genotypes in terms of gene expression rate. The present study, focused on the functional study of those candidate polymorphisms by electrophoretic mobility shift assay (EMSA) and in vitro luciferase expression assays, has revealed that the observed differences in the transcriptional activity of the HSP90AA1 gene as response to heat stress are caused by the presence of a cytosine insertion (rs397514115) and a C to G transversion (rs397514116) at the promoter region. Next, we discovered the presence of epigenetic marks at the promoter and along the gene body founding an allele-specific methylation of the rs397514116 mutation in DNA extracted from blood samples. This regulatory mechanism interacts synergistically to modulate gene expression depending on environmental circumstances. Taking into account the results obtained, it is suggested that the transcription of the HSP90AA1 ovine gene is regulated by a cooperative action of transcription factors (TFs) whose binding sites are polymorphic and where the influence of epigenetic events should be also taken into account.

SMAD4 gene promoter mutations in patients with thyroid tumors

As a key component of the transforming growth factor beta (TGFB) pathway, which regulates the expression of thyroid-specific genes, tumor suppressor SMAD4 is crucial for thyroid development and function. Aberrant expression of SMAD4 in thyroid tumor tissue was reported and mutations affecting the coding region have been detected, but a potential role of mutations in SMAD4 gene regulatory regions remains unexplored. The aim of this study was to analyze SMAD4 gene promoters in thyroid tumors. A total of 76 thyroidectomy specimens were studied, including 42 malignant and 34 benign tumors. The presence of mutations in four SMAD4 gene promoters was analyzed in thyroid tumor tissue and peripheral blood by PCR and DNA sequencing. The expression and intracellular localization of endogenous SMAD4 protein in selected tumor samples was studied by immunostaining and confocal microscopy. Of three novel variants detected, two were within promoter A (-204T/C and -5C/T) and one in promoter D (-180delA). Unlike somatic mutations previously detected in the nearby region, germline mutation -180delA in promoter D doesn't appear to affect SMAD4 expression in the thyroid tumor tissue. However, all newly detected SMAD4 promoter variants affect predicted binding sites of transcription factors involved in cell cycle regulation and should be further characterized functionally. Although not directly involved in carcinogenesis, detected variants may alter SMAD4 transcriptional regulation to some extent. Considering that dosage dependence is of great importance for the role of SMAD4 protein as a tumor suppressor, potential clinical significance of SMAD4 gene promoter mutations is worth further investigation.

Species differential regulation of COX2 can be described by an NFκB-dependent logic AND gate

Cyclooxygenase 2 (COX2), a key regulatory enzyme of the prostaglandin/eicosanoid pathway, is an important target for anti-inflammatory therapy. It is highly induced by pro-inflammatory cytokines in a Nuclear factor kappa B (NFκB)-dependent manner. However, the mechanisms determining the amplitude and dynamics of this important pro-inflammatory event are poorly understood. Furthermore, there is significant difference between human and mouse COX2 expression in response to the inflammatory stimulus tumor necrosis factor alpha (TNFα). Here, we report the presence of a molecular logic AND gate composed of two NFκB response elements (NREs) which controls the expression of human COX2 in a switch-like manner. Combining quantitative kinetic modeling and thermostatistical analysis followed by experimental validation in iterative cycles, we show that the human COX2 expression machinery regulated by NFκB displays features of a logic AND gate. We propose that this provides a digital, noise-filtering mechanism for a tighter control of expression in response to TNFα, such that a threshold level of NFκB activation is required before the promoter becomes active and initiates transcription. This NFκB-regulated AND gate is absent in the mouse COX2 promoter, most likely contributing to its differential graded response in promoter activity and protein expression to TNFα. Our data suggest that the NFκB-regulated AND gate acts as a novel mechanism for controlling the expression of human COX2 to TNFα, and its absence in the mouse COX2 provides the foundation for further studies on understanding species-specific differential gene regulation.

Cancer specificity of promoters of the genes controlling cell proliferation

Violation of proliferation control is a common feature of cancer cells. We put forward the hypothesis that promoters of genes involved in the control of cell proliferation should possess intrinsic cancer specific activity. We cloned promoter regions of CDC6, POLD1, CKS1B, MCM2, and PLK1 genes into pGL3 reporter vector and studied their ability to drive heterologous gene expression in transfected cancer cells of different origin and in normal human fibroblasts. Each promoter was cloned in short (335-800 bp) and long (up to 2.3 kb) variants to cover probable location of core and whole promoter regulatory elements. Cloned promoters were significantly more active in cancer cells than in normal fibroblasts that may indicate their cancer specificity. Both versions of CDC6 promoters were shown to be most active while the activities of others were close to that of BIRC5 gene (survivin) gene promoter. Long and short variants of each cloned promoter demonstrated very similar cancer specificity with the exception of PLK1-long promoter that was substantially more specific than its short variant and other promoters under study. The data indicate that most of the important cis-regulatory transcription elements responsible for intrinsic cancer specificity are located in short variants of the promoters under study. CDC6 short promoter may serve as a promising candidate for transcription targeted cancer gene therapy.

Relating gene expression evolution with CpG content changes

Background

Previous studies have shown that CpG dinucleotides are enriched in a subset of promoters and the CpG content of promoters is positively correlated with gene expression levels. But the relationship between divergence of CpG content and gene expression evolution has not been investigated. Here we calculate the normalized CpG (nCpG) content in DNA regions around transcription start site (TSS) and transcription terminal site (TTS) of genes in nine organisms, and relate them with expression levels measured by RNA-seq.

Results

The nCpG content of TSS shows a bimodal distribution in all organisms except platypus, whereas the nCpG content of TTS only has a single peak. When the nCpG contents are compared between different organisms, we observe a different evolution pattern between TSS and TTS: compared with TTS, TSS exhibits a faster divergence rate between closely related species but are more conserved between distant species. More importantly, we demonstrate the link between gene expression evolution and nCpG content changes: up-/down- regulation of genes in an organism is accompanied by the nCpG content increase/decrease in their TSS and TTS proximal regions.

Conclusions

Our results suggest that gene expression changes between different organisms are correlated with the alterations in normalized CpG contents of promoters. Our analyses provide evidences for the impact of nCpG content on gene expression evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-693) contains supplementary material, which is available to authorized users.

High-throughput functional testing of ENCODE segmentation predictions

The histone modification state of genomic regions is hypothesized to reflect the regulatory activity of the underlying genomic DNA. Based on this hypothesis, the ENCODE Project Consortium measured the status of multiple histone modifications across the genome in several cell types and used these data to segment the genome into regions with different predicted regulatory activities. We measured the cis-regulatory activity of more than 2000 of these predictions in the K562 leukemia cell line. We tested genomic segments predicted to be Enhancers, Weak Enhancers, or Repressed elements in K562 cells, along with other sequences predicted to be Enhancers specific to the H1 human embryonic stem cell line (H1-hESC). Both Enhancer and Weak Enhancer sequences in K562 cells were more active than negative controls, although surprisingly, Weak Enhancer segmentations drove expression higher than did Enhancer segmentations. Lower levels of the covalent histone modifications H3K36me3 and H3K27ac, thought to mark active enhancers and transcribed gene bodies, associate with higher expression and partly explain the higher activity of Weak Enhancers over Enhancer predictions. While DNase I hypersensitivity (HS) is a good predictor of active sequences in our assay, transcription factor (TF) binding models need to be included in order to accurately identify highly expressed sequences. Overall, our results show that a significant fraction (-26%) of the ENCODE enhancer predictions have regulatory activity, suggesting that histone modification states can reflect the cis-regulatory activity of sequences in the genome, but that specific sequence preferences, such as TF-binding sites, are the causal determinants of cis-regulatory activity.

Use of phage display to identify novel mineralocorticoid receptor-interacting proteins

The mineralocorticoid receptor (MR) plays a central role in salt and water homeostasis via the kidney; however, inappropriate activation of the MR in the heart can lead to heart failure. A selective MR modulator that antagonizes MR signaling in the heart but not the kidney would provide the cardiovascular protection of current MR antagonists but allow for normal electrolyte balance. The development of such a pharmaceutical requires an understanding of coregulators and their tissue-selective interactions with the MR, which is currently limited by the small repertoire of MR coregulators described in the literature. To identify potential novel MR coregulators, we used T7 phage display to screen tissue-selective cDNA libraries for MR-interacting proteins. Thirty MR binding peptides were identified, from which three were chosen for further characterization based on their nuclear localization and their interaction with other MR-interacting proteins or, in the case of x-ray repair cross-complementing protein 6, its known status as an androgen receptor coregulator. Eukaryotic elongation factor 1A1, structure-specific recognition protein 1, and x-ray repair cross-complementing protein 6 modulated MR-mediated transcription in a ligand-, cell- and/or promoter-specific manner and colocalized with the MR upon agonist treatment when imaged using immunofluorescence microscopy. These results highlight the utility of phage display for rapid and sensitive screening of MR binding proteins and suggest that eukaryotic elongation factor 1A1, structure-specific recognition protein 1, and x-ray repair cross-complementing protein 6 may be potential MR coactivators whose activity is dependent on the ligand, cellular context, and target gene promoter.

MOF-associated complexes ensure stem cell identity and Xist repression

Histone acetyl transferases (HATs) play distinct roles in many cellular processes and are frequently misregulated in cancers. Here, we study the regulatory potential of MYST1-(MOF)-containing MSL and NSL complexes in mouse embryonic stem cells (ESCs) and neuronal progenitors. We find that both complexes influence transcription by targeting promoters and TSS-distal enhancers. In contrast to flies, the MSL complex is not exclusively enriched on the X chromosome, yet it is crucial for mammalian X chromosome regulation as it specifically regulates Tsix, the major repressor of Xist lncRNA. MSL depletion leads to decreased Tsix expression, reduced REX1 recruitment, and consequently, enhanced accumulation of Xist and variable numbers of inactivated X chromosomes during early differentiation. The NSL complex provides additional, Tsix-independent repression of Xist by maintaining pluripotency. MSL and NSL complexes therefore act synergistically by using distinct pathways to ensure a fail-safe mechanism for the repression of X inactivation in ESCs.

DOI:http://dx.doi.org/10.7554/eLife.02024.001

Gene expression is controlled by a complicated network of mechanisms involving a wide range of enzymes and protein complexes. Many of these mechanisms are identical in males and females, but some are not. Female mammals, for example, carry two X chromosomes, whereas males have one X and one Y chromosome. Since the two X chromosomes in females contain essentially the same set of genes, one of them undergoes silencing to prevent the overproduction of certain proteins. This process, which is called X-inactivation, occurs during different stages of development and it must be tightly controlled.

An enzyme called MOF was originally found in flies in two distinct complexes—the male-specific lethal (MSL) complex, which forms only in males, and the non-specific lethal (NSL) complex, which is ubiquitous in both males and females. These complexes are evolutionary conserved and are also found in mammals. While mammalian MOF is reasonably well understood, the MSL and NSL complexes are not, so Chelmicki, Dündar et al. have used various sequencing techniques, in combination with biochemical experiments, to investigate their roles in embryonic stem cells and neuronal progenitor cells in mice.

These experiments show that MSL and NSL complexes engage in the regulation of thousands of genes. Although the two complexes often show different gene preferences, they often regulate the same cellular processes. The MSL/NSL-dependent regulation of X chromosome inactivation is a prime example of this phenomenon.

The MSL complex reduces the production of an RNA molecule called Xist, which is responsible for the inactivation of one of the two X chromosomes in females. The NSL complex, meanwhile, ensures the production of multiple proteins that are crucial for the development of embryonic stem cells, and are also involved in the repression of X inactivation.

This analysis sheds light on how different complexes can cooperate and complement each other in order to reach the same goal in the cell. The knowledge gained from this study will pave the way towards better understanding of complex processes such as embryonic development, organogenesis and the pathogenesis of disorders like cancer.

DOI:http://dx.doi.org/10.7554/eLife.02024.002

FIREWACh: high-throughput functional detection of transcriptional regulatory modules in mammalian cells

Promoters and enhancers establish precise gene transcription patterns. The development of functional approaches for their identification in mammalian cells has been complicated by the size of these genomes. Here we report a high-throughput functional assay for directly identifying active promoter and enhancer elements called FIREWACh (Functional Identification of Regulatory Elements Within Accessible Chromatin), which we used to simultaneously assess over 80,000 DNA fragments derived from nucleosome-free regions within the chromatin of embryonic stem cells (ESCs) and identify 6,364 active regulatory elements. Many of these represent newly discovered ESC-specific enhancers, showing enriched binding-site motifs for ESC-specific transcription factors including SOX2, POU5F1 (OCT4) and KLF4. The application of FIREWACh to additional cultured cell types will facilitate functional annotation of the genome and expand our view of transcriptional network dynamics.

HMGA proteins as modulators of chromatin structure during transcriptional activation

High mobility group (HMG) proteins are the most abundant non-histone chromatin associated proteins. HMG proteins bind to DNA and nucleosome and alter the structure of chromatin locally and globally. Accessibility to DNA within chromatin is a central factor that affects DNA-dependent nuclear processes, such as transcription, replication, recombination, and repair. HMG proteins associate with different multi-protein complexes to regulate these processes by mediating accessibility to DNA. HMG proteins can be subdivided into three families: HMGA, HMGB, and HMGN. In this review, we will focus on recent advances in understanding the function of HMGA family members, specifically their role in gene transcription regulation during development and cancer.

Genome-wide analysis of promoters: clustering by alignment and analysis of regular patterns

In this paper we perform a genome-wide analysis of H. sapiens promoters. To this aim, we developed and combined two mathematical methods that allow us to (i) classify promoters into groups characterized by specific global structural features, and (ii) recover, in full generality, any regular sequence in the different classes of promoters. One of the main findings of this analysis is that H. sapiens promoters can be classified into three main groups. Two of them are distinguished by the prevalence of weak or strong nucleotides and are characterized by short compositionally biased sequences, while the most frequent regular sequences in the third group are strongly correlated with transposons. Taking advantage of the generality of these mathematical procedures, we have compared the promoter database of H. sapiens with those of other species. We have found that the above-mentioned features characterize also the evolutionary content appearing in mammalian promoters, at variance with ancestral species in the phylogenetic tree, that exhibit a definitely lower level of differentiation among promoters.