Multilineage transcriptional priming and determination of alternate hematopoietic cell fates

Egr2 drives the differentiation of Ly6Chi monocytes into fibrosis-promoting macrophages in metabolic dysfunction-associated steatohepatitis in mice

Metabolic dysfunction-associated steatohepatitis (MASH), previously called non-alcoholic steatohepatitis (NASH), is a growing concern worldwide, with liver fibrosis being a critical determinant of its prognosis. Monocyte-derived macrophages have been implicated in MASH-associated liver fibrosis, yet their precise roles and the underlying differentiation mechanisms remain elusive. In this study, we unveil a key orchestrator of this process: long chain saturated fatty acid-Egr2 pathway. Our findings identify the transcription factor Egr2 as the driving force behind monocyte differentiation into hepatic lipid-associated macrophages (hLAMs) within MASH liver. Notably, Egr2-deficiency reroutes monocyte differentiation towards a macrophage subset resembling resident Kupffer cells, hampering hLAM formation. This shift has a profound impact, suppressing the transition from benign steatosis to liver fibrosis, demonstrating the critical pro-fibrotic role played by hLAMs in MASH pathogenesis. Long-chain saturated fatty acids that accumulate in MASH liver emerge as potent inducers of Egr2 expression in macrophages, a process counteracted by unsaturated fatty acids. Furthermore, oral oleic acid administration effectively reduces hLAMs in MASH mice. In conclusion, our work not only elucidates the intricate interplay between saturated fatty acids, Egr2, and monocyte-derived macrophages but also highlights the therapeutic promise of targeting the saturated fatty acid-Egr2 axis in monocytes for MASH management.

The Hill function is the universal Hopfield barrier for sharpness of input-output responses

The Hill functions, [Formula: see text], have been widely used in biology for over a century but, with the exception of [Formula: see text], they have had no justification other than as a convenient fit to empirical data. Here, we show that they are the universal limit for the sharpness of any input-output response arising from a Markov process model at thermodynamic equilibrium. Models may represent arbitrary molecular complexity, with multiple ligands, internal states, conformations, coregulators, etc, under core assumptions that are detailed in the paper. The model output may be any linear combination of steady-state probabilities, with components other than the chosen input ligand held constant. This formulation generalizes most of the responses in the literature. We use a coarse-graining method in the graph-theoretic linear framework to show that two sharpness measures for input-output responses fall within an effectively bounded region of the positive quadrant, [Formula: see text], for any equilibrium model with [Formula: see text] input binding sites. [Formula: see text] exhibits a cusp which approaches, but never exceeds, the sharpness of [Formula: see text], but the region and the cusp can be exceeded when models are taken away from thermodynamic equilibrium. Such fundamental thermodynamic limits are called Hopfield barriers, and our results provide a biophysical justification for the Hill functions as the universal Hopfield barriers for sharpness. Our results also introduce an object, [Formula: see text], whose structure may be of mathematical interest, and suggest the importance of characterizing Hopfield barriers for other forms of cellular information processing.

Deletion of TECRL promotes skeletal muscle repair by up-regulating EGR2

Myogenic regeneration relies on the proliferation and differentiation of satellite cells. TECRL (trans-2,3-enoyl-CoA reductase like) is an endoplasmic reticulum protein only expressed in cardiac and skeletal muscle. However, its role in myogenesis remains unknown. We show that TECRL expression is increased in response to injury. Satellite cell-specific deletion of TECRL enhances muscle repair by increasing the expression of EGR2 through the activation of the ERK1/2 signaling pathway, which in turn promotes the expression of PAX7. We further show that TECRL deletion led to the upregulation of the histone acetyltransferase general control nonderepressible 5, which enhances the transcription of EGR2 through acetylation. Importantly, we showed that AAV9-mediated TECRL silencing improved muscle repair in mice. These findings shed light on myogenic regeneration and muscle repair.

Accelerated plasma-cell differentiation in Bach2-deficient mouse B cells is caused by altered IRF4 functions

Transcription factors BACH2 and IRF4 are both essential for antibody class-switch recombination (CSR) in activated B lymphocytes, while they oppositely regulate the differentiation of plasma cells (PCs). Here, we investigated how BACH2 and IRF4 interact during CSR and plasma-cell differentiation. We found that BACH2 organizes heterochromatin formation of target gene loci in mouse splenic B cells, including targets of IRF4 activation such as Aicda, an inducer of CSR, and Prdm1, a master plasma-cell regulator. Release of these gene loci from heterochromatin in response to B-cell receptor stimulation was coupled to AKT-mTOR pathway activation. In Bach2-deficient B cells, PC genes' activation depended on IRF4 protein accumulation, without an increase in Irf4 mRNA. Mechanistically, a PU.1-IRF4 heterodimer in activated B cells promoted BACH2 function by inducing gene expression of Bach2 and Pten, a negative regulator of AKT signaling. Elevated AKT activity in Bach2-deficient B cells resulted in IRF4 protein accumulation. Thus, BACH2 and IRF4 mutually modulate the activity of each other, and BACH2 inhibits PC differentiation by both the repression of PC genes and the restriction of IRF4 protein accumulation.

Network switches and their role in circadian clocks

Circadian rhythms are generated by complex interactions among genes and proteins. Self-sustained ∼24 h oscillations require negative feedback loops and sufficiently strong nonlinearities that are the product of molecular and network switches. Here, we review common mechanisms to obtain switch-like behavior, including cooperativity, antagonistic enzymes, multisite phosphorylation, positive feedback, and sequestration. We discuss how network switches play a crucial role as essential components in cellular circadian clocks, serving as integral parts of transcription-translation feedback loops that form the basis of circadian rhythm generation. The design principles of network switches and circadian clocks are illustrated by representative mathematical models that include bistable systems and negative feedback loops combined with Hill functions. This work underscores the importance of negative feedback loops and network switches as essential design principles for biological oscillations, emphasizing how an understanding of theoretical concepts can provide insights into the mechanisms generating biological rhythms.

Gfi1 controls the formation of effector CD8 T cells during chronic infection and cancer

During chronic infections and tumor progression, CD8 T cells gradually lose their effector functions and become exhausted. These exhausted CD8 T cells are heterogeneous and comprised of different subsets, including self-renewing progenitors that give rise to Ly108 - CX3CR1 + effector-like cells. Generation of these effector-like cells is essential for the control of chronic infections and tumors, albeit limited. However, the precise cues and mechanisms directing the formation and maintenance of exhausted effector-like are incompletely understood. Using genetic mouse models challenged with LCMV Clone 13 or syngeneic tumors, we show that the expression of a transcriptional repressor, growth factor independent 1 (Gfi1) is dynamically regulated in exhausted CD8 T cells, which in turn regulates the formation of exhausted effector-like cells. Gfi1 deletion in T cells dysregulates the chromatin accessibility and transcriptomic programs associated with the differentiation of LCMV Clone 13-specific CD8 T cell exhaustion, preventing the formation of effector-like and terminally exhausted cells while maintaining progenitors and a newly identified Ly108 + CX3CR1 + state. These Ly108 + CX3CR1 + cells have a distinct chromatin profile and may represent an alternative target for therapeutic interventions to combat chronic infections and cancer. In sum, we show that Gfi1 is a critical regulator of the formation of exhausted effector-like cells.

The Hill function is the universal Hopfield barrier for sharpness of input-output responses

The Hill functions, ℋ h ( x ) = x h / 1 + x h , have been widely used in biology for over a century but, with the exception of ℋ 1 , they have had no justification other than as a convenient fit to empirical data. Here, we show that they are the universal limit for the sharpness of any input-output response arising from a Markov process model at thermodynamic equilibrium. Models may represent arbitrary molecular complexity, with multiple ligands, internal states, conformations, co-regulators, etc, under core assumptions that are detailed in the paper. The model output may be any linear combination of steady-state probabilities, with components other than the chosen input ligand held constant. This formulation generalises most of the responses in the literature. We use a coarse-graining method in the graph-theoretic linear framework to show that two sharpness measures for input-output responses fall within an effectively bounded region of the positive quadrant, Ω m ⊂ ℝ + 2 , for any equilibrium model with m input binding sites. Ω m exhibits a cusp which approaches, but never exceeds, the sharpness of ℋ m but the region and the cusp can be exceeded when models are taken away from thermodynamic equilibrium. Such fundamental thermodynamic limits are called Hopfield barriers and our results provide a biophysical justification for the Hill functions as the universal Hopfield barriers for sharpness. Our results also introduce an object, Ω m , whose structure may be of mathematical interest, and suggest the importance of characterising Hopfield barriers for other forms of cellular information processing.

Canalizing cell fate by transcriptional repression

Precision in the establishment and maintenance of cellular identities is crucial for the development of multicellular organisms and requires tight regulation of gene expression. While extensive research has focused on understanding cell type-specific gene activation, the complex mechanisms underlying the transcriptional repression of alternative fates are not fully understood. Here, we provide an overview of the repressive mechanisms involved in cell fate regulation. We discuss the molecular machinery responsible for suppressing alternative fates and highlight the crucial role of sequence-specific transcription factors (TFs) in this process. Depletion of these TFs can result in unwanted gene expression and increased cellular plasticity. We suggest that these TFs recruit cell type-specific repressive complexes to their cis-regulatory elements, enabling them to modulate chromatin accessibility in a context-dependent manner. This modulation effectively suppresses master regulators of alternative fate programs and their downstream targets. The modularity and dynamic behavior of these repressive complexes enables a limited number of repressors to canalize and maintain major and minor cell fate decisions at different stages of development.

Regularly spaced tyrosines in EBF1 mediate BRG1 recruitment and formation of nuclear subdiffractive clusters

B lineage priming by pioneer transcription factor EBF1 requires the function of an intrinsically disordered region (IDR). Here, we examine the role of regularly spaced tyrosines in the IDR as potential determinants of IDR function and activity of EBF1. We found that four Y > A mutations in EBF1 reduced the formation of condensates in vitro and subdiffractive clusters in vivo. Notably, Y > A mutant EBF1 was inefficient in promoting B cell differentiation and showed impaired chromatin binding, recruitment of BRG1, and activation of specific target genes. Thus, regularly spaced tyrosines in the IDR contribute to the biophysical and functional properties of EBF1.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

During gene regulation, DNA accessibility is thought to limit the availability of transcription factor (TF) binding sites, while TFs can increase DNA accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events in the modulation of gene expression remain unknown for the vast majority of genes. We utilized deeply sequenced ATAC-Seq data and site-specific knock-in reporter genes to investigate the relationship between the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of Cebpa during macrophage-neutrophil differentiation. While the enhancers upregulate reporter expression during the earliest stages of differentiation, there is little corresponding increase in their total accessibility. Conversely, total accessibility peaks during the last stages of differentiation without any increase in enhancer activity. The accessibility of positions neighboring C/EBP-family TF binding sites, which indicates TF occupancy, does increase significantly during early differentiation, showing that the early upregulation of enhancer activity is driven by TF binding. These results imply that a generalized increase in DNA accessibility is not sufficient, and binding by enhancer-specific TFs is necessary, for the upregulation of gene expression. Additionally, high-coverage ATAC-Seq combined with time-series expression data can infer the sequence of regulatory events at binding-site resolution.

Evolution of hierarchy and irreversibility in theoretical cell differentiation model

The process of cell differentiation in multicellular organisms is characterized by hierarchy and irreversibility in many cases. However, the conditions and selection pressures that give rise to these characteristics remain poorly understood. By using a mathematical model, here we show that the network of differentiation potency (differentiation diagram) becomes necessarily hierarchical and irreversible by increasing the number of terminally differentiated states under certain conditions. The mechanisms generating these characteristics are clarified using geometry in the cell state space. The results demonstrate that the hierarchical organization and irreversibility can manifest independently of direct selection pressures associated with these characteristics, instead they appear to evolve as byproducts of selective forces favoring a diversity of differentiated cell types. The study also provides a new perspective on the structure of gene regulatory networks that produce hierarchical and irreversible differentiation diagrams. These results indicate some constraints on cell differentiation, which are expected to provide a starting point for theoretical discussion of the implicit limits and directions of evolution in multicellular organisms.

A review on gene regulatory network reconstruction algorithms based on single cell RNA sequencing

BACKGROUND

Understanding gene regulatory networks (GRNs) is essential for unraveling the molecular mechanisms governing cellular behavior. With the advent of high-throughput transcriptome measurement technology, researchers have aimed to reverse engineer the biological systems, extracting gene regulatory rules from their outputs, which represented by gene expression data. Bulk RNA sequencing, a widely used method for measuring gene expression, has been employed for GRN reconstruction. However, it falls short in capturing dynamic changes in gene expression at the level of individual cells since it averages gene expression across mixed cell populations.

OBJECTIVE

In this review, we provide an overview of 15 GRN reconstruction tools and discuss their respective strengths and limitations, particularly in the context of single cell RNA sequencing (scRNA-seq).

METHODS

Recent advancements in scRNA-seq break new ground of GRN reconstruction. They offer snapshots of the individual cell transcriptomes and capturing dynamic changes. We emphasize how these technological breakthroughs have enhanced GRN reconstruction.

CONCLUSION

GRN reconstructors can be classified based on their requirement for cellular trajectory, which represents a dynamical cellular process including differentiation, aging, or disease progression. Benchmarking studies support the superiority of GRN reconstructors that do not require trajectory analysis in identifying regulator-target relationships. However, methods equipped with trajectory analysis demonstrate better performance in identifying key regulatory factors. In conclusion, researchers should select a suitable GRN reconstructor based on their specific research objectives.

Interpretable trajectory inference with single-cell Linear Adaptive Negative-binomial Expression (scLANE) testing

The rapid proliferation of trajectory inference methods for single-cell RNA-seq data has allowed researchers to investigate complex biological processes by examining underlying gene expression dynamics. After estimating a latent cell ordering, statistical models are used to determine which genes exhibit changes in expression that are significantly associated with progression through the biological trajectory. While a few techniques for performing trajectory differential expression exist, most rely on the flexibility of generalized additive models in order to account for the inherent nonlinearity of changes in gene expression. As such, the results can be difficult to interpret, and biological conclusions often rest on subjective visual inspections of the most dynamic genes. To address this challenge, we propose scLANE testing, which is built around an interpretable generalized linear model and handles nonlinearity with basis splines chosen empirically for each gene. In addition, extensions to estimating equations and mixed models allow for reliable trajectory testing under complex experimental designs. After validating the accuracy of scLANE under several different simulation scenarios, we apply it to a set of diverse biological datasets and display its ability to provide novel biological information when used downstream of both pseudotime and RNA velocity estimation methods.

Growth Factor Independence-1 (GFI-1) Gene Expression in Hematopoietic Stem Cell Lineage Differentiation in Low Birth Weight Newborns Compared With Normal Birth Weight Newborns at Term Pregnancy

Introduction Low birth weight (LBW), which is a risk factor for noncommunicable diseases throughout life, is a significant public health concern. In addition to regulating myeloid cell differentiation and proliferation, a transcriptional repressor identified as growth factor independence-1 (GFI-1) is essential for hematopoietic stem cell maintenance and self-renewal. The current study was designed to compare the expression of the GFI-1 gene in the differentiation of hematopoietic stem cells in newborns with LBW and those with normal birth weight (NBW). Methods A prospective comparative analytical study was carried out from September 2019 to September 2021 after obtaining Institute Ethical Committee approval at a tertiary care center in north India. The GFI-1 gene expression levels in 50 cord blood samples from women with term gestation and LBW newborns (<2500 grams) were measured using quantitative real-time polymerase chain reaction (RT-PCR) and compared to gene expression levels in 50 cord blood samples from women with term gestation and NBW newborns (≥2500 grams). The data were analyzed using IBM SPSS statistics software version 24.0 (IBM Corp., Armonk, NY). Results The median GFI-1 expression in LBW newborns is 3.1, whereas among NBW newborns it is 9.39. The difference is significant (P <0.001). The level of GFI-1 gene expression in LBW newborns was correlated with their birth weight. The coefficient of correlation was found to be weakly positive (r = 0.223). The birth weight of NBW newborns was correlated to the level of expression of the GFI-1 gene, which was found to be positively correlated (r = 0.332). Conclusion The levels of the GFI-1 gene and newborn birth weight were compared in LBW infants, which were weakly positively correlated. The level of GFI-1 gene expression at birth was compared to the birth weight of NBW newborns, which was positively correlated.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

LucFlow: A method to measure Luciferase reporter expression in single cells

Reporter assays, in which the expression of an inert protein is driven by gene regulatory elements such as promoters and enhancers, are a workhorse for investigating gene regulation. Techniques for measuring reporter gene expression vary from single-cell or single-molecule approaches having low throughput to bulk Luciferase assays that have high throughput. We developed a Luciferase Reporter Assay using Flow-Cytometry (LucFlow), which measures reporter expression in single cells immunostained for Luciferase. We optimized and tested LucFlow with a murine cell line that can be differentiated into neutrophils, into which promoter-reporter and enhancer-promoter-reporter constructs have been integrated in a site-specific manner. The single-cell measurements are comparable to bulk ones but we found that dead cells have no detectable Luciferase protein, so that bulk assays underestimate reporter expression. LucFlow is able to achieve a higher accuracy than bulk methods by excluding dead cells during flow cytometry. Prior to fixation and staining, the samples are spiked with stained cells that can be discriminated during flow cytometry and control for tube-to-tube variation in experimental conditions. Computing fold change relative to control cells allows LucFlow to achieve a high level of precision. LucFlow, therefore, enables the accurate and precise measurement of reporter expression in a high throughput manner.

mTORC2 orchestrates monocytic and granulocytic lineage commitment by an ATF5-mediated pathway

Myeloid hematopoiesis is a finely controlled consecutive developmental process, which is essential to maintain peripheral innate immune homeostasis. Herein, we found that Rictor deletion caused the remarkable reduction of granulocyte-monocyte progenitors (GMPs), monocytes, and macrophages, while the levels of neutrophils were unaffected. Adoptive transfer of Rictor-deleted GMPs or common myeloid progenitors (CMPs) in syngeneic mice showed poor re-constitution of monocytes compared to wild-type GMPs or CMPs. In addition to decreasing the proliferation of CMPs/GMPs, Rictor deletion preferentially inhibited Ly6C+ monocyte differentiation, while enhancing neutrophil differentiation, as determined by colony formation assays. mTORC2 promotes monocyte development by downregulation of the AKT-Foxo4-activating transcription factor 5 (ATF5)-mitochondrial unfolded protein response (mtUPR) pathway. Genetic overexpression of ATF5 or exposure to ethidium bromide significantly rescued monocyte/macrophage differentiation defects of Rictor-deficient myeloid progenitors. Therefore, Rictor is required for CMP/GMP proliferation and acts as an important switch to balance monocytic and granulocytic lineage commitment in bone marrow.

The miR-142 miRNAs: Shaping the naïve immune system

Immunity in a naïve organism is tightly controlled. Adequate proportions of the many immune cell subsets must be produced to mount efficient responses to eventual challenges. In addition, a functioning immune system is highly dynamic at steady state. Mature immune cells must be positioned properly and/or circulate to facilitate the detection of dangers. They must also be poised to promptly react to unusual encounters, while ignoring innocuous germs and self. Numerous regulatory mechanisms act at the molecular level to generate such an exquisite structure, including miRNA-mediated repression of protein synthesis. Notably, the miRNAs from the miR-142 locus are preferentially expressed in hematopoietic cells. Their importance is underscored by the deeply disturbed immune system seen upon inactivation of the locus in mice. In this review, we explore reported roles for the miR-142 miRNAs in the shaping of immunity in vertebrates, discussing in particular their contributions to the generation, migration and survival of hematopoietic cells.

Transcription factor networks link B-lymphocyte development and malignant transformation in leukemia

Rapid advances in genomics have opened unprecedented possibilities to explore the mutational landscapes in malignant diseases, such as B-cell acute lymphoblastic leukemia (B-ALL). This disease is manifested as a severe defect in the production of normal blood cells due to the uncontrolled expansion of transformed B-lymphocyte progenitors in the bone marrow. Even though classical genetics identified translocations of transcription factor-coding genes in B-ALL, the extent of the targeting of regulatory networks in malignant transformation was not evident until the emergence of large-scale genomic analyses. There is now evidence that many B-ALL cases present with mutations in genes that encode transcription factors with critical roles in normal B-lymphocyte development. These include PAX5, IKZF1, EBF1, and TCF3, all of which are targeted by translocations or, more commonly, partial inactivation in cases of B-ALL. Even though there is support for the notion that germline polymorphisms in the PAX5 and IKZF1 genes predispose for B-ALL, the majority of leukemias present with somatic mutations in transcription factor-encoding genes. These genetic aberrations are often found in combination with mutations in genes that encode components of the pre-B-cell receptor or the IL-7/TSLP signaling pathways, all of which are important for early B-cell development. This review provides an overview of our current understanding of the molecular interplay that occurs between transcription factors and signaling events during normal and malignant B-lymphocyte development.

A bistable autoregulatory module in the developing embryo commits cells to binary expression fates

Bistable autoactivation has been proposed as a mechanism for cells to adopt binary fates during embryonic development. However, it is unclear whether the autoactivating modules found within developmental gene regulatory networks are bistable, unless their parameters are quantitatively determined. Here, we combine in vivo live imaging with mathematical modeling to dissect the binary cell fate dynamics of the fruit fly pair-rule gene fushi tarazu (ftz), which is regulated by two known enhancers: the early (non-autoregulating) element and the autoregulatory element. Live imaging of transcription and protein concentration in the blastoderm revealed that binary Ftz fates are achieved as Ftz expression rapidly transitions from being dictated by the early element to the autoregulatory element. Moreover, we discovered that Ftz concentration alone is insufficient to activate the autoregulatory element, and that this element only becomes responsive to Ftz at a prescribed developmental time. Based on these observations, we developed a dynamical systems model and quantitated its kinetic parameters directly from experimental measurements. Our model demonstrated that the ftz autoregulatory module is indeed bistable and that the early element transiently establishes the content of the binary cell fate decision to which the autoregulatory module then commits. Further in silico analysis revealed that the autoregulatory element locks the Ftz fate quickly, within 35 min of exposure to the transient signal of the early element. Overall, our work confirms the widely held hypothesis that autoregulation can establish developmental fates through bistability and, most importantly, provides a framework for the quantitative dissection of cellular decision-making.

Unveiling the complexity of transcription factor networks in hematopoietic stem cells: implications for cell therapy and hematological malignancies

The functionality and longevity of hematopoietic tissue is ensured by a tightly controlled balance between self-renewal, quiescence, and differentiation of hematopoietic stem cells (HSCs) into the many different blood lineages. Cell fate determination in HSCs is influenced by signals from extrinsic factors (e.g., cytokines, irradiation, reactive oxygen species, O2 concentration) that are translated and integrated by intrinsic factors such as Transcription Factors (TFs) to establish specific gene regulatory programs. TFs also play a central role in the establishment and/or maintenance of hematological malignancies, highlighting the need to understand their functions in multiple contexts. TFs bind to specific DNA sequences and interact with each other to form transcriptional complexes that directly or indirectly control the expression of multiple genes. Over the past decades, significant research efforts have unraveled molecular programs that control HSC function. This, in turn, led to the identification of more than 50 TF proteins that influence HSC fate. However, much remains to be learned about how these proteins interact to form molecular networks in combination with cofactors (e.g. epigenetics factors) and how they control differentiation, expansion, and maintenance of cellular identity. Understanding these processes is critical for future applications particularly in the field of cell therapy, as this would allow for manipulation of cell fate and induction of expansion, differentiation, or reprogramming of HSCs using specific cocktails of TFs. Here, we review recent findings that have unraveled the complexity of molecular networks controlled by TFs in HSCs and point towards possible applications to obtain functional HSCs ex vivo for therapeutic purposes including hematological malignancies. Furthermore, we discuss the challenges and prospects for the derivation and expansion of functional adult HSCs in the near future.

Ratiometric sensing of Pnt and Yan transcription factor levels confers ultrasensitivity to photoreceptor fate transitions in Drosophila

Cell state transitions are often triggered by large changes in the concentrations of transcription factors and therefore large differences in their stoichiometric ratios. Whether cells can elicit transitions using modest changes in the ratios of co-expressed factors is unclear. Here, we investigate how cells in the Drosophila eye resolve state transitions by quantifying the expression dynamics of the ETS transcription factors Pnt and Yan. Eye progenitor cells maintain a relatively constant ratio of Pnt/Yan protein, despite expressing both proteins with pulsatile dynamics. A rapid and sustained twofold increase in the Pnt/Yan ratio accompanies transitions to photoreceptor fates. Genetic perturbations that modestly disrupt the Pnt/Yan ratio produce fate transition defects consistent with the hypothesis that transitions are normally driven by a twofold shift in the ratio. A biophysical model based on cooperative Yan-DNA binding coupled with non-cooperative Pnt-DNA binding illustrates how twofold ratio changes could generate ultrasensitive changes in target gene transcription to drive fate transitions. Thus, coupling cell state transitions to the Pnt/Yan ratio sensitizes the system to modest fold-changes, conferring robustness and ultrasensitivity to the developmental program.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

The upregulation of gene expression by enhancers depends upon the interplay between the binding of sequence-specific transcription factors (TFs) and DNA accessibility. DNA accessibility is thought to limit the ability of TFs to bind to their sites, while TFs can increase accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events underlying the modulation of gene expression during cellular differentiation remain unknown for the vast majority of genes. We investigated the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of an important neutrophil gene, Cebpa, during macrophage-neutrophil differentiation. Reporter genes were integrated in a site-specific manner in PUER cells, which are progenitors that can be differentiated into neutrophils or macrophages in vitro by activating the pan-leukocyte TF PU.1. Time series data show that two enhancers upregulate reporter expression during the first 48 hours of neutrophil differentiation. Surprisingly, there is little or no increase in the total accessibility, measured by ATAC-Seq, of the enhancers during the same time period. Conversely, total accessibility peaks 96 hrs after PU.1 activation-consistent with its role as a pioneer-but the enhancers do not upregulate gene expression. Combining deeply sequenced ATAC-Seq data with a new bias-correction method allowed the profiling of accessibility at single-nucleotide resolution and revealed protected regions in the enhancers that match all previously characterized TF binding sites and ChIP-Seq data. Although the accessibility of most positions does not change during early differentiation, that of positions neighboring TF binding sites, an indicator of TF occupancy, did increase significantly. The localized accessibility changes are limited to nucleotides neighboring C/EBP-family TF binding sites, showing that the upregulation of enhancer activity during early differentiation is driven by C/EBP-family TF binding. These results show that increasing the total accessibility of enhancers is not sufficient for upregulating their activity and other events such as TF binding are necessary for upregulation. Also, TF binding can cause upregulation without a perceptible increase in total accessibility. Finally, this study demonstrates the feasibility of comprehensively mapping individual TF binding sites as footprints using high coverage ATAC-Seq and inferring the sequence of events in gene regulation by combining with time-series gene expression data.

In Vitro Human Haematopoietic Stem Cell Expansion and Differentiation

The haematopoietic system plays an essential role in our health and survival. It is comprised of a range of mature blood and immune cell types, including oxygen-carrying erythrocytes, platelet-producing megakaryocytes and infection-fighting myeloid and lymphoid cells. Self-renewing multipotent haematopoietic stem cells (HSCs) and a range of intermediate haematopoietic progenitor cell types differentiate into these mature cell types to continuously support haematopoietic system homeostasis throughout life. This process of haematopoiesis is tightly regulated in vivo and primarily takes place in the bone marrow. Over the years, a range of in vitro culture systems have been developed, either to expand haematopoietic stem and progenitor cells or to differentiate them into the various haematopoietic lineages, based on the use of recombinant cytokines, co-culture systems and/or small molecules. These approaches provide important tractable models to study human haematopoiesis in vitro. Additionally, haematopoietic cell culture systems are being developed and clinical tested as a source of cell products for transplantation and transfusion medicine. This review discusses the in vitro culture protocols for human HSC expansion and differentiation, and summarises the key factors involved in these biological processes.

Airway and parenchyma transcriptomics in a house dust mite model of experimental asthma

Lung transcriptomics studies in asthma have provided valuable information in the whole lung context, however, deciphering the individual contributions of the airway and parenchyma in disease pathogenesis may expedite the development of novel targeted treatment strategies. In this study, we performed transcriptomics on the airway and parenchyma using a house dust mite (HDM)-induced model of experimental asthma that replicates key features of the human disease. HDM exposure increased the expression of 3,255 genes, of which 212 were uniquely increased in the airways, 856 uniquely increased in the parenchyma, and 2187 commonly increased in both compartments. Further interrogation of these genes using a combination of network and transcription factor enrichment analyses identified several transcription factors that regulate airway and/or parenchymal gene expression, including transcription factor EC (TFEC), transcription factor PU.1 (SPI1), H2.0-like homeobox (HLX), metal response element binding transcription factor-1 (MTF1) and E74-like factor 4 (ets domain transcription factor, ELF4) involved in controlling innate immune responses. We next assessed the effects of inhibiting lung SPI1 responses using commercially available DB1976 and DB2313 on key disease outcomes. We found that both compounds had no protective effects on airway inflammation, however DB2313 (8 mg/kg) decreased mucus secreting cell number, and both DB2313 (1 mg/kg) and DB1976 (2.5 mg/kg and 1 mg/kg) reduced small airway collagen deposition. Significantly, both compounds decreased airway hyperresponsiveness. This study demonstrates that SPI1 is important in HDM-induced experimental asthma and that its pharmacological inhibition reduces HDM-induced airway collagen deposition and hyperresponsiveness.

The ENL YEATS epigenetic reader domain critically links MLL-ENL to leukemic stem cell frequency in t(11;19) Leukemia

MLL (KMT2a) translocations are found in ~10% of acute leukemia patients, giving rise to oncogenic MLL-fusion proteins. A common MLL translocation partner is ENL and associated with a poor prognosis in t(11;19) patients. ENL contains a highly conserved N-terminal YEATS domain that binds acetylated histones and interacts with the PAF1c, an epigenetic regulator protein complex essential for MLL-fusion leukemogenesis. Recently, wild-type ENL, and specifically the YEATS domain, was shown to be essential for leukemic cell growth. However, the inclusion and importance of the YEATS domain in MLL-ENL-mediated leukemogenesis remains unexplored. We found the YEATS domain is retained in 84.1% of MLL-ENL patients and crucial for MLL-ENL-mediated leukemogenesis in mouse models. Mechanistically, deletion of the YEATS domain impaired MLL-ENL fusion protein binding and decreased expression of pro-leukemic genes like Eya1 and Meis1. Point mutations that disrupt YEATS domain binding to acetylated histones decreased stem cell frequency and increased MLL-ENL-mediated leukemia latency. Therapeutically, YEATS containing MLL-ENL leukemic cells display increased sensitivity to the YEATS inhibitor SGC-iMLLT compared to control AML cells. Our results demonstrate that the YEATS domain is important for MLL-ENL fusion protein-mediated leukemogenesis and exposes an "Achilles heel" that may be therapeutically targeted for treating t(11;19) patients.

EBF1 primes B-lymphoid enhancers and limits the myeloid bias in murine multipotent progenitors

Hematopoietic stem cells (HSCs) and multipotent progenitors (MPPs) generate all cells of the blood system. Despite their multipotency, MPPs display poorly understood lineage bias. Here, we examine whether lineage-specifying transcription factors, such as the B-lineage determinant EBF1, regulate lineage preference in early progenitors. We detect low-level EBF1 expression in myeloid-biased MPP3 and lymphoid-biased MPP4 cells, coinciding with expression of the myeloid determinant C/EBPα. Hematopoietic deletion of Ebf1 results in enhanced myelopoiesis and reduced HSC repopulation capacity. Ebf1-deficient MPP3 and MPP4 cells exhibit an augmented myeloid differentiation potential and a transcriptome with an enriched C/EBPα signature. Correspondingly, EBF1 binds the Cebpa enhancer, and the deficiency and overexpression of Ebf1 in MPP3 and MPP4 cells lead to an up- and downregulation of Cebpa expression, respectively. In addition, EBF1 primes the chromatin of B-lymphoid enhancers specifically in MPP3 cells. Thus, our study implicates EBF1 in regulating myeloid/lymphoid fate bias in MPPs by constraining C/EBPα-driven myelopoiesis and priming the B-lymphoid fate.

Zebrafish neuromesodermal progenitors undergo a critical state transition in vivo

The transition state model of cell differentiation proposes that a transient window of gene expression stochasticity precedes entry into a differentiated state. Here, we assess this theoretical model in zebrafish neuromesodermal progenitors (NMps) in vivo during late somitogenesis stages. We observed an increase in gene expression variability at the 24 somite stage (24ss) before their differentiation into spinal cord and paraxial mesoderm. Analysis of a published 18ss scRNA-seq dataset showed that the NMp population is noisier than its derivatives. By building in silico composite gene expression maps from image data, we assigned an 'NM index' to in silico NMps based on the expression of neural and mesodermal markers and demonstrated that cell population heterogeneity peaked at 24ss. Further examination revealed cells with gene expression profiles incongruent with their prospective fate. Taken together, our work supports the transition state model within an endogenous cell fate decision making event.

Single-cell RNA sequencing uncovers the nuclear decoy lincRNA PIRAT as a regulator of systemic monocyte immunity during COVID-19

SARS-CoV-2–infected patients often display characteristic changes in the production of immune mediators that trigger life-threatening courses of COVID-19. The underlying molecular mechanisms are not yet fully understood. Here, we used single-cell RNA sequencing to investigate the involvement of the emerging class of long regulatory RNA in COVID-19. Our data reveal that a previously unknown regulatory RNA in the nucleus of immune cells is altered after SARS-CoV-2 infection. The degradation of this RNA removes a natural brake on the production of critical immune mediators that can promote the development of severe COVID-19. We believe that therapeutic intervention in this nuclear RNA circuit could counteract the overproduction of disease-causing immune mediators and protect against severe COVID-19.

The systemic immune response to viral infection is shaped by master transcription factors, such as NF-κB, STAT1, or PU.1. Although long noncoding RNAs (lncRNAs) have been suggested as important regulators of transcription factor activity, their contributions to the systemic immunopathologies observed during SARS-CoV-2 infection have remained unknown. Here, we employed a targeted single-cell RNA sequencing approach to reveal lncRNAs differentially expressed in blood leukocytes during severe COVID-19. Our results uncover the lncRNA PIRAT (PU.1-induced regulator of alarmin transcription) as a major PU.1 feedback-regulator in monocytes, governing the production of the alarmins S100A8/A9, key drivers of COVID-19 pathogenesis. Knockout and transgene expression, combined with chromatin-occupancy profiling, characterized PIRAT as a nuclear decoy RNA, keeping PU.1 from binding to alarmin promoters and promoting its binding to pseudogenes in naïve monocytes. NF-κB–dependent PIRAT down-regulation during COVID-19 consequently releases a transcriptional brake, fueling alarmin production. Alarmin expression is additionally enhanced by the up-regulation of the lncRNA LUCAT1, which promotes NF-κB–dependent gene expression at the expense of targets of the JAK-STAT pathway. Our results suggest a major role of nuclear noncoding RNA networks in systemic antiviral responses to SARS-CoV-2 in humans.

ETS1 governs pathological tissue-remodeling programs in disease-associated fibroblasts

Fibroblasts, the most abundant structural cells, exert homeostatic functions but also drive disease pathogenesis. Single-cell technologies have illuminated the shared characteristics of pathogenic fibroblasts in multiple diseases including autoimmune arthritis, cancer and inflammatory colitis. However, the molecular mechanisms underlying the disease-associated fibroblast phenotypes remain largely unclear. Here, we identify ETS1 as the key transcription factor governing the pathological tissue-remodeling programs in fibroblasts. In arthritis, ETS1 drives polarization toward tissue-destructive fibroblasts by orchestrating hitherto undescribed regulatory elements of the osteoclast differentiation factor receptor activator of nuclear factor-κB ligand (RANKL) as well as matrix metalloproteinases. Fibroblast-specific ETS1 deletion resulted in ameliorated bone and cartilage damage under arthritic conditions without affecting the inflammation level. Cross-tissue fibroblast single-cell data analyses and genetic loss-of-function experiments lent support to the notion that ETS1 defines the perturbation-specific fibroblasts shared among various disease settings. These findings provide a mechanistic basis for pathogenic fibroblast polarization and have important therapeutic implications.

CD11c identifies microbiota and EGR2-dependent MHCII+ serous cavity macrophages with sexually dimorphic fate in mice

The murine serous cavities contain a rare and enigmatic population of short-lived F4/80lo MHCII+ macrophages but what regulates their development, survival, and fate is unclear. Here, we show that mature F4/80lo MHCII+ peritoneal macrophages arise after birth, but that this occurs largely independently of colonization by microbiota. Rather, microbiota specifically regulate development of a subpopulation of CD11c+ cells that express the immunoregulatory cytokine RELM-α, are reliant on the transcription factor EGR2, and develop independently of the growth factor CSF1. Furthermore, we demonstrate that intrinsic expression of RELM-α, a signature marker shared by CD11c+ and CD11c- F4/80lo MHCII+ cavity macrophages, regulates survival and differentiation of these cells in the peritoneal cavity in a sex-specific manner. Thus, we identify a previously unappreciated diversity in serous cavity F4/80lo MHCII+ macrophages that is regulated by microbiota, and describe a novel sex and site-specific function for RELM-α in regulating macrophage endurance that reveals the unique survival challenge presented to monocyte-derived macrophages by the female peritoneal environment.

Early Growth Response Factor 1 in Aging Hematopoietic Stem Cells and Leukemia

Aging is associated with various hematological disorders and a higher risk of myeloproliferative disorders. An aged hematopoietic system can be characterized by decreased immune function and increased myeloid cell production. Hematopoietic stem cells (HSCs) regulate the production of blood cells throughout life. The self-renewal and regenerative potential of HSCs determine the quality and quantity of the peripheral blood cells. External signals from the microenvironment under different conditions determine the fate of the HSCs to proliferate, self-renew, differentiate, or remain quiescent. HSCs respond impromptu to a vast array of extracellular signaling cascades such as cytokines, growth factors, or nutrients, which are crucial in the regulation of HSCs. Early growth response factor 1 (EGR1) is one of the key transcription factors controlling HSC proliferation and their localization in the bone marrow (BM) niche. Downregulation of Egr1 activates and recruits HSCs for their proliferation and differentiation to produce mature blood cells. Increased expression of Egr1 is implicated in immuno-aging of HSCs. However, dysregulation of Egr1 is associated with hematological malignancies such as acute myeloid leukemia (AML), acute lymphoblastic leukemia (ALL), and chronic myelogenous leukemia (CML). Here, we summarize the current understanding of the role of EGR1 in the regulation of HSC functionality and the manifestation of leukemia. We also discuss the alternative strategies to rejuvenate the aged HSCs by targeting EGR1 in different settings.

From developmental to atavistic bet-hedging: How cancer cells pervert the exploitation of random single-cell phenotypic fluctuations

Stochastic gene expression plays a leading developmental role through its contribution to cell differentiation. It is also proposed to promote phenotypic diversification in malignant cells. However, it remains unclear if these two forms of cellular bet-hedging are identical or rather display distinct features. Here we argue that bet-hedging phenomena in cancer cells are more similar to those occurring in unicellular organisms than to those of normal metazoan cells. We further propose that the atavistic bet-hedging strategies in cancer originate from a hijacking of the normal developmental bet-hedging of metazoans. Finally, we discuss the constraints that may shape the atavistic bet-hedging strategies of cancer cells.

Folate Carrier Deficiency Drives Differential Methylation and Enhanced Cellular Potency in the Neural Plate Border

The neural plate border (NPB) of vertebrate embryos segregates from the neural and epidermal regions, and it is comprised of an intermingled group of multipotent progenitor cells. Folate is the precursor of S-adenosylmethionine, the main methyl donor for DNA methylation, and it is critical for embryonic development, including the specification of progenitors which reside in the NPB. Despite the fact that several intersecting signals involved in the specification and territorial restriction of NPB cells are known, the role of epigenetics, particularly DNA methylation, has been a matter of debate. Here, we examined the temporal and spatial distribution of the methyl source and analyzed the abundance of 5mC/5 hmC and their epigenetic writers throughout the segregation of the neural and NPB territories. Reduced representation bisulfite sequencing (RRBS) on Reduced Folate Carrier 1 (RFC1)-deficient embryos leads to the identification of differentially methylated regions (DMRs). In the RFC1-deficient embryos, we identified several DMRs in the Notch1 locus, and the spatiotemporal expression of Notch1 and its downstream target gene Bmp4 were expanded in the NPB. Cell fate analysis on folate deficient embryos revealed a significant increase in the number of cells coexpressing both neural (SOX2) and NPB (PAX7) markers, which may represent an enhancing effect in the cellular potential of those progenitors. Taken together, our findings propose a model where the RFC1 deficiency drives methylation changes in specific genomic regions that are correlated with a dysregulation of pathways involved in early development such as Notch1 and BMP4 signaling. These changes affect the potency of the progenitors residing in the juncture of the neural plate and NPB territories, thus driving them to a primed state.

Quantitative Analysis of the Transcriptome of Two Commonly Used Human Monocytic Cell Lines—THP-1 and Mono Mac 6—Reveals Their Arrest during Early Monocyte/Neutrophil Differentiation

Cell lines of monocyte/macrophage origin are often used as model systems to study monocyte/macrophage biology. A relevant question is how similar these cell lines are to their in vivo counterparts? To address this issue, we performed a detailed analysis of the transcriptome of two commonly used human monocyte/macrophage cell lines, Mono Mac 6 and THP-1. Both of these cell lines originate from leukemic cells with myelo-monocytic characteristics. We found that both Mono Mac 6 and THP-1 represent cells of very immature origin. Their transcriptomes show more similarities to immature neutrophils than cells of the monocyte/macrophage lineage. They express significant levels of N-elastase, proteinase 3, cathepsin G, and azurocidin but very low levels of CD14, ficolin, and complement factor P. All major MHC class II genes are also expressed at low levels. They show high levels of lysozyme and low levels of one of the immunoglobulin Fc receptors, FCGRIIA, which is characteristic of both neutrophils and monocytes. THP-1, but not Mono Mac 6, also expresses the high-affinity receptor for IgG, FCGRIA. Both cell lines lack the expression of the connective tissue components fibronectin, proteoglycan 4, and syndecan 3, which are characteristics of tissue macrophages but are absent in blood monocytes, indicating that they originate from bone marrow precursors and not yolk sac-derived hematopoietic cells. Both of these cell lines seem, therefore, to represent cells arrested during early myelo-monocytic development, at a branch point between neutrophil and monocyte differentiation. Their very immature phenotype indicates that great care should be taken when using these cell lines as models for normal monocyte/macrophage biology.

The journey of neutropoiesis: how complex landscapes in bone marrow guide continuous neutrophil lineage determination

Neutrophils are the most abundant white blood cell, and they differentiate in homeostasis in the bone marrow from hematopoietic stem cells (HSCs) via multiple intermediate progenitor cells into mature cells that enter the circulation. Recent findings support a continuous model of differentiation in the bone marrow of heterogeneous HSCs and progenitor populations. Cell fate decisions at the levels of proliferation and differentiation are enforced through expression of lineage-determining transcription factors and their interactions, which are influenced by intrinsic (intracellular) and extrinsic (extracellular) mechanisms. Neutrophil homeostasis is subjected to positive-feedback loops, stemming from the gut microbiome, as well as negative-feedback loops resulting from the clearance of apoptotic neutrophils by mature macrophages. Finally, the cellular kinetics regarding the replenishing of the mature neutrophil pool is discussed in light of recent contradictory data.

Identification of dynamic driver sets controlling phenotypical landscapes

Controlling phenotypical landscapes is of vital interest to modern biology. This task becomes highly demanding because cellular decisions involve complex networks engaging in crosstalk interactions. Previous work on control theory indicates that small sets of compounds can control single phenotypes. However, a dynamic approach is missing to determine the drivers of the whole network dynamics. By analyzing 35 biologically motivated Boolean networks, we developed a method to identify small sets of compounds sufficient to decide on the entire phenotypical landscape. These compounds do not strictly prefer highly related compounds and show a smaller impact on the stability of the attractor landscape. The dynamic driver sets include many intervention targets and cellular reprogramming drivers in human networks. Finally, by using a new comprehensive model of colorectal cancer, we provide a complete workflow on how to implement our approach to shift from in silico to in vitro guided experiments.

Inflammation in Focus: The Beginning and the End

The inflammation is an important biological response induced by various harmful stimuli, like viruses, bacterial infections, toxins, toxic compounds, tissue injury. During inflammation inflammatory cytokines and reactive oxygen species are produced. Inflammatory cytokines act on various receptors present on the plasma membrane of target cells. To initiate signaling cascade, and activate transcription factors, receptors should be internalized and enter the early endosomes, where the members of the signaling cascade can meet. The further cytoplasmic fate of the receptor plays crucial role in the progression and the course of inflammation. Usually acute inflammation removes injurious stimuli and helps to regain the normal healthy status of the organism. In contrast to this the uncontrolled chronic inflammation—stimulating other than immune cells, inducing transdifferentiation—can provide base of various serious diseases. This paper draws the attention of the long-lasting consequence of chronic inflammation, pointing out that one of the most important step in medication is to identify in time the factors initiating and maintaining inflammation.

Transcriptional Control of Cell Fate Determination in Antigen-Experienced CD8 T Cells

Robust immunity to intracellular infections is mediated by antigen-specific naive CD8 T cells that become activated and differentiate into phenotypically and functionally diverse subsets of effector cells, some of which terminally differentiate and others that give rise to memory cells that provide long-lived protection. This developmental system is an outstanding model with which to elucidate how regulation of chromatin structure and transcriptional control establish gene expression programs that govern cell fate determination, insights from which are likely to be useful for informing the design of immunotherapeutic approaches to engineer durable immunity to infections and tumors. A unifying framework that describes how naive CD8 T cells develop into memory cells is still outstanding. We propose a model that incorporates a common early linear path followed by divergent paths that slowly lose capacity to interconvert and discuss classical and contemporary observations that support these notions, focusing on insights from transcriptional control and chromatin regulation.

Synthetic multistability in mammalian cells

In multicellular organisms, gene regulatory circuits generate thousands of molecularly distinct, mitotically heritable states through the property of multistability. Designing synthetic multistable circuits would provide insight into natural cell fate control circuit architectures and would allow engineering of multicellular programs that require interactions among distinct cell types. We created MultiFate, a naturally inspired, synthetic circuit that supports long-term, controllable, and expandable multistability in mammalian cells. MultiFate uses engineered zinc finger transcription factors that transcriptionally self-activate as homodimers and mutually inhibit one another through heterodimerization. Using a model-based design, we engineered MultiFate circuits that generate as many as seven states, each stable for at least 18 days. MultiFate permits controlled state switching and modulation of state stability through external inputs and can be expanded with additional transcription factors. These results provide a foundation for engineering multicellular behaviors in mammalian cells.

Myeloid Immune Cells CARrying a New Weapon Against Cancer

Chimeric antigen receptor (CAR) engineering for T cells and natural killer cells (NK) are now under clinical evaluation for the treatment of hematologic cancers. Although encouraging clinical results have been reported for hematologic diseases, pre-clinical studies in solid tumors have failed to prove the same effectiveness. Thus, there is a growing interest of the scientific community to find other immune cell candidate to express CAR for the treatment of solid tumors and other diseases. Mononuclear phagocytes may be the most adapted group of cells with potential to overcome the dense barrier imposed by solid tumors. In addition, intrinsic features of these cells, such as migration, phagocytic capability, release of soluble factors and adaptive immunity activation, could be further explored along with gene therapy approaches. Here, we discuss the elements that constitute the tumor microenvironment, the features and advantages of these cell subtypes and the latest studies using CAR-myeloid immune cells in solid tumor models.

Blood stem cell PU.1 upregulation is a consequence of differentiation without fast autoregulation

Transcription factors (TFs) regulate cell fates, and their expression must be tightly regulated. Autoregulation is assumed to regulate many TFs' own expression to control cell fates. Here, we manipulate and quantify the (auto)regulation of PU.1, a TF controlling hematopoietic stem and progenitor cells (HSPCs), and correlate it to their future fates. We generate transgenic mice allowing both inducible activation of PU.1 and noninvasive quantification of endogenous PU.1 protein expression. The quantified HSPC PU.1 dynamics show that PU.1 up-regulation occurs as a consequence of hematopoietic differentiation independently of direct fast autoregulation. In contrast, inflammatory signaling induces fast PU.1 up-regulation, which does not require PU.1 expression or its binding to its own autoregulatory enhancer. However, the increased PU.1 levels induced by inflammatory signaling cannot be sustained via autoregulation after removal of the signaling stimulus. We conclude that PU.1 overexpression induces HSC differentiation before PU.1 up-regulation, only later generating cell types with intrinsically higher PU.1.

Data-driven modeling predicts gene regulatory network dynamics during the differentiation of multipotential hematopoietic progenitors

Cellular differentiation during hematopoiesis is guided by gene regulatory networks (GRNs) comprising transcription factors (TFs) and the effectors of cytokine signaling. Based largely on analyses conducted at steady state, these GRNs are thought to be organized as a hierarchy of bistable switches, with antagonism between Gata1 and PU.1 driving red- and white-blood cell differentiation. Here, we utilize transient gene expression patterns to infer the genetic architecture—the type and strength of regulatory interconnections—and dynamics of a twelve-gene GRN including key TFs and cytokine receptors. We trained gene circuits, dynamical models that learn genetic architecture, on high temporal-resolution gene-expression data from the differentiation of an inducible cell line into erythrocytes and neutrophils. The model is able to predict the consequences of gene knockout, knockdown, and overexpression experiments and the inferred interconnections are largely consistent with prior empirical evidence. The inferred genetic architecture is densely interconnected rather than hierarchical, featuring extensive cross-antagonism between genes from alternative lineages and positive feedback from cytokine receptors. The analysis of the dynamics of gene regulation in the model reveals that PU.1 is one of the last genes to be upregulated in neutrophil conditions and that the upregulation of PU.1 and other neutrophil genes is driven by Cebpa and Gfi1 instead. This model inference is confirmed in an independent single-cell RNA-Seq dataset from mouse bone marrow in which Cebpa and Gfi1 expression precedes the neutrophil-specific upregulation of PU.1 during differentiation. These results demonstrate that full PU.1 upregulation during neutrophil development involves regulatory influences extrinsic to the Gata1-PU.1 bistable switch. Furthermore, although there is extensive cross-antagonism between erythroid and neutrophil genes, it does not have a hierarchical structure. More generally, we show that the combination of high-resolution time series data and data-driven dynamical modeling can uncover the dynamics and causality of developmental events that might otherwise be obscured.

The supply of blood cells is replenished by the maturation of hematopoietic progenitor cells into different cell types. Which cell type a progenitor cell develops into is determined by a complex network of genes whose protein products directly or indirectly regulate each others’ expression and that of downstream genes characteristic of the cell type. We inferred the nature and causality of the regulatory connections in a 12-gene network known to affect the decision between erythrocyte and neutrophil cell fates using a predictive machine-learning approach. Our analysis showed that the overall architecture of the network is densely interconnected and not hierarchical. Furthermore, the model inferred that PU.1, considered a master regulator of all white-blood cell lineages, is upregulated during neutrophil development by two other proteins, Cebpa and Gfi1. We validated this prediction by showing that Cebpa and Gfi1 expression precedes that of PU.1 in single-cell gene expression data from mouse bone marrow. These results revise the architecture of the gene network and the causality of regulatory events guiding hematopoiesis. The results also show that combining machine learning approaches with time course data can help resolve causality during development.

Mechanistic models of blood cell fate decisions in the era of single-cell data

Billions of functionally distinct blood cells emerge from a pool of hematopoietic stem cells in our bodies every day. This progressive differentiation process is hierarchically structured and remarkably robust. We provide an introductory review to mathematical approaches addressing the functional aspects of how lineage choice is potentially implemented on a molecular level. Emerging from studies on the mutual repression of key transcription factors, we illustrate how those simple concepts have been challenged in recent years and subsequently extended. Especially, the analysis of omics data on the single-cell level with computational tools provides descriptive insights on a yet unknown level, while their embedding into a consistent mechanistic and mathematical framework is still incomplete.

Monostability and bistability of biological switches

Cell-fate transition can be modeled by ordinary differential equations (ODEs) which describe the behavior of several molecules in interaction, and for which each stable equilibrium corresponds to a possible phenotype (or 'biological trait'). In this paper, we focus on simple ODE systems modeling two molecules which each negatively (or positively) regulate the other. It is well-known that such models may lead to monostability or multistability, depending on the selected parameters. However, extensive numerical simulations have led systems biologists to conjecture that in the vast majority of cases, there cannot be more than two stable points. Our main result is a proof of this conjecture. More specifically, we provide a criterion ensuring at most bistability, which is indeed satisfied by most commonly used functions. This includes Hill functions, but also a wide family of convex and sigmoid functions. We also determine which parameters lead to monostability, and which lead to bistability, by developing a more general framework encompassing all our results.

Emergent Bistable Switches from the Incoherent Feed-Forward Signaling of a Positive Feedback Loop

Bistability is intrinsically connected to various decision making processes in living systems. The operating principles of a bistable switch, generated from a positive feedback loop, are well understood both in natural and synthetic settings. However, the fate of dynamic modularity of a positive feedback loop is unknown when it is connected to another dynamically modular signaling motif. In order to address this, here we investigate feed-forward signaling of a positive feedback loop to determine the fate of a bistable switch under such signaling. Using the potential energy based high-throughput bifurcation analysis method, we uncover that in addition to the conventional bistability the hybrid motifs generate various emergent bistable switches, namely mushroom and isola switches, which are not produced by the individual motifs. Using random parameter sampling, network perturbation, and phase plane analysis, we establish the design principles of such emergent behaviors. Incoherent feed-forward signaling of a positive feedback loop with distinct regulatory thresholds of the two arms of the feed-forward loop are the key requirements for such emergent behaviors. Our calculations show that the specific types of atypical bistable responses depend on the logic gate configuration of the signals. However, the emergent bistable behaviors of the hybrid networks do not depend on the nature of the positive feedback loop.

The transcription factor EGR2 is indispensable for tissue-specific imprinting of alveolar macrophages in health and tissue repair

Alveolar macrophages are the most abundant macrophages in the healthy lung where they play key roles in homeostasis and immune surveillance against airborne pathogens. Tissue-specific differentiation and survival of alveolar macrophages rely on niche-derived factors, such as granulocyte-macrophage colony-stimulating factor (GM-CSF) and transforming growth factor–β (TGF-β). However, the nature of the downstream molecular pathways that regulate the identity and function of alveolar macrophages and their response to injury remain poorly understood. Here, we identify that the transcription factor EGR2 is an evolutionarily conserved feature of lung alveolar macrophages and show that cell-intrinsic EGR2 is indispensable for the tissue-specific identity of alveolar macrophages. Mechanistically, we show that EGR2 is driven by TGF-β and GM-CSF in a PPAR-γ–dependent manner to control alveolar macrophage differentiation. Functionally, EGR2 was dispensable for the regulation of lipids in the airways but crucial for the effective handling of the respiratory pathogen Streptococcus pneumoniae. Last, we show that EGR2 is required for repopulation of the alveolar niche after sterile, bleomycin-induced lung injury and demonstrate that EGR2-dependent, monocyte-derived alveolar macrophages are vital for effective tissue repair after injury. Collectively, we demonstrate that EGR2 is an indispensable component of the transcriptional network controlling the identity and function of alveolar macrophages in health and disease.

Cyclical fate restriction: a new view of neural crest cell fate specification

Neural crest cells are crucial in development, not least because of their remarkable multipotency. Early findings stimulated two hypotheses for how fate specification and commitment from fully multipotent neural crest cells might occur, progressive fate restriction (PFR) and direct fate restriction, differing in whether partially restricted intermediates were involved. Initially hotly debated, they remain unreconciled, although PFR has become favoured. However, testing of a PFR hypothesis of zebrafish pigment cell development refutes this view. We propose a novel 'cyclical fate restriction' hypothesis, based upon a more dynamic view of transcriptional states, reconciling the experimental evidence underpinning the traditional hypotheses.

An erythroid-to-myeloid cell fate conversion is elicited by LSD1 inactivation

Histone H3 lysine 4 methylation (H3K4Me) is most often associated with chromatin activation, and removing H3K4 methyl groups has been shown to be coincident with gene repression. H3K4Me demethylase KDM1a/LSD1 is a therapeutic target for multiple diseases, including for the potential treatment of β-globinopathies (sickle cell disease and β-thalassemia), because it is a component of γ-globin repressor complexes, and LSD1 inactivation leads to robust induction of the fetal globin genes. The effects of LSD1 inhibition in definitive erythropoiesis are not well characterized, so we examined the consequences of conditional inactivation of Lsd1 in adult red blood cells using a new Gata1creERT2 bacterial artificial chromosome transgene. Erythroid-specific loss of Lsd1 activity in mice led to a block in erythroid progenitor differentiation and to the expansion of granulocyte-monocyte progenitor-like cells, converting hematopoietic differentiation potential from an erythroid fate to a myeloid fate. The analogous phenotype was also observed in human hematopoietic stem and progenitor cells, coincident with the induction of myeloid transcription factors (eg, PU.1 and CEBPα). Finally, blocking the activity of the transcription factor PU.1 or RUNX1 at the same time as LSD1 inhibition rescued myeloid lineage conversion to an erythroid phenotype. These data show that LSD1 promotes erythropoiesis by repressing myeloid cell fate in adult erythroid progenitors and that inhibition of the myeloid-differentiation pathway reverses the lineage switch induced by LSD1 inactivation.