Hematopoietic stem cells (HSCs) are important in cell-based therapies for blood-related disorders. While progress has been made in generating HSCs by directed differentiation of pluripotent stem cells (PSCs), such cultures promote hematopoietic progenitor cells (HPCs) over HSCs. Thus, elucidating markers and factors associated with HSC versus HPC development is imperative. HSCs and HPCs originate from hemogenic endothelium (HE) in a Runx1-dependent process. Here, we characterize a Runx1 enhancer (+110) that distinguishes emerging dorsal aorta HSCs from HPCs. Comparative transcriptomics reveal a 17-gene signature associated with in vivo long-term HSC potential, while single-cell multiome analysis demonstrates a clear epigenetic identity of dorsal aorta HE, preHSCs, and HPCs, providing a resource of stage-specific enhancer activity. Our study demonstrates the power of cell type-specific enhancer-reporter models to dissect cell fate decisions in development and provides new inroads to label and/or perturb HSC versus HPC potential in vivo and in vitro.
Background The development and diversification of sensory proprioceptive neurons, which reside in the dorsal root ganglia (DRG) and express the tropomyosin receptor kinase C (TrkC), depend on the transcription factor (TF) Runx3. Runx3-deficient mice develop severe limb ataxia due to TrkC neuron cell death. Two additional TFs Pou4f1 (also called Brn3a) and Isl1 also play an important role in sensory neuron development. Thus, we aimed to unravel the chromatin state of early-developing TrkC neurons and decipher the Runx3 high-confidence target genes (HCT) and the possible cooperation between Runx3, Brn3a and Isl1 in the regulation of these genes. Methods Runx3 expression is driven by the gene proximal P2 promoter. Transcriptome analysis was conducted by RNA-seq on RNA isolated from heterozygous (P2+/-) vs. homozygous (P2-/-) TrkC neurons and differentially expressed genes (DEGs) were determined. Genome-wide occupancy of Runx3, Brn3a, Isl1 and histone H3 acetylated on lysine 27 (H3K27Ac) was determined using CUT&RUN. The landscape of Transposase-accessible chromatin was analyzed via ATAC-seq. Findings The intersection of Runx3 genomic occupancy-associated genes and DEG data discovered 244 Runx3 HCT. Brn3a and Isl1 were found to bind to numerous genomic loci, some of which overlapped with Runx3. Most genomic regions bound by each of these three TFs or co-bound by them resided in distantly located enhancer regions rather than in gene promoters. In activated and suppressed neuronal Runx3 HCT, Runx3 cooperated mainly with Brn3a to regulate expression through distantly located enhancers. Interestingly, suppression of non-neuronal immune genes was mainly managed via Runx3 without Brn3a. The distribution of ATAC and H3K27Ac marked regions in Runx3 peaks containing at least one RUNX binding site (Runx3_RBS) revealed that while most promoter regions were marked by ATAC, a prominent fraction of intron/intergenic regions occupied by Runx3, Brn3a or Isl1 were unmarked by ATAC and/or H3K27Ac. Conclusions These analyses shed new light on the interplay of Runx3, Brn3a, Isl1, and open chromatin regions in regulating the Runx3 HCT in the early developmental stages of TrkC neurons.
Abstract The development and diversification of proprioceptive neurons, which reside in the dorsal root ganglia (DRG) and express the tropomyosin receptor kinase C (TrkC), depend on the transcription factor (TF) Runx3. As one of the major functions of TrkC neurons is to coordinate limb movements, Runx3-deficient mice develop severe limb ataxia due to TrkC neuron cell death. In TrkC neurons Runx3 expression is driven by the gene proximal P2 promoter. Accordingly, P2-/- mice develop severe ataxia. Heterozygous (P2+/-) vs. homozygous (P2-/-) TrkC neuron transcriptome highlighted a repertoire of differentially expressed genes (DEG). The Intersection of Runx3 genomic occupancy-associated genes and DEG data discovered 244 high-confidence Runx3 target genes (HCT). We also analyzed the genomic occupancy of two additional TFs, Brn3a and Isl1, that play a role in sensory neuron development and determined the landscape of histone H3 acetylated lysine 27 (H3K27Ac) and transpoase-accessible chromatin with sequencing (ATAC-seq)-marked chromatin regions in E11.5 TrkC neurons. These analyses shed new light on the interplay of Runx3, Brn3a, Isl1, and open chromatin regions in regulating the HCT genes in the early developmental stages of TrkC neurons.
Runt domain-related (Runx) transcription factors are essential for early T cell development in mice from uncommitted to committed stages. Single and double Runx knockouts via Cas9 show that target genes responding to Runx activity are not solely controlled by the dominant factor, Runx1. Instead, Runx1 and Runx3 are coexpressed in single cells; bind to highly overlapping genomic sites; and have redundant, collaborative functions regulating genes pivotal for T cell development. Despite stable combined expression levels across pro-T cell development, Runx1 and Runx3 preferentially activate and repress genes that change expression dynamically during lineage commitment, mostly activating T-lineage genes and repressing multipotent progenitor genes. Furthermore, most Runx target genes are sensitive to Runx perturbation only at one stage and often respond to Runx more for expression transitions than for maintenance. Contributing to this highly stage-dependent gene regulation function, Runx1 and Runx3 extensively shift their binding sites during commitment. Functionally distinct Runx occupancy sites associated with stage-specific activation or repression are also distinguished by different patterns of partner factor cobinding. Finally, Runx occupancies change coordinately at numerous clustered sites around positively or negatively regulated targets during commitment. This multisite binding behavior may contribute to a developmental "ratchet" mechanism making commitment irreversible.
The zinc finger transcription factor, Bcl11b, is expressed in T cells and group 2 innate lymphoid cells (ILC2s) among hematopoietic cells. In early T-lineage cells, Bcl11b directly binds and represses the gene encoding the E protein antagonist, Id2, preventing pro-T cells from adopting innate-like fates. In contrast, ILC2s co-express both Bcl11b and Id2. To address this contradiction, we have directly compared Bcl11b action mechanisms in pro-T cells and ILC2s. We found that Bcl11b binding to regions across the genome shows distinct cell type-specific motif preferences. Bcl11b occupies functionally different sites in lineage-specific patterns and controls totally different sets of target genes in these cell types. In addition, Bcl11b bears cell type-specific post-translational modifications and organizes different cell type-specific protein complexes. However, both cell types use the same distal enhancer region to control timing of Bcl11b activation. Therefore, although pro-T cells and ILC2s both need Bcl11b for optimal development and function, Bcl11b works substantially differently in these two cell types.
Mice deficient in the transcription factor Runx3 develop a multitude of immune system defects, including early onset colitis. This paper demonstrates that Runx3 is expressed in colonic mononuclear phagocytes (MNP), including resident macrophages (RM) and dendritic cell subsets (cDC2). Runx3 deletion in MNP causes early onset colitis due to their impaired maturation. Mechanistically, the resulting MNP subset imbalance leads to up-regulation of pro-inflammatory genes as occurs in IL10R-deficient RM. In addition, RM and cDC2 display a marked decrease in expression of anti-inflammatory/TGF β-regulated genes and β-catenin signaling associated genes, respectively. MNP transcriptome and ChIP-seq data analysis suggest that a significant fraction of genes affected by Runx3 loss are direct Runx3 targets. Collectively, Runx3 imposes intestinal immune tolerance by regulating maturation of colonic anti-inflammatory MNP, befitting the identification of RUNX3 as a genome-wide associated risk gene for various immune-related diseases in humans, including gastrointestinal tract diseases such as Crohn's disease and celiac.
Patients with neurofibromatosis type 1 (NF1) are predisposed to develop neurofibromas, but the underlying molecular mechanisms of neurofibromagenesis are not fully understood. We showed dual genetic deletion of Runx1 and Runx3 in Schwann cells (SCs) and SC precursors delayed neurofibromagenesis and prolonged mouse survival. We identified peripheral myelin protein 22 (Pmp22/Gas3) related to neurofibroma initiation. Knockdown of Pmp22 with short hairpin RNAs increased Runx1fl/fl;Runx3fl/fl;Nf1fl/fl;DhhCre tumor-derived sphere numbers and enabled significantly more neurofibroma-like microlesions on transplantation. Conversely, overexpression of Pmp22 in mouse neurofibroma SCs decreased cell proliferation. Mechanistically, RUNX1/3 regulated alternative promoter usage and induced levels of protein expression of Pmp22 to control SC growth. Last, pharmacological inhibition of RUNX/core-binding factor β (CBFB) activity significantly reduced neurofibroma volume in vivo. Thus, we identified a signaling pathway involving RUNX1/3 suppression of Pmp22 in neurofibroma initiation and/or maintenance. Targeting disruption of RUNX/CBFB interaction might provide a novel therapy for patients with neurofibroma.
Maintaining posture requires tight regulation of the position and orientation of numerous spinal components. Yet, surprisingly little is known about this regulatory mechanism, whose failure may result in spinal deformity as in adolescent idiopathic scoliosis. Here, we use genetic mouse models to demonstrate the involvement of proprioception in regulating spine alignment. Null mutants for Runx3 transcription factor, which lack TrkC neurons connecting between proprioceptive mechanoreceptors and spinal cord, developed peripubertal scoliosis not preceded by vertebral dysplasia or muscle asymmetry. Deletion of Runx3 in the peripheral nervous system or specifically in peripheral sensory neurons, or of enhancer elements driving Runx3 expression in proprioceptive neurons, induced a similar phenotype. Egr3 knockout mice, lacking muscle spindles, but not Golgi tendon organs, displayed a less severe phenotype, suggesting that both receptor types may be required for this regulatory mechanism. These findings uncover a central role for the proprioceptive system in maintaining spinal alignment.
The Runx3 transcription factor is essential for development and diversification of the dorsal root ganglia (DRGs) TrkC sensory neurons. In Runx3-deficient mice, developing TrkC neurons fail to extend central and peripheral afferents, leading to cell death and disruption of the stretch reflex circuit, resulting in severe limb ataxia. Despite its central role, the mechanisms underlying the spatiotemporal expression specificities of Runx3 in TrkC neurons were largely unknown. Here we first defined the genomic transcription unit encompassing regulatory elements (REs) that mediate the tissue-specific expression of Runx3. Using transgenic mice expressing BAC reporters spanning the Runx3 locus, we discovered three REs—dubbed R1, R2, and R3—that cross-talk with promoter-2 (P2) to drive TrkC neuron-specific Runx3 transcription. Deletion of single or multiple elements either in the BAC transgenics or by CRISPR/Cas9-mediated endogenous ablation established the REs’ ability to promote and/or repress Runx3 expression in developing sensory neurons. Our analysis reveals that an intricate combinatorial interplay among the three REs governs Runx3 expression in distinct subtypes of TrkC neurons while concomitantly extinguishing its expression in non-TrkC neurons. These findings provide insights into the mechanism regulating cell type-specific expression and subtype diversification of TrkC neurons in developing DRGs.
In this chapter we summarize the pros and cons of the notion that Runx3 is a major tumor suppressor gene (TSG). Inactivation of TSGs in normal cells provides a viability/growth advantage that contributes cell-autonomously to cancer. More than a decade ago it was suggested that RUNX3 is involved in gastric cancer development, a postulate extended later to other epithelial cancers portraying RUNX3 as a major TSG. However, evidence that Runx3 is not expressed in normal gastric and other epithelia has challenged the RUNX3-TSG paradigm. In contrast, RUNX3 is overexpressed in a significant fraction of tumor cells in various human epithelial cancers and its overexpression in pancreatic cancer cells promotes their migration, anchorage-independent growth and metastatic potential. Moreover, recent high-throughput quantitative genome-wide studies on thousands of human samples of various tumors and new investigations of the role of Runx3 in mouse cancer models have unequivocally demonstrated that RUNX3 is not a bona fide cell-autonomous TSG. Importantly, accumulating data demonstrated that RUNX3 functions in control of immunity and inflammation, thereby indirectly influencing epithelial tumor development.
The transcriptional control of lineage commitment to various ILC subsets is incompletely understood. Yokoyama and colleagues show that Runx3 is essential for the normal development of ILC1 and ILC3 cells but not ILC2 cells. Subsets of innate lymphoid cells (ILCs) reside in the mucosa and regulate immune responses to external pathogens. While ILCs can be phenotypically classified into ILC1, ILC2 and ILC3 subsets, the transcriptional control of commitment to each ILC lineage is incompletely understood. Here we report that the transcription factor Runx3 was essential for the normal development of ILC1 and ILC3 cells but not of ILC2 cells. Runx3 controlled the survival of ILC1 cells but not of ILC3 cells. Runx3 was required for expression of the transcription factor RORγt and its downstream target, the transcription factor AHR, in ILC3 cells. The absence of Runx3 in ILCs exacerbated infection with Citrobacter rodentium. Therefore, our data establish Runx3 as a key transcription factor in the lineage-specific differentiation of ILC1 and ILC3 cells.
Natural killer cells belong to the family of innate lymphoid cells comprising the frontline defense against infected and transformed cells. Development and activation of natural killer cells is highly dependent on interleukin-15 signaling. However, very little is known about the transcription program driving this process. The transcription factor Runx3 is highly expressed in natural killer cells, but its function in these cells is largely unknown. We show that loss of Runx3 impaired interleukin-15-dependent accumulation of mature natural killer cells in vivo and under culture conditions and pregnant Runx3(-/-) mice completely lack the unique population of interleukin-15-dependent uterine natural killer cells. Combined chromatin immunoprecipitation sequencing and differential gene expression analysis of wild-type versus Runx3-deficient in vivo activated splenic natural killer cells revealed that Runx3 cooperates with ETS and T-box transcription factors to drive the interleukin-15-mediated transcription program during activation of these cells. Runx3 functions as a nuclear regulator during interleukin-15-dependent activation of natural killer cells by regulating the expression of genes involved in proliferation, maturation, and migration. Similar studies with additional transcription factors will allow the construction of a more detailed transcriptional network that controls natural killer cell development and function.
The transcription factor Runx3 is highly expressed in CD8+ T and NK cytotoxic lymphocytes and is required for their effective activation and proliferation but molecular insights into the transcription program regulated by Runx3 in these cells are still missing. Using Runx3-ChIP-seq and transcriptome analysis of wild type vs. Runx3-/- primary cells we have now identified Runx3-regulated genes in the two cell types at both resting and IL-2-activated states. Runx3-bound genomic regions in both cell types were distantly located relative to gene transcription start sites and were enriched for RUNX and ETS motifs. Bound genomic regions significantly overlapped T-bet and p300-bound enhancer regions in Runx3-expressing Th1 helper cells. Compared to resting cells, IL-2-activated CD8+ T and NK cells contain three times more Runx3-regulated genes that are common to both cell types. Functional annotation of shared CD8+ T and NK Runx3-regulated genes revealed enrichment for immune-associated terms including lymphocyte activation, proliferation, cytotoxicity, migration and cytokine production, highlighting the role of Runx3 in CD8+ T and NK activated cells.
Gastric cancer (GC) is a major cause of cancer mortality. GC studies that aim to identify relevant oncogenes and tumor suppressor genes (TSGs) are essential for devising effective new therapies. A decade ago, RUNX3, a gene that resides on human chromosome 1p36.1, was claimed to be a major TSG in GC. Since then, hundreds of studies involving thousands of GC patients have attempted to verify and extend the RUNX3 TSG paradigm. However, RUNX3 is not recognized as TSG and not listed in the “Cancer Gene Census” website. To be a TSG that protects normal cells against malignancy, the gene must be expressed in the normal tissue from which the cancer arose and its loss or inactivation should contribute to cancer development. This review summarizes compelling body of evidence challenging the RUNX3-TSG paradigm. Studies show unequivocally that RUNX3 is not expressed in normal gastric epithelium and that it fails to fulfill all other premises of a TSG. RUNX3 mutations and 1p36 deletions are not frequent in GC and RUNX3 is not associated with familial GC or with increased risk of GC. Accordingly, Runx3-/- mice do not develop tumors. RUNX3 promoter methylation, which has been reported to be a frequent event in GC, is not relevant to its alleged TSG function, since the gene is already silent in normal gastric epithelium. In sharp contrast, overexpression of RUNX3 was found in several types of human cancers, including GC, and the 1p36.1 region is amplified in B-cell lymphoma. Thus, it is possible that RUNX3 actually promotes cancer development rather than being a TSG. The true targets for GC therapy are discussed below. Those are genes frequently lost or amplified in GC and are well known for their tumor suppressive or oncogenic activity, respectively.
Runx transcription factor family proteins have essential roles during T-cell development by either activating or repressing target genes. For instance, lineage- and stage-specific expression of Cd4 and ThPOK is controlled by a transcriptional silencer embedded in each locus, whose activity requires bindings of Runx complexes. The evolutionarily conserved VWRPY penta-peptide sequences in Runx proteins have been shown to be responsible for repressive function as a platform to recruit Groucho/TLE transcriptional corepressors. However, it remains elusive whether requirement for the VWRPY motif differs among Runx target genes. By examining mice lacking VWRPY motifs in both Runx1 and Runx3 proteins, here, we show a full and partial derepression of Cd4 and ThPOK in CD8-linegae T cells, respectively. Thus, whereas Cd4 silencing completely depends on the VWRPY motif, both VWRPY-dependent and -independent mechanisms operate to repress ThPOK gene. These results indicate that Runx proteins utilize different modes to repress expression of different target genes.
The formation of functional connectivity in the nervous system is governed by axon guidance that instructs nerve growth and branching during development, implying a similarity between neuronal subtypes in terms of nerve extension. We demonstrate the molecular mechanism of another layer of complexity in vertebrates by defining a transcriptional program underlying growth differences between positionally different neurons. The rate of axon extension of the early subset of embryonic dorsal root ganglion sensory neurons is encoded in neurons at different axial levels. This code is determined by a segmental pattern of axial levels of Runx family transcription factor Runx3. Runx3 in turn determines transcription levels of genes encoding cytoskeletal proteins involved in axon extension, including Rock1 and Rock2 which have ongoing activities determining axon growth in early sensory neurons and blocking Rock activity reverses axon extension deficits of Runx3(-/-) neurons. Thus, Runx3 acts to regulate positional differences in axon extension properties apparently without affecting nerve guidance and branching, a principle that could be relevant to other parts of the nervous system.
Runx transcription factor family proteins have essential roles during T-cell development by either activating or repressing target genes. For instance, lineage- and stage-specific expression of Cd4 and ThPOK is controlled by a transcriptional silencer embedded in each locus, whose activity requires bindings of Runx complexes. The evolutionarily conserved VWRPY penta-peptide sequences in Runx proteins have been shown to be responsible for repressive function as a platform to recruit Groucho/TLE transcriptional corepressors. However, it remains elusive whether requirement for the VWRPY motif differs among Runx target genes. By examining mice lacking VWRPY motifs in both Runx1 and Runx3 proteins, here, we show a full and partial derepression of Cd4 and ThPOK in CD8-linegae T cells, respectively. Thus, whereas Cd4 silencing completely depends on the VWRPY motif, both VWRPY-dependent and -independent mechanisms operate to repress ThPOK gene. These results indicate that Runx proteins utilize different modes to repress expression of different target genes.
The Runx3 transcription factor regulates cell fate decisions during embryonic development and in adults. It was previously reported that Runx3 is strongly expressed in embryonic and adult gastrointestinal tract (GIT) epithelium (Ep) and that its loss causes gastric cancer. More than 280 publications have based their research on these findings and concluded that Runx3 is indeed a tumour suppressor (TS). In stark contrast, using various measures, we found that Runx3 expression is undetectable in GIT Ep. Employing a variety of biochemical and genetic techniques, including analysis of Runx3‐GFP and R26LacZ/Runx3Cre or R26tdTomato/Runx3Cre reporter strains, we readily detected Runx3 in GIT‐embedded leukocytes, dorsal root ganglia, skeletal elements and hair follicles. However, none of these approaches revealed detectable Runx3 levels in GIT Ep. Moreover, our analysis of the original Runx3LacZ/LacZ mice used in the previously reported study failed to reproduce the GIT expression of Runx3. The lack of evidence for Runx3 expression in normal GIT Ep creates a serious challenge to the published data and undermines the notion that Runx3 is a TS involved in cancer pathogenesis.