Selective deployment of transcription factor paralogs with submaximal strength facilitates gene regulation in the immune system.

In multicellular organisms, duplicated genes can diverge through tissue-specific gene expression patterns, as exemplified by highly regulated expression of RUNX transcription factor paralogs with apparent functional redundancy. Here we asked what cell-type-specific biologies might be supported by the selective expression of RUNX paralogs during Langerhans cell and inducible regulatory T cell differentiation. We uncovered functional nonequivalence between RUNX paralogs. Selective expression of native paralogs allowed integration of transcription factor activity with extrinsic signals, while non-native paralogs enforced differentiation even in the absence of exogenous inducers. DNA binding affinity was controlled by divergent amino acids within the otherwise highly conserved RUNT domain and evolutionary reconstruction suggested convergence of RUNT domain residues toward submaximal strength. Hence, the selective expression of gene duplicates in specialized cell types can synergize with the acquisition of functional differences to enable appropriate gene expression, lineage choice and differentiation in the mammalian immune system.

Principles of dimer-specific gene regulation revealed by a comprehensive characterization of NF-κB family DNA binding.

The unique DNA-binding properties of distinct NF-κB dimers influence the selective regulation of NF-κB target genes. To more thoroughly investigate these dimer-specific differences, we combined protein-binding microarrays and surface plasmon resonance to evaluate DNA sites recognized by eight different NF-κB dimers. We observed three distinct binding-specificity classes and clarified mechanisms by which dimers might regulate distinct sets of genes. We identified many new nontraditional NF-κB binding site (κB site) sequences and highlight the plasticity of NF-κB dimers in recognizing κB sites with a single consensus half-site. This study provides a database that can be used in efforts to identify NF-κB target sites and uncover gene regulatory circuitry.

Diversification of transcription factor paralogs via noncanonical modularity in C2H2 zinc finger DNA binding.

A major challenge in obtaining a full molecular description of evolutionary adaptation is to characterize how transcription factor (TF) DNA-binding specificity can change. To identify mechanisms of TF diversification, we performed detailed comparisons of yeast C2H2 ZF proteins with identical canonical recognition residues that are expected to bind the same DNA sequences. Unexpectedly, we found that ZF proteins can adapt to recognize new binding sites in a modular fashion whereby binding to common core sites remains unaffected. We identified two distinct mechanisms, conserved across multiple Ascomycota species, by which this molecular adaptation occurred. Our results suggest a route for TF evolution that alleviates negative pleiotropic effects by modularly gaining new binding sites. These findings expand our current understanding of ZF DNA binding and provide evidence for paralogous ZFs utilizing alternate modes of DNA binding to recognize unique sets of noncanonical binding sites.

NextPBM: a platform to study cell-specific transcription factor binding and cooperativity.

High-throughput (HT) in vitro methods for measuring protein-DNA binding have become invaluable for characterizing transcription factor (TF) complexes and modeling gene regulation. However, current methods do not utilize endogenous proteins and, therefore, do not quantify the impact of cell-specific post-translational modifications (PTMs) and cooperative cofactors. We introduce the HT nextPBM (nuclear extract protein-binding microarray) approach to study DNA binding of native cellular TFs that accounts for PTMs and cell-specific cofactors. We integrate immune-depletion and phosphatase treatment steps into our nextPBM pipeline to characterize the impact of cofactors and phosphorylation on TF binding. We analyze binding of PU.1/SPI1 and IRF8 from human monocytes, delineate DNA-sequence determinants for their cooperativity, and show how PU.1 affinity correlates with enhancer status and the presence of cooperative and collaborative cofactors. We describe how nextPBMs, and our accompanying computational framework, can be used to discover cell-specific cofactors, screen for synthetic cooperative DNA elements, and characterize TF cooperativity.

The CLAMP protein links the MSL complex to the X chromosome during Drosophila dosage compensation.

The Drosophila male-specific lethal (MSL) dosage compensation complex increases transcript levels on the single male X chromosome to equal the transcript levels in XX females. However, it is not known how the MSL complex is linked to its DNA recognition elements, the critical first step in dosage compensation. Here, we demonstrate that a previously uncharacterized zinc finger protein, CLAMP (chromatin-linked adaptor for MSL proteins), functions as the first link between the MSL complex and the X chromosome. CLAMP directly binds to the MSL complex DNA recognition elements and is required for the recruitment of the MSL complex. The discovery of CLAMP identifies a key factor required for the chromosome-specific targeting of dosage compensation, providing new insights into how subnuclear domains of coordinate gene regulation are formed within metazoan genomes.

DNA-binding landscape of IRF3, IRF5 and IRF7 dimers: implications for dimer-specific gene regulation.

Transcription factors IRF3, IRF5 and IRF7 (IRF3/5/7) have overlapping, yet distinct, roles in the mammalian response to pathogens. To examine the role that DNA-binding specificity plays in delineating IRF3/5/7-specific gene regulation we used protein-binding microarrays (PBMs) to characterize the DNA binding of IRF3/5/7 homodimers. We identified both common and dimer-specific DNA binding sites, and show that DNA-binding differences can translate into dimer-specific gene regulation. Central to the antiviral response, IRF3/5/7 regulate type I interferon (IFN) genes. We show that IRF3 and IRF7 bind to many interferon-stimulated response element (ISRE)-type sites in the virus-response elements (VREs) of IFN promoters. However, strikingly, IRF5 does not bind the VREs, suggesting evolutionary selection against IRF5 homodimer binding. Mutational analysis reveals a critical specificity-determining residue that inhibits IRF5 binding to the ISRE-variants present in the IFN gene promoters. Integrating PBM and reporter gene data we find that both DNA-binding affinity and affinity-independent mechanisms determine the function of DNA-bound IRF dimers, suggesting that DNA-based allostery plays a role in IRF binding site function. Our results provide new insights into the role and limitations of DNA-binding affinity in delineating IRF3/5/7-specific gene expression.

Expansion of GA Dinucleotide Repeats Increases the Density of CLAMP Binding Sites on the X-Chromosome to Promote Drosophila Dosage Compensation.

Dosage compensation is an essential process that equalizes transcript levels of X-linked genes between sexes by forming a domain of coordinated gene expression. Throughout the evolution of Diptera, many different X-chromosomes acquired the ability to be dosage compensated. Once each newly evolved X-chromosome is targeted for dosage compensation in XY males, its active genes are upregulated two-fold to equalize gene expression with XX females. In Drosophila melanogaster, the CLAMP zinc finger protein links the dosage compensation complex to the X-chromosome. However, the mechanism for X-chromosome identification has remained unknown. Here, we combine biochemical, genomic and evolutionary approaches to reveal that expansion of GA-dinucleotide repeats likely accumulated on the X-chromosome over evolutionary time to increase the density of CLAMP binding sites, thereby driving the evolution of dosage compensation. Overall, we present new insight into how subtle changes in genomic architecture, such as expansions of a simple sequence repeat, promote the evolution of coordinated gene expression.

Noncooperative interactions between transcription factors and clustered DNA binding sites enable graded transcriptional responses to environmental inputs.

A paradigm in transcriptional regulation is that graded increases in transcription factor (TF) concentration are translated into on/off transcriptional responses by cooperative TF binding to adjacent sites. Digital transcriptional responses underlie the definition of anatomical boundaries during development. Here we show that NF-kappaB, a TF controlling inflammation and immunity, is conversely an analog transcriptional regulator that uses clustered binding sites noncooperatively. We observed that increasing concentrations of NF-kappaB are translated into gradual increments in gene transcription. We provide a thermodynamic interpretation of the experimental observations by combining quantitative measurements and a minimal physical model of an NF-kappaB-dependent promoter. We demonstrate that NF-kappaB binds independently to adjacent sites to promote additive RNA Pol II recruitment and graded transcriptional outputs. These findings reveal an alternative mode of operation of clustered TF binding sites, which might function in biological conditions where the transcriptional output is proportional to the strength of an environmental input.

Precise temporal control of the eye regulatory gene Pax6 via enhancer-binding site affinity.

How transcription factors interpret the cis-regulatory logic encoded within enhancers to mediate quantitative changes in spatiotemporally restricted expression patterns during animal development is not well understood. Pax6 is a dosage-sensitive gene essential for eye development. Here, we identify the Prep1 (pKnox1) transcription factor as a critical dose-dependent upstream regulator of Pax6 expression during lens formation. We show that Prep1 activates the Pax6 lens enhancer by binding to two phylogenetically conserved lower-affinity DNA-binding sites. Finally, we describe a mechanism whereby Pax6 levels are determined by transcriptional synergy of Prep1 bound to the two sites, while timing of enhancer activation is determined by binding site affinity.

Protein-DNA binding: complexities and multi-protein codes.

Binding of proteins to particular DNA sites across the genome is a primary determinant of specificity in genome maintenance and gene regulation. DNA-binding specificity is encoded at multiple levels, from the detailed biophysical interactions between proteins and DNA, to the assembly of multi-protein complexes. At each level, variation in the mechanisms used to achieve specificity has led to difficulties in constructing and applying simple models of DNA binding. We review the complexities in protein-DNA binding found at multiple levels and discuss how they confound the idea of simple recognition codes. We discuss the impact of new high-throughput technologies for the characterization of protein-DNA binding, and how these technologies are uncovering new complexities in protein-DNA recognition. Finally, we review the concept of multi-protein recognition codes in which new DNA-binding specificities are achieved by the assembly of multi-protein complexes.

Grhl2 is required in nonneural tissues for neural progenitor survival and forebrain development.

Grainyhead-like genes are part of a highly conserved gene family that play a number of roles in ectoderm development and maintenance in mammals. Here we identify a novel allele of Grhl2, cleft-face 3 (clft3), in a mouse line recovered from an ENU mutagenesis screen for organogenesis defects. Homozygous clft3 mutants have a number of phenotypes in common with other alleles of Grhl2. We note a significant effect of genetic background on the clft3 phenotype. One of these is a reduction in size of the telencephalon where we find abnormal patterns of neural progenitor mitosis and apoptosis in mutant brains. Interestingly, Grhl2 is not expressed in the developing forebrain, suggesting this is a survival factor for neural progenitors exerting a paracrine effect on the neural tissue from the overlying ectoderm where Grhl2 is highly expressed. genesis 53:573-582, 2015. © 2015 Wiley Periodicals, Inc.

Targeting of SUMOylation leads to cBAF complex stabilization and disruption of the SS18::SSX transcriptome in Synovial Sarcoma.

Synovial Sarcoma (SS) is driven by the SS18::SSX fusion oncoprotein and is ultimately refractory to therapeutic approaches. SS18::SSX alters ATP-dependent chromatin remodeling BAF (mammalian SWI/SNF) complexes, leading to the degradation of canonical (cBAF) complex and amplified presence of an SS18::SSX-containing non-canonical BAF (ncBAF or GBAF) that drives an SS-specific transcription program and tumorigenesis. We demonstrate that SS18::SSX activates the SUMOylation program and SSs are sensitive to the small molecule SAE1/2 inhibitor, TAK-981. Mechanistically, TAK-981 de-SUMOylates the cBAF subunit SMARCE1, stabilizing and restoring cBAF on chromatin, shifting away from SS18::SSX-ncBAF-driven transcription, associated with DNA damage and cell death and resulting in tumor inhibition across both human and mouse SS tumor models. TAK-981 synergized with cytotoxic chemotherapy through increased DNA damage, leading to tumor regression. Targeting the SUMOylation pathway in SS restores cBAF complexes and blocks the SS18::SSX-ncBAF transcriptome, identifying a therapeutic vulnerability in SS, positioning the in-clinic TAK-981 to treat SS.

Targeting of SUMOylation leads to cBAF complex stabilization and disruption of the SS18::SSX transcriptome in Synovial Sarcoma.

Synovial Sarcoma (SS) is driven by the SS18::SSX fusion oncoprotein. and is ultimately refractory to therapeutic approaches. SS18::SSX alters ATP-dependent chromatin remodeling BAF (mammalian SWI/SNF) complexes, leading to the degradation of canonical (cBAF) complex and amplified presence of an SS18::SSX-containing non-canonical BAF (ncBAF or GBAF) that drives an SS-specific transcription program and tumorigenesis. We demonstrate that SS18::SSX activates the SUMOylation program and SSs are sensitive to the small molecule SAE1/2 inhibitor, TAK-981. Mechanistically, TAK-981 de-SUMOylates the cBAF subunit SMARCE1, stabilizing and restoring cBAF on chromatin, shifting away from SS18::SSX-ncBAF-driven transcription, associated with DNA damage and cell death and resulting in tumor inhibition across both human and mouse SS tumor models. TAK-981 synergized with cytotoxic chemotherapy through increased DNA damage, leading to tumor regression. Targeting the SUMOylation pathway in SS restores cBAF complexes and blocks the SS18::SSX-ncBAF transcriptome, identifying a therapeutic vulnerability in SS, positioning the in-clinic TAK-981 to treat SS.

Widespread perturbation of ETS factor binding sites in cancer.

Although >90% of somatic mutations reside in non-coding regions, few have been reported as cancer drivers. To predict driver non-coding variants (NCVs), we present a transcription factor (TF)-aware burden test based on a model of coherent TF function in promoters. We apply this test to NCVs from the Pan-Cancer Analysis of Whole Genomes cohort and predict 2555 driver NCVs in the promoters of 813 genes across 20 cancer types. These genes are enriched in cancer-related gene ontologies, essential genes, and genes associated with cancer prognosis. We find that 765 candidate driver NCVs alter transcriptional activity, 510 lead to differential binding of TF-cofactor regulatory complexes, and that they primarily impact the binding of ETS factors. Finally, we show that different NCVs within a promoter often affect transcriptional activity through shared mechanisms. Our integrated computational and experimental approach shows that cancer NCVs are widespread and that ETS factors are commonly disrupted.

Starvation decreases immunity and immune regulatory factor NF-κB in the starlet sea anemone Nematostella vectensis.

Lack of proper nutrition has important consequences for the physiology of all organisms, and nutritional status can affect immunity, based on many studies in terrestrial animals. Here we show a positive correlation between nutrition and immunity in the sea anemone Nematostella vectensis. Gene expression profiling of adult anemones shows downregulation of genes involved in nutrient metabolism, cellular respiration, and immunity in starved animals. Starved adult anemones also have reduced protein levels and activity of immunity transcription factor NF-κB. Starved juvenile anemones have increased sensitivity to bacterial infection and also have lower NF-κB protein levels, as compared to fed controls. Weighted Gene Correlation Network Analysis (WGCNA) is used to identify significantly correlated gene networks that were downregulated with starvation. These experiments demonstrate a correlation between nutrition and immunity in an early diverged marine metazoan, and the results have implications for the survival of marine organisms as they encounter changing environments.

Parallelized multidimensional analytic framework applied to mammary epithelial cells uncovers regulatory principles in EMT.

A proper understanding of disease etiology will require longitudinal systems-scale reconstruction of the multitiered architecture of eukaryotic signaling. Here we combine state-of-the-art data acquisition platforms and bioinformatics tools to devise PAMAF, a workflow that simultaneously examines twelve omics modalities, i.e., protein abundance from whole-cells, nucleus, exosomes, secretome and membrane; N-glycosylation, phosphorylation; metabolites; mRNA, miRNA; and, in parallel, single-cell transcriptomes. We apply PAMAF in an established in vitro model of TGFß-induced epithelial to mesenchymal transition (EMT) to quantify >61,000 molecules from 12 omics and 10 timepoints over 12 days. Bioinformatics analysis of this EMT-ExMap resource allowed us to identify; -topological coupling between omics, -four distinct cell states during EMT, -omics-specific kinetic paths, -stage-specific multi-omics characteristics, -distinct regulatory classes of genes, -ligand-receptor mediated intercellular crosstalk by integrating scRNAseq and subcellular proteomics, and -combinatorial drug targets (e.g., Hedgehog signaling and CAMK-II) to inhibit EMT, which we validate using a 3D mammary duct-on-a-chip platform. Overall, this study provides a resource on TGFß signaling and EMT.

Non-DNA-binding cofactors enhance DNA-binding specificity of a transcriptional regulatory complex.

Recruitment of cofactors to specific DNA sites is integral for specificity in gene regulation. As a model system, we examined how targeting and transcriptional control of the sulfur metabolism genes in Saccharomyces cerevisiae is governed by recruitment of the transcriptional co-activator Met4. We developed genome-scale approaches to measure transcription factor (TF) DNA-binding affinities and cofactor recruitment to >1300 genomic binding site sequences. We report that genes responding to the TF Cbf1 and cofactor Met28 contain a novel 'recruitment motif' (RYAAT), adjacent to Cbf1 binding sites, which enhances the binding of a Met4-Met28-Cbf1 regulatory complex, and that abrogation of this motif significantly reduces gene induction under low-sulfur conditions. Furthermore, we show that correct recognition of this composite motif requires both non-DNA-binding cofactors Met4 and Met28. Finally, we demonstrate that the presence of an RYAAT motif next to a Cbf1 site, rather than Cbf1 binding affinity, specifies Cbf1-dependent sulfur metabolism genes. Our results highlight the need to examine TF/cofactor complexes, as novel specificity can result from cofactors that lack intrinsic DNA-binding specificity.

CASCADE: high-throughput characterization of regulatory complex binding altered by non-coding variants.

Non-coding DNA variants (NCVs) impact gene expression by altering binding sites for regulatory complexes. New high-throughput methods are needed to characterize the impact of NCVs on regulatory complexes. We developed CASCADE (Customizable Approach to Survey Complex Assembly at DNA Elements), an array-based high-throughput method to profile cofactor (COF) recruitment. CASCADE identifies DNA-bound transcription factor-cofactor (TF-COF) complexes in nuclear extracts and quantifies the impact of NCVs on their binding. We demonstrate CASCADE sensitivity in characterizing condition-specific recruitment of COFs p300 and RBBP5 (MLL subunit) to the CXCL10 promoter in lipopolysaccharide (LPS)-stimulated human macrophages and quantify the impact of all possible NCVs. To demonstrate applicability to NCV screens, we profile TF-COF binding to ~1,700 single-nucleotide polymorphism quantitative trait loci (SNP-QTLs) in human macrophages and identify perturbed ETS domain-containing complexes. CASCADE will facilitate high-throughput testing of molecular mechanisms of NCVs for diverse biological applications.

High-Throughput Analysis of the Cell and DNA Site-Specific Binding of Native NF-κB Dimers Using Nuclear Extract Protein-Binding Microarrays (NextPBMs).

Nuclear factor-kappa B (NF-κB) transcription factors coordinate gene expression in response to a broad array of cellular signals. In vertebrates, there are five NF-κB proteins (c-Rel, RelA/p65, RelB, p50, and p52) that can form various dimeric combinations exhibiting both common and dimer-specific DNA-binding specificity. In this chapter, we describe the use of the nuclear extract protein-binding microarray (nextPBM), a high-throughput method to characterize the DNA binding of transcription factors present in cell nuclear extracts. NextPBMs allow for sensitive analysis of the DNA binding of NF-κB dimers and their interactions with cell-specific cofactors.

Comparison of NF-κB from the protists Capsaspora owczarzaki and Acanthoeca spectabilis reveals extensive evolutionary diversification of this transcription factor.

We provide a functional characterization of transcription factor NF-κB in protists and provide information about the evolution and diversification of this biologically important protein. We characterized NF-κB in two protists using phylogenetic, cellular, and biochemical techniques. NF-κB of the holozoan Capsaspora owczarzaki (Co) has an N-terminal DNA-binding domain and a C-terminal Ankyrin repeat (ANK) domain, and its DNA-binding specificity is more similar to metazoan NF-κB proteins than to Rel proteins. Removal of the ANK domain allows Co-NF-κB to enter the nucleus, bind DNA, and activate transcription. However, C-terminal processing of Co-NF-κB is not induced by IκB kinases in human cells. Overexpressed Co-NF-κB localizes to the cytoplasm in Co cells. Co-NF-κB mRNA and DNA-binding levels differ across three Capsaspora life stages. RNA-sequencing and GO analyses identify possible gene targets of Co-NF-κB. Three NF-κB-like proteins from the choanoflagellate Acanthoeca spectabilis (As) contain conserved Rel Homology domain sequences, but lack C-terminal ANK repeats. All three As-NF-κB proteins constitutively enter the nucleus of cells, but differ in their DNA-binding abilities, transcriptional activation activities, and dimerization properties. These results provide a basis for understanding the evolutionary origins of this key transcription factor and could have implications for the origins of regulated immunity in higher taxa.