Wide-ranging functions of E2F4 in transcriptional activation and repression revealed by genome-wide analysis

A Myc-microRNA network promotes exit from quiescence by suppressing the interferon response and cell-cycle arrest genes

The transition of mammalian cells from quiescence to proliferation is accompanied by the differential expression of several microRNAs (miRNAs) and transcription factors. However, the interplay between transcription factors and miRNAs in modulating gene regulatory networks involved in human cell proliferation is largely unknown. Here we show that the miRNA miR-22 promotes proliferation in primary human cells, and through a combination of Argonaute-2 immunoprecipitation and reporter assays, we identified multiple novel targets of miR-22, including several cell-cycle arrest genes that mediate the effects of the tumor-suppressor p53. In addition, we found that miR-22 suppresses interferon gene expression by directly targeting high mobility group box-1 and interferon regulatory factor (IRF)-5, preventing activation of IRF3 and NF-κB, which are activators of interferon genes. The expression of interferon genes is elevated in quiescent cells and their expression is inhibitory for cell proliferation. In addition, we find that miR-22 is activated by the transcription factor Myc when quiescent cells enter proliferation and that miR-22 inhibits the Myc transcriptional repressor MXD4, mediating a feed-forward loop to elevate Myc expression levels. Our results implicate miR-22 in downregulating the anti-proliferative p53 and interferon pathways and reveal a new transcription factor–miRNA network that regulates the transition of primary human cells from quiescence to proliferation.

Epigenetic modulation by TFII-I during embryonic stem cell differentiation

TFII-I transcription factors play an essential role during early vertebrate embryogenesis. Genome-wide mapping studies by ChIP-seq and ChIP-chip revealed that TFII-I primes multiple genomic loci in mouse embryonic stem cells and embryonic tissues. Moreover, many TFII-I-bound regions co-localize with H3K4me3/K27me3 bivalent chromatin within the promoters of lineage-specific genes. This minireview provides a summary of current knowledge regarding the function of TFII-I in epigenetic control of stem cell differentiation.

ChIPBase: a database for decoding the transcriptional regulation of long non-coding RNA and microRNA genes from ChIP-Seq data

Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) represent two classes of important non-coding RNAs in eukaryotes. Although these non-coding RNAs have been implicated in organismal development and in various human diseases, surprisingly little is known about their transcriptional regulation. Recent advances in chromatin immunoprecipitation with next-generation DNA sequencing (ChIP-Seq) have provided methods of detecting transcription factor binding sites (TFBSs) with unprecedented sensitivity. In this study, we describe ChIPBase (http://deepbase.sysu.edu.cn/chipbase/), a novel database that we have developed to facilitate the comprehensive annotation and discovery of transcription factor binding maps and transcriptional regulatory relationships of lncRNAs and miRNAs from ChIP-Seq data. The current release of ChIPBase includes high-throughput sequencing data that were generated by 543 ChIP-Seq experiments in diverse tissues and cell lines from six organisms. By analysing millions of TFBSs, we identified tens of thousands of TF-lncRNA and TF-miRNA regulatory relationships. Furthermore, two web-based servers were developed to annotate and discover transcriptional regulatory relationships of lncRNAs and miRNAs from ChIP-Seq data. In addition, we developed two genome browsers, deepView and genomeView, to provide integrated views of multidimensional data. Moreover, our web implementation supports diverse query types and the exploration of TFs, lncRNAs, miRNAs, gene ontologies and pathways.

Microsatellite polymorphisms in the EGFR, NOTCH4 and E2F4 genes and their association with breast cancer risk

BACKGROUND

The sequences of many human genes that encode proteins involved in cancer contain polymorphic microsatellites. Variations in microsatellite length may constitute risk factors in several human diseases, a possibility that has been little explored in breast cancer. Among the genes that contain polymorphic microsatellites are EGFR, NOTCH4 and E2F4. The length of some of these microsatellites has been associated with breast cancer risk.

PURPOSE AND METHODS

To determine whether the length of the microsatellites (CA)n in EGFR, (CTG)n in NOTCH4 and (AGC)n in E2F4 was associated with breast cancer risk, we genotyped these 3 microsatellites in 212 women with breast cancer and a control group of 308 women from the general population who did not have this disease.

RESULTS AND CONCLUSIONS

The allelic distribution observed for the 3 microsatellites matched that found in other white populations, with the exception of some (AGC)n alleles in E2F4, which have not been described previously. The length of (CA)n in EGFR and (CTG)n in NOTCH4 was not associated with breast cancer (OR=0.99; 95% CI 0.59-1.37; p=0.619 and OR=1.08; 95% CI 0.71-1.65; p=0.725, respectively). Short alleles (<13 repeats) of (AGC)n in E2F4 were less frequent in women with cancer than in the control sample.

Functional analysis of transcription factor binding sites in human promoters

Background

The binding of transcription factors to specific locations in the genome is integral to the orchestration of transcriptional regulation in cells. To characterize transcription factor binding site function on a large scale, we predicted and mutagenized 455 binding sites in human promoters. We carried out functional tests on these sites in four different immortalized human cell lines using transient transfections with a luciferase reporter assay, primarily for the transcription factors CTCF, GABP, GATA2, E2F, STAT, and YY1.

Results

In each cell line, between 36% and 49% of binding sites made a functional contribution to the promoter activity; the overall rate for observing function in any of the cell lines was 70%. Transcription factor binding resulted in transcriptional repression in more than a third of functional sites. When compared with predicted binding sites whose function was not experimentally verified, the functional binding sites had higher conservation and were located closer to transcriptional start sites (TSSs). Among functional sites, repressive sites tended to be located further from TSSs than were activating sites. Our data provide significant insight into the functional characteristics of YY1 binding sites, most notably the detection of distinct activating and repressing classes of YY1 binding sites. Repressing sites were located closer to, and often overlapped with, translational start sites and presented a distinctive variation on the canonical YY1 binding motif.

Conclusions

The genomic properties that we found to associate with functional TF binding sites on promoters -- conservation, TSS proximity, motifs and their variations -- point the way to improved accuracy in future TFBS predictions.

RBF binding to both canonical E2F targets and noncanonical targets depends on functional dE2F/dDP complexes

The retinoblastoma (RB) family of proteins regulate transcription. These proteins lack intrinsic DNA-binding activity but are recruited to specific genomic locations through interactions with sequence-specific DNA-binding factors. The best-known target of RB protein (pRB) is the E2F transcription factor; however, many other chromatin-associated proteins have been described that may allow RB family members to act at additional sites. To gain a perspective on the scale of E2F-dependent and E2F-independent functions, we generated genome-wide binding profiles of RBF1 and dE2F proteins in Drosophila larvae. RBF1 and dE2F2 associate with a large number of binding sites at genes with diverse biological functions. In contrast, dE2F1 was detected at a smaller set of promoters, suggesting that it overrides repression by RBF1/dE2F2 at a specific subset of targets. Approximately 15% of RBF1-bound regions lacked consensus E2F-binding motifs. To test whether RBF1 action at these sites is E2F independent, we examined dDP mutant larvae that lack any functional dE2F/dDP heterodimers. As measured by chromatin immunoprecipitation-microarray analysis (ChIP-chip), ChIP-quantitative PCR (qPCR), and cell fractionation, the stable association of RBF1 with chromatin was eliminated in dDP mutants. This requirement for dDP was seen at classic E2F-regulated promoters and at promoters that lacked canonical E2F-binding sites. These results suggest that E2F/DP complexes are essential for all genomic targeting of RBF1.

E2F7 and E2F8 promote angiogenesis through transcriptional activation of VEGFA in cooperation with HIF1

The E2F family of transcription factors plays an important role in controlling cell-cycle progression. While this is their best-known function, we report here novel functions for the newest members of the E2F family, E2F7 and E2F8 (E2F7/8). We show that simultaneous deletion of E2F7/8 in zebrafish and mice leads to severe vascular defects during embryonic development. Using a panel of transgenic zebrafish with fluorescent-labelled blood vessels, we demonstrate that E2F7/8 are essential for proper formation of blood vessels. Despite their classification as transcriptional repressors, we provide evidence for a molecular mechanism through which E2F7/8 activate the transcription of the vascular endothelial growth factor A (VEGFA), a key factor in guiding angiogenesis. We show that E2F7/8 directly bind and stimulate the VEGFA promoter independent of canonical E2F binding elements. Instead, E2F7/8 form a transcriptional complex with the hypoxia inducible factor 1 (HIF1) to stimulate VEGFA promoter activity. These results uncover an unexpected link between E2F7/8 and the HIF1-VEGFA pathway providing a molecular mechanism by which E2F7/8 control angiogenesis.

E2F mediates enhanced alternative polyadenylation in proliferation

Background

The majority of mammalian genes contain multiple poly(A) sites in their 3' UTRs. Alternative cleavage and polyadenylation are emerging as an important layer of gene regulation as they generate transcript isoforms that differ in their 3' UTRs, thereby modulating genes' response to 3' UTR-mediated regulation. Enhanced cleavage at 3' UTR proximal poly(A) sites resulting in global 3' UTR shortening was recently linked to proliferation and cancer. However, mechanisms that regulate this enhanced alternative polyadenylation are unknown.

Results

Here, we explored, on a transcriptome-wide scale, alternative polyadenylation events associated with cellular proliferation and neoplastic transformation. We applied a deep-sequencing technique for identification and quantification of poly(A) sites to two human cellular models, each examined under proliferative, arrested and transformed states. In both cell systems we observed global 3' UTR shortening associated with proliferation, a link that was markedly stronger than the association with transformation. Furthermore, we found that proliferation is also associated with enhanced cleavage at intronic poly(A) sites. Last, we found that the expression level of the set of genes that encode for 3'-end processing proteins is globally elevated in proliferation, and that E2F transcription factors contribute to this regulation.

Conclusions

Our results comprehensively identify alternative polyadenylation events associated with cellular proliferation and transformation, and demonstrate that the enhanced alternative polyadenylation in proliferative conditions results not only in global 3' UTR shortening but also in enhanced premature cleavage in introns. Our results also indicate that E2F-mediated co-transcriptional regulation of 3'-end processing genes is one of the mechanisms that links enhanced alternative polyadenylation to proliferation.

Rad18 is a transcriptional target of E2F3

The E2F family of transcription factors responds to a variety of intracellular and extracellular signals and, as such, are key regulators of cell growth, differentiation and cell death. The cellular response to DNA damage is a multistep process generally involving the initial detection of DNA damage, propagation of signals via posttranslational modifications (e.g., phosphorylation and ubiquitination) and, finally, the implementation of a response. We have previously reported that E2F3 can be induced by DNA damage, and that it plays an important role in DNA damage-induced apoptosis. Here, we demonstrate that E2F3 knockdown compromises two canonical DNA damage modification events, the ubiquitination of H2AX and PCNA. We find that the defect in these posttranscriptional modifications after E2F3 knockdown is due to reduced expression of important DNA damage responsive ubiquitin ligases. We characterized the regulation of one of these ligases, Rad18, and we demonstrated that E2F3 associates with the Rad18 promoter and directly controls its activity. Furthermore, we find that ectopic expression of Rad18 is sufficient to rescue the PCNA ubiquitination defect resulting from E2F3 knockdown. Our study reveals a novel facet of E2F3's control of the DNA damage response.

Efficient motif search in ranked lists and applications to variable gap motifs

Sequence elements, at all levels—DNA, RNA and protein, play a central role in mediating molecular recognition and thereby molecular regulation and signaling. Studies that focus on measuring and investigating sequence-based recognition make use of statistical and computational tools, including approaches to searching sequence motifs. State-of-the-art motif searching tools are limited in their coverage and ability to address large motif spaces. We develop and present statistical and algorithmic approaches that take as input ranked lists of sequences and return significant motifs. The efficiency of our approach, based on suffix trees, allows searches over motif spaces that are not covered by existing tools. This includes searching variable gap motifs—two half sites with a flexible length gap in between—and searching long motifs over large alphabets. We used our approach to analyze several high-throughput measurement data sets and report some validation results as well as novel suggested motifs and motif refinements. We suggest a refinement of the known estrogen receptor 1 motif in humans, where we observe gaps other than three nucleotides that also serve as significant recognition sites, as well as a variable length motif related to potential tyrosine phosphorylation.

E2F7 represses a network of oscillating cell cycle genes to control S-phase progression

E2F transcription factors are known to be important for timely activation of G₁/S and G₂/M genes required for cell cycle progression, but transcriptional mechanisms for deactivation of cell cycle-regulated genes are unknown. Here, we show that E2F7 is highly expressed during mid to late S-phase, occupies promoters of G₁/S-regulated genes and represses their transcription. ChIP-seq analysis revealed that E2F7 binds preferentially to genomic sites containing the TTCCCGCC motif, which closely resembles the E2F consensus site. We identified 89 target genes that carry E2F7 binding sites close to the transcriptional start site and that are directly repressed by short-term induction of E2F7. Most of these target genes are known to be activated by E2Fs and are involved in DNA replication, metabolism and DNA repair. Importantly, induction of E2F7 during G₀-G₁/S resulted in S-phase arrest and DNA damage, whereas expression of E2F7 during G₂/M failed to disturb cell cycle progression. These findings provide strong evidence that E2F7 directly controls the downswing of oscillating G₁/S genes during S-phase progression.

Cell-type specific and combinatorial usage of diverse transcription factors revealed by genome-wide binding studies in multiple human cells

Cell-type diversity is governed in part by differential gene expression programs mediated by transcription factor (TF) binding. However, there are few systematic studies of the genomic binding of different types of TFs across a wide range of human cell types, especially in relation to gene expression. In the ENCODE Project, we have identified the genomic binding locations across 11 different human cell types of CTCF, RNA Pol II (RNAPII), and MYC, three TFs with diverse roles. Our data and analysis revealed how these factors bind in relation to genomic features and shape gene expression and cell-type specificity. CTCF bound predominantly in intergenic regions while RNAPII and MYC preferentially bound to core promoter regions. CTCF sites were relatively invariant across diverse cell types, while MYC showed the greatest cell-type specificity. MYC and RNAPII co-localized at many of their binding sites and putative target genes. Cell-type specific binding sites, in particular for MYC and RNAPII, were associated with cell-type specific functions. Patterns of binding in relation to gene features were generally conserved across different cell types. RNAPII occupancy was higher over exons than adjacent introns, likely reflecting a link between transcriptional elongation and splicing. TF binding was positively correlated with the expression levels of their putative target genes, but combinatorial binding, in particular of MYC and RNAPII, was even more strongly associated with higher gene expression. These data illuminate how combinatorial binding of transcription factors in diverse cell types is associated with gene expression and cell-type specific biology.