Single-cell mapping of lineage and identity in direct reprogramming

Functional characterization of all CDKN2A missense variants and comparison to in silico models of pathogenicity

Interpretation of variants identified during genetic testing is a significant clinical challenge. In this study, we developed a high-throughput CDKN2A functional assay and characterized all possible CDKN2A missense variants. We found that 17.7% of all missense variants were functionally deleterious. We also used our functional classifications to assess the performance of in silico models that predict the effect of variants, including recently reported models based on machine learning. Notably, we found that all in silico models performed similarly when compared to our functional classifications with accuracies of 39.5-85.4%. Furthermore, while we found that functionally deleterious variants were enriched within ankyrin repeats, we did not identify any residues where all missense variants were functionally deleterious. Our functional classifications are a resource to aid the interpretation of CDKN2A variants and have important implications for the application of variant interpretation guidelines, particularly the use of in silico models for clinical variant interpretation.

Revealing gene function with statistical inference at single-cell resolution

Single-cell and spatial molecular profiling assays have shown large gains in sensitivity, resolution and throughput. Applying these technologies to specimens from human and model organisms promises to comprehensively catalogue cell types, reveal their lineage origins in development and discern their contributions to disease pathogenesis. Moreover, rapidly dropping costs have made well-controlled perturbation experiments and cohort studies widely accessible, illuminating mechanisms that give rise to phenotypes at the scale of the cell, the tissue and the whole organism. Interpreting the coming flood of single-cell data, much of which will be spatially resolved, will place a tremendous burden on existing computational pipelines. However, statistical concepts, models, tools and algorithms can be repurposed to solve problems now arising in genetic and molecular biology studies of development and disease. Here, I review how the questions that recent technological innovations promise to answer can be addressed by the major classes of statistical tools.

Multi-omic lineage tracing predicts the transcriptional, epigenetic and genetic determinants of cancer evolution

Cancer is a highly heterogeneous disease, where phenotypically distinct subpopulations coexist and can be primed to different fates. Both genetic and epigenetic factors may drive cancer evolution, however little is known about whether and how such a process is pre-encoded in cancer clones. Using single-cell multi-omic lineage tracing and phenotypic assays, we investigate the predictive features of either tumour initiation or drug tolerance within the same cancer population. Clones primed to tumour initiation in vivo display two distinct transcriptional states at baseline. Remarkably, these states share a distinctive DNA accessibility profile, highlighting an epigenetic basis for tumour initiation. The drug tolerant niche is also largely pre-encoded, but only partially overlaps the tumour-initiating one and evolves following two genetically and transcriptionally distinct trajectories. Our study highlights coexisting genetic, epigenetic and transcriptional determinants of cancer evolution, unravelling the molecular complexity of pre-encoded tumour phenotypes.

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.

Glioblastoma Neurovascular Progenitor Orchestrates Tumor Cell Type Diversity

Glioblastoma (GBM) is the deadliest form of primary brain tumor with limited treatment options. Recent studies have profiled GBM tumor heterogeneity, revealing numerous axes of variation that explain the molecular and spatial features of the tumor. Here, we seek to bridge descriptive characterization of GBM cell type heterogeneity with the functional role of individual populations within the tumor. Our lens leverages a gene program-centric meta-atlas of published transcriptomic studies to identify commonalities between diverse tumors and cell types in order to decipher the mechanisms that drive them. This approach led to the discovery of a tumor-derived stem cell population with mixed vascular and neural stem cell features, termed a neurovascular progenitor (NVP). Following in situ validation and molecular characterization of NVP cells in GBM patient samples, we characterized their function in vivo. Genetic depletion of NVP cells resulted in altered tumor cell composition, fewer cycling cells, and extended survival, underscoring their critical functional role. Clonal analysis of primary patient tumors in a human organoid tumor transplantation system demonstrated that the NVP has dual potency, generating both neuronal and vascular tumor cells. Although NVP cells comprise a small fraction of the tumor, these clonal analyses demonstrated that they strongly contribute to the total number of cycling cells in the tumor and generate a defined subset of the whole tumor. This study represents a paradigm by which cell type-specific interrogation of tumor populations can be used to study functional heterogeneity and therapeutically targetable vulnerabilities of GBM.

Synthetic DNA barcodes identify singlets in scRNA-seq datasets and evaluate doublet algorithms

Single-cell RNA sequencing (scRNA-seq) datasets contain true single cells, or singlets, in addition to cells that coalesce during the protocol, or doublets. Identifying singlets with high fidelity in scRNA-seq is necessary to avoid false negative and false positive discoveries. Although several methodologies have been proposed, they are typically tested on highly heterogeneous datasets and lack a priori knowledge of true singlets. Here, we leveraged datasets with synthetically introduced DNA barcodes for a hitherto unexplored application: to extract ground-truth singlets. We demonstrated the feasibility of our framework, "singletCode," to evaluate existing doublet detection methods across a range of contexts. We also leveraged our ground-truth singlets to train a proof-of-concept machine learning classifier, which outperformed other doublet detection algorithms. Our integrative framework can identify ground-truth singlets and enable robust doublet detection in non-barcoded datasets.

Multiscale engineering of brain organoids for disease modeling

Brain organoids hold great potential for modeling human brain development and pathogenesis. They recapitulate certain aspects of the transcriptional trajectory, cellular diversity, tissue architecture and functions of the developing brain. In this review, we explore the engineering strategies to control the molecular-, cellular- and tissue-level inputs to achieve high-fidelity brain organoids. We review the application of brain organoids in neural disorder modeling and emerging bioengineering methods to improve data collection and feature extraction at multiscale. The integration of multiscale engineering strategies and analytical methods has significant potential to advance insight into neurological disorders and accelerate drug development.

CellRank 2: unified fate mapping in multiview single-cell data

Single-cell RNA sequencing allows us to model cellular state dynamics and fate decisions using expression similarity or RNA velocity to reconstruct state-change trajectories; however, trajectory inference does not incorporate valuable time point information or utilize additional modalities, whereas methods that address these different data views cannot be combined or do not scale. Here we present CellRank 2, a versatile and scalable framework to study cellular fate using multiview single-cell data of up to millions of cells in a unified fashion. CellRank 2 consistently recovers terminal states and fate probabilities across data modalities in human hematopoiesis and endodermal development. Our framework also allows combining transitions within and across experimental time points, a feature we use to recover genes promoting medullary thymic epithelial cell formation during pharyngeal endoderm development. Moreover, we enable estimating cell-specific transcription and degradation rates from metabolic-labeling data, which we apply to an intestinal organoid system to delineate differentiation trajectories and pinpoint regulatory strategies.

LMNA -Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-derived iPSC Differentiation Support Cell type and Lineage-specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA -Related Dilated Cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C ( LMNA ) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. Molecular mechanisms of disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA -Related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (4 Patient and 8 Control) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for Cardiac Progenitors to Cardiomyocytes (CM) and Epicardium-Derived Cells (EPDC). Data integration and comparative analyses of Patient and Control cells found cell type and lineage differentially expressed genes (DEG) with enrichment to support pathway dysregulation. Top DEG and enriched pathways included: 10 ZNF genes and RNA polymerase II transcription in Pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CM; LMNA and epigenetic regulation and DDIT4 and mTORC1 signaling in EPDC. Top DEG also included: XIST and other X-linked genes, six imprinted genes: SNRPN , PWAR6 , NDN , PEG10 , MEG3 , MEG8 , and enriched gene sets in metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

A temporal developmental map separates human NK cells from noncytotoxic ILCs through clonal and single-cell analysis

ABSTRACT

Natural killer (NK) cells represent the cytotoxic member within the innate lymphoid cell (ILC) family that are important against viral infections and cancer. Although the NK cell emergence from hematopoietic stem and progenitor cells through multiple intermediate stages and the underlying regulatory gene network has been extensively studied in mice, this process is not well characterized in humans. Here, using a temporal in vitro model to reconstruct the developmental trajectory of NK lineage, we identified an ILC-restricted oligopotent stage 3a CD34-CD117+CD161+CD45RA+CD56- progenitor population, that exclusively gave rise to CD56-expressing ILCs in vitro. We also further investigated a previously nonappreciated heterogeneity within the CD56+CD94-NKp44+ subset, phenotypically equivalent to stage 3b population containing both group-1 ILC and RORγt+ ILC3 cells, that could be further separated based on their differential expression of DNAM-1 and CD161 receptors. We confirmed that DNAM-1hi S3b and CD161hiCD117hi ILC3 populations distinctively differed in their expression of effector molecules, cytokine secretion, and cytotoxic activity. Furthermore, analysis of lineage output using DNA-barcode tracing across these stages supported a close developmental relationship between S3b-NK and S4-NK (CD56+CD94+) cells, whereas distant to the ILC3 subset. Cross-referencing gene signatures of culture-derived NK cells and other noncytotoxic ILCs with publicly available data sets validated that these in vitro stages highly resemble transcriptional profiles of respective in vivo ILC counterparts. Finally, by integrating RNA velocity and gene network analysis through single-cell regulatory network inference and clustering we unravel a network of coordinated and highly dynamic regulons driving the cytotoxic NK cell program, as a guide map for future studies on NK cell regulation.

Olfactory neuroblastoma mimics molecular heterogeneity and lineage trajectories of small-cell lung cancer

The olfactory epithelium undergoes neuronal regeneration from basal stem cells and is susceptible to olfactory neuroblastoma (ONB), a rare tumor of unclear origins. Employing alterations in Rb1/Trp53/Myc (RPM), we establish a genetically engineered mouse model of high-grade metastatic ONB exhibiting a NEUROD1+ immature neuronal phenotype. We demonstrate that globose basal cells (GBCs) are a permissive cell of origin for ONB and that ONBs exhibit cell fate heterogeneity that mimics normal GBC developmental trajectories. ASCL1 loss in RPM ONB leads to emergence of non-neuronal histopathologies, including a POU2F3+ microvillar-like state. Similar to small-cell lung cancer (SCLC), mouse and human ONBs exhibit mutually exclusive NEUROD1 and POU2F3-like states, an immune-cold tumor microenvironment, intratumoral cell fate heterogeneity comprising neuronal and non-neuronal lineages, and cell fate plasticity-evidenced by barcode-based lineage tracing and single-cell transcriptomics. Collectively, our findings highlight conserved similarities between ONB and neuroendocrine tumors with significant implications for ONB classification and treatment.

A molecular and cellular perspective on human brain evolution and tempo

The evolution of the modern human brain was accompanied by distinct molecular and cellular specializations, which underpin our diverse cognitive abilities but also increase our susceptibility to neurological diseases. These features, some specific to humans and others shared with related species, manifest during different stages of brain development. In this multi-stage process, neural stem cells proliferate to produce a large and diverse progenitor pool, giving rise to excitatory or inhibitory neurons that integrate into circuits during further maturation. This process unfolds over varying time scales across species and has progressively become slower in the human lineage, with differences in tempo correlating with differences in brain size, cell number and diversity, and connectivity. Here we introduce the terms 'bradychrony' and 'tachycrony' to describe slowed and accelerated developmental tempos, respectively. We review how recent technical advances across disciplines, including advanced engineering of in vitro models, functional comparative genetics and high-throughput single-cell profiling, are leading to a deeper understanding of how specializations of the human brain arise during bradychronic neurodevelopment. Emerging insights point to a central role for genetics, gene-regulatory networks, cellular innovations and developmental tempo, which together contribute to the establishment of human specializations during various stages of neurodevelopment and at different points in evolution.

Clonal tracking in cancer and metastasis

The eradication of many cancers has proven challenging due to the presence of functionally and genetically heterogeneous clones maintained by rare cancer stem cells (CSCs), which contribute to disease progression, treatment refractoriness, and late relapse. The characterization of functional CSC activity has necessitated the development of modern clonal tracking strategies. This review describes viral-based and CRISPR-Cas9-based cellular barcoding, lineage tracing, and imaging-based approaches. DNA-based cellular barcoding technology is emerging as a powerful and robust strategy that has been widely applied to in vitro and in vivo model systems, including patient-derived xenograft models. This review also highlights the potential of these methods for use in the clinical and drug discovery contexts and discusses the important insights gained from such approaches.

Single-cell lineage capture across genomic modalities with CellTag-multi reveals fate-specific gene regulatory changes

Complex gene regulatory mechanisms underlie differentiation and reprogramming. Contemporary single-cell lineage-tracing (scLT) methods use expressed, heritable DNA barcodes to combine cell lineage readout with single-cell transcriptomics. However, reliance on transcriptional profiling limits adaptation to other single-cell assays. With CellTag-multi, we present an approach that enables direct capture of heritable random barcodes expressed as polyadenylated transcripts, in both single-cell RNA sequencing and single-cell Assay for Transposase Accessible Chromatin using sequencing assays, allowing for independent clonal tracking of transcriptional and epigenomic cell states. We validate CellTag-multi to characterize progenitor cell lineage priming during mouse hematopoiesis. Additionally, in direct reprogramming of fibroblasts to endoderm progenitors, we identify core regulatory programs underlying on-target and off-target fates. Furthermore, we reveal the transcription factor Zfp281 as a regulator of reprogramming outcome, biasing cells toward an off-target mesenchymal fate. Our results establish CellTag-multi as a lineage-tracing method compatible with multiple single-cell modalities and demonstrate its utility in revealing fate-specifying gene regulatory changes across diverse paradigms of differentiation and reprogramming.

Multiplexed single-cell lineage tracing of mitotic kinesin inhibitor resistance in glioblastoma

Glioblastoma (GBM) is a deadly brain tumor, and the kinesin motor KIF11 is an attractive therapeutic target with roles in proliferation and invasion. Resistance to KIF11 inhibitors, which has mainly been studied in animal models, presents significant challenges. We use lineage-tracing barcodes and single-cell RNA sequencing to analyze resistance in patient-derived GBM neurospheres treated with ispinesib, a potent KIF11 inhibitor. Similar to GBM progression in patients, untreated cells lose their neural lineage identity and become mesenchymal, which is associated with poor prognosis. Conversely, cells subjected to long-term ispinesib treatment exhibit a proneural phenotype. We generate patient-derived xenografts and show that ispinesib-resistant cells form less aggressive tumors in vivo, even in the absence of drug. Moreover, treatment of human ex vivo GBM slices with ispinesib demonstrates phenotypic alignment with in vitro responses, underscoring the clinical relevance of our findings. Finally, using retrospective lineage tracing, we identify drugs that are synergistic with ispinesib.

Analysis of synthetic cellular barcodes in the genome and transcriptome with BARtab and bartools

Cellular barcoding is a lineage-tracing methodology that couples heritable synthetic barcodes to high-throughput sequencing, enabling the accurate tracing of cell lineages across a range of biological contexts. Recent studies have extended these methods by incorporating lineage information into single-cell or spatial transcriptomics readouts. Leveraging the rich biological information within these datasets requires dedicated computational tools for dataset pre-processing and analysis. Here, we present BARtab, a portable and scalable Nextflow pipeline, and bartools, an open-source R package, designed to provide an integrated end-to-end cellular barcoding analysis toolkit. BARtab and bartools contain methods to simplify the extraction, quality control, analysis, and visualization of lineage barcodes from population-level, single-cell, and spatial transcriptomics experiments. We showcase the utility of our integrated BARtab and bartools workflow via the analysis of exemplar bulk, single-cell, and spatial transcriptomics experiments containing cellular barcoding information.

Novel cell states arise in embryonic cells devoid of key reprogramming factors

The capacity for embryonic cells to differentiate relies on a large-scale reprogramming of the oocyte and sperm nucleus into a transient totipotent state. In zebrafish, this reprogramming step is achieved by the pioneer factors Nanog, Pou5f3, and Sox19b (NPS). Yet, it remains unclear whether cells lacking this reprogramming step are directed towards wild type states or towards novel developmental canals in the Waddington landscape of embryonic development. Here we investigate the developmental fate of embryonic cells mutant for NPS by analyzing their single-cell gene expression profiles. We find that cells lacking the first developmental reprogramming steps can acquire distinct cell states. These states are manifested by gene expression modules that result from a failure of nuclear reprogramming, the persistence of the maternal program, and the activation of somatic compensatory programs. As a result, most mutant cells follow new developmental canals and acquire new mixed cell states in development. In contrast, a group of mutant cells acquire primordial germ cell-like states, suggesting that NPS-dependent reprogramming is dispensable for these cell states. Together, these results demonstrate that developmental reprogramming after fertilization is required to differentiate most canonical developmental programs, and loss of the transient totipotent state canalizes embryonic cells into new developmental states in vivo.

Spatiotemporally resolved colorectal oncogenesis in mini-colons ex vivo

Three-dimensional organoid culture technologies have revolutionized cancer research by allowing for more realistic and scalable reproductions of both tumour and microenvironmental structures1-3. This has enabled better modelling of low-complexity cancer cell behaviours that occur over relatively short periods of time4. However, available organoid systems do not capture the intricate evolutionary process of cancer development in terms of tissue architecture, cell diversity, homeostasis and lifespan. As a consequence, oncogenesis and tumour formation studies are not possible in vitro and instead require the extensive use of animal models, which provide limited spatiotemporal resolution of cellular dynamics and come at a considerable cost in terms of resources and animal lives. Here we developed topobiologically complex mini-colons that are able to undergo tumorigenesis ex vivo by integrating microfabrication, optogenetic and tissue engineering approaches. With this system, tumorigenic transformation can be spatiotemporally controlled by directing oncogenic activation through blue-light exposure, and emergent colon tumours can be tracked in real-time at the single-cell resolution for several weeks without breaking the culture. These induced mini-colons display rich intratumoural and intertumoural diversity and recapitulate key pathophysiological hallmarks displayed by colorectal tumours in vivo. By fine-tuning cell-intrinsic and cell-extrinsic parameters, mini-colons can be used to identify tumorigenic determinants and pharmacological opportunities. As a whole, our study paves the way for cancer initiation research outside living organisms.

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for studying cellular differentiation, but accurately tracking cell fate transitions can be challenging, especially in disease conditions. Here we introduce PhyloVelo, a computational framework that estimates the velocity of transcriptomic dynamics by using monotonically expressed genes (MEGs) or genes with expression patterns that either increase or decrease, but do not cycle, through phylogenetic time. Through integration of scRNA-seq data with lineage information, PhyloVelo identifies MEGs and reconstructs a transcriptomic velocity field. We validate PhyloVelo using simulated data and Caenorhabditis elegans ground truth data, successfully recovering linear, bifurcated and convergent differentiations. Applying PhyloVelo to seven lineage-traced scRNA-seq datasets, generated using CRISPR-Cas9 editing, lentiviral barcoding or immune repertoire profiling, demonstrates its high accuracy and robustness in inferring complex lineage trajectories while outperforming RNA velocity. Additionally, we discovered that MEGs across tissues and organisms share similar functions in translation and ribosome biogenesis.

Lineage Reprogramming: Genetic, Chemical, and Physical Cues for Cell Fate Conversion with a Focus on Neuronal Direct Reprogramming and Pluripotency Reprogramming

Although lineage reprogramming from one cell type to another is becoming a breakthrough technology for cell-based therapy, several limitations remain to be overcome, including the low conversion efficiency and subtype specificity. To address these, many studies have been conducted using genetics, chemistry, physics, and cell biology to control transcriptional networks, signaling cascades, and epigenetic modifications during reprogramming. Here, we summarize recent advances in cellular reprogramming and discuss future directions.

The implications of single-cell RNA-seq analysis in prostate cancer: unraveling tumor heterogeneity, therapeutic implications and pathways towards personalized therapy

In recent years, advancements in single-cell and spatial transcriptomics, which are highly regarded developments in the current era, particularly the emerging integration of single-cell and spatiotemporal transcriptomics, have enabled a detailed molecular comprehension of the complex regulation of cell fate. The insights obtained from these methodologies are anticipated to significantly contribute to the development of personalized medicine. Currently, single-cell technology is less frequently utilized for prostate cancer compared with other types of tumors. Starting from the perspective of RNA sequencing technology, this review outlined the significance of single-cell RNA sequencing (scRNA-seq) in prostate cancer research, encompassing preclinical medicine and clinical applications. We summarize the differences between mouse and human prostate cancer as revealed by scRNA-seq studies, as well as a combination of multi-omics methods involving scRNA-seq to highlight the key molecular targets for the diagnosis, treatment, and drug resistance characteristics of prostate cancer. These studies are expected to provide novel insights for the development of immunotherapy and other innovative treatment strategies for castration-resistant prostate cancer. Furthermore, we explore the potential clinical applications stemming from other single-cell technologies in this review, paving the way for future research in precision medicine.

The Methyltransferases METTL7A and METTL7B Confer Resistance to Thiol-Based Histone Deacetylase Inhibitors

Histone deacetylase inhibitors (HDACi) are part of a growing class of epigenetic therapies used for the treatment of cancer. Although HDACis are effective in the treatment of T-cell lymphomas, treatment of solid tumors with this class of drugs has not been successful. Overexpression of the multidrug resistance protein P-glycoprotein (P-gp), encoded by ABCB1, is known to confer resistance to the HDACi romidepsin in vitro, yet increased ABCB1 expression has not been associated with resistance in patients, suggesting that other mechanisms of resistance arise in the clinic. To identify alternative mechanisms of resistance to romidepsin, we selected MCF-7 breast cancer cells with romidepsin in the presence of the P-gp inhibitor verapamil to reduce the likelihood of P-gp-mediated resistance. The resulting cell line, MCF-7 DpVp300, does not express P-gp and was found to be selectively resistant to romidepsin but not to other HDACis such as belinostat, panobinostat, or vorinostat. RNA-sequencing analysis revealed upregulation of the mRNA coding for the putative methyltransferase, METTL7A, whose paralog, METTL7B, was previously shown to methylate thiol groups on hydrogen sulfide and captopril. As romidepsin has a thiol as the zinc-binding moiety, we hypothesized that METTL7A could inactivate romidepsin and other thiol-based HDACis via methylation of the thiol group. We demonstrate that expression of METTL7A or METTL7B confers resistance to thiol-based HDACis and that both enzymes are capable of methylating thiol-containing HDACis. We thus propose that METTL7A and METTL7B confer resistance to thiol-based HDACis by methylating and inactivating the zinc-binding thiol.

Principles of digital sequencing using unique molecular identifiers

Massively parallel sequencing technologies have long been used in both basic research and clinical routine. The recent introduction of digital sequencing has made previously challenging applications possible by significantly improving sensitivity and specificity to now allow detection of rare sequence variants, even at single molecule level. Digital sequencing utilizes unique molecular identifiers (UMIs) to minimize sequencing-induced errors and quantification biases. Here, we discuss the principles of UMIs and how they are used in digital sequencing. We outline the properties of different UMI types and the consequences of various UMI approaches in relation to experimental protocols and bioinformatics. Finally, we describe how digital sequencing can be applied in specific research fields, focusing on cancer management where it can be used in screening of asymptomatic individuals, diagnosis, treatment prediction, prognostication, monitoring treatment efficacy and early detection of treatment resistance as well as relapse.

Gene-expression memory-based prediction of cell lineages from scRNA-seq datasets

Assigning single cell transcriptomes to cellular lineage trees by lineage tracing has transformed our understanding of differentiation during development, regeneration, and disease. However, lineage tracing is technically demanding, often restricted in time-resolution, and most scRNA-seq datasets are devoid of lineage information. Here we introduce Gene Expression Memory-based Lineage Inference (GEMLI), a computational tool allowing to robustly identify small to medium-sized cell lineages solely from scRNA-seq datasets. GEMLI allows to study heritable gene expression, to discriminate symmetric and asymmetric cell fate decisions and to reconstruct individual multicellular structures from pooled scRNA-seq datasets. In human breast cancer biopsies, GEMLI reveals previously unknown gene expression changes at the onset of cancer invasiveness. The universal applicability of GEMLI allows studying the role of small cell lineages in a wide range of physiological and pathological contexts, notably in vivo. GEMLI is available as an R package on GitHub ( https://github.com/UPSUTER/GEMLI ).

Clonal differences underlie variable responses to sequential and prolonged treatment

Cancer cells exhibit dramatic differences in gene expression at the single-cell level, which can predict whether they become resistant to treatment. Treatment perpetuates this heterogeneity, resulting in a diversity of cell states among resistant clones. However, it remains unclear whether these differences lead to distinct responses when another treatment is applied or the same treatment is continued. In this study, we combined single-cell RNA sequencing with barcoding to track resistant clones through prolonged and sequential treatments. We found that cells within the same clone have similar gene expression states after multiple rounds of treatment. Moreover, we demonstrated that individual clones have distinct and differing fates, including growth, survival, or death, when subjected to a second treatment or when the first treatment is continued. By identifying gene expression states that predict clone survival, this work provides a foundation for selecting optimal therapies that target the most aggressive resistant clones within a tumor. A record of this paper's transparent peer review process is included in the supplemental information.

Protocol to study the inheritance and propagation of non-genetically encoded states using barcode decay lineage tracing

Here, we present a protocol to perform barcode decay lineage tracing followed by single-cell transcriptome analysis (BdLT-Seq). We describe steps for BdLT-Seq experimental design, building barcoded episome reporters, performing episome transfection, and barcode retrieval. We then describe procedures for sequencing library construction while providing options for sample multiplexing and data analysis. This BdLT-Seq technique enables the assessment of clonal evolution in a directional manner while preserving isogeneity, thus allowing the comparison of non-genetic molecular features between isogenic cell lineages. For complete details on the use and execution of this protocol, please refer to Shlyakhtina et al. (2023).1.

Ex vivo expansion potential of murine hematopoietic stem cells is a rare property only partially predicted by phenotype

The scarcity of hematopoietic stem cells (HSCs) restricts their use in both clinical settings and experimental research. Here, we examined a recently developed method for expanding rigorously purified murine HSCs ex vivo. After 3 weeks of culture, only 0.1% of cells exhibited the input HSC phenotype, but these accounted for almost all functional long-term HSC activity. Input HSCs displayed varying potential for ex vivo self-renewal, with alternative outcomes revealed by single-cell multimodal RNA and ATAC sequencing profiling. While most HSC progeny offered only transient in vivo reconstitution, these cells efficiently rescued mice from lethal myeloablation. The amplification of functional HSC activity allowed for long-term multilineage engraftment in unconditioned hosts that associated with a return of HSCs to quiescence. Thereby, our findings identify several key considerations for ex vivo HSC expansion, with major implications also for assessment of normal HSC activity.

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.

Cell fate decision in erythropoiesis: Insights from multiomics studies

Every second, the body produces 2 million red blood cells through a process called erythropoiesis. Erythropoiesis is hierarchical in that it results from a series of cell fate decisions whereby hematopoietic stem cells progress toward the erythroid lineage. Single-cell transcriptomic and proteomic approaches have revolutionized the way we understand erythropoiesis, revealing it to be a gradual process that underlies a progressive restriction of fate potential driven by quantitative changes in lineage-specifying transcription factors. Despite these major advances, we still know very little about what cell fate decision entails at the molecular level. Novel approaches that simultaneously measure additional properties in single cells, including chromatin accessibility, transcription factor binding, and/or cell surface proteins are being developed at a fast pace, providing the means to exciting new advances in the near future. In this review, we briefly summarize the main findings obtained from single-cell studies of erythropoiesis, highlight outstanding questions, and suggest recent technological advances to address them.

Retrospective identification of cell-intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of OCT4, KLF4, SOX2, and MYC. Usually only a rare subset of cells reprogram, and the molecular characteristics of this subset remain unknown. We apply retrospective clone tracing to identify and characterize the rare human fibroblasts primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis increased the reprogramming efficiency. We provide evidence for a unified model in which cells can move into and out of the primed state over time, explaining how reprogramming appears deterministic at short timescales and stochastic at long timescales. Furthermore, inhibiting the activity of LSD1 enlarged the pool of cells that were primed for reprogramming. Thus, even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Tracing back primed resistance in cancer via sister cells

Exploring non-genetic evolution of cell states during cancer treatments has become attainable by recent advances in lineage-tracing methods. However, transcriptional changes that drive cells into resistant fates may be subtle, necessitating high resolution analysis. Here, we present ReSisTrace that uses shared transcriptomic features of sister cells to predict the states priming treatment resistance. Applying ReSisTrace in ovarian cancer cells perturbed with olaparib, carboplatin or natural killer (NK) cells reveals pre-resistant phenotypes defined by proteostatic and mRNA surveillance features, reflecting traits enriched in the upcoming subclonal selection. Furthermore, we show that DNA repair deficiency renders cells susceptible to both DNA damaging agents and NK killing in a context-dependent manner. Finally, we leverage the obtained pre-resistance profiles to predict and validate small molecules driving cells to sensitive states prior to treatment. In summary, ReSisTrace resolves pre-existing transcriptional features of treatment vulnerability, facilitating both molecular patient stratification and discovery of synergistic pre-sensitizing therapies.

Extracting, filtering and simulating cellular barcodes using CellBarcode tools

Identifying true DNA cellular barcodes among polymerase chain reaction and sequencing errors is challenging. Current tools are restricted in the diversity of barcode types supported or the analysis strategies implemented. As such, there is a need for more versatile and efficient tools for barcode extraction, as well as for tools to investigate which factors impact barcode detection and which filtering strategies to best apply. Here we introduce the package CellBarcode and its barcode simulation kit, CellBarcodeSim, that allows efficient and versatile barcode extraction and filtering for a range of barcode types from bulk or single-cell sequencing data using a variety of filtering strategies. Using the barcode simulation kit and biological data, we explore the technical and biological factors influencing barcode identification and provide a decision tree on how to optimize barcode identification for different barcode settings. We believe that CellBarcode and CellBarcodeSim have the capability to enhance the reproducibility and interpretation of barcode results across studies.

deMULTIplex2: robust sample demultiplexing for scRNA-seq

Sample multiplexing enables pooled analysis during single-cell RNA sequencing workflows, thereby increasing throughput and reducing batch effects. A challenge for all multiplexing techniques is to link sample-specific barcodes with cell-specific barcodes, then demultiplex sample identity post-sequencing. However, existing demultiplexing tools fail under many real-world conditions where barcode cross-contamination is an issue. We therefore developed deMULTIplex2, an algorithm inspired by a mechanistic model of barcode cross-contamination. deMULTIplex2 employs generalized linear models and expectation-maximization to probabilistically determine the sample identity of each cell. Benchmarking reveals superior performance across various experimental conditions, particularly on large or noisy datasets with unbalanced sample compositions.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

Lineage Tracing by Single-Cell Transcriptomics Decoding Dynamics of Lineage Commitment

Tracing the fate of individual cells and their progeny is necessary and significant for stem cell research and cancer research. Changes in lineage-specific transcription factor levels during lineage commitment are gradual and continuous. Development of single-cell sequencing technology allows many different states of cells to be sequenced at an unprecedented resolution, and it has been proved that single-cell transcriptomics meets lineage tracing. Here, we introduce a detailed protocol for the lineage tracing by single-cell transcriptomics to clarify the dynamics of lineage commitment.

DNA Barcoding Technology for Lineage Recording and Tracing to Resolve Cell Fate Determination

In various biological contexts, cells receive signals and stimuli that prompt them to change their current state, leading to transitions into a future state. This change underlies the processes of development, tissue maintenance, immune response, and the pathogenesis of various diseases. Following the path of cells from their initial identity to their current state reveals how cells adapt to their surroundings and undergo transformations to attain adjusted cellular states. DNA-based molecular barcoding technology enables the documentation of a phylogenetic tree and the deterministic events of cell lineages, providing the mechanisms and timing of cell lineage commitment that can either promote homeostasis or lead to cellular dysregulation. This review comprehensively presents recently emerging molecular recording technologies that utilize CRISPR/Cas systems, base editing, recombination, and innate variable sequences in the genome. Detailing their underlying principles, applications, and constraints paves the way for the lineage tracing of every cell within complex biological systems, encompassing the hidden steps and intermediate states of organism development and disease progression.

Insights and applications of direct neuronal reprogramming

Direct neuronal reprogramming converts somatic cells of a defined lineage into induced neuronal cells without going through a pluripotent intermediate. This approach not only provides access to the otherwise largely inaccessible cells of the brain for neuronal disease modeling, but also holds great promise for ultimately enabling neuronal cell replacement without the use of transplantation. To improve efficiency and specificity of direct neuronal reprogramming, much of the current efforts aim to understand the mechanisms that safeguard cell identities and how the reprogramming cells overcome the barriers resisting fate changes. Here, we review recent discoveries into the mechanisms by which the donor cell program is silenced, and new cell identities are established. We also discuss advancements that have been made toward fine-tuning the output of these reprogramming systems to generate specific types of neuronal cells. Finally, we highlight the benefit of using direct neuronal reprogramming to study age-related disorders and the potential of in vivo direct reprogramming in regenerative medicine.

A mouse model with high clonal barcode diversity for joint lineage, transcriptomic, and epigenomic profiling in single cells

Cellular lineage histories and their molecular states encode fundamental principles of tissue development and homeostasis. Current lineage-recording mouse models have insufficient barcode diversity and single-cell lineage coverage for profiling tissues composed of millions of cells. Here, we developed DARLIN, an inducible Cas9 barcoding mouse line that utilizes terminal deoxynucleotidyl transferase (TdT) and 30 CRISPR target sites. DARLIN is inducible, generates massive lineage barcodes across tissues, and enables the detection of edited barcodes in ∼70% of profiled single cells. Using DARLIN, we examined fate bias within developing hematopoietic stem cells (HSCs) and revealed unique features of HSC migration. Additionally, we established a protocol for joint transcriptomic and epigenomic single-cell measurements with DARLIN and found that cellular clonal memory is associated with genome-wide DNA methylation rather than gene expression or chromatin accessibility. DARLIN will enable the high-resolution study of lineage relationships and their molecular signatures in diverse tissues and physiological contexts.

Disrupting cellular memory to overcome drug resistance

Gene expression states persist for varying lengths of time at the single-cell level, a phenomenon known as gene expression memory. When cells switch states, losing memory of their prior state, this transition can occur in the absence of genetic changes. However, we lack robust methods to find regulators of memory or track state switching. Here, we develop a lineage tracing-based technique to quantify memory and identify cells that switch states. Applied to melanoma cells without therapy, we quantify long-lived fluctuations in gene expression that are predictive of later resistance to targeted therapy. We also identify the PI3K and TGF-β pathways as state switching modulators. We propose a pretreatment model, first applying a PI3K inhibitor to modulate gene expression states, then applying targeted therapy, which leads to less resistance than targeted therapy alone. Together, we present a method for finding modulators of gene expression memory and their associated cell fates.

Expressed Barcoding Enables High-Resolution Tracking of the Evolution of Drug Tolerance

For a majority of patients with non-small cell lung cancer with EGFR mutations, treatment with EGFR inhibitors (EGFRi) induces a clinical response. Despite this initial reduction in tumor size, residual disease persists that leads to disease relapse. Elucidating the preexisting biological differences between sensitive cells and surviving drug-tolerant persister cells and deciphering how drug-tolerant cells evolve in response to treatment could help identify strategies to improve the efficacy of EGFRi. In this study, we tracked the origins and clonal evolution of drug-tolerant cells at a high resolution by using an expressed barcoding system coupled with single-cell RNA sequencing. This platform enabled longitudinal profiling of gene expression and drug sensitivity in response to EGFRi across a large number of clones. Drug-tolerant cells had higher expression of key survival pathways such as YAP and EMT at baseline and could also differentially adapt their gene expression following EGFRi treatment compared with sensitive cells. In addition, drug combinations targeting common downstream components (MAPK) or orthogonal factors (chemotherapy) showed greater efficacy than EGFRi alone, which is attributable to broader targeting of the heterogeneous EGFRi-tolerance mechanisms present in tumors. Overall, this approach facilitates thorough examination of clonal evolution in response to therapy that could inform the development of improved diagnostic approaches and treatment strategies for targeting drug-tolerant cells.

SIGNIFICANCE

The evolution and heterogeneity of EGFR inhibitor tolerance are identified in a large number of clones at enhanced cellular and temporal resolution using an expressed barcode technology coupled with single-cell RNA sequencing.

Efficient generation of functional neurons from mouse embryonic stem cells via neurogenin-2 expression

The production of induced neuronal (iN) cells from human embryonic stem cells (ESCs) and induced pluripotent stem cells by the forced expression of proneural transcription factors is rapid, efficient and reproducible. The ability to generate large numbers of human neurons in such a robust manner enables large-scale studies of human neural differentiation and neuropsychiatric diseases. Surprisingly, similar transcription factor-based approaches for converting mouse ESCs into iN cells have been challenging, primarily because of low cell survival. Here, we provide a detailed approach for the efficient and reproducible generation of functional iN cells from mouse ESC cultures by the genetically induced expression of neurogenin-2. The resulting iN cells display mature pre- and postsynaptic specializations and form synaptic networks. Our method provides the basis for studying neuronal development and enables the direct comparison of cellular phenotypes in mouse and human neurons generated in an equivalent way. The procedure requires 14 d and can be carried out by users with expertise in stem cell culture.

PC3T: a signature-driven predictor of chemical compounds for cellular transition

Cellular transitions hold great promise in translational medicine research. However, therapeutic applications are limited by the low efficiency and safety concerns of using transcription factors. Small molecules provide a temporal and highly tunable approach to overcome these issues. Here, we present PC3T, a computational framework to enrich molecules that induce desired cellular transitions, and PC3T was able to consistently enrich small molecules that had been experimentally validated in both bulk and single-cell datasets. We then predicted small molecule reprogramming of fibroblasts into hepatic progenitor-like cells (HPLCs). The converted cells exhibited epithelial cell-like morphology and HPLC-like gene expression pattern. Hepatic functions were also observed, such as glycogen storage and lipid accumulation. Finally, we collected and manually curated a cell state transition resource containing 224 time-course gene expression datasets and 153 cell types. Our framework, together with the data resource, is freely available at http://pc3t.idrug.net.cn/ . We believe that PC3T is a powerful tool to promote chemical-induced cell state transitions.

Inferring single-cell transcriptomic dynamics with structured latent gene expression dynamics

Gene expression dynamics provide directional information for trajectory inference from single-cell RNA sequencing data. Traditional approaches compute RNA velocity using strict modeling assumptions about transcription and splicing of RNA. This can fail in scenarios where multiple lineages have distinct gene dynamics or where rates of transcription and splicing are time dependent. We present "LatentVelo," an approach to compute a low-dimensional representation of gene dynamics with deep learning. LatentVelo embeds cells into a latent space with a variational autoencoder and models differentiation dynamics on this "dynamics-based" latent space with neural ordinary differential equations. LatentVelo infers a latent regulatory state that controls the dynamics of an individual cell to model multiple lineages. LatentVelo can predict latent trajectories, describing the inferred developmental path for individual cells rather than just local RNA velocity vectors. The dynamics-based embedding batch corrects cell states and velocities, outperforming comparable autoencoder batch correction methods that do not consider gene expression dynamics.

Single-cell epigenome analysis identifies molecular events controlling direct conversion of human fibroblasts to pancreatic ductal-like cells

Cell fate can be reprogrammed by ectopic expression of lineage-specific transcription factors (TFs). However, the exact cell state transitions during transdifferentiation are still poorly understood. Here, we have generated pancreatic exocrine cells of ductal epithelial identity from human fibroblasts using a set of six TFs. We mapped the molecular determinants of lineage dynamics using a factor-indexing method based on single-nuclei multiome sequencing (FI-snMultiome-seq) that enables dissecting the role of each individual TF and pool of TFs in cell fate conversion. We show that transition from mesenchymal fibroblast identity to epithelial pancreatic exocrine fate involves two deterministic steps: an endodermal progenitor state defined by activation of HHEX with FOXA2 and SOX17 and a temporal GATA4 activation essential for the maintenance of pancreatic cell fate program. Collectively, our data suggest that transdifferentiation-although being considered a direct cell fate conversion method-occurs through transient progenitor states orchestrated by stepwise activation of distinct TFs.

Multiplexed transcriptomic profiling of the fate of human CAR T cells in vivo via genetic barcoding with shielded small nucleotides

The design of chimeric antigen receptor (CAR) T cells would benefit from knowledge of the fate of the cells in vivo. This requires the permanent labelling of CAR T cell products and their pooling in the same microenvironment. Here, we report a cell-barcoding method for the multiplexed longitudinal profiling of cells in vivo using single-cell RNA sequencing (scRNA-seq). The method, which we named shielded-small-nucleotide-based scRNA-seq (SSN-seq), is compatible with both 3' and 5' single-cell profiling, and enables the recording of cell identity, from cell infusion to isolation, by leveraging the ubiquitous Pol III U6 promoters to robustly express small-RNA barcodes modified with direct-capture sequences. By using SSN-seq to track the dynamics of the states of CAR T cells in a tumour-rechallenge mouse model of leukaemia, we found that a combination of cytokines and small-molecule inhibitors that are used in the ex vivo manufacturing of CAR T cells promotes the in vivo expansion of persistent populations of CD4+ memory T cells. By facilitating the probing of cell-state dynamics in vivo, SSN-seq may aid the development of adoptive cell therapies.

FOXC2 promotes vasculogenic mimicry and resistance to anti-angiogenic therapy

Vasculogenic mimicry (VM) describes the formation of pseudo blood vessels constructed of tumor cells that have acquired endothelial-like properties. VM channels endow the tumor with a tumor-derived vascular system that directly connects to host blood vessels, and their presence is generally associated with poor patient prognosis. Here we show that the transcription factor, Foxc2, promotes VM in diverse solid tumor types by driving ectopic expression of endothelial genes in tumor cells, a process that is stimulated by hypoxia. VM-proficient tumors are resistant to anti-angiogenic therapy, and suppression of Foxc2 augments response. This work establishes co-option of an embryonic endothelial transcription factor by tumor cells as a key mechanism driving VM proclivity and motivates the search for VM-inhibitory agents that could form the basis of combination therapies with anti-angiogenics.

Cell-type specification in the retina: Recent discoveries from transcriptomic approaches

Understanding the formation of the complex nervous system hinges on decoding the mechanism that specifies a vast array of neuronal types, each endowed with a unique morphology, physiology, and connectivity. As a pivotal step towards addressing this problem, seminal work has been devoted to characterizing distinct neuronal types. In recent years, high-throughput, single-cell transcriptomic methods have enabled a rapid inventory of cell types in various regions of the nervous system, with the retina exhibiting complete molecular characterization across many vertebrate species. This invaluable resource has furnished a fresh perspective for investigating the molecular principles of cell-type specification, thereby advancing our understanding of retinal development. Accordingly, this review focuses on the most recent transcriptomic characterizations of retinal cells, with a particular focus on amacrine cells and retinal ganglion cells. These investigations have unearthed new insights into their cell-type specification.

Comprehensive spatiotemporal mapping of single-cell lineages in developing mouse brain by CRISPR-based barcoding

A fundamental interest in developmental neuroscience lies in the ability to map the complete single-cell lineages within the brain. To this end, we developed a CRISPR editing-based lineage-specific tracing (CREST) method for clonal tracing in Cre mice. We then used two complementary strategies based on CREST to map single-cell lineages in developing mouse ventral midbrain (vMB). By applying snapshotting CREST (snapCREST), we constructed a spatiotemporal lineage landscape of developing vMB and identified six progenitor archetypes that could represent the principal clonal fates of individual vMB progenitors and three distinct clonal lineages in the floor plate that specified glutamatergic, dopaminergic or both neurons. We further created pandaCREST (progenitor and derivative associating CREST) to associate the transcriptomes of progenitor cells in vivo with their differentiation potentials. We identified multiple origins of dopaminergic neurons and demonstrated that a transcriptome-defined progenitor type comprises heterogeneous progenitors, each with distinct clonal fates and molecular signatures. Therefore, the CREST method and strategies allow comprehensive single-cell lineage analysis that could offer new insights into the molecular programs underlying neural specification.

Diverse clonal fates emerge upon drug treatment of homogeneous cancer cells

Even among genetically identical cancer cells, resistance to therapy frequently emerges from a small subset of those cells1-7. Molecular differences in rare individual cells in the initial population enable certain cells to become resistant to therapy7-9; however, comparatively little is known about the variability in the resistance outcomes. Here we develop and apply FateMap, a framework that combines DNA barcoding with single-cell RNA sequencing, to reveal the fates of hundreds of thousands of clones exposed to anti-cancer therapies. We show that resistant clones emerging from single-cell-derived cancer cells adopt molecularly, morphologically and functionally distinct resistant types. These resistant types are largely predetermined by molecular differences between cells before drug addition and not by extrinsic factors. Changes in the dose and type of drug can switch the resistant type of an initial cell, resulting in the generation and elimination of certain resistant types. Samples from patients show evidence for the existence of these resistant types in a clinical context. We observed diversity in resistant types across several single-cell-derived cancer cell lines and cell types treated with a variety of drugs. The diversity of resistant types as a result of the variability in intrinsic cell states may be a generic feature of responses to external cues.

NF-kB's contribution to B cell fate decisions

NF-κB signaling is essential to an effective innate and adaptive immune response. Many immune-specific functional and developmental outcomes depend in large on NF-κB. The formidable task of sorting out the mechanisms behind the regulation and outcome of NF-κB signaling remains an important area of immunology research. Here we briefly discuss the role of NF-κB in regulating cell fate decisions at various times in the path of B cell development, activation, and the generation of long-term humoral immunity.