Search

Systemic inflammation (SI) plays a detrimental role in various conditions with high mortality rates1–4. SI manifests an acute hyperinflammation followed by long-lasting immunosuppression, increasing patients’ risks for secondary infections and impaired clinical outcomes5–7. Due to the extensive heterogeneity in SI etiology, the mechanisms governing these states are incompletely understood. Here, we characterized acute and late effects of lipopolysaccharide (LPS)-induced SI (LPS-SI8) on blood monocytes and bone marrow (BM) cells of healthy volunteers. Like clinical SI, LPS administration elicited a profound but transient acute response. Single-cell transcriptomic analysis of acute LPS-SI unveiled loss of BM monocytes and appearance of an inflammatory monocyte-like (i-Mono’s) population, expressing gene programs similar to early-stage sepsis patients9. In the ensuing late phase of LPS-SI, we observed reduced expression of interferon type I (IFN-I) responsive genes in monocytes and profound attenuation of in vivo response to a second LPS challenge. Furthermore, late LPS-SI led to impaired myelopoiesis with a loss of intermediate and non-classical monocytes. In accordance, we show compromised myelopoiesis also occurs in late-stage sepsis. Finally, IFNβ treatment reversed LPS-induced immunosuppression in monocytes. Our results reveal long-lasting effects of SI on myelopoiesis and substantiate the importance of IFN-I in the pathophysiology of SI-induced immunosuppression.

Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors of varied genetic makeup, associated with high cardiovascular morbidity and a variable risk of malignancy. The source of the transcriptional heterogeneity of the disease and the underlying biological processes determining the outcome in PCPG remains largely unclear. We focused on PCPG tumors with germline SDHB and RET mutations, representing distinct prognostic groups with worse or better prognoses, respectively. We applied single-nuclei RNA sequencing (snRNA-seq) on tissue samples from 11 patients and found high patient-to-patient transcriptome heterogeneity of the neuroendocrine tumor cells. The tumor microenvironment also showed heterogeneous profiles mainly contributed by macrophages of the immune cell clusters and Schwann cells of the stroma. Performing non-negative matrix factorization we identified common transcriptional programs active in RET and SDHB as well as distinct modules? including neuronal development, hormone synthesis and secretion, and DNA replication. Comparison of the SDHB and RET transcriptomes with that of developmental stages of adrenal gland formation suggests different developmental stages at which PC and PG tumors appear to be arrested.

Macrophage (Mϕ) repolarization from a pro-tumor, immunosuppressive phenotype towards an anti-tumor, pro-inflammatory state represents a promising therapeutic strategy in patients with cancer1. Successful reprogramming of Mϕ in a clinical setting has not been documented. Here, we traced the evolution at single-cell resolution of a diffuse midline glioma (DMG) with H3K27M mutation which metastasized in the abdomen after placing a ventriculoperitoneal (VP) shunt and exploited this information for therapeutic decision-making. The primary tumor showed a complex cellular and genomic landscape characterized by heterogeneous cancer cells and a tumor-supportive immune microenvironment. Upon metastasis, malignant cells populated the peritoneum and triggered a massive anti-inflammatory immune cell infiltration and expansion of myeloid-derived suppressor cells (MDSCs) in peripheral blood. Hydroxychloroquine adjuvant treatment started to overcome the immunosuppressive milieu, which resulted in a decrease in peritoneal cancer cells and pro-tumor innate immune cells. Importantly, an emergence of anti-tumor, pro-inflammatory macrophages and cytotoxic T-cells was observed, accompanied by the activation of monocytes in the blood. Our study advocates the employment of single-cell technologies to better understand and inspire therapeutic regimens in patients with cancer.

Significance During infection, extracellular “labile” heme, released from damaged red blood or parenchymal cells, acts as prototypical alarmin stimulating myeloid cells. A characteristic hallmark of myeloid cell activation is the development of trained immunity, specified as long-lasting adaptations based on transcriptional and epigenetic modifications. In vivo, this is maintained by the rerouting of hematopoiesis. We found that heme is a previously unrecognized trained immunity inducer promoting resistance to bacterial infection in mice. This goes along with extensive long-lasting epigenetic memory in hematopoietic stem cells provoking drastic changes in the transcription factor–binding landscape of myeloid progenitor cells. Given the critical role of heme during infections, we propose that trained immunity is a more general component of innate immunity than previously suggested., Trained immunity defines long-lasting adaptations of innate immunity based on transcriptional and epigenetic modifications of myeloid cells and their bone marrow progenitors [M. Divangahi et al., Nat. Immunol. 22, 2–6 (2021)]. Innate immune cells, however, do not exclusively differentiate between foreign and self but also react to host-derived molecules referred to as alarmins. Extracellular “labile” heme, released during infections, is a bona fide alarmin promoting myeloid cell activation [M. P. Soares, M. T. Bozza, Curr. Opin. Immunol. 38, 94–100 (2016)]. Here, we report that labile heme is a previously unrecognized inducer of trained immunity that confers long-term regulation of lineage specification of hematopoietic stem cells and progenitor cells. In contrast to previous reports on trained immunity, essentially mediated by pathogen-associated molecular patterns, heme training depends on spleen tyrosine kinase signal transduction pathway acting upstream of c-Jun N-terminal kinases. Heme training promotes resistance to sepsis, is associated with the expansion of self-renewing hematopoetic stem cells primed toward myelopoiesis and to the occurrence of a specific myeloid cell population. This is potentially evoked by sustained activity of Nfix, Runx1, and Nfe2l2 and dissociation of the transcriptional repressor Bach2. Previously reported trained immunity inducers are, however, infrequently present in the host, whereas heme abundantly occurs during noninfectious and infectious disease. This difference might explain the vanishing protection exerted by heme training in sepsis over time with sustained long-term myeloid adaptations. Hence, we propose that trained immunity is an integral component of innate immunity with distinct functional differences on infectious disease outcome depending on its induction by pathogenic or endogenous molecules.

Clusters of enhancers, referred as to super-enhancers (SEs), control the expression of cell identity genes. The organisation of these clusters, and how they are remodelled upon developmental transitions remain poorly understood. Here, we report the existence of two types of enhancer units within SEs typified by distinctive CpG methylation dynamics in embryonic stem cells (ESCs). We find that these units are either prone for decommissioning or remain constitutively active in epiblast stem cells (EpiSCs), as further established in the peri-implantation epiblast in vivo. Mechanistically, we show a pivotal role for ESRRB in regulating the activity of ESC-specific enhancer units and propose that the developmentally regulated silencing of ESRRB triggers the selective inactivation of these units within SEs. Our study provides insights into the molecular events that follow the loss of ESRRB binding, and offers a mechanism by which the naive pluripotency transcriptional programme can be partially reset upon embryo implantation., Clusters of enhancers, called super-enhancers (SEs), control the expression of cell identity genes. Here, the authors identify two types of enhancer units within SEs in ESCs, characterised by distinctive CpG methylation dynamics, and find that ESRRB regulates the enhancer units decommissioned at exit from naive pluripotency.

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.

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.

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.

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.

Contains fulltext : 171940.pdf (Publisher’s version ) (Open Access) We have developed the Weighted Gene Expression Tool and database (WeGET, http://weget.cmbi.umcn.nl) for the prediction of new genes of a molecular system by correlated gene expression. WeGET utilizes a compendium of 465 human and 560 murine gene expression datasets that have been collected from multiple tissues under a wide range of experimental conditions. It exploits this abundance of expression data by assigning a high weight to datasets in which the known genes of a molecular system are harmoniously up- and down-regulated. WeGET ranks new candidate genes by calculating their weighted co-expression with that system. A weighted rank is calculated for human genes and their mouse orthologs. Then, an integrated gene rank and p-value is computed using a rank-order statistic. We applied our method to predict novel genes that have a high degree of co-expression with Gene Ontology terms and pathways from KEGG and Reactome. For each query set we provide a list of predicted novel genes, computed weights for transcription datasets used and cell and tissue types that contributed to the final predictions. The performance for each query set is assessed by 10-fold cross-validation. Finally, users can use the WeGET to predict novel genes that co-express with a custom query set.

Contains fulltext : 165879.pdf (Publisher’s version ) (Open Access) Innate immune memory is the phenomenon whereby innate immune cells such as monocytes or macrophages undergo functional reprogramming after exposure to microbial components such as lipopolysaccharide (LPS). We apply an integrated epigenomic approach to characterize the molecular events involved in LPS-induced tolerance in a time-dependent manner. Mechanistically, LPS-treated monocytes fail to accumulate active histone marks at promoter and enhancers of genes in the lipid metabolism and phagocytic pathways. Transcriptional inactivity in response to a second LPS exposure in tolerized macrophages is accompanied by failure to deposit active histone marks at promoters of tolerized genes. In contrast, beta-glucan partially reverses the LPS-induced tolerance in vitro. Importantly, ex vivo beta-glucan treatment of monocytes from volunteers with experimental endotoxemia re-instates their capacity for cytokine production. Tolerance is reversed at the level of distal element histone modification and transcriptional reactivation of otherwise unresponsive genes. VIDEO ABSTRACT.

Contains fulltext : 152194.PDF (Publisher’s version ) (Open Access) MOTIVATION: Genome-scale metabolic networks can be modeled in a constraint-based fashion. Reaction stoichiometry combined with flux capacity constraints determine the space of allowable reaction rates. This space is often large and a central challenge in metabolic modeling is finding the biologically most relevant flux distributions. A widely used method is flux balance analysis (FBA), which optimizes a biologically relevant objective such as growth or ATP production. Although FBA has proven to be highly useful for predicting growth and byproduct secretion, it cannot predict the intracellular fluxes under all environmental conditions. Therefore, alternative strategies have been developed to select flux distributions that are in agreement with experimental "omics" data, or by incorporating experimental flux measurements. The latter, unfortunately can only be applied to a limited set of reactions and is currently not feasible at the genome-scale. On the other hand, it has been observed that micro-organisms favor a suboptimal growth rate, possibly in exchange for a more "flexible" metabolic network. Instead of dedicating the internal network state to an optimal growth rate in one condition, a suboptimal growth rate is used, that allows for an easier switch to other nutrient sources. A small decrease in growth rate is exchanged for a relatively large gain in metabolic capability to adapt to changing environmental conditions. RESULTS: Here, we propose Maximum Metabolic Flexibility (MMF) a computational method that utilizes this observation to find the most probable intracellular flux distributions. By mapping measured flux data from central metabolism to the genome-scale models of Escherichia coli and Saccharomyces cerevisiae we show that i) indeed, most of the measured fluxes agree with a high adaptability of the network, ii) this result can be used to further reduce the space of feasible solutions iii) this reduced space improves the quantitative predictions made by FBA and contains a significantly larger fraction of the measured fluxes compared to the flux space that was reduced by a uniform sampling approach and iv) MMF can be used to select reactions in the network that contribute most to the steady-state flux space. Constraining the selected reactions improves the quantitative predictions of FBA considerably more than adding an equal amount of flux constraints, selected using a more naive approach. Our method can be applied to any cell type without requiring prior information. AVAILABILITY: MMF is freely available as a MATLAB plugin at: http://cs.ru.nl/~wmegchel/mmf.

Synthetic Lethal (SL) genetic interactions play a key role in various types of biological research, ranging from understanding genotype-phenotype relationships to identifying drug-targets against cancer. Despite recent advances in empirical measuring SL interactions in human cells, the human genetic interaction map is far from complete. Here, we present a novel approach to predict this map by exploiting patterns in cancer genome evolution. First, we show that empirically determined SL interactions are reflected in various gene presence, absence, and duplication patterns in hundreds of cancer genomes. The most evident pattern that we discovered is that when one member of an SL interaction gene pair is lost, the other gene tends not to be lost, i.e. the absence of co-loss. This observation is in line with expectation, because the loss of an SL interacting pair will be lethal to the cancer cell. SL interactions are also reflected in gene expression profiles, such as an under representation of cases where the genes in an SL pair are both under expressed, and an over representation of cases where one gene of an SL pair is under expressed, while the other one is over expressed. We integrated the various previously unknown cancer genome patterns and the gene expression patterns into a computational model to identify SL pairs. This simple, genome-wide model achieves a high prediction power (AUC = 0.75) for known genetic interactions. It allows us to present for the first time a comprehensive genome-wide list of SL interactions with a high estimated prediction precision, covering up to 591,000 gene pairs. This unique list can potentially be used in various application areas ranging from biotechnology to medical genetics.

Contains fulltext : 149260.pdf (Publisher’s version ) (Closed access) Synthetic dosage lethality (SDL) denotes a genetic interaction between two genes whereby the underexpression of gene A combined with the overexpression of gene B is lethal. SDLs offer a promising way to kill cancer cells by inhibiting the activity of SDL partners of activated oncogenes in tumors, which are often difficult to target directly. As experimental genome-wide SDL screens are still scarce, here we introduce a network-level computational modeling framework that quantitatively predicts human SDLs in metabolism. For each enzyme pair (A, B) we systematically knock out the flux through A combined with a stepwise flux increase through B and search for pairs that reduce cellular growth more than when either enzyme is perturbed individually. The predictive signal of the emerging network of 12,000 SDLs is demonstrated in five different ways. (i) It can be successfully used to predict gene essentiality in shRNA cancer cell line screens. Moving to clinical tumors, we show that (ii) SDLs are significantly underrepresented in tumors. Furthermore, breast cancer tumors with SDLs active (iii) have smaller sizes and (iv) result in increased patient survival, indicating that activation of SDLs increases cancer vulnerability. Finally, (v) patient survival improves when multiple SDLs are present, pointing to a cumulative effect. This study lays the basis for quantitative identification of cancer SDLs in a model-based mechanistic manner. The approach presented can be used to identify SDLs in species and cell types in which "omics" data necessary for data-driven identification are missing.