Integrated multi-omics analysis of PBX1 in mouse adult neural stem- and progenitor cells identifies a transcriptional module that functionally links PBX1 to TCF3/4.

Developmental transcription factors act in networks, but how these networks achieve cell- and tissue specificity is still poorly understood. Here, we explored pre-B cell leukemia homeobox 1 (PBX1) in adult neurogenesis combining genomic, transcriptomic, and proteomic approaches. ChIP-seq analysis uncovered PBX1 binding to numerous genomic sites. Integration of PBX1 ChIP-seq with ATAC-seq data predicted interaction partners, which were subsequently validated by mass spectrometry. Whole transcriptome spatial RNA analysis revealed shared expression dynamics of Pbx1 and interacting factors. Among these were class I bHLH proteins TCF3 and TCF4. RNA-seq following Pbx1, Tcf3 or Tcf4 knockdown identified proliferation- and differentiation associated genes as shared targets, while sphere formation assays following knockdown argued for functional cooperativity of PBX1 and TCF3 in progenitor cell proliferation. Notably, while physiological PBX1-TCF interaction has not yet been described, chromosomal translocation resulting in genomic TCF3::PBX1 fusion characterizes a subtype of acute lymphoblastic leukemia. Introducing Pbx1 into Nalm6 cells, a pre-B cell line expressing TCF3 but lacking PBX1, upregulated the leukemogenic genes BLK and NOTCH3, arguing that functional PBX1-TCF cooperativity likely extends to hematopoiesis. Our study hence uncovers a transcriptional module orchestrating the balance between progenitor cell proliferation and differentiation in adult neurogenesis with potential implications for leukemia etiology.

Combinatorial action of NF-Y and TALE at embryonic enhancers defines distinct gene expression programs during zygotic genome activation in zebrafish.

Animal embryogenesis is initiated by maternal factors, but zygotic genome activation (ZGA) shifts regulatory control to the embryo during blastula stages. ZGA is thought to be mediated by maternally provided transcription factors (TFs), but few such TFs have been identified in vertebrates. Here we report that NF-Y and TALE TFs bind zebrafish genomic elements associated with developmental control genes already at ZGA. In particular, co-regulation by NF-Y and TALE is associated with broadly acting genes involved in transcriptional control, while regulation by either NF-Y or TALE defines genes in specific developmental processes, such that NF-Y controls a cilia gene expression program while TALE controls expression of hox genes. We also demonstrate that NF-Y and TALE-occupied genomic elements function as enhancers during embryogenesis. We conclude that combinatorial use of NF-Y and TALE at developmental enhancers permits the establishment of distinct gene expression programs at zebrafish ZGA.

p63 exerts spatio-temporal control of palatal epithelial cell fate to prevent cleft palate.

Cleft palate is a common congenital disorder that affects up to 1 in 2500 live births and results in considerable morbidity to affected individuals and their families. The aetiology of cleft palate is complex with both genetic and environmental factors implicated. Mutations in the transcription factor p63 are one of the major individual causes of cleft palate; however, the gene regulatory networks in which p63 functions remain only partially characterized. Our findings demonstrate that p63 functions as an essential regulatory molecule in the spatio-temporal control of palatal epithelial cell fate to ensure appropriate fusion of the palatal shelves. Initially, p63 induces periderm formation and controls its subsequent maintenance to prevent premature adhesion between adhesion-competent, intra-oral epithelia. Subsequently, TGFß3-induced down-regulation of p63 in the medial edge epithelia of the palatal shelves is a pre-requisite for palatal fusion by facilitating periderm migration from, and reducing the proliferative potential of, the midline epithelial seam thereby preventing cleft palate.

Regulation of the BMP Signaling-Responsive Transcriptional Network in the Drosophila Embryo.

The BMP signaling pathway has a conserved role in dorsal-ventral axis patterning during embryonic development. In Drosophila, graded BMP signaling is transduced by the Mad transcription factor and opposed by the Brinker repressor. In this study, using the Drosophila embryo as a model, we combine RNA-seq with Mad and Brinker ChIP-seq to decipher the BMP-responsive transcriptional network underpinning differentiation of the dorsal ectoderm during dorsal-ventral axis patterning. We identify multiple new BMP target genes, including positive and negative regulators of EGF signaling. Manipulation of EGF signaling levels by loss- and gain-of-function studies reveals that EGF signaling negatively regulates embryonic BMP-responsive transcription. Therefore, the BMP gene network has a self-regulating property in that it establishes a balance between its activity and that of the antagonistic EGF signaling pathway to facilitate correct patterning. In terms of BMP-dependent transcription, we identify key roles for the Zelda and Zerknüllt transcription factors in establishing the resulting expression domain, and find widespread binding of insulator proteins to the Mad and Brinker-bound genomic regions. Analysis of embryos lacking the BEAF-32 insulator protein shows reduced transcription of a peak BMP target gene and a reduction in the number of amnioserosa cells, the fate specified by peak BMP signaling. We incorporate our findings into a model for Mad-dependent activation, and discuss its relevance to BMP signal interpretation in vertebrates.

eIF2ß is critical for eIF5-mediated GDP-dissociation inhibitor activity and translational control.

In protein synthesis translation factor eIF2 binds initiator tRNA to ribosomes and facilitates start codon selection. eIF2 GDP/GTP status is regulated by eIF5 (GAP and GDI functions) and eIF2B (GEF and GDF activities), while eIF2α phosphorylation in response to diverse signals is a major point of translational control. Here we characterize a growth suppressor mutation in eIF2ß that prevents eIF5 GDI and alters cellular responses to reduced eIF2B activity, including control of GCN4 translation. By monitoring the binding of fluorescent nucleotides and initiator tRNA to purified eIF2 we show that the eIF2ß mutation does not affect intrinsic eIF2 affinities for these ligands, neither does it interfere with eIF2 binding to 43S pre-initiation complex components. Instead we show that the eIF2ß mutation prevents eIF5 GDI stabilizing nucleotide binding to eIF2, thereby altering the off-rate of GDP from eIF2•GDP/eIF5 complexes. This enables cells to grow with reduced eIF2B GEF activity but impairs activation of GCN4 targets in response to amino acid starvation. These findings provide support for the importance of eIF5 GDI activity in vivo and demonstrate that eIF2ß acts in concert with eIF5 to prevent premature release of GDP from eIF2γ and thereby ensure tight control of protein synthesis initiation.

Genome-wide occupancy links Hoxa2 to Wnt-ß-catenin signaling in mouse embryonic development.

The regulation of gene expression is central to developmental programs and largely depends on the binding of sequence-specific transcription factors with cis-regulatory elements in the genome. Hox transcription factors specify the spatial coordinates of the body axis in all animals with bilateral symmetry, but a detailed knowledge of their molecular function in instructing cell fates is lacking. Here, we used chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) to identify Hoxa2 genomic locations in a time and space when it is actively instructing embryonic development in mouse. Our data reveals that Hoxa2 has large genome coverage and potentially regulates thousands of genes. Sequence analysis of Hoxa2-bound regions identifies high occurrence of two main classes of motifs, corresponding to Hox and Pbx-Hox recognition sequences. Examination of the binding targets of Hoxa2 faithfully captures the processes regulated by Hoxa2 during embryonic development; in addition, it uncovers a large cluster of potential targets involved in the Wnt-signaling pathway. In vivo examination of canonical Wnt-ß-catenin signaling reveals activity specifically in Hoxa2 domain of expression, and this is undetectable in Hoxa2 mutant embryos. The comprehensive mapping of Hoxa2-binding sites provides a framework to study Hox regulatory networks in vertebrate developmental processes.

Aspergillus fumigatus Can Display Persistence to the Fungicidal Drug Voriconazole.

Aspergillus fumigatus is a filamentous fungus that can infect the lungs of patients with immunosuppression and/or underlying lung diseases. The mortality associated with chronic and invasive aspergillosis infections remain very high, despite availability of antifungal treatments. In the last decade, there has been a worrisome emergence and spread of resistance to the first-line antifungals, the azoles. The mortality caused by resistant isolates is even higher, and patient management is complicated as the therapeutic options are reduced. Nevertheless, treatment failure is also common in patients infected with azole-susceptible isolates, which can be due to several non-mutually exclusive reasons, such as poor drug absorption. In addition, the phenomena of tolerance or persistence, where susceptible pathogens can survive the action of an antimicrobial for extended periods, have been associated with treatment failure in bacterial infections, and their occurrence in fungal infections already proposed. Here, we demonstrate that some isolates of A. fumigatus display persistence to voriconazole. A subpopulation of the persister isolates can survive for extended periods and even grow at low rates in the presence of supra-MIC of voriconazole and seemingly other azoles. Persistence cannot be eradicated with adjuvant drugs or antifungal combinations and seemed to reduce the efficacy of treatment for certain individuals in a Galleria mellonella model of infection. Furthermore, persistence implies a distinct transcriptional profile, demonstrating that it is an active response. We propose that azole persistence might be a relevant and underestimated factor that could influence the outcome of infection in human aspergillosis. IMPORTANCE The phenomena of antibacterial tolerance and persistence, where pathogenic microbes can survive for extended periods in the presence of cidal drug concentrations, have received significant attention in the last decade. Several mechanisms of action have been elucidated, and their relevance for treatment failure in bacterial infections demonstrated. In contrast, our knowledge of antifungal tolerance and, in particular, persistence is still very limited. In this study, we have characterized the response of the prominent fungal pathogen Aspergillus fumigatus to the first-line therapy antifungal voriconazole. We comprehensively show that some isolates display persistence to this fungicidal antifungal and propose various potential mechanisms of action. In addition, using an alternative model of infection, we provide initial evidence to suggest that persistence may cause treatment failure in some individuals. Therefore, we propose that azole persistence is an important factor to consider and further investigate in A. fumigatus.

Elucidation of the ELK1 target gene network reveals a role in the coordinate regulation of core components of the gene regulation machinery.

Transcription factors play an important role in orchestrating the activation of specific networks of genes through targeting their proximal promoter and distal enhancer regions. However, it is unclear how the specificity of downstream responses is maintained by individual members of transcription-factor families and, in most cases, what their target repertoire is. We have used ChIP-chip analysis to identify the target genes of the ETS-domain transcription factor ELK1. Two distinct modes of ELK1 target gene selection are identified; the first involves redundant promoter binding with other ETS-domain family members; the second occurs through combinatorial binding with a second transcription factor SRF, which specifies a unique group of target genes. One of the most prominent groups of genes forming the ELK1 target network includes classes involved in core gene expression control, namely, components of the basal transcriptional machinery, the spliceosome and the ribosome. Amongst the set of genes encoding the basal transcription machinery components, are a functionally linked subset of GTFs and TAFs. Our study, therefore, reveals an unsuspected level of coordinate regulation of components of the core gene expression control machinery and also identifies two different modes of promoter targeting through binding with a second transcription factor or redundant binding with other ETS-domain family members.

Proneural bHLH and Brn proteins coregulate a neurogenic program through cooperative binding to a conserved DNA motif.

Proneural proteins play a central role in vertebrate neurogenesis, but little is known of the genes that they regulate and of the factors that interact with proneural proteins to activate a neurogenic program. Here, we demonstrate that the proneural protein Mash1 and the POU proteins Brn1 and Brn2 interact on the promoter of the Notch ligand Delta1 and synergistically activate Delta1 transcription, a key step in neurogenesis. Overexpression experiments in vivo indicate that Brn2, like Mash1, regulates additional aspects of neurogenesis, including the division of progenitors and the differentiation and migration of neurons. We identify by in silico screening a number of additional candidate target genes, which are recognized by Mash1 and Brn proteins through a DNA-binding motif similar to that found in the Delta1 gene and present a broad range of activities. We thus propose that Mash1 synergizes with Brn factors to regulate multiple steps of neurogenesis.

ID1 promotes expansion and survival of primary erythroid cells and is a target of JAK2V617F-STAT5 signaling.

The discovery of JAK2V617F as an acquired mutation in the majority of patients with myeloproliferative disorders (MPDs) and the key role of the JAK2-STAT5 signaling cascade in normal hematopoiesis has focused attention on the downstream transcriptional targets of STAT5. Despite evidence of its vital role in normal erythropoiesis and its ability to recapitulate many of the features of myeloid malignancies, including the MPDs, few functionally validated targets of STAT5 have been described. Here we used a combination of comparative genomics and chromatin immunoprecipitation assays to identify ID1 as a novel target of the JAK2-STAT5 signaling axis in erythroid cells. STAT5 binds and transactivates a downstream enhancer of ID1, and ID1 expression levels correlate with the JAK2V617F mutation in both retrovirally transfected fetal liver cells and polycythemia vera patients. Knockdown and overexpression studies in a well-characterized erythroid differentiation assay from primary murine fetal liver cells demonstrated a survival-promoting action of ID1. This hitherto unrecognized function implicates ID1 in the expansion of erythroblasts during terminal differentiation and suggests that ID1 plays an important role in the pathogenesis of polycythemia vera. Furthermore, our findings contribute to an increasing body of evidence implicating ID proteins in a wider range of cellular functions than initially appreciated.

The transcriptional program controlled by the stem cell leukemia gene Scl/Tal1 during early embryonic hematopoietic development.

The basic helix-loop-helix transcription factor Scl/Tal1 controls the development and subsequent differentiation of hematopoietic stem cells (HSCs). However, because few Scl target genes have been validated to date, the underlying mechanisms have remained largely unknown. In this study, we have used ChIP-Seq technology (coupling chromatin immunoprecipitation with deep sequencing) to generate a genome-wide catalog of Scl-binding events in a stem/progenitor cell line, followed by validation using primary fetal liver cells and comprehensive transgenic mouse assays. Transgenic analysis provided in vivo validation of multiple new direct Scl target genes and allowed us to reconstruct an in vivo validated network consisting of 17 factors and their respective regulatory elements. By coupling ChIP-Seq in model cell lines with in vivo transgenic validation and sophisticated bioinformatic analysis, we have identified a widely applicable strategy for the reconstruction of stem cell regulatory networks in which biologic material is otherwise limiting. Moreover, in addition to revealing multiple previously unrecognized links to known HSC regulators, as well as novel links to genes not previously implicated in HSC function, comprehensive transgenic analysis of regulatory elements provided substantial new insights into the transcriptional control of several important hematopoietic regulators, including Cbfa2t3h/Eto2, Cebpe, Nfe2, Zfpm1/Fog1, Erg, Mafk, Gfi1b, and Myb.

Expression of the leukemia oncogene Lmo2 is controlled by an array of tissue-specific elements dispersed over 100 kb and bound by Tal1/Lmo2, Ets, and Gata factors.

The Lmo2 gene encodes a transcriptional cofactor critical for the development of hematopoietic stem cells. Ectopic LMO2 expression causes leukemia in T-cell acute lymphoblastic leukemia (T-ALL) patients and severe combined immunodeficiency patients undergoing retroviral gene therapy. Tightly controlled Lmo2 expression is therefore essential, yet no comprehensive analysis of Lmo2 regulation has been published so far. By comparative genomics, we identified 17 highly conserved noncoding elements, 9 of which revealed specific acetylation marks in chromatin-immunoprecipitation and microarray (ChIP-chip) assays performed across 250 kb of the Lmo2 locus in 11 cell types covering different stages of hematopoietic differentiation. All candidate regulatory regions were tested in transgenic mice. An extended LMO2 proximal promoter fragment displayed strong endothelial activity, while the distal promoter showed weak forebrain activity. Eight of the 15 distal candidate elements functioned as enhancers, which together recapitulated the full expression pattern of Lmo2, directing expression to endothelium, hematopoietic cells, tail, and forebrain. Interestingly, distinct combinations of specific distal regulatory elements were required to extend endothelial activity of the LMO2 promoter to yolk sac or fetal liver hematopoietic cells. Finally, Sfpi1/Pu.1, Fli1, Gata2, Tal1/Scl, and Lmo2 were shown to bind to and transactivate Lmo2 hematopoietic enhancers, thus identifying key upstream regulators and positioning Lmo2 within hematopoietic regulatory networks.

Overlapping promoter targeting by Elk-1 and other divergent ETS-domain transcription factor family members.

ETS-domain transcription factors play important roles in controlling gene expression in a variety of different contexts; however, these proteins bind to very similar sites and it is unclear how in vivo specificity is achieved. In silico analysis is unlikely to reveal specific targets for individual family members and direct experimental approaches are therefore required. Here, we take advantage of an inducible dominant-negative expression system to identify a group of novel target genes for the ETS-domain transcription factor Elk-1. Elk-1 is thought to mainly function through cooperation with a second transcription factor SRF, but the targets we identify are largely SRF-independent. Furthermore, we demonstrate that there is a high degree of overlapping, cell type-specific, target gene binding by Elk-1 and other ETS-domain transcription factors. Our results are therefore consistent with the notion that there is a high degree of functional redundancy in target gene regulation by ETS-domain transcription factors in addition to the specific target gene regulation that can be dictated through heterotypic interactions exemplified by the Elk-1-SRF complex.

ORegAnno: an open-access community-driven resource for regulatory annotation.

ORegAnno is an open-source, open-access database and literature curation system for community-based annotation of experimentally identified DNA regulatory regions, transcription factor binding sites and regulatory variants. The current release comprises 30 145 records curated from 922 publications and describing regulatory sequences for over 3853 genes and 465 transcription factors from 19 species. A new feature called the 'publication queue' allows users to input relevant papers from scientific literature as targets for annotation. The queue contains 4438 gene regulation papers entered by experts and another 54 351 identified by text-mining methods. Users can enter or 'check out' papers from the queue for manual curation using a series of user-friendly annotation pages. A typical record entry consists of species, sequence type, sequence, target gene, binding factor, experimental outcome and one or more lines of experimental evidence. An evidence ontology was developed to describe and categorize these experiments. Records are cross-referenced to Ensembl or Entrez gene identifiers, PubMed and dbSNP and can be visualized in the Ensembl or UCSC genome browsers. All data are freely available through search pages, XML data dumps or web services at: http://www.oreganno.org.

Hypoxia regulates GR function through multiple mechanisms involving microRNAs 103 and 107.

Glucocorticoids (Gcs) potently inhibit inflammation, and regulate liver energy metabolism, often acting in a hypoxic environment. We now show hypoxic conditions open a specific GR cistrome, and prevent access of GR to part of the normoxic GR cistrome. Motif analysis identified enrichment of KLF4 binding sites beneath those peaks of GR binding exclusive to normoxia, implicating KLF4 as a pioneer, or co-factor under these conditions. Hypoxia reduced KLF4 expression, however, knockdown of KLF4 did not impair GR recruitment. KLF4 is a known target of microRNAs 103 and 107, both of which are induced by hypoxia. Expression of mimics to either microRNA103, or microRNA107 inhibited GR transactivation of normoxic target genes, thereby replicating the hypoxic effect. Therefore, studies in hypoxia reveal that microRNAs 103 and 107 are potent regulators of GR function. We have now identified a new pathway linking hypoxia through microRNAs 103 and 107 to regulation of GR function.

Vulnerability of drug-resistant EML4-ALK rearranged lung cancer to transcriptional inhibition.

A subset of lung adenocarcinomas is driven by the EML4-ALK translocation. Even though ALK inhibitors in the clinic lead to excellent initial responses, acquired resistance to these inhibitors due to on-target mutations or parallel pathway alterations is a major clinical challenge. Exploring these mechanisms of resistance, we found that EML4-ALK cells parental or resistant to crizotinib, ceritinib or alectinib are remarkably sensitive to inhibition of CDK7/12 with THZ1 and CDK9 with alvocidib or dinaciclib. These compounds robustly induce apoptosis through transcriptional inhibition and downregulation of anti-apoptotic genes. Importantly, alvocidib reduced tumour progression in xenograft mouse models. In summary, our study takes advantage of the transcriptional addiction hypothesis to propose a new treatment strategy for a subset of patients with acquired resistance to first-, second- and third-generation ALK inhibitors.

A living biobank of ovarian cancer ex vivo models reveals profound mitotic heterogeneity.

High-grade serous ovarian carcinoma is characterised by TP53 mutation and extensive chromosome instability (CIN). Because our understanding of CIN mechanisms is based largely on analysing established cell lines, we developed a workflow for generating ex vivo cultures from patient biopsies to provide models that support interrogation of CIN mechanisms in cells not extensively cultured in vitro. Here, we describe a "living biobank" of ovarian cancer models with extensive replicative capacity, derived from both ascites and solid biopsies. Fifteen models are characterised by p53 profiling, exome sequencing and transcriptomics, and karyotyped using single-cell whole-genome sequencing. Time-lapse microscopy reveals catastrophic and highly heterogeneous mitoses, suggesting that analysis of established cell lines probably underestimates mitotic dysfunction in advanced human cancers. Drug profiling reveals cisplatin sensitivities consistent with patient responses, demonstrating that this workflow has potential to generate personalized avatars with advantages over current pre-clinical models and the potential to guide clinical decision making.

CoMoDis: composite motif discovery in mammalian genomes.

Specificity of mammalian gene regulatory regions is achieved to a large extent through the combinatorial binding of sets of transcription factors to distinct binding sites, discrete combinations of which are often referred to as regulatory modules. Identification and subsequent characterization of gene regulatory modules will be a key step in assembling transcriptional regulatory networks from gene expression profiling data, with the ultimate goal of unravelling the regulatory codes that govern gene expression in various cell types. Here we describe the new bioinformatics tool, Composite Motif Discovery (CoMoDis), which streamlines computational identification of novel regulatory modules starting from a single seed motif. Seed motifs represent binding sites conserved across mammalian species. CoMoDis facilitates novel motif discovery by automating the extraction of DNA sequences flanking seed motifs and streamlining downstream motif discovery using a variety of tools, including several that utilize phylogenetic conservation criteria. CoMoDis is available at http://hscl.cimr.cam.ac.uk/CoMoDis_portal.html.

TFBScluster web server for the identification of mammalian composite regulatory elements.

Identification of transcriptional regulatory elements represents a critical step in our ability to reconstruct transcriptional regulatory networks from gene expression profiling datasets. To facilitate computational identification of candidate gene regulatory elements from whole genome sequences, we have developed the TFBScluster web server that integrates several tools for the genome-wide identification and subsequent characterization of transcription factor binding site clusters that are conserved in multiple mammalian species. Either the human or mouse genomes can be used as the reference sequence with direct links from the search results to the ENSEMBL and UCSC genome browsers. Moreover, TFBScluster provides seamless integration of transcription factor binding site searches with genome annotation and gene expression profiling data, to allow prioritising computational predictions for subsequent experimental validation. TFBScluster is publicly available at http://hscl.cimr.cam.ac.uk/TFBScluster_genome_portal.html.

Identification of the REST regulon reveals extensive transposable element-mediated binding site duplication.

The genome-wide mapping of gene-regulatory motifs remains a major goal that will facilitate the modelling of gene-regulatory networks and their evolution. The repressor element 1 is a long, conserved transcription factor-binding site which recruits the transcriptional repressor REST to numerous neuron-specific target genes. REST plays important roles in multiple biological processes and disease states. To map RE1 sites and target genes, we created a position specific scoring matrix representing the RE1 and used it to search the human and mouse genomes. We identified 1301 and 997 RE1s inhuman and mouse genomes, respectively, of which >40% are novel. By employing an ontological analysis we show that REST target genes are significantly enriched in a number of functional classes. Taking the novel REST target gene CACNA1A as an experimental model, we show that it can be regulated by multiple RE1s of different binding affinities, which are only partially conserved between human and mouse. A novel BLAST methodology indicated that many RE1s belong to closely related families. Most of these sequences are associated with transposable elements, leading us to propose that transposon-mediated duplication and insertion of RE1s has led to the acquisition of novel target genes by REST during evolution.