Glucocorticoid signaling induces transcriptional memory and universally reversible chromatin changes

Fibroblasts as an in vitro model of circadian genetic and genomic studies

Bipolar disorder (BD) is a heritable disorder characterized by shifts in mood that manifest in manic or depressive episodes. Clinical studies have identified abnormalities of the circadian system in BD patients as a hallmark of underlying pathophysiology. Fibroblasts are a well-established in vitro model for measuring circadian patterns. We set out to examine the underlying genetic architecture of circadian rhythm in fibroblasts, with the goal to assess its contribution to the polygenic nature of BD disease risk. We collected, from primary cell lines of 6 healthy individuals, temporal genomic features over a 48 h period from transcriptomic data (RNA-seq) and open chromatin data (ATAC-seq). The RNA-seq data showed that only a limited number of genes, primarily the known core clock genes such as ARNTL, CRY1, PER3, NR1D2 and TEF display circadian patterns of expression consistently across cell cultures. The ATAC-seq data identified that distinct transcription factor families, like those with the basic helix-loop-helix motif, were associated with regions that were increasing in accessibility over time. Whereas known glucocorticoid receptor target motifs were identified in those regions that were decreasing in accessibility. Further evaluation of these regions using stratified linkage disequilibrium score regression analysis failed to identify a significant presence of them in the known genetic architecture of BD, and other psychiatric disorders or neurobehavioral traits in which the circadian rhythm is affected. In this study, we characterize the biological pathways that are activated in this in vitro circadian model, evaluating the relevance of these processes in the context of the genetic architecture of BD and other disorders, highlighting its limitations and future applications for circadian genomic studies.

Detailing organelle division and segregation in Plasmodium falciparum

The malaria causing parasite, Plasmodium falciparum, replicates through a tightly orchestrated process termed schizogony, where approximately 32 daughter parasites are formed in a single infected red blood cell and thousands of daughter cells in mosquito or liver stages. One-per-cell organelles, such as the mitochondrion and apicoplast, need to be properly divided and segregated to ensure a complete set of organelles per daughter parasites. Although this is highly essential, details about the processes and mechanisms involved remain unknown. We developed a new reporter parasite line that allows visualization of the mitochondrion in blood and mosquito stages. Using high-resolution 3D-imaging, we found that the mitochondrion orients in a cartwheel structure, prior to stepwise, non-geometric division during the last stage of schizogony. Analysis of focused ion beam scanning electron microscopy (FIB-SEM) data confirmed these mitochondrial division stages. Furthermore, these data allowed us to elucidate apicoplast division steps, highlighted its close association with the mitochondrion, and showed putative roles of the centriolar plaques (CPs) in apicoplast segregation. These observations form the foundation for a new detailed mechanistic model of mitochondrial and apicoplast division and segregation during P. falciparum schizogony and pave the way for future studies into the proteins and protein complexes involved in organelle division and segregation.

Human cortical neurogenesis is altered via glucocorticoid-mediated regulation of ZBTB16 expression

Glucocorticoids are important for proper organ maturation, and their levels are tightly regulated during development. Here, we use human cerebral organoids and mice to study the cell-type-specific effects of glucocorticoids on neurogenesis. We show that glucocorticoids increase a specific type of basal progenitors (co-expressing PAX6 and EOMES) that has been shown to contribute to cortical expansion in gyrified species. This effect is mediated via the transcription factor ZBTB16 and leads to increased production of neurons. A phenome-wide Mendelian randomization analysis of an enhancer variant that moderates glucocorticoid-induced ZBTB16 levels reveals causal relationships with higher educational attainment and altered brain structure. The relationship with postnatal cognition is also supported by data from a prospective pregnancy cohort study. This work provides a cellular and molecular pathway for the effects of glucocorticoids on human neurogenesis that relates to lasting postnatal phenotypes.

Oncogenic enhancers prime quiescent metastatic cells to escape NK immune surveillance by eliciting transcriptional memory

Metastasis arises from disseminated tumour cells (DTCs) that are characterized by intrinsic phenotypic plasticity and the capability of seeding to secondary organs. DTCs can remain latent for years before giving rise to symptomatic overt metastasis. In this context, DTCs fluctuate between a quiescent and proliferative state in response to systemic and microenvironmental signals including immune-mediated surveillance. Despite its relevance, how intrinsic mechanisms sustain DTCs plasticity has not been addressed. By interrogating the epigenetic state of metastatic cells, we find that tumour progression is coupled with the activation of oncogenic enhancers that are organized in variable interconnected chromatin domains. This spatial chromatin context leads to the activation of a robust transcriptional response upon repeated exposure to retinoic acid (RA). We show that this adaptive mechanism sustains the quiescence of DTCs through the activation of the master regulator SOX9. Finally, we determine that RA-stimulated transcriptional memory increases the fitness of metastatic cells by supporting the escape of quiescent DTCs from NK-mediated immune surveillance. Overall, these findings highlight the contribution of oncogenic enhancers in establishing transcriptional memories as an adaptive mechanism to reinforce cancer dormancy and immune escape, thus amenable for therapeutic intervention.

Fibroblasts as an in vitro model of circadian genetic and genomic studies: A temporal analysis

Bipolar disorder (BD) is a heritable disorder characterized by shifts in mood that manifest in manic or depressive episodes. Clinical studies have identified abnormalities of the circadian system in BD patients as a hallmark of underlying pathophysiology. Fibroblasts are a well-established in vitro model for measuring circadian patterns. We set out to examine the underlying genetic architecture of circadian rhythm in fibroblasts, with the goal to assess its contribution to the polygenic nature of BD disease risk. We collected, from primary cell lines of 6 healthy individuals, temporal genomic features over a 48 hour period from transcriptomic data (RNA-seq) and open chromatin data (ATAC-seq). The RNA-seq data showed that only a limited number of genes, primarily the known core clock genes such as ARNTL, CRY1, PER3, NR1D2 and TEF display circadian patterns of expression consistently across cell cultures. The ATAC-seq data identified that distinct transcription factor families, like those with the basic helix-loop-helix motif, were associated with regions that were increasing in accessibility over time. Whereas known glucocorticoid receptor target motifs were identified in those regions that were decreasing in accessibility. Further evaluation of these regions using stratified linkage disequilibrium score regression (sLDSC) analysis failed to identify a significant presence of them in the known genetic architecture of BD, and other psychiatric disorders or neurobehavioral traits in which the circadian rhythm is affected. In this study, we characterize the biological pathways that are activated in this in vitro circadian model, evaluating the relevance of these processes in the context of the genetic architecture of BD and other disorders, highlighting its limitations and future applications for circadian genomic studies.

High-throughput screening of glucocorticoid-induced enhancer activity reveals mechanisms of stress-related psychiatric disorders

Exposure to stressful life events increases the risk for psychiatric disorders. Mechanistic insight into the genetic factors moderating the impact of stress can increase our understanding of disease processes. Here, we test 3,662 single nucleotide polymorphisms (SNPs) from preselected expression quantitative trait loci in massively parallel reporter assays to identify genetic variants that modulate the activity of regulatory elements sensitive to glucocorticoids, important mediators of the stress response. Of the tested SNP sequences, 547 were located in glucocorticoid-responsive regulatory elements of which 233 showed allele-dependent activity. Transcripts regulated by these functional variants were enriched for those differentially expressed in psychiatric disorders in the postmortem brain. Phenome-wide Mendelian randomization analysis in 4,439 phenotypes revealed potentially causal associations specifically in neurobehavioral traits, including major depression and other psychiatric disorders. Finally, a functional gene score derived from these variants was significantly associated with differences in the physiological stress response, suggesting that these variants may alter disease risk by moderating the individual set point of the stress response.

Quantification of transcriptome changes to investigate the role of glucocorticoid receptor-RNA binding during dexamethasone treatment

OBJECTIVES

The glucocorticoid receptor (GR) is a well-studied, ligand-activated transcription factor and a common target of anti-inflammatory treatments. Recently, several studies have drawn attention the effects of binding of GR to RNA rather than DNA and the potential implications of this activity for GR function. The objective of our study was to further characterize the relationship between GR function and RNA binding by measuring changes in the glucocorticoid-driven transcriptome in the presence of a GR mutant that exhibited reduced RNA affinity.

DATA DESCRIPTION

GR was activated in three cell lines containing GR constructs (GR-HaloTag). One of the cell lines contained a wild-type GR-HaloTag. Another contained GR-HaloTag with a mutation that reduced RNA affinity and slightly reduced DNA affinity. The third cell line contained GR-HaloTag with a mutation that only slightly reduced DNA affinity. All three cell lines were treated with dexamethasone, a GR agonist. RNA-seq samples were collected every hour for 3 h. Moreover, transcriptome quantification was accomplished via labeling of RNAs transcribed in the final hour of dexamethasone treatment using 4-thiouridine. These labeled RNAs were then purified and sequenced. This data set is the first of its kind for GR and contains valuable insights into the function of RNA binding by GR.

m 6 A mRNA Methylation Regulates Early Pancreatic β-Cell Differentiation

N 6 -methyladenosine (m 6 A) is the most abundant chemical modification in mRNA, and plays important roles in human and mouse embryonic stem cell pluripotency, maintenance, and differentiation. We have recently reported, for the first time, the role of m 6 A in the postnatal control of β-cell function in physiological states and in Type 1 and 2 Diabetes. However, the precise mechanisms by which m 6 A acts to regulate the development of human and mouse β-cells are unexplored. Here, we show that the m 6 A landscape is dynamic during human pancreas development, and that METTL14, one of the m 6 A writer complex proteins, is essential for the early differentiation of both human and mouse β-cells.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

Inducible nitric oxide synthase (iNOS) is necessary for GBP-mediated T. gondii restriction in murine macrophages via vacuole nitration and intravacuolar network collapse

Toxoplasma gondii is an obligate intracellular, protozoan pathogen of rodents and humans. T. gondii's ability to grow within cells and evade cell-autonomous immunity depends on the integrity of the parasitophorous vacuole (PV). Interferon-inducible guanylate binding proteins (GBPs) are central mediators of T. gondii clearance, however, the precise mechanism linking GBP recruitment to the PV and T. gondii restriction is not clear. This knowledge gap is linked to heterogenous GBP-targeting across a population of vacuoles and the lack of tools to selectively purify the intact PV. To identify mediators of parasite clearance associated with GBP2-positive vacuoles, we employed a novel protein discovery tool automated spatially targeted optical micro proteomics (autoSTOMP). This approach identified inducible nitric oxide synthetase (iNOS) enriched at levels similar to the GBPs in infected bone marrow-derived myeloid cells. iNOS expression on myeloid cells was necessary for mice to control T. gondii growth in vivo and survive acute infection. T. gondii infection of IFNγ-primed macrophage was sufficient to robustly induce iNOS expression. iNOS restricted T. gondii infection through nitric oxide synthesis rather than arginine depletion, leading to robust and selective nitration of the PV. Optimal parasite restriction by iNOS and vacuole nitration depended on the chromosome 3 GBPs. Notably, GBP2 recruitment and ruffling of the PV membrane occurred in iNOS knockouts, however, these vacuoles contained dividing parasites. iNOS activity was necessary for the collapse of the intravacuolar network of nanotubular membranes which connects parasites to each other and the host cytosol. Based on these data we conclude reactive nitrogen species generated by iNOS cooperate with the chromosome 3 GBPs to target distinct biology of the PV that are necessary for optimal parasite clearance in murine myeloid cells.

RNA binding by the glucocorticoid receptor attenuates dexamethasone-induced gene activation

The glucocorticoid receptor (GR) is a ligand-activated transcription factor that regulates a suite of genes through direct binding of GR to specific DNA promoter elements. GR also interacts with RNA, but the function of this RNA-binding activity remains elusive. Current models speculate that RNA could repress the transcriptional activity of GR. To investigate the function of the GR-RNA interaction on GR's transcriptional activity, we generated cells that stably express a mutant of GR with reduced RNA binding affinity and treated the cells with the GR agonist dexamethasone. Changes in the dexamethasone-driven transcriptome were quantified using 4-thiouridine labeling of RNAs followed by high-throughput sequencing. We find that while many genes are unaffected, GR-RNA binding is repressive for specific subsets of genes in both dexamethasone-dependent and independent contexts. Genes that are dexamethasone-dependent are activated directly by chromatin-bound GR, suggesting a competition-based repression mechanism in which increasing local concentrations of RNA may compete with DNA for binding to GR at sites of transcription. Unexpectedly, genes that are dexamethasone-independent instead display a localization to specific chromosomal regions, which points to changes in chromatin accessibility or architecture. These results show that RNA binding plays a fundamental role in regulating GR function and highlights potential functions for transcription factor-RNA interactions.

How to tame your genes: mechanisms of inflammatory gene repression by glucocorticoids

Glucocorticoids (GCs) are widely used therapeutic agents to treat a broad range of inflammatory conditions. Their functional effects are elicited by binding to the glucocorticoid receptor (GR), which regulates transcription of distinct gene networks in response to ligand. However, the mechanisms governing various aspects of undesired side effects versus beneficial immunomodulation upon GR activation remain complex and incompletely understood. In this review, we discuss emerging models of inflammatory gene regulation by GR, highlighting GR's regulatory specificity conferred by context-dependent changes in chromatin architecture and transcription factor or co-regulator dynamics. GR controls both gene activation and repression, with the repression mechanism being central to favorable clinical outcomes. We describe current knowledge about 3D genome organization and its role in spatiotemporal transcriptional control by GR. Looking beyond, we summarize the evidence for dynamics in gene regulation by GR through cooperative convergence of epigenetic modifications, transcription factor crosstalk, molecular condensate formation and chromatin looping. Further characterizing these genomic events will reframe our understanding of mechanisms of transcriptional repression by GR.