Search

Chronic hepatitis B virus (HBV) infection affects 300 million patients worldwide1,2, in whom virus-specific CD8 T cells by still ill-defined mechanisms lose their function and cannot eliminate HBV-infected hepatocytes3–7. Here we demonstrate that a liver immune rheostat renders virus-specific CD8 T cells refractory to activation and leads to their loss of effector functions. In preclinical models of persistent infection with hepatotropic viruses such as HBV, dysfunctional virus-specific CXCR6+CD8 T cells accumulated in the liver and, as a characteristic hallmark, showed enhanced transcriptional activity of cAMP-responsive element modulator (CREM) distinct from T cell exhaustion. In patients with chronic hepatitis B, circulating and intrahepatic HBV-specific CXCR6+CD8 T cells with enhanced CREMexpression and transcriptional activity were detected at a frequency of 12–22% of HBV-specific CD8 T cells. Knocking out the inhibitory CREM/ICERisoform in T cells, however, failed to rescue T cell immunity. This indicates that CREM activity was a consequence, rather than the cause, of loss in T cell function, further supported by the observation of enhanced phosphorylation of protein kinase A (PKA) which is upstream of CREM. Indeed, we found that enhanced cAMP–PKA-signalling from increased T cell adenylyl cyclase activity augmented CREM activity and curbed T cell activation and effector function in persistent hepatic infection. Mechanistically, CD8 T cells recognizing their antigen on hepatocytes established close and extensive contact with liver sinusoidal endothelial cells, thereby enhancing adenylyl cyclase–cAMP–PKA signalling in T cells. In these hepatic CD8 T cells, which recognize their antigen on hepatocytes, phosphorylation of key signalling kinases of the T cell receptor signalling pathway was impaired, which rendered them refractory to activation. Thus, close contact with liver sinusoidal endothelial cells curbs the activation and effector function of HBV-specific CD8 T cells that target hepatocytes expressing viral antigens by means of the adenylyl cyclase–cAMP–PKA axis in an immune rheostat-like fashion.

Deregulation of the EVI1proto-oncogene by the GATA2distal hematopoietic enhancer (G2DHE) is a key event in high-risk acute myeloid leukemia carrying 3q21q26 aberrations (3q-AML). Upon chromosomal rearrangement, G2DHEacquires characteristics of a super-enhancer and causes overexpression of EVI1at 3q26.2. However, the transcription factor (TF) complex of G2DHEremains poorly characterized. The aim of this study was to unravel key components of G2DHE-bound TFs involved in the deregulation of EVI1. We have identified several CEBPA and RUNX1 binding sites to be enriched and critical for G2DHEfunction in 3q-AML cells. Using ChIP-SICAP (ChIP followed by selective isolation of chromatin-associated proteins), a panel of chromatin interactors of RUNX1 and CEBPA were detected in 3q-AML, including PARP1 and IKZF1. PARP1 inhibition (PARPi) caused a reduction of EVI1expression and a decrease in EVI1–G2DHEinteraction frequency, highlighting the involvement of PARP1 in oncogenic super-enhancer formation. Furthermore, 3q-AML cells were highly sensitive to PARPi and displayed morphological changes with higher rates of differentiation and apoptosis as well as depletion of CD34 + cells. In summary, integrative analysis of the 3q-AML super-enhancer complex identified CEBPA and RUNX1 associated proteins and nominated PARP1 as a potential new therapeutic target in EVI1+ 3q-AML.

Neuroblastoma is a pediatric tumor of the developing sympathetic nervous system. However, the cellular origin of neuroblastoma has yet to be defined. Here we studied the single-cell transcriptomes of neuroblastomas and normal human developing adrenal glands at various stages of embryonic and fetal development. We defined normal differentiation trajectories from Schwann cell precursors over intermediate states to neuroblasts or chromaffin cells and showed that neuroblastomas transcriptionally resemble normal fetal adrenal neuroblasts. Importantly, neuroblastomas with varying clinical phenotypes matched different temporal states along normal neuroblast differentiation trajectories, with the degree of differentiation corresponding to clinical prognosis. Our work highlights the roles of oncogenic MYCN and loss of TFAP2B in blocking differentiation and may provide the basis for designing therapeutic interventions to overcome differentiation blocks.

Acute myeloid leukemia is characterized by the accumulation of clonal myeloid blast cells unable to differentiate into mature leukocytes. Chemotherapy induces remission in the majority of patients, but relapse rates are high and lead to poor clinical outcomes. Because this is primarily caused by chemotherapy-resistant leukemic stem cells (LSCs), it is essential to eradicate LSCs to improve patient survival. LSCs have predominantly been studied at the transcript level, thus information about posttranscriptionally regulated genes and associated networks is lacking. Here, we extend our previous report on LSC proteomes to healthy age-matched hematopoietic stem and progenitor cells (HSPCs) and correlate the proteomes to the corresponding transcriptomes. By comparing LSCs to leukemic blasts and healthy HSPCs, we validate candidate LSC markers and highlight novel and potentially targetable proteins that are absent or only lowly expressed in HSPCs. In addition, our data provide strong evidence that LSCs harbor a characteristic energy metabolism, adhesion molecule composition, as well as RNA-processing properties. Furthermore, correlating proteome and transcript data of the same individual samples highlights the strength of proteome analyses, which are particularly potent in detecting alterations in metabolic pathways. In summary, our study provides a comprehensive proteomic and transcriptomic characterization of functionally validated LSCs, blasts, and healthy HSPCs, representing a valuable resource helping to design LSC-directed therapies.

Acute myeloid leukemia is characterized by the accumulation of clonal myeloid blast cells unable to differentiate into mature leukocytes. Chemotherapy induces remission in the majority of patients, but relapse rates are high and lead to poor clinical outcomes. Because this is primarily caused by chemotherapy-resistant leukemic stem cells (LSCs), it is essential to eradicate LSCs to improve patient survival. LSCs have predominantly been studied at the transcript level, thus information about posttranscriptionally regulated genes and associated networks is lacking. Here, we extend our previous report on LSC proteomes to healthy age-matched hematopoietic stem and progenitor cells (HSPCs) and correlate the proteomes to the corresponding transcriptomes. By comparing LSCs to leukemic blasts and healthy HSPCs, we validate candidate LSC markers and highlight novel and potentially targetable proteins that are absent or only lowly expressed in HSPCs. In addition, our data provide strong evidence that LSCs harbor a characteristic energy metabolism, adhesion molecule composition, as well as RNA-processing properties. Furthermore, correlating proteome and transcript data of the same individual samples highlights the strength of proteome analyses, which are particularly potent in detecting alterations in metabolic pathways. In summary, our study provides a comprehensive proteomic and transcriptomic characterization of functionally validated LSCs, blasts, and healthy HSPCs, representing a valuable resource helping to design LSC-directed therapies.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1–4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5′ region of TP53, in the 3′ untranslated regions of NFKBIZand TOB1, focal deletions in BRD4and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

The branched-chain amino acid (BCAA) pathway and high levels of BCAA transaminase 1 (BCAT1) have recently been associated with aggressiveness in several cancer entities. However, the mechanistic role of BCAT1 in this process remains largely uncertain. Here, by performing high-resolution proteomic analysis of human acute myeloid leukaemia (AML) stem-cell and non-stem-cell populations, we find the BCAA pathway enriched and BCAT1 protein and transcripts overexpressed in leukaemia stem cells. We show that BCAT1, which transfers α-amino groups from BCAAs to α-ketoglutarate (αKG), is a critical regulator of intracellular αKG homeostasis. Further to its role in the tricarboxylic acid cycle, αKG is an essential cofactor for αKG-dependent dioxygenases such as Egl-9 family hypoxia inducible factor 1 (EGLN1) and the ten-eleven translocation (TET) family of DNA demethylases. Knockdown of BCAT1 in leukaemia cells caused accumulation of αKG, leading to EGLN1-mediated HIF1α protein degradation. This resulted in a growth and survival defect and abrogated leukaemia-initiating potential. By contrast, overexpression of BCAT1 in leukaemia cells decreased intracellular αKG levels and caused DNA hypermethylation through altered TET activity. AML with high levels of BCAT1 (BCAT1high) displayed a DNA hypermethylation phenotype similar to cases carrying a mutant isocitrate dehydrogenase (IDHmut), in which TET2 is inhibited by the oncometabolite 2-hydroxyglutarate. High levels of BCAT1 strongly correlate with shorter overall survival in IDHWTTET2WT, but not IDHmutor TET2mutAML. Gene sets characteristic for IDHmutAML were enriched in samples from patients with an IDHWTTET2WTBCAT1highstatus. BCAT1highAML showed robust enrichment for leukaemia stem-cell signatures, and paired sample analysis showed a significant increase in BCAT1 levels upon disease relapse. In summary, by limiting intracellular αKG, BCAT1 links BCAA catabolism to HIF1α stability and regulation of the epigenomic landscape, mimicking the effects of IDH mutations. Our results suggest the BCAA–BCAT1–αKG pathway as a therapeutic target to compromise leukaemia stem-cell function in patients with IDHWTTET2WTAML.

Regulatory T cells (Tregcells) perform two distinct functions: they maintain self-tolerance, and they support organ homeostasis by differentiating into specialized tissue Tregcells. We found that epigenetic modifications defined the molecular characteristics of tissue Tregcells. Tagmentation-based whole-genome bisulfite sequencing revealed more than 11,000 regions that were methylated differentially in pairwise comparisons of tissue Tregcell populations and lymphoid T cells. Similarities in the epigenetic landscape led to the identification of a common tissue Tregcell population that was present in many organs and was characterized by gain and loss of DNA methylation that included many gene sites associated with the TH2 subset of helper T cells, such as the gene encoding cytokine IL-33 receptor ST2, as well as the production of tissue-regenerative factors. Furthermore, the ST2-expressing population was dependent on the transcriptional regulator BATF and could be expanded by IL-33. Thus, tissue Tregcells integrate multiple waves of epigenetic reprogramming that define their tissue-restricted specialization.

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+T cells, highlighting the importance of host genome organization in HIV-1 infection.

Epigenetic alterations during cellular differentiation are a key molecular mechanism which both instructs and reinforces the process of lineage commitment. Within the haematopoietic system, progressive changes in the DNA methylome of haematopoietic stem cells (HSCs) are essential for the effective production of mature blood cells. Inhibition or loss of function of the cellular DNA methylation machinery has been shown to lead to a severe perturbation in blood production and is also an important driver of malignant transformation. HSCs constitute a very rare cell population in the bone marrow, capable of life-long self-renewal and multi-lineage differentiation. The low abundance of HSCs has been a major technological barrier to the global analysis of the CpG methylation status within both HSCs and their immediate progeny, the multipotent progenitors (MPPs). Within this Extra View article, we review the current understanding of how the DNA methylome regulates normal and malignant hematopoiesis. We also discuss the current methodologies that are available for interrogating the DNA methylation status of HSCs and MPPs and describe a new data set that was generated using tagmentation-based whole genome bisulfite sequencing (TWGBS) in order to comprehensively map methylated cytosines using the limited amount of genomic DNA that can be harvested from rare cell populations. Extended analysis of this data set clearly demonstrates the added value of genome-wide sequencing of methylated cytosines and identifies novel important cis-acting regulatory regions that are dynamically remodeled during the first steps of haematopoietic differentiation.

In this study, we present integrated quantitative proteome, transcriptome, and methylome analyses of hematopoietic stem cells (HSCs) and four multipotent progenitor (MPP) populations. From the characterization of more than 6,000 proteins, 27,000 transcripts, and 15,000 differentially methylated regions (DMRs), we identified coordinated changes associated with early differentiation steps. DMRs show continuous gain or loss of methylation during differentiation, and the overall change in DNA methylation correlates inversely with gene expression at key loci. Our data reveal the differential expression landscape of 493 transcription factors and 682 lncRNAs and highlight specific expression clusters operating in HSCs. We also found an unexpectedly dynamic pattern of transcript isoform regulation, suggesting a critical regulatory role during HSC differentiation, and a cell cycle/DNA repair signature associated with multipotency in MPP2 cells. This study provides a comprehensive genome-wide resource for the functional exploration of molecular, cellular, and epigenetic regulation at the top of the hematopoietic hierarchy.

Myelodysplastic syndromes (MDSs) are a heterogeneous group of myeloid neoplasms with defects in hematopoietic stem and progenitor cells (HSPCs) and possibly the HSPC niche. Here, we show that patient-derived mesenchymal stromal cells (MDS MSCs) display a disturbed differentiation program and are essential for the propagation of MDS-initiating Lin−CD34+CD38−stem cells in orthotopic xenografts. Overproduction of niche factors such as CDH2 (N-Cadherin), IGFBP2, VEGFA, and LIF is associated with the ability of MDS MSCs to enhance MDS expansion. These factors represent putative therapeutic targets in order to disrupt critical hematopoietic-stromal interactions in MDS. Finally, healthy MSCs adopt MDS MSC-like molecular features when exposed to hematopoietic MDS cells, indicative of an instructive remodeling of the microenvironment. Therefore, this patient-derived xenograft model provides functional and molecular evidence that MDS is a complex disease that involves both the hematopoietic and stromal compartments. The resulting deregulated expression of niche factors may well also be a feature of other hematopoietic malignancies.

This protocol explains how to use the online integrated pipeline 'peak-motifs' (http://rsat.ulb.ac.be/rsat/) to predict motifs and binding sites in full-size peak sets obtained by chromatin immunoprecipitation–sequencing (ChIP-seq) or related technologies. The workflow combines four time- and memory-efficient motif discovery algorithms to extract significant motifs from the sequences. Discovered motifs are compared with databases of known motifs to identify potentially bound transcription factors. Sequences are scanned to predict transcription factor binding sites and analyze their enrichment and positional distribution relative to peak centers. Peaks and binding sites are exported as BED tracks that can be uploaded into the University of California Santa Cruz (UCSC) genome browser for visualization in the genomic context. This protocol is illustrated with the analysis of a set of 6,000 peaks (8 Mb in total) bound by the Drosophilatranscription factor Krüppel. The complete workflow is achieved in about 25 min of computational time on the Regulatory Sequence Analysis Tools (RSAT) Web server. This protocol can be followed in about 1 h.

We report the direct reprogramming of both adult human fibroblasts and blood cells into induced neural plate border stem cells (iNBSCs) by ectopic expression of four neural transcription factors. Self-renewing, clonal iNBSCs can be robustly expanded in defined media while retaining multilineage differentiation potential. They generate functional cell types of neural crest and CNS lineages and could be used to model a human pain syndrome via gene editing of SCN9A in iNBSCs. NBSCs can also be derived from human pluripotent stem cells and share functional and molecular features with NBSCs isolated from embryonic day 8.5 (E8.5) mouse neural folds. Single-cell RNA sequencing identified the anterior hindbrain as the origin of mouse NBSCs, with human iNBSCs sharing a similar regional identity. In summary, we identify embryonic NBSCs and report their generation by direct reprogramming in human, which may facilitate insights into neural development and provide a neural stem cell source for applications in regenerative medicine.

In Extended Data Fig. 1a of this Letter, the flow cytometry plot depicting the surface phenotype of AML sample DD08 was a duplicate of the plot for AML sample DD06. Supplementary Data 4 has been added to the Supplementary Information of the original Letter to clarify the proteome data acquisition and presentation. The original Letter has been corrected online.

The branched chain amino acid (BCAA) pathway and high levels of BCAA transaminase 1 (BCAT1) have recently been associated with aggressiveness in several cancer entities. However, the mechanistic role of BCAT1 in this process remains largely uncertain. By performing high-resolution proteomic analysis of human acute myeloid leukaemia (AML) stem cell (LSC) and non-LSC populations, we found the BCAA pathway enriched and BCAT1 protein and transcripts overexpressed in LSCs. We show that BCAT1, which transfers α-amino groups from BCAAs to α-ketoglutarate (αKG), is a critical regulator of intracellular αKG homeostasis. Next to its role in the tricarboxylic acid (TCA) cycle, αKG is an essential co-factor for αKG-dependent dioxygenases such as EGLN1 and the TET family of DNA demethylases. Knockdown of BCAT1 in leukaemia cells caused accumulation of αKG leading to EGLN1-mediated HIF1α protein degradation. This resulted in a growth and survival defect and abrogated leukaemia-initiating potential. In contrast, overexpression (OE) of BCAT1 in leukaemia cells decreased intracellular αKG levels and caused DNA hypermethylation via altered TET activity. BCAT1highAMLs displayed a DNA hypermethylation phenotype similar to cases carrying mutant isocitrate dehydrogenase (IDHmut), in which TET2 is inhibited by the oncometabolite 2-hydroxyglutarate. High levels of BCAT1 strongly correlate with shorter overall survival in IDHwtTET2wt, but not IDHmutor TET2mutAMLs. Gene sets characteristic for IDHmutAMLs were enriched in IDHwtTETwtBCAT1highpatient samples. BCAT1highAMLs showed robust enrichment for LSC signatures, and paired sample analysis revealed a significant increase in BCAT1 levels upon disease relapse. In summary, by limiting intracellular αKG, BCAT1 links BCAA catabolism to HIF1α stability and regulation of the epigenomic landscape mimicking the effects of IDH mutations. In essence, our data strongly suggest that suppression of αKG-dependent dioxygenases constitutes a common feature of AML and, next to mutations in IDH and TET2, we identify the BCAA-BCAT1-αKG axis as alternative, non-genetic mechanism in LSCs to control the activity of these enzymes. Strategies to increase αKG levels, i.e. by inhibition of BCAT1 in order to compromise cancer stem cell function may lead to improved clinical outcomes of patients suffering from AML or other tumour entities expressing high levels of BCAT1 in the absence of IDH or TET mutations.

This corrects the article DOI: 10.1038/ni.3799

No relevant conflicts of interest to declare.

Acute Myeloid Leukemia (AML) is a heterogeneous group of myeloid malignancies. The classification of AML subtypes is based on recurrent genetic abnormalities and typical histopathological features, which impact on the patient’s prognosis. AML is also the model disease for the cancer stem cell model with leukemia stem cells (LSCs) residing at the top of a hierarchical organization. LSCs have self-renewal activity and generate leukemic progeny, which make up the majority of leukemic cells. Because LSCs can be quiescent and reside in specific niches in the bone marrow, rendering them resistant to conventional chemotherapy approaches, they are considered the source of relapse. Hence, further strategies to eradicate LSCs are pivotal to improve patient outcomes of this dismal disease.