ß-Glucan Reverses the Epigenetic State of LPS-Induced Immunological Tolerance.

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, ß-glucan partially reverses the LPS-induced tolerance in vitro. Importantly, ex vivo ß-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.

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease.

Many common variants have been associated with hematological traits, but identification of causal genes and pathways has proven challenging. We performed a genome-wide association analysis in the UK Biobank and INTERVAL studies, testing 29.5 million genetic variants for association with 36 red cell, white cell, and platelet properties in 173,480 European-ancestry participants. This effort yielded hundreds of low frequency (<5%) and rare (<1%) variants with a strong impact on blood cell phenotypes. Our data highlight general properties of the allelic architecture of complex traits, including the proportion of the heritable component of each blood trait explained by the polygenic signal across different genome regulatory domains. Finally, through Mendelian randomization, we provide evidence of shared genetic pathways linking blood cell indices with complex pathologies, including autoimmune diseases, schizophrenia, and coronary heart disease and evidence suggesting previously reported population associations between blood cell indices and cardiovascular disease may be non-causal.

Molecular mechanism of anaerobic ammonium oxidation.

Two distinct microbial processes, denitrification and anaerobic ammonium oxidation (anammox), are responsible for the release of fixed nitrogen as dinitrogen gas (N(2)) to the atmosphere. Denitrification has been studied for over 100 years and its intermediates and enzymes are well known. Even though anammox is a key biogeochemical process of equal importance, its molecular mechanism is unknown, but it was proposed to proceed through hydrazine (N(2)H(4)). Here we show that N(2)H(4) is produced from the anammox substrates ammonium and nitrite and that nitric oxide (NO) is the direct precursor of N(2)H(4). We resolved the genes and proteins central to anammox metabolism and purified the key enzymes that catalyse N(2)H(4) synthesis and its oxidation to N(2). These results present a new biochemical reaction forging an N-N bond and fill a lacuna in our understanding of the biochemical synthesis of the N(2) in the atmosphere. Furthermore, they reinforce the role of nitric oxide in the evolution of the nitrogen cycle.

Nitrite-driven anaerobic methane oxidation by oxygenic bacteria.

Only three biological pathways are known to produce oxygen: photosynthesis, chlorate respiration and the detoxification of reactive oxygen species. Here we present evidence for a fourth pathway, possibly of considerable geochemical and evolutionary importance. The pathway was discovered after metagenomic sequencing of an enrichment culture that couples anaerobic oxidation of methane with the reduction of nitrite to dinitrogen. The complete genome of the dominant bacterium, named 'Candidatus Methylomirabilis oxyfera', was assembled. This apparently anaerobic, denitrifying bacterium encoded, transcribed and expressed the well-established aerobic pathway for methane oxidation, whereas it lacked known genes for dinitrogen production. Subsequent isotopic labelling indicated that 'M. oxyfera' bypassed the denitrification intermediate nitrous oxide by the conversion of two nitric oxide molecules to dinitrogen and oxygen, which was used to oxidize methane. These results extend our understanding of hydrocarbon degradation under anoxic conditions and explain the biochemical mechanism of a poorly understood freshwater methane sink. Because nitrogen oxides were already present on early Earth, our finding opens up the possibility that oxygen was available to microbial metabolism before the evolution of oxygenic photosynthesis.

Streptococcus pneumoniae folate biosynthesis responds to environmental CO2 levels.

Although carbon dioxide (CO2) is known to be essential for Streptococcus pneumoniae growth, it is poorly understood how this respiratory tract pathogen adapts to the large changes in environmental CO2 levels it encounters during transmission, host colonization, and disease. To identify the molecular mechanisms that facilitate pneumococcal growth under CO2-poor conditions, we generated a random S. pneumoniae R6 mariner transposon mutant library representing mutations in 1,538 different genes and exposed it to CO2-poor ambient air. With Tn-seq, we found mutations in two genes that were involved in S. pneumoniae adaptation to changes in CO2 availability. The gene pca, encoding pneumococcal carbonic anhydrase (PCA), was absolutely essential for S. pneumoniae growth under CO2-poor conditions. PCA catalyzes the reversible hydration of endogenous CO2 to bicarbonate (HCO3(-)) and was previously demonstrated to facilitate HCO3(-)-dependent fatty acid biosynthesis. The gene folC that encodes the dihydrofolate/folylpolyglutamate synthase was required at the initial phase of bacterial growth under CO2-poor culture conditions. FolC compensated for the growth-phase-dependent decrease in S. pneumoniae intracellular long-chain (n > 3) polyglutamyl folate levels, which was most pronounced under CO2-poor growth conditions. In conclusion, S. pneumoniae adaptation to changes in CO2 availability involves the retention of endogenous CO2 and the preservation of intracellular long-chain polyglutamyl folate pools.

ChIP-Seq of ERalpha and RNA polymerase II defines genes differentially responding to ligands.

We used ChIP-Seq to map ERalpha-binding sites and to profile changes in RNA polymerase II (RNAPII) occupancy in MCF-7 cells in response to estradiol (E2), tamoxifen or fulvestrant. We identify 10 205 high confidence ERalpha-binding sites in response to E2 of which 68% contain an estrogen response element (ERE) and only 7% contain a FOXA1 motif. Remarkably, 596 genes change significantly in RNAPII occupancy (59% up and 41% down) already after 1 h of E2 exposure. Although promoter proximal enrichment of RNAPII (PPEP) occurs frequently in MCF-7 cells (17%), it is only observed on a minority of E2-regulated genes (4%). Tamoxifen and fulvestrant partially reduce ERalpha DNA binding and prevent RNAPII loading on the promoter and coding body on E2-upregulated genes. Both ligands act differently on E2-downregulated genes: tamoxifen acts as an agonist thus downregulating these genes, whereas fulvestrant antagonizes E2-induced repression and often increases RNAPII occupancy. Furthermore, our data identify genes preferentially regulated by tamoxifen but not by E2 or fulvestrant. Thus (partial) antagonist loaded ERalpha acts mechanistically different on E2-activated and E2-repressed genes.

The metagenome of the marine anammox bacterium 'Candidatus Scalindua profunda' illustrates the versatility of this globally important nitrogen cycle bacterium.

Anaerobic ammonium-oxidizing (anammox) bacteria are responsible for a significant portion of the loss of fixed nitrogen from the oceans, making them important players in the global nitrogen cycle. To date, marine anammox bacteria found in marine water columns and sediments worldwide belong almost exclusively to the 'Candidatus Scalindua' species, but the molecular basis of their metabolism and competitive fitness is presently unknown. We applied community sequencing of a marine anammox enrichment culture dominated by 'Candidatus Scalindua profunda' to construct a genome assembly, which was subsequently used to analyse the most abundant gene transcripts and proteins. In the S. profunda assembly, 4756 genes were annotated, and only about half of them showed the highest identity to the only other anammox bacterium of which a metagenome assembly had been constructed so far, the freshwater 'Candidatus Kuenenia stuttgartiensis'. In total, 2016 genes of S. profunda could not be matched to the K. stuttgartiensis metagenome assembly at all, and a similar number of genes in K.stuttgartiensis could not be found in S. profunda. Most of these genes did not have a known function but 98 expressed genes could be attributed to oligopeptide transport, amino acid metabolism, use of organic acids and electron transport. On the basis of the S. profunda metagenome, and environmental metagenome data, we observed pronounced differences in the gene organization and expression of important anammox enzymes, such as hydrazine synthase (HzsAB), nitrite reductase (NirS) and inorganic nitrogen transport proteins. Adaptations of Scalindua to the substrate limitation of the ocean may include highly expressed ammonium, nitrite and oligopeptide transport systems and pathways for the transport, oxidation, and assimilation of small organic compounds that may allow a more versatile lifestyle contributing to the competitive fitness of Scalindua in the marine realm.

Crosstalk between c-Jun and TAp73alpha/beta contributes to the apoptosis-survival balance.

The p53-family member p73 plays a role in various cellular signaling pathways during development and growth control and it can have tumor suppressor properties. Several isoforms of p73 exist with considerable differences in their function. Whereas the functions of the N-terminal isoforms (TA and ΔNp73) and their opposing pro- and antiapoptotic roles have become evident, the functional differences of the distinct C-terminal splice forms of TAp73 have remained unclear. Here, we characterized the global genomic binding sites for TAp73α and TAp73ß by chromatin immunoprecipitation sequencing as well as the transcriptional responses by performing RNA sequencing. We identified a specific p73 consensus binding motif and found a strong enrichment of AP1 motifs in close proximity to binding sites for TAp73α. These AP1 motif-containing target genes are selectively upregulated by TAp73α, while their mRNA expression is repressed upon TAp73ß induction. We show that their expression is dependent on endogenous c-Jun and that recruitment of c-Jun to the respective AP1 sites was impaired upon TAp73ß expression, in part due to downregulation of c-Jun. Several of these AP1-site containing TAp73α-induced genes impinge on apoptosis induction, suggesting an underlying molecular mechanism for the observed functional differences between TAp73α and TAp73ß.

H2A.Z demarcates intergenic regions of the plasmodium falciparum epigenome that are dynamically marked by H3K9ac and H3K4me3.

Epigenetic regulatory mechanisms and their enzymes are promising targets for malaria therapeutic intervention; however, the epigenetic component of gene expression in P. falciparum is poorly understood. Dynamic or stable association of epigenetic marks with genomic features provides important clues about their function and helps to understand how histone variants/modifications are used for indexing the Plasmodium epigenome. We describe a novel, linear amplification method for next-generation sequencing (NGS) that allows unbiased analysis of the extremely AT-rich Plasmodium genome. We used this method for high resolution, genome-wide analysis of a histone H2A variant, H2A.Z and two histone H3 marks throughout parasite intraerythrocytic development. Unlike in other organisms, H2A.Z is a constant, ubiquitous feature of euchromatic intergenic regions throughout the intraerythrocytic cycle. The almost perfect colocalisation of H2A.Z with H3K9ac and H3K4me3 suggests that these marks are preferentially deposited on H2A.Z-containing nucleosomes. By performing RNA-seq on 8 time-points, we show that acetylation of H3K9 at promoter regions correlates very well with the transcriptional status whereas H3K4me3 appears to have stage-specific regulation, being low at early stages, peaking at trophozoite stage, but does not closely follow changes in gene expression. Our improved NGS library preparation procedure provides a foundation to exploit the malaria epigenome in detail. Furthermore, our findings place H2A.Z at the cradle of P. falciparum epigenetic regulation by stably defining intergenic regions and providing a platform for dynamic assembly of epigenetic and other transcription related complexes.

Dynamic histone H3 epigenome marking during the intraerythrocytic cycle of Plasmodium falciparum.

Epigenome profiling has led to the paradigm that promoters of active genes are decorated with H3K4me3 and H3K9ac marks. To explore the epigenome of Plasmodium falciparum asexual stages, we performed MS analysis of histone modifications and found a general preponderance of H3/H4 acetylation and H3K4me3. ChIP-on-chip profiling of H3, H3K4me3, H3K9me3, and H3K9ac from asynchronous parasites revealed an extensively euchromatic epigenome with heterochromatin restricted to variant surface antigen gene families (VSA) and a number of genes hitherto unlinked to VSA. Remarkably, the vast majority of the genome shows an unexpected pattern of enrichment of H3K4me3 and H3K9ac. Analysis of synchronized parasites revealed significant developmental stage specificity of the epigenome. In rings, H3K4me3 and H3K9ac are homogenous across the genes marking active and inactive genes equally, whereas in schizonts, they are enriched at the 5' end of active genes. This study reveals an unforeseen and unique plasticity in the use of the epigenetic marks and implies the presence of distinct epigenetic pathways in gene silencing/activation throughout the erythrocytic cycle.

Single cell transcriptomic analysis of the immune cell compartment in the human small intestine and in Celiac disease.

Celiac disease is an autoimmune disorder in which ingestion of dietary gluten triggers an immune reaction in the small intestine leading to destruction of the lining epithelium. Current treatment focusses on lifelong adherence to a gluten-free diet. Gluten-specific CD4+ T cells and cytotoxic intraepithelial CD8+ T cells have been proposed to be central in disease pathogenesis. Here we use unbiased single-cell RNA-sequencing and explore the heterogeneity of CD45+ immune cells in the human small intestine. We show altered myeloid cell transcriptomes present in active celiac lesions. CD4+ and CD8+ T cells transcriptomes show extensive changes and we define a natural intraepithelial lymphocyte population that is reduced in celiac disease. We show that the immune landscape in Celiac patients on a gluten-free diet is only partially restored compared to control samples. Altogether, we provide a single cell transcriptomic resource that can inform the immune landscape of the small intestine during Celiac disease.

The novel p53 target gene IRF2BP2 participates in cell survival during the p53 stress response.

The tumor suppressor p53 contributes to the cellular fate after genotoxic insults, mainly through the regulation of target genes, thereby allowing e.g. repair mechanisms resulting in cell survival or inducing apoptosis. Unresolved so far is the issue, which exact mechanisms lead to one or the other cellular outcome. Here, we describe the interferon regulatory factor-2-binding protein-2 (IRF2BP2) as a new direct target gene of p53, influencing the p53-mediated cellular decision. We show that upregulation of IRF2BP2 after treatment with actinomycin D (Act.D) is dependent on functional p53 in different cell lines. This occurs in parallel with the down-regulation of the interacting partner of IRF2BP2, the interferon regulatory factor-2 (IRF2), which is known to positively influence cell growth. Analyzing the molecular functions of IRF2BP2, it appears to be able to impede on the p53-mediated transactivation of the p21- and the Bax-gene. We show here that overexpressed IRF2BP2 has an impact on the cellular stress response after Act.D treatment and that it diminishes the induction of apoptosis after doxorubicin treatment. Furthermore, the knockdown of IRF2BP2 leads to an upregulation of p21 and faster induction of apoptosis after doxorubicin as well as Act.D treatment.

Blueprint of human thymopoiesis reveals molecular mechanisms of stage-specific TCR enhancer activation.

Cell differentiation is accompanied by epigenetic changes leading to precise lineage definition and cell identity. Here we present a comprehensive resource of epigenomic data of human T cell precursors along with an integrative analysis of other hematopoietic populations. Although T cell commitment is accompanied by large scale epigenetic changes, we observed that the majority of distal regulatory elements are constitutively unmethylated throughout T cell differentiation, irrespective of their activation status. Among these, the TCRA gene enhancer (Eα) is in an open and unmethylated chromatin structure well before activation. Integrative analyses revealed that the HOXA5-9 transcription factors repress the Eα enhancer at early stages of T cell differentiation, while their decommission is required for TCRA locus activation and enforced αß T lineage differentiation. Remarkably, the HOXA-mediated repression of Eα is paralleled by the ectopic expression of homeodomain-related oncogenes in T cell acute lymphoblastic leukemia. These results highlight an analogous enhancer repression mechanism at play in normal and cancer conditions, but imposing distinct developmental constraints.

The acute myeloid leukemia associated AML1-ETO fusion protein alters the transcriptome and cellular progression in a single-oncogene expressing in vitro induced pluripotent stem cell based granulocyte differentiation model.

Acute myeloid leukemia (AML) is characterized by recurrent mutations that affect normal hematopoiesis. The analysis of human AMLs has mostly been performed using end-point materials, such as cell lines and patient derived AMLs that also carry additional contributing mutations. The molecular effects of a single oncogenic hit, such as expression of the AML associated oncoprotein AML1-ETO on hematopoietic development and transformation into a (pre-) leukemic state still needs further investigation. Here we describe the development and characterization of an induced pluripotent stem cell (iPSC) system that allows in vitro differentiation towards different mature myeloid cell types such as monocytes and granulocytes. During in vitro differentiation we expressed the AML1-ETO fusion protein and examined the effects of the oncoprotein on differentiation and the underlying alterations in the gene program at 8 different time points. Our analysis revealed that AML1-ETO as a single oncogenic hit in a non-mutated background blocks granulocytic differentiation, deregulates the gene program via altering the acetylome of the differentiating granulocytic cells, and induces t(8;21) AML associated leukemic characteristics. Together, these results reveal that inducible oncogene expression during in vitro differentiation of iPS cells provides a valuable platform for analysis of aberrant regulation in disease.

Chromatin-Based Classification of Genetically Heterogeneous AMLs into Two Distinct Subtypes with Diverse Stemness Phenotypes.

Global investigation of histone marks in acute myeloid leukemia (AML) remains limited. Analyses of 38 AML samples through integrated transcriptional and chromatin mark analysis exposes 2 major subtypes. One subtype is dominated by patients with NPM1 mutations or MLL-fusion genes, shows activation of the regulatory pathways involving HOX-family genes as targets, and displays high self-renewal capacity and stemness. The second subtype is enriched for RUNX1 or spliceosome mutations, suggesting potential interplay between the 2 aberrations, and mainly depends on IRF family regulators. Cellular consequences in prognosis predict a relatively worse outcome for the first subtype. Our integrated profiling establishes a rich resource to probe AML subtypes on the basis of expression and chromatin data.

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.

Multi-omics profiling reveals a distinctive epigenome signature for high-risk acute promyelocytic leukemia.

Epigenomic alterations have been associated with both pathogenesis and progression of cancer. Here, we analyzed the epigenome of two high-risk APL (hrAPL) patients and compared it to non-high-risk APL cases. Despite the lack of common genetic signatures, we found that human hrAPL blasts from patients with extremely poor prognosis display specific patterns of histone H3 acetylation, specifically hyperacetylation at a common set of enhancer regions. In addition, unique profiles of the repressive marks H3K27me3 and DNA methylation were exposed in high-risk APLs. Epigenetic comparison with low/intermediate-risk APLs and AMLs revealed hrAPL-specific patterns of histone acetylation and DNA methylation, suggesting these could be further developed into markers for clinical identification. The epigenetic drug MC2884, a newly generated general HAT/EZH2 inhibitor, induces apoptosis of high-risk APL blasts and reshapes their epigenomes by targeting both active and repressive marks. Together, our analysis uncovers distinctive epigenome signatures of hrAPL patients, and provides proof of concept for use of epigenome profiling coupled to epigenetic drugs to 'personalize' precision medicine.

Dynamics of Transcription Regulation in Human Bone Marrow Myeloid Differentiation to Mature Blood Neutrophils.

Neutrophils are short-lived blood cells that play a critical role in host defense against infections. To better comprehend neutrophil functions and their regulation, we provide a complete epigenetic overview, assessing important functional features of their differentiation stages from bone marrow-residing progenitors to mature circulating cells. Integration of chromatin modifications, methylation, and transcriptome dynamics reveals an enforced regulation of differentiation, for cellular functions such as release of proteases, respiratory burst, cell cycle regulation, and apoptosis. We observe an early establishment of the cytotoxic capability, while the signaling components that activate these antimicrobial mechanisms are transcribed at later stages, outside the bone marrow, thus preventing toxic effects in the bone marrow niche. Altogether, these data reveal how the developmental dynamics of the chromatin landscape orchestrate the daily production of a large number of neutrophils required for innate host defense and provide a comprehensive overview of differentiating human neutrophils.

Combined HAT/EZH2 modulation leads to cancer-selective cell death.

Epigenetic alterations have been associated with both pathogenesis and progression of cancer. By screening of library compounds, we identified a novel hybrid epi-drug MC2884, a HAT/EZH2 inhibitor, able to induce bona fide cancer-selective cell death in both solid and hematological cancers in vitro, ex vivo and in vivo xenograft models. Anticancer action was due to an epigenome modulation by H3K27me3, H3K27ac, H3K9/14ac decrease, and to caspase-dependent apoptosis induction. MC2884 triggered mitochondrial pathway apoptosis by up-regulation of cleaved-BID, and strong down-regulation of BCL2. Even aggressive models of cancer, such as p53-/- or TET2-/- cells, responded to MC2884, suggesting MC2884 therapeutic potential also for the therapy of TP53 or TET2-deficient human cancers. MC2884 induced massive apoptosis in ex vivo human primary leukemia blasts with poor prognosis in vivo, by targeting BCL2 expression. MC2884-treatment reduced acetylation of the BCL2 promoter at higher level than combined p300 and EZH2 inhibition. This suggests a key role for BCL-2 reduction in potentiating responsiveness, also in combination therapy with BCL2 inhibitors. Finally, we identified both the mechanism of MC2884 action as well as a potential therapeutic scheme of its use. Altogether, this provides proof of concept for the use of epi-drugs coupled with epigenome analyses to 'personalize' precision medicine.