Increased bleeding and thrombosis in myeloproliferative neoplasms mediated through altered expression of inherited platelet disorder genes.

An altered thrombo-hemorrhagic profile has long been observed in patients with myeloproliferative neoplasms (MPNs). We hypothesized that this observed clinical phenotype may result from altered expression of genes known to harbor genetic variants in bleeding, thrombotic, or platelet disorders. Here, we identify 32 genes from a clinically validated gene panel that were also significantly differentially expressed in platelets from MPN patients as opposed to healthy donors. This work begins to unravel previously unclear mechanisms underlying an important clinical reality in MPNs. Knowledge of altered platelet gene expression in MPN thrombosis/bleeding diathesis opens opportunities to advance clinical care by: (1) enabling risk stratification, in particular, for patients undergoing invasive procedures, and (2) facilitating tailoring of treatment strategies for those at highest risk, for example, in the form of antifibrinolytics, desmopressin or platelet transfusions (not current routine practice). Marker genes identified in this work may also enable prioritization of candidates in future MPN mechanistic as well as outcome studies.

Non-alcoholic fatty liver disease is characterised by a reduced polyunsaturated fatty acid transport via free fatty acids and high-density lipoproteins (HDL).

BACKGROUND AND OBJECTIVES: Non-alcoholic fatty liver disease (NAFLD) develops due to impaired hepatic lipid fluxes and is a risk factor for chronic liver disease and atherosclerosis. Lipidomic studies consistently reported characteristic hepatic/VLDL "lipid signatures" in NAFLD; whole plasma traits are more debated. Surprisingly, the HDL lipid composition by mass spectrometry has not been characterised across the NAFLD spectrum, despite HDL being a possible source of hepatic lipids delivered from peripheral tissues alongside free fatty acids (FFA). This study characterises the HDL lipidomic signature in NAFLD, and its correlation with metabolic and liver disease markers. METHODS: We used liquid chromatography-mass spectrometry to determine the whole serum and HDL lipidomic profile in 89 biopsy-proven NAFLD patients and 20 sex and age-matched controls. RESULTS: In the whole serum of NAFLD versus controls, we report a depletion in polyunsaturated (PUFA) phospholipids (PL) and FFA; with PUFA PL being also lower in HDL, and negatively correlated with BMI, insulin resistance, triglycerides, and hepatocyte ballooning. In the HDL of the NAFLD group we also describe higher saturated ceramides, which positively correlate with insulin resistance and transaminases. CONCLUSION: NAFLD features lower serum lipid species containing polyunsaturated fatty acids; the most affected lipid fractions are FFA and (HDL) phospholipids; our data suggest a possible defect in the transfer of PUFA from peripheral tissues to the liver in NAFLD. Mechanistic studies are required to explore the biological implications of our findings addressing if HDL composition can influence liver metabolism and damage, thus contributing to NAFLD pathophysiology.

Somatostatin Receptor PET/MR Imaging of Inflammation in Patients With Large Vessel Vasculitis and Atherosclerosis.

BACKGROUND: Assessing inflammatory disease activity in large vessel vasculitis (LVV) can be challenging by conventional measures. OBJECTIVES: We aimed to investigate somatostatin receptor 2 (SST2) as a novel inflammation-specific molecular imaging target in LVV. METHODS: In a prospective, observational cohort study, in vivo arterial SST2 expression was assessed by positron emission tomography/magnetic resonance imaging (PET/MRI) using 68Ga-DOTATATE and 18F-FET-ßAG-TOCA. Ex vivo mapping of the imaging target was performed using immunofluorescence microscopy; imaging mass cytometry; and bulk, single-cell, and single-nucleus RNA sequencing. RESULTS: Sixty-one participants (LVV: n = 27; recent atherosclerotic myocardial infarction of ≤2 weeks: n = 25; control subjects with an oncologic indication for imaging: n = 9) were included. Index vessel SST2 maximum tissue-to-blood ratio was 61.8% (P < 0.0001) higher in active/grumbling LVV than inactive LVV and 34.6% (P = 0.0002) higher than myocardial infarction, with good diagnostic accuracy (area under the curve: ≥0.86; P < 0.001 for both). Arterial SST2 signal was not elevated in any of the control subjects. SST2 PET/MRI was generally consistent with 18F-fluorodeoxyglucose PET/computed tomography imaging in LVV patients with contemporaneous clinical scans but with very low background signal in the brain and heart, allowing for unimpeded assessment of nearby coronary, myocardial, and intracranial artery involvement. Clinically effective treatment for LVV was associated with a 0.49 ± 0.24 (standard error of the mean [SEM]) (P = 0.04; 22.3%) reduction in the SST2 maximum tissue-to-blood ratio after 9.3 ± 3.2 months. SST2 expression was localized to macrophages, pericytes, and perivascular adipocytes in vasculitis specimens, with specific receptor binding confirmed by autoradiography. SSTR2-expressing macrophages coexpressed proinflammatory markers. CONCLUSIONS: SST2 PET/MRI holds major promise for diagnosis and therapeutic monitoring in LVV. (PET Imaging of Giant Cell and Takayasu Arteritis [PITA], NCT04071691; Residual Inflammation and Plaque Progression Long-Term Evaluation [RIPPLE], NCT04073810).

Transcriptional, epigenetic and metabolic signatures in cardiometabolic syndrome defined by extreme phenotypes.

BACKGROUND: This work is aimed at improving the understanding of cardiometabolic syndrome pathophysiology and its relationship with thrombosis by generating a multi-omic disease signature. METHODS/RESULTS: We combined classic plasma biochemistry and plasma biomarkers with the transcriptional and epigenetic characterisation of cell types involved in thrombosis, obtained from two extreme phenotype groups (morbidly obese and lipodystrophy) and lean individuals to identify the molecular mechanisms at play, highlighting patterns of abnormal activation in innate immune phagocytic cells. Our analyses showed that extreme phenotype groups could be distinguished from lean individuals, and from each other, across all data layers. The characterisation of the same obese group, 6 months after bariatric surgery, revealed the loss of the abnormal activation of innate immune cells previously observed. However, rather than reverting to the gene expression landscape of lean individuals, this occurred via the establishment of novel gene expression landscapes. NETosis and its control mechanisms emerge amongst the pathways that show an improvement after surgical intervention. CONCLUSIONS: We showed that the morbidly obese and lipodystrophy groups, despite some differences, shared a common cardiometabolic syndrome signature. We also showed that this could be used to discriminate, amongst the normal population, those individuals with a higher likelihood of presenting with the disease, even when not displaying the classic features.

Unique molecular and functional features of extramedullary hematopoietic stem and progenitor cell reservoirs in humans.

Rare hematopoietic stem and progenitor cell (HSPC) pools outside the bone marrow (BM) contribute to blood production in stress and disease but remain ill-defined. Although nonmobilized peripheral blood (PB) is routinely sampled for clinical management, the diagnosis and monitoring potential of PB HSPCs remain untapped, as no healthy PB HSPC baseline has been reported. Here we comprehensively delineate human extramedullary HSPC compartments comparing spleen, PB, and mobilized PB to BM using single-cell RNA-sequencing and/or functional assays. We uncovered HSPC features shared by extramedullary tissues and others unique to PB. First, in contrast to actively dividing BM HSPCs, we found no evidence of substantial ongoing hematopoiesis in extramedullary tissues at steady state but report increased splenic HSPC proliferative output during stress erythropoiesis. Second, extramedullary hematopoietic stem cells/multipotent progenitors (HSCs/MPPs) from spleen, PB, and mobilized PB share a common transcriptional signature and increased abundance of lineage-primed subsets compared with BM. Third, healthy PB HSPCs display a unique bias toward erythroid-megakaryocytic differentiation. At the HSC/MPP level, this is functionally imparted by a subset of phenotypic CD71+ HSCs/MPPs, exclusively producing erythrocytes and megakaryocytes, highly abundant in PB but rare in other adult tissues. Finally, the unique erythroid-megakaryocytic-skewing of PB is perturbed with age in essential thrombocythemia and ß-thalassemia. Collectively, we identify extramedullary lineage-primed HSPC reservoirs that are nonproliferative in situ and report involvement of splenic HSPCs during demand-adapted hematopoiesis. Our data also establish aberrant composition and function of circulating HSPCs as potential clinical indicators of BM dysfunction.

Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease.

Neutrophils play fundamental roles in innate immune response, shape adaptive immunity, and are a potentially causal cell type underpinning genetic associations with immune system traits and diseases. Here, we profile the binding of myeloid master regulator PU.1 in primary neutrophils across nearly a hundred volunteers. We show that variants associated with differential PU.1 binding underlie genetically-driven differences in cell count and susceptibility to autoimmune and inflammatory diseases. We integrate these results with other multi-individual genomic readouts, revealing coordinated effects of PU.1 binding variants on the local chromatin state, enhancer-promoter contacts and downstream gene expression, and providing a functional interpretation for 27 genes underlying immune traits. Collectively, these results demonstrate the functional role of PU.1 and its target enhancers in neutrophil transcriptional control and immune disease susceptibility.

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.

Cell type-specific novel long non-coding RNA and circular RNA in the BLUEPRINT hematopoietic transcriptomes atlas.

Transcriptional profiling of hematopoietic cell subpopulations has helped to characterize the developmental stages of the hematopoietic system and the molecular bases of malignant and non-malignant blood diseases. Previously, only the genes targeted by expression microarrays could be profiled genome-wide. High-throughput RNA sequencing, however, encompasses a broader repertoire of RNA molecules, without restriction to previously annotated genes. We analyzed the BLUEPRINT consortium RNA-sequencing data for mature hematopoietic cell types. The data comprised 90 total RNA-sequencing samples, each composed of one of 27 cell types, and 32 small RNA-sequencing samples, each composed of one of 11 cell types. We estimated gene and isoform expression levels for each cell type using existing annotations from Ensembl. We then used guided transcriptome assembly to discover unannotated transcripts. We identified hundreds of novel non-coding RNA genes and showed that the majority have cell type-dependent expression. We also characterized the expression of circular RNA and found that these are also cell type-specific. These analyses refine the active transcriptional landscape of mature hematopoietic cells, highlight abundant genes and transcriptional isoforms for each blood cell type, and provide a valuable resource for researchers of hematologic development and diseases. Finally, we made the data accessible via a web-based interface: https://blueprint.haem.cam.ac.uk/bloodatlas/.

Whole-genome sequencing of patients with rare diseases in a national health system.

Most patients with rare diseases do not receive a molecular diagnosis and the aetiological variants and causative genes for more than half such disorders remain to be discovered1. Here we used whole-genome sequencing (WGS) in a national health system to streamline diagnosis and to discover unknown aetiological variants in the coding and non-coding regions of the genome. We generated WGS data for 13,037 participants, of whom 9,802 had a rare disease, and provided a genetic diagnosis to 1,138 of the 7,065 extensively phenotyped participants. We identified 95 Mendelian associations between genes and rare diseases, of which 11 have been discovered since 2015 and at least 79 are confirmed to be aetiological. By generating WGS data of UK Biobank participants2, we found that rare alleles can explain the presence of some individuals in the tails of a quantitative trait for red blood cells. Finally, we identified four novel non-coding variants that cause disease through the disruption of transcription of ARPC1B, GATA1, LRBA and MPL. Our study demonstrates a synergy by using WGS for diagnosis and aetiological discovery in routine healthcare.

Sphingolipid dysregulation due to lack of functional KDSR impairs proplatelet formation causing thrombocytopenia.

Sphingolipids are fundamental to membrane trafficking, apoptosis, and cell differentiation and proliferation. KDSR or 3-keto-dihydrosphingosine reductase is an essential enzyme for de novo sphingolipid synthesis, and pathogenic mutations in KDSR result in the severe skin disorder erythrokeratodermia variabilis et progressiva-4 Four of the eight reported cases also had thrombocytopenia but the underlying mechanism has remained unexplored. Here we expand upon the phenotypic spectrum of KDSR deficiency with studies in two siblings with novel compound heterozygous variants associated with thrombocytopenia, anemia, and minimal skin involvement. We report a novel phenotype of progressive juvenile myelofibrosis in the propositus, with spontaneous recovery of anemia and thrombocytopenia in the first decade of life. Examination of bone marrow biopsies showed megakaryocyte hyperproliferation and dysplasia. Megakaryocytes obtained by culture of CD34+ stem cells confirmed hyperproliferation and showed reduced proplatelet formation. The effect of KDSR insufficiency on the sphingolipid profile was unknown, and was explored in vivo and in vitro by a broad metabolomics screen that indicated activation of an in vivo compensatory pathway that leads to normalization of downstream metabolites such as ceramide. Differentiation of propositus-derived induced pluripotent stem cells to megakaryocytes followed by expression of functional KDSR showed correction of the aberrant cellular and biochemical phenotypes, corroborating the critical role of KDSR in proplatelet formation. Finally, Kdsr depletion in zebrafish recapitulated the thrombocytopenia and showed biochemical changes similar to those observed in the affected siblings. These studies support an important role for sphingolipids as regulators of cytoskeletal organization during megakaryopoiesis and proplatelet formation.

Histone Acetyltransferase KAT2A Stabilizes Pluripotency with Control of Transcriptional Heterogeneity.

Cell fate transitions in mammalian stem cell systems have often been associated with transcriptional heterogeneity; however, existing data have failed to establish a functional or mechanistic link between the two phenomena. Experiments in unicellular organisms support the notion that transcriptional heterogeneity can be used to facilitate adaptability to environmental changes and have identified conserved chromatin-associated factors that modulate levels of transcriptional noise. Herein, we show destabilization of pluripotency-associated gene regulatory networks through increased transcriptional heterogeneity of mouse embryonic stem cells in which paradigmatic histone acetyl-transferase, and candidate noise modulator, Kat2a (yeast orthologue Gcn5), have been inhibited. Functionally, network destabilization associates with reduced pluripotency and accelerated mesendodermal differentiation, with increased probability of transitions into lineage commitment. Thus, we show evidence of a relationship between transcriptional heterogeneity and cell fate transitions through manipulation of the histone acetylation landscape of mouse embryonic stem cells, suggesting a general principle that could be exploited in other normal and malignant stem cell fate transitions. Stem Cells 2018;36:1828-11.

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.

SNP in human ARHGEF3 promoter is associated with DNase hypersensitivity, transcript level and platelet function, and Arhgef3 KO mice have increased mean platelet volume.

Genome-wide association studies have identified a genetic variant at 3p14.3 (SNP rs1354034) that strongly associates with platelet number and mean platelet volume in humans. While originally proposed to be intronic, analysis of mRNA expression in primary human hematopoietic subpopulations reveals that this SNP is located directly upstream of the predominantly expressed ARHGEF3 isoform in megakaryocytes (MK). We found that ARHGEF3, which encodes a Rho guanine exchange factor, is dramatically upregulated during both human and murine MK maturation. We show that the SNP (rs1354034) is located in a DNase I hypersensitive region in human MKs and is an expression quantitative locus (eQTL) associated with ARHGEF3 expression level in human platelets, suggesting that it may be the causal SNP that accounts for the variations observed in human platelet traits and ARHGEF3 expression. In vitro human platelet activation assays revealed that rs1354034 is highly correlated with human platelet activation by ADP. In order to test whether ARHGEF3 plays a role in MK development and/or platelet function, we developed an Arhgef3 KO/LacZ reporter mouse model. Reflecting changes in gene expression, LacZ expression increases during MK maturation in these mice. Although Arhgef3 KO mice have significantly larger platelets, loss of Arhgef3 does not affect baseline MK or platelets nor does it affect platelet function or platelet recovery in response to antibody-mediated platelet depletion compared to littermate controls. In summary, our data suggest that modulation of ARHGEF3 gene expression in humans with a promoter-localized SNP plays a role in human MKs and human platelet function-a finding resulting from the biological follow-up of human genetic studies. Arhgef3 KO mice partially recapitulate the human phenotype.

Megakaryocytes in Myeloproliferative Neoplasms Have Unique Somatic Mutations.

Myeloproliferative neoplasms (MPNs) are a group of related clonal hemopoietic stem cell disorders associated with hyperproliferation of myeloid cells. They are driven by mutations in the hemopoietic stem cell, most notably JAK2V617F, CALR, and MPL. Clinically, they have the propensity to progress to myelofibrosis and transform to acute myeloid leukemia. Megakaryocytic hyperplasia with abnormal features are characteristic, and it is thought that these cells stimulate and drive fibrotic progression. The biological defects underpinning this remain to be explained. In this study we examined the megakaryocyte genome in 12 patients with MPNs to determine whether there are somatic variants and whether there is any association with marrow fibrosis. We performed targeted next-generation sequencing for 120 genes associated with myeloid neoplasms on megakaryocytes isolated from aspirated bone marrow. Ten of the 12 patients had genomic defects in megakaryocytes that were not present in nonmegakaryocytic hemopoietic marrow cells from the same patient. The greatest allelic burden was in patients with increased reticulin deposition. The megakaryocyte-unique mutations were predominantly in genes that regulate chromatin remodeling, chromosome alignment, and stability. These findings show that genomic abnormalities are present in megakaryocytes in MPNs and that these appear to be associated with progression to bone marrow fibrosis.

Detection of Atherosclerotic Inflammation by 68Ga-DOTATATE PET Compared to [18F]FDG PET Imaging.

BACKGROUND: Inflammation drives atherosclerotic plaque rupture. Although inflammation can be measured using fluorine-18-labeled fluorodeoxyglucose positron emission tomography ([18F]FDG PET), [18F]FDG lacks cell specificity, and coronary imaging is unreliable because of myocardial spillover. OBJECTIVES: This study tested the efficacy of gallium-68-labeled DOTATATE (68Ga-DOTATATE), a somatostatin receptor subtype-2 (SST2)-binding PET tracer, for imaging atherosclerotic inflammation. METHODS: We confirmed 68Ga-DOTATATE binding in macrophages and excised carotid plaques. 68Ga-DOTATATE PET imaging was compared to [18F]FDG PET imaging in 42 patients with atherosclerosis. RESULTS: Target SSTR2 gene expression occurred exclusively in "proinflammatory" M1 macrophages, specific 68Ga-DOTATATE ligand binding to SST2 receptors occurred in CD68-positive macrophage-rich carotid plaque regions, and carotid SSTR2 mRNA was highly correlated with in vivo 68Ga-DOTATATE PET signals (r = 0.89; 95% confidence interval [CI]: 0.28 to 0.99; p = 0.02). 68Ga-DOTATATE mean of maximum tissue-to-blood ratios (mTBRmax) correctly identified culprit versus nonculprit arteries in patients with acute coronary syndrome (median difference: 0.69; interquartile range [IQR]: 0.22 to 1.15; p = 0.008) and transient ischemic attack/stroke (median difference: 0.13; IQR: 0.07 to 0.32; p = 0.003). 68Ga-DOTATATE mTBRmax predicted high-risk coronary computed tomography features (receiver operating characteristics area under the curve [ROC AUC]: 0.86; 95% CI: 0.80 to 0.92; p < 0.0001), and correlated with Framingham risk score (r = 0.53; 95% CI: 0.32 to 0.69; p <0.0001) and [18F]FDG uptake (r = 0.73; 95% CI: 0.64 to 0.81; p < 0.0001). [18F]FDG mTBRmax differentiated culprit from nonculprit carotid lesions (median difference: 0.12; IQR: 0.0 to 0.23; p = 0.008) and high-risk from lower-risk coronary arteries (ROC AUC: 0.76; 95% CI: 0.62 to 0.91; p = 0.002); however, myocardial [18F]FDG spillover rendered coronary [18F]FDG scans uninterpretable in 27 patients (64%). Coronary 68Ga-DOTATATE PET scans were readable in all patients. CONCLUSIONS: We validated 68Ga-DOTATATE PET as a novel marker of atherosclerotic inflammation and confirmed that 68Ga-DOTATATE offers superior coronary imaging, excellent macrophage specificity, and better power to discriminate high-risk versus low-risk coronary lesions than [18F]FDG. (Vascular Inflammation Imaging Using Somatostatin Receptor Positron Emission Tomography [VISION]; NCT02021188).

CSB ablation induced apoptosis is mediated by increased endoplasmic reticulum stress response.

The DNA repair protein Cockayne syndrome group B (CSB) has been recently identified as a promising anticancer target. Suppression, by antisense technology, of this protein causes devastating effects on tumor cells viability, through a massive induction of apoptosis, while being non-toxic to non-transformed cells. To gain insights into the mechanisms underlying the pro-apoptotic effects observed after CSB ablation, global gene expression patterns were determined, to identify genes that were significantly differentially regulated as a function of CSB expression. Our findings revealed that response to endoplasmic reticulum stress and response to unfolded proteins were ranked top amongst the cellular processes affected by CSB suppression. The major components of the endoplasmic reticulum stress-mediated apoptosis pathway, including pro-apoptotic factors downstream of the ATF3-CHOP cascade, were dramatically up-regulated. Altogether our findings add new pieces to the understanding of CSB mechanisms of action and to the molecular basis of CS syndrome.

DNA Methylation Dynamics of Human Hematopoietic Stem Cell Differentiation.

Hematopoietic stem cells give rise to all blood cells in a differentiation process that involves widespread epigenome remodeling. Here we present genome-wide reference maps of the associated DNA methylation dynamics. We used a meta-epigenomic approach that combines DNA methylation profiles across many small pools of cells and performed single-cell methylome sequencing to assess cell-to-cell heterogeneity. The resulting dataset identified characteristic differences between HSCs derived from fetal liver, cord blood, bone marrow, and peripheral blood. We also observed lineage-specific DNA methylation between myeloid and lymphoid progenitors, characterized immature multi-lymphoid progenitors, and detected progressive DNA methylation differences in maturing megakaryocytes. We linked these patterns to gene expression, histone modifications, and chromatin accessibility, and we used machine learning to derive a model of human hematopoietic differentiation directly from DNA methylation data. Our results contribute to a better understanding of human hematopoietic stem cell differentiation and provide a framework for studying blood-linked diseases.

Distinct Trends of DNA Methylation Patterning in the Innate and Adaptive Immune Systems.

DNA methylation and the localization and post-translational modification of nucleosomes are interdependent factors that contribute to the generation of distinct phenotypes from genetically identical cells. With 112 whole-genome bisulfite sequencing datasets from the BLUEPRINT Epigenome Project, we analyzed the global development of DNA methylation patterns during lineage commitment and maturation of a range of immune system effector cells and the cancers that arise from them. We show clear trends in methylation patterns that are distinct in the innate and adaptive arms of the human immune system, both globally and in relation to consistently positioned nucleosomes. Most notable are a progressive loss of methylation in developing lymphocytes and the consistent occurrence of non-CG methylation in specific cell types. Cancer samples from the two lineages are further polarized, suggesting the involvement of distinct lineage-specific epigenetic mechanisms. We anticipate broad utility for this resource as a basis for further comparative epigenetic analyses.

eFORGE: A Tool for Identifying Cell Type-Specific Signal in Epigenomic Data.

Epigenome-wide association studies (EWAS) provide an alternative approach for studying human disease through consideration of non-genetic variants such as altered DNA methylation. To advance the complex interpretation of EWAS, we developed eFORGE (http://eforge.cs.ucl.ac.uk/), a new standalone and web-based tool for the analysis and interpretation of EWAS data. eFORGE determines the cell type-specific regulatory component of a set of EWAS-identified differentially methylated positions. This is achieved by detecting enrichment of overlap with DNase I hypersensitive sites across 454 samples (tissues, primary cell types, and cell lines) from the ENCODE, Roadmap Epigenomics, and BLUEPRINT projects. Application of eFORGE to 20 publicly available EWAS datasets identified disease-relevant cell types for several common diseases, a stem cell-like signature in cancer, and demonstrated the ability to detect cell-composition effects for EWAS performed on heterogeneous tissues. Our approach bridges the gap between large-scale epigenomics data and EWAS-derived target selection to yield insight into disease etiology.

Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters.

Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases.