Clonal Heterogeneity and Tumor Evolution: Past, Present, and the Future

From genetic mosaicism to tumorigenesis through indirect genetic effects

Genetic mosaicism has long been linked to aging, and several hypotheses have been proposed to explain the potential connections between mosaicism and susceptibility to cancer. It has been proposed that mosaicism may disrupt tissue homeostasis by affecting intercellular communications and releasing microenvironmental constraints within tissues. The underlying mechanisms driving these tissue-level influences remain unidentified, however. Here, we present an evolutionary perspective on the interplay between mosaicism and cancer, suggesting that the tissue-level impacts of genetic mosaicism can be attributed to Indirect Genetic Effects (IGEs). IGEs can increase the level of cellular stochasticity and phenotypic instability among adjacent cells, thereby elevating the risk of cancer development within the tissue. Moreover, as cells experience phenotypic changes in response to challenging microenvironmental conditions, these changes can initiate a cascade of nongenetic alterations, referred to as Indirect non-Genetic Effects (InGEs), which in turn catalyze IGEs among surrounding cells. We argue that incorporating both InGEs and IGEs into our understanding of the process of oncogenic transformation could trigger a major paradigm shift in cancer research with far-reaching implications for practical applications.

FPR3 reprograms glycolytic metabolism and stemness in gastric cancer via calcium-NFATc1 pathway

Aerobic glycolysis accelerates tumor proliferation and progression, and inhibitors or drugs targeting abnormal cancer metabolism have been developing. Cancer stem-like cells (CSCs) significantly contribute to tumor initiation, metastasis, therapy resistance, and recurrence. Formyl peptide receptor 3 (FPR3), a member of FPR family, involves in inflammation, tissue repair, and angiogenesis. However, studies in exploring the regulatory mechanisms of aerobic glycolysis and CSCs by FPR3 in gastric cancer (GC) remain unknown. Here, we demonstrated that overexpressed FPR3 suppressed glycolytic capacity and stemness of tumor cells, then inhibited GC cells proliferation. Mechanistically, FPR3 impeded cytoplasmic calcium ion flux and hindered nuclear factor of activated T cells 1 (NFATc1) nuclear translocation, leading to the transcriptional inactivation of NFATc1-binding neurogenic locus notch homolog protein 3 (NOTCH3) promoter, subsequently obstructing NOTCH3 expression and the AKT/mTORC1 signaling pathway, and ultimately downregulating glycolysis. Additionally, NFATc1 directly binds to the sex determining region Y-box 2 (SOX2) promoter and modifies stemness in GC. In conclusion, our work illustrated that FPR3 played a negative role in GC progression by modulating NFATc1-mediated glycolysis and stemness in a calcium-dependent manner, providing potential insights into cancer therapy.

Beyond typical histology of BAP1-inactivated melanocytoma

BAP1-inactivated melanocytoma (BIM) is a novel subgroup of melanocytic neoplasm listed in the 5th edition of WHO classification of skin tumor. BIM is characterized by two molecular alterations, including a mitogenic driver mutation (usually BRAF gene) and the loss of function of BAP1, a tumor suppressor gene located on chromosome 3p21, which encodes for BRCA1-associated protein (BAP1). The latter represents a nuclear-localized deubiquitinase involved in several cellular processes including cell cycle regulation, chromatin remodeling, DNA damage response, differentiation, senescence and cell death. BIMs are histologically characterized by a population of large epithelioid melanocytes with well-demarcated cytoplasmic borders and copious eosinophilic cytoplasm, demonstrating loss of BAP1 nuclear expression by immunohistochemistry. Recently, we have published a series of 50 cases, extending the morphological spectrum of the neoplasm and highlighting some new microscopic features. In the current article, we focus on some new histological features, attempting to explain and link them to certain mechanisms of tumor development, including senescence, endoreplication, endocycling, asymmetric cytokinesis, entosis and others. In light of the morphological and molecular findings observed in BIM, we postulated that this entity unmasks a fine mechanism of tumor in which both clonal/stochastic and hierarchical model can be unified.

GWAS advancements to investigate disease associations and biological mechanisms

Genome-wide association studies (GWAS) have been instrumental in elucidating the genetic architecture of various traits and diseases. Despite the success of GWAS, inherent limitations such as identifying rare and ultra-rare variants, the potential for spurious associations, and in pinpointing causative agents can undermine diagnostic capabilities. This review provides an overview of GWAS and highlights recent advances in genetics that employ a range of methodologies, including Whole Genome Sequencing (WGS), Mendelian Randomization (MR), the Pangenome's high-quality T2T-CHM13 panel, and the Human BioMolecular Atlas Program (HuBMAP), as potential enablers of current and future GWAS research. State of the literature demonstrate the capabilities of these techniques in enhancing the statistical power of GWAS. WGS, with its comprehensive approach, captures the entire genome, surpassing the capabilities of the traditional GWAS technique focused on predefined Single Nucleotide Polymorphism (SNP) sites. The Pangenome's T2T-CHM13 panel, with its holistic approach, aids in the analysis of regions with high sequence identity, such as segmental duplications (SDs). Mendelian Randomization has advanced causative inference, improving clinical diagnostics and facilitating definitive conclusions. Furthermore, spatial biology techniques like HuBMAP, enable 3D molecular mapping of tissues at single-cell resolution, offering insights into pathology of complex traits. This study aims to elucidate and advocate for the increased application of these technologies, highlighting their potential to shape the future of GWAS research.

Diffusion kurtosis imaging-based habitat analysis identifies high-risk molecular subtypes and heterogeneity matching in diffuse gliomas

OBJECTIVE

High-risk types of diffuse gliomas in adults include isocitrate dehydrogenase (IDH) wild-type glioblastomas and grade 4 astrocytomas. Achieving noninvasive prediction of high-risk molecular subtypes of gliomas is important for personalized and precise diagnosis and treatment.

METHODS

We retrospectively collected data from 116 patients diagnosed with adult diffuse gliomas. Multiple high-risk molecular markers were tested, and various habitat models and whole-tumor models were constructed based on preoperative routine and diffusion kurtosis imaging (DKI) sequences to predict high-risk molecular subtypes of gliomas. Feature selection and model construction utilized Least absolute shrinkage and selection operator (LASSO) and support vector machine (SVM). Finally, the Wilcoxon rank-sum test was employed to explore the correlation between habitat quantitative features (intra-tumor heterogeneity score，ITH score) and heterogeneity, as well as high-risk molecular subtypes.

RESULTS

The results showed that the habitat analysis model based on DKI performed remarkably well (with AUC values reaching 0.977 and 0.902 in the training and test sets, respectively). The model's performance was further enhanced when combined with clinical variables. (The AUC values were 0.994 and 0.920, respectively.) Additionally, we found a close correlation between ITH score and heterogeneity, with statistically significant differences observed between high-risk and non-high-risk molecular subtypes.

INTERPRETATION

The habitat model based on DKI is an ideal means for preoperatively predicting high-risk molecular subtypes of gliomas, holding significant value for noninvasively alerting malignant gliomas and those with malignant transformation potential.

Microenvironment shapes small-cell lung cancer neuroendocrine states and presents therapeutic opportunities

Small-cell lung cancer (SCLC) is the most fatal form of lung cancer. Intratumoral heterogeneity, marked by neuroendocrine (NE) and non-neuroendocrine (non-NE) cell states, defines SCLC, but the cell-extrinsic drivers of SCLC plasticity are poorly understood. To map the landscape of SCLC tumor microenvironment (TME), we apply spatially resolved transcriptomics and quantitative mass spectrometry-based proteomics to metastatic SCLC tumors obtained via rapid autopsy. The phenotype and overall composition of non-malignant cells in the TME exhibit substantial variability, closely mirroring the tumor phenotype, suggesting TME-driven reprogramming of NE cell states. We identify cancer-associated fibroblasts (CAFs) as a crucial element of SCLC TME heterogeneity, contributing to immune exclusion, and predicting exceptionally poor prognosis. Our work provides a comprehensive map of SCLC tumor and TME ecosystems, emphasizing their pivotal role in SCLC's adaptable nature, opening possibilities for reprogramming the TME-tumor communications that shape SCLC tumor states.

ScRNA-seq revealed the tumor microenvironment heterogeneity related to the occurrence and metastasis in upper urinary tract urothelial carcinoma

Metastasis is the greatest clinical challenge for UTUCs, which may have distinct molecular and cellular characteristics from earlier cancers. Herein, we provide single-cell transcriptome profiles of UTUC para cancer normal tissue, primary tumor lesions, and lymphatic metastases to explore possible mechanisms associated with UTUC occurrence and metastasis. From 28,315 cells obtained from normal and tumor tissues of 3 high-grade UTUC patients, we revealed the origin of UTUC tumor cells and the homology between metastatic and primary tumor cells. Unlike the immunomicroenvironment suppression of other tumors, we found no immunosuppression in the tumor microenvironment of UTUC. Moreover, it is imperative to note that stromal cells are pivotal in the advancement of UTUC. This comprehensive single-cell exploration enhances our comprehension of the molecular and cellular dynamics of metastatic UTUCs and discloses promising diagnostic and therapeutic targets in cancer-microenvironment interactions.

Integration of Pan-Cancer Cell Line and Single-Cell Transcriptomic Profiles Enables Inference of Therapeutic Vulnerabilities in Heterogeneous Tumors

Single-cell RNA sequencing (scRNA-seq) greatly advanced the understanding of intratumoral heterogeneity by identifying distinct cancer cell subpopulations. However, translating biological differences into treatment strategies is challenging due to a lack of tools to facilitate efficient drug discovery that tackles heterogeneous tumors. Developing such approaches requires accurate prediction of drug response at the single-cell level to offer therapeutic options to specific cell subpopulations. Here, we developed a transparent computational framework (nicknamed scIDUC) to predict therapeutic efficacies on an individual cell basis by integrating single-cell transcriptomic profiles with large, data-rich pan-cancer cell line screening data sets. This method achieved high accuracy in separating cells into their correct cellular drug response statuses. In three distinct prospective tests covering different diseases (rhabdomyosarcoma, pancreatic ductal adenocarcinoma, and castration-resistant prostate cancer), the predicted results using scIDUC were accurate and mirrored biological expectations. In the first two tests, the framework identified drugs for cell subpopulations that were resistant to standard-of-care (SOC) therapies due to intrinsic resistance or tumor microenvironmental effects, and the results showed high consistency with experimental findings from the original studies. In the third test using newly generated SOC therapy-resistant cell lines, scIDUC identified efficacious drugs for the resistant line, and the predictions were validated with in vitro experiments. Together, this study demonstrates the potential of scIDUC to quickly translate scRNA-seq data into drug responses for individual cells, displaying the potential as a tool to improve the treatment of heterogenous tumors.

SIGNIFICANCE

A versatile method that infers cell-level drug response in scRNA-seq data facilitates the development of therapeutic strategies to target heterogeneous subpopulations within a tumor and address issues such as treatment failure and resistance.

Single-cell analysis of tumor microenvironment and cell adhesion reveals that interleukin-1 beta promotes cancer cell proliferation in breast cancer

BACKGROUND

Triple-negative breast cancer (TNBC), which is so called because of the lack of estrogen receptors (ER), progesterone receptors (PR), and human epidermal growth factor receptor 2 (HER2) receptors on the cancer cells, accounts for 10%-15% of all breast cancers. The heterogeneity of the tumor microenvironment is high. However, the role of plasma cells controlling the tumor migration progression in TNBC is still not fully understood.

METHODS

We analyzed single-cell RNA sequencing data from five HER2 positive, 12 ER positive/PR positive, and nine TNBC samples. The potential targets were validated by immunohistochemistry.

RESULTS

Plasma cells were enriched in TNBC samples, which was consistent with validation using data from The Cancer Genome Atlas. Cell communication analysis revealed that plasma cells interact with T cells through the intercellular adhesion molecule 2-integrin-aLb2 complex, and then release interleukin 1 beta (IL1B), as verified by immunohistochemistry, ultimately promoting tumor growth.

CONCLUSION

Our results revealed the role of plasma cells in TNBC and identified IL1B as a new prognostic marker for TNBC.

Whole-tumoral metabolic heterogeneity in 18F-FDG PET/CT is a novel prognostic marker for neuroblastoma

BACKGROUND

Neuroblastoma (NB) is a highly heterogeneous tumor, and more than half of newly diagnosed NB are associated with extensive metastases. Accurately characterizing the heterogeneity of whole-body tumor lesions remains clinical challenge. This study aims to quantify whole-tumoral metabolic heterogeneity (WMH) derived from whole-body tumor lesions, and investigate the prognostic value of WMH in NB.

METHODS

We retrospectively enrolled 95 newly diagnosed pediatric NB patients in our department. Traditional semi-quantitative PET/CT parameters including the maximum standardized uptake value (SUVmax), the mean standardized uptake value (SUVmean), the peak standardized uptake value (SUVpeak), metabolic tumor volume (MTV) and total lesion glycolysis (TLG) were measured. These PET/CT parameters were expressed as PSUVmax, PSUVmean, PSUVpeak, PMTV, PTLG for primary tumor, WSUVmax, WSUVmean, WSUVpeak, WMTV, WTLG for whole-body tumor lesions. The metabolic heterogeneity was quantified using the areas under the curve of the cumulative SUV-volume histogram index (AUC-CSH index). Intra-tumoral metabolic heterogeneity (IMH) and WMH were extracted from primary tumor and whole-body tumor lesions, respectively. The outcome endpoints were overall survival (OS) and progression-free survival (PFS). Survival analysis was performed utilizing the univariate and multivariate Cox proportional hazards regression. The optimal cut-off values for metabolic parameters were obtained by receiver operating characteristic curve (ROC).

RESULTS

During follow up, 27 (28.4%) patients died, 21 (22.1%) patients relapsed and 47 (49.5%) patients remained progression-free survival, with a median follow-up of 35.0 months. In survival analysis, WMTV and WTLG were independent indicators of PFS, and WMH was an independent risk factor of PFS and OS. However, IMH only showed association with PFS and OS. In addition to metabolic parameters, the International Neuroblastoma Staging System (INSS) was identified as an independent risk factor for PFS, and neuron-specific enolase (NSE) served as an independent predictor of OS.

CONCLUSION

WMH was an independent risk factor for PFS and OS, suggesting its potential as a novel prognostic marker for newly diagnosed NB patients.

Integrative molecular and spatial analysis reveals evolutionary dynamics and tumor-immune interplay of in situ and invasive acral melanoma

In acral melanoma (AM), progression from in situ (AMis) to invasive AM (iAM) leads to significantly reduced survival. However, evolutionary dynamics during this process remain elusive. Here, we report integrative molecular and spatial characterization of 147 AMs using genomics, bulk and single-cell transcriptomics, and spatial transcriptomics and proteomics. Vertical invasion from AMis to iAM displays an early and monoclonal seeding pattern. The subsequent regional expansion of iAM exhibits two distinct patterns, clonal expansion and subclonal diversification. Notably, molecular subtyping reveals an aggressive iAM subset featured with subclonal diversification, increased epithelial-mesenchymal transition (EMT), and spatial enrichment of APOE+/CD163+ macrophages. In vitro and ex vivo experiments further demonstrate that APOE+CD163+ macrophages promote tumor EMT via IGF1-IGF1R interaction. Adnexal involvement can predict AMis with higher invasive potential whereas APOE and CD163 serve as prognostic biomarkers for iAM. Altogether, our results provide implications for the early detection and treatment of AM.

TMO-Net: an explainable pretrained multi-omics model for multi-task learning in oncology

Cancer is a complex disease composing systemic alterations in multiple scales. In this study, we develop the Tumor Multi-Omics pre-trained Network (TMO-Net) that integrates multi-omics pan-cancer datasets for model pre-training, facilitating cross-omics interactions and enabling joint representation learning and incomplete omics inference. This model enhances multi-omics sample representation and empowers various downstream oncology tasks with incomplete multi-omics datasets. By employing interpretable learning, we characterize the contributions of distinct omics features to clinical outcomes. The TMO-Net model serves as a versatile framework for cross-modal multi-omics learning in oncology, paving the way for tumor omics-specific foundation models.

Suppressing mitochondrial inner membrane protein (IMMT) inhibits the proliferation of breast cancer cells through mitochondrial remodeling and metabolic regulation

Metabolic reprogramming is widely recognized as a hallmark of malignant tumors, and the targeting of metabolism has emerged as an appealing approach for cancer treatment. Mitochondria, as pivotal organelles, play a crucial role in the metabolic regulation of tumor cells, and their morphological and functional alterations are intricately linked to the biological characteristics of tumors. As a key regulatory subunit of mitochondria, mitochondrial inner membrane protein (IMMT), plays a vital role in degenerative diseases, but its role in tumor is almost unknown. The objective of this research was to investigate the roles that IMMT play in the development and progression of breast cancer (BC), as well as to elucidate the underlying biological mechanisms that drive these effects. In this study, it was confirmed that the expression of IMMT in BC tissues was significantly higher than that in normal tissues. The analysis of The Cancer Genome Atlas (TCGA) database revealed that IMMT can serve as an independent prognostic factor for BC patients. Additionally, verification in clinical specimens of BC demonstrated a positive association between high IMMT expression and larger tumor size (> 2 cm), Ki-67 expression (> 15%), and HER-2 status. Furthermore, in vitro experiments have substantiated that the suppression of IMMT expression resulted in a reduction in cell proliferation and alterations in mitochondrial cristae, concomitant with the liberation of cytochrome c, but it did not elicit mitochondrial apoptosis. Through Gene Set Enrichment Analysis (GSEA) analysis, we have predicted the associated metabolic genes and discovered that IMMT potentially modulates the advancement of BC through its interaction with 16 metabolic-related genes, and the changes in glycolysis related pathways have been validated in BC cell lines after IMMT inhibition. Consequently, this investigation furnishes compelling evidence supporting the classification of IMMT as prognostic marker in BC, and underscoring its prospective utility as a novel target for metabolic therapy.

Characterization of EpCAM in thyroid cancer biology by three-dimensional spheroids in vitro model

BACKGROUND

Thyroid cancer (TC) is the most common endocrine malignancy. Nowadays, undifferentiated thyroid cancers (UTCs) are still lethal, mostly due to the insurgence of therapy resistance and disease relapse. These events are believed to be caused by a subpopulation of cancer cells with stem-like phenotype and specific tumor-initiating abilities, known as tumor-initiating cells (TICs). A comprehensive understanding of how to isolate and target these cells is necessary. Here we provide insights into the role that the protein Epithelial Cell Adhesion Molecule (EpCAM), a known TICs marker for other solid tumors, may have in TC biology, thus considering EpCAM a potential marker of thyroid TICs in UTCs.

METHODS

The characterization of EpCAM was accomplished through Western Blot and Immunofluorescence on patient-derived tissue samples, adherent cell cultures, and 3D sphere cultures of poorly differentiated thyroid cancer (PDTC) and anaplastic thyroid cancer (ATC) cell lines. The frequency of tumor cells with putative tumor-initiating ability within the 3D cultures was assessed through extreme limiting dilution analysis (ELDA). EpCAM proteolytic cleavages were studied through treatments with different cleavages' inhibitors. To evaluate the involvement of EpCAM in inducing drug resistance, Vemurafenib (PLX-4032) treatments were assessed through MTT assay.

RESULTS

Variable EpCAM expression pattern was observed in TC tissue samples, with increased cleavage in the more UTC. We demonstrated that EpCAM is subjected to an intense cleavage process in ATC-derived 3D tumor spheres and that the 3D model faithfully mimics what was observed in patient's samples. We also proved that the integrity of the protein appears to be crucial for the generation of 3D spheres, and its expression and cleavage in a 3D system could contribute to drug resistance in thyroid TICs.

CONCLUSIONS

Our data provide novel information on the role of EpCAM expression and cleavage in the biology of thyroid TICs, and our 3D model reflects the variability of EpCAM cleavage observed in tissue samples. EpCAM evaluation could play a role in clinical decisions regarding patient therapy since its expression and cleavage may have a fundamental role in the switch to a drug-resistant phenotype of UTC cells.

Subpopulation commensalism promotes Rac1-dependent invasion of single cells via laminin-332

Phenotypic heterogeneity poses a significant hurdle for cancer treatment but is under-characterized in the context of tumor invasion. Amidst the range of phenotypic heterogeneity across solid tumor types, collectively invading cells and single cells have been extensively characterized as independent modes of invasion, but their intercellular interactions have rarely been explored. Here, we isolated collectively invading cells and single cells from the heterogeneous 4T1 cell line and observed extensive transcriptional and epigenetic diversity across these subpopulations. By integrating these datasets, we identified laminin-332 as a protein complex exclusively secreted by collectively invading cells. Live-cell imaging revealed that laminin-332 derived from collectively invading cells increased the velocity and directionality of single cells. Despite collectively invading and single cells having similar expression of the integrin α6β4 dimer, single cells demonstrated higher Rac1 activation upon laminin-332 binding to integrin α6β4. This mechanism suggests a novel commensal relationship between collectively invading and single cells, wherein collectively invading cells promote the invasive potential of single cells through a laminin-332/Rac1 axis.

Strength of selection in lung tumors correlates with clinical features better than tumor mutation burden

Single nucleotide substitutions are the most common type of somatic mutations in cancer genome. The goal of this study was to use publicly available somatic mutation data to quantify negative and positive selection in individual lung tumors and test how strength of directional and absolute selection is associated with clinical features. The analysis found a significant variation in strength of selection (both negative and positive) among tumors, with median selection tending to be negative even though tumors with strong positive selection also exist. Strength of selection estimated as the density of missense mutations relative to the density of silent mutations showed only a weak correlation with tumor mutation burden. In the "all histology together" analysis we found that absolute strength of selection was strongly correlated with all clinically relevant features analyzed. In histology-stratified analysis selection was strongest in small cell lung cancer. Selection in adenocarcinoma was somewhat higher compared to squamous cell carcinoma. The study suggests that somatic mutation- based quantifying of directional and absolute selection in individual tumors can be a useful biomarker of tumor aggressiveness.

Research autopsy programmes in oncology: shared experience from 14 centres across the world

While there is a great clinical need to understand the biology of metastatic cancer in order to treat it more effectively, research is hampered by limited sample availability. Research autopsy programmes can crucially advance the field through synchronous, extensive, and high-volume sample collection. However, it remains an underused strategy in translational research. Via an extensive questionnaire, we collected information on the study design, enrolment strategy, study conduct, sample and data management, and challenges and opportunities of research autopsy programmes in oncology worldwide. Fourteen programmes participated in this study. Eight programmes operated 24 h/7 days, resulting in a lower median postmortem interval (time between death and start of the autopsy, 4 h) compared with those operating during working hours (9 h). Most programmes (n = 10) succeeded in collecting all samples within a median of 12 h after death. A large number of tumour sites were sampled during each autopsy (median 15.5 per patient). The median number of samples collected per patient was 58, including different processing methods for tumour samples but also non-tumour tissues and liquid biopsies. Unique biological insights derived from these samples included metastatic progression, treatment resistance, disease heterogeneity, tumour dormancy, interactions with the tumour micro-environment, and tumour representation in liquid biopsies. Tumour patient-derived xenograft (PDX) or organoid (PDO) models were additionally established, allowing for drug discovery and treatment sensitivity assays. Apart from the opportunities and achievements, we also present the challenges related with postmortem sample collections and strategies to overcome them, based on the shared experience of these 14 programmes. Through this work, we hope to increase the transparency of postmortem tissue donation, to encourage and aid the creation of new programmes, and to foster collaborations on these unique sample collections. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Epigenetic regulation of breast cancer metastasis

Breast cancer is the most frequently diagnosed malignancy and the second leading cause of cancer-related mortality among women worldwide. Recurrent metastasis is associated with poor patient outcomes and poses a significant challenge in breast cancer therapies. Cancer cells adapting to a new tissue microenvironment is the key event in distant metastasis development, where the disseminating tumor cells are likely to acquire genetic and epigenetic alterations during the process of metastatic colonization. Despite several decades of research in this field, the exact mechanisms governing metastasis are not fully understood. However, emerging body of evidence indicates that in addition to genetic changes, epigenetic reprogramming of cancer cells and the metastatic niche are paramount toward successful metastasis. Here, we review and discuss the latest knowledge about the salient attributes of metastasis and epigenetic regulation in breast cancer and crucial research domains that need further investigation.

Chromatin as an old and new anticancer target

Recent genome-wide analyses identified chromatin modifiers as one of the most frequently mutated classes of genes across all cancers. However, chemotherapies developed for cancers involving DNA damage remain the standard of care for chromatin-deranged malignancies. In this review we address this conundrum by establishing the concept of 'chromatin damage': the non-genetic damage to protein-DNA interactions induced by certain small molecules. We highlight anthracyclines, a class of chemotherapeutic agents ubiquitously applied in oncology, as an example of overlooked chromatin-targeting agents. We discuss our current understanding of this phenomenon and explore emerging chromatin-damaging agents as a basis for further studies to maximize their impact in modern cancer treatment.

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.

Development of new techniques and clinical applications of liquid biopsy in lung cancer management

Lung cancer is an exceedingly malignant tumor reported as having the highest morbidity and mortality of any cancer worldwide, thus posing a great threat to global health. Despite the growing demand for precision medicine, current methods for early clinical detection, treatment and prognosis monitoring in lung cancer are hampered by certain bottlenecks. Studies have found that during the formation and development of a tumor, molecular substances carrying tumor-related genetic information can be released into body fluids. Liquid biopsy (LB), a method for detecting these tumor-related markers in body fluids, maybe a way to make progress in these bottlenecks. In recent years, LB technology has undergone rapid advancements. Therefore, this review will provide information on technical updates to LB and its potential clinical applications, evaluate its effectiveness for specific applications, discuss the existing limitations of LB, and present a look forward to possible future clinical applications. Specifically, this paper will introduce technical updates from the prospectives of engineering breakthroughs in the detection of membrane-based LB biomarkers and other improvements in sequencing technology. Additionally, it will summarize the latest applications of liquid biopsy for the early detection, diagnosis, treatment, and prognosis of lung cancer. We will present the interconnectedness of clinical and laboratory issues and the interplay of technology and application in LB today.

In silico study of heterogeneous tumour-derived organoid response to CAR T-cell therapy

Chimeric antigen receptor (CAR) T-cell therapy is a promising immunotherapy for treating cancers. This method consists in modifying the patients' T-cells to directly target antigen-presenting cancer cells. One of the barriers to the development of this type of therapies, is target antigen heterogeneity. It is thought that intratumour heterogeneity is one of the leading determinants of therapeutic resistance and treatment failure. While understanding antigen heterogeneity is important for effective therapeutics, a good therapy strategy could enhance the therapy efficiency. In this work we introduce an agent-based model (ABM), built upon a previous ABM, to rationalise the outcomes of different CAR T-cells therapies strategies over heterogeneous tumour-derived organoids. We found that one dose of CAR T-cell therapy should be expected to reduce the tumour size as well as its growth rate, however it may not be enough to completely eliminate it. Moreover, the amount of free CAR T-cells (i.e. CAR T-cells that did not kill any cancer cell) increases as we increase the dosage, and so does the risk of side effects. We tested different strategies to enhance smaller dosages, such as enhancing the CAR T-cells long-term persistence and multiple dosing. For both approaches an appropriate dosimetry strategy is necessary to produce "effective yet safe" therapeutic results. Moreover, an interesting emergent phenomenon results from the simulations, namely the formation of a shield-like structure of cells with low antigen expression. This shield turns out to protect cells with high antigen expression. Finally we tested a multi-antigen recognition therapy to overcome antigen escape and heterogeneity. Our studies suggest that larger dosages can completely eliminate the organoid, however the multi-antigen recognition increases the risk of side effects. Therefore, an appropriate small dosages dosimetry strategy is necessary to improve the outcomes. Based on our results, it is clear that a proper therapeutic strategy could enhance the therapies outcomes. In that direction, our computational approach provides a framework to model treatment combinations in different scenarios and to explore the characteristics of successful and unsuccessful treatments.

Themis: advancing precision oncology through comprehensive molecular subtyping and optimization

Recent advances in tumor molecular subtyping have revolutionized precision oncology, offering novel avenues for patient-specific treatment strategies. However, a comprehensive and independent comparison of these subtyping methodologies remains unexplored. This study introduces 'Themis' (Tumor HEterogeneity analysis on Molecular subtypIng System), an evaluation platform that encapsulates a few representative tumor molecular subtyping methods, including Stemness, Anoikis, Metabolism, and pathway-based classifications, utilizing 38 test datasets curated from The Cancer Genome Atlas (TCGA) and significant studies. Our self-designed quantitative analysis uncovers the relative strengths, limitations, and applicability of each method in different clinical contexts. Crucially, Themis serves as a vital tool in identifying the most appropriate subtyping methods for specific clinical scenarios. It also guides fine-tuning existing subtyping methods to achieve more accurate phenotype-associated results. To demonstrate the practical utility, we apply Themis to a breast cancer dataset, showcasing its efficacy in selecting the most suitable subtyping methods for personalized medicine in various clinical scenarios. This study bridges a crucial gap in cancer research and lays a foundation for future advancements in individualized cancer therapy and patient management.

Dual targeting negative enrichment strategy for highly sensitive and purity detection of CTCs

Circulating tumor cells (CTCs) have significant clinical value in early tumor detection, dynamic monitoring and immunotherapy. CTC detection stands out as a leading non-invasive approach for tumor diagnostics and therapeutics. However, the high heterogeneity of CTCs and the occurrence of epithelial-mesenchymal transition (EMT) during metastasis pose challenges to methods relying on EpCAM-positive enrichment. To address these limitations, a method based on negative enrichment of CTCs using specific leukocyte targets has been developed. In this study, aiming to overcome the low purity associated with immunomagnetic beads targeting solely the leukocyte common antigen CD45, we introduced CD66b-modified immunomagnetic beads. CD66b, a specific target for neutrophils with abundant residues, was chosen as a complementary approach. The process involved initial collection of nucleated cells from whole blood samples using density gradient centrifugation. Subsequently, magnetically labeled leukocytes were removed by magnetic field, enabling the capture of CTCs with higher sensitivity and purity while retaining their activity. Finally, we selected 20 clinical blood samples from patients with various cancers to validate the effectiveness of this strategy, providing a new generalized tool for the clinical detection of CTCs.

Oncogenic composite mutations can be predicted by co-mutations and their chromosomal location

Genetic heterogeneity in tumors can show a remarkable selectivity when two or more independent genetic events occur in the same gene. This phenomenon, called composite mutation, points toward a selective pressure, which frequently causes therapy resistance to mutation-specific drugs. Since composite mutations have been described to occur in sub-clonal populations, they are not always captured through biopsy sampling. Here, we provide a proof of concept to predict composite mutations to anticipate which patients might be at risk for sub-clonally driven therapy resistance. We found that composite mutations occur in 5% of cancer patients, mostly affecting the PIK3CA, EGFR, BRAF, and KRAS genes, which are common precision medicine targets. Furthermore, we found a strong and significant relationship between the frequencies of composite mutations with commonly co-occurring mutations in a non-composite context. We also found that co-mutations are significantly enriched on the same chromosome. These observations were independently confirmed using cell line data. Finally, we show the feasibility of predicting compositive mutations based on their co-mutations (AUC 0.62, 0.81, 0.82, and 0.91 for EGFR, PIK3CA, KRAS, and BRAF, respectively). This prediction model could help to stratify patients who are at risk of developing therapy resistance-causing mutations.

Nomograms Predict PFS and OS for SCLC Patients After Standardized Treatment: A Real-World Study

Purpose

This study aims to investigate the process of small cell lung cancer (SCLC) patients from achieving optimal efficacy to experiencing disease progression until death. It examines the predictive value of the treatment response on progression free survival (PFS) and overall survival (OS) of SCLC patients.

Patients and Methods

We conducted a retrospective analysis on 136 SCLC patients diagnosed from 1992 to 2018. Important prognostic factors were identified to construct nomogram models. The predictive performance of the models was evaluated using the receiver operating characteristic curves and calibration curves. Survival differences between groups were compared using Kaplan-Meier survival curves. Subsequently, an independent cohort consisting of 106 SCLC patients diagnosed from 2014 to 2021 was used for validation.

Results

We constructed two nomograms to predict first-line PFS (PFS1) and OS of SCLC. The area under the receiver operating characteristic curves for the PFS1 nomogram predicting PFS at 3-, 6-, and 12-months were 0.919 (95% CI: 0.867-0.970), 0.908 (95% CI: 0.860-0.956) and 0.878 (95% CI: 0.798-0.958), and for the OS nomogram predicting OS at 6-, 12-, and 24-months were 0.814 (95% CI: 0.736-0.892), 0.819 (95% CI: 0.749-0.889) and 0.809 (95% CI: 0.678-0.941), indicating those two models with a high discriminative ability. The calibration curves demonstrated the models had a high degree of consistency between predicted and observed values. According to the risk scores, patients were divided into high-risk and low-risk groups, showing a significant difference in survival rate. And these findings were validated in another independent validation cohort.

Conclusion

Based on the patients' treatment response after standardized treatment, we developed and validated two nomogram models to predict PFS1 and OS of SCLC. The models demonstrated good accuracy, reliability and clinical applicability by validating in an independent cohort.

Targeting Tumor Heterogeneity by Breaking a Stem Cell and Epithelial Niche Interaction Loop

Tumor heterogeneity, the presence of multiple distinct subpopulations of cancer cells between patients or among the same tumors, poses a major challenge to current targeted therapies. The way these different subpopulations interact among themselves and the stromal niche environment, and how such interactions affect cancer stem cell behavior has remained largely unknown. Here, it is shown that an FGF-BMP7-INHBA signaling positive feedback loop integrates interactions among different cell populations, including mammary gland stem cells, luminal epithelial and stromal fibroblast niche components not only in organ regeneration but also, with certain modifications, in cancer progression. The reciprocal dependence of basal stem cells and luminal epithelium is based on basal-derived BMP7 and luminal-derived INHBA, which promote their respective expansion, and is regulated by stromal-epithelial FGF signaling. Targeting this interaction loop, for example, by reducing the function of one or more of its components, inhibits organ regeneration and breast cancer progression. The results have profound implications for overcoming drug resistance because of tumor heterogeneity in future targeted therapies.

Distinct molecular profiles drive multifaceted characteristics of colorectal cancer metastatic seeds

Metastasis of primary tumors remains a challenge for early diagnosis and prevention. The cellular properties and molecular drivers of metastatically competent clones within primary tumors remain unclear. Here, we generated 10-16 single cell-derived lines from each of three colorectal cancer (CRC) tumors to identify and characterize metastatic seeds. We found that intrinsic factors conferred clones with distinct metastatic potential and cellular communication capabilities, determining organ-specific metastasis. Poorly differentiated or highly metastatic clones, rather than drug-resistant clones, exhibited poor clinical prognostic impact. Personalized genetic alterations, instead of mutation burden, determined the occurrence of metastatic potential during clonal evolution. Additionally, we developed a gene signature for capturing metastatic potential of primary CRC tumors and demonstrated a strategy for identifying metastatic drivers using isogenic clones with distinct metastatic potential in primary tumors. This study provides insight into the origin and mechanisms of metastasis and will help develop potential anti-metastatic therapeutic targets for CRC patients.

Identification and classification of glioma subtypes based on RNA-binding proteins

BACKGROUND

Glioma is a common and aggressive primary malignant cancer known for its high morbidity, mortality, and recurrence rates. Despite this, treatment options for glioma are currently restricted. The dysregulation of RBPs has been linked to the advancement of several types of cancer, but their precise role in glioma evolution is still not fully understood. This study sought to investigate how RBPs may impact the development and prognosis of glioma, with potential implications for prognosis and therapy.

METHODS

RNA-seq profiles of glioma and corresponding clinical data from the CGGA database were initially collected for analysis. Unsupervised clustering was utilized to identify crucial tumor subtypes in glioma development. Subsequent time-series analysis and MS model were employed to track the progression of these identified subtypes. RBPs playing a significant role in glioma progression were then pinpointed using WGCNA and Lasso Cox regression models. Functional analysis of these key RBP-related genes was conducted through GSEA. Additionally, the CIBERSORT algorithm was utilized to estimate immune infiltrating cells, while the STRING database was consulted to uncover potential mechanisms of the identified biomarkers.

RESULTS

Six tumor subgroups were identified and found to be highly homogeneous within each subgroup. The progression stages of these tumor subgroups were determined using time-series analysis and a MS model. Through WGCNA, Lasso Cox, and multivariate Cox regression analysis, it was confirmed that BCLAF1 is correlated with survival in glioma patients and is closely linked to glioma progression. Functional annotation suggests that BCLAF1 may impact glioma progression by influencing RNA splicing, which in turn affects the cell cycle, Wnt signaling pathway, and other cancer development pathways.

CONCLUSIONS

The study initially identified six subtypes of glioma progression and assessed their malignancy ranking. Furthermore, it was determined that BCLAF1 could serve as an RBP-related prognostic marker, offering significant implications for the clinical diagnosis and personalized treatment of glioma.

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.

Realizing the promise of Project Optimus: Challenges and emerging opportunities for dose optimization in oncology drug development

Project Optimus is a US Food and Drug Administration Oncology Center of Excellence initiative aimed at reforming the dose selection and optimization paradigm in oncology drug development. This project seeks to bring together pharmaceutical companies, international regulatory agencies, academic institutions, patient advocates, and other stakeholders. Although there is much promise in this initiative, there are several challenges that need to be addressed, including multidimensionality of the dose optimization problem in oncology, the heterogeneity of cancer and patients, importance of evaluating long-term tolerability beyond dose-limiting toxicities, and the lack of reliable biomarkers for long-term efficacy. Through the lens of Totality of Evidence and with the mindset of model-informed drug development, we offer insights into dose optimization by building a quantitative knowledge base integrating diverse sources of data and leveraging quantitative modeling tools to build evidence for drug dosage considering exposure, disease biology, efficacy, toxicity, and patient factors. We believe that rational dose optimization can be achieved in oncology drug development, improving patient outcomes by maximizing therapeutic benefit while minimizing toxicity.

The immunosuppressive landscape in tumor microenvironment

Recent advances in cancer immunotherapy, especially immune checkpoint inhibitors (ICIs), have revolutionized the clinical outcome of many cancer patients. Despite the fact that impressive progress has been made in recent decades, the response rate remains unsatisfactory, and many patients do not benefit from ICIs. Herein, we summarized advanced studies and the latest insights on immune inhibitory factors in the tumor microenvironment. Our in-depth discussion and updated landscape of tumor immunosuppressive microenvironment may provide new strategies for reversing tumor immune evasion, enhancing the efficacy of ICIs therapy, and ultimately achieving a better clinical outcome.

Immunologic tumor microenvironment modulators for turning cold tumors hot

Tumors can be classified into distinct immunophenotypes based on the presence and arrangement of cytotoxic immune cells within the tumor microenvironment (TME). Hot tumors, characterized by heightened immune activity and responsiveness to immune checkpoint inhibitors (ICIs), stand in stark contrast to cold tumors, which lack immune infiltration and remain resistant to therapy. To overcome immune evasion mechanisms employed by tumor cells, novel immunologic modulators have emerged, particularly ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) and programmed cell death protein 1/programmed death-ligand 1(PD-1/PD-L1). These agents disrupt inhibitory signals and reactivate the immune system, transforming cold tumors into hot ones and promoting effective antitumor responses. However, challenges persist, including primary resistance to immunotherapy, autoimmune side effects, and tumor response heterogeneity. Addressing these challenges requires innovative strategies, deeper mechanistic insights, and a combination of immune interventions to enhance the effectiveness of immunotherapies. In the landscape of cancer medicine, where immune cold tumors represent a formidable hurdle, understanding the TME and harnessing its potential to reprogram the immune response is paramount. This review sheds light on current advancements and future directions in the quest for more effective and safer cancer treatment strategies, offering hope for patients with immune-resistant tumors.

MOB-CBAM: A dual-channel attention-based deep learning generalizable model for breast cancer molecular subtypes prediction using mammograms

BACKGROUND AND OBJECTIVE

Deep Learning models have emerged as a significant tool in generating efficient solutions for complex problems including cancer detection, as they can analyze large amounts of data with high efficiency and performance. Recent medical studies highlight the significance of molecular subtype detection in breast cancer, aiding the development of personalized treatment plans as different subtypes of cancer respond better to different therapies.

METHODS

In this work, we propose a novel lightweight dual-channel attention-based deep learning model MOB-CBAM that utilizes the backbone of MobileNet-V3 architecture with a Convolutional Block Attention Module to make highly accurate and precise predictions about breast cancer. We used the CMMD mammogram dataset to evaluate the proposed model in our study. Nine distinct data subsets were created from the original dataset to perform coarse and fine-grained predictions, enabling it to identify masses, calcifications, benign, malignant tumors and molecular subtypes of cancer, including Luminal A, Luminal B, HER-2 Positive, and Triple Negative. The pipeline incorporates several image pre-processing techniques, including filtering, enhancement, and normalization, for enhancing the model's generalization ability.

RESULTS

While identifying benign versus malignant tumors, i.e., coarse-grained classification, the MOB-CBAM model produced exceptional results with 99 % accuracy, precision, recall, and F1-score values of 0.99 and MCC of 0.98. In terms of fine-grained classification, the MOB-CBAM model has proven to be highly efficient in accurately identifying mass with (benign/malignant) and calcification with (benign/malignant) classification tasks with an impressive accuracy rate of 98 %. We have also cross-validated the efficiency of the proposed MOB-CBAM deep learning architecture on two datasets: MIAS and CBIS-DDSM. On the MIAS dataset, an accuracy of 97 % was reported for the task of classifying benign, malignant, and normal images, while on the CBIS-DDSM dataset, an accuracy of 98 % was achieved for the classification of mass with either benign or malignant, and calcification with benign and malignant tumors.

CONCLUSION

This study presents lightweight MOB-CBAM, a novel deep learning framework, to address breast cancer diagnosis and subtype prediction. The model's innovative incorporation of the CBAM enhances precise predictions. The extensive evaluation of the CMMD dataset and cross-validation on other datasets affirm the model's efficacy.

RNA expression of 6 genes from metastatic mucosal gastric cancer serves as the global prognostic marker for gastric cancer with functional validation

BACKGROUND

Molecular analysis of advanced tumors can increase tumor heterogeneity and selection bias. We developed a robust prognostic signature for gastric cancer by comparing RNA expression between very rare early gastric cancers invading only mucosal layer (mEGCs) with lymph node metastasis (Npos) and those without metastasis (Nneg).

METHODS

Out of 1003 mEGCs, all Npos were matched to Nneg using propensity scores. Machine learning approach comparing Npos and Nneg was used to develop prognostic signature. The function and robustness of prognostic signature was validated using cell lines and external datasets.

RESULTS

Extensive machine learning with cross-validation identified the prognostic classifier consisting of four overexpressed genes (HDAC5, NPM1, DTX3, and PPP3R1) and two downregulated genes (MED12 and TP53), and enabled us to develop the risk score predicting poor prognosis. Cell lines engineered to high-risk score showed increased invasion, migration, and resistance to 5-FU and Oxaliplatin but maintained sensitivity to an HDAC inhibitor. Mouse models after tail vein injection of cell lines with high-risk score revealed increased metastasis. In three external cohorts, our risk score was identified as the independent prognostic factor for overall and recurrence-free survival.

CONCLUSION

The risk score from the 6-gene classifier can successfully predict the prognosis of gastric cancer.

Pathologic Assessment and Staging of Multiple Non-Small Cell Lung Carcinomas: A Paradigm Shift with the Emerging Role of Molecular Methods

Non-small cell lung carcinomas (NSCLCs) commonly present as 2 or more separate tumors. Biologically, this encompasses 2 distinct processes: separate primary lung carcinomas (SPLCs), representing independently arising tumors, and intrapulmonary metastases (IPMs), representing intrapulmonary spread of a single tumor. The advent of computed tomography imaging has substantially increased the detection of multifocal NSCLCs. The strategies and approaches for distinguishing between SPLCs and IPMs have evolved significantly over the years. Recently, genomic sequencing of somatic mutations has been widely adopted to identify targetable alterations in NSCLC. These molecular techniques have enabled pathologists to reliably discern clonal relationships among multiple NSCLCs in clinical practice. However, a standardized approach to evaluating and staging multiple NSCLCs using molecular methods is still lacking. Here, we reviewed the historical context and provided an update on the growing applications of genomic testing as a clinically relevant benchmark for determining clonal relationships in multiple NSCLCs, a practice we have designated "comparative molecular profiling." We examined the strengths and limitations of the morphology-based distinction of SPLCs vs IPMs and highlighted pivotal clinical and pathologic insights that have emerged from studying multiple NSCLCs using genomic approaches as a gold standard. Lastly, we suggest a practical approach for evaluating multiple NSCLCs in the clinical setting, considering the varying availability of molecular techniques.

Unveiling the Dynamics behind Glioblastoma Multiforme Single-Cell Data Heterogeneity

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

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 .

Modeling tumors as species-rich ecological communities

Many advanced cancers resist therapeutic intervention. This process is fundamentally related to intra-tumor heterogeneity: multiple cell populations, each with different mutational and phenotypic signatures, coexist within a tumor and its metastatic nodes. Like species in an ecosystem, many cancer cell populations are intertwined in a complex network of ecological interactions. Most mathematical models of tumor ecology, however, cannot account for such phenotypic diversity nor are able to predict its consequences. Here we propose that the Generalized Lotka-Volterra model (GLV), a standard tool to describe complex, species-rich ecological communities, provides a suitable framework to describe the ecology of heterogeneous tumors. We develop a GLV model of tumor growth and discuss how its emerging properties, such as outgrowth and multistability, provide a new understanding of the disease. Additionally, we discuss potential extensions of the model and their application to three active areas of cancer research, namely phenotypic plasticity, the cancer-immune interplay and the resistance of metastatic tumors to treatment. Our work outlines a set of questions and a tentative road map for further research in cancer ecology.

Circulating tumour DNA detects somatic variants contributing to spatial and temporal intra-tumoural heterogeneity in head and neck squamous cell carcinoma

Background

As circulating tumour DNA (ctDNA) liquid biopsy analysis is increasingly incorporated into modern oncological practice, establishing the impact of genomic intra-tumoural heterogeneity (ITH) upon data output is paramount. Despite advances in other cancer types the evidence base in head and neck squamous cell carcinoma (HNSCC) remains poor. We sought to investigate the utility of ctDNA to detect ITH in HNSCC.

Methods

In a pilot cohort of 9 treatment-naïve HNSCC patients, DNA from two intra-tumoural sites (core and margin) was whole-exome sequenced. A 9-gene panel was designed to perform targeted sequencing on pre-treatment plasma cell-free DNA and selected post-treatment samples.

Results

Rates of genomic ITH among the 9 patients was high. COSMIC variants from 19 TCGA HNSCC genes demonstrated an 86.9% heterogeneity rate (present in one tumour sub-site only). Across all patients, cell-free DNA (ctDNA) identified 12.9% (range 7.5-19.8%) of tumour-specific variants, of which 55.6% were specific to a single tumour sub-site only. CtDNA identified 79.0% (range: 55.6-90.9%) of high-frequency variants (tumour VAF>5%). Analysis of ctDNA in serial post-treatment blood samples in patients who suffered recurrence demonstrated dynamic changes in both tumour-specific and acquired variants that predicted recurrence ahead of clinical detection.

Conclusion

We demonstrate that a ctDNA liquid biopsy identified spatial genomic ITH in HNSCC and reliably detected high-frequency driver mutations. Serial sampling allowed post-treatment surveillance and early identification of treatment failure.

scGAL: unmask tumor clonal substructure by jointly analyzing independent single-cell copy number and scRNA-seq data

BACKGROUND

Accurately deciphering clonal copy number substructure can provide insights into the evolutionary mechanism of cancer, and clustering single-cell copy number profiles has become an effective means to unmask intra-tumor heterogeneity (ITH). However, copy numbers inferred from single-cell DNA sequencing (scDNA-seq) data are error-prone due to technically confounding factors such as amplification bias and allele-dropout, and this makes it difficult to precisely identify the ITH.

RESULTS

We introduce a hybrid model called scGAL to infer clonal copy number substructure. It combines an autoencoder with a generative adversarial network to jointly analyze independent single-cell copy number profiles and gene expression data from same cell line. Under an adversarial learning framework, scGAL exploits complementary information from gene expression data to relieve the effects of noise in copy number data, and learns latent representations of scDNA-seq cells for accurate inference of the ITH. Evaluation results on three real cancer datasets suggest scGAL is able to accurately infer clonal architecture and surpasses other similar methods. In addition, assessment of scGAL on various simulated datasets demonstrates its high robustness against the changes of data size and distribution. scGAL can be accessed at: https://github.com/zhyu-lab/scgal .

CONCLUSIONS

Joint analysis of independent single-cell copy number and gene expression data from a same cell line can effectively exploit complementary information from individual omics, and thus gives more refined indication of clonal copy number substructure.

A novel multidrug-resistant cell line from a Chinese patient with pancreatic ductal adenocarcinoma

Chemotherapy resistance poses clinical challenges in pancreatic cancer treatment. Developing cell lines resistant to chemotherapy is crucial for investigating drug resistance mechanisms and identifying alternative treatment pathways. The genetic and biological attributes of pancreatic cancer depend on its aetiology, racial demographics and anatomical origin, underscoring the need for models that comprehensively represent these characteristics. Here, we introduce PDAC-X2, a pancreatic cancer cell line derived from Chinese patients. We conducted a comprehensive analysis encompassing the immune phenotype, biology, genetics, molecular characteristics and tumorigenicity of the cell line. PDAC-X2 cells displayed epithelial morphology and expressed cell markers (CK7 and CK19) alongside other markers (E-cadherin, Vimentin, Ki-67, CEA and CA19-9). The population doubling time averaged around 69 h. In vivo, PDAC-X2 cells consistently maintained their tumorigenicity, achieving a 100% tumour formation rate. Characterised by a predominantly tetraploid karyotype, this cell line exhibited a complex genetic markup. Notably, PDAC-X2 cells demonstrated resistance to multiple drugs, including gemcitabine, paclitaxel, 5-fluorouracil and oxaliplatin. In conclusion, PDAC-X2 presents an invaluable preclinical model. Its utility lies in facilitating the study of drug resistance mechanisms and the exploration of alternative therapeutic approaches aimed at enhancing the prognosis of this tumour type.

A patient-specific lung cancer assembloid model with heterogeneous tumor microenvironments

Cancer models play critical roles in basic cancer research and precision medicine. However, current in vitro cancer models are limited by their inability to mimic the three-dimensional architecture and heterogeneous tumor microenvironments (TME) of in vivo tumors. Here, we develop an innovative patient-specific lung cancer assembloid (LCA) model by using droplet microfluidic technology based on a microinjection strategy. This method enables precise manipulation of clinical microsamples and rapid generation of LCAs with good intra-batch consistency in size and cell composition by evenly encapsulating patient tumor-derived TME cells and lung cancer organoids inside microgels. LCAs recapitulate the inter- and intratumoral heterogeneity, TME cellular diversity, and genomic and transcriptomic landscape of their parental tumors. LCA model could reconstruct the functional heterogeneity of cancer-associated fibroblasts and reflect the influence of TME on drug responses compared to cancer organoids. Notably, LCAs accurately replicate the clinical outcomes of patients, suggesting the potential of the LCA model to predict personalized treatments. Collectively, our studies provide a valuable method for precisely fabricating cancer assembloids and a promising LCA model for cancer research and personalized medicine.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.

Adjuvant PD-1 Checkpoint Inhibition in Early Cutaneous Melanoma: Immunological Mode of Action and the Role of Ultraviolet Radiation

Melanoma ranks as the fifth most common solid cancer in adults worldwide and is responsible for a significant proportion of skin-tumor-related deaths. The advent of immune checkpoint inhibition with anti-programmed death protein-1 (PD-1) antibodies has revolutionized the adjuvant treatment of high-risk, completely resected stage III/IV melanoma. However, not all patients benefit equally. Current strategies for improving outcomes involve adjuvant treatment in earlier disease stages (IIB/C) as well as perioperative treatment approaches. Interfering with T-cell exhaustion to counteract cancer immune evasion and the immunogenic nature of melanoma is key for anti-PD-1 effectiveness. Yet, the biological rationale for the efficacy of adjuvant treatment in clinically tumor-free patients remains to be fully elucidated. High-dose intermittent sun exposure (sunburn) is a well-known primary risk factor for melanomagenesis. Also, ultraviolet radiation (UVR)-induced immunosuppression may impair anti-cancer immune surveillance. In this review, we summarize the current knowledge about adjuvant anti-PD-1 blockade, including a characterization of the main cell types most likely responsible for its efficacy. In conclusion, we propose that local and systemic immunosuppression, to some extent UVR-mediated, can be restored by adjuvant anti-PD-1 therapy, consequently boosting anti-melanoma immune surveillance and the elimination of residual melanoma cell clones.

Leveraging cross-source heterogeneity to improve the performance of bulk gene expression deconvolution

A main limitation of bulk transcriptomic technologies is that individual measurements normally contain contributions from multiple cell populations, impeding the identification of cellular heterogeneity within diseased tissues. To extract cellular insights from existing large cohorts of bulk transcriptomic data, we present CSsingle, a novel method designed to accurately deconvolve bulk data into a predefined set of cell types using a scRNA-seq reference. Through comprehensive benchmark evaluations and analyses using diverse real data sets, we reveal the systematic bias inherent in existing methods, stemming from differences in cell size or library size. Our extensive experiments demonstrate that CSsingle exhibits superior accuracy and robustness compared to leading methods, particularly when dealing with bulk mixtures originating from cell types of markedly different cell sizes, as well as when handling bulk and single-cell reference data obtained from diverse sources. Our work provides an efficient and robust methodology for the integrated analysis of bulk and scRNA-seq data, facilitating various biological and clinical studies.

Pharmacogenomic profiling of intra-tumor heterogeneity using a large organoid biobank of liver cancer

Inter- and intra-tumor heterogeneity is a major hurdle in primary liver cancer (PLC) precision therapy. Here, we establish a PLC biobank, consisting of 399 tumor organoids derived from 144 patients, which recapitulates histopathology and genomic landscape of parental tumors, and is reliable for drug sensitivity screening, as evidenced by both in vivo models and patient response. Integrative analysis dissects PLC heterogeneity, regarding genomic/transcriptomic characteristics and sensitivity to seven clinically relevant drugs, as well as clinical associations. Pharmacogenomic analysis identifies and validates multi-gene expression signatures predicting drug response for better patient stratification. Furthermore, we reveal c-Jun as a major mediator of lenvatinib resistance through JNK and β-catenin signaling. A compound (PKUF-01) comprising moieties of lenvatinib and veratramine (c-Jun inhibitor) is synthesized and screened, exhibiting a marked synergistic effect. Together, our study characterizes the landscape of PLC heterogeneity, develops predictive biomarker panels, and identifies a lenvatinib-resistant mechanism for combination therapy.

Morphological diversity of cancer cells predicts prognosis across tumor types

BACKGROUND

Intratumor heterogeneity drives disease progression and treatment resistance, which can lead to poor patient outcomes. Here, we present a computational approach for quantification of cancer cell diversity in routine hematoxylin-eosin-stained histopathology images.

METHODS

We analyzed publicly available digitized whole-slide hematoxylin-eosin images for 2000 patients. Four tumor types were included: lung, head and neck, colon, and rectal cancers, representing major histology subtypes (adenocarcinomas and squamous cell carcinomas). We performed single-cell analysis on hematoxylin-eosin images and trained a deep convolutional autoencoder to automatically learn feature representations of individual cancer nuclei. We then computed features of intranuclear variability and internuclear diversity to quantify tumor heterogeneity. Finally, we used these features to build a machine-learning model to predict patient prognosis.

RESULTS

A total of 68 million cancer cells were segmented and analyzed for nuclear image features. We discovered multiple morphological subtypes of cancer cells (range = 15-20) that co-exist within the same tumor, each with distinct phenotypic characteristics. Moreover, we showed that a higher morphological diversity is associated with chromosome instability and genomic aneuploidy. A machine-learning model based on morphological diversity demonstrated independent prognostic values across tumor types (hazard ratio range = 1.62-3.23, P < .035) in validation cohorts and further improved prognostication when combined with clinical risk factors.

CONCLUSIONS

Our study provides a practical approach for quantifying intratumor heterogeneity based on routine histopathology images. The cancer cell diversity score can be used to refine risk stratification and inform personalized treatment strategies.

Spatial transcriptomics reveals discrete tumour microenvironments and autocrine loops within ovarian cancer subclones

High-grade serous ovarian carcinoma (HGSOC) is genetically unstable and characterised by the presence of subclones with distinct genotypes. Intratumoural heterogeneity is linked to recurrence, chemotherapy resistance, and poor prognosis. Here, we use spatial transcriptomics to identify HGSOC subclones and study their association with infiltrating cell populations. Visium spatial transcriptomics reveals multiple tumour subclones with different copy number alterations present within individual tumour sections. These subclones differentially express various ligands and receptors and are predicted to differentially associate with different stromal and immune cell populations. In one sample, CosMx single molecule imaging reveals subclones differentially associating with immune cell populations, fibroblasts, and endothelial cells. Cell-to-cell communication analysis identifies subclone-specific signalling to stromal and immune cells and multiple subclone-specific autocrine loops. Our study highlights the high degree of subclonal heterogeneity in HGSOC and suggests that subclone-specific ligand and receptor expression patterns likely modulate how HGSOC cells interact with their local microenvironment.

MoS2-Based Sensor Array for Accurate Identification of Cancer Cells with Ensemble-Modified Aptamers

In this work, we present an array-based chemical nose sensor that utilizes a set of ensemble-modified aptamer (EMAmer) probes to sense subtle physicochemical changes on the cell surface for cancer cell identification. The EMAmer probes are engineered by domain-selective incorporation of different types and/or copies of positively charged functional groups into DNA scaffolds, and their differential interactions with cancer cells can be transduced through competitive adsorption of fluorophore-labeled EMAmer probes loaded on MoS2 nanosheets. We demonstrate that this MoS2-EMAmer-based sensor array enables rapid and effective discrimination among six types of cancer cells and their mixtures with a concentration of 104 cells within 60 min, achieving a 94.4% accuracy in identifying blinded unknown cell samples. The established MoS2-EMAmer sensing platform is anticipated to show significant promise in the advancement of cancer diagnostics.