The forkhead transcription factor FOXK2 premarks lineage-specific genes in human embryonic stem cells for activation during differentiation.

Enhancers play important roles in controlling gene expression in a choreographed spatial and temporal manner during development. However, it is unclear how these regulatory regions are established during differentiation. Here we investigated the genome-wide binding profile of the forkhead transcription factor FOXK2 in human embryonic stem cells (ESCs) and downstream cell types. This transcription factor is bound to thousands of regulatory regions in human ESCs, and binding at many sites is maintained as cells differentiate to mesendodermal and neural precursor cell (NPC) types, alongside the emergence of new binding regions. FOXK2 binding is generally associated with active histone marks in any given cell type. Furthermore newly acquired, or retained FOXK2 binding regions show elevated levels of activating histone marks following differentiation to NPCs. In keeping with this association with activating marks, we demonstrate a role for FOXK transcription factors in gene activation during NPC differentiation. FOXK2 occupancy in ESCs is therefore an early mark for delineating the regulatory regions, which become activated in later lineages.

Complexities in the role of acetylation dynamics in modifying inducible gene activation parameters.

High levels of histone acetylation are associated with the regulatory elements of active genes, suggesting a link between acetylation and gene activation. We revisited this model, in the context of EGF-inducible gene expression and found that rather than a simple unifying model, there are two broad classes of genes; one in which high lysine acetylation activity is required for efficient gene activation, and a second group where the opposite occurs and high acetylation activity is inhibitory. We examined the latter class in more detail using EGR2 as a model gene and found that lysine acetylation levels are critical for several activation parameters, including the timing of expression onset, and overall amplitudes of the transcriptional response. In contrast, DUSP1 responds in the canonical manner and its transcriptional activity is promoted by acetylation. Single cell approaches demonstrate heterogenous activation kinetics of a given gene in response to EGF stimulation. Acetylation levels modify these heterogenous patterns and influence both allele activation frequencies and overall expression profile parameters. Our data therefore point to a complex interplay between acetylation equilibria and target gene induction where acetylation level thresholds are an important determinant of transcriptional induction dynamics that are sensed in a gene-specific manner.

Identification of a primitive intestinal transcription factor network shared between esophageal adenocarcinoma and its precancerous precursor state.

Esophageal adenocarcinoma (EAC) is one of the most frequent causes of cancer death, and yet compared to other common cancers, we know relatively little about the molecular composition of this tumor type. To further our understanding of this cancer, we have used open chromatin profiling to decipher the transcriptional regulatory networks that are operational in EAC. We have uncovered a transcription factor network that is usually found in primitive intestinal cells during embryonic development, centered on HNF4A and GATA6. These transcription factors work together to control the EAC transcriptome. We show that this network is activated in Barrett's esophagus, the putative precursor state to EAC, thereby providing novel molecular evidence in support of stepwise malignant transition. Furthermore, we show that HNF4A alone is sufficient to drive chromatin opening and activation of a Barrett's-like chromatin signature when expressed in normal human epithelial cells. Collectively, these data provide a new way to categorize EAC at a genome scale and implicate HNF4A activation as a potential pivotal event in its malignant transition from healthy cells.

Classifying cells with Scasat, a single-cell ATAC-seq analysis tool.

ATAC-seq is a recently developed method to identify the areas of open chromatin in a cell. These regions usually correspond to active regulatory elements and their location profile is unique to a given cell type. When done at single-cell resolution, ATAC-seq provides an insight into the cell-to-cell variability that emerges from otherwise identical DNA sequences by identifying the variability in the genomic location of open chromatin sites in each of the cells. This paper presents Scasat (single-cell ATAC-seq analysis tool), a complete pipeline to process scATAC-seq data with simple steps. Scasat treats the data as binary and applies statistical methods that are especially suitable for binary data. The pipeline is developed in a Jupyter notebook environment that holds the executable code along with the necessary description and results. It is robust, flexible, interactive and easy to extend. Within Scasat we developed a novel differential accessibility analysis method based on information gain to identify the peaks that are unique to a cell. The results from Scasat showed that open chromatin locations corresponding to potential regulatory elements can account for cellular heterogeneity and can identify regulatory regions that separates cells from a complex population.

Open chromatin profiling identifies AP1 as a transcriptional regulator in oesophageal adenocarcinoma.

Oesophageal adenocarcinoma (OAC) is one of the ten most prevalent forms of cancer and is showing a rapid increase in incidence and yet exhibits poor survival rates. Compared to many other common cancers, the molecular changes that occur in this disease are relatively poorly understood. However, genes encoding chromatin remodeling enzymes are frequently mutated in OAC. This is consistent with the emerging concept that cancer cells exhibit reprogramming of their chromatin environment which leads to subsequent changes in their transcriptional profile. Here, we have used ATAC-seq to interrogate the chromatin changes that occur in OAC using both cell lines and patient-derived material. We demonstrate that there are substantial changes in the regulatory chromatin environment in the cancer cells and using this data we have uncovered an important role for ETS and AP1 transcription factors in driving the changes in gene expression found in OAC cells.

A competitive transcription factor binding mechanism determines the timing of late cell cycle-dependent gene expression.

Transcriptional control is exerted by the antagonistic activities of activator and repressor proteins. In Saccharomyces cerevisiae, transcription factor complexes containing the MADS box protein Mcm1p are key regulators of cell cycle-dependent transcription at both the G2/M and M/G1 transitions. The homeodomain repressor protein Yox1p acts in a complex with Mcm1p to control the timing of gene expression. Here, we show that Yox1p interacts with Mcm1p through a motif located N terminally to its homeodomain. Yox1p functions as a transcriptional repressor by competing with the forkhead transcription activator protein Fkh2p for binding to Mcm1p through protein-protein interactions at promoters of a subset of Mcm1p-regulated genes. Importantly, this competition is not through binding the same DNA site that is commonly observed. Thus, this study describes a different mechanism for determining the timing of cell cycle-dependent gene expression that involves competition between short peptide motifs in repressor and activator proteins for interaction with a common binding partner.

Genome-wide binding studies reveal DNA binding specificity mechanisms and functional interplay amongst Forkhead transcription factors.

Transcription factors belonging to the same transcription factor families contain very similar DNA binding domains and hence have the potential to bind to related DNA sequences. However, subtle differences in binding specificities can be detected in vitro with the potential to direct specific responses in vivo. Here, we have examined the binding properties of three Forkhead (FOX) transcription factors, FOXK2, FOXO3 and FOXJ3 in vivo. Extensive overlap in chromatin binding is observed, although underlying differential DNA binding specificity can dictate the recruitment of FOXK2 and FOXJ3 to chromatin. However, functionally, FOXO3-dependent gene regulation is generally mediated not through uniquely bound regions but through regions occupied by both FOXK2 and FOXO3 where both factors play a regulatory role. Our data point to a model whereby FOX transcription factors control gene expression through dynamically binding and generating partial occupancy of the same site rather than mutually exclusive binding derived by stable binding of individual FOX proteins.

Screen for multi-SUMO-binding proteins reveals a multi-SIM-binding mechanism for recruitment of the transcriptional regulator ZMYM2 to chromatin.

Protein SUMOylation has emerged as an important regulatory event, particularly in nuclear processes such as transcriptional control and DNA repair. In this context, small ubiquitin-like modifier (SUMO) often provides a binding platform for the recruitment of proteins via their SUMO-interacting motifs (SIMs). Recent discoveries point to an important role for multivalent SUMO binding through multiple SIMs in the binding partner as exemplified by poly-SUMOylation acting as a binding platform for ubiquitin E3 ligases such as ring finger protein 4. Here, we have investigated whether other types of protein are recruited through multivalent SUMO interactions. We have identified dozens of proteins that bind to multi-SUMO platforms, thereby uncovering a complex potential regulatory network. Multi-SUMO binding is mediated through multi-SIM modules, and the functional importance of these interactions is demonstrated for the transcriptional corepressor ZMYM2/ZNF198 where its multi-SUMO-binding activity is required for its recruitment to chromatin.

Jun-Mediated Changes in Cell Adhesion Contribute to Mouse Embryonic Stem Cell Exit from Ground State Pluripotency.

Embryonic stem cells (ESC) are able to give rise to any somatic cell type. A lot is known about how ESC pluripotency is maintained, but comparatively less is known about how differentiation is promoted. Cell fate decisions are regulated by interactions between signaling and transcriptional networks. Recent studies have shown that the overexpression or downregulation of the transcription factor Jun can affect the ESC fate. Here we have focussed on the role of the Jun in the exit of mouse ESCs from ground state pluripotency and the onset of early differentiation. Transcriptomic analysis of differentiating ESCs reveals that Jun is required to upregulate a programme of genes associated with cell adhesion as ESCs exit the pluripotent ground state. Several of these Jun-regulated genes are shown to be required for efficient adhesion. Importantly this adhesion is required for the timely regulated exit of ESCs from ground state pluripotency and the onset of early differentiation events. Stem Cells 2016;34:1213-1224.

JNK-associated Leucine Zipper Protein Functions as a Docking Platform for Polo-like Kinase 1 and Regulation of the Associating Transcription Factor Forkhead Box Protein K1.

JLP (JNK-associated leucine zipper protein) is a scaffolding protein that interacts with various signaling proteins associated with coordinated regulation of cellular process such as endocytosis, motility, neurite outgrowth, cell proliferation, and apoptosis. Here we identified PLK1 (Polo-like kinase 1) as a novel interaction partner of JLP through mass spectrometric approaches. Our results indicate that JLP is phospho-primed by PLK1 on Thr-351, which is recognized by the Polo box domain of PLK1 leading to phosphorylation of JLP at additional sites. Stable isotope labeling by amino acids in cell culture and quantitative LC-MS/MS analysis was performed to identify PLK1-dependent JLP-interacting proteins. Treatment of cells with the PLK1 kinase inhibitor BI2536 suppressed binding of the Forkhead box protein K1 (FOXK1) transcriptional repressor to JLP. JLP was found to interact with PLK1 and FOXK1 during mitosis. Moreover, knockdown of PLK1 affected the interaction between JLP and FOXK1. FOXK1 is a known transcriptional repressor of the CDK inhibitor p21/WAF1, and knockdown of JLP resulted in increased FOXK1 protein levels and a reduction of p21 transcript levels. Our results suggest a novel mechanism by which FOXK1 protein levels and activity are regulated by associating with JLP and PLK1.

MAP kinase-mediated c-fos regulation relies on a histone acetylation relay switch.

Gene activation is often associated with high levels of histone acetylation. Enhanced acetylation levels can promote the recruitment of further chromatin modifying complexes or the basal transcription machinery. Here, we have studied MAP kinase-mediated upregulation of c-fos and uncover a role for histone acetylation in promoting the recruitment of a second transcription factor, NFI. MAP kinase signaling to Elk-1 enhances the net histone acetylase activity associated with the c-fos promoter, which leads to changes in the acetylation state and structure of a promoter-proximal nucleosome, which allows NFI binding. Binding of NFI provides a permissive state for the recruitment of basal machinery and subsequent promoter activation. Our results provide insights into how MAP kinase signaling promotes inducible gene expression; phosphorylation of recipient transcription factors (primary effectors) triggers a HAT relay switch, which facilitates the recruitment of additional transcription factors (secondary effectors) through alteration of the local nucleosomal structure.

The forkhead transcription factor FOXK2 acts as a chromatin targeting factor for the BAP1-containing histone deubiquitinase complex.

There are numerous forkhead transcription factors in mammalian cells but we know little about the molecular functions of the majority of these. FOXK2 is a ubiquitously expressed family member suggesting an important function across multiple cell types. Here, we show that FOXK2 binds to the SIN3A and PR-DUB complexes. The PR-DUB complex contains the important tumour suppressor protein, the deubiquitinase BAP1. FOXK2 recruits BAP1 to DNA, promotes local histone deubiquitination and causes changes in target gene activity. Our results therefore provide an important link between BAP1 and the transcription factor FOXK2 and demonstrate how BAP1 can be recruited to specific regulatory loci.

Deregulation of the FOXM1 target gene network and its coregulatory partners in oesophageal adenocarcinoma.

BACKGROUND: Survival rates for oesophageal adenocarcinoma (OAC) remain disappointingly poor and current conventional treatment modalities have minimal impact on long-term survival. This is partly due to a lack of understanding of the molecular changes that occur in this disease. Previous studies have indicated that the transcription factor FOXM1 is commonly upregulated in this cancer type but the impact of this overexpression on gene expression in the context of OAC is largely unknown. FOXM1 does not function alone but works alongside the antagonistically-functioning co-regulatory MMB and DREAM complexes. METHODS: To establish how FOXM1 affects gene expression in OAC we have identified the FOXM1 target gene network in OAC-derived cells using ChIP-seq and determined the expression of both its coregulatory partners and members of this target gene network in OAC by digital transcript counting using the Nanostring gene expression assay. RESULTS: We find co-upregulation of FOXM1 with its target gene network in OAC. Furthermore, we find changes in the expression of its coregulatory partners, including co-upregulation of LIN9 and, surprisingly, reduced expression of LIN54. Mechanistically, we identify LIN9 as the direct binding partner for FOXM1 in the MMB complex. In the context of OAC, both coregulator (eg LIN54) and target gene (eg UHRF1) expression levels are predictive of disease stage. CONCLUSIONS: Together our data demonstrate that there are global changes to the FOXM1 regulatory network in OAC and the expression of components of this network help predict cancer prognosis.

PARP1 orchestrates variant histone exchange in signal-mediated transcriptional activation.

Transcriptional activation is accompanied by multiple molecular events that remodel the local chromatin environment in promoter regions. These molecular events are often orchestrated in response to the activation of signalling pathways, as exemplified by the response of immediate early genes such as FOS to ERK MAP kinase signalling. Here, we demonstrate that inducible NFI recruitment permits PARP1 binding to the FOS promoter by a mutually reinforcing loop. PARP1 and its poly(ADP-ribosyl)ation activity are required for maintaining FOS activation kinetics. We also show that the histone variant H2A.Z associates with the FOS promoter and acts in a transcription-suppressive manner. However, in response to ERK pathway signalling, H2A.Z is replaced by H2A; PARP1 activity is required to promote this exchange. Thus, our work has revealed an additional facet of PARP1 function in promoting dynamic remodelling of promoter-associated nucleosomes to allow transcriptional activation in response to cellular signalling.

ELK1 uses different DNA binding modes to regulate functionally distinct classes of target genes.

Eukaryotic transcription factors are grouped into families and, due to their similar DNA binding domains, often have the potential to bind to the same genomic regions. This can lead to redundancy at the level of DNA binding, and mechanisms are required to generate specific functional outcomes that enable distinct gene expression programmes to be controlled by a particular transcription factor. Here we used ChIP-seq to uncover two distinct binding modes for the ETS transcription factor ELK1. In one mode, other ETS transcription factors can bind regulatory regions in a redundant fashion; in the second, ELK1 binds in a unique fashion to another set of genomic targets. Each binding mode is associated with different binding site features and also distinct regulatory outcomes. Furthermore, the type of binding mode also determines the control of functionally distinct subclasses of genes and hence the phenotypic response elicited. This is demonstrated for the unique binding mode where a novel role for ELK1 in controlling cell migration is revealed. We have therefore uncovered an unexpected link between the type of binding mode employed by a transcription factor, the subsequent gene regulatory mechanisms used, and the functional categories of target genes controlled.

A genome-wide RNAi screen reveals MAP kinase phosphatases as key ERK pathway regulators during embryonic stem cell differentiation.

Embryonic stem cells and induced pluripotent stem cells represent potentially important therapeutic agents in regenerative medicine. Complex interlinked transcriptional and signaling networks control the fate of these cells towards maintenance of pluripotency or differentiation. In this study we have focused on how mouse embryonic stem cells begin to differentiate and lose pluripotency and, in particular, the role that the ERK MAP kinase and GSK3 signaling pathways play in this process. Through a genome-wide siRNA screen we have identified more than 400 genes involved in loss of pluripotency and promoting the onset of differentiation. These genes were functionally associated with the ERK and/or GSK3 pathways, providing an important resource for studying the roles of these pathways in controlling escape from the pluripotent ground state. More detailed analysis identified MAP kinase phosphatases as a focal point of regulation and demonstrated an important role for these enzymes in controlling ERK activation kinetics and subsequently determining early embryonic stem cell fate decisions.

MBD5 and MBD6 interact with the human PR-DUB complex through their methyl-CpG-binding domain.

MBD5 and MBD6 are two members of the methyl-CpG-binding domain (MBD) family of proteins that are poorly characterized. Studies performed thus far have failed to show binding of the MBD5 and MBD6 MBD to methylated DNA. Here, we show that both MBD5 and MBD6 interact with the mammalian PR-DUB Polycomb protein complex in a mutually exclusive manner. Strikingly, the MBD of MBD5 and MBD6 is both necessary and sufficient to mediate this interaction. Chromatin immunoprecipitation analyses reveal that MBD6 and FOXK2/PR-DUB share a subset of genomic target genes, suggesting a functional interaction in vivo. Finally, we show that MBD6, but not MBD5, is recruited to sites of DNA damage in a PR-DUB independent manner. Our study thus implies a shared function for MBD5 and MBD6 through an interaction with PR-DUB, as well as an MBD6-specific recruitment to sites of DNA damage.

WDR5, ASH2L, and RBBP5 control the efficiency of FOS transcript processing.

H3K4 trimethylation is strongly associated with active transcription. The deposition of this mark is catalyzed by SET-domain methyltransferases, which consist of a subcomplex containing WDR5, ASH2L, and RBBP5 (the WAR subcomplex); a catalytic SET-domain protein; and additional complexspecific subunits. The ERK MAPK pathway also plays an important role in gene regulation via phosphorylation of transcription factors, co-regulators, or histone modifier complexes. However, the potential interactions between these two pathways remain largely unexplored. We investigated their potential interplay in terms of the regulation of the immediate early gene (IEG) regulatory network. We found that depletion of components of the WAR subcomplex led to increased levels of unspliced transcripts of IEGs that did not necessarily reflect changes in their mature transcripts. This occurs in a manner independent from changes in the H3K4me3 levels at the promoter region. We focused on FOS and found that the depletion of WAR subcomplex components affected the efficiency of FOS transcript processing. Our findings show a new aspect of WAR subcomplex function in coordinating active transcription with efficient pre-mRNA processing.

Nuclear localization of heparanase 2 (Hpa2) attenuates breast carcinoma growth and metastasis.

Unlike the intense research effort devoted to exploring the significance of heparanase in cancer, very little attention was given to Hpa2, a close homolog of heparanase. Here, we explored the role of Hpa2 in breast cancer. Unexpectedly, we found that patients endowed with high levels of Hpa2 exhibited a higher incidence of tumor metastasis and survived less than patients with low levels of Hpa2. Immunohistochemical examination revealed that in normal breast tissue, Hpa2 localizes primarily in the cell nucleus. In striking contrast, in breast carcinoma, Hpa2 expression is not only decreased but also loses its nuclear localization and appears diffuse in the cell cytoplasm. Importantly, breast cancer patients in which nuclear localization of Hpa2 is retained exhibited reduced lymph-node metastasis, suggesting that nuclear localization of Hpa2 plays a protective role in breast cancer progression. To examine this possibility, we engineered a gene construct that directs Hpa2 to the cell nucleus (Hpa2-Nuc). Notably, overexpression of Hpa2 in breast carcinoma cells resulted in bigger tumors, whereas targeting Hpa2 to the cell nucleus attenuated tumor growth and tumor metastasis. RNAseq analysis was performed to reveal differentially expressed genes (DEG) in Hpa2-Nuc tumors vs. control. The analysis revealed, among others, decreased expression of genes associated with the hallmark of Kras, beta-catenin, and TNF-alpha (via NFkB) signaling. Our results imply that nuclear localization of Hpa2 prominently regulates gene transcription, resulting in attenuation of breast tumorigenesis. Thus, nuclear Hpa2 may be used as a predictive parameter in personalized medicine for breast cancer patients.