Melanocyte lineage dynamics in development, growth and disease.

Melanocytes evolved to produce the melanin that gives colour to our hair, eyes and skin. The melanocyte lineage also gives rise to melanoma, the most lethal form of skin cancer. The melanocyte lineage differentiates from neural crest cells during development, and most melanocytes reside in the skin and hair, where they are replenished by melanocyte stem cells. Because the molecular mechanisms necessary for melanocyte specification, migration, proliferation and differentiation are co-opted during melanoma initiation and progression, studying melanocyte development is directly relevant to human disease. Here, through the lens of advances in cellular omic and genomic technologies, we review the latest findings in melanocyte development and differentiation, and how these developmental pathways become dysregulated in disease.

Hepatic stellate cells in zone 1 engage in capillarization rather than myofibroblast formation in murine liver fibrosis.

The combination of lineage tracing and immunohistochemistry has helped to identify subpopulations and fate of hepatic stellate cells (HSC) in murine liver. HSC are sinusoidal pericytes that act as myofibroblast precursors after liver injury. Single cell RNA sequencing approaches have recently helped to differentiate central and portal HSC. A specific Cre line to lineage trace portal HSC has not yet been described. We used three Cre lines (Lrat-Cre, PDGFRß-CreERT2 and SMMHC-CreERT2) known to label mesenchymal cells including HSC in combination with a tdTomato-expressing reporter. All three Cre lines labeled populations of HSC as well as smooth muscle cells (SMC). Using the SMMHC-CreERT2, we identified a subtype of HSC in the periportal area of the hepatic lobule (termed zone 1-HSC). We lineage traced tdTomato-expressing zone 1-HSC over 1 year, described fibrotic behavior in two fibrosis models and investigated their possible role during fibrosis. This HSC subtype resides in zone 1 under healthy conditions; however, zonation is disrupted in preclinical models of liver fibrosis (CCl4 and MASH). Zone 1-HSC do not transform into αSMA-expressing myofibroblasts. Rather, they participate in sinusoidal capillarization. We describe a novel subtype of HSC restricted to zone 1 under physiological conditions and its possible function after liver injury. In contrast to the accepted notion, this HSC subtype does not transform into αSMA-positive myofibroblasts; rather, zone 1-HSC adopt properties of capillary pericytes, thereby participating in sinusoidal capillarization.

Nanos2 marks precursors of somatic lineages and is required for germline formation in the sea anemone Nematostella vectensis.

In animals, stem cell populations of varying potency facilitate regeneration and tissue homeostasis. Notably, germline stem cells in both vertebrates and invertebrates express highly conserved RNA binding proteins, such as nanos, vasa, and piwi. In highly regenerative animals, these genes are also expressed in somatic stem cells, which led to the proposal that they had an ancestral role in all stem cells. In cnidarians, multi- and pluripotent interstitial stem cells have only been identified in hydrozoans. Therefore, it is currently unclear if cnidarian stem cell systems share a common evolutionary origin. We, therefore, aimed to characterize conserved stem cell marker genes in the sea anemone Nematostella vectensis. Through transgenic reporter genes and single-cell transcriptomics, we identify cell populations expressing the germline-associated markers piwi1 and nanos2 in the soma and germline, and gene knockout shows that Nanos2 is indispensable for germline formation. This suggests that nanos and piwi genes have a conserved role in somatic and germline stem cells in cnidarians.

A pre-vertebrate endodermal origin of calcitonin-producing neuroendocrine cells.

Vertebrate calcitonin-producing cells (C-cells) are neuroendocrine cells that secrete the small peptide hormone calcitonin in response to elevated blood calcium levels. Whereas mouse C-cells reside within the thyroid gland and derive from pharyngeal endoderm, avian C-cells are located within ultimobranchial glands and have been reported to derive from the neural crest. We use a comparative cell lineage tracing approach in a range of vertebrate model systems to resolve the ancestral embryonic origin of vertebrate C-cells. We find, contrary to previous studies, that chick C-cells derive from pharyngeal endoderm, with neural crest-derived cells instead contributing to connective tissue intimately associated with C-cells in the ultimobranchial gland. This endodermal origin of C-cells is conserved in a ray-finned bony fish (zebrafish) and a cartilaginous fish (the little skate, Leucoraja erinacea). Furthermore, we discover putative C-cell homologs within the endodermally-derived pharyngeal epithelium of the ascidian Ciona intestinalis and the amphioxus Branchiostoma lanceolatum, two invertebrate chordates that lack neural crest cells. Our findings point to a conserved endodermal origin of C-cells across vertebrates and to a pre-vertebrate origin of this cell type along the chordate stem.

Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes.

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs in NEPC, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

Reconstruction of single-cell lineage trajectories and identification of diversity in fates during the epithelial-to-mesenchymal transition.

Exploring the complexity of the epithelial-to-mesenchymal transition (EMT) unveils a diversity of potential cell fates; however, the exact timing and mechanisms by which early cell states diverge into distinct EMT trajectories remain unclear. Studying these EMT trajectories through single-cell RNA sequencing is challenging due to the necessity of sacrificing cells for each measurement. In this study, we employed optimal-transport analysis to reconstruct the past trajectories of different cell fates during TGF-beta-induced EMT in the MCF10A cell line. Our analysis revealed three distinct trajectories leading to low EMT, partial EMT, and high EMT states. Cells along the partial EMT trajectory showed substantial variations in the EMT signature and exhibited pronounced stemness. Throughout this EMT trajectory, we observed a consistent downregulation of the EED and EZH2 genes. This finding was validated by recent inhibitor screens of EMT regulators and CRISPR screen studies. Moreover, we applied our analysis of early-phase differential gene expression to gene sets associated with stemness and proliferation, pinpointing ITGB4, LAMA3, and LAMB3 as genes differentially expressed in the initial stages of the partial versus high EMT trajectories. We also found that CENPF, CKS1B, and MKI67 showed significant upregulation in the high EMT trajectory. While the first group of genes aligns with findings from previous studies, our work uniquely pinpoints the precise timing of these upregulations. Finally, the identification of the latter group of genes sheds light on potential cell cycle targets for modulating EMT trajectories.

Combinatorial regulatory states define cell fate diversity during embryogenesis.

Cell fate specification occurs along invariant species-specific trajectories that define the animal body plan. This process is controlled by gene regulatory networks that regulate the expression of the limited set of transcription factors encoded in animal genomes. Here we globally assess the spatial expression of ~90% of expressed transcription factors during sea urchin development from embryo to larva to determine the activity of gene regulatory networks and their regulatory states during cell fate specification. We show that >200 embryonically expressed transcription factors together define >70 cell fates that recapitulate the morphological and functional organization of this organism. Most cell fate-specific regulatory states consist of ~15-40 transcription factors with similarity particularly among functionally related cell types regardless of developmental origin. Temporally, regulatory states change continuously during development, indicating that progressive changes in regulatory circuit activity determine cell fate specification. We conclude that the combinatorial expression of transcription factors provides molecular definitions that suffice for the unique specification of cell states in time and space during embryogenesis.

A pathologically expanded, clonal lineage of IL-21-producing CD4+ T cells drives inflammatory neuropathy.

Inflammatory neuropathies, which include chronic inflammatory demyelinating polyneuropathy (CIDP) and Guillain Barré syndrome (GBS), result from autoimmune destruction of the PNS and are characterized by progressive weakness and sensory loss. CD4+ T cells play a key role in the autoimmune destruction of the PNS. Yet, key properties of pathogenic CD4+ T cells remain incompletely understood. Here, we used paired single-cell RNA-Seq (scRNA-Seq) and single-cell T cell receptor-sequencing (scTCR-Seq) of peripheral nerves from an inflammatory neuropathy mouse model to identify IL-21-expressing CD4+ T cells that were clonally expanded and multifunctional. These IL-21-expressing CD4+ T cells consisted of 2 transcriptionally distinct expanded cell populations, which expressed genes associated with T follicular helper (Tfh) and T peripheral helper (Tph) cell subsets. Remarkably, TCR clonotypes were shared between these 2 IL-21-expressing cell populations, suggesting a common lineage differentiation pathway. Finally, we demonstrated that IL-21 receptor-KO (IL-21R-KO) mice were protected from neuropathy development and had decreased immune infiltration into peripheral nerves. IL-21 signaling upregulated CXCR6, a chemokine receptor that promotes CD4+ T cell localization in peripheral nerves. Together, these findings point to IL-21 signaling, Tfh/Tph differentiation, and CXCR6-mediated cellular localization as potential therapeutic targets in inflammatory neuropathies.

Cellular proliferation biases clonal lineage tracing and trajectory inference.

MOTIVATION: Lineage tracing and trajectory inference from single-cell RNA-sequencing data hold tremendous potential for uncovering the genetic programs driving development and disease. Single cell datasets are thought to provide an unbiased view on the diverse cellular architecture of tissues. Sampling bias, however, can skew single cell datasets away from the cellular composition they are meant to represent. RESULTS: We demonstrate a novel form of sampling bias, caused by a statistical phenomenon related to repeated sampling from a growing, heterogeneous population. Relative growth rates of cells influence the probability that they will be sampled in clones observed across multiple time points. We support our probabilistic derivations with a simulation study and an analysis of a real time-course of T-cell development. We find that this bias can impact fate probability predictions, and we explore how to develop trajectory inference methods which are robust to this bias. AVAILABILITY AND IMPLEMENTATION: Source code for the simulated datasets and to create the figures in this manuscript is freely available in python at https://github.com/rbonhamcarter/simulate-clones. A python implementation of the extension of the LineageOT method is freely available at https://github.com/rbonhamcarter/LineageOT/tree/multi-time-clones.

Simulation of CRISPR-Cas9 editing on evolving barcode and accuracy of lineage tracing.

We designed a simulation program that mimics the CRISPR-Cas9 editing on evolving barcode and double strand break repair procedure along with cell divisions. Emerging barcode mutations tend to build upon previously existing mutations, occurring sequentially with each generation. This process results in a unique mutation profile in each cell. We sample the barcodes in leaf cells and reconstruct the lineage, comparing it to the original lineage tree to test algorithm accuracy under different parameter settings. Our computational simulations validate the reasonable assumptions deduced from experimental observations, emphasizing that factors such as sampling size, barcode length, multiple barcodes, indel probabilities, and Cas9 activity are critical for accurate and successful lineage tracing. Among the many factors we found that sampling size and indel probabilities are two major ones that affect lineage tracing accuracy. Large segment deletions in early generations could greatly impact lineage accuracy. These simulation results offer insightful recommendations for enhancing the design and analysis of Cas9-mediated molecular barcodes in actual experiments.

Spectrin mediates 3D-specific matrix stress-relaxation response in neural stem cell lineage commitment.

While extracellular matrix (ECM) stress relaxation is increasingly appreciated to regulate stem cell fate commitment and other behaviors, much remains unknown about how cells process stress-relaxation cues in tissue-like three-dimensional (3D) geometries versus traditional 2D cell culture. Here, we develop an oligonucleotide-crosslinked hyaluronic acid-based ECM platform with tunable stress relaxation properties capable of use in either 2D or 3D. Strikingly, stress relaxation favors neural stem cell (NSC) neurogenesis in 3D but suppresses it in 2D. RNA sequencing and functional studies implicate the membrane-associated protein spectrin as a key 3D-specific transducer of stress-relaxation cues. Confining stress drives spectrin's recruitment to the F-actin cytoskeleton, where it mechanically reinforces the cortex and potentiates mechanotransductive signaling. Increased spectrin expression is also accompanied by increased expression of the transcription factor EGR1, which we previously showed mediates NSC stiffness-dependent lineage commitment in 3D. Our work highlights spectrin as an important molecular sensor and transducer of 3D stress-relaxation cues.

Redox heterogeneity in mouse embryonic stem cells individualizes cell fate decisions.

Pluripotent embryonic stem cells (ESCs) can develop into any cell type in the body. Yet, the regulatory mechanisms that govern cell fate decisions during embryogenesis remain largely unknown. We now demonstrate that mouse ESCs (mESCs) display large natural variations in mitochondrial reactive oxygen species (mitoROS) levels that individualize their nuclear redox state, H3K4me3 landscape, and cell fate. While mESCs with high mitoROS levels (mitoROSHIGH) differentiate toward mesendoderm and form the primitive streak during gastrulation, mESCs, which generate less ROS, choose the alternative neuroectodermal fate. Temporal studies demonstrated that mesendodermal (ME) specification of mitoROSHIGH mESCs is mediated by a Nrf2-controlled switch in the nuclear redox state, triggered by the accumulation of redox-sensitive H3K4me3 marks, and executed by a hitherto unknown ROS-dependent activation process of the Wnt signaling pathway. In summary, our study explains how ESC heterogeneity is generated and used by individual cells to decide between distinct cellular fates.

A single-cell transcriptomic map of the developing Atoh1 lineage identifies neural fate decisions and neuronal diversity in the hindbrain.

Proneural transcription factors establish molecular cascades to orchestrate neuronal diversity. One such transcription factor, Atonal homolog 1 (Atoh1), gives rise to cerebellar excitatory neurons and over 30 distinct nuclei in the brainstem critical for hearing, breathing, and balance. Although Atoh1 lineage neurons have been qualitatively described, the transcriptional programs that drive their fate decisions and the full extent of their diversity remain unknown. Here, we analyzed single-cell RNA sequencing and ATOH1 DNA binding in Atoh1 lineage neurons of the developing mouse hindbrain. This high-resolution dataset identified markers for specific brainstem nuclei and demonstrated that transcriptionally heterogeneous progenitors require ATOH1 for proper migration. Moreover, we identified a sizable population of proliferating unipolar brush cell progenitors in the mouse Atoh1 lineage, previously described in humans as the origin of one medulloblastoma subtype. Collectively, our data provide insights into the developing mouse hindbrain and markers for functional assessment of understudied neuronal populations.

A developmental mechanism to regulate alternative polyadenylation in an adult stem cell lineage.

Alternative cleavage and polyadenylation (APA) often results in production of mRNA isoforms with either longer or shorter 3' UTRs from the same genetic locus, potentially impacting mRNA translation, localization, and stability. Developmentally regulated APA can thus make major contributions to cell type-specific gene expression programs as cells differentiate. During Drosophila spermatogenesis, ∼500 genes undergo APA when proliferating spermatogonia differentiate into spermatocytes, producing transcripts with shortened 3' UTRs, leading to profound stage-specific changes in the proteins expressed. The molecular mechanisms that specify usage of upstream polyadenylation sites in spermatocytes are thus key to understanding the changes in cell state. Here, we show that upregulation of PCF11 and Cbc, the two components of cleavage factor II (CFII), orchestrates APA during Drosophila spermatogenesis. Knockdown of PCF11 or cbc in spermatocytes caused dysregulation of APA, with many transcripts normally cleaved at a proximal site in spermatocytes now cleaved at their distal site, as in spermatogonia. Forced overexpression of CFII components in spermatogonia switched cleavage of some transcripts to the proximal site normally used in spermatocytes. Our findings reveal a developmental mechanism where changes in expression of specific cleavage factors can direct cell type-specific APA at selected genes.

Clonal analysis of fetal hematopoietic stem/progenitor cells reveals how post-transplantation capabilities are distributed.

It has been proposed that adult hematopoiesis is sustained by multipotent progenitors (MPPs) specified during embryogenesis. Adult-like hematopoietic stem cell (HSC) and MPP immunophenotypes are present in the fetus, but knowledge of their functional capacity is incomplete. We found that fetal MPP populations were functionally similar to adult cells, albeit with some differences in lymphoid output. Clonal assessment revealed that lineage biases arose from differences in patterns of single-/bi-lineage differentiation. Long-term (LT)- and short-term (ST)-HSC populations were distinguished from MPPs according to capacity for clonal multilineage differentiation. We discovered that a large cohort of long-term repopulating units (LT-RUs) resides within the ST-HSC population; a significant portion of these were labeled using Flt3-cre. This finding has two implications: (1) use of the CD150+ LT-HSC immunophenotype alone will significantly underestimate the size and diversity of the LT-RU pool and (2) LT-RUs in the ST-HSC population have the attributes required to persist into adulthood.

Dissecting reactive astrocyte responses: lineage tracing and morphology-based clustering.

Brain damage triggers diverse cellular and molecular events, with astrocytes playing a crucial role in activating local neuroprotective and reparative signaling within damaged neuronal circuits. Here, we investigated reactive astrocytes using a multidimensional approach to categorize their responses into different subtypes based on morphology. This approach utilized the StarTrack lineage tracer, single-cell imaging reconstruction and multivariate data analysis. Our findings identified three profiles of reactive astrocyte responses, categorized by their effects on cell size- and shape- related morphological parameters: "moderate", "strong," and "very strong". We also examined the heterogeneity of astrocyte reactivity, focusing on spatial and clonal distribution. Our research revealed a notable enrichment of protoplasmic and fibrous astrocytes within the "strong" and "very strong" response subtypes. Overall, our study contributes to a better understanding of astrocyte heterogeneity in response to an injury. By characterizing the diverse reactive responses among astrocyte subpopulations, we provide insights that could guide future research aimed at identifying novel therapeutic targets to mitigate brain damage and promote neural repair.

Emerging role of glutathione peroxidase 4 in myeloid cell lineage development and acute myeloid leukemia.

Phospholipid Hydroperoxide Gluthatione Peroxidase also called Glutathione Peroxidase 4 is one of the 25 described human selenoproteins. It plays an essential role in eliminating toxic lipid hydroxy peroxides, thus inhibiting ferroptosis and favoring cell survival. GPX4 is differentially expressed according to myeloid differentiation stage, exhibiting lower expression in hematopoietic stem cells and polymorphonuclear leucocytes, while harboring higher level of expression in common myeloid progenitors and monocytes. In addition, GPX4 is highly expressed in most of acute myeloid leukemia (AML) subtypes compared to normal hematopoietic stem cells. High GPX4 expression is consistently correlated to poor prognosis in patients suffering AML. However, the role of GPX4 in the development of the myeloid lineage and in the initiation and progression of myeloid leukemia remains poorly explored. Given its essential role in the detoxification of lipid hydroperoxides, and its overexpression in most of myeloid malignancies, GPX4 inhibition has emerged as a promising therapeutic strategy to specifically trigger ferroptosis and eradicate myeloid leukemia cells. In this review, we describe the most recent advances concerning the role of GPX4 and, more generally ferroptosis in the myeloid lineage and in the emergence of AML. We also discuss the therapeutic interest and limitations of GPX4 inhibition alone or in combination with other drugs as innovative therapies to treat AML patients.

Tumour vasculature at single-cell resolution.

Tumours can obtain nutrients and oxygen required to progress and metastasize through the blood supply1. Inducing angiogenesis involves the sprouting of established vessel beds and their maturation into an organized network2,3. Here we generate a comprehensive atlas of tumour vasculature at single-cell resolution, encompassing approximately 200,000 cells from 372 donors representing 31 cancer types. Trajectory inference suggested that tumour angiogenesis was initiated from venous endothelial cells and extended towards arterial endothelial cells. As neovascularization elongates (through angiogenic stages SI, SII and SIII), APLN+ tip cells at the SI stage (APLN+ TipSI) advanced to TipSIII cells with increased Notch signalling. Meanwhile, stalk cells, following tip cells, transitioned from high chemokine expression to elevated TEK (also known as Tie2) expression. Moreover, APLN+ TipSI cells not only were associated with disease progression and poor prognosis but also hold promise for predicting response to anti-VEGF therapy. Lymphatic endothelial cells demonstrated two distinct differentiation lineages: one responsible for lymphangiogenesis and the other involved in antigen presentation. In pericytes, endoplasmic reticulum stress was associated with the proangiogenic BASP1+ matrix-producing pericytes. Furthermore, intercellular communication analysis showed that neovascular endothelial cells could shape an immunosuppressive microenvironment conducive to angiogenesis. This study depicts the complexity of tumour vasculature and has potential clinical significance for anti-angiogenic therapy.

Transcriptional programs of Pitx2 and Tfap2a/Tfap2b controlling lineage specification of mandibular epithelium during tooth initiation.

How the dorsal-ventral axis of the vertebrate jaw, particularly the position of tooth initiation site, is established remains a critical and unresolved question. Tooth development starts with the formation of the dental lamina, a localized thickened strip within the maxillary and mandibular epithelium. To identify transcriptional regulatory networks (TRN) controlling the specification of dental lamina from the naïve mandibular epithelium, we utilized Laser Microdissection coupled low-input RNA-seq (LMD-RNA-seq) to profile gene expression of different domains of the mandibular epithelium along the dorsal-ventral axis. We comprehensively identified transcription factors (TFs) and signaling pathways that are differentially expressed along mandibular epithelial domains (including the dental lamina). Specifically, we found that the TFs Sox2 and Tfap2 (Tfap2a/Tfap2b) formed complimentary expression domains along the dorsal-ventral axis of the mandibular epithelium. Interestingly, both classic and novel dental lamina specific TFs-such as Pitx2, Ascl5 and Zfp536-were found to localize near the Sox2:Tfap2a/Tfap2b interface. To explore the functional significance of these domain specific TFs, we next examined loss-of-function mouse models of these domain specific TFs, including the dental lamina specific TF, Pitx2, and the ventral surface ectoderm specific TFs Tfap2a and Tfap2b. We found that disruption of domain specific TFs leads to an upregulation and expansion of the alternative domain's TRN. The importance of this cross-repression is evident by the ectopic expansion of Pitx2 and Sox2 positive dental lamina structure in Tfap2a/Tfap2b ectodermal double knockouts and the emergence of an ectopic tooth in the ventral surface ectoderm. Finally, we uncovered an unappreciated interface of mesenchymal SHH and WNT signaling pathways, at the site of tooth initiation, that were established by the epithelial domain specific TFs including Pitx2 and Tfap2a/Tfap2b. These results uncover a previously unknown molecular mechanism involving cross-repression of domain specific TFs including Pitx2 and Tfap2a/Tfap2b in patterning the dorsal-ventral axis of the mouse mandible, specifically the regulation of tooth initiation site.

Characterizing Rare Dormant and Cycling Lineages Using the Watermelon System.

Barcode-based lineage tracing approaches enable molecular characterization of clonal cell families. Barcodes that are expressed as mRNA can be used to deconvolve lineage identity from single-cell RNA sequencing transcriptional data. Here we describe the Watermelon system, which facilitates the simultaneous tracing of lineage, transcriptional, and proliferative state at a single cell level.