Transcriptomic profiling of osteoarthritis synovial macrophages reveals a tolerized phenotype compounded by a weak corticosteroid response.

OBJECTIVES: It is well-known that long-term osteoarthritis prognosis is not improved by corticosteroid treatments. Here we investigate what could underlie this phenomenon by measuring the short term corticosteroid response of OA-Mf. METHODS: We determined the genome-wide transcriptomic response to corticosteroids of end-stage osteoarthritic joint synovial macrophages (OA-Mf). This was compared with LPS-tolerized and ß-glucan-trained circulating blood monocyte-derived macrophage models. RESULTS: Upon corticosteroid stimulation, the trained and tolerized macrophages significantly alter the abundance of 201 and 257 RNA transcripts, respectively. By contrast, by the same criteria, OA-Mf have a very restricted corticosteroid response of only 12 RNA transcripts. Furthermore, while metalloproteinases 1, -2, -3 and -10 expression clearly distinguish OA-Mf from both the tolerized and trained macrophage models, OA-Mf Interleukin 1 (IL1), chemokine (CXCL) and cytokine (CCL) family member profiles resemble the tolerized macrophage model, with the exception that OA-Mf show high levels of CCL20. CONCLUSION: Terminal osteoarthritis joints therefore harbor macrophages with an inflammatory state that closely resembles the tolerized macrophage state and this is compounded by a weak corticosteroid response capacity that may explain the lack of positive long-term effects of corticosteroid treatment for osteoarthritis patients.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene.

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

Quantification of absolute transcription factor binding affinities in the native chromatin context using BANC-seq.

Transcription factor binding across the genome is regulated by DNA sequence and chromatin features. However, it is not yet possible to quantify the impact of chromatin context on transcription factor binding affinities. Here, we report a method called binding affinities to native chromatin by sequencing (BANC-seq) to determine absolute apparent binding affinities of transcription factors to native DNA across the genome. In BANC-seq, a concentration range of a tagged transcription factor is added to isolated nuclei. Concentration-dependent binding is then measured per sample to quantify apparent binding affinities across the genome. BANC-seq adds a quantitative dimension to transcription factor biology, which enables stratification of genomic targets based on transcription factor concentration and prediction of transcription factor binding sites under non-physiological conditions, such as disease-associated overexpression of (onco)genes. Notably, whereas consensus DNA binding motifs for transcription factors are important to establish high-affinity binding sites, these motifs are not always strictly required to generate nanomolar-affinity interactions in the genome.

The gene regulation knowledge commons: the action area of GREEKC.

As computational modeling becomes more essential to analyze and understand biological regulatory mechanisms, governance of the many databases and knowledge bases that support this domain is crucial to guarantee reliability and interoperability of resources. To address this, the COST Action Gene Regulation Ensemble Effort for the Knowledge Commons (GREEKC, CA15205, www.greekc.org) organized nine workshops in a four-year period, starting September 2016. The workshops brought together a wide range of experts from all over the world working on various steps in the knowledge management process that focuses on understanding gene regulatory mechanisms. The discussions between ontologists, curators, text miners, biologists, bioinformaticians, philosophers and computational scientists spawned a host of activities aimed to standardize and update existing knowledge management workflows and involve end-users in the process of designing the Gene Regulation Knowledge Commons (GRKC). Here the GREEKC consortium describes its main achievements in improving this GRKC.

Retinoic acid signaling drives differentiation toward the absorptive lineage in colorectal cancer.

Retinoic acid (RA) signaling is an important and conserved pathway that regulates cellular proliferation and differentiation. Furthermore, perturbed RA signaling is implicated in cancer initiation and progression. However, the mechanisms by which RA signaling contributes to homeostasis, malignant transformation, and disease progression in the intestine remain incompletely understood. Here, we report, in agreement with previous findings, that activation of the Retinoic Acid Receptor and the Retinoid X Receptor results in enhanced transcription of enterocyte-specific genes in mouse small intestinal organoids. Conversely, inhibition of this pathway results in reduced expression of genes associated with the absorptive lineage. Strikingly, this latter effect is conserved in a human organoid model for colorectal cancer (CRC) progression. We further show that RXR motif accessibility depends on progression state of CRC organoids. Finally, we show that reduced RXR target gene expression correlates with worse CRC prognosis, implying RA signaling as a putative therapeutic target in CRC.

A GO catalogue of human DNA-binding transcription factors.

To control gene transcription, DNA-binding transcription factors recognise specific sequence motifs in gene regulatory regions. A complete and reliable GO annotation of all DNA-binding transcription factors is key to investigating the delicate balance of gene regulation in response to environmental and developmental stimuli. The need for such information is demonstrated by the many lists of transcription factors that have been produced over the past decade. The COST Action Gene Regulation Ensemble Effort for the Knowledge Commons (GREEKC) Consortium brought together experts in the field of transcription with the aim of providing high quality and interoperable gene regulatory data. The Gene Ontology (GO) Consortium provides strict definitions for gene product function, including factors that regulate transcription. The collaboration between the GREEKC and GO Consortia has enabled the application of those definitions to produce a new curated catalogue of over 1400 human DNA-binding transcription factors, that can be accessed at https://www.ebi.ac.uk/QuickGO/targetset/dbTF. This catalogue has facilitated an improvement in the GO annotation of human DNA-binding transcription factors and led to the GO annotation of almost sixty thousand DNA-binding transcription factors in over a hundred species. Thus, this work will aid researchers investigating the regulation of transcription in both biomedical and basic science.

Formalization of gene regulation knowledge using ontologies and gene ontology causal activity models.

Gene regulation computational research requires handling and integrating large amounts of heterogeneous data. The Gene Ontology has demonstrated that ontologies play a fundamental role in biological data interoperability and integration. Ontologies help to express data and knowledge in a machine processable way, which enables complex querying and advanced exploitation of distributed data. Contributing to improve data interoperability in gene regulation is a major objective of the GREEKC Consortium, which aims to develop a standardized gene regulation knowledge commons. GREEKC proposes the use of ontologies and semantic tools for developing interoperable gene regulation knowledge models, which should support data annotation. In this work, we study how such knowledge models can be generated from cartoons of gene regulation scenarios. The proposed method consists of generating descriptions in natural language of the cartoons; extracting the entities from the texts; finding those entities in existing ontologies to reuse as much content as possible, especially from well known and maintained ontologies such as the Gene Ontology, the Sequence Ontology, the Relations Ontology and ChEBI; and implementation of the knowledge models. The models have been implemented using Protégé, a general ontology editor, and Noctua, the tool developed by the Gene Ontology Consortium for the development of causal activity models to capture more comprehensive annotations of genes and link their activities in a causal framework for Gene Ontology Annotations. We applied the method to two gene regulation scenarios and illustrate how to apply the models generated to support the annotation of data from research articles.

Integration of transcription coregulator complexes with sequence-specific DNA-binding factor interactomes.

The domain of transcription regulation has been notoriously difficult to annotate in the Gene Ontology, partly because of the intricacies of gene regulation which involve molecular interactions with DNA as well as amongst protein complexes. The molecular function 'transcription coregulator activity' is a part of the biological process 'regulation of transcription, DNA-templated' that occurs in the cellular component 'chromatin'. It can mechanistically link sequence-specific DNA-binding transcription factor (dbTF) regulatory DNA target sites to coactivator and corepressor target sites through the molecular function 'cis-regulatory region sequence-specific DNA binding'. Many questions arise about transcription coregulators (coTF). Here, we asked how many unannotated, putative coregulators can be identified in protein complexes? Therefore, we mined the CORUM and hu.MAP protein complex databases with known and strongly presumed human transcription coregulators. In addition, we trawled the BioGRID and IntAct molecular interaction databases for interactors of the known 1457 human dbTFs annotated by the GREEKC and GO consortia. This yielded 1093 putative transcription factor coregulator complex subunits, of which 954 interact directly with a dbTF. This substantially expands the set of coTFs that could be annotated to 'transcription coregulator activity' and sets the stage for renewed annotation and wet-lab research efforts. To this end, we devised a prioritisation score based on existing GO annotations of already curated transcription coregulators as well as interactome representation. Since all the proteins that we mined are parts of protein complexes, we propose to concomitantly engage in annotation of the putative transcription coregulator-containing complexes in the Complex Portal database.

Sequence Ontology terminology for gene regulation.

The Sequence Ontology (SO) is a structured, controlled vocabulary that provides terms and definitions for genomic annotation. The Gene Regulation Ensemble Effort for the Knowledge Commons (GREEKC) initiative has gathered input from many groups of researchers, including the SO, the Gene Ontology (GO), and gene regulation experts, with the goal of curating information about how gene expression is regulated at the molecular level. Here we discuss recent updates to the SO reflecting current knowledge. We have developed more accurate human-readable terms (also known as classes), including new definitions, and relationships related to the expression of genes. New findings continue to give us insight into the biology of gene regulation, including the order of events, and participants in those events. These updates to the SO support logical reasoning with the current understanding of gene expression regulation at the molecular level.

Gene Ontology representation for transcription factor functions.

Transcription plays a central role in defining the identity and functionalities of cells, as well as in their responses to changes in the cellular environment. The Gene Ontology (GO) provides a rigorously defined set of concepts that describe the functions of gene products. A GO annotation is a statement about the function of a particular gene product, represented as an association between a gene product and the biological concept a GO term defines. Critically, each GO annotation is based on traceable scientific evidence. Here, we describe the different GO terms that are associated with proteins involved in transcription and its regulation, focusing on the standard of evidence required to support these associations. This article is intended to help users of GO annotations understand how to interpret the annotations and can contribute to the consistency of GO annotations. We distinguish between three classes of activities involved in transcription or directly regulating it - general transcription factors, DNA-binding transcription factors, and transcription co-regulators.

Fibronectin-functionalization of 3D collagen networks supports immune tolerance and inflammation suppression in human monocyte-derived macrophages.

The extracellular matrix (ECM) is dynamically reorganized during wound healing. Concomitantly, recruited monocytes differentiate into macrophages. However, the role of the wound's ECM during this transition remain to be fully understood. Fibronectin is a multifunctional glycoprotein present in early wound ECM with a potential immunomodulatory role during monocyte-to-macrophage differentiation. Hence, to investigate the impact of fibronectin during this differentiation step, 3D fibrillar collagen type I networks with or without fibronectin-functionalization were engineered with defined topology (fibril and pore diameter: 0.8 µm; 7 µm) and amount of adsorbed fibronectin (0.15 µg per µg collagen). Primary, human monocytes were then differentiated into macrophages inside these networks. The immunological imprinting of the resulting macrophages was monitored by means of the expression of FABP4, CLEC4E, SLC2A6, and SOD2 which discriminate naïve and tolerized macrophages, as well pro-inflammatory (M1) and anti-inflammatory (M2) macrophage polarization. The analyses indicate that fibronectin-functionalization of collagen I networks induces macrophage tolerance rather than M1 or M2 macrophage phenotypes. This finding was confirmed by release profiles of pro- and anti-inflammatory cytokines such as IL6, IL8, CXCL10, and IL10. Nevertheless, upon LPS challenge, immune suppression by fibronectin was overridden since these macrophages could then deploy an efficient immune response. Our results therefore provide new perspectives in biomaterial science of wound healing scaffolds and the design of instructive materials for human monocyte-derived cells.

Spatial patterns of CTCF sites define the anatomy of TADs and their boundaries.

BACKGROUND: Topologically associating domains (TADs) are genomic regions of self-interaction. Additionally, it is known that TAD boundaries are enriched in CTCF binding sites. In turn, CTCF sites are known to be asymmetric, whereby the convergent configuration of a pair of CTCF sites leads to the formation of a chromatin loop in vivo. However, to date, it has been unclear how to reconcile TAD structure with CTCF-based chromatin loops. RESULTS: We approach this problem by analysing CTCF binding site strengths and classifying clusters of CTCF sites along the genome on the basis of their relative orientation. Analysis of CTCF site orientation classes as a function of their spatial distribution along the human genome reveals that convergent CTCF site clusters are depleted while divergent CTCF clusters are enriched in the 5- to 100-kb range. We then analyse the distribution of CTCF binding sites as a function of TAD boundary conservation across seven primary human blood cell types. This reveals divergent CTCF site enrichment at TAD boundaries. Furthermore, convergent arrays of CTCF sites separate the left and right sections of TADs that harbour internal CTCF sites, resulting in unequal TAD 'halves'. CONCLUSIONS: The orientation-based CTCF binding site cluster classification that we present reconciles TAD boundaries and CTCF site clusters in a mechanistically elegant fashion. This model suggests that the emergent structure of nuclear chromatin in the form of TADs relies on the obligate alternation of divergent and convergent CTCF site clusters that occur at different length scales along the genome.

The alarmin S100A9 hampers osteoclast differentiation from human circulating precursors by reducing the expression of RANK.

The alarmin S100A8/A9 is implicated in sterile inflammation-induced bone resorption and has been shown to increase the bone-resorptive capacity of mature osteoclasts. Here, we investigated the effects of S100A9 on osteoclast differentiation from human CD14+ circulating precursors. Hereto, human CD14+ monocytes were isolated and differentiated toward osteoclasts with M-CSF and receptor activator of NF-κB (RANK) ligand (RANKL) in the presence or absence of S100A9. Tartrate-resistant acid phosphatase staining showed that exposure to S100A9 during monocyte-to-osteoclast differentiation strongly decreased the numbers of multinucleated osteoclasts. This was underlined by a decreased resorption of a hydroxyapatite-like coating. The thus differentiated cells showed a high mRNA and protein production of proinflammatory factors after 16 h of exposure. In contrast, at d 4, the cells showed a decreased production of the osteoclast-promoting protein TNF-α. Interestingly, S100A9 exposure during the first 16 h of culture only was sufficient to reduce osteoclastogenesis. Using fluorescently labeled RANKL, we showed that, within this time frame, S100A9 inhibited the M-CSF-mediated induction of RANK. Chromatin immunoprecipitation showed that this was associated with changes in various histone marks at the epigenetic level. This S100A9-induced reduction in RANK was in part recovered by blocking TNF-α but not IL-1. Together, our data show that S100A9 impedes monocyte-to-osteoclast differentiation, probably via a reduction in RANK expression.-Di Ceglie, I., Blom, A. B., Davar, R., Logie, C., Martens, J. H. A., Habibi, E., Böttcher, L.-M., Roth, J., Vogl, T., Goodyear, C. S., van der Kraan, P. M., van Lent, P. L., van den Bosch, M. H. The alarmin S100A9 hampers osteoclast differentiation from human circulating precursors by reducing the expression of RANK.

Epigenetic modulation of a hardwired 3D chromatin landscape in two naive states of pluripotency.

The mechanisms underlying enhancer activation and the extent to which enhancer-promoter rewiring contributes to spatiotemporal gene expression are not well understood. Using integrative and time-resolved analyses we show that the extensive transcriptome and epigenome resetting during the conversion between 'serum' and '2i' states of mouse embryonic stem cells (ESCs) takes place with minimal enhancer-promoter rewiring that becomes more evident in primed-state pluripotency. Instead, differential gene expression is strongly linked to enhancer activation via H3K27ac. Conditional depletion of transcription factors and allele-specific enhancer analysis reveal an essential role for Esrrb in H3K27 acetylation and activation of 2i-specific enhancers. Restoration of a polymorphic ESRRB motif using CRISPR-Cas9 in a hybrid ESC line restores ESRRB binding and enhancer H3K27ac in an allele-specific manner but has no effect on chromatin interactions. Our study shows that enhancer activation in serum- and 2i-ESCs is largely driven by transcription factor binding and epigenetic marking in a hardwired network of chromatin interactions.

Publisher Correction: Epigenetic modulation of a hardwired 3D chromatin landscape in two naive states of pluripotency.

In the version of the article originally published, extra lines were displayed in Fig. 7. Fig. 7a contained a solid black line that extended into panel b, and Fig. 7c contained two extra scale bars on the left. These have been removed from the figure. The errors have been corrected in the HTML and PDF versions of the article.

Underreported and unknown student harassment at the Faculty of Science.

Reports of sexual harassment at medical faculties throughout the world, including the Radboud University, raised the question how prevalent this is at the Faculty of Science. We performed a survey among students to assess their experiences with harassment. This questionnaire consisted of questions from the EGERA survey, a questionnaire held among staff of multiple European Universities. We found that 9% of the respondents had observed or experienced harassment at the Faculty. Hardly any of these cases were reported to one of the institutional services. Moreover, most students did not now any of the provided services. We therefore suggest raising awareness on harassment and to make students more familiar with the trust person.

Extensive epigenomic integration of the glucocorticoid response in primary human monocytes and in vitro derived macrophages.

Glucocorticoid receptor is a transcription factor that is ubiquitously expressed. Glucocorticoids are circadian steroids that regulate a wide range of bodily functions, including immunity. Here we report that synthetic glucocorticoids affect 1035 mRNAs in isolated healthy human blood monocytes but only 165 in the respective six day-old monocyte-derived macrophages. The majority of the glucocorticoid response in monocytes concerns genes that are dynamic upon monocyte to macrophage differentiation, whereby macrophage-like mRNA levels are often reached in monocytes within four hours of treatment. Concomitantly, over 5000 chromosomal H3K27ac regions undergo remodelling, of which 60% involve increased H3K27ac signal. We find that chromosomal glucocorticoid receptor binding sites correlate with positive but not with negative local epigenomic effects. To investigate further we assigned our data to topologically associating domains (TADs). This shows that about 10% of macrophage TADs harbour at least one GR binding site and that half of all the glucocorticoid-induced H3K27ac regions are confined to these TADs. Our analyses are therefore consistent with the notion that TADs naturally accommodate information from sets of distal glucocorticoid response elements.

Mutant p63 Affects Epidermal Cell Identity through Rewiring the Enhancer Landscape.

Transcription factor p63 is a key regulator of epidermal keratinocyte proliferation and differentiation. Mutations in the p63 DNA-binding domain are associated with ectrodactyly, ectodermal dysplasia, and cleft lip/palate (EEC) syndrome. However, the underlying molecular mechanism of these mutations remains unclear. Here, we characterized the transcriptome and epigenome of p63 mutant keratinocytes derived from EEC patients. The transcriptome of p63 mutant keratinocytes deviated from the normal epidermal cell identity. Epigenomic analyses showed an altered enhancer landscape in p63 mutant keratinocytes contributed by loss of p63-bound active enhancers and unexpected gain of enhancers. The gained enhancers were frequently bound by deregulated transcription factors such as RUNX1. Reversing RUNX1 overexpression partially rescued deregulated gene expression and the altered enhancer landscape. Our findings identify a disease mechanism whereby mutant p63 rewires the enhancer landscape and affects epidermal cell identity, consolidating the pivotal role of p63 in controlling the enhancer landscape of epidermal keratinocytes.

Dendritic Cells Actively Limit Interleukin-10 Production Under Inflammatory Conditions via DC-SCRIPT and Dual-Specificity Phosphatase 4.

Dendritic cell (DC)-based immunotherapy makes use of the DC's ability to direct the adaptive immune response toward activation or inhibition. DCs perform this immune orchestration in part by secretion of selected cytokines. The most potent anti-inflammatory cytokine interleukin-10 (IL-10) is under tight regulation, as it needs to be predominantly expressed during the resolution phase of the immune response. Currently it is not clear whether there is active suppression of IL-10 by DCs at the initial pro-inflammatory stage of the immune response. Previously, knockdown of the DC-specific transcription factor DC-SCRIPT has been demonstrated to mediate an extensive increase in IL-10 production upon encounter with pro-inflammatory immune stimuli. Here, we explored how DC-SCRIPT contributes to IL-10 suppression under pro-inflammatory conditions by applying chromatin immunoprecipitation sequencing analysis of DC-SCRIPT and the epigenetic marks H3K4me3 and H3K27ac in human DCs. The data showed binding of DC-SCRIPT to a GA-rich motif at H3K27ac-marked genomic enhancers that associated with genes encoding MAPK dual-specificity phosphatases (DUSPs). Functional studies revealed that upon knockdown of DC-SCRIPT, human DCs express much less DUSP4 and exhibit increased phosphorylation of the three major MAPKs (ERK, JNK, and p38). Enhanced ERK signaling in DC-SCRIPT-knockdown-DCs led to higher production of IL-10, which was reverted by rescuing DUSP4 expression. Finally, DC-SCRIPT-knockdown-DCs induced less IFN-γ and increased IL-10 production in naïve T cells, indicative for a more anti-inflammatory phenotype. In conclusion, we have delineated a new mechanism by which DC-SCRIPT allows DCs to limit IL-10 production under inflammatory conditions and potentiate pro-inflammatory Th1 responses. These insights may be exploited to improve DC-based immunotherapies.