Role of the chromatin landscape and sequence in determining cell type-specific genomic glucocorticoid receptor binding and gene regulation.

The genomic loci bound by the glucocorticoid receptor (GR), a hormone-activated transcription factor, show little overlap between cell types. To study the role of chromatin and sequence in specifying where GR binds, we used Bayesian modeling within the universe of accessible chromatin. Taken together, our results uncovered that although GR preferentially binds accessible chromatin, its binding is biased against accessible chromatin located at promoter regions. This bias can only be explained partially by the presence of fewer GR recognition sequences, arguing for the existence of additional mechanisms that interfere with GR binding at promoters. Therefore, we tested the role of H3K9ac, the chromatin feature with the strongest negative association with GR binding, but found that this correlation does not reflect a causative link. Finally, we find a higher percentage of promoter-proximal GR binding for genes regulated by GR across cell types than for cell type-specific target genes. Given that GR almost exclusively binds accessible chromatin, we propose that cell type-specific regulation by GR preferentially occurs via distal enhancers, whose chromatin accessibility is typically cell type-specific, whereas ubiquitous target gene regulation is more likely to result from binding to promoter regions, which are often accessible regardless of cell type examined.

Incoherent feed-forward regulatory logic underpinning glucocorticoid receptor action.

The complexity and specificity of metazoan transcription are determined by combinatorial control of the composition and activity of regulatory complexes. To investigate the basis of this specificity, we focused on the glucocorticoid receptor (GR), a single regulatory factor that integrates multiple signals to give rise to many distinct patterns of expression. We measured the expression of a set of genes, each directly GR-regulated, but by different mechanisms in two cell lines. We varied ligand (dose, chemistry, and duration of treatment), GR (expression level and functionality), and a non-GR regulatory factor that commonly interacts with GR. Our study revealed distinct expression patterns within this set of genes, but all could be modeled by an incoherent feed-forward regulatory logic. Cellular signals, operating on GR and other factors within regulatory complexes, may define and modulate the kinetics and strength of the activating or inhibitory paths of the regulatory logic. Thus, characterizing systems behavior by perturbing single or multiple signals can reveal general principles of regulation, providing an approach to the dissection and deconvolution of combinatorial control.

A naturally occurring 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.

The Mediator subunit MDT-15 confers metabolic adaptation to ingested material.

In eukaryotes, RNA polymerase II (Pol(II)) dependent gene expression requires accessory factors termed transcriptional coregulators. One coregulator that universally contributes to Pol(II)-dependent transcription is the Mediator, a multisubunit complex that is targeted by many transcriptional regulatory factors. For example, the Caenorhabditis elegans Mediator subunit MDT-15 confers the regulatory actions of the sterol response element binding protein SBP-1 and the nuclear hormone receptor NHR-49 on fatty acid metabolism. Here, we demonstrate that MDT-15 displays a broader spectrum of activities, and that it integrates metabolic responses to materials ingested by C. elegans. Depletion of MDT-15 protein or mutation of the mdt-15 gene abrogated induction of specific detoxification genes in response to certain xenobiotics or heavy metals, rendering these animals hypersensitive to toxin exposure. Intriguingly, MDT-15 appeared to selectively affect stress responses related to ingestion, as MDT-15 functional defects did not abrogate other stress responses, e.g., thermotolerance. Together with our previous finding that MDT-15:NHR-49 regulatory complexes coordinate a sector of the fasting response, we propose a model whereby MDT-15 integrates several transcriptional regulatory pathways to monitor both the availability and quality of ingested materials, including nutrients and xenobiotic compounds.

Conservation analysis predicts in vivo occupancy of glucocorticoid receptor-binding sequences at glucocorticoid-induced genes.

The glucocorticoid receptor (GR) interacts with specific GR-binding sequences (GBSs) at glucocorticoid response elements (GREs) to orchestrate transcriptional networks. Although the sequences of the GBSs are highly variable among different GREs, the precise sequence within an individual GRE is highly conserved. In this study, we examined whether sequence conservation of sites resembling GBSs is sufficient to predict GR occupancy of GREs at genes responsive to glucocorticoids. Indeed, we found that the level of conservation of these sites at genes up-regulated by glucocorticoids in mouse C3H10T1/2 mesenchymal stem-like cells correlated directly with the extent of occupancy by GR. In striking contrast, we failed to observe GR occupancy of GBSs at genes repressed by glucocorticoids, despite the occurrence of these sites at a frequency similar to that of the induced genes. Thus, GR occupancy of the GBS motif correlates with induction but not repression, and GBS conservation alone is sufficient to predict GR occupancy and GRE function at induced genes.

Systematic quantification of HDR and NHEJ reveals effects of locus, nuclease, and cell type on genome-editing.

Precise genome-editing relies on the repair of sequence-specific nuclease-induced DNA nicking or double-strand breaks (DSBs) by homology-directed repair (HDR). However, nonhomologous end-joining (NHEJ), an error-prone repair, acts concurrently, reducing the rate of high-fidelity edits. The identification of genome-editing conditions that favor HDR over NHEJ has been hindered by the lack of a simple method to measure HDR and NHEJ directly and simultaneously at endogenous loci. To overcome this challenge, we developed a novel, rapid, digital PCR-based assay that can simultaneously detect one HDR or NHEJ event out of 1,000 copies of the genome. Using this assay, we systematically monitored genome-editing outcomes of CRISPR-associated protein 9 (Cas9), Cas9 nickases, catalytically dead Cas9 fused to FokI, and transcription activator-like effector nuclease at three disease-associated endogenous gene loci in HEK293T cells, HeLa cells, and human induced pluripotent stem cells. Although it is widely thought that NHEJ generally occurs more often than HDR, we found that more HDR than NHEJ was induced under multiple conditions. Surprisingly, the HDR/NHEJ ratios were highly dependent on gene locus, nuclease platform, and cell type. The new assay system, and our findings based on it, will enable mechanistic studies of genome-editing and help improve genome-editing technology.

Opposing chemokine gradients control human thymocyte migration in situ.

The ordered migration of thymocytes from the cortex to the medulla is critical for the appropriate selection of the mature T cell repertoire. Most studies of thymocyte migration rely on mouse models, but we know relatively little about how human thymocytes find their appropriate anatomical niches within the thymus. Moreover, the signals that retain CD4+CD8+ double-positive (DP) thymocytes in the cortex and prevent them from entering the medulla prior to positive selection have not been identified in mice or humans. Here, we examined the intrathymic migration of human thymocytes in both mouse and human thymic stroma and found that human thymocyte subsets localized appropriately to the cortex on mouse thymic stroma and that MHC-dependent interactions between human thymocytes and mouse stroma could maintain the activation and motility of DP cells. We also showed that CXCR4 was required to retain human DP thymocytes in the cortex, whereas CCR7 promoted migration of mature human thymocytes to the medulla. Thus, 2 opposing chemokine gradients control the migration of thymocytes from the cortex to the medulla. These findings point to significant interspecies conservation in thymocyte-stroma interactions and provide the first evidence that chemokines not only attract mature thymocytes to the medulla, but also play an active role in retaining DP thymocytes in the cortex prior to positive selection.

Codon usage determines the mutational robustness, evolutionary capacity, and virulence of an RNA virus.

RNA viruses exist as dynamic and diverse populations shaped by constant mutation and selection. Yet little is known about how the mutant spectrum contributes to virus evolvability and pathogenesis. Because several codon choices are available for a given amino acid, a central question concerns whether viral sequences have evolved to optimize not only the protein coding consensus, but also the DNA/RNA sequences accessible through mutation. Here we directly test this hypothesis by comparing wild-type poliovirus to synthetic viruses carrying re-engineered capsid sequences with hundreds of synonymous mutations. Strikingly, such rewiring of the population's mutant network reduced its robustness and attenuated the virus in an animal model of infection. We conclude that the position of a virus in sequence space defines its mutant spectrum, evolutionary trajectory, and pathogenicity. This organizing principle for RNA virus populations confers tolerance to mutations and facilitates replication and spread within the dynamic host environment.

Occupancy of chromatin organizers in the Epstein-Barr virus genome.

The human CCCTC-binding factor, CTCF, regulates transcription of the double-stranded DNA genomes of herpesviruses. The architectural complex cohesin and RNA Polymerase II also contribute to this organization. We profiled the occupancy of CTCF, cohesin, and RNA Polymerase II on the episomal genome of the Epstein-Barr virus in a cell culture model of latent infection. CTCF colocalizes with cohesin but not RNA Polymerase II. CTCF and cohesin bind specific sequences throughout the genome that are found not just proximal to the regulatory elements of latent genes, but also near lytic genes. In addition to tracking with known transcripts, RNA Polymerase II appears at two unannotated positions, one of which lies within the latent origin of replication. The widespread occupancy profile of each protein reveals binding near or at a myriad of regulatory elements and suggests context-dependent functions.