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.

Loss of TET2 in hematopoietic cells leads to DNA hypermethylation of active enhancers and induction of leukemogenesis.

DNA methylation is tightly regulated throughout mammalian development, and altered DNA methylation patterns are a general hallmark of cancer. The methylcytosine dioxygenase TET2 is frequently mutated in hematological disorders, including acute myeloid leukemia (AML), and has been suggested to protect CG dinucleotide (CpG) islands and promoters from aberrant DNA methylation. In this study, we present a novel Tet2-dependent leukemia mouse model that closely recapitulates gene expression profiles and hallmarks of human AML1-ETO-induced AML. Using this model, we show that the primary effect of Tet2 loss in preleukemic hematopoietic cells is progressive and widespread DNA hypermethylation affecting up to 25% of active enhancer elements. In contrast, CpG island and promoter methylation does not change in a Tet2-dependent manner but increases relative to population doublings. We confirmed this specific enhancer hypermethylation phenotype in human AML patients with TET2 mutations. Analysis of immediate gene expression changes reveals rapid deregulation of a large number of genes implicated in tumorigenesis, including many down-regulated tumor suppressor genes. Hence, we propose that TET2 prevents leukemic transformation by protecting enhancers from aberrant DNA methylation and that it is the combined silencing of several tumor suppressor genes in TET2 mutated hematopoietic cells that contributes to increased stem cell proliferation and leukemogenesis.

Transforming activities of the NUP98-KMT2A fusion gene associated with myelodysplasia and acute myeloid leukemia.

Inv(11)(p15q23), found in myelodysplastic syndromes and acute myeloid leukemia, leads to expression of a fusion protein consisting of the N-terminal of nucleoporin 98 (NUP98) and the majority of the lysine methyltransferase 2A (KMT2A). To explore the transforming potential of this fusion we established inducible iNUP98-KMT2A transgenic mice. After a median latency of 80 weeks, over 90% of these mice developed signs of disease, with anemia and reduced bone marrow cellularity, increased white blood cell numbers, extramedullary hematopoiesis, and multilineage dysplasia. Additionally, induction of iNUP98-KMT2A led to elevated lineage marker-negative Sca-1+ c-Kit+ cell numbers in the bone marrow, which outcompeted wildtype cells in repopulation assays. Six iNUP98-KMT2A mice developed transplantable acute myeloid leukemia with leukemic blasts infiltrating multiple organs. Notably, as reported for patients, iNUP98-KMT2A leukemic blasts did not express increased levels of the HoxA-B-C gene cluster, and in contrast to KMT2A-AF9 leukemic cells, the cells were resistant to pharmacological targeting of menin and BET family proteins by MI-2-2 or JQ1, respectively. Expression of iNUP98-KMT2A in mouse embryonic fibroblasts led to an accumulation of cells in G1 phase, and abrogated replicative senescence. In bone marrow-derived hematopoietic progenitors, iNUP98-KMT2A expression similarly resulted in increased cell numbers in the G1 phase of the cell cycle, with aberrant gene expression of Sirt1, Tert, Rbl2, Twist1, Vim, and Prkcd, mimicking that seen in mouse embryonic fibroblasts. In summary, we demonstrate that iNUP98-KMT2A has in vivo transforming activity and interferes with cell cycle progression rather than primarily blocking differentiation.

Peak-valley-peak pattern of histone modifications delineates active regulatory elements and their directionality.

Formation of nucleosome free region (NFR) accompanied by specific histone modifications at flanking nucleosomes is an important prerequisite for enhancer and promoter activity. Due to this process, active regulatory elements often exhibit a distinct shape of histone signal in the form of a peak-valley-peak (PVP) pattern. However, different features of PVP patterns and their robustness in predicting active regulatory elements have never been systematically analyzed. Here, we present PARE, a novel computational method that systematically analyzes the H3K4me1 or H3K4me3 PVP patterns to predict NFRs. We show that NFRs predicted by H3K4me1 and me3 patterns are associated with active enhancers and promoters, respectively. Furthermore, asymmetry in the height of peaks flanking the central valley can predict the directionality of stable transcription at promoters. Using PARE on ChIP-seq histone modifications from four ENCODE cell lines and four hematopoietic differentiation stages, we identified several enhancers whose regulatory activity is stage specific and correlates positively with the expression of proximal genes in a particular stage. In conclusion, our results demonstrate that PVP patterns delineate both the histone modification landscape and the transcriptional activities governed by active enhancers and promoters, and therefore can be used for their prediction. PARE is freely available at http://servers.binf.ku.dk/pare.

Autophagy is required for stem cell mobilization by G-CSF.

Granulocyte colony-stimulating factor (G-CSF) is widely used clinically to prevent neutropenia after cytotoxic chemotherapy and to mobilize hematopoietic stem cells (HSCs) for transplantation. Autophagy, a process of cytoplasmic component recycling, maintains cellular homeostasis and protects the cell during periods of metabolic stress or nutrient deprivation. We have observed that G-CSF activates autophagy in neutrophils and HSCs from both mouse and human donors. Furthermore, G-CSF-induced neutrophil and HSC mobilization is impaired in the absence of autophagy. In contrast, autophagy is dispensable for direct HSC mobilization in response to the CXCR4 antagonist AMD3100. Altogether, these data demonstrate an important role for G-CSF in invoking autophagy within hematopoietic and myeloid cells and suggest that this pathway is critical for ensuring cell survival in response to clinically relevant cytokine-induced stress. These findings have direct relevance to HSC transplantation and the increasing clinical use of agents that modulate autophagy.

C/EBPα is required for long-term self-renewal and lineage priming of hematopoietic stem cells and for the maintenance of epigenetic configurations in multipotent progenitors.

Transcription factors are key regulators of hematopoietic stem cells (HSCs) and act through their ability to bind DNA and impact on gene transcription. Their functions are interpreted in the complex landscape of chromatin, but current knowledge on how this is achieved is very limited. C/EBPα is an important transcriptional regulator of hematopoiesis, but its potential functions in HSCs have remained elusive. Here we report that C/EBPα serves to protect adult HSCs from apoptosis and to maintain their quiescent state. Consequently, deletion of Cebpa is associated with loss of self-renewal and HSC exhaustion. By combining gene expression analysis with genome-wide assessment of C/EBPα binding and epigenetic configurations, we show that C/EBPα acts to modulate the epigenetic states of genes belonging to molecular pathways important for HSC function. Moreover, our data suggest that C/EBPα acts as a priming factor at the HSC level where it actively promotes myeloid differentiation and counteracts lymphoid lineage choice. Taken together, our results show that C/EBPα is a key regulator of HSC biology, which influences the epigenetic landscape of HSCs in order to balance different cell fate options.

Modification of T cell responses by stem cell mobilization requires direct signaling of the T cell by G-CSF and IL-10.

The majority of allogeneic stem cell transplants are currently undertaken using G-CSF mobilized peripheral blood stem cells. G-CSF has diverse biological effects on a broad range of cells and IL-10 is a key regulator of many of these effects. Using mixed radiation chimeras in which the hematopoietic or nonhematopoietic compartments were wild-type, IL-10(-/-), G-CSFR(-/-), or combinations thereof we demonstrated that the attenuation of alloreactive T cell responses after G-CSF mobilization required direct signaling of the T cell by both G-CSF and IL-10. IL-10 was generated principally by radio-resistant tissue, and was not required to be produced by T cells. G-CSF mobilization significantly modulated the transcription profile of CD4(+)CD25(+) regulatory T cells, promoted their expansion in the donor and recipient and their depletion significantly increased graft-versus-host disease (GVHD). In contrast, stem cell mobilization with the CXCR4 antagonist AMD3100 did not alter the donor T cell's ability to induce acute GVHD. These studies provide an explanation for the effects of G-CSF on T cell function and demonstrate that IL-10 is required to license regulatory function but T cell production of IL-10 is not itself required for the attenuation GVHD. Although administration of CXCR4 antagonists is an efficient means of stem cell mobilization, this fails to evoke the immunomodulatory effects seen during G-CSF mobilization. These data provide a compelling rationale for considering the immunological benefits of G-CSF in selecting mobilization protocols for allogeneic stem cell transplantation.

Amplification of pico-scale DNA mediated by bacterial carrier DNA for small-cell-number transcription factor ChIP-seq.

BACKGROUND: Chromatin-Immunoprecipitation coupled with deep sequencing (ChIP-seq) is used to map transcription factor occupancy and generate epigenetic profiles genome-wide. The requirement of nano-scale ChIP DNA for generation of sequencing libraries has impeded ChIP-seq on in vivo tissues of low cell numbers. RESULTS: We describe a robust, simple and scalable methodology for ChIP-seq of low-abundant cell populations, verified down to 10,000 cells. By employing non-mammalian genome mapping bacterial carrier DNA during amplification, we reliably amplify down to 50 pg of ChIP DNA from transcription factor (CEBPA) and histone mark (H3K4me3) ChIP. We further demonstrate that genomic profiles are highly resilient to changes in carrier DNA to ChIP DNA ratios. CONCLUSIONS: This represents a significant advance compared to existing technologies, which involve either complex steps of pre-selection for nucleosome-containing chromatin or pre-amplification of precipitated DNA, making them prone to introduce experimental biases.

Combining Homologous Recombination-Deficient Testing and Functional RAD51 Analysis Enhances the Prediction of Poly(ADP-Ribose) Polymerase Inhibitor Sensitivity.

PURPOSE: To meet the urgent need for accessible homologous recombination-deficient (HRD) test options, we validated a laboratory-developed test (LDT) and a functional RAD51 assay to assess patients with ovarian cancer and predict the clinical benefit of poly(ADP-ribose) polymerase inhibitor therapy. METHODS: Optimization of the LDT cutoff and validation on the basis of samples from 91 patients enrolled in the ENGOT-ov24/NSGO-AVANOVA1&2 trial (ClinicalTrials.gov identifier: NCT02354131), previously subjected to commercial CDx HRD testing (CDx). RAD51 foci analysis was performed and tumors with ≥five foci/nucleus were classified as RAD51-positive (homologous recombination-proficient). RESULTS: The optimal LDT cutoff is 54. Comparing CDx genome instability score and LDT HRD scores show a Spearman's correlation of rho = 0.764 (P < .0001). Cross-tabulation analysis shows that the sensitivity of the LDT HRD score is 86% and of the LDT HRD status is 91.8% (Fisher's exact test P < .001). Survival analysis on progression-free survival (PFS) of LDT-assessed patients show a Cox regression P < .05. RAD51 assays show a correlation between low RAD51 foci detection (<20% RAD51+ cells) and significantly prolonged PFS (P < .001). CONCLUSION: The robust concordance between the open standard LDT and the CDx, especially the correlation with PFS, warrants future validation and implementation of the open standard LDT for HRD testing in diagnostic settings.

Transcriptional characterization of human megakaryocyte polyploidization and lineage commitment.

BACKGROUND: Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs. OBJECTIVES: To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis. RESULTS: We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state. CONCLUSIONS: Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC.

ERG Controls B Cell Development by Promoting Igh V-to-DJ Recombination.

B cell development depends on the coordinated expression and cooperation of several transcription factors. Here we show that the transcription factor ETS-related gene (ERG) is crucial for normal B cell development and that its deletion results in a substantial loss of bone marrow B cell progenitors and peripheral B cells, as well as a skewing of splenic B cell populations. We find that ERG-deficient B lineage cells exhibit an early developmental block at the pre-B cell stage and proliferate less. The cells fail to express the immunoglobulin heavy chain due to inefficient V-to-DJ recombination, and cells that undergo recombination display a strong bias against incorporation of distal V gene segments. Furthermore, antisense transcription at PAX5-activated intergenic repeat (PAIR) elements, located in the distal region of the Igh locus, depends on ERG. These findings show that ERG serves as a critical regulator of B cell development by ensuring efficient and balanced V-to-DJ recombination.

A Gain-of-Function p53-Mutant Oncogene Promotes Cell Fate Plasticity and Myeloid Leukemia through the Pluripotency Factor FOXH1.

Mutations in the TP53 tumor suppressor gene are common in many cancer types, including the acute myeloid leukemia (AML) subtype known as complex karyotype AML (CK-AML). Here, we identify a gain-of-function (GOF) Trp53 mutation that accelerates CK-AML initiation beyond p53 loss and, surprisingly, is required for disease maintenance. The Trp53R172H mutation (TP53R175H in humans) exhibits a neomorphic function by promoting aberrant self-renewal in leukemic cells, a phenotype that is present in hematopoietic stem and progenitor cells (HSPC) even prior to their transformation. We identify FOXH1 as a critical mediator of mutant p53 function that binds to and regulates stem cell-associated genes and transcriptional programs. Our results identify a context where mutant p53 acts as a bona fide oncogene that contributes to the pathogenesis of CK-AML and suggests a common biological theme for TP53 GOF in cancer. SIGNIFICANCE: Our study demonstrates how a GOF p53 mutant can hijack an embryonic transcription factor to promote aberrant self-renewal. In this context, mutant Trp53 functions as an oncogene to both initiate and sustain myeloid leukemia and suggests a potential convergent activity of mutant Trp53 across cancer types.This article is highlighted in the In This Issue feature, p. 813.

Genome-wide analysis of differential transcriptional and epigenetic variability across human immune cell types.

BACKGROUND: A healthy immune system requires immune cells that adapt rapidly to environmental challenges. This phenotypic plasticity can be mediated by transcriptional and epigenetic variability. RESULTS: We apply a novel analytical approach to measure and compare transcriptional and epigenetic variability genome-wide across CD14+CD16- monocytes, CD66b+CD16+ neutrophils, and CD4+CD45RA+ naïve T cells from the same 125 healthy individuals. We discover substantially increased variability in neutrophils compared to monocytes and T cells. In neutrophils, genes with hypervariable expression are found to be implicated in key immune pathways and are associated with cellular properties and environmental exposure. We also observe increased sex-specific gene expression differences in neutrophils. Neutrophil-specific DNA methylation hypervariable sites are enriched at dynamic chromatin regions and active enhancers. CONCLUSIONS: Our data highlight the importance of transcriptional and epigenetic variability for the key role of neutrophils as the first responders to inflammatory stimuli. We provide a resource to enable further functional studies into the plasticity of immune cells, which can be accessed from: http://blueprint-dev.bioinfo.cnio.es/WP10/hypervariability .

A computational method for identification of vaccine targets from protein regions of conserved human leukocyte antigen binding.

BACKGROUND: Computational methods for T cell-based vaccine target discovery focus on selection of highly conserved peptides identified across pathogen variants, followed by prediction of their binding of human leukocyte antigen molecules. However, experimental studies have shown that T cells often target diverse regions in highly variable viral pathogens and this diversity may need to be addressed through redefinition of suitable peptide targets. METHODS: We have developed a method for antigen assessment and target selection for polyvalent vaccines, with which we identified immune epitopes from variable regions, where all variants bind HLA. These regions, although variable, can thus be considered stable in terms of HLA binding and represent valuable vaccine targets. RESULTS: We applied this method to predict CD8+ T-cell targets in influenza A H7N9 hemagglutinin and significantly increased the number of potential vaccine targets compared to the number of targets discovered using the traditional approach where low-frequency peptides are excluded. CONCLUSIONS: We developed a webserver with an intuitive visualization scheme for summarizing the T cell-based antigenic potential of any given protein or proteome using human leukocyte antigen binding predictions and made a web-accessible software implementation freely available at http://met-hilab.cbs.dtu.dk/blockcons/.

Telomerase inhibition effectively targets mouse and human AML stem cells and delays relapse following chemotherapy.

Acute myeloid leukemia (AML) is an aggressive and lethal blood cancer maintained by rare populations of leukemia stem cells (LSCs). Selective targeting of LSCs is a promising approach for treating AML and preventing relapse following chemotherapy, and developing such therapeutic modalities is a key priority. Here, we show that targeting telomerase activity eradicates AML LSCs. Genetic deletion of the telomerase subunit Terc in a retroviral mouse AML model induces cell-cycle arrest and apoptosis of LSCs, and depletion of telomerase-deficient LSCs is partially rescued by p53 knockdown. Murine Terc(-/-) LSCs express a specific gene expression signature that can be identified in human AML patient cohorts and is positively correlated with patient survival following chemotherapy. In xenografts of primary human AML, genetic or pharmacological inhibition of telomerase targets LSCs, impairs leukemia progression, and delays relapse following chemotherapy. Altogether, these results establish telomerase inhibition as an effective strategy for eliminating AML LSCs.

Lack of the p42 form of C/EBPα leads to spontaneous immortalization and lineage infidelity of committed myeloid progenitors.

Acute myeloid leukemia (AML) develops via a multistep process involving several genetic and epigenetic events, which ultimately leads to the formation of a heterogeneous population of malignant cells, of which only a small subpopulation termed the leukemia initiating cell (LIC) is able to sustain the leukemia. The identity of the LIC is highly diverse and ranges from populations resembling hematopoietic stem cells or multipotent progenitors (MPPs) to more committed myeloid progenitors, and the question still remains whether this is a direct consequence of which cells are targets of the final transforming events. In this study, we use premalignant cells from a Cebpa mutant AML model, in which the LIC population resembles granulocyte-macrophage progenitors (GMPs), to show that premalignant GMPs undergo spontaneous immortalization with a high clonal frequency when cultured in vitro, suggesting that these cells constitute the target of the final transforming events. Furthermore, we show that premalignant GMPs are characterized by a distinct T cell gene expression signature correlating with an increased potential for differentiation toward the T cell lineage. These findings have implications for our understanding of the transcriptional wiring in premalignant myeloid progenitors and how this contributes to the development of AML.