Lineage Decisions Research Articles

The fate of secretory cells during intestinal homeostasis, regeneration, and tumor formation is regulated by Tcf4.

The single-layer epithelium of the gastrointestinal tract is a dynamically renewing tissue that ensures nutrient absorption, secretory and barrier functions and is involved in immune responses. The basis for this homeostatic renewal is the Wnt signaling pathway. Blocking this pathway can lead to epithelial damage, while its abnormal activation can result in the development of intestinal tumors. In this study, we investigated the dynamics of intestinal epithelial cells and tumorigenesis using a conditional mouse model. Using single-cell and bulk RNA sequencing and histological analysis, we elucidated the cellular responses following the loss of specific cell types. We focused on the fate of cells in the lower parts of the intestinal crypts and the development of colon adenomas. By partially inactivating the transcription factor Tcf4, a key effector of the Wnt signaling pathway, we analyzed the regeneration of isolated hyperproliferative foci (crypts). Our results suggest that the damaged epithelium is not restored by a specific regeneration program associated with oncofetal gene production, but rather by a standard homeostatic renewal pathway. Moreover, disruption of Tcf4 in secretory progenitors resulted in a significant shift in the cell lineage from Paneth cells to goblet cells, characterized by morphological changes and loss of Paneth cell-specific genes. We also found that hyperactivation of the Wnt signaling pathway in colonic adenomas correlated with the upregulation of genes typical of Paneth cells in the intestine, followed by the emergence of secretory tumor cells producing the Wnt3 ligand. The absence of Tcf4 led to a phenotypic shift of the tumor cells towards goblet cells. Our study presents a new model of epithelial regeneration based on the genetically driven partial elimination of intestinal crypts. We highlight the critical role of Tcf4 in the control of cell lineage decisions in the intestinal epithelium and colon tumors.

Heterogeneous stiffness of the bone marrow microenvironment regulates the fate decision of haematopoietic stem and progenitor cells.

The bone marrow (BM) niches are the complex microenvironments that surround cells, providing various external stimuli to regulate a range of haematopoietic stem cell (HSC) behaviours. Recently, it has been proposed that the fate decision of HSCs is often correlated with significantly altered biophysical signals of BM niches. To thoroughly elucidate the effect of mechanical microenvironments on cell fates, we constructed 2D and 3D cell culture hydrogels using polyacrylamide to replicate the mechanical properties of heterogeneous sub-niches, including the inherent rigidity of marrow adipose tissue (2 kPa), perivascular tissue (8 kPa) and endosteum region (35 kPa) in BM. Our observations suggest that HSCs can respond to the mechanical heterogeneity of the BM microenvironment, exhibiting diversity in cell mechanics, haematopoietic pool maintenance and differentiated lineages. Hydrogels with higher stiffness promote the preservation of long-term repopulating HSCs (LT-HSCs), while those with lower stiffness support multi-potent progenitors (MPPs) viability in vitro. Furthermore, we established a comprehensive transcriptional profile of haematopoietic subpopulations to reflect the multipotency of haematopoietic stem and progenitor cells (HSPCs) that are modulated by niche-like stiffness. Our findings demonstrate that HSPCs exhibit completely distinct downstream differentiated preferences within hydrogel systems of varying stiffness. This highlights the crucial role of tissue-specific mechanical properties in HSC lineage decisions, which may provide innovative solutions to clinical challenges.

Effects of miR-29c on proliferation and adipogenic differentiation of porcine bone marrow mesenchymal stromal cells

ABSTRACT microRNAs (miRNAs), a subclass of noncoding short RNAs, direct cells fate decisions that are important for cell proliferation and cell lineage decisions. Adipogenic differentiation contributes greatly to the development of white adipose tissue, involving of highly organized regulation by miRNAs. In the present study, we screened and identified 78 differently expressed miRNAs of porcine BMSCs during adipogenic differentiation. Of which, the role of miR-29c in regulating the proliferation and adipogenic differentiation was proved and detailed. Specifically, over-expression miR-29c inhibits the proliferation and adipogenic differentiation of BMSCs, which were reversed upon miR-29c inhibitor. Interference of IGF1 inhibits the proliferation and adipogenic differentiation of BMSCs. Mechanistically, miR-29c regulates the proliferation and adipogenic differentiation of BMSCs by targeting IGF1 and further regulating the MAPK pathway and the PI3K-AKT-mTOR pathway, respectively. In conclusion, we highlight the important role of miR-29c in regulating proliferation and adipogenic differentiation of BMSCs.

Smad4 is essential for epiblast scaling and morphogenesis after implantation, but nonessential before implantation.

Bone morphogenic protein (BMP) signaling plays an essential and highly conserved role in embryo axial patterning in animal species. However, in mammalian embryos, which develop inside the mother, early development includes a preimplantation stage, which does not occur in externally developing embryos. During preimplantation, the epiblast is segregated from extra-embryonic lineages that enable implantation and development in utero. Yet, the requirement for BMP signaling is imprecisely defined in mouse early embryos. Here, we show that, in contrast to previous reports, BMP signaling (SMAD1/5/9 phosphorylation) is not detectable until implantation when it is detected in the primitive endoderm - an extra-embryonic lineage. Moreover, preimplantation development appears to be normal following deletion of maternal and zygotic Smad4, an essential effector of canonical BMP signaling. In fact, mice lacking maternal Smad4 are viable. Finally, we uncover a new requirement for zygotic Smad4 in epiblast scaling and cavitation immediately after implantation, via a mechanism involving FGFR/ERK attenuation. Altogether, our results demonstrate no role for BMP4/SMAD4 in the first lineage decisions during mouse development. Rather, multi-pathway signaling among embryonic and extra-embryonic cell types drives epiblast morphogenesis postimplantation.

Branching topology of the human embryo transcriptome revealed by Entropy Sort Feature Weighting.

Analysis of single cell transcriptomics (scRNA-seq) data is typically performed after subsetting to highly variable genes (HVGs). Here, we show that Entropy Sorting provides an alternative mathematical framework for feature selection. On synthetic datasets, continuous Entropy Sort Feature Weighting (cESFW) outperforms HVG selection in distinguishing cell-state-specific genes. We apply cESFW to six merged scRNA-seq datasets spanning human early embryo development. Without smoothing or augmenting the raw counts matrices, cESFW generates a high-resolution embedding displaying coherent developmental progression from eight-cell to post-implantation stages and delineating 15 distinct cell states. The embedding highlights sequential lineage decisions during blastocyst development, while unsupervised clustering identifies branch point populations obscured in previous analyses. The first branching region, where morula cells become specified for inner cell mass or trophectoderm, includes cells previously asserted to lack a developmental trajectory. We quantify the relatedness of different pluripotent stem cell cultures to distinct embryo cell types and identify marker genes of naïve and primed pluripotency. Finally, by revealing genes with dynamic lineage-specific expression, we provide markers for staging progression from morula to blastocyst.

Feedback between stochastic gene networks and population dynamics enables cellular decision-making.

Phenotypic selection occurs when genetically identical cells are subject to different reproductive abilities due to cellular noise. Such noise arises from fluctuations in reactions synthesizing proteins and plays a crucial role in how cells make decisions and respond to stress or drugs. We propose a general stochastic agent-based model for growing populations capturing the feedback between gene expression and cell division dynamics. We devise a finite state projection approach to analyze gene expression and division distributions and infer selection from single-cell data in mother machines and lineage trees. We use the theory to quantify selection in multi-stable gene expression networks and elucidate that the trade-off between phenotypic switching and selection enables robust decision-making essential for synthetic circuits and developmental lineage decisions. Using live-cell data, we demonstrate that combining theory and inference provides quantitative insights into bet-hedging-like response to DNA damage and adaptation during antibiotic exposure in Escherichia coli.

Investigating the basis of lineage decisions and developmental trajectories in the dorsal spinal cord through pseudotime analyses.

Dorsal interneurons (dIs) in the spinal cord encode the perception of touch, pain, heat, itchiness and proprioception. Previous studies using genetic strategies in animal models have revealed important insights into dI development, but the molecular details of how dIs arise as distinct populations of neurons remain incomplete. We have developed a resource to investigate dI fate specification by combining a single-cell RNA-Seq atlas of mouse embryonic stem cell-derived dIs with pseudotime analyses. To validate this in silico resource as a useful tool, we used it to first identify genes that are candidates for directing the transition states that lead to distinct dI lineage trajectories, and then validated them using in situ hybridization analyses in the developing mouse spinal cord in vivo. We have also identified an endpoint of the dI5 lineage trajectory and found that dIs become more transcriptionally homogeneous during terminal differentiation. This study introduces a valuable tool for further discovery about the timing of gene expression during dI differentiation and demonstrates its utility in clarifying dI lineage relationships.

A geometrical model of cell fate specification in the mouse blastocyst.

The lineage decision that generates the epiblast and primitive endoderm from the inner cell mass (ICM) is a paradigm for cell fate specification. Recent mathematics has formalized Waddington's landscape metaphor and proven that lineage decisions in detailed gene network models must conform to a small list of low-dimensional stereotypic changes called bifurcations. The most plausible bifurcation for the ICM is the so-called heteroclinic flip that we define and elaborate here. Our re-analysis of recent data suggests that there is sufficient cell movement in the ICM so the FGF signal, which drives the lineage decision, can be treated as spatially uniform. We thus extend the bifurcation model for a single cell to the entire ICM by means of a self-consistently defined time-dependent FGF signal. This model is consistent with available data and we propose additional dynamic experiments to test it further. This demonstrates that simplified, quantitative and intuitively transparent descriptions are possible when attention is shifted from specific genes to lineages. The flip bifurcation is a very plausible model for any situation where the embryo needs control over the relative proportions of two fates by a morphogen feedback.

Thyroid hormone action in adult neurogliogenic niches: the known and unknown.

Over the last decades, thyroid hormones (THs) signaling has been established as a key signaling cue for the proper maintenance of brain functions in adult mammals, including humans. One of the most fascinating roles of THs in the mature mammalian brain is their ability to regulate adult neurogliogenic processes. In this respect, THs control the generation of new neuronal and glial progenitors from neural stem cells (NSCs) as well as their final differentiation and maturation programs. In this review, we summarize current knowledge on the cellular organization of adult rodent neurogliogenic niches encompassing well-established niches in the subventricular zone (SVZ) lining the lateral ventricles, the hippocampal subgranular zone (SGZ), and the hypothalamus, but also less characterized niches in the striatum and the cerebral cortex. We then discuss critical questions regarding how THs availability is regulated in the respective niches in rodents and larger mammals as well as how modulating THs availability in those niches interferes with lineage decision and progression at the molecular, cellular, and functional levels. Based on those alterations, we explore the novel therapeutic avenues aiming at harnessing THs regulatory influences on neurogliogenic output to stimulate repair processes by influencing the generation of either new neurons (i.e. Alzheimer's, Parkinson's diseases), oligodendrocytes (multiple sclerosis) or both (stroke). Finally, we point out future challenges, which will shape research in this exciting field in the upcoming years.

TBX3 reciprocally controls key trophoblast lineage decisions in villi during human placenta development in the first trimester

Human trophoblastic lineage development is intertwined with placental development and pregnancy outcomes, but the regulatory mechanisms underpinning this process remain inadequately understood. In this study, based on single-nuclei RNA sequencing (snRNA-seq) analysis of the human early maternal-fetal interface, we compared the gene expression pattern of trophoblast at different developmental stages. Our findings reveal a predominant upregulation of TBX3 during the transition from villous cytotrophoblast (VCT) to syncytiotrophoblast (SCT), but downregulation of TBX3 as VCT progresses into extravillous trophoblast cells (EVT). Immunofluorescence analysis verified the primary expression of TBX3 in SCT, partial expression in MKi67-positive VCT, and absence in HLA-G-positive EVT, consistent with our snRNA-seq results. Using immortalized trophoblastic cell lines (BeWo and HTR8/SVneo) and human primary trophoblast stem cells (hTSCs), we observed that TBX3 knockdown impedes SCT formation through RAS-MAPK signaling, while TBX3 overexpression disrupts the cytoskeleton structure of EVT and hinders EVT differentiation by suppressing FAK signaling.In conclusion, our study suggests that the spatiotemporal expression of TBX3 plays a critical role in regulating trophoblastic lineage development via distinct signaling pathways. This underscores TBX3 as a key determinant during hemochorial placental development.

DeepVelo: deep learning extends RNA velocity to multi-lineage systems with cell-specific kinetics

Existing RNA velocity estimation methods strongly rely on predefined dynamics and cell-agnostic constant transcriptional kinetic rates, assumptions often violated in complex and heterogeneous single-cell RNA sequencing (scRNA-seq) data. Using a graph convolution network, DeepVelo overcomes these limitations by generalizing RNA velocity to cell populations containing time-dependent kinetics and multiple lineages. DeepVelo infers time-varying cellular rates of transcription, splicing, and degradation, recovers each cell’s stage in the differentiation process, and detects functionally relevant driver genes regulating these processes. Application to various developmental and pathogenic processes demonstrates DeepVelo’s capacity to study complex differentiation and lineage decision events in heterogeneous scRNA-seq data.

BCOR Loss Confers Increased Stemness and Partially Rescues RUNX1-Deficient Phenotypes in Human Hematopoietic Stem and Progenitor Cells

Introduction: Loss of function mutations in BCL-6 corepressor (BCOR) commonly occur in de novo acute myeloid leukemia (AML) and germline hematopoietic malignancies, most notably in those arising from RUNX1 familial platelet disorder (RUNX1-FPD). Second only to loss of heterozygosity of RUNX1, mutations in BCOR are the most commonly acquired mutations upon progression to myeloid malignancy in RUNX1-FPD. Additionally, RUNX1 and BCOR mutations commonly co-occur in de novo AML, and both are classified as “adverse risk” by the European LeukemiaNet if found at AML diagnosis. However, the role of BCOR in leukemogenesis is poorly understood, and there are limited studies investigating co-loss of BCOR and RUNX1 in hematopoietic stem and progenitor cells (HSPCs). Methods: We developed a model of BCOR loss in primary human CD34+ HSPCs using CRISPR/Cas9 and AAV6-mediated homology directed repair to knock-in a fluorescent reporter expression cassette into the BCOR locus. We used this system to investigate how BCOR knockout (KO) perturbs normal HSPC function. We combined this approach with a previously developed analogous system to disrupt the RUNX1 locus to investigatehow BCOR loss in the setting of RUNX1 KO modulates the effects of RUNX1 loss on HSPC function. Results: Here, we found that BCOR loss increased capacity for serial replating and conferred a competitive advantage over control HSPCs. In serial replating assays, BCOR-deficient HSPCs had a 5-fold increase in colony formation after one passage with a striking majority of colonies being erythroid. In liquid differentiation assays, BCOR-deficient HSPCs consistently had a 2-fold increase in CD34+ cells over a 21 day time course. When co-cultured with control HSPCs in a 1:1 ratio in stem expansion media (base media + SCF, FLT3-L, TPO, IL-6, UM-171), BCOR-deficient HSPCs expanded to comprise 70% of the culture at 14 days. This competitive advantage was recapitulated in vivo. At the 12-week timepoint of a competitive transplantation assay in NSG mice using HSPCs from 3 different donors, BCOR KO cells (marked by BFP) on average made up 75% of the engrafted human CD45+ population, while control cells (expressing GFP from the AAVS1 safe harbor locus) made up only 21%. Additionally, BCOR loss partially rescued RUNX1-induced stem cell dysfunction in vitro. Compared to RUNX1 deficiency alone, HSPCs deficient in both RUNX1 and BCOR had a 4-fold increase in colony formation after two passages in serial replating assays. Our lab has previously shown that RUNX1 loss ablates erythroid colony formation and skews towards monocytic colonies. Despite BCOR loss alone having a strong skew towards erythroid colonies, BCOR loss in RUNX1-deficient HSPCs did not induce erythroid colony formation. This suggests that in double mutant HSPCs, RUNX1 loss drives lineage decisions while BCOR loss drives the increase in serial replating potential. Ongoing efforts are aimed at characterizing combinatorial effects of BCOR and RUNX1 loss on hematopoiesis and stem cell function through in vitro liquid differentiation assays and in vivo transplantations, as well as investigating the molecular underpinnings of these phenotypes through RNA-seq. Conclusion: In summary, we have established a model of BCOR loss in primary human HSPCs to evaluate not only the role of BCOR in hematopoiesis, but also to characterize how BCOR is involved in leukemic transformation when combined with RUNX1. We show that BCOR deficiency enhances several attributes associated with stem cell function, including increased potential for serial replating, increased CD34 positive population, and increased engraftment compared to control HSPCs. When combined with RUNX1 loss, BCOR loss increases the serial replating potential of RUNX1-deficient cells, suggesting some rescue of RUNX1 KO-induced stem cell dysfunction. Ultimately, this model will enable us to elucidate the molecular mechanisms underlying increased stemness and competitive advantage in BCOR-deficient HSPCs, interrogate RUNX1-FPD disease progression associated with BCOR mutation, and validate potential therapeutics.

Myelo-Erythroid Lineage Segregation Is Regulated By the GATA-2 Interactome

Hematopoiesis arises from hematopoietic stem cell differentiating into downstream progenitors through a multistep process of sequential lineage potential restriction leading to the final production of >10 hematopoietic cell type.The regulation of hematopoietic lineage development involves synergistic and antagonistic interactions between transcription factors (TF) with lineage-restricted expression patterns, both through protein-protein interaction and direct transcriptional activation.Eosinophils, basophils and mast cells are key myeloid cell types involved in allergic inflammation and have been proposed to arise from a granulocyte-macrophage progenitor (GMP), that also generates neutrophils and macrophages, with GATA-2 up-regulation as a critical determinant. However, recent work identified a distinct basophil/eosinophil/mast cells progenitor (BEMP) that arises from multi-potent erythroid-primed multipotent progenitors (EMPP) where GATA-2 is already expressed ( Drissen et al., 2016). Here, we show that the GATA-2 interactome, rather than GATA-2 expression, controls myeloid lineage choice in EMPPs. To investigate the molecular mechanism instructing lineage decision in EMPP, we performed gene-expression profiling (RNA-seq) comparing its downstream progeny, early restricted erythroid/megakaryocyte progenitor (preMegE) and BEMP. This analysis identified differential expressed genes, including antagonist gene-expression pattern between the LIM-only protein 2 ( Lmo2) and protein 4 ( Lmo4), both previously shown as GATA-2 interactor. While Lmo2 function in hematopoiesis have been exhaustively studied, the role of Lmo4 remains unknown. In colony forming unit (CFU)-assay, overexpression of Lmo4 in EMPP mediated by lentiviral delivery, resulted in an increased number of mast-colonies associated with a decreased of erythroid/megakaryocyte colonies, while the total colony number was unchanged. In addition, transplantation of hematopoietic lineage-negative cells transduced with Lmo4/empty(Ctrl) lentiviral particles, showed an increased frequency of basophils and eosinophils associated with a decrease of erythroid cells in Lmo4 overexpressing cells compared to Ctrl, in both bone marrow and spleen of recipient mice, 8 weeks after transplantation. Analysis of transduced Lineage-negative Kit + (LK) compartment of recipient mice, also revealed an increased frequency of BEMP associated with a decreased frequency of preMegE in Lmo4overexpressing cells compared to Ctrl. Taken together, these data suggested that Lmo4 overexpression instruct lineage decision of EMPP toward basophil/eosinophil/mast lineages.Over-expression followed by co-immunoprecipitation in HEK-293T cell line revealed that LMO4 compete with LMO2 binding to GATA-2, suggesting a mutually exclusive transcriptional complex. More importantly, lack of interaction between LMO4 and GATA-2 mediated by a point mutation in the zinc-finger 1 domain (ZnF1) of GATA-2 (GATA2-G320D), frequently identified in AML/AEL, resulted in an absence of mast lineage potential of EMPP compared to GATA2 +/+. Moreover, over-expression of Lmo4 in EMPP from GATA2 G320D/G320D mice compared to GATA2 +/+, failed to instruct myeloid fate in vitro. Overall, our data identified a specific molecular mechanism of hematopoietic cell fate decision driven by protein-protein interaction with GATA-2. Perturbation of GATA-2 interactome through over-expression of the cofactor Lmo4 or specific point mutation in GATA2 ZnF1 domain results in lineage fate switch of EMPP toward the erythroid/megakaryocyte or basophil/eosinophil/mast lineages.

Functional Interrogation of Clinically Associated Venetoclax-Sensitive Genes Identifies an Essential Role of Transcription Factor ZNF740 in Therapeutic Response

Venetoclax, a targeted inhibitor of BCL-2, has shown therapeutic benefits in acute myeloid leukemia (AML). However, the emergence of drug resistance remains a significant clinical challenge, emphasizing the need for deeper understanding of its mechanisms. While comprehensive transcriptome analyses have identified many genes associated with venetoclax responsiveness, it remains crucial to distinguish between genes that play a direct role and those that are simply associated. In this study, we focused on genes strongly linked to a positive venetoclax response in primary AML samples using the BeatAML2.0 dataset. We identified ~1,400 top genes whose expression inversely correlated with venetoclax AUC. We theorized that if any of these genes were essential for venetoclax sensitivity, their absence would lead to resistance. To explore their significance, we developed a focused CRISPR knockout library of these ~1,400 genes, with each gene represented by 7 guide RNAs, and an additional 250 guide RNAs as negative controls. We then performed CRISPR screens on AML cell lines, using either DMSO or venetoclax as treatments. Interestingly, the majority of the gene knockouts did not lead to resistance, suggesting that these genes likely have an associative rather than a direct impact on venetoclax sensitivity. However, our screen did identify several genes that, when depleted, consistently resulted in resistance across different AML cell models. One of the top hits was the zinc finger transcription factor ZNF740, whose expression is strongly correlated with venetoclax AUC ( R = -0.52, p = 2.71e-26). For verification, we used CRISPR approach to knock out ZNF740 in Cas9-expressing OCI-AML2 and MOLM-13 cell lines. We then performed IC 50 and competitive proliferation assays under venetoclax treatment. We found that AML cells without ZNF740 exhibited a markedly higher IC 50 and demonstrated enhanced survival under venetoclax exposure. Conversely, enforced overexpression of ZNF740 rendered hypersensitivity to venetoclax. Intriguingly, our transcriptomic study showed that ZNF740 deletion upregulated monocyte differentiation signature genes, a change previously linked to venetoclax resistance. Collectively, our findings indicate that ZNF740 serves as a newly identified, clinically relevant transcriptional modulator of venetoclax responsiveness. We are now investigating the molecular targets of ZNF740, particularly its role in guiding AML lineage decisions and its impact on venetoclax response. Gaining insight into ZNF740's function and mechanism could pave the way for innovative strategies to counteract venetoclax resistance observed in clinical settings.

MMP13 promotes the osteogenic potential of BMP9 by enhancing Wnt/β-catenin signaling via HIF-1α upregulation in mouse embryonic fibroblasts.

Bone morphogenetic protein 9 (BMP9) has been validated as one of the most potent osteoinduction factors, but its underlying mechanism remains unclear. As a member of the matrix metalloproteinase (MMP) family, MMP13 may be involved in regulating the lineage-specific differentiation of mouse embryonic fibroblasts (MEFs). The goal of this study was to determine whether MMP13 regulates the osteoinduction potential of BMP9 in MEFs, which are multipotent progenitor cells widely used for stem cell biology research. In vitro and in vivo experiments showed that BMP9-induced osteogenic markers and/or bone were enhanced by exogenous MMP13 in MEFs, but were reduced by MMP13 knockdown or inhibition. The expression of hypoxia inducible factor 1 alpha (HIF-1α) was induced by BMP9, which was enhanced by MMP13. The protein expression of β-catenin and phosphorylation level of glycogen synthase kinase-3 beta (GSK-3β) were increased by BMP9 in MEFs, as was the translocation of β-catenin from the cytoplasm to the nucleus; all these effects of BMP9 were enhanced by MMP13. Furthermore, the MMP13 effects of increasing BMP9-induced β-catenin protein expression and GSK-3β phosphorylation level were partially reversed by HIF-1α knockdown. These results suggest that MMP13 can enhance the osteoinduction potential of BMP9, which may be mediated, at least in part, through the HIF-1α/β-catenin axis. Our findings demonstrate a novel role of MMP13 in the lineage decision of progenitor cells and provide a promising strategy to speed up bone regeneration.

Recognition and reconstruction of cell differentiation patterns with deep learning.

Cell lineage decisions occur in three-dimensional spatial patterns that are difficult to identify by eye. There is an ongoing effort to replicate such patterns using mathematical modeling. One approach uses long ranging cell-cell communication to replicate common spatial arrangements like checkerboard and engulfing patterns. In this model, the cell-cell communication has been implemented as a signal that disperses throughout the tissue. On the other hand, machine learning models have been developed for pattern recognition and pattern reconstruction tasks. We combined synthetic data generated by the mathematical model with spatial summary statistics and deep learning algorithms to recognize and reconstruct cell fate patterns in organoids of mouse embryonic stem cells. Application of Moran's index and pair correlation functions for in vitro and synthetic data from the model showed local clustering and radial segregation. To assess the patterns as a whole, a graph neural network was developed and trained on synthetic data from the model. Application to in vitro data predicted a low signal dispersion value. To test this result, we implemented a multilayer perceptron for the prediction of a given cell fate based on the fates of the neighboring cells. The results show a 70% accuracy of cell fate imputation based on the nine nearest neighbors of a cell. Overall, our approach combines deep learning with mathematical modeling to link cell fate patterns with potential underlying mechanisms.

Polycomb repressive complex 2 regulates basal cell fate during adult olfactory neurogenesis

Adult neurogenesis occurs in the mammalian olfactory epithelium to maintain populations of neurons that are vulnerable to injury yet essential for olfaction. Multipotent olfactory basal stem cells are activated by damage, although mechanisms regulating lineage decisions are not understood. Using mouse lesion models, we focused on defining the role of Polycomb repressive complexes (PRCs) in olfactory neurogenesis. PRC2 has a well-established role in developing tissues, orchestrating transcriptional programs via chromatin modification. PRC2 proteins are expressed in olfactory globose basal cells (GBCs) and nascent neurons. Conditional PRC2 loss perturbs lesion-induced neuron production, accompanied by altered histone modifications and misexpression of lineage-specific transcription factors in GBCs. De-repression of Sox9 in PRC2-mutant GBCs is accompanied by increased Bowman's gland production, defining an unrecognized role for PRC2 in regulating gland versus neuron cell fate. Our findings support a model for PRC2-dependent mechanisms promoting sensory neuronal differentiation in an adult neurogenic niche.

Silica Replicas Derived from Mammalian Cells as an Innovative Approach to Physically Direct Cell Lineage Decisions of Mesenchymal Stem Cells.

By means of a "live-cell" template strategy, silica replicas displaying the same morphology and topography of the mammalian cells used as templates are fabricated. The replicas are used as substrates to direct the differentiation of mesenchymal stem cells (MSCs) to predefined cell lineages. Upregulation of specific genes shows how the silica replica-based substrates have the ability to induce the molecular characteristics of the mature cell types from which they have been derived from. Thus, MSCs cultured in the presence of silica replicas of human osteoblasts (HObs) differentiate into HObs-like cells, as shown by the upregulation of specific osteogenic genes. Likewise, when MSCs are incubated with silica replicas derived from human chondrocytes, an enhanced expression of chondrogenic markers is observed. Importantly, the effects of the silica replicas are cell type-specific since the incubation of MSCs with HObs silica replicas does not result in upregulation of chondrogenic markers and vice versa. What is more, for both cases, the differentiation rate is enhanced when the silica replicas are used in combination with growth factors, suggesting a potential synergistic effect. These results demonstrate the potential of this innovative method as an efficient and cheap approach with the potential to eliminate, or at least reduce, the use of biochemically soluble compounds in stem cells research.

Lineage tracing identifies heterogeneous hepatoblast contribution to cell lineages and postembryonic organ growth dynamics.

To meet the physiological demands of the body, organs need to establish a functional tissue architecture and adequate size as the embryo develops to adulthood. In the liver, uni- and bipotent progenitor differentiation into hepatocytes and biliary epithelial cells (BECs), and their relative proportions, comprise the functional architecture. Yet, the contribution of individual liver progenitors at the organ level to both fates, and their specific proportion, is unresolved. Combining mathematical modelling with organ-wide, multispectral FRaeppli-NLS lineage tracing in zebrafish, we demonstrate that a precise BEC-to-hepatocyte ratio is established (i) fast, (ii) solely by heterogeneous lineage decisions from uni- and bipotent progenitors, and (iii) independent of subsequent cell type-specific proliferation. Extending lineage tracing to adulthood determined that embryonic cells undergo spatially heterogeneous three-dimensional growth associated with distinct environments. Strikingly, giant clusters comprising almost half a ventral lobe suggest lobe-specific dominant-like growth behaviours. We show substantial hepatocyte polyploidy in juveniles representing another hallmark of postembryonic liver growth. Our findings uncover heterogeneous progenitor contributions to tissue architecture-defining cell type proportions and postembryonic organ growth as key mechanisms forming the adult liver.

Harnessing spatial transcriptomics to study fate decisions and lineage specification in T cell development

Abstract T cells are critical components of the adaptive immune system; they develop from bone-marrow derived precursors that migrate to the thymus. There, precursors undergo a sequence of lineage specification and commitment steps that give rise to multiple T cell lineages serving essential functions in defenses against infection and immune homeostasis. The diversity of thymocytes, and their complex interactions with the thymic stroma, including epithelial cells and other hematopoietic compartments, have hampered our understanding of T cell development. For greater insight into intrathymic developmental trajectories and their transcriptional control, we and others have utilized single cell RNA- and ATAC-sequencing (scRNA-seq and scATAC-seq). However, single-cell approaches do not assess contextual cell-cell interactions. To address this limitation, we used in-situ transcriptomics, specifically the 10× Genomics “Visium” technology, which assays in-situ gene expression in 55-micron spots spanning whole mouse thymus lobes. We will present initial results of these experiments, including mapping of thymocyte differentiation stages, location of developmental bifurcations, and the spatial progression of lineage specific gene expression. Briefly, we use statistical models informed by scRNA-seq data to infer the location, relative to the thymic architecture, of lineage decisions. Additionally, we combine positional information of thymocytes and stromal cells to infer cell-cell interactions that direct lineage decisions. Collectively, this analysis sheds light on the interactions that guide development of the diverse T cell family from a common precursor and characterizes the journey of these cells to reach maturity.