Chromatin accessibility pre-determines glucocorticoid receptor binding patterns

Endogenous glucocorticoids are required for normal macrophage activation and gastric Helicobacter pylori immunity

Glucocorticoids are steroid hormones well known for their potent anti-inflammatory effects. However, their immunomodulatory properties are multifaceted. Increasing evidence suggests that glucocorticoid signaling promotes effective immunity and that disruption of glucocorticoid signaling impairs immune function. In this study, we conditionally deleted the glucocorticoid receptor (GR) in the myeloid lineage using the LysM-Cre driver (myGRKO). We examined the impact on macrophage activation and gastric immune responses to Helicobacter pylori, the best-known risk factor of gastric cancer. Our results indicate that, compared with wild type (WT), glucocorticoid receptor knockout (GRKO) macrophages exhibited higher expression of proinflammatory genes in steroid-free conditions. However, when challenged in vivo, GRKO macrophages exhibited aberrant chromatin landscapes and impaired proinflammatory gene expression profiles. Moreover, gastric colonization with H. pylori revealed impaired gastric immune responses and reduced T cell recruitment in myGRKO mice. As a result, myGRKO mice were protected from atrophic gastritis and pyloric metaplasia development. These results demonstrate a dual role for glucocorticoid signaling in preparing macrophages to respond to bacterial infection but limiting their pathogenic activation. In addition, our results support that macrophages are critical for gastric H. pylori immunity.NEW & NOTEWORTHY Signaling by endogenous glucocorticoids primes macrophages toward more robust responses to pathogens. Disruption of glucocorticoid signaling caused dysregulation of the chromatin landscape, blunted proinflammatory gene activation upon bacterial challenge, and impaired the gastric inflammatory response to Helicobacter pylori infection.

Epigenetic adaptation drives monocyte differentiation into microglia-like cells upon engraftment into the retina

The identification of specific markers for microglia has been a long-standing challenge. Recently, markers such as P2ry12, TMEM119, and Fcrls have been proposed as microglia-specific and widely used to explore microglial functions within various central nervous system (CNS) contexts. The specificity of these markers was based on the assumption that circulating monocytes retain their distinct signatures even after infiltrating the CNS. However, recent findings reveal that infiltrating monocytes can adopt microglia-like characteristics while maintaining a pro-inflammatory profile upon permanent engraftment in the CNS.In this study, we utilize bone marrow chimeras, single-cell RNA sequencing, ATAC-seq, flow cytometry, and immunohistochemistry to demonstrate that engrafted monocytes acquire expression of established microglia markers-P2ry12, TMEM119, Fcrls-and the pan-myeloid marker Iba1, which has been commonly mischaracterized as microglia-specific. These changes are accompanied by alterations in chromatin accessibility and shifts in chromatin binding motifs that are indicative of microglial identity. Moreover, we show that engrafted monocytes dynamically regulate the expression of CX3CR1, CCR2, Ly6C, and transcription factors PU.1, CTCF, RUNX, AP-1, CEBP, and IRF2, all of which are crucial for shaping microglial identity. This study is the first to illustrate that engrafted monocytes in the retina undergo both epigenetic and transcriptional changes, enabling them to express microglia-like signatures. These findings highlight the need for future research to account for these changes when assessing the roles of monocytes and microglia in CNS pathology.

bHLH transcription factors cooperate with chromatin remodelers to regulate cell fate decisions during Arabidopsis stomatal development

The development of multicellular organisms requires coordinated changes in gene expression that are often mediated by the interaction between transcription factors (TFs) and their corresponding cis-regulatory elements (CREs). During development and differentiation, the accessibility of CREs is dynamically modulated by the epigenome. How the epigenome, CREs, and TFs together exert control over cell fate commitment remains to be fully understood. In the Arabidopsis leaf epidermis, meristemoids undergo a series of stereotyped cell divisions, then switch fate to commit to stomatal differentiation. Newly created or reanalyzed scRNA-seq and ChIP-seq data confirm that stomatal development involves distinctive phases of transcriptional regulation and that differentially regulated genes are bound by the stomatal basic helix-loop-helix (bHLH) TFs. Targets of the bHLHs often reside in repressive chromatin before activation. MNase-seq evidence further suggests that the repressive state can be overcome and remodeled upon activation by specific stomatal bHLHs. We propose that chromatin remodeling is mediated through the recruitment of a set of physical interactors that we identified through proximity labeling-the ATPase-dependent chromatin remodeling SWI/SNF complex and the histone acetyltransferase HAC1. The bHLHs and chromatin remodelers localize to overlapping genomic regions in a hierarchical order. Furthermore, plants with stage-specific knockdown of the SWI/SNF components or HAC1 fail to activate specific bHLH targets and display stomatal development defects. Together, these data converge on a model for how stomatal TFs and epigenetic machinery cooperatively regulate transcription and chromatin remodeling during progressive fate specification.

Multiome Perturb-seq unlocks scalable discovery of integrated perturbation effects on the transcriptome and epigenome

Single-cell CRISPR screens link genetic perturbations to transcriptional states, but high-throughput methods connecting these induced changes to their regulatory foundations are limited. Here we introduce Multiome Perturb-seq, extending single-cell CRISPR screens to simultaneously measure perturbation-induced changes in gene expression and chromatin accessibility. We apply Multiome Perturb-seq in a CRISPRi screen of 13 chromatin remodelers in human RPE-1 cells, achieving efficient assignment of sgRNA identities to single nuclei via an improved method for capturing barcode transcripts from nuclear RNA. We organize expression and accessibility measurements into coherent programs describing the integrated effects of perturbations on cell state, finding that ARID1A and SUZ12 knockdowns induce programs enriched for developmental features. Pseudotime analysis of perturbations connects accessibility changes to changes in gene expression, highlighting the value of multimodal profiling. Overall, our method provides a scalable and simply implemented system to dissect the regulatory logic underpinning cell state.

A female-biased gene expression signature of dominance in cooperatively breeding meerkats

Dominance is a primary determinant of social dynamics and resource access in social animals. Recent studies show that dominance is also reflected in the gene regulatory profiles of peripheral immune cells. However, the strength and direction of this relationship differs across the species and sex combinations investigated, potentially due to variation in the predictors and energetic consequences of dominance status. Here, we investigated the association between social status and gene expression in the blood of wild meerkats (Suricata suricatta; n = 113 individuals), including in response to lipopolysaccharide, Gardiquimod (an agonist of TLR7, which detects single-stranded RNA in vivo) and glucocorticoid stimulation. Meerkats are cooperatively breeding social carnivores in which breeding females physically outcompete other females to suppress reproduction, resulting in high reproductive skew. They therefore present an opportunity to disentangle the effects of social dominance from those of sex per se. We identify a sex-specific signature of dominance, including 1045 differentially expressed genes in females but none in males. Dominant females exhibit elevated activity in innate immune pathways and a larger fold-change response to LPS challenge. Based on these results and a preliminary comparison to other mammals, we speculate that the gene regulatory signature of social status in the immune system depends on the determinants and energetic costs of social dominance, such that it is most pronounced in hierarchies where physical competition is important and reproductive skew is large. Such a pattern has the potential to mediate life history trade-offs between investment in reproduction versus somatic maintenance.

The Glucocorticoid Receptor: Isoforms, Functions, and Contribution to Glucocorticoid Sensitivity

Glucocorticoids exert pleiotropic effects on all tissues to regulate cellular and metabolic homeostasis. Synthetic forms are used therapeutically in a wide range of conditions for their anti-inflammatory benefits, at the cost of dose and duration-dependent side effects. Significant variability occurs between tissues, disease states, and individuals with regard to both the beneficial and deleterious effects. The glucocorticoid receptor (GR) is the site of action for these hormones and a vast body of work has been conducted understanding its function. Traditionally, it was thought that the anti-inflammatory benefits of glucocorticoids were mediated by transrepression of pro-inflammatory transcription factors, while the adverse metabolic effects resulted from direct transactivation. This canonical understanding of the GR function has been brought into question over the past 2 decades with advances in the resolution of scientific techniques, and the discovery of multiple isoforms of the receptor present in most tissues. Here we review the structure and function of the GR, the nature of the receptor isoforms, and the contribution of the receptor to glucocorticoid sensitivity, or resistance in health and disease.

Genome-Wide Association Study for Weight Loss at the End of Dry-Curing of Hams Produced from Purebred Heavy Pigs

Dissecting the genetics of production traits in livestock is of outmost importance, both to understand biological mechanisms underlying those traits and to facilitate the design of selection programs incorporating that information. For the pig industry, traits related to curing are key for protected designation of origin productions. In particular, appropriate ham weight loss after dry-curing ensures high quality of the final product and avoids economic losses. In this study, we analyzed data (N = 410) of ham weight loss after approximately 20 months of dry-curing. The animals used for ham production were purebred pigs belonging to a commercial line. A genome-wide association study (GWAS) of 29,844 SNP markers revealed the polygenic nature of the trait: 221 loci explaining a small percentage of the variance (0.3-1.65%) were identified on almost all Sus scrofa chromosomes. Post-GWAS analyses revealed 32 windows located within regulatory regions and 94 windows located in intronic regions of specific genes. In total, 30 candidate genes encoding receptors and enzymes associated with ham weight loss (MTHFD1L, DUSP8), proteolysis (SPARCL1, MYH8), drip loss (TNNI2), growth (CDCA3, LSP1, CSMD1, AP2A2, TSPAN4), and fat metabolism (AGPAT4, IGF2R, PTDSS2, HRAS, TALDO1, BRSK2, TNNI2, SYT8, GTF2I, GTF2IRD1, LPCAT3, ATN1, GNB3, CMIP, SORCS2, CCSER1, SPP1) were detected.

Glucocorticoid receptor action in prostate cancer: the role of transcription factor crosstalk

Prostate cancer is one of the most prevalent malignancies and is primarily driven by aberrant androgen receptor (AR) signaling. While AR-targeted therapies form the cornerstone of prostate cancer treatment, they often inadvertently activate compensatory pathways, leading to therapy resistance. This resistance is frequently mediated through changes in transcription factor (TF) crosstalk, reshaping gene regulatory programs and ultimately weakening treatment efficacy. Consequently, investigating TF interactions has become crucial for understanding the mechanisms driving therapy-resistant cancers. Recent evidence has highlighted the crosstalk between the glucocorticoid receptor (GR) and AR, demonstrating that GR can induce prostate cancer therapy resistance by replacing the inactivated AR, thereby becoming a driver of the disease. In addition to this oncogenic role, GR has also been shown to act as a tumor suppressor in prostate cancer. Owing to this dual role and the widespread use of glucocorticoids as adjuvant therapy, it is essential to understand GR's actions across different stages of prostate cancer development. In this review, we explore the current knowledge of GR in prostate cancer, with a specific focus on its crosstalk with other TFs. GR can directly and indirectly interact with a variety of TFs, and these interactions vary significantly depending on the type of prostate cancer cells. By highlighting these crosstalk interactions, we aim to provide insights that can guide the research and development of new GR-targeted therapies to mitigate its harmful effects in prostate cancer.

Genome-wide nucleosome and transcription factor responses to genetic perturbations reveal chromatin-mediated mechanisms of transcriptional regulation

Epigenetic mechanisms contribute to gene regulation by altering chromatin accessibility through changes in transcription factor (TF) and nucleosome occupancy throughout the genome. Despite numerous studies focusing on changes in gene expression, the intricate chromatin-mediated regulatory code remains largely unexplored on a comprehensive scale. We address this by employing a factor-agnostic, reverse-genetics approach that uses MNase-seq to capture genome-wide TF and nucleosome occupancies in response to the individual deletion of 201 transcriptional regulators in Saccharomyces cerevisiae, thereby assaying nearly one million mutant-gene interactions. We develop a principled approach to identify and quantify chromatin changes genome-wide, observing differences in TF and nucleosome occupancy that recapitulate well-established pathways identified by gene expression data. We also discover distinct chromatin signatures associated with the up- and downregulation of genes, and use these signatures to reveal regulatory mechanisms previously unexplored in expression-based studies. Finally, we demonstrate that chromatin features are predictive of transcriptional activity and leverage these features to reconstruct chromatin-based transcriptional regulatory networks. Overall, these results illustrate the power of an approach combining genetic perturbation with high-resolution epigenomic profiling; the latter enables a close examination of the interplay between TFs and nucleosomes genome-wide, providing a deeper, more mechanistic understanding of the complex relationship between chromatin organization and transcription.

Epigenetic Mechanisms Modulated by Glucocorticoids With a Focus on Cushing Syndrome

In Cushing syndrome (CS), prolonged exposure to high cortisol levels results in a wide range of devastating effects causing multisystem morbidity. Despite the efficacy of treatment leading to disease remission and clinical improvement, hypercortisolism-induced complications may persist. Since glucocorticoids use the epigenetic machinery as a mechanism of action to modulate gene expression, the persistence of some comorbidities may be mediated by hypercortisolism-induced long-lasting epigenetic changes. Additionally, glucocorticoids influence microRNA expression, which is an important epigenetic regulator as it modulates gene expression without changing the DNA sequence. Evidence suggests that chronically elevated glucocorticoid levels may induce aberrant microRNA expression which may impact several cellular processes resulting in cardiometabolic disorders. The present article reviews the evidence on epigenetic changes induced by (long-term) glucocorticoid exposure. Key aspects of some glucocorticoid-target genes and their implications in the context of CS are described. Lastly, the effects of epigenetic drugs influencing glucocorticoid effects are discussed for their ability to be potentially used as adjunctive therapy in CS.

Nuclear receptor corepressors non-canonically drive glucocorticoid receptor-dependent activation of hepatic gluconeogenesis

Nuclear receptor corepressors (NCoRs) function in multiprotein complexes containing histone deacetylase 3 (HDAC3) to alter transcriptional output primarily through repressive chromatin remodelling at target loci1-5. In the liver, loss of HDAC3 causes a marked hepatosteatosis largely because of de-repression of genes involved in lipid metabolism6,7; however, the individual roles and contribution of other complex members to hepatic and systemic metabolic regulation are unclear. Here we show that adult loss of both NCoR1 and NCoR2 (double knockout (KO)) in hepatocytes phenocopied the hepatomegalic fatty liver phenotype of HDAC3 KO. In addition, double KO livers exhibited a dramatic reduction in glycogen storage and gluconeogenic gene expression that was not observed with hepatic KO of individual NCoRs or HDAC3, resulting in profound fasting hypoglycaemia. This surprising HDAC3-independent activation function of NCoR1 and NCoR2 is due to an unexpected loss of chromatin accessibility on deletion of NCoRs that prevented glucocorticoid receptor binding and stimulatory effect on gluconeogenic genes. These studies reveal an unanticipated, non-canonical activation function of NCoRs that is required for metabolic health.

Glucocorticoid receptor signaling: intricacies and therapeutic opportunities

The glucocorticoid receptor (GR) is a major nuclear receptor (NR) drug target for the treatment of inflammatory disorders and several cancers. Despite the effectiveness of GR ligands, their systemic action triggers a plethora of side effects, limiting long-term use. Here, we discuss new concepts of and insights into GR mechanisms of action to assist in the identification of routes toward enhanced therapeutic benefits. We zoom in on the communication between different GR domains and how this is influenced by different ligands. We detail findings on the interaction between GR and chromatin, and highlight how condensate formation and coregulator confinement can perturb GR transcriptional responses. Last, we discuss the potential of novel ligands and the therapeutic exploitation of crosstalk with other NRs.

Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network

BACKGROUND

The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice.

METHODS

FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription.

RESULTS

FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency.

CONCLUSIONS

Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

Context is key: glucocorticoid receptor and corticosteroid therapeutics in outcomes after traumatic brain injury

Traumatic brain injury (TBI) is a global health burden, and survivors suffer functional and psychiatric consequences that can persist long after injury. TBI induces a physiological stress response by activating the hypothalamic-pituitary-adrenal (HPA) axis, but the effects of injury on the stress response become more complex in the long term. Clinical and experimental evidence suggests long lasting dysfunction of the stress response after TBI. Additionally, pre- and post-injury stress both have negative impacts on outcome following TBI. This bidirectional relationship between stress and injury impedes recovery and exacerbates TBI-induced psychiatric and cognitive dysfunction. Previous clinical and experimental studies have explored the use of synthetic glucocorticoids as a therapeutic for stress-related TBI outcomes, but these have yielded mixed results. Furthermore, long-term steroid treatment is associated with multiple negative side effects. There is a pressing need for alternative approaches that improve stress functionality after TBI. Glucocorticoid receptor (GR) has been identified as a fundamental link between stress and immune responses, and preclinical evidence suggests GR plays an important role in microglia-mediated outcomes after TBI and other neuroinflammatory conditions. In this review, we will summarize GR-mediated stress dysfunction after TBI, highlighting the role of microglia. We will discuss recent studies which target microglial GR in the context of stress and injury, and we suggest that cell-specific GR interventions may be a promising strategy for long-term TBI pathophysiology.

Epigenetic profiling reveals key genes and cis-regulatory networks specific to human parathyroids

In all terrestrial vertebrates, the parathyroid glands are critical regulators of calcium homeostasis and the sole source of parathyroid hormone (PTH). Hyperparathyroidism and hypoparathyroidism are clinically important disorders affecting multiple organs. However, our knowledge regarding regulatory mechanisms governing the parathyroids has remained limited. Here, we present the comprehensive maps of the chromatin landscape of the human parathyroid glands, identifying active regulatory elements and chromatin interactions. These data allow us to define regulatory circuits and previously unidentified genes that play crucial roles in parathyroid biology. We experimentally validate candidate parathyroid-specific enhancers and demonstrate their integration with GWAS SNPs for parathyroid-related diseases and traits. For instance, we observe reduced activity of a parathyroid-specific enhancer of the Calcium Sensing Receptor gene, which contains a risk allele associated with higher PTH levels compared to the wildtype allele. Our datasets provide a valuable resource for unraveling the mechanisms governing parathyroid gland regulation in health and disease.

JMnorm: a novel joint multi-feature normalization method for integrative and comparative epigenomics

Combinatorial patterns of epigenetic features reflect transcriptional states and functions of genomic regions. While many epigenetic features have correlated relationships, most existing data normalization approaches analyze each feature independently. Such strategies may distort relationships between functionally correlated epigenetic features and hinder biological interpretation. We present a novel approach named JMnorm that simultaneously normalizes multiple epigenetic features across cell types, species, and experimental conditions by leveraging information from partially correlated epigenetic features. We demonstrate that JMnorm-normalized data can better preserve cross-epigenetic-feature correlations across different cell types and enhance consistency between biological replicates than data normalized by other methods. Additionally, we show that JMnorm-normalized data can consistently improve the performance of various downstream analyses, which include candidate cis-regulatory element clustering, cross-cell-type gene expression prediction, detection of transcription factor binding and changes upon perturbations. These findings suggest that JMnorm effectively minimizes technical noise while preserving true biologically significant relationships between epigenetic datasets. We anticipate that JMnorm will enhance integrative and comparative epigenomics.

Endogenous glucocorticoids are required for normal macrophage activation and gastric Helicobacter pylori immunity

Glucocorticoids are steroid hormones well-known for their potent anti-inflammatory effects. However, their immunomodulatory properties are multifaceted. Increasing evidence suggests that glucocorticoid signaling promotes effective immunity and that disruption of glucocorticoid signaling impairs immune function. In this study, we conditionally deleted the glucocorticoid receptor (GR) in the myeloid lineage using the LysM-Cre driver (myGRKO). We examined the impact on macrophage activation and gastric immune responses to Helicobacter pylori , the best-known risk factor of gastric cancer. Our results indicate that compared to WT, GRKO macrophages exhibited higher expression of proinflammatory genes in steroid-free conditions. However, when challenged in vivo, GRKO macrophages exhibited aberrant chromatin landscapes and impaired proinflammatory gene expression profiles. Moreover, gastric colonization with Helicobacter revealed impaired gastric immune responses and reduced T cell recruitment in myGRKO mice. As a result, myGRKO mice were protected from atrophic gastritis and pyloric metaplasia development. These results demonstrate a dual role for glucocorticoid signaling in preparing macrophages to respond to bacterial infection but limiting their pathogenic activation. In addition, our results support that macrophages are critical for gastric anti- Helicobacter immunity.

Estrogen Receptor Chromatin Profiling by CUT&RUN

Gonadal steroid hormones, namely, testosterone, progesterone, and estrogens, influence the physiological state of an organism through the regulation of gene transcription. Steroid hormones activate nuclear hormone receptor (HR), transcription factors (TFs), which bind DNA in a tissue- and cell type-specific manner to influence cellular function. Identifying the genomic binding sites of HRs is essential to understanding mechanisms of hormone signaling across tissues and disease contexts. Traditionally, chromatin immunoprecipitation followed by sequencing (ChIP-seq) has been used to map the genomic binding of HRs in cancer cell lines and large tissues. However, ChIP-seq lacks the sensitivity to detect TF binding in small numbers of cells, such as genetically defined neuronal subtypes in the brain. Cleavage Under Targets & Release Under Nuclease (CUT&RUN) resolves most of the technical limitations of ChIP-seq, enabling the detection of protein-DNA interactions with as few as 100-1000 cells. In this chapter, we provide a stepwise CUT&RUN protocol for detecting and analyzing the genome-wide binding of estrogen receptor α (ERα) in mouse brain tissue. The steps described here can be used to identify the genomic binding sites of most TFs in the brain.

Glucocorticoids, their uses, sexual dimorphisms, and diseases: new concepts, mechanisms, and discoveries

The normal stress response in humans is governed by the hypothalamic-pituitary-adrenal (HPA) axis through heightened mechanisms during stress, raising blood levels of the glucocorticoid hormone cortisol. Glucocorticoids are quintessential compounds that balance the proper functioning of numerous systems in the mammalian body. They are also generated synthetically and are the preeminent therapy for inflammatory diseases. They act by binding to the nuclear receptor transcription factor glucocorticoid receptor (GR), which has two main isoforms (GRα and GRβ). Our classical understanding of glucocorticoid signaling is from the GRα isoform, which binds the hormone, whereas GRβ has no known ligands. With glucocorticoids being involved in many physiological and cellular processes, even small disruptions in their release via the HPA axis, or changes in GR isoform expression, can have dire ramifications on health. Long-term chronic glucocorticoid therapy can lead to a glucocorticoid-resistant state, and we deliberate how this impacts disease treatment. Chronic glucocorticoid treatment can lead to noticeable side effects such as weight gain, adiposity, diabetes, and others that we discuss in detail. There are sexually dimorphic responses to glucocorticoids, and women tend to have a more hyperresponsive HPA axis than men. This review summarizes our understanding of glucocorticoids and critically analyzes the GR isoforms and their beneficial and deleterious mechanisms and the sexual differences that cause a dichotomy in responses. We also discuss the future of glucocorticoid therapy and propose a new concept of dual GR isoform agonist and postulate why activating both isoforms may prevent glucocorticoid resistance.

NET formation is a default epigenetic program controlled by PAD4 in apoptotic neutrophils

Neutrophil extracellular traps (NETs) not only counteract bacterial and fungal pathogens but can also promote thrombosis, autoimmunity, and sterile inflammation. The presence of citrullinated histones, generated by the peptidylarginine deiminase 4 (PAD4), is synonymous with NETosis and is considered independent of apoptosis. Mitochondrial- and death receptor-mediated apoptosis promote gasdermin E (GSDME)-dependent calcium mobilization and membrane permeabilization leading to histone H3 citrullination (H3Cit), nuclear DNA extrusion, and cytoplast formation. H3Cit is concentrated at the promoter in bone marrow neutrophils and redistributes in a coordinated process from promoter to intergenic and intronic regions during apoptosis. Loss of GSDME prevents nuclear and plasma membrane disruption of apoptotic neutrophils but prolongs early apoptosis-induced cellular changes to the chromatin and cytoplasmic granules. Apoptotic signaling engages PAD4 in neutrophils, establishing a cellular state that is primed for NETosis, but that occurs only upon membrane disruption by GSDME, thereby redefining the end of life for neutrophils.

Thyroid hormone-regulated chromatin landscape and transcriptional sensitivity of the pituitary gland

Thyroid hormone (3,5,3'-triiodothyronine, T3) is a key regulator of pituitary gland function. The response to T3 is thought to hinge crucially on interactions of nuclear T3 receptors with enhancers but these sites in pituitary chromatin remain surprisingly obscure. Here, we investigate genome-wide receptor binding in mice using tagged endogenous thyroid hormone receptor β (TRβ) and analyze T3-regulated open chromatin using an anterior pituitary-specific Cre driver (Thrbb2Cre). Strikingly, T3 regulates histone modifications and chromatin opening primarily at sites that maintain TRβ binding regardless of T3 levels rather than at sites where T3 abolishes or induces de novo binding. These sites associate more frequently with T3-activated than T3-suppressed genes. TRβ-deficiency blunts T3-regulated gene expression, indicating that TRβ confers transcriptional sensitivity. We propose a model of gene activation in which poised receptor-enhancer complexes facilitate adjustable responses to T3 fluctuations, suggesting a genomic basis for T3-dependent pituitary function or pituitary dysfunction in thyroid disorders.

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.

E47 as a novel glucocorticoid-dependent gene mediating lipid metabolism in patients with endogenous glucocorticoid excess

Purpose

E47 has been identified as a modulating transcription factor of glucocorticoid receptor target genes, its loss protecting mice from metabolic adverse effects of glucocorticoids. We aimed to analyze the role of E47 in patients with endogenous glucocorticoid excess [Cushing's syndrome (CS)] and its association with disorders of lipid and glucose metabolism.

Methods

This is a prospective cohort study including 120 female patients with CS (ACTH-dependent = 79; ACTH-independent = 41) and 26 healthy female controls. Morning whole blood samples after an overnight fast were used to determine E47 mRNA expression levels in patients with overt CS before and 6-12 months after curative surgery. Expression levels were correlated with the clinical phenotype of the patients. Control subjects underwent ACTH stimulation tests and dexamethasone suppression tests to analyze short-term regulation of E47.

Results

E47 gene expression showed significant differences in patient cohorts with overt CS vs. patients in remission (p = 0.0474) and in direct intraindividual comparisons pre- vs. post-surgery (p = 0.0353). ACTH stimulation of controls resulted in a significant decrease of E47 mRNA expression 30 min after i.v. injection compared to baseline measurements. Administration of 1 mg of dexamethasone overnight in controls did not change E47 mRNA expression. E47 gene expression showed a positive correlation with total serum cholesterol (p = 0.0036), low-density lipoprotein cholesterol (p = 0.0157), and waist-arm ratio (p = 0.0138) in patients with CS in remission.

Conclusion

E47 is a GC-dependent gene that is upregulated in GC excess potentially aiming at reducing metabolic glucocorticoid side effects such as dyslipidemia.

Exploring the impact of dexamethasone on gene regulation in myeloma cells

Among glucocorticoids (GCs), dexamethasone (Dex) is widely used in treatment of multiple myelomas. However, despite a definite benefit, all patients relapse. Moreover, the molecular basis of glucocorticoid efficacy remains elusive. To determine genomic response to Dex in myeloma cells, we generated bulk and single-cell multi-omics data and high-resolution contact maps of active enhancers and target genes. We show that a minority of glucocorticoid receptor-binding sites are associated with enhancer activity gains, increased interaction loops, and transcriptional activity. We identified and characterized a predominant enhancer enriched in cohesin (RAD21) and more accessible upon Dex exposure. Analysis of four gene-specific networks revealed the importance of the CTCF-cohesin couple and the synchronization of regulatory sequence openings for efficient transcription in response to Dex. Notably, these epigenomic changes are associated with cell-to-cell transcriptional heterogeneity, in particular, lineage-specific genes. As consequences, BCL2L11-encoding BIM critical for Dex-induced apoptosis and CXCR4 protective from chemotherapy-induced apoptosis are rather up-regulated in different cells. In summary, our work provides new insights into the molecular mechanisms involved in Dex escape.

Cellular allostatic load is linked to increased energy expenditure and accelerated biological aging

Stress triggers anticipatory physiological responses that promote survival, a phenomenon termed allostasis. However, the chronic activation of energy-dependent allostatic responses results in allostatic load, a dysregulated state that predicts functional decline, accelerates aging, and increases mortality in humans. The energetic cost and cellular basis for the damaging effects of allostatic load have not been defined. Here, by longitudinally profiling three unrelated primary human fibroblast lines across their lifespan, we find that chronic glucocorticoid exposure increases cellular energy expenditure by ∼60%, along with a metabolic shift from glycolysis to mitochondrial oxidative phosphorylation (OxPhos). This state of stress-induced hypermetabolism is linked to mtDNA instability, non-linearly affects age-related cytokines secretion, and accelerates cellular aging based on DNA methylation clocks, telomere shortening rate, and reduced lifespan. Pharmacologically normalizing OxPhos activity while further increasing energy expenditure exacerbates the accelerated aging phenotype, pointing to total energy expenditure as a potential driver of aging dynamics. Together, our findings define bioenergetic and multi-omic recalibrations of stress adaptation, underscoring increased energy expenditure and accelerated cellular aging as interrelated features of cellular allostatic load.

Cell Fate Programming by Transcription Factors and Epigenetic Machinery in Stomatal Development

The development of multi-cellular organisms requires coordinated changes in gene expression that are often mediated by the interaction between transcription factors (TFs) and their corresponding cis-regulatory elements (CREs). During development and differentiation, the accessibility of CREs is dynamically modulated by the epigenome. How the epigenome, CREs and TFs together exert control over cell fate commitment remains to be fully understood. In the Arabidopsis leaf epidermis, meristemoids undergo a series of stereotyped cell divisions, then switch fate to commit to stomatal differentiation. Newly created or reanalyzed scRNA-seq and ChIP-seq data confirm that stomatal development involves distinctive phases of transcriptional regulation and that differentially regulated genes are bound by the stomatal basic-helix-loop-helix (bHLH) TFs. Targets of the bHLHs often reside in repressive chromatin before activation. MNase-seq evidence further suggests that the repressive state can be overcome and remodeled upon activation by specific stomatal bHLHs. We propose that chromatin remodeling is mediated through the recruitment of a set of physical interactors that we identified through proximity labeling - the ATPase-dependent chromatin remodeling SWI/SNF complex and the histone acetyltransferase HAC1. The bHLHs and chromatin remodelers localize to overlapping genomic regions in a hierarchical order. Furthermore, plants with stage-specific knock-down of the SWI/SNF components or HAC1 fail to activate specific bHLH targets and display stomatal development defects. Together these data converge on a model for how stomatal TFs and epigenetic machinery cooperatively regulate transcription and chromatin remodeling during progressive fate specification.

Crosstalk between glucocorticoid and mineralocorticoid receptors boosts glucocorticoid-induced killing of multiple myeloma cells

The glucocorticoid receptor (GR) is a crucial drug target in multiple myeloma as its activation with glucocorticoids effectively triggers myeloma cell death. However, as high-dose glucocorticoids are also associated with deleterious side effects, novel approaches are urgently needed to improve GR action in myeloma. Here, we reveal a functional crosstalk between GR and the mineralocorticoid receptor (MR) that plays a role in improved myeloma cell killing. We show that the GR agonist dexamethasone (Dex) downregulates MR levels in a GR-dependent way in myeloma cells. Co-treatment of Dex with the MR antagonist spironolactone (Spi) enhances Dex-induced cell killing in primary, newly diagnosed GC-sensitive myeloma cells. In a relapsed GC-resistant setting, Spi alone induces distinct myeloma cell killing. On a mechanistic level, we find that a GR-MR crosstalk likely arises from an endogenous interaction between GR and MR in myeloma cells. Quantitative dimerization assays show that Spi reduces Dex-induced GR-MR heterodimerization and completely abolishes Dex-induced MR-MR homodimerization, while leaving GR-GR homodimerization intact. Unbiased transcriptomics analyses reveal that c-myc and many of its target genes are downregulated most by combined Dex-Spi treatment. Proteomics analyses further identify that several metabolic hallmarks are modulated most by this combination treatment. Finally, we identified a subset of Dex-Spi downregulated genes and proteins that may predict prognosis in the CoMMpass myeloma patient cohort. Our study demonstrates that GR-MR crosstalk is therapeutically relevant in myeloma as it provides novel strategies for glucocorticoid-based dose-reduction.

H3K36 methylation maintains cell identity by regulating opposing lineage programmes

The epigenetic mechanisms that maintain differentiated cell states remain incompletely understood. Here we employed histone mutants to uncover a crucial role for H3K36 methylation in the maintenance of cell identities across diverse developmental contexts. Focusing on the experimental induction of pluripotency, we show that H3K36M-mediated depletion of H3K36 methylation endows fibroblasts with a plastic state poised to acquire pluripotency in nearly all cells. At a cellular level, H3K36M facilitates epithelial plasticity by rendering fibroblasts insensitive to TGFβ signals. At a molecular level, H3K36M enables the decommissioning of mesenchymal enhancers and the parallel activation of epithelial/stem cell enhancers. This enhancer rewiring is Tet dependent and redirects Sox2 from promiscuous somatic to pluripotency targets. Our findings reveal a previously unappreciated dual role for H3K36 methylation in the maintenance of cell identity by integrating a crucial developmental pathway into sustained expression of cell-type-specific programmes, and by opposing the expression of alternative lineage programmes through enhancer methylation.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

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.

Chromatin accessibility dynamics insight into crosstalk between regulatory landscapes in poplar responses to multiple treatments

Perennial trees develop and coordinate endogenous response signaling pathways, including their crosstalk and convergence, to cope with various environmental stresses which occur simultaneously in most cases. These processes are involved in gene transcriptional regulations that depend on dynamic interactions between regulatory proteins and corresponding chromatin regions, but the mechanisms remain poorly understood in trees. In this study, we detected chromatin regulatory landscapes of poplar under abscisic acid, methyl jasmonate, salicylic acid and sodium chloride (NaCl) treatment, through integrating ATAC-seq and RNA-seq data. Our results showed that the degree of chromatin accessibility for a given gene is closely related to its expression level. However, unlike the gene expression that shows treatment-specific response patterns, changes in chromatin accessibility exhibit high similarities under these treatments. We further proposed and experimentally validated that a homologous gene copy of RESPONSIVE TO DESICCATION 26 mediates the crosstalk between jasmonic acid and NaCl signaling pathways by directly regulating the stress-responsive genes and that circadian clock-related transcription factors like REVEILLE8 play a central role in response of poplar to these treatments. Overall, our study provides a chromatin insight into the molecular mechanism of transcription regulatory networks in response to different environmental stresses and raises the key roles of the circadian clock of poplar to adapt to adverse environments.

A survey on algorithms to characterize transcription factor binding sites

Transcription factors (TFs) are key regulatory proteins that control the transcriptional rate of cells by binding short DNA sequences called transcription factor binding sites (TFBS) or motifs. Identifying and characterizing TFBS is fundamental to understanding the regulatory mechanisms governing the transcriptional state of cells. During the last decades, several experimental methods have been developed to recover DNA sequences containing TFBS. In parallel, computational methods have been proposed to discover and identify TFBS motifs based on these DNA sequences. This is one of the most widely investigated problems in bioinformatics and is referred to as the motif discovery problem. In this manuscript, we review classical and novel experimental and computational methods developed to discover and characterize TFBS motifs in DNA sequences, highlighting their advantages and drawbacks. We also discuss open challenges and future perspectives that could fill the remaining gaps in the field.

Developing OCHROdb, a comprehensive quality checked database of open chromatin regions from sequencing data

International consortia, including ENCODE, Roadmap Epigenomics, Genomics of Gene Regulation and Blueprint Epigenome have made large-scale datasets of open chromatin regions publicly available. While these datasets are extremely useful for studying mechanisms of gene regulation in disease and cell development, they only identify open chromatin regions in individual samples. A uniform comparison of accessibility of the same regulatory sites across multiple samples is necessary to correlate open chromatin accessibility and expression of target genes across matched cell types. Additionally, although replicate samples are available for majority of cell types, a comprehensive replication-based quality checking of individual regulatory sites is still lacking. We have integrated 828 DNase-I hypersensitive sequencing samples, which we have uniformly processed and then clustered their regulatory regions across all samples. We checked the quality of open-chromatin regions using our replication test. This has resulted in a comprehensive, quality-checked database of Open CHROmatin (OCHROdb) regions for 194 unique human cell types and cell lines which can serve as a reference for gene regulatory studies involving open chromatin. We have made this resource publicly available: users can download the whole database, or query it for their genomic regions of interest and visualize the results in an interactive genome browser.

PGE2 alters chromatin through H2A.Z-variant enhancer nucleosome modification to promote hematopoietic stem cell fate

Prostaglandin E2 (PGE2) and 16,16-dimethyl-PGE2 (dmPGE2) are important regulators of hematopoietic stem and progenitor cell (HSPC) fate and offer potential to enhance stem cell therapies [C. Cutler et al. Blood 122, 3074-3081(2013); W. Goessling et al. Cell Stem Cell 8, 445-458 (2011); W. Goessling et al. Cell 136, 1136-1147 (2009)]. Here, we report that PGE2-induced changes in chromatin at enhancer regions through histone-variant H2A.Z permit acute inflammatory gene induction to promote HSPC fate. We found that dmPGE2-inducible enhancers retain MNase-accessible, H2A.Z-variant nucleosomes permissive of CREB transcription factor (TF) binding. CREB binding to enhancer nucleosomes following dmPGE2 stimulation is concomitant with deposition of histone acetyltransferases p300 and Tip60 on chromatin. Subsequent H2A.Z acetylation improves chromatin accessibility at stimuli-responsive enhancers. Our findings support a model where histone-variant nucleosomes retained within inducible enhancers facilitate TF binding. Histone-variant acetylation by TF-associated nucleosome remodelers creates the accessible nucleosome landscape required for immediate enhancer activation and gene induction. Our work provides a mechanism through which inflammatory mediators, such as dmPGE2, lead to acute transcriptional changes and modify HSPC behavior to improve stem cell transplantation.

Glucocorticoid activation of anti-inflammatory macrophages protects against insulin resistance

Insulin resistance (IR) during obesity is linked to adipose tissue macrophage (ATM)-driven inflammation of adipose tissue. Whether anti-inflammatory glucocorticoids (GCs) at physiological levels modulate IR is unclear. Here, we report that deletion of the GC receptor (GR) in myeloid cells, including macrophages in mice, aggravates obesity-related IR by enhancing adipose tissue inflammation due to decreased anti-inflammatory ATM leading to exaggerated adipose tissue lipolysis and severe hepatic steatosis. In contrast, GR deletion in Kupffer cells alone does not alter IR. Co-culture experiments show that the absence of GR in macrophages directly causes reduced phospho-AKT and glucose uptake in adipocytes, suggesting an important function of GR in ATM. GR-deficient macrophages are refractory to alternative ATM-inducing IL-4 signaling, due to reduced STAT6 chromatin loading and diminished anti-inflammatory enhancer activation. We demonstrate that GR has an important function in macrophages during obesity by limiting adipose tissue inflammation and lipolysis to promote insulin sensitivity.

Evaluation of MYBL1 as the master regulator for pachytene spermatocyte genes dysregulated in interspecific hybrid dzo

Misregulation of spermatogenesis transcription factors (TF) in hybrids can lead to misexpression, which is a mechanism for hybrid male sterility (HMS). We used dzo (male offspring of Bos taurus ♂ × Bos grunniens ♀) in bovines to investigate the relationship of the key TF with HMS via RNA sequencing and assay for transposase-accessible chromatin with high-throughput sequencing analyses. RNA sequencing showed that the widespread misexpression in dzo was associated with spermatogenesis-related genes and somatic or progenitor genes. The transition from leptotene or zygotene spermatocytes to pachytene spermatocytes may be the key stage for meiosis arrest in dzo. The analysis of TF-binding motif enrichment revealed that the male meiosis-specific master TF MYB proto-oncogene like 1 (MYBL1, known as A-MYB) motif was enriched on the promoters of downregulated pachytene spermatocyte genes in dzo. Assay for transposase-accessible chromatin with high-throughput sequencing revealed that TF-binding sites for MYBL1, nuclear transcription factor Y, and regulatory factor X were enriched in the low-chromatin accessibility region of dzo. The target genes of the MYBL1-binding motif were associated with meiosis-specific genes and significantly downregulated in dzo testis. The transcription factor MYBL1 may be the candidate master regulator for pachytene spermatocyte genes dysregulated in interspecific HMS dzo. This study reported that a few upstream TF regulation changes might exert a cascading effect downstream in a regulatory network as a mechanism for HMS.

Remembering through the genome: the role of chromatin states in brain functions and diseases

Chromatin is the physical substrate of the genome that carries the DNA sequence and ensures its proper functions and regulation in the cell nucleus. While a lot is known about the dynamics of chromatin during programmed cellular processes such as development, the role of chromatin in experience-dependent functions remains not well defined. Accumulating evidence suggests that in brain cells, environmental stimuli can trigger long-lasting changes in chromatin structure and tri-dimensional (3D) organization that can influence future transcriptional programs. This review describes recent findings suggesting that chromatin plays an important role in cellular memory, particularly in the maintenance of traces of prior activity in the brain. Inspired by findings in immune and epithelial cells, we discuss the underlying mechanisms and the implications for experience-dependent transcriptional regulation in health and disease. We conclude by presenting a holistic view of chromatin as potential molecular substrate for the integration and assimilation of environmental information that may constitute a conceptual basis for future research.

Glucocorticoid activation by HSD11B1 limits T cell-driven interferon signaling and response to PD-1 blockade in melanoma

BACKGROUND

Immune responses against tumors are subject to negative feedback regulation. Immune checkpoint inhibitors (ICIs) blocking Programmed cell death protein 1 (PD-1), a receptor expressed on T cells, or its ligand PD-L1 have significantly improved the treatment of cancer, in particular malignant melanoma. Nevertheless, responses and durability are variables, suggesting that additional critical negative feedback mechanisms exist and need to be targeted to improve therapeutic efficacy.

METHODS

We used different syngeneic melanoma mouse models and performed PD-1 blockade to identify novel mechanisms of negative immune regulation. Genetic gain-of-function and loss-of-function approaches as well as small molecule inhibitor applications were used for target validation in our melanoma models. We analyzed mouse melanoma tissues from treated and untreated mice by RNA-seq, immunofluorescence and flow cytometry to detect changes in pathway activities and immune cell composition of the tumor microenvironment. We analyzed tissue sections of patients with melanoma by immunohistochemistry as well as publicly available single-cell RNA-seq data and correlated target expression with clinical responses to ICIs.

RESULTS

Here, we identified 11-beta-hydroxysteroid dehydrogenase-1 (HSD11B1), an enzyme that converts inert glucocorticoids into active forms in tissues, as negative feedback mechanism in response to T cell immunotherapies. Glucocorticoids are potent suppressors of immune responses. HSD11B1 was expressed in different cellular compartments of melanomas, most notably myeloid cells but also T cells and melanoma cells. Enforced expression of HSD11B1 in mouse melanomas limited the efficacy of PD-1 blockade, whereas small molecule HSD11B1 inhibitors improved responses in a CD8+ T cell-dependent manner. Mechanistically, HSD11B1 inhibition in combination with PD-1 blockade augmented the production of interferon-γ by T cells. Interferon pathway activation correlated with sensitivity to PD-1 blockade linked to anti-proliferative effects on melanoma cells. Furthermore, high levels of HSD11B1, predominantly expressed by tumor-associated macrophages, were associated with poor responses to ICI therapy in two independent cohorts of patients with advanced melanomas analyzed by different methods (scRNA-seq, immunohistochemistry).

CONCLUSION

As HSD11B1 inhibitors are in the focus of drug development for metabolic diseases, our data suggest a drug repurposing strategy combining HSD11B1 inhibitors with ICIs to improve melanoma immunotherapy. Furthermore, our work also delineated potential caveats emphasizing the need for careful patient stratification.

Regulatory and coding sequences of TRNP1 co-evolve with brain size and cortical folding in mammals

Brain size and cortical folding have increased and decreased recurrently during mammalian evolution. Identifying genetic elements whose sequence or functional properties co-evolve with these traits can provide unique information on evolutionary and developmental mechanisms. A good candidate for such a comparative approach is TRNP1, as it controls proliferation of neural progenitors in mice and ferrets. Here, we investigate the contribution of both regulatory and coding sequences of TRNP1 to brain size and cortical folding in over 30 mammals. We find that the rate of TRNP1 protein evolution (ω) significantly correlates with brain size, slightly less with cortical folding and much less with body size. This brain correlation is stronger than for >95% of random control proteins. This co-evolution is likely affecting TRNP1 activity, as we find that TRNP1 from species with larger brains and more cortical folding induce higher proliferation rates in neural stem cells. Furthermore, we compare the activity of putative cis-regulatory elements (CREs) of TRNP1 in a massively parallel reporter assay and identify one CRE that likely co-evolves with cortical folding in Old World monkeys and apes. Our analyses indicate that coding and regulatory changes that increased TRNP1 activity were positively selected either as a cause or a consequence of increases in brain size and cortical folding. They also provide an example how phylogenetic approaches can inform biological mechanisms, especially when combined with molecular phenotypes across several species.

Postoperative systemic inflammatory dysregulation and corticosteroids: a narrative review

In some patients, the inflammatory-immune response to surgical injury progresses to a harmful, dysregulated state. We posit that postoperative systemic inflammatory dysregulation forms part of a pathophysiological response to surgical injury that places patients at increased risk of complications and subsequently prolongs hospital stay. In this narrative review, we have outlined the evolution, measurement and prediction of postoperative systemic inflammatory dysregulation, distinguishing it from a healthy and self-limiting host response. We reviewed the actions of glucocorticoids and the potential for heterogeneous responses to peri-operative corticosteroid supplementation. We have then appraised the evidence highlighting the safety of corticosteroid supplementation, and the potential benefits of high/repeated doses to reduce the risks of major complications and death. Finally, we addressed how clinical trials in the future should target patients at higher risk of peri-operative inflammatory complications, whereby corticosteroid regimes should be tailored to modify not only the a priori risk, but also further adjusted in response to markers of an evolving pathophysiological response.

Genomic glucocorticoid action in embryonic mouse neural stem cells

Prenatal exposure to synthetic glucocorticoids (sGCs) reprograms brain development and predisposes the developing fetus towards potential adverse neurodevelopmental outcomes. Using a mouse model of sGC administration, previous studies show that these changes are accompanied by sexually dimorphic alterations in the transcriptome of neural stem and progenitor cells (NSPCs) derived from the embryonic telencephalon. Because cell type-specific gene expression profiles tightly regulate cell fate decisions and are controlled by a flexible landscape of chromatin domains upon which transcription factors and enhancer elements act, we multiplexed data from four genome-wide assays: RNA-seq, ATAC-seq (assay for transposase accessible chromatin followed by genome wide sequencing), dual cross-linking ChIP-seq (chromatin immunoprecipitation followed by genome wide sequencing), and microarray gene expression to identify novel relationships between gene regulation, chromatin structure, and genomic glucocorticoid receptor (GR) action in NSPCs. These data reveal that GR binds preferentially to predetermined regions of accessible chromatin to influence gene programming and cell fate decisions. In addition, we identify SOX2 as a transcription factor that impacts the genomic response of select GR target genes to sGCs (i.e., dexamethasone) in NSPCs.

Enhancer Clusters Drive Type I Interferon-Induced TRAIL Overexpression in Cancer, and Its Intracellular Protein Accumulation Fails to Induce Apoptosis

Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) is a cytokine produced and secreted by immune cells in response to an infection, often in response to interferon (IFN) stimulation. In cancer, it has also been shown that IFN stimulates the production of TRAIL, and it has been proposed that this TRAIL can induce apoptosis in an autocrine or paracrine manner in different cancer cells. Yet, the mechanism mediating TRAIL upregulation and the implications of TRAIL as an apoptotic molecule in cancer cells are still poorly understood. We show here that in certain cancer cells, TRAIL is upregulated by enhancer clusters, potent genomic regulatory regions containing densely packed enhancers that have combinatorial and additive activity and that are usually found to be associated with cancer-promoting genes. Moreover, we found that TRAIL upregulation by IFNα is mediated by these enhancer clusters in breast and lung cancer cells. Surprisingly, IFNα stimulation leads to the intracellular accumulation of TRAIL protein in these cancer cells. Consequently, this TRAIL is not capable of inducing apoptosis. Our study provides novel insights into the mechanism behind the interferon-mediated upregulation of TRAIL and its protein accumulation in cancer cells. Further investigation is required to understand the role of intracellular TRAIL or depict the mechanisms mediating its apoptosis impairment in cancer cells.

Transcriptional glucocorticoid effects in the brain: Finding the relevant target genes

Glucocorticoids are powerful modulators of brain function. They act via mineralocorticoid and glucocorticoid receptors (MR and GR). These are best understood as transcription factors. Although many glucocorticoid effects depend on the modulation of gene transcription, it is a major challenge to link gene expression to function given the large-scale, apparently pleiotropic genomic responses. The extensive sets of MR and GR target genes are highly specific per cell type, and the brain contains many different (neuronal and non-neuronal) cell types. Next to the set "trait" of cellular context, the "state" of other active signaling pathways will affect MR and GR transcriptional activity. Here, we discuss receptor specificity and contextual factors that determine the transcriptional outcome of MR/GR signaling, experimental possibilities offered by single-cell transcriptomics approaches, and reflect on how to make sense of lists of target genes in relation to understanding the functional effects of steroid receptor activation.

Human zygotic genome activation is initiated from paternal genome

Although parental genomes undergo extensive epigenetic reprogramming to be equalized after fertilization, whether they play different roles in human zygotic genome activation (ZGA) remains unknown. Here, we mapped parental transcriptomes by using human parthenogenetic (PG) and androgenetic (AG) embryos during ZGA. Our data show that human ZGA is launched at the 8-cell stage in AG and bi-parental embryos, but at the morula stage in PG embryos. In contrast, mouse ZGA occurs at the same stage in PG and AG embryos. Mechanistically, primate-specific ZNF675 with AG-specific expression plays a role in human ZGA initiated from paternal genome at the 8-cell stage. AG-specifically expressed LSM1 is also critical for human maternal RNA degradation (MRD) and ZGA. The allelic expressions of ZNF675 and LSM1 are associated with their allelically epigenetic states. Notably, the paternally specific expressions of ZNF675 and LSM1 are also observed in diploid embryos. Collectively, human ZGA is initiated from paternal genome.

Quality assessment and refinement of chromatin accessibility data using a sequence-based predictive model

Chromatin accessibility assays are central to the genome-wide identification of gene regulatory elements associated with transcriptional regulation. However, the data have highly variable quality arising from several biological and technical factors. To surmount this problem, we developed a sequence-based machine learning method to evaluate and refine chromatin accessibility data. Our framework, gapped k-mer SVM quality check (gkmQC), provides the quality metrics for a sample based on the prediction accuracy of the trained models. We tested 886 DNase-seq samples from the ENCODE/Roadmap projects to demonstrate that gkmQC can effectively identify "high-quality" (HQ) samples with low conventional quality scores owing to marginal read depths. Peaks identified in HQ samples are more accurately aligned at functional regulatory elements, show greater enrichment of regulatory elements harboring functional variants, and explain greater heritability of phenotypes from their relevant tissues. Moreover, gkmQC can optimize the peak-calling threshold to identify additional peaks, especially for rare cell types in single-cell chromatin accessibility data.

A multi-omics longitudinal aging dataset in primary human fibroblasts with mitochondrial perturbations

Aging is a process of progressive change. To develop biological models of aging, longitudinal datasets with high temporal resolution are needed. Here we report a multi-omics longitudinal dataset for cultured primary human fibroblasts measured across their replicative lifespans. Fibroblasts were sourced from both healthy donors (n = 6) and individuals with lifespan-shortening mitochondrial disease (n = 3). The dataset includes cytological, bioenergetic, DNA methylation, gene expression, secreted proteins, mitochondrial DNA copy number and mutations, cell-free DNA, telomere length, and whole-genome sequencing data. This dataset enables the bridging of mechanistic processes of aging as outlined by the "hallmarks of aging", with the descriptive characterization of aging such as epigenetic age clocks. Here we focus on bridging the gap for the hallmark mitochondrial metabolism. Our dataset includes measurement of healthy cells, and cells subjected to over a dozen experimental manipulations targeting oxidative phosphorylation (OxPhos), glycolysis, and glucocorticoid signaling, among others. These experiments provide opportunities to test how cellular energetics affect the biology of cellular aging. All data are publicly available at our webtool: https://columbia-picard.shinyapps.io/shinyapp-Lifespan_Study/.

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.

Chromatin accessibility: methods, mechanisms, and biological insights

Access to DNA is a prerequisite to the execution of essential cellular processes that include transcription, replication, chromosomal segregation, and DNA repair. How the proteins that regulate these processes function in the context of chromatin and its dynamic architectures is an intensive field of study. Over the past decade, genome-wide assays and new imaging approaches have enabled a greater understanding of how access to the genome is regulated by nucleosomes and associated proteins. Additional mechanisms that may control DNA accessibility in vivo include chromatin compaction and phase separation - processes that are beginning to be understood. Here, we review the ongoing development of accessibility measurements, we summarize the different molecular and structural mechanisms that shape the accessibility landscape, and we detail the many important biological functions that are linked to chromatin accessibility.