Zebrafish Expression Ontology of Gene Sets (ZEOGS): a tool to analyze enrichment of zebrafish anatomical terms in large gene sets

The zebrafish (Danio rerio) is an established model organism for developmental and biomedical research. It is frequently used for high-throughput functional genomics experiments, such as genome-wide gene expression measurements, to systematically analyze molecular mechanisms. However, the use of whole embryos or larvae in such experiments leads to a loss of the spatial information. To address this problem, we have developed a tool called Zebrafish Expression Ontology of Gene Sets (ZEOGS) to assess the enrichment of anatomical terms in large gene sets. ZEOGS uses gene expression pattern data from several sources: first, in situ hybridization experiments from the Zebrafish Model Organism Database (ZFIN); second, it uses the Zebrafish Anatomical Ontology, a controlled vocabulary that describes connected anatomical structures; and third, the available connections between expression patterns and anatomical terms contained in ZFIN. Upon input of a gene set, ZEOGS determines which anatomical structures are overrepresented in the input gene set. ZEOGS allows one for the first time to look at groups of genes and to describe them in terms of shared anatomical structures. To establish ZEOGS, we first tested it on random gene selections and on two public microarray datasets with known tissue-specific gene expression changes. These tests showed that ZEOGS could reliably identify the tissues affected, whereas only very few enriched terms to none were found in the random gene sets. Next we applied ZEOGS to microarray datasets of 24 and 72 h postfertilization zebrafish embryos treated with beclomethasone, a potent glucocorticoid. This analysis resulted in the identification of several anatomical terms related to glucocorticoid-responsive tissues, some of which were stage-specific. Our studies highlight the ability of ZEOGS to extract spatial information from datasets derived from whole embryos, indicating that ZEOGS could be a useful tool to automatically analyze gene expression pattern features of any large zebrafish gene set

Identification and characterization of BATF3 as a context-specific coactivator of the glucocorticoid receptor

The ability of the glucocorticoid receptor (GR) to regulate the transcriptional output of genes relies on its interactions with transcriptional coregulators. However, which coregulators are required for GR-dependent activation is context-dependent and can be influenced by the sequence of the DNA bound by GR and by the nature of the GR isoform responsible for the regulation of a gene. Here, we screened for GR-interacting proteins for which the interaction signal differed between two GR isoforms GRalpha and GRgamma. These isoforms diverge by a single amino acid insertion in a domain, the lever arm, which adopts DNA sequence-specific conformations. We identify Basic Leucine Zipper ATF-Like Transcription Factor 3 (BATF3), an AP-1 family transcription factor, as a GR coregulator whose interaction with GR is modulated by the lever arm. Further, a combination of experiments uncovered that BATF3 acts as a gene-specific coactivator of GR whose coactivator potency is influenced by the sequence of the GR binding site. Together, our findings suggest that GR isoform and the sequence of GR binding site influence the interaction of GR with BATF3, which might direct the assembly of gene-specific regulatory complexes to fine-tune the expression of individual GR target genes

A naturally occuring insertion of a single amino acid rewires transcriptional regulation by glucocorticoid receptor isoforms

In addition to guiding proteins to defined genomic loci, DNA can act as an allosteric ligand that influences protein structure and activity. Here we compared genome-wide binding, transcriptional regulation, and, using NMR, the conformation of two glucocorticoid receptor (GR) isoforms that differ by a single amino acid insertion in the lever arm, a domain that adopts DNA sequence-specific conformations. We show that these isoforms differentially regulate gene expression levels through two mechanisms: differential DNA binding and altered communication between GR domains. Our studies suggest a versatile role for DNA in both modulating GR activity and also in directing the use of GR isoforms. We propose that the lever arm is a "fulcrum" for bidirectional allosteric signaling, conferring conformational changes in the DNA reading head that influence DNA sequence selectivity, as well as conferring changes in the dimerization domain that connect functionally with remote regulatory surfaces, thereby influencing which genes are regulated and the magnitude of their regulation

Expanding the repertoire of glucocorticoid receptor target genes by engineering genomic response elements

The glucocorticoid receptor (GR), a hormone-activated transcription factor, binds to a myriad of genomic binding sites yet seems to regulate a much smaller number of genes. Genome-wide analysis of GR binding and gene regulation has shown that the likelihood of GR-dependent regulation increases with decreased distance of its binding to the transcriptional start site of a gene. To test if we can adopt this knowledge to expand the repertoire of GR target genes, we used CRISPR/Cas-mediated homology-directed repair to add a single GR-binding site directly upstream of the transcriptional start site of each of four genes. To our surprise, we found that the addition of a single GR-binding site can be enough to convert a gene into a GR target. The gain of GR-dependent regulation was observed for two of four genes analyzed and coincided with acquired GR binding at the introduced binding site. However, the gene-specific gain of GR-dependent regulation could not be explained by obvious differences in chromatin accessibility between converted genes and their non-converted counterparts. Furthermore, by introducing GR-binding sequences with different nucleotide compositions, we show that activation can be facilitated by distinct sequences without obvious differences in activity between the GR-binding sequence variants we tested. The approach to use genome engineering to build genomic response elements facilitates the generation of cell lines with tailored repertoires of GR-responsive genes and a framework to test and refine our understanding of the cis-regulatory logic of gene regulation by testing if engineered response elements behave as predicted

Pluripotency reprogramming by competent and incompetent POU factors uncovers temporal dependency for Oct4 and Sox2

Oct4, along with Sox2 and Klf4 (SK), can induce pluripotency but structurally similar factors like Oct6 cannot. To decode why Oct4 has this unique ability, we compare Oct4-binding, accessibility patterns and transcriptional waves with Oct6 and an Oct4 mutant defective in the dimerization with Sox2 (Oct4(defSox2)). We find that initial silencing of the somatic program proceeds indistinguishably with or without Oct4. Oct6 mitigates the mesenchymal-to-epithelial transition and derails reprogramming. These effects are a consequence of differences in genome-wide binding, as the early binding profile of Oct4(defSox2) resembles Oct4, whilst Oct6 does not bind pluripotency enhancers. Nevertheless, in the Oct6-SK condition many otherwise Oct4-bound locations become accessible but chromatin opening is compromised when Oct4(defSox2) occupies these sites. We find that Sox2 predominantly facilitates chromatin opening, whilst Oct4 serves an accessory role. Formation of Oct4/Sox2 heterodimers is essential for pluripotency establishment; however, reliance on Oct4/Sox2 heterodimers declines during pluripotency maintenance

Dynamic regulation of glucocorticoid signalling in health and disease

Activation of the glucocorticoid receptor (GR) by endogenous and synthetic glucocorticoids regulates hundreds of genes to control regulatory networks in development, metabolism, cognition and inflammation. Elucidation of the mechanisms that regulate glucocorticoid action has highlighted the dynamic nature of hormone signalling and provides novel insights into genomic glucocorticoid actions. The major factors that regulate GR function include chromatin structure, epigenetics, genetic variation and the pattern of glucocorticoid hormone secretion. We review our current understanding of the mechanisms that contribute to GR signalling and how these contribute to glucocorticoid sensitivity, resistance and side effects

Insights into Hox Protein Function from a Large Scale Combinatorial Analysis of Protein Domains

Protein function is encoded within protein sequence and protein domains. However, how protein domains cooperate within a protein to modulate overall activity and how this impacts functional diversification at the molecular and organism levels remains largely unaddressed. Focusing on three domains of the central class Drosophila Hox transcription factor AbdominalA (AbdA), we used combinatorial domain mutations and most known AbdA developmental functions as biological readouts to investigate how protein domains collectively shape protein activity. The results uncover redundancy, interactivity, and multifunctionality of protein domains as salient features underlying overall AbdA protein activity, providing means to apprehend functional diversity and accounting for the robustness of Hox-controlled developmental programs. Importantly, the results highlight context-dependency in protein domain usage and interaction, allowing major modifications in domains to be tolerated without general functional loss. The non-pleoitropic effect of domain mutation suggests that protein modification may contribute more broadly to molecular changes underlying morphological diversification during evolution, so far thought to rely largely on modification in gene cis-regulatory sequences

The Elongator Complex Interacts with PCNA and Modulates Transcriptional Silencing and Sensitivity to DNA Damage Agents

Histone chaperones CAF-1 and Asf1 function to deposit newly synthesized histones onto replicating DNA to promote nucleosome formation in a proliferating cell nuclear antigen (PCNA) dependent process. The DNA replication- or DNA repair-coupled nucleosome assembly pathways are important for maintenance of transcriptional gene silencing and genome stability. However, how these pathways are regulated is not well understood. Here we report an interaction between the Elongator histone acetyltransferase and the proliferating cell nuclear antigen. Cells lacking Elp3 (K-acetyltransferase Kat9), the catalytic subunit of the six-subunit Elongator complex, partially lose silencing of reporter genes at the chromosome VIIL telomere and at the HMR locus, and are sensitive to the DNA replication inhibitor hydroxyurea (HU) and the damaging agent methyl methanesulfonate (MMS). Like deletion of the ELP3, mutation of each of the four other subunits of the Elongator complex as well as mutations in Elp3 that compromise the formation of the Elongator complex also result in loss of silencing and increased HU sensitivity. Moreover, Elp3 is required for S-phase progression in the presence of HU. Epistasis analysis indicates that the elp3Δ mutant, which itself is sensitive to MMS, exacerbates the MMS sensitivity of cells lacking histone chaperones Asf1, CAF-1 and the H3 lysine 56 acetyltransferase Rtt109. The elp3Δ mutant has allele specific genetic interactions with mutations in POL30 that encodes PCNA and PCNA binds to the Elongator complex both in vivo and in vitro. Together, these results uncover a novel role for the intact Elongator complex in transcriptional silencing and maintenance of genome stability, and it does so in a pathway linked to the DNA replication and DNA repair protein PCNA

Low-Resolution Molecular Models Reveal the Oligomeric State of the PPAR and the Conformational Organization of Its Domains in Solution

The peroxisome proliferator-activated receptors (PPARs) regulate genes involved in lipid and carbohydrate metabolism, and are targets of drugs approved for human use. Whereas the crystallographic structure of the complex of full length PPARγ and RXRα is known, structural alterations induced by heterodimer formation and DNA contacts are not well understood. Herein, we report a small-angle X-ray scattering analysis of the oligomeric state of hPPARγ alone and in the presence of retinoid X receptor (RXR). The results reveal that, in contrast with other studied nuclear receptors, which predominantly form dimers in solution, hPPARγ remains in the monomeric form by itself but forms heterodimers with hRXRα. The low-resolution models of hPPARγ/RXRα complexes predict significant changes in opening angle between heterodimerization partners (LBD) and extended and asymmetric shape of the dimer (LBD-DBD) as compared with X-ray structure of the full-length receptor bound to DNA. These differences between our SAXS models and the high-resolution crystallographic structure might suggest that there are different conformations of functional heterodimer complex in solution. Accordingly, hydrogen/deuterium exchange experiments reveal that the heterodimer binding to DNA promotes more compact and less solvent-accessible conformation of the receptor complex

Widespread Site-Dependent Buffering of Human Regulatory Polymorphism

The average individual is expected to harbor thousands of variants within non-coding genomic regions involved in gene regulation. However, it is currently not possible to interpret reliably the functional consequences of genetic variation within any given transcription factor recognition sequence. To address this, we comprehensively analyzed heritable genome-wide binding patterns of a major sequence-specific regulator (CTCF) in relation to genetic variability in binding site sequences across a multi-generational pedigree. We localized and quantified CTCF occupancy by ChIP-seq in 12 related and unrelated individuals spanning three generations, followed by comprehensive targeted resequencing of the entire CTCF–binding landscape across all individuals. We identified hundreds of variants with reproducible quantitative effects on CTCF occupancy (both positive and negative). While these effects paralleled protein–DNA recognition energetics when averaged, they were extensively buffered by striking local context dependencies. In the significant majority of cases buffering was complete, resulting in silent variants spanning every position within the DNA recognition interface irrespective of level of binding energy or evolutionary constraint. The prevalence of complex partial or complete buffering effects severely constrained the ability to predict reliably the impact of variation within any given binding site instance. Surprisingly, 40% of variants that increased CTCF occupancy occurred at positions of human–chimp divergence, challenging the expectation that the vast majority of functional regulatory variants should be deleterious. Our results suggest that, even in the presence of “perfect” genetic information afforded by resequencing and parallel studies in multiple related individuals, genomic site-specific prediction of the consequences of individual variation in regulatory DNA will require systematic coupling with empirical functional genomic measurements