Genome-wide Analysis Identifies Nuclear Factor 1C as a Novel Transcription Factor and Potential Therapeutic Target in Small Cell Lung Cancer

BACKGROUND

Recent insights regarding mechanisms mediating stemness, heterogeneity, and metastatic potential of lung cancers have yet to be fully translated to effective regimens for the treatment of these malignancies. This study sought to identify novel targets for lung cancer therapy.

METHODS

Transcriptomes and DNA methylomes of 14 SCLC and 10 NSCLC lines were compared to normal human small airway epithelial cells (SAEC) and induced pluripotent stem cell (iPSC) clones derived from SAEC. SCLC lines, lung iPSC (Lu-iPSC), and SAEC were further evaluated by DNase I hypersensitivity (DHS-seq). Changes in chromatin accessibility and depths of transcription factor (TF) footprints were quantified using Bivariate analysis of Genomic Footprint. Standard techniques were used to examine growth and tumorigenencity as well as changes in transcriptomes and glucose metabolism of SCLC cells following Nuclear Factor 1C (NFIC) knockdown, and to examine NFIC expression in SCLC cells following exposure to BET inhibitors.

RESULTS

Significant commonality of transcriptomes and DNA methylomes was observed between Lu-iPSC and SCLC; however, this analysis was uninformative regarding pathways unique to lung cancer. Linking results of DNase-seq to RNA-seq enabled identification of networks not previously associated with SCLC. When combined with footprint depth, NFIC, a transcription factor not previously associated with SCLC, had the highest score of occupancy at open chromatin sites. Knockdown of NFIC impaired glucose metabolism, decreased stemness, and inhibited growth of SCLC cells in-vitro and in-vivo. ChIP-seq analysis identified numerous sites occupied by Bromodomain-containing protein 4 (BRD4) in the NFIC promoter region. Knock-down of BRD4 or treatment with Bromodomain and extra-terminal domain (BET) inhibitors (BETi) markedly reduced NFIC expression in SCLC cells and SCLC PDX models. Approximately 8% of genes downregulated by BETi treatment were repressed by NFIC knockdown in SCLC, while 34% of genes repressed following NFIC knockdown were also downregulated in SCLC cells following BETi treatment.

CONCLUSIONS

NFIC is a key TF and possible mediator of transcriptional regulation by BET family proteins in SCLC. Our findings highlight the potential of genome-wide chromatin accessibility analysis for elucidating mechanisms of pulmonary carcinogenesis and identifying novel targets for lung cancer therapy.

Novel biosensor for high-throughput detection of progesterone receptor-interacting endocrine disruptors

Progesterone receptor (PR)-interacting compounds in the environment are associated with serious health hazards. However, methods for their detection in environmental samples are cumbersome. We report a sensitive activity-based biosensor for rapid and reliable screening of progesterone receptor (PR)-interacting endocrine disrupting chemicals (EDCs). The biosensor is a cell line which expresses nuclear mCherry-NF1 and a green fluorescent protein (GFP)-tagged chimera of glucocorticoid receptor (GR) N terminus fused to the ligand binding domain (LBD) of PR (GFP-GR-PR). As this LBD is shared by the PRA and PRB, the biosensor reports on the activation of both PR isoforms. This GFP-GR-PR chimera is cytoplasmic in the absence of hormone and translocates rapidly to the nucleus in response to PR agonists or antagonists in concentration- and time-dependent manner. In live cells, presence of nuclear NF1 label eliminates cell fixation and nuclear staining resulting in efficient screening. The assay can be used in screens for novel PR ligands and PR-interacting contaminants in environmental samples. A limited screen of river water samples indicated a widespread, low-level contamination with PR-interacting contaminants in all tested samples.

The mineralocorticoid receptor forms higher order oligomers upon DNA binding

The prevailing model of steroid hormone nuclear receptor function assumes ligand-induced homodimer formation followed by binding to DNA hormone response elements (HREs). This model has been challenged by evidence showing that the glucocorticoid receptor (GR) forms tetramers upon ligand and DNA binding, which then drive receptor-mediated gene transactivation and transrepression. GR and the closely-related mineralocorticoid receptors (MR) interact to transduce corticosteroid hormone signaling, but whether they share the same quaternary arrangement is unknown. Here, we used a fluorescence imaging technique, Number & Brightness, to study oligomerization in a cell system allowing real-time analysis of receptor-DNA interactions. Agonist-bound MR forms tetramers in the nucleoplasm and higher order oligomers upon binding to HREs. Antagonists form intermediate-size quaternary arrangements, suggesting that large oligomers are essential for function. Divergence between MR and GR quaternary structure is driven by different functionality of known and new multimerization interfaces, which does not preclude formation of heteromers. Thus, influencing oligomerization may be important to selectively modulate corticosteroid signaling.

Transcription Factor Dynamics: One Molecule at a Time

Cells must tightly regulate their gene expression programs and yet rapidly respond to acute biochemical and biophysical cues within their environment. This information is transmitted to the nucleus through various signaling cascades, culminating in the activation or repression of target genes. Transcription factors (TFs) are key mediators of these signals, binding to specific regulatory elements within chromatin. While live-cell imaging has conclusively proven that TF-chromatin interactions are highly dynamic, how such transient interactions can have long-term impacts on developmental trajectories and disease progression is still largely unclear. In this review, we summarize our current understanding of the dynamic nature of TF functions, starting with a historical overview of early live-cell experiments. We highlight key factors that govern TF dynamics and how TF dynamics, in turn, affect downstream transcriptional bursting. Finally, we conclude with open challenges and emerging technologies that will further our understanding of transcriptional regulation.

Long-range gene regulation in hormone-dependent cancer

The human genome is organized into multiple structural layers, ranging from chromosome territories to progressively smaller substructures, such as topologically associating domains (TADs) and chromatin loops. These substructures, collectively referred to as long-range chromatin interactions (LRIs), have a significant role in regulating gene expression. TADs are regions of the genome that harbour groups of genes and regulatory elements that frequently interact with each other and are insulated from other regions, thereby preventing widespread uncontrolled DNA contacts. Chromatin loops formed within TADs through enhancer and promoter interactions are elastic, allowing transcriptional heterogeneity and stochasticity. Over the past decade, it has become evident that the 3D genome structure, also referred to as the chromatin architecture, is central to many transcriptional cellular decisions. In this Review, we delve into the intricate relationship between steroid receptors and LRIs, discussing how steroid receptors interact with and modulate these chromatin interactions. Genetic alterations in the many processes involved in organizing the nuclear architecture are often associated with the development of hormone-dependent cancers. A better understanding of the interplay between architectural proteins and hormone regulatory networks can ultimately be exploited to develop improved approaches for cancer treatment.

SMYD3 represses tumor-intrinsic interferon response in HPV-negative squamous cell carcinoma of the head and neck

Cancers often display immune escape, but the mechanisms are incompletely understood. Herein, we identify SMYD3 as a mediator of immune escape in human papilloma virus (HPV)-negative head and neck squamous cell carcinoma (HNSCC), an aggressive disease with poor response to immunotherapy with pembrolizumab. SMYD3 depletion induces upregulation of multiple type I interferon (IFN) response and antigen presentation machinery genes in HNSCC cells. Mechanistically, SMYD3 binds to and regulates the transcription of UHRF1, encoding for a reader of H3K9me3, which binds to H3K9me3-enriched promoters of key immune-related genes, recruits DNMT1, and silences their expression. SMYD3 further maintains the repression of immune-related genes through intragenic deposition of H4K20me3. In vivo, Smyd3 depletion induces influx of CD8+ T cells and increases sensitivity to anti-programmed death 1 (PD-1) therapy. SMYD3 overexpression is associated with decreased CD8 T cell infiltration and poor response to neoadjuvant pembrolizumab. These data support combining SMYD3 depletion strategies with checkpoint blockade to overcome anti-PD-1 resistance in HPV-negative HNSCC.

The Glucocorticoid Receptor is Required for Efficient Aldosterone-Induced Transcription by the Mineralocorticoid Receptor

The glucocorticoid and mineralocorticoid receptors (GR and MR, respectively) have distinct, yet overlapping physiological and pathophysiological functions. There are indications that both receptors interact functionally and physically, but the precise role of this interdependence is poorly understood. Here, we analyzed the impact of GR co-expression on MR genome-wide chromatin binding and transcriptional responses to aldosterone and glucocorticoids, both physiological ligands of this receptor. Our data show that GR co-expression alters MR genome-wide binding to consensus DNA sequences in a locus- and ligand-specific way. MR binding to consensus DNA sequences is affected by GR. Transcriptional responses of MR in the absence of GR are weak and show poor correlation with chromatin binding. In contrast, co-expression of GR potentiates MR-mediated transcription, particularly in response to aldosterone. Finally, single-molecule tracking of MR suggests that the presence of GR contributes to productive binding of MR/aldosterone complexes to chromatin. Together, our data indicate that co-expression of GR potentiates aldosterone-mediated MR transcriptional activity, even in the absence of glucocorticoids.

Dynamic switching of transcriptional regulators between two distinct low-mobility chromatin states

How chromatin dynamics relate to transcriptional activity remains poorly understood. Using single-molecule tracking, coupled with machine learning, we show that histone H2B and multiple chromatin-bound transcriptional regulators display two distinct low-mobility states. Ligand activation results in a marked increase in the propensity of steroid receptors to bind in the lowest-mobility state. Mutational analysis revealed that interactions with chromatin in the lowest-mobility state require an intact DNA binding domain and oligomerization domains. These states are not spatially separated as previously believed, but individual H2B and bound-TF molecules can dynamically switch between them on time scales of seconds. Single bound-TF molecules with different mobilities exhibit different dwell time distributions, suggesting that the mobility of TFs is intimately coupled with their binding dynamics. Together, our results identify two unique and distinct low-mobility states that appear to represent common pathways for transcription activation in mammalian cells.

The mineralocorticoid receptor forms higher order oligomers upon DNA binding

The prevailing model of steroid hormone nuclear receptor function assumes ligand-induced homodimer formation followed by binding to DNA hormone response elements (HREs). This model has been challenged by evidence showing that the glucocorticoid receptor (GR) forms tetramers upon ligand and DNA binding, which then drive receptor-mediated gene transactivation and transrepression. GR and the closely-related mineralocorticoid receptors (MR) interact to transduce corticosteroid hormone signaling, but whether they share the same quaternary arrangement is unknown. Here, we used a fluorescence imaging technique, Number & Brightness, to study oligomerization in a cell system allowing real-time analysis of receptor-DNA interactions. Agonist-bound MR forms tetramers in the nucleoplasm and higher order oligomers upon binding to HREs. Antagonists form intermediate quaternary arrangements, suggesting that large oligomers are essential for function. Divergence between MR and GR quaternary structure is driven by different functionality of known and new multimerization interfaces, which does not preclude formation of heteromers. Thus, influencing oligomerization may be important to selectively modulate corticosteroid signaling.

Juxtaposition of intensive agriculture, vulnerable aquifers, and mixed chemical/microbial exposures in private-well tapwater in northeast Iowa

In the United States and globally, contaminant exposure in unregulated private-well point-of-use tapwater (TW) is a recognized public-health data gap and an obstacle to both risk-management and homeowner decision making. To help address the lack of data on broad contaminant exposures in private-well TW from hydrologically-vulnerable (alluvial, karst) aquifers in agriculturally-intensive landscapes, samples were collected in 2018-2019 from 47 northeast Iowa farms and analyzed for 35 inorganics, 437 unique organics, 5 in vitro bioassays, and 11 microbial assays. Twenty-six inorganics and 51 organics, dominated by pesticides and related transformation products (35 herbicide-, 5 insecticide-, and 2 fungicide-related), were observed in TW. Heterotrophic bacteria detections were near ubiquitous (94 % of the samples), with detection of total coliform bacteria in 28 % of the samples and growth on at least one putative-pathogen selective media across all TW samples. Health-based hazard index screening levels were exceeded frequently in private-well TW and attributed primarily to inorganics (nitrate, uranium). Results support incorporation of residential treatment systems to protect against contaminant exposure and the need for increased monitoring of rural private-well homes. Continued assessment of unmonitored and unregulated private-supply TW is needed to model contaminant exposures and human-health risks.

Extrachromosomal DNA Amplification Contributes to Small Cell Lung Cancer Heterogeneity and Is Associated with Worse Outcomes

Small-cell lung cancer (SCLC) is an aggressive neuroendocrine lung cancer. Oncogenic MYC amplifications drive SCLC heterogeneity, but the genetic mechanisms of MYC amplification and phenotypic plasticity, characterized by neuroendocrine and nonneuroendocrine cell states, are not known. Here, we integrate whole-genome sequencing, long-range optical mapping, single-cell DNA sequencing, and fluorescence in situ hybridization to find extrachromosomal DNA (ecDNA) as a primary source of SCLC oncogene amplifications and driver fusions. ecDNAs bring to proximity enhancer elements and oncogenes, creating SCLC transcription-amplifying units, driving exceptionally high MYC gene dosage. We demonstrate that cell-free nucleosome profiling can noninvasively detect ecDNA amplifications in plasma, facilitating its genome-wide interrogation in SCLC and other cancers. Altogether, our work provides the first comprehensive map of SCLC ecDNA and describes a new mechanism that governs MYC-driven SCLC heterogeneity. ecDNA-enabled transcriptional flexibility may explain the significantly worse survival outcomes of SCLC harboring complex ecDNA amplifications.

SIGNIFICANCE

MYC drives SCLC progression, but the genetic basis of MYC-driven SCLC evolution is unknown. Using SCLC as a paradigm, we report how ecDNA amplifications function as MYC-amplifying units, fostering tumor plasticity and a high degree of tumor heterogeneity. This article is highlighted in the In This Issue feature, p. 799.

HIV vaccine candidate efficacy in female macaques mediated by cAMP-dependent efferocytosis and V2-specific ADCC

The development of an effective vaccine to protect against HIV acquisition will be greatly bolstered by in-depth understanding of the innate and adaptive responses to vaccination. We report here that the efficacy of DNA/ALVAC/gp120/alum vaccines, based on V2-specific antibodies mediating apoptosis of infected cells (V2-ADCC), is complemented by efferocytosis, a cyclic AMP (cAMP)-dependent antiphlogistic engulfment of apoptotic cells by CD14+ monocytes. Central to vaccine efficacy is the engagement of the CCL2/CCR2 axis and tolerogenic dendritic cells producing IL-10 (DC-10). Epigenetic reprogramming in CD14+ cells of the cyclic AMP/CREB pathway and increased systemic levels of miRNA-139-5p, a negative regulator of expression of the cAMP-specific phosphodiesterase PDE4D, correlated with vaccine efficacy. These data posit that efferocytosis, through the prompt and effective removal of apoptotic infected cells, contributes to vaccine efficacy by decreasing inflammation and maintaining tissue homeostasis.

Identification of a novel GR-ARID1a-P53BP1 protein complex involved in DNA damage repair and cell cycle regulation

ARID1a (BAF250), a component of human SWI/SNF chromatin remodeling complexes, is frequently mutated across numerous cancers, and its loss of function has been putatively linked to glucocorticoid resistance. Here, we interrogate the impact of siRNA knockdown of ARID1a compared to a functional interference approach in the HeLa human cervical cancer cell line. We report that ARID1a knockdown resulted in a significant global decrease in chromatin accessibility in ATAC-Seq analysis, as well as affecting a subset of genome-wide GR binding sites determined by analyzing GR ChIP-Seq data. Interestingly, the specific effects on gene expression were limited to a relatively small subset of glucocorticoid-regulated genes, notably those involved in cell cycle regulation and DNA repair. The vast majority of glucocorticoid-regulated genes were largely unaffected by ARID1a knockdown or functional interference, consistent with a more specific role for ARID1a in glucocorticoid function than previously speculated. Using liquid chromatography-mass spectrometry, we have identified a chromatin-associated protein complex comprising GR, ARID1a, and several DNA damage repair proteins including P53 binding protein 1 (P53BP1), Poly(ADP-Ribose) Polymerase 1 (PARP1), DNA damage-binding protein 1 (DDB1), DNA mismatch repair protein MSH6 and splicing factor proline and glutamine-rich protein (SFPQ), as well as the histone acetyltransferase KAT7, an epigenetic regulator of steroid-dependent transcription, DNA damage repair and cell cycle regulation. Not only was this protein complex ablated with both ARID1a knockdown and functional interference, but spontaneously arising DNA damage was also found to accumulate in a manner consistent with impaired DNA damage repair mechanisms. Recovery from dexamethasone-dependent cell cycle arrest was also significantly impaired. Taken together, our data demonstrate that although glucocorticoids can still promote cell cycle arrest in the absence of ARID1a, the purpose of this arrest to allow time for DNA damage repair is hindered.

Transcription factors TCF-1 and GATA3 are key factors for the epigenetic priming of early innate lymphoid progenitors toward distinct cell fates

The differentiation of innate lymphoid cells (ILCs) from hematopoietic stem cells needs to go through several multipotent progenitor stages. However, it remains unclear whether the fates of multipotent progenitors are predefined by epigenetic states. Here, we report the identification of distinct accessible chromatin regions in all lymphoid progenitors (ALPs), EILPs, and ILC precursors (ILCPs). Single-cell MNase-seq analyses revealed that EILPs contained distinct subpopulations epigenetically primed toward either dendritic cell lineages or ILC lineages. We found that TCF-1 and GATA3 co-bound to the lineage-defining sites for ILCs (LDS-Is), whereas PU.1 binding was enriched in the LDSs for alternative dendritic cells (LDS-As). TCF-1 and GATA3 were indispensable for the epigenetic priming of LDSs at the EILP stage. Our results suggest that the multipotency of progenitor cells is defined by the existence of a heterogeneous population of cells epigenetically primed for distinct downstream lineages, which are regulated by key transcription factors.

The glucocorticoid receptor associates with the cohesin loader NIPBL to promote long-range gene regulation

The cohesin complex is central to chromatin looping, but mechanisms by which these long-range chromatin interactions are formed and persist remain unclear. We demonstrate that interactions between a transcription factor (TF) and the cohesin loader NIPBL regulate enhancer-dependent gene activity. Using mass spectrometry, genome mapping, and single-molecule tracking methods, we demonstrate that the glucocorticoid (GC) receptor (GR) interacts with NIPBL and the cohesin complex at the chromatin level, promoting loop extrusion and long-range gene regulation. Real-time single-molecule experiments show that loss of cohesin markedly diminishes the concentration of TF molecules at specific nuclear confinement sites, increasing TF local concentration and promoting gene regulation. Last, patient-derived acute myeloid leukemia cells harboring cohesin mutations exhibit a reduced response to GCs, suggesting that the GR-NIPBL-cohesin interaction is defective in these patients, resulting in poor response to GC treatment.

Endocrine disruptors of sex hormone activities

Sex hormones, such as androgens, estrogens and progestins are naturally occurring compounds that tightly regulate endocrine systems in a variety of living organisms. Uncontrolled environmental exposure to these hormones or their biological and synthetic mimetics has been widely documented. Furthermore, water contaminants penetrate soil to affect flora, fauna and ultimately humans. Because endocrine systems evolved to respond to very small changes in hormone levels, the low levels found in the environment cannot be ignored. The combined actions of sex hormones with glucocorticoids and other nuclear receptors disruptors creates additional level of complexity including the newly described "dynamic assisted loading" mechanism. We reviewed the extensive literature pertaining to world-wide detection of these disruptors and created a detailed Table on the development and current status of methods used for their analysis.

Nuclear pore protein NUP210 depletion suppresses metastasis through heterochromatin-mediated disruption of tumor cell mechanical response

Mechanical signals from the extracellular microenvironment have been implicated in tumor and metastatic progression. Here, we identify nucleoporin NUP210 as a metastasis susceptibility gene for human estrogen receptor positive (ER+) breast cancer and a cellular mechanosensor. Nup210 depletion suppresses lung metastasis in mouse models of breast cancer. Mechanistically, NUP210 interacts with LINC complex protein SUN2 which connects the nucleus to the cytoskeleton. In addition, the NUP210/SUN2 complex interacts with chromatin via the short isoform of BRD4 and histone H3.1/H3.2 at the nuclear periphery. In Nup210 knockout cells, mechanosensitive genes accumulate H3K27me3 heterochromatin modification, mediated by the polycomb repressive complex 2 and differentially reposition within the nucleus. Transcriptional repression in Nup210 knockout cells results in defective mechanotransduction and focal adhesion necessary for their metastatic capacity. Our study provides an important role of nuclear pore protein in cellular mechanosensation and metastasis.

The involvement of nuclear pore proteins in cellular mechanosensing and metastasis is unclear. Here the authors identify that nuclear pore protein NUP210 promotes metastasis through the interaction with mechanotransducer LINC complex protein and chromatin to regulate mechanosensitive genes.

Activation of Crtc2/Creb1 in skeletal muscle enhances weight loss during intermittent fasting

The Creb-Regulated Transcriptional Coactivator (Crtc) family of transcriptional coregulators drive Creb1-mediated transcription effects on metabolism in many tissues, but the in vivo effects of Crtc2/Creb1 transcription on skeletal muscle metabolism are not known. Skeletal muscle-specific overexpression of Crtc2 (Crtc2 mice) induced greater mitochondrial activity, metabolic flux capacity for both carbohydrates and fats, improved glucose tolerance and insulin sensitivity, and increased oxidative capacity, supported by upregulation of key metabolic genes. Crtc2 overexpression led to greater weight loss during alternate day fasting (ADF), selective loss of fat rather than lean mass, maintenance of higher energy expenditure during the fast and reduced binge-eating during the feeding period. ADF downregulated most of the mitochondrial electron transport genes, and other regulators of mitochondrial function, that were substantially reversed by Crtc2-driven transcription. Glucocorticoids acted with AMPK to drive atrophy and mitophagy, which was reversed by Crtc2/Creb1 signaling. Crtc2/Creb1-mediated signaling coordinates metabolic adaptations in skeletal muscle that explain how Crtc2/Creb1 regulates metabolism and weight loss.

Corticosterone pattern-dependent glucocorticoid receptor binding and transcriptional regulation within the liver

Ultradian glucocorticoid rhythms are highly conserved across mammalian species, however, their functional significance is not yet fully understood. Here we demonstrate that pulsatile corticosterone replacement in adrenalectomised rats induces a dynamic pattern of glucocorticoid receptor (GR) binding at ~3,000 genomic sites in liver at the pulse peak, subsequently not found during the pulse nadir. In contrast, constant corticosterone replacement induced prolonged binding at the majority of these sites. Additionally, each pattern further induced markedly different transcriptional responses. During pulsatile treatment, intragenic occupancy by active RNA polymerase II exhibited pulsatile dynamics with transient changes in enrichment, either decreased or increased depending on the gene, which mostly returned to baseline during the inter-pulse interval. In contrast, constant corticosterone exposure induced prolonged effects on RNA polymerase II occupancy at the majority of gene targets, thus acting as a sustained regulatory signal for both transactivation and repression of glucocorticoid target genes. The nett effect of these differences were consequently seen in the liver transcriptome as RNA-seq analysis indicated that despite the same overall amount of corticosterone infused, twice the number of transcripts were regulated by constant corticosterone infusion, when compared to pulsatile. Target genes that were found to be differentially regulated in a pattern-dependent manner were enriched in functional pathways including carbohydrate, cholesterol, glucose and fat metabolism as well as inflammation, suggesting a functional role for dysregulated glucocorticoid rhythms in the development of metabolic dysfunction.

Adrenal glucocorticoid hormones are released in a characteristic ultradian rhythm that becomes dysregulated during chronic stress, disease, or synthetic corticosteroid treatment. Metabolic dysfunction is a comorbidity associated with all these conditions, but the role that altered glucocorticoid dynamics play is unknown. As the liver is a major site of glucocorticoid action on metabolic homeostasis regulated by the glucocorticoid receptor, we have assessed how different patterns of hormone replacement in adrenalectomised rats differentially regulate gene pathways involved in type II diabetes, cirrhosis, and fatty liver development, via altering the pattern of glucocorticoid receptor binding to regulatory sites. We believe our findings have important implications for therapies that can reproduce the endogenous glucocorticoid rhythm and thus minimize adverse metabolic side-effects in patients.

Power-law behavior of transcription factor dynamics at the single-molecule level implies a continuum affinity model

Single-molecule tracking (SMT) allows the study of transcription factor (TF) dynamics in the nucleus, giving important information regarding the diffusion and binding behavior of these proteins in the nuclear environment. Dwell time distributions obtained by SMT for most TFs appear to follow bi-exponential behavior. This has been ascribed to two discrete populations of TFs-one non-specifically bound to chromatin and another specifically bound to target sites, as implied by decades of biochemical studies. However, emerging studies suggest alternate models for dwell-time distributions, indicating the existence of more than two populations of TFs (multi-exponential distribution), or even the absence of discrete states altogether (power-law distribution). Here, we present an analytical pipeline to evaluate which model best explains SMT data. We find that a broad spectrum of TFs (including glucocorticoid receptor, oestrogen receptor, FOXA1, CTCF) follow a power-law distribution of dwell-times, blurring the temporal line between non-specific and specific binding, suggesting that productive binding may involve longer binding events than previously believed. From these observations, we propose a continuum of affinities model to explain TF dynamics, that is consistent with complex interactions of TFs with multiple nuclear domains as well as binding and searching on the chromatin template.

Mapping multiple endocrine disrupting activities in Virginia rivers using effect-based assays

Water sources are frequently contaminated with natural and anthropogenic substances having known or suspected endocrine disrupting activities; however, these activities are not routinely measured and monitored. Phenotypic bioassays are a promising new approach for detection and quantitation of endocrine disrupting chemicals (EDCs). We developed cell lines expressing fluorescent chimeric constructs capable of detecting environmental contaminants which interact with multiple nuclear receptors. Using these assays, we tested water samples collected in the summers of 2016, 2017 and 2018 from two major Virginia rivers. Samples were concentrated 200× and screened for contaminants interacting with the androgen (AR), glucocorticoid (GR), aryl hydrocarbon (AhR) and thyroid receptors. Among 45 tested sites, over 70% had AR activity and 60% had AhR activity. Many sites were also positive for GR and TRβ activation (22% and 42%, respectively). Multiple sites were positive for more than one type of contaminants, indicating presence of complex mixtures. These activities may negatively impact river ecosystems and consequently human health.

Mechanical Regulation of Transcription: Recent Advances

Mechanotransduction is the ability of a cell to sense mechanical cues from its microenvironment and convert them into biochemical signals to elicit adaptive transcriptional and other cellular responses. Here, we describe recent advances in the field of mechanical regulation of transcription, highlight mechanical regulation of the epigenome as a key novel aspect of mechanotransduction, and describe recent technological advances that could further elucidate the link between mechanical stimuli and gene expression. In this review, we emphasize the importance of mechanotransduction as one of the governing principles of cancer progression, underscoring the need to conduct further studies of the molecular mechanisms involved in sensing mechanical cues and coordinating transcriptional responses.

Novel Assay for Detection of Progesterone Receptor-Interacting Endocrine Disruptors

The presence of progesterone receptor (PR)-interacting compounds in the environment may have serious health consequences for humans and wildlife, but the methods for their detection and monitoring are limited. Here we report the development and testing of a cell line expressing a chimeric construct containing ligand-binding domain of progesterone receptor and green fluorescent protein-tagged domain of the glucocorticoid receptor (GFP-GR-PR) under tetracycline regulation. Unlike the constitutively nuclear PR, this chimera is cytoplasmic in the absence of the ligand and translocates to the nucleus in response to the hormone or its analogues. The GFP-GR-PR chimera maintains specificity for binding to progesterone and does not cross-react with GR-activating hormones. A concentration- and time-dependent translocation in response to progesterone confirmed picomolar sensitivity for detecting PR ligands. Importantly, the assay can detect both agonist and antagonist activities and thus can be used for screening environmental samples for contamination with endocrine disruptors and for drug development. Using this approach, we screened water samples collected at 23 sites along 2 major rivers in Virginia: Mattaponi and Rappahannock Rivers. We detected a low, but reproducible PR-binding activity in 34.8 % of the sites tested. The calculated progesterone equivalent concentration (EQ) in some of these sites reached ~ 0.8 ng/L. The assay provides an effect-based approach for screening PR-interacting endocrine disrupting chemicals regardless of whether they exert agonist or antagonist activities. Either one could be seriously disruptive for the health of humans and wildlife.

An intrinsically disordered region-mediated confinement state contributes to the dynamics and function of transcription factors

Transcription factors (TFs) regulate gene expression by binding to specific consensus motifs within the local chromatin context. The mechanisms by which TFs navigate the nuclear environment as they search for binding sites remain unclear. Here, we used single-molecule tracking and machine-learning-based classification to directly measure the nuclear mobility of the glucocorticoid receptor (GR) in live cells. We revealed two distinct and dynamic low-mobility populations. One accounts for specific binding to chromatin, while the other represents a confinement state that requires an intrinsically disordered region (IDR), implicated in liquid-liquid condensate subdomains. Further analysis showed that the dwell times of both subpopulations follow a power-law distribution, consistent with a broad distribution of affinities on the GR cistrome and interactome. Together, our data link IDRs with a confinement state that is functionally distinct from specific chromatin binding and modulates the transcriptional output by increasing the local concentration of TFs at specific sites.

Meta-analysis of Chromatin Programming by Steroid Receptors

Transcription factors (TFs) must access chromatin to bind to their response elements and regulate gene expression. A widely accepted model proposes that only a special subset of TFs, pioneer factors, can associate with condensed chromatin and initiate chromatin opening. We previously reported that steroid receptors (SRs), not considered pioneer factors, can assist the binding of an archetypal pioneer, the forkhead box protein 1 (FOXA1), at a subset of receptor-activated enhancers. These findings have been challenged recently, with the suggestion that newly acquired data fully support the prevailing pioneer model. Here, we reexamine our results and confirm the original conclusions. We also analyze and discuss a number of available datasets relevant to chromatin penetration by SRs and find a general consensus supporting our original observations. Hence, we propose that chromatin opening at some sites can be initiated by SRs, with a parallel recruitment of factors often treated as having a unique pioneer function. This Matters Arising paper is in response to Glont et al. (2019), published in Cell Reports.

Single-Cell Resolution and Quantitation of Targeted Glucocorticoid Delivery in the Thymus

Glucocorticoids are lipid-soluble hormones that signal via the glucocorticoid receptor (GR), a ligand-dependent transcription factor. Circulating glucocorticoids derive from the adrenals, but it is now apparent that paracrine glucocorticoid signaling occurs in multiple tissues. Effective local glucocorticoid concentrations and whether glucocorticoid delivery can be targeted to specific cell subsets are unknown. We use fluorescence detection of chromatin-associated GRs as biosensors of ligand binding and observe signals corresponding to steroid concentrations over physiological ranges in vitro and in vivo. In the thymus, where thymic epithelial cell (TEC)-synthesized glucocorticoids antagonize negative selection, we find that CD4⁺CD8⁺TCR^hi cells, a small subset responding to self-antigens and undergoing selection, are specific targets of TEC-derived glucocorticoids and are exposed to 3-fold higher levels than other cells. These results demonstrate and quantitate targeted delivery of paracrine glucocorticoids. This approach may be used to assess in situ nuclear receptor signaling in a variety of physiological and pathological contexts.

OR12-03 Mineralocorticoid and Glucocorticoid Receptors Adopt Distinct Quaternary Structures and Can Form Heteromultimers That Affect Chromatin-Binding Profiles

The mineralocorticoid and glucocorticoid receptors (MR and GR) are evolutionarily related nuclear receptors with high sequence conservation and a shared hormone response element (HRE). Both receptors are activated by glucocorticoids, but MR can be selectively activated by aldosterone. Using the imaging technique Number & Brightness (N&B) it has recently been proposed that liganded GR dimers form tetramers upon binding to HREs in live cells. We now show that agonist-bound MR adopts a tetrameric organization in the nucleoplasm and forms complexes with an average of 7 receptor units upon binding an HRE. Interestingly, MR antagonists eplerenone and spironolactone induced intermediate oligomerization arrangements, strongly suggesting that higher order oligomerization is essential for receptor activity. Site-directed mutagenesis and deletion analysis suggest that the N-terminus of MR is a main determinant of higher order oligomerization. Both with corticosterone and aldosterone, GR can incorporate into MR complexes partially displacing MR monomers. Genome-wide chromatin binding studies suggest that the presence of GR in the same cells profoundly change MR interaction with distinct sets of enhancers in a ligand-dependent way, contributing to receptor-specific signaling. Certain genes respond to only one receptor while others respond to both receptors. The interaction of these two closely related receptors has important implications for the mechanisms for glucocorticoid signaling and transcription factors in general.

Pilot study of global endocrine disrupting activity in Iowa public drinking water utilities using cell-based assays

Some anthropogenic substances in drinking water are known or suspected endocrine disrupting compounds (EDCs), but EDCs are not routinely measured. We conducted a pilot study of 10 public drinking water utilities in Iowa, where common contaminants (e.g., pesticides) are suspected EDCs. Raw (untreated) and finished (treated) drinking water samples were collected in spring and fall and concentrated using solid phase extraction. We assessed multiple endocrine disrupting activities using novel mammalian cell-based assays that express nuclear steroid receptors (aryl hydrocarbon [AhR], androgenic [AR], thyroid [TR], estrogenic [ER] and glucocorticoid [GR]). We quantified each receptor's activation relative to negative controls and compared activity by season and utility/sample characteristics. Among 62 samples, 69% had AhR, 52% AR, 3% TR, 2% ER, and 0% GR activity. AhR and AR activities were detected more frequently in spring (p =0 .002 and < 0.001, respectively). AR activity was more common in samples of raw water (p =0 .02) and from surface water utilities (p =0 .05), especially in fall (p =0 .03). Multivariable analyses suggested spring season, surface water, and nitrate and disinfection byproduct concentrations as determinants of bioactivity. Our results demonstrate that AR and AhR activities are commonly found in Iowa drinking water, and that their detection varies by season and utility/sample characteristics. Screening EDCs with cell-based bioassays holds promise for characterizing population exposure to diverse EDCs mixtures.

Transcriptional Bursting and Co-bursting Regulation by Steroid Hormone Release Pattern and Transcription Factor Mobility

Genes are transcribed in a discontinuous pattern referred to as RNA bursting, but the mechanisms regulating this process are unclear. Although many physiological signals, including glucocorticoid hormones, are pulsatile, the effects of transient stimulation on bursting are unknown. Here we characterize RNA synthesis from single-copy glucocorticoid receptor (GR)-regulated transcription sites (TSs) under pulsed (ultradian) and constant hormone stimulation. In contrast to constant stimulation, pulsed stimulation induces restricted bursting centered around the hormonal pulse. Moreover, we demonstrate that transcription factor (TF) nuclear mobility determines burst duration, whereas its bound fraction determines burst frequency. Using 3D tracking of TSs, we directly correlate TF binding and RNA synthesis at a specific promoter. Finally, we uncover a striking co-bursting pattern between TSs located at proximal and distal positions in the nucleus. Together, our data reveal a dynamic interplay between TF mobility and RNA bursting that is responsive to stimuli strength, type, modality, and duration.

Limited Chemical Structural Diversity Found to Modulate Thyroid Hormone Receptor in the Tox21 Chemical Library

BACKGROUND

Thyroid hormone receptors (TRs) are critical endocrine receptors that regulate a multitude of processes in adult and developing organisms, and thyroid hormone disruption is of high concern for neurodevelopmental and reproductive toxicities in particular. To date, only a small number of chemical classes have been identified as possible TR modulators, and the receptors appear highly selective with respect to the ligand structural diversity. Thus, the question of whether TRs are an important screening target for protection of human and wildlife health remains.

OBJECTIVE

Our goal was to evaluate the hypothesis that there is limited structural diversity among environmentally relevant chemicals capable of modulating TR activity via the collaborative interagency Tox21 project.

METHODS

We screened the Tox21 chemical library (8,305 unique structures) in a quantitative high-throughput, cell-based reporter gene assay for TR agonist or antagonist activity. Active compounds were further characterized using additional orthogonal assays, including mammalian one-hybrid assays, coactivator recruitment assays, and a high-throughput, fluorescent imaging, nuclear receptor translocation assay.

RESULTS

Known agonist reference chemicals were readily identified in the TR transactivation assay, but only a single novel, direct agonist was found, the pharmaceutical betamipron. Indirect activation of TR through activation of its heterodimer partner, the retinoid-X-receptor (RXR), was also readily detected by confirmation in an RXR agonist assay. Identifying antagonists with high confidence was a challenge with the presence of significant confounding cytotoxicity and other, non-TR-specific mechanisms common to the transactivation assays. Only three pharmaceuticals-mefenamic acid, diclazuril, and risarestat-were confirmed as antagonists.

DISCUSSION

The results support limited structural diversity for direct ligand effects on TR and imply that other potential target sites in the thyroid hormone axis should be a greater priority for bioactivity screening for thyroid axis disruptors. https://doi.org/10.1289/EHP5314.

Glucocorticoid receptor quaternary structure drives chromatin occupancy and transcriptional outcome

Most transcription factors, including nuclear receptors, are widely modeled as binding regulatory elements as monomers, homodimers, or heterodimers. Recent findings in live cells show that the glucocorticoid receptor NR3C1 (also known as GR) forms tetramers on enhancers, owing to an allosteric alteration induced by DNA binding, and suggest that higher oligomerization states are important for the gene regulatory responses of GR. By using a variant (GRtetra) that mimics this allosteric transition, we performed genome-wide studies using a GR knockout cell line with reintroduced wild-type GR or reintroduced GRtetra. GRtetra acts as a super receptor by binding to response elements not accessible to the wild-type receptor and both induces and represses more genes than GRwt. These results argue that DNA binding induces a structural transition to the tetrameric state, forming a transient higher-order structure that drives both the activating and repressive actions of glucocorticoids.

Conventional and pioneer modes of glucocorticoid receptor interaction with enhancer chromatin in vivo

Glucocorticoid hormone plays a major role in metabolism and disease. The hormone-bound glucocorticoid receptor (GR) binds to a specific set of enhancers in different cell types, resulting in unique patterns of gene expression. We have addressed the role of chromatin structure in GR binding by mapping nucleosome positions in mouse adenocarcinoma cells. Before hormone treatment, GR-enhancers exist in one of three chromatin states: (i) Nucleosome-depleted enhancers that are DNase I-hypersensitive, associated with the Brg1 chromatin remodeler and flanked by nucleosomes incorporating histone H2A.Z. (ii) Nucleosomal enhancers that are DNase I-hypersensitive, marked by H2A.Z and associated with Brg1. (iii) Nucleosomal enhancers that are inaccessible to DNase I, incorporate little or no H2A.Z and lack Brg1. Hormone-induced GR binding results in nucleosome shifts at all types of GR-enhancer, coinciding with increased recruitment of Brg1. We propose that nucleosome-depleted GR-enhancers are formed and maintained by other transcription factors which recruit Brg1 whereas, at nucleosomal enhancers, GR behaves like a pioneer factor, interacting with nucleosomal sites and recruiting Brg1 to remodel the chromatin.

High Quality ATAC-Seq Data Recovered from Cryopreserved Breast Cell Lines and Tissue

DNA accessibility to transcription regulators varies between cells and modulates gene expression patterns. Several “open” chromatin profiling methods that provide valuable insight into the activity of these regulatory regions have been developed. However, their application to clinical samples has been limited despite the discovery that the Analysis of Transposase-Accessible Chromatin followed by sequencing (ATAC-seq) method can be performed using fewer cells than other techniques. Obtaining fresh rather than stored samples and a lack of adequate optimization and quality controls are major barriers to ATAC’s clinical implementation. Here, we describe an optimized ATAC protocol in which we varied nuclear preparation conditions and transposase concentrations and applied rigorous quality control measures before testing fresh, flash frozen, and cryopreserved breast cells and tissue. We obtained high quality data from small cell number. Furthermore, the genomic distribution of sequencing reads, their enrichment at transcription start sites, and transcription factor footprint analyses were similar between cryopreserved and fresh samples. This updated method is applicable to clinical samples, including cells from fine needle aspiration and tissues obtained via core needle biopsy or surgery. Chromatin accessibility analysis using patient samples will greatly expand the range of translational research and personalized medicine by identification of clinically-relevant epigenetic features.

Chromatin reprogramming in breast cancer

Reprogramming of the chromatin landscape is a critical component to the transcriptional response in breast cancer. Effects of sex hormones such as estrogens and progesterone have been well described to have a critical impact on breast cancer proliferation. However, the complex network of the chromatin landscape, enhancer regions and mode of function of steroid receptors (SRs) and other transcription factors (TFs), is an intricate web of signaling and functional processes that is still largely misunderstood at the mechanistic level. In this review, we describe what is currently known about the dynamic interplay between TFs with chromatin and the reprogramming of enhancer elements. Emphasis has been placed on characterizing the different modes of action of TFs in regulating enhancer activity, specifically, how different SRs target enhancer regions to reprogram chromatin in breast cancer cells. In addition, we discuss current techniques employed to study enhancer function at a genome-wide level. Further, we have noted recent advances in live cell imaging technology. These single-cell approaches enable the coupling of population-based assays with real-time studies to address many unsolved questions about SRs and chromatin dynamics in breast cancer.

Dual Roles for Ikaros in Regulation of Macrophage Chromatin State and Inflammatory Gene Expression

Macrophage activation by bacterial LPS leads to induction of a complex inflammatory gene program dependent on numerous transcription factor families. The transcription factor Ikaros has been shown to play a critical role in lymphoid cell development and differentiation; however, its function in myeloid cells and innate immune responses is less appreciated. Using comprehensive genomic analysis of Ikaros-dependent transcription, DNA binding, and chromatin accessibility, we describe unexpected dual repressor and activator functions for Ikaros in the LPS response of murine macrophages. Consistent with the described function of Ikaros as transcriptional repressor, Ikzf1^-/- macrophages showed enhanced induction for select responses. In contrast, we observed a dramatic defect in expression of many delayed LPS response genes, and chromatin immunoprecipitation sequencing analyses support a key role for Ikaros in sustained NF-κB chromatin binding. Decreased Ikaros expression in Ikzf1^+/- mice and human cells dampens these Ikaros-enhanced inflammatory responses, highlighting the importance of quantitative control of Ikaros protein level for its activator function. In the absence of Ikaros, a constitutively open chromatin state was coincident with dysregulation of LPS-induced chromatin remodeling, gene expression, and cytokine responses. Together, our data suggest a central role for Ikaros in coordinating the complex macrophage transcriptional program in response to pathogen challenge.

Bivariate Genomic Footprinting Detects Changes in Transcription Factor Activity

In response to activating signals, transcription factors (TFs) bind DNA and regulate gene expression. TF binding can be measured by protection of the bound sequence from DNase digestion (i.e., footprint). Here, we report that 80% of TF binding motifs do not show a measurable footprint, partly because of a variable cleavage pattern within the motif sequence. To more faithfully portray the effect of TFs on chromatin, we developed an algorithm that captures two TF-dependent effects on chromatin accessibility: footprinting and motif-flanking accessibility. The algorithm, termed bivariate genomic footprinting (BaGFoot), efficiently detects TF activity. BaGFoot is robust to different accessibility assays (DNase-seq, ATAC-seq), all examined peak-calling programs, and a variety of cut bias correction approaches. BaGFoot reliably predicts TF binding and provides valuable information regarding the TFs affecting chromatin accessibility in various biological systems and following various biological events, including in cases where an absolute footprint cannot be determined.

The Three Ds of Transcription Activation by Glucagon: Direct, Delayed, and Dynamic

Upon lowered blood glucose occurring during fasting, glucagon is secreted from pancreatic islets, exerting various metabolic effects to normalize glucose levels. A considerable portion of these effects is mediated by glucagon-activated transcription factors (TFs) in liver. Glucagon directly activates several TFs via immediate cyclic adenosine monophosphate (cAMP)- and calcium-dependent signaling events. Among these TFs, cAMP response element-binding protein (CREB) is a major factor. CREB recruits histone-modifying enzymes and cooperates with other TFs on the chromatin template to increase the rate of gene transcription. In addition to direct signal transduction, the transcriptional effects of glucagon are also influenced by dynamic TF cross talk. Specifically, assisted loading of one TF by a companion TF leads to increased binding and activity. Lastly, transcriptional regulation by glucagon is also exerted by TF cascades by which a primary TF induces the gene expression of secondary TFs that bring about their activity a few hours after the initial glucagon signal. This mechanism of a delayed response may be instrumental in establishing the temporal organization of the fasting response by which distinct metabolic events separate early from prolonged fasting. In this mini-review, we summarize recent advances and critical discoveries in glucagon-dependent gene regulation with a focus on direct TF activation, dynamic TF cross talk, and TF cascades.

Dynamic enhancer function in the chromatin context

Enhancers serve as critical regulatory elements in higher eukaryotic cells. The characterization of enhancer function has evolved primarily from genome-wide methodologies, including chromatin immunoprecipitation (ChIP-seq), DNase-I hypersensitivity (DNase-seq), digital genomic footprinting (DGF), and the chromosome conformation capture techniques (3C, 4C, and Hi-C). These population-based assays average signals across millions of cells and lead to enhancer models characterized by static and sequential binding. More recently, fluorescent microscopy techniques, including fluorescence recovery after photobleaching, fluorescence correlation spectroscopy, and single molecule tracking (SMT), reveal a highly dynamic binding behavior for these factors in live cells. Furthermore, a refined analysis of genomic footprinting suggests that many transcription factors leave minimal or no footprints in chromatin, even when present and active in a given cell type. In this study, we review the implications of these new approaches for an accurate understanding of enhancer function in real time. In vivo SMT, in particular, has recently evolved as a promising methodology to probe enhancer function in live cells. Integration of findings from the many approaches now employed in the study of enhancer function suggest a highly dynamic view for the action of enhancer activating factors, viewed on a time scale of milliseconds to seconds, rather than minutes to hours. WIREs Syst Biol Med 2018, 10:e1390. doi: 10.1002/wsbm.1390 This article is categorized under: Analytical and Computational Methods > Computational Methods Laboratory Methods and Technologies > Genetic/Genomic Methods Laboratory Methods and Technologies > Imaging.

Synergistic gene expression during the acute phase response is characterized by transcription factor assisted loading

The cytokines interleukin 1β and 6 (IL-1β, IL-6) mediate the acute phase response (APR). In liver, they regulate the secretion of acute phase proteins. Using RNA-seq in primary hepatocytes, we show that these cytokines regulate transcription in a bifurcated manner, leading to both synergistic and antagonistic gene expression. By mapping changes in enhancer landscape and transcription factor occupancy (using ChIP-seq), we show that synergistic gene induction is achieved by assisted loading of STAT3 on chromatin by NF-κB. With IL-6 treatment alone, STAT3 does not efficiently bind 20% of its coordinated binding sites. In the presence of IL-1β, NF-κB is activated, binds a subset of enhancers and primes their activity, as evidenced by increasing H3K27ac. This facilitates STAT3 binding and synergistic gene expression. Our findings reveal an enhancer-specific crosstalk whereby NF-κB enables STAT3 binding at some enhancers while perturbing it at others. This model reconciles seemingly contradictory reports of NF-κB-STAT3 crosstalk.

The cytokines IL-1β and IL-6 mediate the systemic acute phase response (APR). Here, the authors provide evidence that these cytokines lead to both synergistic and antagonistic gene expression during APR; synergistic induction occurs by assisted loading of STAT3 on chromatin by NF-κB.

Abstract 4991: Dynamic reprogramming of the chromatin landscape in cancer: Studies in real time

Reprogramming of the chromatin landscape by sequence-specific transcription factors (TF) can influence in the prevalence and progression of several diseases, including cancers. Steroid receptors are a class of hormone-regulated TF that have a major impact on the progression of cancers. The link between the androgen receptor (AR) action in prostate cancer and the glucocorticoid (GR) and estrogen receptor (ER) action in breast cancer is well-known. Furthermore, other TFs such as the pioneer factor FoxA1 has been extensively linked to AR and ER action in these cancers. Although the action of these TF have been widely studied, the knowledge is almost exclusively based on population-based assays such as ChIP-seq and DHS-seq. Hence, the dynamic binding behavior of these factors in living cell remains unclear. For example, the classical pioneer factor model suggests long term static binding events between pioneer TFs, such as FoxA1, and chromatin. However, little is known on the dynamic action of FoxA1 in living cells.

We have examined the behavior of steroid receptors and cofactors, including FoxA1, at the single-molecule level in living cells. These factors were fused with HaloTag, and labeled with the bright and stable fluorophore JF549. Single-molecule tracking (SMT) of these proteins (using HILO microscopy) revealed a highly dynamic binding behavior of TFs. This supports transient rather than stable TF chromatin interactions. Specifically, two distinct binding populations have been observed for all factors tested; fast or slow stops. Fast stops represent genomic scanning of TF, while slow stops, with residence time of 6-14 sec, represent functional binding events at specific response elements. The vast majority of molecules at any given time are either diffusing or exhibit fast stops, while only a small percentage exhibit slow stops. However, hormone activation of AR, GR or ER results in an increase in the percentage of molecules with slow stops compared to the unstimulated state. Thus, our results indicate that on a single-molecule level, only a small proportion TFs are functionally bound at any given time.

In addition to affirming the general model that many TFs are highly dynamic in their chromatin binding activity, SMT experiments in live cells revealed a highly dynamic interaction of FoxA1 with chromatin in vivo. Unexpectedly, we also observed by ChIP-seq and DHS-seq analysis that at subset of genomic sites, the role of pioneer can be reversed, with the steroid receptors serving to enhance binding of FoxA1. We propose a general model, dynamic assisted loading, wherein interactions between TFs and pioneer factors are highly dynamic. Our findings show that the action of key TFs in cancer development are more complex and dynamic in real time than previously reported.

Citation Format: Ville Paakinaho, Diego M. Presman, Erin E. Swinstead, Tina B. Miranda, David A. Ball, Tatiana S. Karpova, Gordon L. Hager. Dynamic reprogramming of the chromatin landscape in cancer: Studies in real time [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4991. doi:10.1158/1538-7445.AM2017-4991

Abstract 3867: Epigenetic siRNA and chemical screens identify SETD8 inhibition as a therapeutic strategy to reactivate p53 in high-risk neuroblastoma

Neuroblastoma (NB) is considered a failure of sympathoadrenal differentiation. High-risk neuroblastoma (HR-NB) is an aggressive pediatric tumor accounting for 15% of all pediatric oncology deaths. Less than 50% of HR-NB patients have long-term survival, despite intense multimodality treatment. Given the paucity of druggable mutations and findings that epigenetic drivers contribute to NB tumorigenesis, we undertook a chromatin-focused siRNA screen to uncover epigenetic regulators critical for survival of high-risk NBs. Of the 400 genes analyzed, high-content Opera imaging identified 53 genes whose loss of expression led to significant decreases in NB cell number with 16 also inducing differentiation. A screen with 21 epigenetic compounds in 8 NB cell lines and 2 non-transformed cell lines prioritized those siRNA hits with active tool compounds in the drug development pipeline. This revealed UNC0379 (targets SETD8) inhibited NB cell growth and identified SETD8 as an important and druggable NB target. SETD8 is the H4K20me1 methyltransferase which regulates DNA replication, chromosome condensation and gene expression. Analysis of primary NB revealed that high expression of SETD8 is associated with poor prognosis in NB (R2 platform ex. Kocak; p=1.4e-07). Levels of SETD8 were not significantly different between Stage 4 MYCN-amp compared to MYCN-WT tumors but high SETD8 levels were only associated with poor prognosis in the Stage 4 MYCN-WT(p=0.03). To understand SETD8-mechanism of action, we performed RNA-seq transcriptome analyses after genetic or pharmacological inhibition of SETD8. Ingenuity Pathway Analysis revealed that SETD8 ablation rescued p53 pro-apoptotic and cell-cycle arrest functions by activating the canonical p53 pathway. Functional studies showed SETD8 methylates p53 (K382) leading to its inactivation. Levels of p53K382me1 are higher in MYCN-WT NB cell lines compared to those with MYCN-amp. Less than 2% of NB tumors have p53 mutations but multiple mechanisms have been identified in MYCN-amp NB that functionally inactivate p53. This study identified that SETD8 inactivates p53 in NB and may be an important mechanism to inactivate p53 in MYCN-WT HR-NB. This subgroup represents 60-70% of HR-NB tumors. SETD8 inhibition led to increases in caspase-dependent cell death only in p53-WT but not -mutant or -null NB cells. Genetic rescue experiments confirmed that SETD8-induced cell death is p53 dependent and p53K382 is important for this activity. Our in vivo xenograft NB models, showed that genetic or pharmacologic (UNC0379) inhibition of SETD8 confers a significant survival advantage. This work identifies that SETD8 is a novel therapeutic target and its inhibition may be especially relevant for the subset of high-risk NB tumors with wildtype MYCN. This is the first in vivo preclinical study showing that targeting SETD8 inhibits tumor growth.

Citation Format: Veronica Veschi, Zhihui Liu, Ty C. Voss, Laurent Ozbun, Berkley Gryder, Chunhua Yan, Ying Hu, Anqi Ma, Jian Jin, Sharlyn J. Mazur, Norris Lam, Barbara K. Souza, Giuseppe Giannini, Gordon L. Hager, Cheryl H. Arrowsmith, Javed Khan, Ettore Appella, Carol Thiele. Epigenetic siRNA and chemical screens identify SETD8 inhibition as a therapeutic strategy to reactivate p53 in high-risk neuroblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3867. doi:10.1158/1538-7445.AM2017-3867

Quantifying transcription factor binding dynamics at the single-molecule level in live cells

Progressive, technological achievements in the quantitative fluorescence microscopy field are allowing researches from many different areas to start unraveling the dynamic intricacies of biological processes inside living cells. From super-resolution microscopy techniques to tracking of individual proteins, fluorescence microscopy is changing our perspective on how the cell works. Fortunately, a growing number of research groups are exploring single-molecule studies in living cells. However, no clear consensus exists on several key aspects of the technique such as image acquisition conditions, or analysis of the obtained data. Here, we describe a detailed approach to perform single-molecule tracking (SMT) of transcription factors in living cells to obtain key binding characteristics, namely their residence time and bound fractions. We discuss different types of fluorophores, labeling density, microscope, cameras, data acquisition, and data analysis. Using the glucocorticoid receptor as a model transcription factor, we compared alternate tags (GFP, mEOS, HaloTag, SNAP-tag, CLIP-tag) for potential multicolor applications. We also examine different methods to extract the dissociation rates and compare them with simulated data. Finally, we discuss several challenges that this exciting technique still faces.

Abstract 5510: Pro inflammatory cytokines induce synergistic gene expression by assisted loading of STAT3 at a subset of primed enhancers

During inflammation, the liver responds to pro-inflammatory signals by producing and secreting acute phase reactants (APR) involved in the inflammatory process. The notable increase in the level of these proteins is due to transcriptional induction that is jointly mediated by the pro-inflammatory cytokines IL-1β, TNFα and IL-6. Although critical to restoring homeostasis during infection or trauma, the hepatic response to these cytokines can be deleterious in pathologies involving chronic inflammation such as cancer.

We set out to decipher how pro-inflammatory cytokines cooperate to elicit a transcriptional response in hepatocytes. We examined the transcriptome of primary mouse hepatocytes treated with different combinations of cytokines and found several gene expression patterns showing a complex cross-talk between cytokine signaling pathways. Specifically, we found that a dual treatment of IL-1β and IL-6 results in a synergistic expression of 164 genes, many of which are APR genes. TNFα had an effect similar to IL-1β, albeit weaker. A sequential rather than a combined treatment attenuated the synergistic effect, suggesting direct cooperation between the two cytokines in inducing transcription (as opposed to a cascade of transcription). Assaying enhancer activity by H3K27ac ChIP-seq revealed the dynamic nature of cytokine-responsive enhancers, with hundreds of enhancers changing their activity level following either treatment. Importantly, enhancers becoming more active upon dual IL-1/IL-6 treatment were enriched with binding sequences for STAT3, CEBP, AP1 and NFkB. Because IL-6 is the major pro-inflammatory cytokine and its transcriptional effects are mediated by the STAT3 transcription factor (TF), we assayed STAT3 binding patterns by ChIP-seq. As expected, STAT3 binding was elevated following IL-6 treatment genome-wide. Surprisingly, in 351 enhancers, STAT3 binding was markedly increased in the dual treatment compared to IL-6 treatment alone, suggesting a supporting role for IL-1β in assisting STAT3 loading at a subset of enhancers. This sub-population of enhancers had lower motif scores compared to unassisted STAT3 binding sites, indicating optimal STAT3 binding at these enhancers might require some facilitation by other factors. Accordingly, these enhancers were less active at the untreated state, but increased in activity following IL-1β treatment, an effect we termed ‘IL-1β priming’.

Our results are consistent with a scenario whereby IL-1β-activated TFs (C/EBPβ, AP1 and/or NFκB) assist STAT3 loading by binding a subset of enhancers, activating them thereby augmenting STAT3 binding. This unexplored mechanism of transcriptional regulation by cytokines may be targeted to alleviate the deleterious effects of pro-tumorigenic inflammatory programs.

Citation Format: Ido Goldstein, Songjoon Baek, Myong-Hee Sung, Gordon L. Hager. Pro inflammatory cytokines induce synergistic gene expression by assisted loading of STAT3 at a subset of primed enhancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5510. doi:10.1158/1538-7445.AM2017-5510

Single-molecule analysis of steroid receptor and cofactor action in living cells

Population-based assays have been employed extensively to investigate the interactions of transcription factors (TFs) with chromatin and are often interpreted in terms of static and sequential binding. However, fluorescence microscopy techniques reveal a more dynamic binding behaviour of TFs in live cells. Here we analyse the strengths and limitations of in vivo single-molecule tracking and performed a comprehensive analysis on the intranuclear dwell times of four steroid receptors and a number of known cofactors. While the absolute residence times estimates can depend on imaging acquisition parameters due to sampling bias, our results indicate that only a small proportion of factors are specifically bound to chromatin at any given time. Interestingly, the glucocorticoid receptor and its cofactors affect each other’s dwell times in an asymmetric manner. Overall, our data indicate transient rather than stable TF-cofactors chromatin interactions at response elements at the single-molecule level.

Transcription factors (TFs) are thought to regulate gene expression by stably binding to target DNA elements. Here, the authors use single-molecule tracking to analyse the dynamic behaviour of steroid receptors TFs and show that most specific interactions with chromatin are transient and dynamic.

Identifying environmental chemicals as agonists of the androgen receptor by using a quantitative high-throughput screening platform

The androgen receptor (AR, NR3C4) is a nuclear receptor whose main function is acting as a transcription factor regulating gene expression for male sexual development and maintaining accessory sexual organ function. It is also a necessary component of female fertility by affecting the functionality of ovarian follicles and ovulation. Pathological processes involving AR include Kennedy's disease and Klinefelter's syndrome, as well as prostate, ovarian, and testicular cancer. Strict regulation of sex hormone signaling is required for normal reproductive organ development and function. Therefore, testing small molecules for their ability to modulate AR is a first step in identifying potential endocrine disruptors. We screened the Tox21 10K compound library in a quantitative high-throughput format to identify activators of AR using two reporter gene cell lines, AR β-lactamase (AR-bla) and AR-luciferase (AR-luc). Seventy-five compounds identified through the primary assay were characterized as potential agonists or inactives through confirmation screens and secondary assays. Biochemical binding and AR nuclear translocation assays were performed to confirm direct binding and activation of AR from these compounds. The top seventeen compounds identified were found to bind to AR, and sixteen of them translocated AR from the cytoplasm into the nucleus. Five potentially novel or not well-characterized AR agonists were discovered through primary and follow-up studies. We have identified multiple AR activators, including known AR agonists such as testosterone, as well as novel/not well-known compounds such as prulifloxacin. The information gained from the current study can be directly used to prioritize compounds for further in-depth toxicological evaluations, as well as their potential to disrupt the endocrine system via AR activation.

Specific genomic cues regulate Cajal body assembly

The assembly of specialized sub-nuclear microenvironments known as nuclear bodies (NBs) is important for promoting efficient nuclear function. In particular, the Cajal body (CB), a prominent NB that facilitates spliceosomal snRNP biogenesis, assembles in response to genomic cues. Here, we detail the factors that regulate CB assembly and structural maintenance. These include the importance of transcription at nucleating gene loci, the grouping of these genes on human chromosomes 1, 6 and 17, as well as cell cycle and biochemical regulation of CB protein function. We also speculate on the correlation between CB formation and RNA splicing levels in neurons and cancer. The timing and location of these specific molecular events is critical to CB assembly and its contribution to genome function. However, further work is required to explore the emerging biophysical characteristics of CB assembly and the impact upon subsequent genome reorganization.

Genome-Wide Identification of Basic Helix-Loop-Helix and NF-1 Motifs Underlying GR Binding Sites in Male Rat Hippocampus

Glucocorticoids regulate hippocampal function in part by modulating gene expression through the glucocorticoid receptor (GR). GR binding is highly cell type specific, directed to accessible chromatin regions established during tissue differentiation. Distinct classes of GR binding sites are dependent on the activity of additional signal-activated transcription factors that prime chromatin toward context-specific organization. We hypothesized a stress context dependency for GR binding in hippocampus as a consequence of rapidly induced stress mediators priming chromatin accessibility. Using chromatin immunoprecipitation sequencing to interrogate GR binding, we found no effect of restraint stress context on GR binding, although analysis of sequences underlying GR binding sites revealed mechanistic detail for hippocampal GR function. We note enrichment of GR binding sites proximal to genes linked to structural and organizational roles, an absence of major tethering partners for GRs, and little or no evidence for binding at negative glucocorticoid response elements. A basic helix-loop-helix motif closely resembling a NeuroD1 or Olig2 binding site was found underlying a subset of GR binding sites and is proposed as a candidate lineage-determining transcription factor directing hippocampal chromatin access for GRs. Of our GR binding sites, 54% additionally contained half-sites for nuclear factor (NF)-1 that we propose as a collaborative or general transcription factor involved in hippocampal GR function. Our findings imply a dose-dependent and context-independent action of GRs in the hippocampus. Alterations in the expression or activity of NF-1/basic helix-loop-helix factors may play an as yet undetermined role in glucocorticoid-related disease susceptibility and outcome by altering GR access to hippocampal binding sites.