Emergence of a High-Plasticity Cell State during Lung Cancer Evolution

Senescent fibroblasts in the tumor stroma rewire lung cancer metabolism and plasticity

Senescence has been demonstrated to either inhibit or promote tumorigenesis. Resolving this paradox requires spatial mapping and functional characterization of senescent cells in the native tumor niche. Here, we identified senescent p16 Ink4a + cancer-associated fibroblasts with a secretory phenotype that promotes fatty acid uptake and utilization by aggressive lung adenocarcinoma driven by Kras and p53 mutations. Furthermore, rewiring of lung cancer metabolism by p16 Ink4a + cancer-associated fibroblasts also altered tumor cell identity to a highly plastic/dedifferentiated state associated with progression in murine and human LUAD. Our ex vivo senolytic screening platform identified XL888, a HSP90 inhibitor, that cleared p16 Ink4a + cancer-associated fibroblasts in vivo. XL888 administration after establishment of advanced lung adenocarcinoma significantly reduced tumor burden concurrent with the loss of plastic tumor cells. Our study identified a druggable component of the tumor stroma that fulfills the metabolic requirement of tumor cells to acquire a more aggressive phenotype.

KLF4 is an epigenetically modulated, context-dependent tumor suppressor

The epigenetic layer of regulation has become increasingly relevant in the research focused on tumor suppressors. KLF4 is a well-described zinc-finger transcription factor, mainly known for its role in the acquisition of cell pluripotency. Here we report and describe the most relevant epigenetic regulation mechanisms that affect KLF4 expression in tumors. CpG island methylation emerges as the most common mechanism in several tumors including lung adenocarcinoma, hepatocellular carcinoma, non-Hodgkin lymphomas, among others. Further layers of regulation represented by histone methylation and acetylation and by non-coding RNAs are described. Overall, KLF4 emerges as a crucial target in the fight against cancer.

LncRNAs and asymmetric cell division: the epigenetic mechanisms

Asymmetric cell division (ACD) plays a pivotal role in development, tissue homeostasis, and stem cell maintenance. Emerging evidence suggests that long non-coding RNAs (lncRNAs) are key regulators of ACD, orchestrating the intricate molecular machinery that governs cell fate determination. This review summarizes current literature to elucidate the diverse roles of lncRNAs in modulating ACD across various biological contexts. The regulatory mechanisms of asymmetric cell division mediated by lncRNAs, including their interactions with protein effectors, epigenetic regulation, and subcellular localization are explored. Additionally, we discuss the implications of dysregulated lncRNAs in mediating ACD that lead to tumorigenesis. By integrating findings from diverse experimental models and cell types, this review provides insights into the multifaceted roles of lncRNAs in governing asymmetric cell division, shedding light on fundamental biological processes. Further research in this area may lead to the development of novel therapies targeting dysregulated lncRNAs to restore proper cell division and function. The knowledge of lncRNAs regulating ACD could potentially revolutionize the field of regenerative medicine and cancer therapy by targeting specific lncRNAs involved in ACD. By unraveling the complex interactions between lncRNAs and cellular processes, the potential novel opportunities for precision medicine approaches may be uncovered.

Early-stage lung cancer is driven by a transitional cell state dependent on a KRAS-ITGA3-SRC axis

Glycine-12 mutations in the GTPase KRAS (KRASG12) are an initiating event for development of lung adenocarcinoma (LUAD). KRASG12 mutations promote cell-intrinsic rewiring of alveolar type-II progenitor (AT2) cells, but to what extent such changes interplay with lung homeostasis and cell fate pathways is unclear. Here, we generated single-cell RNA-seq (scRNA-seq) profiles from AT2-mesenchyme organoid co-cultures, mice, and stage-IA LUAD patients, identifying conserved regulators of AT2 transcriptional dynamics and defining the impact of KRASG12D mutation with temporal resolution. In AT2WT organoids, we found a transient injury/plasticity state preceding AT2 self-renewal and AT1 differentiation. Early-stage AT2KRAS cells exhibited perturbed gene expression dynamics, most notably retention of the injury/plasticity state. The injury state in AT2KRAS cells of patients, mice, and organoids was distinguishable from AT2WT states via altered receptor expression, including co-expression of ITGA3 and SRC. The combination of clinically relevant KRASG12D and SRC inhibitors impaired AT2KRAS organoid growth. Together, our data show that an injury/plasticity state essential for lung repair is co-opted during AT2 self-renewal and LUAD initiation, suggesting that early-stage LUAD may be susceptible to interventions that target specifically the oncogenic nature of this cell state.

Aging limits stemness and tumorigenesis in the lung by reprogramming iron homeostasis

Aging is associated with a decline in the number and fitness of adult stem cells 1-4 . Aging-associated loss of stemness is posited to suppress tumorigenesis 5,6 , but this hypothesis has not been tested in vivo . Here, using physiologically aged autochthonous genetically engineered mouse models and primary cells 7,8 , we demonstrate aging suppresses lung cancer initiation and progression by degrading stemness of the alveolar cell of origin. This phenotype is underpinned by aging-associated induction of the transcription factor NUPR1 and its downstream target lipocalin-2 in the cell of origin in mice and humans, leading to a functional iron insufficiency in the aged cells. Genetic inactivation of the NUPR1-lipocalin-2 axis or iron supplementation rescue stemness and promote tumorigenic potential of aged alveolar cells. Conversely, targeting the NUPR1- lipocalin-2 axis is detrimental to young alveolar cells via induction of ferroptosis. We find that aging-associated DNA hypomethylation at specific enhancer sites associates with elevated NUPR1 expression, which is recapitulated in young alveolar cells by inhibition of DNA methylation. We uncover that aging drives a functional iron insufficiency, which leads to loss of stemness and tumorigenesis, but promotes resistance to ferroptosis. These findings have significant implications for the therapeutic modulation of cellular iron homeostasis in regenerative medicine and in cancer prevention. Furthermore, our findings are consistent with a model whereby most human cancers initiate in young individuals, revealing a critical window for such cancer prevention efforts.

Smoking-induced CCNA2 expression promotes lung adenocarcinoma tumorigenesis by boosting AT2/AT2-like cell differentiation

Lung adenocarcinoma (LUAD), a type of non-small cell lung cancer (NSCLC), originates from not only bronchial epithelial cells but also alveolar type 2 (AT2) cells, which could differentiate into AT2-like cells. AT2-like cells function as cancer stem cells (CSCs) of LUAD tumorigenesis to give rise to adenocarcinoma. However, the mechanism underlying AT2 cell differentiation into AT2-like cells in LUAD remains unknown. We analyze genes differentially expressed and genes with significantly different survival curves in LUAD, and the combination of these two analyses yields 147 differential genes, in which 14 differentially expressed genes were enriched in cell cycle pathway. We next analyze the protein levels of these genes in LUAD and find that Cyclin-A2 (CCNA2) is closely associated with LUAD tumorigenesis. Unexpectedly, high CCNA2 expression in LUAD is restrictedly associated with smoking and independent of other driver mutations. Single-cell sequencing analyses reveal that CCNA2 is predominantly involved in AT2-like cell differentiation, while inhibition of CCNA2 significantly reverses smoking-induced AT2-like cell differentiation. Mechanistically, CCNA2 binding to CDK2 phosphorylates the AXIN1 complex, which in turn induces ubiquitination-dependent degradation of β-catenin and inhibits the WNT signaling pathway, thereby failing AT2 cell maintenance. These results uncover smoking-induced CCNA2 overexpression and subsequent WNT/β-catenin signaling inactivation as a hitherto uncharacterized mechanism controlling AT2 cell differentiation and LUAD tumorigenesis.

Tumor initiation and early tumorigenesis: molecular mechanisms and interventional targets

Tumorigenesis is a multistep process, with oncogenic mutations in a normal cell conferring clonal advantage as the initial event. However, despite pervasive somatic mutations and clonal expansion in normal tissues, their transformation into cancer remains a rare event, indicating the presence of additional driver events for progression to an irreversible, highly heterogeneous, and invasive lesion. Recently, researchers are emphasizing the mechanisms of environmental tumor risk factors and epigenetic alterations that are profoundly influencing early clonal expansion and malignant evolution, independently of inducing mutations. Additionally, clonal evolution in tumorigenesis reflects a multifaceted interplay between cell-intrinsic identities and various cell-extrinsic factors that exert selective pressures to either restrain uncontrolled proliferation or allow specific clones to progress into tumors. However, the mechanisms by which driver events induce both intrinsic cellular competency and remodel environmental stress to facilitate malignant transformation are not fully understood. In this review, we summarize the genetic, epigenetic, and external driver events, and their effects on the co-evolution of the transformed cells and their ecosystem during tumor initiation and early malignant evolution. A deeper understanding of the earliest molecular events holds promise for translational applications, predicting individuals at high-risk of tumor and developing strategies to intercept malignant transformation.

ALCAM-mediated cDC1 CD8 T cells interactions are suppressed in advanced lung tumors

Conventional type 1 dendritic cells (cDC1) are critical regulators of anti-tumoral T-cell responses. The structure and abundance of intercellular contacts between cDC1 and CD8 T cells in cancer tissues is important to determine the outcome of the T-cell response. However, the molecular determinants controlling the stability of cDC1-CD8 interactions during cancer progression remain poorly investigated. Here, we generated a genetic model of non-small cell lung cancer crossed to a fluorescent cDC1 reporter (KP-XCR1venus) to allow the detection of cDC1-CD8T cell clusters in tumor tissues across tumor stages. We found that cDC1-CD8 clusters are abundant and productive at the early stages of tumor development but progressively diminish in advanced tumors. Transcriptional profiling and flow cytometry identified the adhesion molecule ALCAM/CD166 (Activated Leukocyte Cell Adhesion Molecule, ligand of CD6) as highly expressed by lung cDC1 and significantly downregulated in advanced tumors. Analysis of human datasets indicated that ALCAM is downregulated in non-small cell lung cancer and its expression correlates to better prognosis. Mechanistically, triggering ALCAM on lung cDC1 induces cytoskeletal remodeling and contact formation whereas its blockade prevents T-cell activation. Together, our results indicate that ALCAM is important to stabilize cDC1-CD8 interactions at early tumor stages, while its loss in advanced tumors contributes to immune evasion.

A developmental constraint model of cancer cell states and tumor heterogeneity

Cancer is a disease that stems from a fundamental liability inherent to multicellular life forms in which an individual cell is capable of reneging on the interests of the collective organism. Although cancer is commonly described as an evolutionary process, a less appreciated aspect of tumorigenesis may be the constraints imposed by the organism's developmental programs. Recent work from single-cell transcriptomic analyses across a range of cancer types has revealed the recurrence, plasticity, and co-option of distinct cellular states among cancer cell populations. Here, we note that across diverse cancer types, the observed cell states are proximate within the developmental hierarchy of the cell of origin. We thus posit a model by which cancer cell states are directly constrained by the organism's "developmental map." According to this model, a population of cancer cells traverses the developmental map, thereby generating a heterogeneous set of states whose interactions underpin emergent tumor behavior.

Phenotypic noise and plasticity in cancer evolution

Non-genetic alterations can produce changes in a cell's phenotype. In cancer, these phenomena can influence a cell's fitness by conferring access to heritable, beneficial phenotypes. Herein, we argue that current discussions of 'phenotypic plasticity' in cancer evolution ignore a salient feature of the original definition: namely, that it occurs in response to an environmental change. We suggest 'phenotypic noise' be used to distinguish non-genetic changes in phenotype that occur independently from the environment. We discuss the conceptual and methodological techniques used to identify these phenomena during cancer evolution. We propose that the distinction will guide efforts to define mechanisms of phenotype change, accelerate translational work to manipulate phenotypes through treatment, and, ultimately, improve patient outcomes.

NANOS1 restricts oral cancer cell motility and TGF-ß signaling

Oral squamous cell carcinoma (OSCC) is the most frequent type of cancer of the head and neck area accounting for approx. 377,000 new cancer cases every year. The epithelial-to-mesenchymal transition (EMT) program plays an important role in OSCC progression and metastasis therefore contributing to a poor prognosis in patients with advanced disease. Transforming growth factor beta (TGF-ß) is a powerful inducer of EMT thereby increasing cancer cell aggressiveness. Here, we aimed at identifying RNA-binding proteins (RBPs) that affect TGF-ß-induced EMT. To this end we treated oral cancer cells with TGF-ß and identified a total of 643 significantly deregulated protein-coding genes in response to TGF-ß. Of note, 19 genes encoded RBPs with NANOS1 being the most downregulated RBP. Subsequent cellular studies demonstrated a strong inhibitory effect of NANOS1 on migration and invasion of SAS oral cancer cells. Further mechanistic studies revealed an interaction of NANOS1 with the TGF-ß receptor 1 (TGFBR1) mRNA, leading to increased decay of this transcript and a reduced TGFBR1 protein expression, thereby preventing downstream TGF-ß/SMAD signaling. In summary, we identified NANOS1 as negative regulator of TGF-ß signaling in oral cancer cells.

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.

Integrative Single Cell Atlas Revealed Intratumoral Heterogeneity Generation from an Adaptive Epigenetic Cell State in Human Bladder Urothelial Carcinoma

Intratumor heterogeneity (ITH) of bladder cancer (BLCA) contributes to therapy resistance and immune evasion affecting clinical prognosis. The molecular and cellular mechanisms contributing to BLCA ITH generation remain elusive. It is found that a TM4SF1-positive cancer subpopulation (TPCS) can generate ITH in BLCA, evidenced by integrative single cell atlas analysis. Extensive profiling of the epigenome and transcriptome of all stages of BLCA revealed their evolutionary trajectories. Distinct ancestor cells gave rise to low-grade noninvasive and high-grade invasive BLCA. Epigenome reprograming led to transcriptional heterogeneity in BLCA. During early oncogenesis, epithelial-to-mesenchymal transition generated TPCS. TPCS has stem-cell-like properties and exhibited transcriptional plasticity, priming the development of transcriptionally heterogeneous descendent cell lineages. Moreover, TPCS prevalence in tumor is associated with advanced stage cancer and poor prognosis. The results of this study suggested that bladder cancer interacts with its environment by acquiring a stem cell-like epigenomic landscape, which might generate ITH without additional genetic diversification.

Stem-cell states converge in multistage cutaneous squamous cell carcinoma development

Stem cells play a critical role in cancer development by contributing to cell heterogeneity, lineage plasticity, and drug resistance. We created gene expression networks from hundreds of mouse tissue samples (both normal and tumor) and integrated these with lineage tracing and single-cell RNA-seq, to identify convergence of cell states in premalignant tumor cells expressing markers of lineage plasticity and drug resistance. Two of these cell states representing multilineage plasticity or proliferation were inversely correlated, suggesting a mutually exclusive relationship. Treatment of carcinomas in vivo with chemotherapy repressed the proliferative state and activated multilineage plasticity whereas inhibition of differentiation repressed plasticity and potentiated responses to cell cycle inhibitors. Manipulation of this cell state transition point may provide a source of potential combinatorial targets for cancer therapy.

Decoding the principle of cell-fate determination for its reverse control

Understanding and manipulating cell fate determination is pivotal in biology. Cell fate is determined by intricate and nonlinear interactions among molecules, making mathematical model-based quantitative analysis indispensable for its elucidation. Nevertheless, obtaining the essential dynamic experimental data for model development has been a significant obstacle. However, recent advancements in large-scale omics data technology are providing the necessary foundation for developing such models. Based on accumulated experimental evidence, we can postulate that cell fate is governed by a limited number of core regulatory circuits. Following this concept, we present a conceptual control framework that leverages single-cell RNA-seq data for dynamic molecular regulatory network modeling, aiming to identify and manipulate core regulatory circuits and their master regulators to drive desired cellular state transitions. We illustrate the proposed framework by applying it to the reversion of lung cancer cell states, although it is more broadly applicable to understanding and controlling a wide range of cell-fate determination processes.

Extracellular vesicles in non-small cell lung cancer stemness and clinical applications

Non-small cell lung carcinoma (NSCLC) accounts for 85% of lung cancers, the leading cause of cancer associated deaths in the US and worldwide. Within NSCLC tumors, there is a subpopulation of cancer cells termed cancer stem cells (CSCs) which exhibit stem-like properties that drive NSCLC progression, metastasis, relapse, and therapeutic resistance. Extracellular vesicles (EVs) are membrane-bound nanoparticles secreted by cells that carry vital messages for short- and long-range intercellular communication. Numerous studies have implicated NSCLC CSC-derived EVs in the factors associated with NSCLC lethality. In this review, we have discussed mechanisms of EV-directed cross-talk between CSCs and cells of the tumor microenvironment that promote stemness, tumor progression and metastasis in NSCLC. The mechanistic studies discussed herein have provided insights for developing novel NSCLC diagnostic and prognostic biomarkers and strategies to therapeutically target the NSCLC CSC niche.

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.

A novel bioinformatic approach reveals cooperation between Cancer/Testis genes in basal-like breast tumors

Breast cancer is the most prevalent type of cancer in women worldwide. Within breast tumors, the basal-like subtype has the worst prognosis, prompting the need for new tools to understand, detect, and treat these tumors. Certain germline-restricted genes show aberrant expression in tumors and are known as Cancer/Testis genes; their misexpression has diagnostic and therapeutic applications. Here we designed a new bioinformatic approach to examine Cancer/Testis gene misexpression in breast tumors. We identify several new markers in Luminal and HER-2 positive tumors, some of which predict response to chemotherapy. We then use machine learning to identify the two Cancer/Testis genes most associated with basal-like breast tumors: HORMAD1 and CT83. We show that these genes are expressed by tumor cells and not by the microenvironment, and that they are not expressed by normal breast progenitors; in other words, their activation occurs de novo. We find these genes are epigenetically repressed by DNA methylation, and that their activation upon DNA demethylation is irreversible, providing a memory of past epigenetic disturbances. Simultaneous expression of both genes in breast cells in vitro has a synergistic effect that increases stemness and activates a transcriptional profile also observed in double-positive tumors. Therefore, we reveal a functional cooperation between Cancer/Testis genes in basal breast tumors; these findings have consequences for the understanding, diagnosis, and therapy of the breast tumors with the worst outcomes.

Advances in molecular targeted drugs in combination with CAR-T cell therapy for hematologic malignancies

Molecular targeted drugs and chimeric antigen receptor (CAR) T cell therapy represent specific biological treatments that have significantly improved the efficacy of treating hematologic malignancies. However, they face challenges such as drug resistance and recurrence after treatment. Combining molecular targeted drugs and CAR-T cells could regulate immunity, improve tumor microenvironment (TME), promote cell apoptosis, and enhance sensitivity to tumor cell killing. This approach might provide a dual coordinated attack on cancer cells, effectively eliminating minimal residual disease and overcoming therapy resistance. Moreover, molecular targeted drugs can directly or indirectly enhance the anti-tumor effect of CAR-T cells by inducing tumor target antigen expression, reversing CAR-T cell exhaustion, and reducing CAR-T cell associated toxic side effects. Therefore, combining molecular targeted drugs with CAR-T cells is a promising and novel tactic for treating hematologic malignancies. In this review article, we focus on analyzing the mechanism of therapy resistance and its reversal of CAR-T cell therapy resistance, as well as the synergistic mechanism, safety, and future challenges in CAR-T cell therapy in combination with molecular targeted drugs. We aim to explore the benefits of this combination therapy for patients with hematologic malignancies and provide a rationale for subsequent clinical studies.

Single-cell and single-nucleus RNA-sequencing from paired normal-adenocarcinoma lung samples provide both common and discordant biological insights

Whether single-cell RNA-sequencing (scRNA-seq) captures the same biological information as single-nucleus RNA-sequencing (snRNA-seq) remains uncertain and likely to be context-dependent. Herein, a head-to-head comparison was performed in matched normal-adenocarcinoma human lung samples to assess biological insights derived from scRNA-seq versus snRNA-seq and better understand the cellular transition that occurs from normal to tumoral tissue. Here, the transcriptome of 160,621 cells/nuclei was obtained. In non-tumor lung, cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former (81.5%) and epithelial cells (69.9%) in the later. Similar results were observed in adenocarcinomas, in addition to an overall increase in cell type heterogeneity and a greater prevalence of copy number variants in cells of epithelial origin, which suggests malignant assignment. The cell type transition that occurs from normal lung tissue to adenocarcinoma was not always concordant whether cells or nuclei were examined. As expected, large differential expression of the whole-cell and nuclear transcriptome was observed, but cell-type specific changes of paired normal and tumor lung samples revealed a set of common genes in the cells and nuclei involved in cancer-related pathways. In addition, we showed that the ligand-receptor interactome landscape of lung adenocarcinoma was largely different whether cells or nuclei were evaluated. Immune cell depletion in fresh specimens partly mitigated the difference in cell type composition observed between cells and nuclei. However, the extra manipulations affected cell viability and amplified the transcriptional signatures associated with stress responses. In conclusion, research applications focussing on mapping the immune landscape of lung adenocarcinoma benefit from scRNA-seq in fresh samples, whereas snRNA-seq of frozen samples provide a low-cost alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content.

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.

Cell-intrinsic and microenvironmental determinants of metastatic colonization

Cancer metastasis is a biologically complex process that remains a major challenge in the oncology clinic, accounting for nearly all of the mortality associated with malignant neoplasms. To establish metastatic growths, carcinoma cells must disseminate from the primary tumour, survive in unfamiliar tissue microenvironments, re-activate programs of proliferation, and escape innate and adaptive immunosurveillance. The entire process is extremely inefficient and can occur over protracted timescales, yielding only a vanishingly small number of carcinoma cells that are able to complete all of the required steps. Here we review both the cancer-cell-intrinsic mechanisms and microenvironmental interactions that enable metastatic colonization. In particular, we highlight recent work on the behaviour of already-disseminated tumour cells, since meaningful progress in treating metastatic disease will clearly require a better understanding of the cells that spawn metastases, which generally have disseminated by the time of initial diagnosis.

Adeno-to-squamous transition drives resistance to KRAS inhibition in LKB1 mutant lung cancer

KRASG12C inhibitors (adagrasib and sotorasib) have shown clinical promise in targeting KRASG12C-mutated lung cancers; however, most patients eventually develop resistance. In lung patients with adenocarcinoma with KRASG12C and STK11/LKB1 co-mutations, we find an enrichment of the squamous cell carcinoma gene signature in pre-treatment biopsies correlates with a poor response to adagrasib. Studies of Lkb1-deficient KRASG12C and KrasG12D lung cancer mouse models and organoids treated with KRAS inhibitors reveal tumors invoke a lineage plasticity program, adeno-to-squamous transition (AST), that enables resistance to KRAS inhibition. Transcriptomic and epigenomic analyses reveal ΔNp63 drives AST and modulates response to KRAS inhibition. We identify an intermediate high-plastic cell state marked by expression of an AST plasticity signature and Krt6a. Notably, expression of the AST plasticity signature and KRT6A at baseline correlates with poor adagrasib responses. These data indicate the role of AST in KRAS inhibitor resistance and provide predictive biomarkers for KRAS-targeted therapies in lung cancer.

An atlas of epithelial cell states and plasticity in lung adenocarcinoma

Understanding the cellular processes that underlie early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies1. Here we studied 246,102 single epithelial cells from 16 early-stage LUADs and 47 matched normal lung samples. Epithelial cells comprised diverse normal and cancer cell states, and diversity among cancer cells was strongly linked to LUAD-specific oncogenic drivers. KRAS mutant cancer cells showed distinct transcriptional features, reduced differentiation and low levels of aneuploidy. Non-malignant areas surrounding human LUAD samples were enriched with alveolar intermediate cells that displayed elevated KRT8 expression (termed KRT8+ alveolar intermediate cells (KACs) here), reduced differentiation, increased plasticity and driver KRAS mutations. Expression profiles of KACs were enriched in lung precancer cells and in LUAD cells and signified poor survival. In mice exposed to tobacco carcinogen, KACs emerged before lung tumours and persisted for months after cessation of carcinogen exposure. Moreover, they acquired Kras mutations and conveyed sensitivity to targeted KRAS inhibition in KAC-enriched organoids derived from alveolar type 2 (AT2) cells. Last, lineage-labelling of AT2 cells or KRT8+ cells following carcinogen exposure showed that KACs are possible intermediates in AT2-to-tumour cell transformation. This study provides new insights into epithelial cell states at the root of LUAD development, and such states could harbour potential targets for prevention or intervention.

Cellular plasticity in reprogramming, rejuvenation and tumorigenesis: a pioneer TF perspective

The multistep process of in vivo reprogramming, mediated by the transcription factors (TFs) Oct4, Sox2, Klf4, and c-Myc (OSKM), holds great promise for the development of rejuvenating and regenerative strategies. However, most of the approaches developed so far are accompanied by a persistent risk of tumorigenicity. Here, we review the groundbreaking effects of in vivo reprogramming with a particular focus on rejuvenation and regeneration. We discuss how the activity of pioneer TFs generates cellular plasticity that may be critical for inducing not only reprogramming and regeneration, but also cancer initiation. Finally, we highlight how a better understanding of the uncoupled control of cellular identity, plasticity, and aging during reprogramming might pave the way to the development of rejuvenating/regenerating strategies in a nontumorigenic manner.

Lineage and ecology define liver tumor evolution in response to treatment

A tumor ecosystem constantly evolves over time in the face of immune predation or therapeutic intervention, resulting in treatment failure and tumor progression. Here, we present a single-cell transcriptome-based strategy to determine the evolution of longitudinal tumor biopsies from liver cancer patients by measuring cellular lineage and ecology. We construct a lineage and ecological score as joint dynamics of tumor cells and their microenvironments. Tumors may be classified into four main states in the lineage-ecological space, which are associated with clinical outcomes. Analysis of longitudinal samples reveals the evolutionary trajectory of tumors in response to treatment. We validate the lineage-ecology-based scoring system in predicting clinical outcomes using bulk transcriptomic data of additional cohorts of 716 liver cancer patients. Our study provides a framework for monitoring tumor evolution in response to therapeutic intervention.

Artificial Intelligence-Based Treatment Decisions: A New Era for NSCLC

Non-small cell lung cancer (NSCLC) is the leading cause of cancer-related mortality among women and men, in developed countries, despite the public health interventions including tobacco-free campaigns, screening and early detection methods, recent therapeutic advances, and ongoing intense research on novel antineoplastic modalities. Targeting oncogenic driver mutations and immune checkpoint inhibition has indeed revolutionized NSCLC treatment, yet there still remains the unmet need for robust and standardized predictive biomarkers to accurately inform clinical decisions. Artificial intelligence (AI) represents the computer-based science concerned with large datasets for complex problem-solving. Its concept has brought a paradigm shift in oncology considering its immense potential for improved diagnosis, treatment guidance, and prognosis. In this review, we present the current state of AI-driven applications on NSCLC management, with a particular focus on radiomics and pathomics, and critically discuss both the existing limitations and future directions in this field. The thoracic oncology community should not be discouraged by the likely long road of AI implementation into daily clinical practice, as its transformative impact on personalized treatment approaches is undeniable.

A cellular and spatial atlas of TP53 -associated tissue remodeling in lung adenocarcinoma

TP53 is the most frequently mutated gene across many cancers and is associated with shorter survival in lung adenocarcinoma (LUAD). To define how TP53 mutations affect the LUAD tumor microenvironment (TME), we constructed a multi-omic cellular and spatial tumor atlas of 23 treatment-naïve human lung tumors. We found that TP53 -mutant ( TP53 mut ) malignant cells lose alveolar identity and upregulate highly proliferative and entropic gene expression programs consistently across resectable LUAD patient tumors, genetically engineered mouse models, and cell lines harboring a wide spectrum of TP53 mutations. We further identified a multicellular tumor niche composed of SPP1 + macrophages and collagen-expressing fibroblasts that coincides with hypoxic, pro-metastatic expression programs in TP53 mut tumors. Spatially correlated angiostatic and immune checkpoint interactions, including CD274 - PDCD1 and PVR - TIGIT , are also enriched in TP53 mut LUAD tumors, which may influence response to checkpoint blockade therapy. Our methodology can be further applied to investigate mutation-specific TME changes in other cancers.

Functional analysis of cell plasticity using single-cell technologies

Metazoan organisms are heterocellular systems composed of hundreds of different cell types, which arise from an isogenic genome through differentiation. Cellular 'plasticity' further enables cells to alter their fate in response to exogenous cues and is involved in a variety of processes, such as wound healing, infection, and cancer. Recent advances in cellular model systems, high-dimensional single-cell technologies, and lineage tracing have sparked a renaissance in plasticity research. Here, we discuss the definition of cell plasticity, evaluate state-of-the-art model systems and techniques to study cell-fate dynamics, and explore the application of single-cell technologies to obtain functional insights into cell plasticity in healthy and diseased tissues. The integration of advanced biomimetic model systems, single-cell technologies, and high-throughput perturbation studies is enabling a new era of research into non-genetic plasticity in metazoan systems.

Alveolar Differentiation Drives Resistance to KRAS Inhibition in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD), commonly driven by KRAS mutations, is responsible for 7% of all cancer mortality. The first allele-specific KRAS inhibitors were recently approved in LUAD, but the clinical benefit is limited by intrinsic and acquired resistance. LUAD predominantly arises from alveolar type 2 (AT2) cells, which function as facultative alveolar stem cells by self-renewing and replacing alveolar type 1 (AT1) cells. Using genetically engineered mouse models, patient-derived xenografts, and patient samples, we found inhibition of KRAS promotes transition to a quiescent AT1-like cancer cell state in LUAD tumors. Similarly, suppressing Kras induced AT1 differentiation of wild-type AT2 cells upon lung injury. The AT1-like LUAD cells exhibited high growth and differentiation potential upon treatment cessation, whereas ablation of the AT1-like cells robustly improved treatment response to KRAS inhibitors. Our results uncover an unexpected role for KRAS in promoting intratumoral heterogeneity and suggest that targeting alveolar differentiation may augment KRAS-targeted therapies in LUAD.

SIGNIFICANCE

Treatment resistance limits response to KRAS inhibitors in LUAD patients. We find LUAD residual disease following KRAS targeting is composed of AT1-like cancer cells with the capacity to reignite tumorigenesis. Targeting the AT1-like cells augments responses to KRAS inhibition, elucidating a therapeutic strategy to overcome resistance to KRAS-targeted therapy. This article is featured in Selected Articles from This Issue, p. 201.

Exosomal lncCRLA is predictive for the evolvement and development of lung adenocarcinoma

Lung adenocarcinoma, the most common histological subtype of non-small cell lung cancer, exhibits heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Our team uncovers that lncRNA related to chemotherapy resistance in lung adenocarcinoma (lncCRLA) is preferentially expressed in lung adenocarcinoma cells with the mesenchymal phenotype. lncCRLA can not enhance chemotherapy resistance in lung adenocarcinoma due to its binding to RIPK1 in exosomes, which is released into intercellular media and transferred by exosomes from mesenchymal-like to epithelial-like cells. However, plasmatic lncCRLA corresponding to tissue lncCRLA functions as a preferred biomarker to reflect the response to chemotherapy and disease progression of lung adenocarcinoma. Through single-cell sequencing, RNA-Mutect technique and spatial transcriptomics, a handful of hybrid EMT cells with elevated lncCRLA are characterized as the origin of lung adenocarcinoma, which are indiscriminated from hybrid EMT cells by the in-depth sequencing. Plasmatic lncCRLA is properly predictive for the preinvasive lesion of lung adenocarcinoma that would evolve to invasive lesion. That notion is confirmed by a brand-new transgenic mouse model in which EMT is tracked by Cre and Dre system. Dasatinib is potential to hinder the spontaneous progression from preinvasive to invasive lesion of lung adenocarcinoma. Together, plasmatic lncCRLA is defined as a brand-new circulating biomarker to predict the occurrence and evolvement of lung adenocarcinoma, a light for early detection of lung adenocarcinoma.

Solid tumor growth depends on an intricate equilibrium of malignant cell states

Control of cell identity and number is central to tissue function, yet principles governing organization of malignant cells in tumor tissues remain poorly understood. Using mathematical modeling and candidate-based analysis, we discover primary and metastatic pancreatic ductal adenocarcinoma (PDAC) organize in a stereotypic pattern whereby PDAC cells responding to WNT signals (WNT-R) neighbor WNT-secreting cancer cells (WNT-S). Leveraging lineage-tracing, we reveal the WNT-R state is transient and gives rise to the WNT-S state that is highly stable and committed to organizing malignant tissue. We further show that a subset of WNT-S cells expressing the Notch ligand DLL1 form a functional niche for WNT-R cells. Genetic inactivation of WNT secretion or Notch pathway components, or cytoablation of the WNT-S state disrupts PDAC tissue organization, suppressing tumor growth and metastasis. This work indicates PDAC growth depends on an intricately controlled equilibrium of functionally distinct cancer cell states, uncovering a fundamental principle governing solid tumor growth and revealing new opportunities for therapeutic intervention.

Temporal evolution reveals bifurcated lineages in aggressive neuroendocrine small cell prostate cancer trans-differentiation

Trans-differentiation from an adenocarcinoma to a small cell neuroendocrine state is associated with therapy resistance in multiple cancer types. To gain insight into the underlying molecular events of the trans-differentiation, we perform a multi-omics time course analysis of a pan-small cell neuroendocrine cancer model (termed PARCB), a forward genetic transformation using human prostate basal cells and identify a shared developmental, arc-like, and entropy-high trajectory among all transformation model replicates. Further mapping with single cell resolution reveals two distinct lineages defined by mutually exclusive expression of ASCL1 or ASCL2. Temporal regulation by groups of transcription factors across developmental stages reveals that cellular reprogramming precedes the induction of neuronal programs. TFAP4 and ASCL1/2 feedback are identified as potential regulators of ASCL1 and ASCL2 expression. Our study provides temporal transcriptional patterns and uncovers pan-tissue parallels between prostate and lung cancers, as well as connections to normal neuroendocrine cell states.

A Mutation-driven oncofetal regression fuels phenotypic plasticity in colorectal cancer

Targeting cancer stem cells (CSCs) is crucial for effective cancer treatment 1 . However, the molecular mechanisms underlying resistance to LGR5 + CSCs depletion in colorectal cancer (CRC) 2,3 remain largely elusive. Here, we unveil the existence of a primitive cell state dubbed the oncofetal (OnF) state, which works in tandem with the LGR5 + stem cells (SCs) to fuel tumor evolution in CRC. OnF cells emerge early during intestinal tumorigenesis and exhibit features of lineage plasticity. Normally suppressed by the Retinoid X Receptor (RXR) in mature SCs, the OnF program is triggered by genetic deletion of the gatekeeper APC. We demonstrate that diminished RXR activity unlocks an epigenetic circuity governed by the cooperative action of YAP and AP1, leading to OnF reprogramming. This high-plasticity state is inherently resistant to conventional chemotherapies and its adoption by LGR5 + CSCs enables them to enter a drug-tolerant state. Furthermore, through phenotypic tracing and ablation experiments, we uncover a functional redundancy between the OnF and stem cell (SC) states and show that targeting both cellular states is essential for sustained tumor regression in vivo . Collectively, these findings establish a mechanistic foundation for developing effective combination therapies with enduring impact on CRC treatment.

An optimized protocol for the generation and monitoring of conditional orthotopic lung cancer in the KP mouse model using an adeno-associated virus vector compatible with biosafety level 1

BACKGROUND

The inducible Kras/p53 lung adenocarcinoma mouse model, which faithfully recapitulates human disease, is routinely initiated by the intratracheal instillation of a virus-based Cre recombinase delivery system. Handling virus-based delivery systems requires elevated biosafety levels, e.g., biosafety level 2 (BSL-2). However, in experimental animal research facilities, following exposure to viral vectors in a BSL-2 environment, rodents may not be reclassified to BSL-1 according to standard practice, preventing access to small animal micro-computed tomography (micro-CT) scanners that are typically housed in general access areas such as BSL-1 rooms. Therefore, our goal was to adapt the protocol so that the Cre-induced KP mouse model could be handled under BSL-1 conditions during the entire procedure.

RESULTS

The Kras-Lox-STOP-Lox-G12D/p53 flox/flox (KP)-based lung adenocarcinoma mouse model was activated by intratracheal instillation of either an adenoviral-based or a gutless, adeno-associated viral-based Cre delivery system. Tumor growth was monitored over time by micro-CT. We have successfully substituted the virus-based Cre delivery system with a commercially available, gutless, adeno-associated, Cre-expressing vector that allows the KP mouse model to be handled and imaged in a BSL-1 facility. By optimizing the anesthesia protocol and switching to a microscope-guided vector instillation procedure, productivity was increased and procedure-related complications were significantly reduced. In addition, repeated micro-CT analysis of individual animals allowed us to monitor tumor growth longitudinally, dramatically reducing the number of animals required per experiment. Finally, we documented the evolution of tumor volume for different doses, which revealed that individual tumor nodules induced by low-titer AAV-Cre transductions can be monitored over time by micro-CT.

CONCLUSION

Modifications to the anesthesia and instillation protocols increased the productivity of the original KP protocol. In addition, the switch to a gutless, adeno-associated, Cre-expressing vector allowed longitudinal monitoring of tumor growth under BSL-1 conditions, significantly reducing the number of animals required for an experiment, in line with the 3R principles.

Cellular dynamics as a marker of normal-to-cancer transition in human cells

Normal-to-cancer (NTC) transition is known to be closely associated to cell´s biomechanical properties which are dependent on the dynamics of the intracellular medium. This study probes different human cancer cells (breast, prostate and lung), concomitantly to their healthy counterparts, aiming at characterising the dynamical profile of water in distinct cellular locations, for each type of cell, and how it changes between normal and cancer states. An increased plasticity of the cytomatrix is observed upon normal-to-malignant transformation, the lung carcinoma cells displaying the highest flexibility followed by prostate and breast cancers. Also, lung cells show a distinct behaviour relative to breast and prostate, with a higher influence from hydration water motions and localised fast rotations upon NTC transformation. Quasielastic neutron scattering techniques allowed to accurately distinguish the different dynamical processes taking place within these highly heterogeneous cellular systems. The results thus obtained suggest that intracellular water dynamics may be regarded as a specific reporter of the cellular conditions-either healthy or malignant.

Activation of KrasG12D in Subset of Alveolar Type II Cells Enhances Cellular Plasticity in Lung Adenocarcinoma

We have previously identified alveolar type II cell as the cell-of-origin of KrasG12D-induced lung adenocarcinoma using cell lineage-specific inducible Cre mouse models. Using gain-of-function and loss-of-function genetic models, we discovered that active Notch signaling and low Sox2 levels dictate the ability of type II cells to proliferate and progress into lung adenocarcinoma upon KrasG12D activation. Here, we examine the phenotype of type II cells after Kras activation and find evidence for proliferation of cells that coexpress type I and type II markers. Three-dimensional organoid culture and transplantation studies determine that these dual-positive cells are highly plastic and tumor initiating in vivo. RNA sequencing analysis reveals that these dual-positive cells are enriched in Ras/MAPK, EGFR, and Notch pathways. Furthermore, the proliferation of these cells requires active Notch signaling and is inhibited by genetic/chemical Sox2 upregulation. Our findings could provide new therapeutic strategies to target KRAS-activated lung adenocarcinomas.

SIGNIFICANCE

Identification of progenitor like tumor-initiating cells in KRAS-mutant lung adenocarcinoma may allow development of novel targeted therapeutics.

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.

Pan-cancer spatially resolved single-cell analysis reveals the crosstalk between cancer-associated fibroblasts and tumor microenvironment

Cancer-associated fibroblasts (CAFs) are a heterogeneous cell population that plays a crucial role in remodeling the tumor microenvironment (TME). Here, through the integrated analysis of spatial and single-cell transcriptomics data across six common cancer types, we identified four distinct functional subgroups of CAFs and described their spatial distribution characteristics. Additionally, the analysis of single-cell RNA sequencing (scRNA-seq) data from three additional common cancer types and two newly generated scRNA-seq datasets of rare cancer types, namely epithelial-myoepithelial carcinoma (EMC) and mucoepidermoid carcinoma (MEC), expanded our understanding of CAF heterogeneity. Cell-cell interaction analysis conducted within the spatial context highlighted the pivotal roles of matrix CAFs (mCAFs) in tumor angiogenesis and inflammatory CAFs (iCAFs) in shaping the immunosuppressive microenvironment. In patients with breast cancer (BRCA) undergoing anti-PD-1 immunotherapy, iCAFs demonstrated heightened capacity in facilitating cancer cell proliferation, promoting epithelial-mesenchymal transition (EMT), and contributing to the establishment of an immunosuppressive microenvironment. Furthermore, a scoring system based on iCAFs showed a significant correlation with immune therapy response in melanoma patients. Lastly, we provided a web interface ( https://chenxisd.shinyapps.io/pancaf/ ) for the research community to investigate CAFs in the context of pan-cancer.

β-catenin/CBP activation of mTORC1 signaling promotes partial epithelial-mesenchymal states in head and neck cancer

Head and neck cancers, which include oral squamous cell carcinoma (OSCC) as a major subsite, exhibit cellular plasticity that includes features of an epithelial-mesenchymal transition (EMT), referred to as partial-EMT (p-EMT). To identify molecular mechanisms contributing to OSCC plasticity, we performed a multiphase analysis of single cell RNA sequencing (scRNAseq) data from human OSCC. This included a multiresolution characterization of cancer cell subgroups to identify pathways and cell states that are heterogeneously represented, followed by casual inference analysis to elucidate activating and inhibitory relationships between these pathways and cell states. This approach revealed signaling networks associated with hierarchical cell state transitions, which notably included an association between β-catenin-driven CREB-binding protein (CBP) activity and mTORC1 signaling. This network was associated with subpopulations of cancer cells that were enriched for markers of the p-EMT state and poor patient survival. Functional analyses revealed that β-catenin/CBP induced mTORC1 activity in part through the transcriptional regulation of a raptor-interacting protein, chaperonin containing TCP1 subunit 5 (CCT5). Inhibition of β-catenin-CBP activity through the use of the orally active small molecule, E7386, reduced the expression of CCT5 and mTORC1 activity in vitro, and inhibited p-EMT-associated markers and tumor development in a murine model of OSCC. Our study highlights the use of multiresolution network analyses of scRNAseq data to identify targetable signals for therapeutic benefit, thus defining an underappreciated association between β-catenin/CBP and mTORC1 signaling in head and neck cancer plasticity.

Development of two mouse strains conditionally expressing bright luciferases with distinct emission spectra as new tools for in vivo imaging

In vivo bioluminescence imaging (BLI) has been an invaluable noninvasive method to visualize molecular and cellular behaviors in laboratory animals. Bioluminescent reporter mice harboring luciferases for general use have been limited to a classical luciferase, Luc2, from Photinus pyralis, and have been extremely powerful for various in vivo studies. However, applicability of reporter mice for in vivo BLI could be further accelerated by increasing light intensity through the use of other luciferases and/or by improving the biodistribution of their substrates in the animal body. Here we created two Cre-dependent reporter mice incorporating luciferases oFluc derived from Pyrocoeli matsumurai and Akaluc, both of which had been reported previously to be brighter than Luc2 when using appropriate substrates; we then tested their bioluminescence in neural tissues and other organs in living mice. When expressed throughout the body, both luciferases emitted an intense yellow (oFluc) or far-red (Akaluc) light easily visible to the naked eye. oFluc and Akaluc were similarly bright in the pancreas for in vivo BLI; however, Akaluc was superior to oFluc for brain imaging, because its substrate, AkaLumine-HCl, was distributed to the brain more efficiently than the oFluc substrate, D-luciferin. We also demonstrated that the lights produced by oFluc and Akaluc were sufficiently spectrally distinct from each other for dual-color imaging in a single living mouse. Taken together, these novel bioluminescent reporter mice are an ideal source of cells with bright bioluminescence and may facilitate in vivo BLI of various tissues/organs for preclinical and biomedical research in combination with a wide variety of Cre-driver mice.

Alveolar differentiation drives resistance to KRAS inhibition in lung adenocarcinoma

Lung adenocarcinoma (LUAD), commonly driven by KRAS mutations, is responsible for 7% of all cancer mortality. The first allele-specific KRAS inhibitors were recently approved in LUAD, but clinical benefit is limited by intrinsic and acquired resistance. LUAD predominantly arises from alveolar type 2 (AT2) cells, which function as facultative alveolar stem cells by self-renewing and replacing alveolar type 1 (AT1) cells. Using genetically engineered mouse models, patient-derived xenografts, and patient samples we found inhibition of KRAS promotes transition to a quiescent AT1-like cancer cell state in LUAD tumors. Similarly, suppressing Kras induced AT1 differentiation of wild-type AT2 cells upon lung injury. The AT1-like LUAD cells exhibited high growth and differentiation potential upon treatment cessation, whereas ablation of the AT1-like cells robustly improved treatment response to KRAS inhibitors. Our results uncover an unexpected role for KRAS in promoting intra-tumoral heterogeneity and suggest targeting alveolar differentiation may augment KRAS-targeted therapies in LUAD.

Significance

Treatment resistance limits response to KRAS inhibitors in LUAD patients. We find LUAD residual disease following KRAS targeting is composed of AT1-like cancer cells with the capacity to reignite tumorigenesis. Targeting the AT1-like cells augments responses to KRAS inhibition, elucidating a therapeutic strategy to overcome resistance to KRAS-targeted therapy.

The circadian clock circuitry deconvolutes colorectal cancer and lung adenocarcinoma heterogeneity in a dynamic time-related framework

Increasing evidence imputes cancer progression and resistance to therapy to intra-tumor molecular heterogeneity set off by cancer cell plasticity. Re-activation of developmental programs strictly linked to epithelial-to-mesenchymal transition and gaining of stem cells properties are crucial in this setting. Many biological processes involved in cancer onset and progression show rhythmic fluctuations driven by the circadian clock circuitry. Novel cancer patient stratification tools taking into account the temporal dimension of these biological processes are definitely needed. Lung cancer and colorectal cancer (CRC) are the leading causes of cancer death worldwide. Here, by developing an innovative computational approach we named Phase-Finder, we show that the molecular heterogeneity characterizing the two deadliest cancers, CRC and lung adenocarcinoma (LUAD), rather than a merely stochastic event is the readout of specific cancer molecular states which correlate with time-qualified patterns of gene expression. We performed time-course transcriptome analysis of CRC and LUAD cell lines and upon computing circadian genes expression-based correlation matrices we derived pseudo-time points to infer time-qualified patterns in the transcriptomic analysis of real-world data (RWD) from large cohorts of CRC and LUAD patients. Our temporal classification of CRC and LUAD cohorts was able to effectively render time-specific patterns in cancer phenotype switching determining dynamical distribution of molecular subtypes impacting patient prognosis.

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.

Spatial transcriptomics delineates molecular features and cellular plasticity in lung adenocarcinoma progression

Indolent (lepidic) and aggressive (micropapillary, solid, and poorly differentiated acinar) histologic subtypes often coexist within a tumor tissue of lung adenocarcinoma (LUAD), but the molecular features associated with different subtypes and their transitions remain elusive. Here, we combine spatial transcriptomics and multiplex immunohistochemistry to elucidate molecular characteristics and cellular plasticity of distinct histologic subtypes of LUAD. We delineate transcriptional reprogramming and dynamic cell signaling that determine subtype progression, especially hypoxia-induced regulatory network. Different histologic subtypes exhibit heterogeneity in dedifferentiation states. Additionally, our results show that macrophages are the most abundant cell type in LUAD, and identify different tumor-associated macrophage subpopulations that are unique to each histologic subtype, which might contribute to an immunosuppressive microenvironment. Our results provide a systematic landscape of molecular profiles that drive LUAD subtype progression, and demonstrate potentially novel therapeutic strategies and targets for invasive lung adenocarcinoma.

A novel necroptosis related gene signature and regulatory network for overall survival prediction in lung adenocarcinoma

We downloaded the mRNA expression profiles of patients with LUAD and corresponding clinical data from The Cancer Genome Atlas (TCGA) database and used the Least Absolute Shrinkage and Selection Operator Cox regression model to construct a multigene signature in the TCGA cohort, which was validated with patient data from the GEO cohort. Results showed differences in the expression levels of 120 necroptosis-related genes between normal and tumor tissues. An eight-gene signature (CYLD, FADD, H2AX, RBCK1, PPIA, PPID, VDAC1, and VDAC2) was constructed through univariate Cox regression, and patients were divided into two risk groups. The overall survival of patients in the high-risk group was significantly lower than of the patients in the low-risk group in the TCGA and GEO cohorts, indicating that the signature has a good predictive effect. The time-ROC curves revealed that the signature had a reliable predictive role in both the TCGA and GEO cohorts. Enrichment analysis showed that differential genes in the risk subgroups were associated with tumor immunity and antitumor drug sensitivity. We then constructed an mRNA-miRNA-lncRNA regulatory network, which identified lncRNA AL590666. 2/let-7c-5p/PPIA as a regulatory axis for LUAD. Real-time quantitative PCR (RT-qPCR) was used to validate the expression of the 8-gene signature. In conclusion, necroptosis-related genes are important factors for predicting the prognosis of LUAD and potential therapeutic targets.

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.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

Single-cell data analysis of malignant epithelial cell heterogeneity in lung adenocarcinoma for patient classification and prognosis prediction

Lung cancer is one of the leading causes of cancer-related death. Most advanced lung adenocarcinoma (LUAD) patients have poor survival because of drug resistance and relapse. Neglecting intratumoral heterogeneity might be one of the reasons for treatment insensitivity, while single-cell RNA sequencing (scRNA-seq) technologies can provide transcriptome information at the single-cell level. Herein, we combined scRNA-seq and bulk RNA-seq data of LUAD and identified a novel cluster of malignant epithelial cells - KRT81+ malignant epithelial cells - associated with worse prognoses. Further analysis revealed that the hypoxia and EMT pathways of these cells were activated to predispose them to differentiate into metastatic lung adenocarcinoma cells. Finally, we also studied the role of these tumor cells in the immune microenvironment and their role in the classification and prognosis prediction of lung adenocarcinoma patients.

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.