Single-cell lineages reveal the rates, routes, and drivers of metastasis in cancer xenografts

Detecting branching rate heterogeneity in multifurcating trees with applications in lineage tracing data

Understanding cellular birth rate differences is crucial for predicting cancer progression and interpreting tumor-derived genetic data. Lineage tracing experiments enable detailed reconstruction of cellular genealogies, offering new opportunities to measure branching rate heterogeneity. However, the lineage tracing process can introduce complex tree features that complicate this effort. Here, we examine tree characteristics in lineage tracing-derived genealogies and find that editing window placement leads to multifurcations at a tree's root or tips. We propose several ways in which existing tree topology-based metrics can be extended to test for rate heterogeneity on trees even in the presence of lineage-tracing associated distortions. Although these methods vary in power and robustness, a test based on the J 1 statistic effectively detects branching rate heterogeneity in simulated lineage tracing data. Tests based on other common statistics ( ŝ and the Sackin index) show interior performance to J 1 . We apply our validated methods to xenograft experimental data and find widespread rate heterogeneity across multiple study systems. Our results demonstrate the potential of tree topology statistics in analyzing lineage tracing data, and highlight the challenges associated with adapting phylogenetic methods to these systems.

Reconstructing axial progenitor field dynamics in mouse stem cell-derived embryoids

Embryogenesis requires substantial coordination to translate genetic programs to the collective behavior of differentiating cells, but understanding how cellular decisions control tissue morphology remains conceptually and technically challenging. Here, we combine continuous Cas9-based molecular recording with a mouse embryonic stem cell-based model of the embryonic trunk to build single-cell phylogenies that describe the behavior of transient, multipotent neuro-mesodermal progenitors (NMPs) as they commit into neural and somitic cell types. We find that NMPs show subtle transcriptional signatures related to their recent differentiation and contribute to downstream lineages through a surprisingly broad distribution of individual fate outcomes. Although decision-making can be heavily influenced by environmental cues to induce morphological phenotypes, axial progenitors intrinsically mature over developmental time to favor the neural lineage. Using these data, we present an experimental and analytical framework for exploring the non-homeostatic dynamics of transient progenitor populations as they shape complex tissues during critical developmental windows.

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.

Heterogeneity generating capacity in tumorigenesis and cancer therapeutics

Cells of multicellular organisms generate heterogeneity in a controlled and transient fashion during embryogenesis, which can be reactivated in pathologies such as cancer. Although genomic heterogeneity is an important part of tumorigenesis, continuous generation of phenotypic heterogeneity is central for the adaptation of cancer cells to the challenges of tumorigenesis and response to therapy. Here I discuss the capacity of generating heterogeneity, hereafter called cell hetness, in cancer cells both as the activation of hetness oncogenes and inactivation of hetness tumor suppressor genes, which increase the generation of heterogeneity, ultimately producing an increase in adaptability and cell fitness. Transcriptomic high hetness states in therapy-tolerant cell states denote its importance in cancer resistance to therapy. The definition of the concept of hetness will allow the understanding of its origins, its control during embryogenesis, its loss of control in tumorigenesis and cancer therapeutics and its active targeting.

Assessment of the impact of residual tumors at different sites post-neoadjuvant chemotherapy on prognosis in breast cancer patients and development of a disease-free survival prediction model

Background

Residual disease after neoadjuvant chemotherapy (NAC) in breast cancer patients predicts worse outcomes than pathological complete response. Differing prognostic impacts based on the anatomical site of residual tumors are not well studied.

Objectives

The study aims to assess disease-free survival (DFS) in breast cancer patients with different residual tumor sites following NAC and to develop a nomogram for predicting 1- to 3-year DFS in these patients.

Design

A retrospective cohort study.

Methods

Retrospective analysis of 953 lymph node-positive breast cancer patients with residual disease post-NAC. Patients were categorized into three groups: residual disease in breast (RDB), residual disease in lymph nodes (RDN), and residual disease in both (RDBN). DFS compared among groups. Patients were divided into a training set and a validation set in a 7:3 ratio. Prognostic factors for DFS were analyzed to develop a nomogram prediction model.

Results

RDB patients had superior 3-year DFS of 94.6% versus 85.2% for RDN and 81.8% for RDBN (p < 0.0001). Clinical T stage, N stage, molecular subtype, and postoperative pN stage were independently associated with DFS on both univariate and multivariate analyses. Nomogram integrating clinical tumor-node-metastasis (TNM) stage, molecular subtype, pathological response demonstrated good discrimination (C-index 0.748 training, 0.796 validation cohort), and calibration.

Conclusion

The location of residual disease has prognostic implications, with nodal residuals predicting poorer DFS. The validated nomogram enables personalized DFS prediction to guide treatment decisions.

How CRISPR Is Revolutionizing the Generation of New Models for Cancer Research

Cancers arise through acquisition of mutations in genes that regulate core biological processes like cell proliferation and cell death. Decades of cancer research have led to the identification of genes and mutations causally involved in disease development and evolution, yet defining their precise function across different cancer types and how they influence therapy responses has been challenging. Mouse models have helped define the in vivo function of cancer-associated alterations, and genome-editing approaches using CRISPR have dramatically accelerated the pace at which these models are developed and studied. Here, we highlight how CRISPR technologies have impacted the development and use of mouse models for cancer research and discuss the many ways in which these rapidly evolving platforms will continue to transform our understanding of this disease.

Lung endothelium exploits susceptible tumor cell states to instruct metastatic latency

In metastasis, cancer cells travel around the circulation to colonize distant sites. Due to the rarity of these events, the immediate fates of metastasizing tumor cells (mTCs) are poorly understood while the role of the endothelium as a dissemination interface remains elusive. Using a newly developed combinatorial mTC enrichment approach, we provide a transcriptional blueprint of the early colonization process. Following their arrest at the metastatic site, mTCs were found to either proliferate intravascularly or extravasate, thereby establishing metastatic latency. Endothelial-derived angiocrine Wnt factors drive this bifurcation, instructing mTCs to follow the extravasation-latency route. Surprisingly, mTC responsiveness towards niche-derived Wnt was established at the epigenetic level, which predetermined tumor cell behavior. Whereas hypomethylation enabled high Wnt activity leading to metastatic latency, methylated mTCs exhibited low activity and proliferated intravascularly. Collectively the data identify the predetermined methylation status of disseminated tumor cells as a key regulator of mTC behavior in the metastatic niche.

Decoding the interplay between genetic and non-genetic drivers of metastasis

Metastasis is a multistep process by which cancer cells break away from their original location and spread to distant organs, and is responsible for the vast majority of cancer-related deaths. Preventing early metastatic dissemination would revolutionize the ability to fight cancer. Unfortunately, the relatively poor understanding of the molecular underpinnings of metastasis has hampered the development of effective anti-metastatic drugs. Although it is now accepted that disseminating tumour cells need to acquire multiple competencies to face the many obstacles they encounter before reaching their metastatic site(s), whether these competencies are acquired through an accumulation of metastasis-specific genetic alterations and/or non-genetic events is often debated. Here we review a growing body of literature highlighting the importance of both genetic and non-genetic reprogramming events during the metastatic cascade, and discuss how genetic and non-genetic processes act in concert to confer metastatic competencies. We also describe how recent technological advances, and in particular the advent of single-cell multi-omics and barcoding approaches, will help to better elucidate the cross-talk between genetic and non-genetic mechanisms of metastasis and ultimately inform innovative paths for the early detection and interception of this lethal process.

Genetic Tools for Cell Lineage Tracing and Profiling Developmental Trajectories in the Skin

The epidermis is the body's first line of protection against dehydration and pathogens, continually regenerating the outermost protective skin layers throughout life. During both embryonic development and wound healing, epidermal stem and progenitor cells must respond to external stimuli and insults to build, maintain, and repair the cutaneous barrier. Recent advances in CRISPR-based methods for cell lineage tracing have remarkably expanded the potential for experiments that track stem and progenitor cell proliferation and differentiation over the course of tissue and even organismal development. Additional tools for DNA-based recording of cellular signaling cues promise to deepen our understanding of the mechanisms driving normal skin morphogenesis and response to stressors as well as the dysregulation of cell proliferation and differentiation in skin diseases and cancer. In this review, we highlight cutting-edge methods for cell lineage tracing, including in organoids and model organisms, and explore how cutaneous biology researchers might leverage these techniques to elucidate the developmental programs that support the regenerative capacity and plasticity of the skin.

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.

NestedBD: Bayesian inference of phylogenetic trees from single-cell copy number profiles under a birth-death model

Copy number aberrations (CNAs) are ubiquitous in many types of cancer. Inferring CNAs from cancer genomic data could help shed light on the initiation, progression, and potential treatment of cancer. While such data have traditionally been available via "bulk sequencing," the more recently introduced techniques for single-cell DNA sequencing (scDNAseq) provide the type of data that makes CNA inference possible at the single-cell resolution. We introduce a new birth-death evolutionary model of CNAs and a Bayesian method, NestedBD, for the inference of evolutionary trees (topologies and branch lengths with relative mutation rates) from single-cell data. We evaluated NestedBD's performance using simulated data sets, benchmarking its accuracy against traditional phylogenetic tools as well as state-of-the-art methods. The results show that NestedBD infers more accurate topologies and branch lengths, and that the birth-death model can improve the accuracy of copy number estimation. And when applied to biological data sets, NestedBD infers plausible evolutionary histories of two colorectal cancer samples. NestedBD is available at https://github.com/Androstane/NestedBD .

Single-Cell Profiling Reveals Immune-Based Mechanisms Underlying Tumor Radiosensitization by a Novel Mn Porphyrin Clinical Candidate, MnTnBuOE-2-PyP5+ (BMX-001)

Manganese porphyrins reportedly exhibit synergic effects when combined with irradiation. However, an in-depth understanding of intratumoral heterogeneity and immune pathways, as affected by Mn porphyrins, remains limited. Here, we explored the mechanisms underlying immunomodulation of a clinical candidate, MnTnBuOE-2-PyP5+ (BMX-001, MnBuOE), using single-cell analysis in a murine carcinoma model. Mice bearing 4T1 tumors were divided into four groups: control, MnBuOE, radiotherapy (RT), and combined MnBuOE and radiotherapy (MnBuOE/RT). In epithelial cells, the epithelial-mesenchymal transition, TNF-α signaling via NF-кB, angiogenesis, and hypoxia-related genes were significantly downregulated in the MnBuOE/RT group compared with the RT group. All subtypes of cancer-associated fibroblasts (CAFs) were clearly reduced in MnBuOE and MnBuOE/RT. Inhibitory receptor-ligand interactions, in which epithelial cells and CAFs interacted with CD8+ T cells, were significantly lower in the MnBuOE/RT group than in the RT group. Trajectory analysis showed that dendritic cells maturation-associated markers were increased in MnBuOE/RT. M1 macrophages were significantly increased in the MnBuOE/RT group compared with the RT group, whereas myeloid-derived suppressor cells were decreased. CellChat analysis showed that the number of cell-cell communications was the lowest in the MnBuOE/RT group. Our study is the first to provide evidence for the combined radiotherapy with a novel Mn porphyrin clinical candidate, BMX-001, from the perspective of each cell type within the tumor microenvironment.

Punctuational evolution is pervasive in distal site metastatic colonization

The evolution of metastasis represents a lethal stage of cancer progression. Yet, the evolutionary kinetics of metastatic disease remain unresolved. Here, using single cell CRISPR-Cas9 lineage tracing data, we show that in metastatic disease, gradual molecular evolution is punctuated by episodes of rapid evolutionary change associated with lineage divergence. By measuring punctuational effects across the metastatic cascade, we show that punctuational effects contribute more to the molecular diversity at distal site metastases compared to the paired primary tumor, suggesting qualitatively different modes of evolution may drive primary and metastatic tumor progression. This is the first empirical evidence for distinct patterns of molecular evolution at early and late stages of metastasis and demonstrates the complex interplay of cell intrinsic and extrinsic factors that shape lethal cancer.

Developmental progression of DNA double-strand break repair deciphered by a single-allele resolution mutation classifier

DNA double-strand breaks (DSBs) are repaired by a hierarchically regulated network of pathways. Factors influencing the choice of particular repair pathways, however remain poorly characterized. Here we develop an Integrated Classification Pipeline (ICP) to decompose and categorize CRISPR/Cas9 generated mutations on genomic target sites in complex multicellular insects. The ICP outputs graphic rank ordered classifications of mutant alleles to visualize discriminating DSB repair fingerprints generated from different target sites and alternative inheritance patterns of CRISPR components. We uncover highly reproducible lineage-specific mutation fingerprints in individual organisms and a developmental progression wherein Microhomology-Mediated End-Joining (MMEJ) or Insertion events predominate during early rapid mitotic cell cycles, switching to distinct subsets of Non-Homologous End-Joining (NHEJ) alleles, and then to Homology-Directed Repair (HDR)-based gene conversion. These repair signatures enable marker-free tracking of specific mutations in dynamic populations, including NHEJ and HDR events within the same samples, for in-depth analysis of diverse gene editing events.

Joint analysis of mutational and transcriptional landscapes in human cancer reveals key perturbations during cancer evolution

BACKGROUND

Tumors are able to acquire new capabilities, including traits such as drug resistance and metastasis that are associated with unfavorable clinical outcomes. Single-cell technologies have made it possible to study both mutational and transcriptomic profiles, but as most studies have been conducted on model systems, little is known about cancer evolution in human patients. Hence, a better understanding of cancer evolution could have important implications for treatment strategies.

RESULTS

Here, we analyze cancer evolution and clonal selection by jointly considering mutational and transcriptomic profiles of single cells acquired from tumor biopsies from 49 lung cancer samples and 51 samples with chronic myeloid leukemia. Comparing the two profiles, we find that each clone is associated with a preferred transcriptional state. For metastasis and drug resistance, we find that the number of mutations affecting related genes increases as the clone evolves, while changes in gene expression profiles are limited. Surprisingly, we find that mutations affecting ligand-receptor interactions with the tumor microenvironment frequently emerge as clones acquire drug resistance.

CONCLUSIONS

Our results show that lung cancer and chronic myeloid leukemia maintain a high clonal and transcriptional diversity, and we find little evidence in favor of clonal sweeps. This suggests that for these cancers selection based solely on growth rate is unlikely to be the dominating driving force during cancer evolution.

Mutations, Bottlenecks, and Clonal Sweeps: How Environmental Carcinogens and Genomic Changes Shape Clonal Evolution during Tumor Progression

The transition from a single, initiated cell to a full-blown malignant tumor involves significant genomic evolution. Exposure to carcinogens-whether directly mutagenic or not-can drive progression toward malignancy, as can stochastic acquisition of cancer-promoting genetic events. Mouse models using both carcinogens and germline genetic manipulations have enabled precise inquiry into the evolutionary dynamics that take place as a tumor progresses from benign to malignant to metastatic stages. Tumor progression is characterized by changes in somatic point mutations and copy-number alterations, even though any single tumor can itself have a high or low burden of genomic alterations. Further, lineage-tracing, single-cell analyses and CRISPR barcoding have revealed the distinct clonal dynamics within benign and malignant tumors. Application of these tools in a range of mouse models can shed unique light on the patterns of clonal evolution that take place in both mouse and human tumors.

Past, present, and future of CRISPR genome editing technologies

Genome editing has been a transformative force in the life sciences and human medicine, offering unprecedented opportunities to dissect complex biological processes and treat the underlying causes of many genetic diseases. CRISPR-based technologies, with their remarkable efficiency and easy programmability, stand at the forefront of this revolution. In this Review, we discuss the current state of CRISPR gene editing technologies in both research and therapy, highlighting limitations that constrain them and the technological innovations that have been developed in recent years to address them. Additionally, we examine and summarize the current landscape of gene editing applications in the context of human health and therapeutics. Finally, we outline potential future developments that could shape gene editing technologies and their applications in the coming years.

In vivo CRISPR knockout screen identifies p47 as a suppressor of HER2+ breast cancer metastasis by regulating NEMO trafficking and autophagy flux

Autophagy is a conserved cellular process, and its dysfunction is implicated in cancer and other diseases. Here, we employ an in vivo CRISPR screen targeting genes implicated in the regulation of autophagy to identify the Nsfl1c gene encoding p47 as a suppressor of human epidermal growth factor receptor 2 (HER2)+ breast cancer metastasis. p47 ablation specifically increases metastasis without promoting primary mammary tumor growth. Analysis of human breast cancer patient databases and tissue samples indicates a correlation of lower p47 expression levels with metastasis and decreased survival. Mechanistic studies show that p47 functions in the repair of lysosomal damage for autophagy flux and in the endosomal trafficking of nuclear factor κB essential modulator for lysosomal degradation to promote metastasis. Our results demonstrate a role and mechanisms of p47 in the regulation of breast cancer metastasis. They highlight the potential to exploit p47 as a suppressor of metastasis through multiple pathways in HER2+ breast cancer cells.

Exact and efficient phylodynamic simulation from arbitrarily large populations

Many biological studies involve inferring the genealogical history of a sample of individuals from a large population and interpreting the reconstructed tree. Such an ascertained tree typically represents only a small part of a comprehensive population tree and is distorted by survivorship and sampling biases. Inferring evolutionary parameters from ascertained trees requires modeling both the underlying population dynamics and the ascertainment process. A crucial component of this phylodynamic modeling involves tree simulation, which is used to benchmark probabilistic inference methods. To simulate an ascertained tree, one must first simulate the full population tree and then prune unobserved lineages. Consequently, the computational cost is determined not by the size of the final simulated tree, but by the size of the population tree in which it is embedded. In most biological scenarios, simulations of the entire population are prohibitively expensive due to computational demands placed on lineages without sampled descendants. Here, we address this challenge by proving that, for any partially ascertained process from a general multi-type birth-death-mutation-sampling (BDMS) model, there exists an equivalent pure birth process (i.e., no death) with mutation and complete sampling. The final trees generated under these processes have exactly the same distribution. Leveraging this property, we propose a highly efficient algorithm for simulating trees under a general BDMS model. Our algorithm scales linearly with the size of the final simulated tree and is independent of the population size, enabling simulations from extremely large populations beyond the reach of current methods but essential for various biological applications. We anticipate that this unprecedented speedup will significantly advance the development of novel inference methods that require extensive training data.

A novel batch-effect correction method for scRNA-seq data based on Adversarial Information Factorization

Single-cell RNA sequencing (scRNA-seq) technology produces an unprecedented resolution at the level of a unique cell, raising great hopes in medicine. Nevertheless, scRNA-seq data suffer from high variations due to the experimental conditions, called batch effects, preventing any aggregated downstream analysis. Adversarial Information Factorization provides a robust batch-effect correction method that does not rely on prior knowledge of the cell types nor a specific normalization strategy while being adapted to any downstream analysis task. It compares to and even outperforms state-of-the-art methods in several scenarios: low signal-to-noise ratio, batch-specific cell types with few cells, and a multi-batches dataset with imbalanced batches and batch-specific cell types. Moreover, it best preserves the relative gene expression between cell types, yielding superior differential expression analysis results. Finally, in a more complex setting of a Leukemia cohort, our method preserved most of the underlying biological information for each patient while aligning the batches, improving the clustering metrics in the aggregated dataset.

Single-cell lineage tracing with endogenous markers

Resolving lineage relationships between cells in an organism provides key insights into the fate of individual cells and drives a fundamental understanding of the process of development and disease. A recent rapid increase in experimental and computational advances for detecting naturally occurring somatic nuclear and mitochondrial mutation at single-cell resolution has expanded lineage tracing from model organisms to humans. This review discusses the advantages and challenges of experimental and computational techniques for cell lineage tracing using somatic mutation as endogenous DNA barcodes to decipher the relationships between cells during development and tumour evolution. We outlook the advantages of spatial clonal evolution analysis and single-cell lineage tracing using endogenous genetic markers.

Data Storage Using DNA

The exponential growth of global data has outpaced the storage capacities of current technologies, necessitating innovative storage strategies. DNA, as a natural medium for preserving genetic information, has emerged as a highly promising candidate for next-generation storage medium. Storing data in DNA offers several advantages, including ultrahigh physical density and exceptional durability. Facilitated by significant advancements in various technologies, such as DNA synthesis, DNA sequencing, and DNA nanotechnology, remarkable progress has been made in the field of DNA data storage over the past decade. However, several challenges still need to be addressed to realize practical applications of DNA data storage. In this review, the processes and strategies of in vitro DNA data storage are first introduced, highlighting recent advancements. Next, a brief overview of in vivo DNA data storage is provided, with a focus on the various writing strategies developed to date. At last, the challenges encountered in each step of DNA data storage are summarized and promising techniques are discussed that hold great promise in overcoming these obstacles.

Lineage Plasticity: The New Cancer Hallmark on the Block

Plasticity refers to the ability of cells to adopt a spectrum of states or phenotypes. In cancer, it is a critical contributor to tumor initiation, progression, invasiveness, and therapy resistance, and it has recently been recognized as an emerging cancer hallmark. Plasticity can occur as a result of cell-intrinsic factors (e.g., genetic, transcriptional, or epigenetic fluctuations), or through cell-extrinsic cues (e.g., signaling from components of the tumor microenvironment or selective pressure from therapy). Over the past decade, technological advances, analysis of patient samples, and studies in mouse model systems have led to a deeper understanding of how such plastic states come about. In this review, we discuss: (i) the definition of plasticity; (ii) methods to measure and quantify plasticity; (iii) the clinical relevance of plasticity; and (iv) therapeutic hypotheses to modulate plasticity in the clinic.

Charting Single Cell Lineage Dynamics and Mutation Networks via Homing CRISPR

Single cell lineage tracing, essential for unraveling cellular dynamics in disease evolution is critical for developing targeted therapies. CRISPR-Cas9, known for inducing permanent and cumulative mutations, is a cornerstone in lineage tracing. The novel homing guide RNA (hgRNA) technology enhances this by enabling dynamic retargeting and facilitating ongoing genetic modifications. Charting these mutations, especially through successive hgRNA edits, poses a significant challenge. Our solution, LINEMAP, is a computational framework designed to trace and map these mutations with precision. LINEMAP meticulously discerns mutation alleles at single-cell resolution and maps their complex interrelationships through a mutation evolution network. By utilizing a Markov Process model, we can predict mutation transition probabilities, revealing potential mutational routes and pathways. Our reconstruction algorithm, anchored in the Markov model's attributes, reconstructs cellular lineage pathways, shedding light on the cell's evolutionary journey to the minutiae of single-cell division. Our findings reveal an intricate network of mutation evolution paired with a predictive Markov model, advancing our capability to reconstruct single-cell lineage via hgRNA. This has substantial implications for advancing our understanding of biological mechanisms and propelling medical research forward.

KRT17 From Keratinocytes With High Glucose Stimulation Inhibit Dermal Fibroblasts Migration Through Integrin α11

Objective

To investigate the effects of overexpressed keratin 17 (KRT17) on the biology of human dermal fibroblasts (HDFs) and to explore the mechanism of KRT17 in diabetic wound healing.

Methods

KRT17 expression was tested in diabetic keratinocytes, animal models, and patient skin tissues (Huazhong University of Science and Technology Ethics Committee, [2022] No. 3110). Subsequently, HDFs were stimulated with different concentrations of KRT17 in vitro. Changes in the proliferation and migration of HDFs were observed. Then, identification of KRT17-induced changes in dermal fibroblast of RNA sequencing-based transcriptome analysis was performed.

Results

KRT17 expression was upregulated under pathological conditions. In vitro stimulation of HDFs with different concentrations of KRT17 inhibited cell migration. RNA-seq data showed that enriched GO terms were extracellular matrix components and their regulation. KEGG analysis revealed that the highest number of enriched genes was PI3K-Akt, in which integrin alpha-11 (ITGA11) mRNA, a key molecule that regulates cell migration, was significantly downregulated. Decreased ITGA11 expression was observed after stimulation of HDFs with KRT17 in vitro.

Conclusion

Increased expression of KRT17 in diabetic pathological surroundings inhibits fibroblast migration by downregulating the expression of ITGA11. Thus, KRT17 may be a molecular target for the treatment of diabetic wounds.

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.

Alone you go faster, together you go farther

The metastatic process is an extraordinarily complex step-by-step procedure, characterized by many analogies with migratory patterns of humans or animals across our planet. The ongoing interrogation of circulating tumor cells (CTCs), caught in the act of spreading from one location to another, is revealing distinct behaviors including biological, physical, and mechanical features that impact on their likelihood to form metastasis. In this viewpoint, I will discuss some of these findings and provide a perspective on the metastatic journey, open questions and opportunities to exploit some of the most recent discoveries for the development of antimetastasis medicines.

ConvexML: Scalable and accurate inference of single-cell chronograms from CRISPR/Cas9 lineage tracing data

CRISPR/Cas9 gene editing technology has enabled lineage tracing for thousands of cells in vivo. However, most of the analysis of CRISPR/Cas9 lineage tracing data has so far been limited to the reconstruction of single-cell tree topologies, which depict lineage relationships between cells, but not the amount of time that has passed between ancestral cell states and the present. Time-resolved trees, known as chronograms, would allow one to study the evolutionary dynamics of cell populations at an unprecedented level of resolution. Indeed, time-resolved trees would reveal the timing of events on the tree, the relative fitness of subclones, and the dynamics underlying phenotypic changes in the cell population - among other important applications. In this work, we introduce the first scalable and accurate method to refine any given single-cell tree topology into a single-cell chronogram by estimating its branch lengths. To do this, we leverage a statistical model of CRISPR/Cas9 cutting with missing data, paired with a conservative version of maximum parsimony that reconstructs only the ancestral states that we are confident about. As part of our method, we propose a novel approach to represent and handle missing data - specifically, double-resection events - which greatly simplifies and speeds up branch length estimation without compromising quality. All this leads to a convex maximum likelihood estimation (MLE) problem that can be readily solved in seconds with off-the-shelf convex optimization solvers. To stabilize estimates in low-information regimes, we propose a simple penalized version of MLE using a minimum branch length and pseudocounts. We benchmark our method using simulations and show that it performs well on several tasks, outperforming more naive baselines. Our method, which we name 'ConvexML', is available through the cassiopeia open source Python package.

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.

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.

Application of single-cell sequencing to the research of tumor microenvironment

Single-cell sequencing is a technique for detecting and analyzing genomes, transcriptomes, and epigenomes at the single-cell level, which can detect cellular heterogeneity lost in conventional sequencing hybrid samples, and it has revolutionized our understanding of the genetic heterogeneity and complexity of tumor progression. Moreover, the tumor microenvironment (TME) plays a crucial role in the formation, development and response to treatment of tumors. The application of single-cell sequencing has ushered in a new age for the TME analysis, revealing not only the blueprint of the pan-cancer immune microenvironment, but also the heterogeneity and differentiation routes of immune cells, as well as predicting tumor prognosis. Thus, the combination of single-cell sequencing and the TME analysis provides a unique opportunity to unravel the molecular mechanisms underlying tumor development and progression. In this review, we summarize the recent advances in single-cell sequencing and the TME analysis, highlighting their potential applications in cancer research and clinical translation.

RAI2 acts as a tumor suppressor with functional significance in gastric cancer

Metastasis of gastric cancer (GC) is one of the major causes of death among GC patients. GC metastasis involves numerous biological processes, yet the specific molecular biological mechanisms have not been elucidated. Here, we report a novel tumor suppressor, retinoic acid-induced 2 (RAI2), which is located in the Xp22 region of the chromosome and plays a role in inhibiting GC growth and invasion. In this study, integrated analysis of The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) datasets and immunohistochemistry staining data suggested that RAI2 expression in GC samples was low. Moreover, the immune infiltration analysis indicated that low expression of RAI2 in GC was associated with a higher intensity of tumor-infiltrating lymphocytes (TILs) and an abundance of Programmed death ligand 1 (PD-L1) expression. Gene set enrichment analysis (GSEA) analysis further revealed that RAI2 regulated some pathways including the GAP junction, focal adhesion and ECM receptor interaction pathway, immune regulation, PI3K-Akt signaling, MAPK signaling, cell cycle, and DNA replication. Furthermore, the knockdown of RAI2 promoted GC cell proliferation, migration, and invasion in vitro. Taken together, these results suggest that the tumor suppressor RAI2 could be a potential target for the development of anti-cancer strategies in GC.

KRT17 from skin cells with high glucose stimulation promotes keratinocytes proliferation and migration

Impaired diabetic wound healing is an important issue in diabetic complications. Proliferation and migration of keratinocytes are major processes of skin wound repair after injury. However, hyperkeratosis can affect the speed of wound healing. Based on the results of preliminary experiments on increased KRT17 expression after high glucose stimulation of human skin tissue cells, a cell model of human immortalized keratinocyte (HaCaT) stimulation with different concentrations of KRT17 was established in vitro, and the promotion in cell proliferation and migration were discovered. KRT17 silencing promoted diabetic wound healing in the db/db diabetic wound model. Transcriptome sequencing (RNA-seq) was performed on HaCaT cells after KRT17 stimulation, and analysis showed significant enrichment in the PI3K-AKT signaling pathway, in which the regulation of cell c-MYB mRNA, a key molecule regulating cell proliferation and migration, was significantly upregulated. In vitro assays showed increased c-MYB expression and enhanced pAKT activity after HaCaT cell stimulation by KRT17. We speculate that KRT17 is upregulated under high glucose and promotes keratinocyte proliferation and migration caused hyperkeratosis, through the c-MYB/PI3K-AKT pathway, contributing to delayed wound healing.

Single-cell RNA sequencing reveals the cellular and molecular characteristics of high-grade and metastatic bladder cancer

PURPOSE

Metastatic bladder cancer (BC) has the highest somatic mutation frequency and recurrence rate of all tumors. However, the cellular and molecular characteristics of BC remain unclear.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) on the samples of paracancerous normal tissue (PNT), primary tumor (PT) and lymph node metastasis (LNM). The proportions and gene expression profiles of different cell types in the tumor microenvironment (TME) were investigated.

RESULTS

In total, 50,158 cells were classified into six populations. Malignant cells of PT and LNM exhibited large mutant DNA fragments, while the cell phenotypes and gene expression profiles differed during differentiation. Metastasis was associated with a poorer prognosis than PT. Tumor-associated stromal cells and inhibitory immune cells were the main cell populations in PT and LNM. Cell-cell communication analysis revealed the roles of signaling pathways of inflammatory cancer-associated fibroblast (iCAF) and tumor-associated macrophage (TAM) in exhaustion of T cells. In addition, iCAF may recruit TAM to promote formation of the TME earlier than the differentiation of tumor cells.

CONCLUSION

This study through scRNA-seq enhanced our understanding of new features about the cellular and molecular similarities and differences of high-grade and metastatic bladder cancer, which might provide potential therapeutic targets in future treatment.

Lymph node metastasis in cancer progression: molecular mechanisms, clinical significance and therapeutic interventions

Lymph nodes (LNs) are important hubs for metastatic cell arrest and growth, immune modulation, and secondary dissemination to distant sites through a series of mechanisms, and it has been proved that lymph node metastasis (LNM) is an essential prognostic indicator in many different types of cancer. Therefore, it is important for oncologists to understand the mechanisms of tumor cells to metastasize to LNs, as well as how LNM affects the prognosis and therapy of patients with cancer in order to provide patients with accurate disease assessment and effective treatment strategies. In recent years, with the updates in both basic and clinical studies on LNM and the application of advanced medical technologies, much progress has been made in the understanding of the mechanisms of LNM and the strategies for diagnosis and treatment of LNM. In this review, current knowledge of the anatomical and physiological characteristics of LNs, as well as the molecular mechanisms of LNM, are described. The clinical significance of LNM in different anatomical sites is summarized, including the roles of LNM playing in staging, prognostic prediction, and treatment selection for patients with various types of cancers. And the novel exploration and academic disputes of strategies for recognition, diagnosis, and therapeutic interventions of metastatic LNs are also discussed.

Microfluidic design in single-cell sequencing and application to cancer precision medicine

Single-cell sequencing (SCS) is a crucial tool to reveal the genetic and functional heterogeneity of tumors, providing unique insights into the clonal evolution, microenvironment, drug resistance, and metastatic progression of cancers. Microfluidics is a critical component of many SCS technologies and workflows, conferring advantages in throughput, economy, and automation. Here, we review the current landscape of microfluidic architectures and sequencing techniques for single-cell omics analysis and highlight how these have enabled recent applications in oncology research. We also discuss the challenges and the promise of microfluidics-based single-cell analysis in the future of precision oncology.

Decoding Cancer Evolution: Integrating Genetic and Non-Genetic Insights

The development of cancer begins with cells transitioning from their multicellular nature to a state akin to unicellular organisms. This shift leads to a breakdown in the crucial regulators inherent to multicellularity, resulting in the emergence of diverse cancer cell subpopulations that have enhanced adaptability. The presence of different cell subpopulations within a tumour, known as intratumoural heterogeneity (ITH), poses challenges for cancer treatment. In this review, we delve into the dynamics of the shift from multicellularity to unicellularity during cancer onset and progression. We highlight the role of genetic and non-genetic factors, as well as tumour microenvironment, in promoting ITH and cancer evolution. Additionally, we shed light on the latest advancements in omics technologies that allow for in-depth analysis of tumours at the single-cell level and their spatial organization within the tissue. Obtaining such detailed information is crucial for deepening our understanding of the diverse evolutionary paths of cancer, allowing for the development of effective therapies targeting the key drivers of cancer evolution.

Deciphering the landscape of transcriptional heterogeneity across cancer

By integrating scRNA-seq datasets across 77 studies and 24 cancer types, in Nature, Gavish et al. uncover recurrent patterns of gene expression that explain a significant proportion of transcriptomic heterogeneity observed in cancer and explore their functional significance.

Machine learning based lineage tree reconstruction improved with knowledge of higher level relationships between cells and genomic barcodes

Tracking cells as they divide and progress through differentiation is a fundamental step in understanding many biological processes, such as the development of organisms and progression of diseases. In this study, we investigate a machine learning approach to reconstruct lineage trees in experimental systems based on mutating synthetic genomic barcodes. We refine previously proposed methodology by embedding information of higher level relationships between cells and single-cell barcode values into a feature space. We test performance of the algorithm on shallow trees (up to 100 cells) and deep trees (up to 10 000 cells). Our proposed algorithm can improve tree reconstruction accuracy in comparison to reconstructions based on a maximum parsimony method, but this comes at a higher computational time requirement.

Clinical and translational relevance of intratumor heterogeneity

Intratumor heterogeneity (ITH) is a driver of tumor evolution and a main cause of therapeutic resistance. Despite its importance, measures of ITH are still not incorporated into clinical practice. Consequently, standard treatment is frequently ineffective for patients with heterogeneous tumors as changes to treatment regimens are made only after recurrence and disease progression. More effective combination therapies require a mechanistic understanding of ITH and ways to assess it in clinical samples. The growth of technologies enabling the spatially intact analysis of tumors at the single-cell level and the development of sophisticated preclinical models give us hope that ITH will not simply be used as a predictor of a poor outcome but will guide treatment decisions from diagnosis through treatment.

Spectral top-down recovery of latent tree models

Modeling the distribution of high-dimensional data by a latent tree graphical model is a prevalent approach in multiple scientific domains. A common task is to infer the underlying tree structure, given only observations of its terminal nodes. Many algorithms for tree recovery are computationally intensive, which limits their applicability to trees of moderate size. For large trees, a common approach, termed divide-and-conquer, is to recover the tree structure in two steps. First, separately recover the structure of multiple, possibly random subsets of the terminal nodes. Second, merge the resulting subtrees to form a full tree. Here, we develop spectral top-down recovery (STDR), a deterministic divide-and-conquer approach to infer large latent tree models. Unlike previous methods, STDR partitions the terminal nodes in a non random way, based on the Fiedler vector of a suitable Laplacian matrix related to the observed nodes. We prove that under certain conditions, this partitioning is consistent with the tree structure. This, in turn, leads to a significantly simpler merging procedure of the small subtrees. We prove that STDR is statistically consistent and bound the number of samples required to accurately recover the tree with high probability. Using simulated data from several common tree models in phylogenetics, we demonstrate that STDR has a significant advantage in terms of runtime, with improved or similar accuracy.

Phylogenetic inference from single-cell RNA-seq data

Tumors are comprised of subpopulations of cancer cells that harbor distinct genetic profiles and phenotypes that evolve over time and during treatment. By reconstructing the course of cancer evolution, we can understand the acquisition of the malignant properties that drive tumor progression. Unfortunately, recovering the evolutionary relationships of individual cancer cells linked to their phenotypes remains a difficult challenge. To address this need, we have developed PhylinSic, a method that reconstructs the phylogenetic relationships among cells linked to their gene expression profiles from single cell RNA-sequencing (scRNA-Seq) data. This method calls nucleotide bases using a probabilistic smoothing approach and then estimates a phylogenetic tree using a Bayesian modeling algorithm. We showed that PhylinSic identified evolutionary relationships underpinning drug selection and metastasis and was sensitive enough to identify subclones from genetic drift. We found that breast cancer tumors resistant to chemotherapies harbored multiple genetic lineages that independently acquired high K-Ras and β-catenin, suggesting that therapeutic strategies may need to control multiple lineages to be durable. These results demonstrated that PhylinSic can reconstruct evolution and link the genotypes and phenotypes of cells across monophyletic tumors using scRNA-Seq.

Improving cellular phylogenies through the integrated use of mutation order and optimality principles

The study of tumor evolution is being revolutionalized by single-cell sequencing technologies that survey the somatic variation of cancer cells. In these endeavors, reliable inference of the evolutionary relationship of single cells is a key step. However, single-cell sequences contain many errors and missing bases, which necessitate advancing standard molecular phylogenetics approaches for applications in analyzing these datasets. We have developed a computational approach that integratively applies standard phylogenetic optimality principles and patterns of co-occurrence of sequence variations to produce more expansive and accurate cellular phylogenies from single-cell sequence datasets. We found the new approach to also perform well for CRISPR/Cas9 genome editing datasets, suggesting that it can be useful for various applications. We apply the new approach to some empirical datasets to showcase its use for reconstructing recurrent mutations and mutational reversals as well as for phylodynamics analysis to infer metastatic cell migrations between tumors.

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.

Single-cell methylation sequencing data reveal succinct metastatic migration histories and tumor progression models

Recent studies exploring the impact of methylation in tumor evolution suggest that although the methylation status of many of the CpG sites are preserved across distinct lineages, others are altered as the cancer progresses. Because changes in methylation status of a CpG site may be retained in mitosis, they could be used to infer the progression history of a tumor via single-cell lineage tree reconstruction. In this work, we introduce the first principled distance-based computational method, Sgootr, for inferring a tumor's single-cell methylation lineage tree and for jointly identifying lineage-informative CpG sites that harbor changes in methylation status that are retained along the lineage. We apply Sgootr on single-cell bisulfite-treated whole-genome sequencing data of multiregionally sampled tumor cells from nine metastatic colorectal cancer patients, as well as multiregionally sampled single-cell reduced-representation bisulfite sequencing data from a glioblastoma patient. We show that the tumor lineages constructed reveal a simple model underlying tumor progression and metastatic seeding. A comparison of Sgootr against alternative approaches shows that Sgootr can construct lineage trees with fewer migration events and with more in concordance with the sequential-progression model of tumor evolution, with a running time a fraction of that used in prior studies. Lineage-informative CpG sites identified by Sgootr are in inter-CpG island (CGI) regions, as opposed to intra-CGIs, which have been the main regions of interest in genomic methylation-related analyses.

Cellular barcoding tracks heterogeneous clones through selective pressures and phenotypic transitions

Genomic DNA barcoding has emerged as a sensitive and flexible tool to measure the fates of clonal subpopulations within a heterogeneous cancer cell population. Coupling cellular barcoding with single-cell transcriptomics permits the longitudinal analysis of molecular mechanisms with detailed clone-level resolution. Numerous recent studies have employed these tools to track clonal cell states in cancer progression and treatment response. With these new technologies comes the opportunity to examine longstanding questions about the origins and contributions of tumor cell heterogeneity and the roles of selection and phenotypic plasticity in disease progression and treatment.

Immune Resistance Mechanisms and the Road to Personalized Immunotherapy

What does the future of cancer immunotherapy look like and how do we get there? Find out where we've been and where we're headed in A Report on Resistance: The Road to personalized immunotherapy.

Applications of single-cell RNA sequencing in drug discovery and development

Single-cell technologies, particularly single-cell RNA sequencing (scRNA-seq) methods, together with associated computational tools and the growing availability of public data resources, are transforming drug discovery and development. New opportunities are emerging in target identification owing to improved disease understanding through cell subtyping, and highly multiplexed functional genomics screens incorporating scRNA-seq are enhancing target credentialling and prioritization. ScRNA-seq is also aiding the selection of relevant preclinical disease models and providing new insights into drug mechanisms of action. In clinical development, scRNA-seq can inform decision-making via improved biomarker identification for patient stratification and more precise monitoring of drug response and disease progression. Here, we illustrate how scRNA-seq methods are being applied in key steps in drug discovery and development, and discuss ongoing challenges for their implementation in the pharmaceutical industry.

Insights into the metastatic cascade through research autopsies

Metastasis is a complex process and the leading cause of cancer-related death globally. Recent studies have demonstrated that genomic sequencing data from paired primary and metastatic tumours can be used to trace the evolutionary origins of cells responsible for metastasis. This approach has yielded new insights into the genomic alterations that engender metastatic potential, and the mechanisms by which cancer spreads. Given that the reliability of these approaches is contingent upon how representative the samples are of primary and metastatic tumour heterogeneity, we review insights from studies that have reconstructed the evolution of metastasis within the context of their cohorts and designs. We discuss the role of research autopsies in achieving the comprehensive sampling necessary to advance the current understanding of metastasis.

Single-cell technologies in multiple myeloma: new insights into disease pathogenesis and translational implications

Multiple myeloma (MM) is a hematological malignancy characterized by clonal proliferation of plasma cells. Although therapeutic advances have been made to improve clinical outcomes and to prolong patients' survival in the past two decades, MM remains largely incurable. Single-cell sequencing (SCS) is a powerful method to dissect the cellular and molecular landscape at single-cell resolution, instead of providing averaged results. The application of single-cell technologies promises to address outstanding questions in myeloma biology and has revolutionized our understanding of the inter- and intra-tumor heterogeneity, tumor microenvironment, and mechanisms of therapeutic resistance in MM. In this review, we summarize the recently developed SCS methodologies and latest MM research progress achieved by single-cell profiling, including information regarding the cancer and immune cell landscapes, tumor heterogeneities, underlying mechanisms and biomarkers associated with therapeutic response and resistance. We also discuss future directions of applying transformative SCS approaches with contribution to clinical translation.