Critical Modulation of Hematopoietic Lineage Fate by Hepatic Leukemia Factor

CUX1 regulates human hematopoietic stem cell chromatin accessibility via the BAF complex

CUX1 is a homeodomain-containing transcription factor that is essential for the development and differentiation of multiple tissues. CUX1 is recurrently mutated or deleted in cancer, particularly in myeloid malignancies. However, the mechanism by which CUX1 regulates gene expression and differentiation remains poorly understood, creating a barrier to understanding the tumor-suppressive functions of CUX1. Here, we demonstrate that CUX1 directs the BAF chromatin remodeling complex to DNA to increase chromatin accessibility in hematopoietic cells. CUX1 preferentially regulates lineage-specific enhancers, and CUX1 target genes are predictive of cell fate in vivo. These data indicate that CUX1 regulates hematopoietic lineage commitment and homeostasis via pioneer factor activity, and CUX1 deficiency disrupts these processes in stem and progenitor cells, facilitating transformation.

Lineage-tracing hematopoietic stem cell origins in vivo to efficiently make human HLF+ HOXA+ hematopoietic progenitors from pluripotent stem cells

The developmental origin of blood-forming hematopoietic stem cells (HSCs) is a longstanding question. Here, our non-invasive genetic lineage tracing in mouse embryos pinpoints that artery endothelial cells generate HSCs. Arteries are transiently competent to generate HSCs for 2.5 days (∼E8.5-E11) but subsequently cease, delimiting a narrow time frame for HSC formation in vivo. Guided by the arterial origins of blood, we efficiently and rapidly differentiate human pluripotent stem cells (hPSCs) into posterior primitive streak, lateral mesoderm, artery endothelium, hemogenic endothelium, and >90% pure hematopoietic progenitors within 10 days. hPSC-derived hematopoietic progenitors generate T, B, NK, erythroid, and myeloid cells in vitro and, critically, express hallmark HSC transcription factors HLF and HOXA5-HOXA10, which were previously challenging to upregulate. We differentiated hPSCs into highly enriched HLF+ HOXA+ hematopoietic progenitors with near-stoichiometric efficiency by blocking formation of unwanted lineages at each differentiation step. hPSC-derived HLF+ HOXA+ hematopoietic progenitors could avail both basic research and cellular therapies.

Transcriptional regulation of genetic variants in the SLC40A1 promoter

Solute carrier 40A1 (SLC40A1) encodes ferroportin, which is the only known transmembrane protein that exports elemental iron from mammalian cells and is essential for iron homeostasis. Mutations in SLC40A1 are associated with iron-overload disorders. In addition to ferroportin diseases, SLC40A1 expression is downregulated in various cancer types. Despite the clinical significance of the SLC40A1 transporter, only a few studies have investigated genetic variants in SLC40A1. The present study was performed to identify genetic variations in the SLC40A1 promoter and functionally characterize each variant using in vitro assays. We investigated four haplotypes and five variants in the SLC40A1 promoter. We observed that haplotype 3 (H3) had significantly lower promoter activity than H1, whereas the activity of H4 was significantly higher than that of H1. Luciferase activity of H2 was comparable to that of H1. In addition, four variants of SLC40A1, c.-1355G>C, c.-662C>T, c.-98G>C, and c.-8C>G, showed significantly increased luciferase activity compared to the wild type (WT), whereas c.-750G>A showed significantly decreased luciferase activity compared to the WT. Three transcription factors, cAMP response element-binding protein-1 (CREB-1), chicken ovalbumin upstream promoter transcription factor 1, and hepatic leukemia factor (HLF), were predicted to bind to the promoter regions of SLC40A1 near c.-662C>T, c.-98G>C, and c.-8C>G, respectively. Among these, CREB-1 and HLF bound more strongly to the variant sequences than to the WT and functioned as activators of SLC40A1 transcription. Collectively, our findings indicate that the two SLC40A1 promoter haplotypes affect SLC40A1 transcription, which is regulated by CREB-1 and HLF.

Clonal hematopoiesis related TET2 loss-of-function impedes IL1β-mediated epigenetic reprogramming in hematopoietic stem and progenitor cells

Clonal hematopoiesis (CH) is defined as a single hematopoietic stem/progenitor cell (HSPC) gaining selective advantage over a broader range of HSPCs. When linked to somatic mutations in myeloid malignancy-associated genes, such as TET2-mediated clonal hematopoiesis of indeterminate potential or CHIP, it represents increased risk for hematological malignancies and cardiovascular disease. IL1β is elevated in patients with CHIP, however, its effect is not well understood. Here we show that IL1β promotes expansion of pro-inflammatory monocytes/macrophages, coinciding with a failure in the demethylation of lymphoid and erythroid lineage associated enhancers and transcription factor binding sites, in a mouse model of CHIP with hematopoietic-cell-specific deletion of Tet2. DNA-methylation is significantly lost in wild type HSPCs upon IL1β administration, which is resisted by Tet2-deficient HSPCs, and thus IL1β enhances the self-renewing ability of Tet2-deficient HSPCs by upregulating genes associated with self-renewal and by resisting demethylation of transcription factor binding sites related to terminal differentiation. Using aged mouse models and human progenitors, we demonstrate that targeting IL1 signaling could represent an early intervention strategy in preleukemic disorders. In summary, our results show that Tet2 is an important mediator of an IL1β-promoted epigenetic program to maintain the fine balance between self-renewal and lineage differentiation during hematopoiesis.

Neutral sphingomyelinase blockade enhances hematopoietic stem cell fitness through an integrated stress response

Hematopoietic stem and progenitor cell (HSPC) transplantation serves as a curative therapy for many benign and malignant hematopoietic disorders and as a platform for gene therapy. However, growing needs for ex vivo manipulation of HSPC-graft products are limited by barriers in maintaining critical self-renewal and quiescence properties. The role of sphingolipid metabolism in safeguarding these essential cellular properties has been recently recognized, but not yet widely explored. Here, we demonstrate that pharmacologic and genetic inhibition of neutral sphingomyelinase 2 (nSMase-2) leads to sustained improvements in long-term competitive transplantation efficiency after ex vivo culture. Mechanistically, nSMase-2 blockade activates a canonical integrated stress response (ISR) and promotes metabolic quiescence in human and murine HSPCs. These adaptations result in part from disruption in sphingolipid metabolism that impairs the release of nSMase-2-dependent extracellular vesicles (EVs). The aggregate findings link EV trafficking and the ISR as a regulatory dyad guarding HSPC homeostasis and long-term fitness. Translationally, transient nSMase-2 inhibition enables ex vivo graft manipulation with enhanced HSPC potency.

HLF promotes ovarian cancer progression and chemoresistance via regulating Hippo signaling pathway

Hepatic leukemia factor (HLF) is aberrantly expressed in human malignancies. However, the role of HLF in the regulation of ovarian cancer (OC) remains unknown. Herein, we reported that HLF expression was upregulated in OC tissues and ovarian cancer stem cells (CSCs). Functional studies have revealed that HLF regulates OC cell stemness, proliferation, and metastasis. Mechanistically, HLF transcriptionally activated Yes-associated protein 1 (YAP1) expression and subsequently modulated the Hippo signaling pathway. Moreover, we found that miR-520e directly targeted HLF 3'-UTR in OC cells. miR-520e expression was negatively correlated with HLF and YAP1 expression in OC tissues. The combined immunohistochemical (IHC) panels exhibited a better prognostic value for OC patients than any of these components alone. Importantly, the HLF/YAP1 axis determines the response of OC cells to carboplatin treatment and HLF depletion or the YAP1 inhibitor verteporfin abrogated carboplatin resistance. Analysis of patient-derived xenografts (PDXs) further suggested that HLF might predict carboplatin benefits in OC patients. In conclusion, these findings suggest a crucial role of the miR-520e/HLF/YAP1 axis in OC progression and chemoresistance, suggesting potential therapeutic targets for OC.

In vivo screening characterizes chromatin factor functions during normal and malignant hematopoiesis

Cellular differentiation requires extensive alterations in chromatin structure and function, which is elicited by the coordinated action of chromatin and transcription factors. By contrast with transcription factors, the roles of chromatin factors in differentiation have not been systematically characterized. Here, we combine bulk ex vivo and single-cell in vivo CRISPR screens to characterize the role of chromatin factor families in hematopoiesis. We uncover marked lineage specificities for 142 chromatin factors, revealing functional diversity among related chromatin factors (i.e. barrier-to-autointegration factor subcomplexes) as well as shared roles for unrelated repressive complexes that restrain excessive myeloid differentiation. Using epigenetic profiling, we identify functional interactions between lineage-determining transcription factors and several chromatin factors that explain their lineage dependencies. Studying chromatin factor functions in leukemia, we show that leukemia cells engage homeostatic chromatin factor functions to block differentiation, generating specific chromatin factor-transcription factor interactions that might be therapeutically targeted. Together, our work elucidates the lineage-determining properties of chromatin factors across normal and malignant hematopoiesis.

Quantification of hematopoietic stem and progenitor cells by targeted DNA methylation analysis

Hematopoietic stem and progenitor cells (HSPCs) are quantified in daily clinical practice by flow cytometry. In this study, we provide proof of concept that HSPCs can also be estimated by targeted DNA methylation (DNAm) analysis. The DNAm levels at three individual CG dinucleotides (CpG sites) in the genes MYO1D, STK17A, and SP140 correlated with CD34⁺ cell numbers in mobilized peripheral blood and with blast counts in leukemia. In the future, such epigenetic biomarkers can support the evaluation of stem cell mobilization, HSPC harvesting, or blast count in leukemia.

Temporal multimodal single-cell profiling of native hematopoiesis illuminates altered differentiation trajectories with age

Aging negatively affects hematopoiesis, with consequences for immunity and acquired blood cell disorders. Although impairments in hematopoietic stem cell (HSC) function contribute to this, the in vivo dynamics of such changes remain obscure. Here, we integrate extensive longitudinal functional assessments of HSC-specific lineage tracing with single-cell transcriptome and epitope profiling. In contrast to recent suggestions from single-cell RNA sequencing alone, our data favor a defined structure of HSC/progenitor differentiation that deviates substantially from HSC-derived hematopoiesis following transplantation. Native age-dependent attrition in HSC differentiation manifests as drastically reduced lymphoid output through an early lymphoid-primed progenitor (MPP Ly-I). While in vitro activation fails to rescue lymphoid differentiation from most aged HSCs, robust lymphopoiesis can be achieved by culturing elevated numbers of candidate HSCs. Therefore, our data position rare chronologically aged HSC clones, fully competent at producing lymphoid offspring, as a prime target for approaches aimed to improve lymphopoiesis in the elderly.

m6A RNA methylation-mediated upregulation of HLF promotes intrahepatic cholangiocarcinoma progression by regulating the FZD4/β-catenin signaling pathway

Hepatic leukemia factor (HLF) is aberrantly expressed in human malignancies. However, its role in regulating intrahepatic cholangiocarcinoma (ICC) remains unclear. This study aimed to define the role of HLF in ICC progression. Here, we showed that HLF expression is upregulated in ICC and predicts the poor prognosis in patients. Mechanistically, HLF activation in ICC is mediated by METTL3-dependent m6A methylation of the HLF mRNA. As per the results from the loss- or gain-of-function experiments, HLF promoted the self-renewal, tumorigenicity, proliferation and metastasis of ICC cells. RNA-seq and CUT&Tag analyses showed that frizzled-4 (FZD4) and forkhead box Q1 (FOXQ1) are target genes of HLF. Moreover, FOXQ1 transcriptionally activates METTL3 expression, forming a positive feedback loop, which subsequently activates WNT/β-catenin signaling and downstream tumor stemness. Furthermore, HLF expression was positively correlated with METTL3, IGF2BP3, FZD4 and FOXQ1 expression in ICC tissues in a large ICC cohort. The combined IHC panels exhibited a better prognostic value for patients with ICC than any of these components alone. In conclusions, these findings demonstrated that the METTL3/HLF/FOXQ1 regulatory circuit drove FZD4-mediated WNT/β-catenin activation in ICC progression, suggesting that targeting this axis could be novel therapeutic strategy for ICC.

Cell-Type-Specific Gene Regulatory Networks of Pro-Inflammatory and Pro-Resolving Lipid Mediator Biosynthesis in the Immune System

Lipid mediators are important regulators in inflammatory responses, and their biosynthetic pathways are targeted by commonly used anti-inflammatory drugs. Switching from pro-inflammatory lipid mediators (PIMs) to specialized pro-resolving (SPMs) is a critical step toward acute inflammation resolution and preventing chronic inflammation. Although the biosynthetic pathways and enzymes for PIMs and SPMs have now been largely identified, the actual transcriptional profiles underlying the immune cell type-specific transcriptional profiles of these mediators are still unknown. Using the Atlas of Inflammation Resolution, we created a large network of gene regulatory interactions linked to the biosynthesis of SPMs and PIMs. By mapping single-cell sequencing data, we identified cell type-specific gene regulatory networks of the lipid mediator biosynthesis. Using machine learning approaches combined with network features, we identified cell clusters of similar transcriptional regulation and demonstrated how specific immune cell activation affects PIM and SPM profiles. We found substantial differences in regulatory networks in related cells, accounting for network-based preprocessing in functional single-cell analyses. Our results not only provide further insight into the gene regulation of lipid mediators in the immune response but also shed light on the contribution of selected cell types in their biosynthesis.

A transcriptional constraint mechanism limits the homeostatic response to activity deprivation in mammalian neocortex

Healthy neuronal networks rely on homeostatic plasticity to maintain stable firing rates despite changing synaptic drive. These mechanisms, however, can themselves be destabilizing if activated inappropriately or excessively. For example, prolonged activity deprivation can lead to rebound hyperactivity and seizures. While many forms of homeostasis have been described, whether and how the magnitude of homeostatic plasticity is constrained remains unknown. Here, we uncover negative regulation of cortical network homeostasis by the PARbZIP family of transcription factors. In cortical slice cultures made from knockout mice lacking all three of these factors, the network response to prolonged activity withdrawal measured with calcium imaging is much stronger, while baseline activity is unchanged. Whole-cell recordings reveal an exaggerated increase in the frequency of miniature excitatory synaptic currents reflecting enhanced upregulation of recurrent excitatory synaptic transmission. Genetic analyses reveal that two of the factors, Hlf and Tef, are critical for constraining plasticity and for preventing life-threatening seizures. These data indicate that transcriptional activation is not only required for many forms of homeostatic plasticity but is also involved in restraint of the response to activity deprivation.

Nuclear factor I-C overexpression promotes monocytic development and cell survival in acute myeloid leukemia

Nuclear factor I-C (NFIC) belongs to a family of NFI transcription factors that binds to DNA through CAATT-boxes and are involved in cellular differentiation and stem cell maintenance. Here we show NFIC protein is significantly overexpressed in 69% of acute myeloid leukemia patients. Examination of the functional consequences of NFIC overexpression in HSPCs showed that this protein promoted monocytic differentiation. Single-cell RNA sequencing analysis further demonstrated that NFIC overexpressing monocytes had increased expression of growth and survival genes. In contrast, depletion of NFIC through shRNA decreased cell growth, increased cell cycle arrest and apoptosis in AML cell lines and AML patient blasts. Further, in AML cell lines (THP-1), bulk RNA sequencing of NFIC knockdown led to downregulation of genes involved in cell survival and oncogenic signaling pathways including mixed lineage leukemia-1 (MLL-1). Lastly, we show that NFIC knockdown in an ex vivo mouse MLL::AF9 pre-leukemic stem cell model, decreased their growth and colony formation and increased expression of myeloid differentiation markers Gr1 and Mac1. Collectively, our results suggest that NFIC is an important transcription factor in myeloid differentiation as well as AML cell survival and is a potential therapeutic target in AML.

Clonal Hematopoiesis: Role in Hematologic and Non-Hematologic Malignancies

Hematopoietic stem cells (HSCs) ensure the coordinated and balanced production of all hematopoietic cell types throughout life. Aging is associated with a gradual decline of the self-renewal and regenerative potential of HSCs and with the development of clonal hematopoiesis. Clonal hematopoiesis of indeterminate potential (CHIP) defines the clonal expansion of genetically variant hematopoietic cells bearing one or more gene mutations and/or structural variants (such as copy number alterations). CHIP increases exponentially with age and is associated with cancers, including hematologic neoplasia, cardiovascular and other diseases. The presence of CHIP consistently increases the risk of hematologic malignancy, particularly in individuals who have CHIP in association with peripheral blood cytopenia.

Hepatic leukemia factor-expressing paraxial mesoderm cells contribute to the developing brain vasculature

Recent genetic lineage tracing studies reveal heterogeneous origins of vascular endothelial cells and pericytes in the developing brain vasculature, despite classical experimental evidence for a mesodermal origin. Here we provide evidence through a genetic lineage tracing experiment that cephalic paraxial mesodermal cells give rise to endothelial cells and pericytes in the developing mouse brain. We show that Hepatic leukemia factor (Hlf) is transiently expressed by cephalic paraxial mesenchyme at embryonic day (E) 8.0-9.0 and the genetically-marked E8.0 Hlf-expressing cells mainly contribute to the developing brain vasculature. Interestingly, the genetically-marked E10.5 Hlf-expressing cells, which have been previously reported to contain embryonic hematopoietic stem cells, fail to contribute to the vascular cells. Combined, our genetic lineage tracing data demonstrate that a transient expression of Hlf marks a cephalic paraxial mesenchyme contributing to the developing brain vasculature.

Conservation of dynamic characteristics of transcriptional regulatory elements in periodic biological processes

Background

Cell and circadian cycles control a large fraction of cell and organismal physiology by regulating large periodic transcriptional programs that encompass anywhere from 15 to 80% of the genome despite performing distinct functions. In each case, these large periodic transcriptional programs are controlled by gene regulatory networks (GRNs), and it has been shown through genetics and chromosome mapping approaches in model systems that at the core of these GRNs are small sets of genes that drive the transcript dynamics of the GRNs. However, it is unlikely that we have identified all of these core genes, even in model organisms. Moreover, large periodic transcriptional programs controlling a variety of processes certainly exist in important non-model organisms where genetic approaches to identifying networks are expensive, time-consuming, or intractable. Ideally, the core network components could be identified using data-driven approaches on the transcriptome dynamics data already available.

Results

This study shows that a unified set of quantified dynamic features of high-throughput time series gene expression data are more prominent in the core transcriptional regulators of cell and circadian cycles than in their outputs, in multiple organism, even in the presence of external periodic stimuli. Additionally, we observe that the power to discriminate between core and non-core genes is largely insensitive to the particular choice of quantification of these features.

Conclusions

There are practical applications of the approach presented in this study for network inference, since the result is a ranking of genes that is enriched for core regulatory elements driving a periodic phenotype. In this way, the method provides a prioritization of follow-up genetic experiments. Furthermore, these findings reveal something unexpected—that there are shared dynamic features of the transcript abundance of core components of unrelated GRNs that control disparate periodic phenotypes.

Identification of a Four-Gene Signature Associated with the Prognosis Prediction of Lung Adenocarcinoma Based on Integrated Bioinformatics Analysis

Lung adenocarcinoma (LUAD) is often diagnosed at an advanced stage, so it is necessary to identify potential biomarkers for the early diagnosis and prognosis of LUAD. In our study, a gene co-expression network was constructed using weighted gene co-expression network analysis (WGCNA) in order to obtain the key modules and genes correlated with LUAD prognosis. Four hub genes (HLF, CHRDL1, SELENBP1, and TMEM163) were screened out using least absolute shrinkage and selection operator (LASSO)–Cox regression analysis; then, a prognostic model was established for predicting overall survival (OS) based on these four hub genes..Furthermore, the prognostic values of this four-gene signature were verified in four validation sets (GSE26939, GSE31210, GSE72094, and TCGA-LUAD) as well as in the GEPIA database. To assess the prognostic values of hub genes, receiver operating characteristic (ROC) curves were constructed and a nomogram was created. We found that a higher expression of four hub genes was associated with a lower risk of patient death. In a training set, it was demonstrated that this four-gene signature was a better prognostic factor than clinical factors such as age and stage of disease. Moreover, our results revealed that these four genes were suppressor factors of LUAD and that their high expression was associated with a lower risk of death. In summary, we demonstrated that this four-gene signature could be a potential prognostic factor for LUAD patients. These findings provide a theoretical basis for exploring potential biomarkers for LUAD prognosis prediction in the future.

HLF expression defines the human hematopoietic stem cell state

Hematopoietic stem cells (HSCs) sustain blood cell homeostasis throughout life and can regenerate all blood lineages after transplantation. Despite this clear functional definition, highly enriched isolation of human HSCs can currently only be achieved through combinatorial assessment of multiple surface antigens. Although several transgenic HSC reporter mouse strains have been described, no analogous approach to prospectively isolate human HSCs has been reported. To identify genes with the most selective expression in human HSCs, we profiled population and single-cell transcriptomes of unexpanded and ex vivo cultured cord blood-derived hematopoietic stem and progenitor cells as well as peripheral blood, adult bone marrow, and fetal liver. On the basis of these analyses, we propose the master transcription factor HLF (hepatic leukemia factor) as one of the most specific HSC marker genes. To directly track its expression in human hematopoietic cells, we developed a genomic HLF reporter strategy, capable of selectively labeling the most immature blood cells on the basis of a single engineered parameter. Most importantly, HLF-expressing cells comprise all stem cell activity in culture and in vivo during serial transplantation. Taken together, these results experimentally establish HLF as a defining gene of the human HSC state and outline a new approach to continuously mark these cells with high fidelity.

FLT3-ITD allelic ratio and HLF expression predict FLT3 inhibitor efficacy in adult AML

FLT3 internal tandem duplication (FLT3-ITD) is a frequent mutation in acute myeloid leukemia (AML) and remains a strong prognostic factor due to high rate of disease recurrence. Several FLT3-targeted agents have been developed, but determinants of variable responses to these agents remain understudied. Here, we investigated the role FLT3-ITD allelic ratio (ITD-AR), ITD length, and associated gene expression signatures on FLT3 inhibitor response in adult AML. We performed fragment analysis, ex vivo drug testing, and next generation sequencing (RNA, exome) to 119 samples from 87 AML patients and 13 healthy bone marrow controls. We found that ex vivo response to FLT3 inhibitors is significantly associated with ITD-AR, but not with ITD length. Interestingly, we found that the HLF gene is overexpressed in FLT3-ITD⁺ AML and associated with ITD-AR. The retrospective analysis of AML patients treated with FLT3 inhibitor sorafenib showed that patients with high HLF expression and ITD-AR had better clinical response to therapy compared to those with low ITD-AR and HLF expression. Thus, our findings suggest that FLT3 ITD-AR together with increased HLF expression play a role in variable FLT3 inhibitor responses observed in FLT3-ITD⁺ AML patients.

Characterization of a genomic region 8 kb downstream of GFI1B associated with myeloproliferative neoplasms

A genomic locus 8 kb downstream of the transcription factor GFI1B (Growth Factor Independence 1B) predisposes to clonal hematopoiesis and myeloproliferative neoplasms. One of the most significantly associated polymorphisms in this region is rs621940-G. GFI1B auto-represses GFI1B, and altered GFI1B expression contributes to myeloid neoplasms. We studied whether rs621940-G affects GFI1B expression and growth of immature cells. GFI1B ChIP-seq showed clear binding to the rs621940 locus. Preferential binding of various hematopoietic transcription factors to either the rs621940-C or -G allele was observed, but GFI1B showed no preference. In gene reporter assays the rs621940 region inhibited GFI1B promoter activity with the G-allele having less suppressive effects compared to the C-allele. However, CRISPR-Cas9 mediated deletion of the locus in K562 cells did not alter GFI1B expression nor auto-repression. In healthy peripheral blood mononuclear cells GFI1B expression did not differ consistently between the rs621940 alleles. Long range and targeted deep sequencing did not detect consistent effects of rs621940-G on allelic GFI1B expression either. Finally, we observed that myeloid colony formation was not significantly affected by either rs621940 allele in 193 healthy donors. Together, these findings show no evidence that rs621940 or its locus affect GFI1B expression, auto-repression or growth of immature myeloid cells.

Discrete regulatory modules instruct hematopoietic lineage commitment and differentiation

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. Here we perform dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit the stem and progenitor stage. Collectively, these data provide insights into the temporally regulated synergy of the cis- and the trans-regulatory components underlying hematopoietic lineage commitment and differentiation.

Hlf Expression Marks Early Emergence of Hematopoietic Stem Cell Precursors With Adult Repopulating Potential and Fate

In the aorta-gonad-mesonephros (AGM) region of mouse embryos, pre-hematopoietic stem cells (pre-HSCs) are generated from rare and specialized hemogenic endothelial cells (HECs) via endothelial-to-hematopoietic transition, followed by maturation into bona fide hematopoietic stem cells (HSCs). As HECs also generate a lot of hematopoietic progenitors not fated to HSCs, powerful tools that are pre-HSC/HSC-specific become urgently critical. Here, using the gene knockin strategy, we firstly developed an Hlf-tdTomato reporter mouse model and detected Hlf-tdTomato expression exclusively in the hematopoietic cells including part of the immunophenotypic CD45^– and CD45⁺ pre-HSCs in the embryonic day (E) 10.5 AGM region. By in vitro co-culture together with long-term transplantation assay stringent for HSC precursor identification, we further revealed that unlike the CD45^– counterpart in which both Hlf-tdTomato-positive and negative sub-populations harbored HSC competence, the CD45⁺ E10.5 pre-HSCs existed exclusively in Hlf-tdTomato-positive cells. The result indicates that the cells should gain the expression of Hlf prior to or together with CD45 to give rise to functional HSCs. Furthermore, we constructed a novel Hlf-CreER mouse model and performed time-restricted genetic lineage tracing by a single dose induction at E9.5. We observed the labeling in E11.5 AGM precursors and their contribution to the immunophenotypic HSCs in fetal liver (FL). Importantly, these Hlf-labeled early cells contributed to and retained the size of the HSC pool in the bone marrow (BM), which continuously differentiated to maintain a balanced and long-term multi-lineage hematopoiesis in the adult. Therefore, we provided another valuable mouse model to specifically trace the fate of emerging HSCs during development.

Single-cell ATAC-seq reveals GATA2-dependent priming defect in myeloid and a maturation bottleneck in lymphoid lineages

Germline heterozygous mutations in GATA2 are associated with a syndrome characterized by cytopenias, atypical infections, and increased risk of hematologic malignancies. Here, we generated a zebrafish mutant of gata2b that recapitulated the myelomonocytopenia and B-cell lymphopenia of GATA2 deficiency syndrome. Using single-cell assay for transposase accessible chromatin with sequencing of marrow cells, we showed that loss of gata2b led to contrasting alterations in chromosome accessibility in early myeloid and lymphoid progenitors, associated with defects in gene expression. Within the myeloid lineage in gata2b mutant zebrafish, we identified an attenuated myeloid differentiation with reduced transcriptional priming and skewing away from the monocytic program. In contrast, in early lymphoid progenitors, gata2b loss led to accumulation of B-lymphoid transcription factor accessibility coupled with increased expression of the B-cell lineage-specification program. However, gata2b mutant zebrafish had incomplete B-cell lymphopoiesis with loss of lineage-specific transcription factor accessibility in differentiating B cells, in the context of aberrantly reduced oxidative metabolic pathways. Our results establish that transcriptional events in early progenitors driven by Gata2 are required to complete normal differentiation.

Deciphering transcriptional and functional heterogeneity in hematopoiesis with single-cell genomics

PURPOSE OF REVIEW

Single-cell genomic approaches have uncovered cell fate biases and heterogeneity within hematopoietic subpopulations. However, standard single-cell transcriptomics suffers from high sampling noise, which particularly skews the distribution of lowly expressed genes, such as transcription factors (TFs). This might preclude the identification of rare transcripts that define cell identity and demarcate cell fate biases. Moreover, these studies need to go hand in hand with relevant functional assays to ensure that observed gene expression changes represent biologically meaningful alterations.

RECENT FINDINGS

Single-cell lineage tracing and functional validation studies have uncovered cell fate bias within transcriptionally distinct hematopoietic stem and progenitor subpopulations. Novel markers identified using these strategies have been proposed to prospectively isolate functionally distinct subpopulations, including long-term hematopoietic stem cells for ex vivo applications. Furthermore, the continuous nature of hematopoiesis has prompted the study of the relationship between stochastic transcriptional noise in hematopoietic TFs and cell fate determination.

SUMMARY

An understanding of the limitations of single-cell genomic approaches and follow-up functional assays is critical to discern the technical and biological contribution of noise in hematopoietic heterogeneity, to identify rare gene expression states, and to uncover functionally distinct subpopulations within hematopoiesis.

SUPPLEMENTARY VIDEO

http://links.lww.com/COH/A23.

Discovery and validation of an expression signature for recurrence prediction in high-risk diffuse-type gastric cancer

BACKGROUND

Diffuse type gastric cancer (DGC), represented by low sensitivity to chemotherapy and poor prognosis, is a heterogenous malignancy in which patient subsets exhibit diverse oncological risk-profiles. This study aimed to develop molecular biomarkers for robust prognostic risk-stratification and improve survival outcomes in patients with diffuse type gastric cancer (DGC).

METHODS

We undertook a systematic and comprehensive discovery and validation effort to identify recurrence prediction biomarkers by analyzing genome-wide transcriptomic profiling data from 157 patients with DGC, followed by their validation in 254 patients from 2 clinical cohorts.

RESULTS

Genome-wide transcriptomic profiling identified a 7-gene panel for robust prediction of recurrence in DGC patients (AUC = 0.91), which was successfully validated in an independent dataset (AUC = 0.86). Examination of 180 specimens from a training cohort allowed us to establish a gene-based risk prediction model (AUC = 0.78; 95% CI 0.71-0.84), which was subsequently validated in an independent cohort of 74 GC patients (AUC = 0.83; 95% CI 0.72-0.90). The Kaplan-Meier analyses exhibited a consistently superior performance of our risk-prediction model in the identification of high- and low-risk patient subgroups, which was significantly improved when we combined our gene signature with the tumor stage in both clinical cohorts (AUC of 0.83 in the training cohort and 0.89 in the validation cohort). Finally, for an easier clinical translation, we established a nomogram that robustly predicted prognosis in patients with DGC.

CONCLUSIONS

Our novel transcriptomic signature for risk-stratification and identification of high-risk patients with recurrence could serve as an important clinical decision-making tool in patients with DGC.

High Hepatic leukemia factor expression indicates a favorable survival in glioma patients

Hepatic leukemia factor (HLF) is an oncogenic transcript factor, but its role in gliomas is unclear.With the open-access data from the Cancer Genome Atlatls (TCGA), HLF expression was compared between normal and glioma tissues and its correlation to patient survival, age, gender, race, and tumor grade was analyzed. Multivariate Cox regression was adopted to explore the independent risk factors for patient survival. Survivals between high and low HLF expression, and high and low model predicted risk subgroups were compared. 1, 2, 3, and 5-year patient survival were predicted with the Cox regression model. Gene set enrichment analysis (GSEA) was performed to predict the potential function of HLF.Expression and clinical data of 5 normal brain samples and 655 glioma samples were obtained from TCGA. HLF expression was downregulated in gliomas than normal brain tissue (P = .007), and negatively related to patient age and advancing tumor grade (P < .001). HLF was a protective factor for patient survival (OR = 0.81, 95%CI 0.67-0.99, P = .035). Patients' survivals were poorer in low HLF expression subgroups and the Cox regression model predicted high-risk subgroups (P < .001). The accuracy of the model in predicting 1, 2, 3, and 5-year patient survival was 0.864, 0.895, 0.907, and 0.893, respectively. GSEA revealed HLF mainly took part in regulating tumor cell metabolism and cell cycle.HLF was downregulated in gliomas than normal tissue, negatively related to patient age and tumor grade, and was an independent protective factor for glioma patients.

Transcriptome Response of Atlantic Salmon (Salmo salar) to a New Piscine Orthomyxovirus

Pilchard orthomyxovirus (POMV) is an emerging pathogen of concern to the salmon industry in Australia. To explore the molecular events that underpin POMV infection, we challenged Atlantic salmon (Salmo salar) post-smolts in seawater via cohabitation. Tissue samples of the head kidney and liver were collected from moribund and surviving individuals and analyzed using transcriptome sequencing. Viral loads were higher in the head kidney compared to the liver, yet the liver presented more upregulated genes. Fish infected with POMV showed a strong innate immune response that included the upregulation of pathogen recognition receptors such as RIG-I and Toll-like receptors as well as the induction of interferon-stimulated genes (MX, ISG15). Moribund fish also presented a dramatic induction of pro-inflammatory cytokines, contributing to severe tissue damage and morbidity. An induction of major histocompatibility complex (MHC) class I genes (B2M) and markers of T cell-mediated immunity (CD8-alpha, CD8-beta, Perforin-1, Granzyme-A) was observed in both moribund fish and survivors. In addition, differential connectivity analysis showed that three key regulators (RELA/p65, PRDM1, and HLF) related to cell-mediated immunity had significant differences in connectivity in "clinically healthy" versus "clinically affected" or moribund fish. Collectively, our results show that T cell-mediated immunity plays a central role in the response of Atlantic salmon to the infection with POMV.

ARID1a Associates with Lymphoid-Restricted Transcription Factors and Has an Essential Role in T Cell Development

Maturation of lymphoid cells is controlled by the action of stage and lineage-restricted transcription factors working in concert with the general transcription and chromatin remodeling machinery to regulate gene expression. To better understand this functional interplay, we used Biotin Identification in human embryonic kidney cells to identify proximity interaction partners for GATA3, TCF7 (TCF1), SPI1, HLF, IKZF1, PAX5, ID1, and ID2. The proximity interaction partners shared among the lineage-restricted transcription factors included ARID1a, a BRG1-associated factor complex component. CUT&RUN analysis revealed that ARID1a shared binding with TCF7 and GATA3 at a substantial number of putative regulatory elements in mouse T cell progenitors. In support of an important function for ARID1a in lymphocyte development, deletion of Arid1a in early lymphoid progenitors in mice resulted in a pronounced developmental arrest in early T cell development with a reduction of CD4⁺CD8⁺ cells and a 20-fold reduction in thymic cellularity. Exploring gene expression patterns in DN3 cells from Wt and Arid1a-deficient mice suggested that the developmental block resided in the DN3a to DN3b transition, indicating a deficiency in β-selection. Our work highlights the critical importance of functional interactions between stage and lineage-restricted factors and the basic transcription machinery during lymphocyte differentiation.

Targeting the oncogenic activity of TCF3-HLF in leukemia

The oncogenic fusion transcription factor TCF3-HLF identifies an aggressive subtype of acute lymphoblastic leukemia. TCF3-HLF imposes a malignant program by activation of lineage-specific oncogenic enhancers. Among critical cofactors of the TCF3-HLF complex we identified EP300, which functional inhibition results in potent anti-leukemic activity by interference with the specific gene expression.

Downregulation of the circadian rhythm regulator HLF promotes multiple-organ distant metastases in non-small cell lung cancer through PPAR/NF-κb signaling

Non-small cell lung cancer (NSCLC) is the leading cause of cancer-related death due to its early recurrence and widespread metastatic potential. Accumulating studies have reported that dysregulation of circadian rhythms-associated regulators is implicated in the recurrence and metastasis of NSCLC. Therefore, identification of metastasis-associated circadian rhythm genes is clinically necessary. Here we report that the circadian gene hepatic leukemia factor (HLF), which was dramatically reduced in early-relapsed NSCLC tissues, was significantly correlated with early progression and distant metastasis in NSCLC patients. Upregulating HLF inhibited, while silencing HLF promoted lung colonization, as well as metastasis of NSCLC cells to bone, liver and brain in vivo. Importantly, downexpression of HLF promoted anaerobic metabolism to support anchorage-independent growth of NSCLC cells under low nutritional condition by activating NF-κB/p65 signaling through disrupting translocation of PPARα and PPARγ. Further investigations revealed that both genetic deletion and methylation contribute to downexpression of HLF in NSCLC tissues. In conclusion, our results shed light on a plausible mechanism by which HLF inhibits distant metastasis in NSCLC, suggesting that HLF may serve as a novel target for clinical intervention in NSCLC.

Bayesian integrative analysis of epigenomic and transcriptomic data identifies Alzheimer's disease candidate genes and networks

Biomedical research studies have generated large multi-omic datasets to study complex diseases like Alzheimer’s disease (AD). An important aim of these studies is the identification of candidate genes that demonstrate congruent disease-related alterations across the different data types measured by the study. We developed a new method to detect such candidate genes in large multi-omic case-control studies that measure multiple data types in the same set of samples. The method is based on a gene-centric integrative coefficient quantifying to what degree consistent differences are observed in the different data types. For statistical inference, a Bayesian hierarchical model is used to study the distribution of the integrative coefficient. The model employs a conditional autoregressive prior to integrate a functional gene network and to share information between genes known to be functionally related. We applied the method to an AD dataset consisting of histone acetylation, DNA methylation, and RNA transcription data from human cortical tissue samples of 233 subjects, and we detected 816 genes with consistent differences between persons with AD and controls. The findings were validated in protein data and in RNA transcription data from two independent AD studies. Finally, we found three subnetworks of jointly dysregulated genes within the functional gene network which capture three distinct biological processes: myeloid cell differentiation, protein phosphorylation and synaptic signaling. Further investigation of the myeloid network indicated an upregulation of this network in early stages of AD prior to accumulation of hyperphosphorylated tau and suggested that increased CSF1 transcription in astrocytes may contribute to microglial activation in AD. Thus, we developed a method that integrates multiple data types and external knowledge of gene function to detect candidate genes, applied the method to an AD dataset, and identified several disease-related genes and processes demonstrating the usefulness of the integrative approach.

Recent technological advances have led to a new generation of studies that interrogate multiple molecular levels in the same target tissue of a set of subjects, generating complex multi-omic datasets with which to study disease mechanism. These datasets of genetic, epigenomic, transcriptomic, and other data have the potential to reveal novel biological insights; however, integrative analyses remain challenging and require new computational methods. We developed an integrative Bayesian approach to detect genes with consistent differences between case and control samples across multiple data types. The method further integrates prior knowledge about gene function in the form of a gene functional similarity network to improve statistical inference by sharing information between related genes. We applied our method to an Alzheimer’s disease dataset of epigenomic and transcriptomic data and detected and then validated several novel and known candidate genes as well as three major disease-related biological processes. One of these processes reflected microglial activation and included the cytokine CSF1. Single-nucleus data revealed that CSF1 was primarily upregulated in astrocytes, implicating the involvement of this cell type in microglial activation. Hence, we demonstrated that integrative analysis approaches to multi-omic datasets can improve candidate gene detection and thereby generate new insights into complex diseases.

The Leukemogenic TCF3-HLF Complex Rewires Enhancers Driving Cellular Identity and Self-Renewal Conferring EP300 Vulnerability

The chimeric transcription factor TCF3-HLF defines an incurable acute lymphoblastic leukemia subtype. Here we decipher the regulome of endogenous TCF3-HLF and dissect its essential transcriptional components and targets by functional genomics. We demonstrate that TCF3-HLF recruits HLF binding sites at hematopoietic stem cell/myeloid lineage associated (super-) enhancers to drive lineage identity and self-renewal. Among direct targets, hijacking an HLF binding site in a MYC enhancer cluster by TCF3-HLF activates a conserved MYC-driven transformation program crucial for leukemia propagation in vivo. TCF3-HLF pioneers the cooperation with ERG and recruits histone acetyltransferase p300 (EP300), conferring susceptibility to EP300 inhibition. Our study provides a framework for targeting driving transcriptional dependencies in this fatal leukemia.

Oncofetal HLF transactivates c-Jun to promote hepatocellular carcinoma development and sorafenib resistance

BACKGROUND AND AIMS

The unique expression pattern makes oncofetal proteins ideal diagnostic biomarkers and therapeutic targets in cancer. However, few oncofetal proteins have been identified and entered clinical practice.

METHODS

Fetal liver, adult liver and hepatocellular carcinoma (HCC) tissues were employed to assess the expression of hepatic leukaemia factor (HLF). The impact of HLF on HCC onset and progression was investigated both in vivo and in vitro. The association between HLF and patient prognosis was determined in patient cohorts. The correlation between HLF expression and sorafenib benefits in HCC was further evaluated in patient cohorts and patient-derived xenografts (PDXs).

RESULTS

HLF is a novel oncofetal protein which is reactivated in HCC by SOX2 and OCT4. Functional studies revealed that HLF transactivates c-Jun to promote tumour initiating cell (TIC) generation and enhances TIC-like properties of hepatoma cells, thus driving HCC initiation and progression. Consistently, our clinical investigations elucidated the association between HLF and patient prognosis and unravelled the close correlation between HLF levels and c-Jun expression in patient HCCs. Importantly, HLF/c-Jun axis determines the responses of hepatoma cells to sorafenib treatment, and interference of HLF abrogated c-Jun activation and enhanced sorafenib response. Analysis of patient cohorts and PDXs further suggests that HLF/c-Jun axis might serve as a biomarker for sorafenib benefits in HCC patients.

CONCLUSIONS

Our findings uncovered HLF as a novel oncofetal protein and revealed the crucial role of the HLF/c-Jun axis in HCC development and sorafenib response, rendering HLF as an optimal target for the prevention and intervention of HCC.

Transcription factor Oct1 protects against hematopoietic stress and promotes acute myeloid leukemia

A better understanding of the development and progression of acute myelogenous leukemia (AML) is necessary to improve patient outcome. Here we define roles for the transcription factor Oct1/Pou2f1 in AML and normal hematopoiesis. Inappropriate reactivation of the CDX2 gene is widely observed in leukemia patients and in leukemia mouse models. We show that Oct1 associates with the CDX2 promoter in both normal and AML primary patient samples, but recruits the histone demethylase Jmjd1a/Kdm3a to remove the repressive H3K9me2 mark only in malignant specimens. The CpG DNA immediately adjacent to the Oct1 binding site within the CDX2 promoter exhibits variable DNA methylation in healthy control blood and bone marrow samples, but complete demethylation in AML samples. In MLL-AF9-driven mouse models, partial loss of Oct1 protects from myeloid leukemia. Complete Oct1 loss completely suppresses leukemia but results in lethality from bone marrow failure. Loss of Oct1 in normal hematopoietic transplants results in superficially normal long-term reconstitution; however, animals become acutely sensitive to 5-fluorouracil, indicating that Oct1 is dispensable for normal hematopoiesis but protects blood progenitor cells against external chemotoxic stress. These findings elucidate a novel and important role for Oct1 in AML.

Mutant CEBPA directly drives the expression of the targetable tumor-promoting factor CD73 in AML

The key myeloid transcription factor (TF), CEBPA, is frequently mutated in acute myeloid leukemia (AML), but the direct molecular effects of this leukemic driver mutation remain elusive. To investigate CEBPA mutant AML, we performed microscale, in vivo chromatin immunoprecipitation sequencing and identified a set of aberrantly activated enhancers, exclusively occupied by the leukemia-associated CEBPA-p30 isoform. Comparing gene expression changes in human CEBPA mutant AML and the corresponding Cebpa ^Lp30 mouse model, we identified Nt5e, encoding CD73, as a cross-species AML gene with an upstream leukemic enhancer physically and functionally linked to the gene. Increased expression of CD73, mediated by the CEBPA-p30 isoform, sustained leukemic growth via the CD73/A2AR axis. Notably, targeting of this pathway enhanced survival of AML-transplanted mice. Our data thus indicate a first-in-class link between a cancer driver mutation in a TF and a druggable, direct transcriptional target.

Hlf marks the developmental pathway for hematopoietic stem cells but not for erythro-myeloid progenitors

Hematopoietic stem cells (HSCs) and HSC-independent progenitors are generated from hemogenic endothelium. Yokomizo et al. show that Hlf expression distinguishes nascent HSCs from HSC-independent progenitors. HSC specification, regulated by the Evi-1/Hlf axis, is activated only within Hlf⁺ nascent hematopoietic clusters.

Before the emergence of hematopoietic stem cells (HSCs), lineage-restricted progenitors, such as erythro-myeloid progenitors (EMPs), are detected in the embryo or in pluripotent stem cell cultures in vitro. Although both HSCs and EMPs are derived from hemogenic endothelium, it remains unclear how and when these two developmental programs are segregated during ontogeny. Here, we show that hepatic leukemia factor (Hlf) expression specifically marks a developmental continuum between HSC precursors and HSCs. Using the Hlf-tdTomato reporter mouse, we found that Hlf is expressed in intra-aortic hematopoietic clusters and fetal liver HSCs. In contrast, EMPs and yolk sac hematopoietic clusters before embryonic day 9.5 do not express Hlf. HSC specification, regulated by the Evi-1/Hlf axis, is activated only within Hlf⁺ nascent hematopoietic clusters. These results strongly suggest that HSCs and EMPs are generated from distinct cohorts of hemogenic endothelium. Selective induction of the Hlf⁺ lineage pathway may lead to the in vitro generation of HSCs from pluripotent stem cells.

Hepatic leukemia factor is a novel leukemic stem cell regulator in DNMT3A, NPM1, and FLT3-ITD triple-mutated AML

FLT3, DNMT3A, and NPM1 are the most frequently mutated genes in cytogenetically normal acute myeloid leukemia (AML), but little is known about how these mutations synergize upon cooccurrence. Here we show that triple-mutated AML is characterized by high leukemia stem cell (LSC) frequency, an aberrant leukemia-specific GPR56 ^highCD34^low immunophenotype, and synergistic upregulation of Hepatic Leukemia Factor (HLF). Cell sorting based on the LSC marker GPR56 allowed isolation of triple-mutated from DNMT3A/NPM1 double-mutated subclones. Moreover, in DNMT3A R882-mutated patients, CpG hypomethylation at the HLF transcription start site correlated with high HLF mRNA expression, which was itself associated with poor survival. Loss of HLF via CRISPR/Cas9 significantly reduced the CD34⁺GPR56⁺ LSC compartment of primary human triple-mutated AML cells in serial xenotransplantation assays. HLF knockout cells were more actively cycling when freshly harvested from mice, but rapidly exhausted when reintroduced in culture. RNA sequencing of primary human triple-mutated AML cells after shRNA-mediated HLF knockdown revealed the NOTCH target Hairy and Enhancer of Split 1 (HES1) and the cyclin-dependent kinase inhibitor CDKN1C/p57 as novel targets of HLF, potentially mediating these effects. Overall, our data establish HLF as a novel LSC regulator in this genetically defined high-risk AML subgroup.

Hepatic Leukemia Factor Maintains Quiescence of Hematopoietic Stem Cells and Protects the Stem Cell Pool during Regeneration

The transcription factor hepatic leukemia factor (HLF) is strongly expressed in hematopoietic stem cells (HSCs) and is thought to influence both HSC self-renewal and leukemogenesis. However, the physiological role of HLF in hematopoiesis and HSC function is unclear. Here, we report that mice lacking Hlf are viable with essentially normal hematopoietic parameters, including an intact HSC pool during steady-state hematopoiesis. In contrast, when challenged through transplantation, Hlf-deficient HSCs showed an impaired ability to reconstitute hematopoiesis and became gradually exhausted upon serial transplantation. Transcriptional profiling of Hlf-deficient HSCs revealed changes associated with enhanced cellular activation, and cell-cycle analysis demonstrated a significant reduction of quiescent HSCs. Accordingly, toxic insults targeting dividing cells completely eradicated the HSC pool in Hlf-deficient mice. In summary, our findings point to HLF as a critical regulator of HSC quiescence and as an essential factor for maintaining the HSC pool during regeneration.