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

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.

Impact of Clonal Heterogeneity in Multiple Myeloma

Multiple myeloma is characterized by a highly heterogeneous disease distribution within the bone marrow-containing skeletal system. In this review, we introduce the molecular mechanisms underlying clonal heterogeneity and the spatio-temporal evolution of myeloma. We discuss the clinical impact of clonal heterogeneity, which is thought to be one of the biggest obstacles to overcome therapy resistance and to achieve cure.

A randomized, non-comparative phase 2 study of neoadjuvant immune-checkpoint blockade in retroperitoneal dedifferentiated liposarcoma and extremity/truncal undifferentiated pleomorphic sarcoma

Based on the demonstrated clinical activity of immune-checkpoint blockade (ICB) in advanced dedifferentiated liposarcoma (DDLPS) and undifferentiated pleomorphic sarcoma (UPS), we conducted a randomized, non-comparative phase 2 trial ( NCT03307616 ) of neoadjuvant nivolumab or nivolumab/ipilimumab in patients with resectable retroperitoneal DDLPS (n = 17) and extremity/truncal UPS (+ concurrent nivolumab/radiation therapy; n = 10). The primary end point of pathologic response (percent hyalinization) was a median of 8.8% in DDLPS and 89% in UPS. Secondary end points were the changes in immune infiltrate, radiographic response, 12- and 24-month relapse-free survival and overall survival. Lower densities of regulatory T cells before treatment were associated with a major pathologic response (hyalinization > 30%). Tumor infiltration by B cells was increased following neoadjuvant treatment and was associated with overall survival in DDLPS. B cell infiltration was associated with higher densities of regulatory T cells before treatment, which was lost upon ICB treatment. Our data demonstrate that neoadjuvant ICB is associated with complex immune changes within the tumor microenvironment in DDLPS and UPS and that neoadjuvant ICB with concurrent radiotherapy has significant efficacy in UPS.

Changes in chromatin accessibility and transcriptional landscape induced by HDAC inhibitors in TP53 mutated patient-derived colon cancer organoids

Here we present the generation and characterization of patient-derived organoids (PDOs) from colorectal cancer patients. PDOs derived from two patients with TP53 mutations were tested with two different HDAC inhibitors (SAHA and NKL54). Cell death induction, transcriptome, and chromatin accessibility changes were analyzed. HDACIs promote the upregulation of low expressed genes and the downregulation of highly expressed genes. A similar differential effect is observed at the level of chromatin accessibility. Only SAHA is a potent inducer of cell death, which is characterized by the upregulation of BH3-only genes BIK and BMF. Up-regulation of BIK is associated with increased accessibility in an intronic region that has enhancer properties. SAHA, but not NKL54, also causes downregulation of BCL2L1 and decreases chromatin accessibility in three distinct regions of the BCL2L1 locus. Both inhibitors upregulate the expression of innate immunity genes and members of the MHC family. In summary, our exploratory study indicates a mechanism of action for SAHA and demonstrate the low efficacy of NKL54 as a single agent for apoptosis induction, using two PDOs. These observations need to be validated in a larger cohort of PDOs.

The impact of rare cancer and early-line treatments on the benefit of comprehensive genome profiling-based precision oncology

BACKGROUND

Comprehensive genome profiling (CGP) serves as a guide for suitable genomically matched therapies for patients with cancer. However, little is known about the impact of the timing and types of cancer on the therapeutic benefit of CGP.

MATERIALS AND METHODS

A single hospital-based pan-cancer prospective study (TOP-GEAR; UMIN000011141) was conducted to examine the benefit of CGP with respect to the timing and types of cancer. Patients with advanced solid tumors (>30 types) who either progressed with or without standard treatments were genotyped using a single CGP test. The subjects were followed up for a median duration of 590 days to examine therapeutic response, using progression-free survival (PFS), PFS ratio, and factors associated with therapeutic response.

RESULTS

Among the 507 patients, 62 (12.2%) received matched therapies with an overall response rate (ORR) of 32.3%. The PFS ratios (≥1.3) were observed in 46.3% (19/41) of the evaluated patients. The proportion of subjects receiving such therapies in the rare cancer cohort was lower than that in the non-rare cancer cohort (9.6% and 17.4%, respectively; P = 0.010). However, ORR of the rare cancer patients was higher than that in the non-rare cancer cohort (43.8% and 20.0%, respectively; P = 0.046). Moreover, ORR of matched therapies in the first or second line after receiving the CGP test was higher than that in the third or later lines (62.5% and 21.7%, respectively; P = 0.003). Rare cancer and early-line treatment were significantly and independently associated with ORR of matched therapies in multivariable analysis (P = 0.017 and 0.004, respectively).

CONCLUSION

Patients with rare cancer preferentially benefited from tumor mutation profiling by increasing the chances of therapeutic response to matched therapies. Early-line treatments after profiling increase the therapeutic benefit, irrespective of tumor types.

Pursuing dynamics of minimal residual leukemic subclones in relapsed and refractory acute myeloid leukemia during conventional therapy

BACKGROUND

Acute myeloid leukemia (AML) is characterized by clonal heterogeneity, leading to frequent relapses and drug resistance despite intensive clinical therapy. Although AML's clonal architecture has been addressed in many studies, practical monitoring of dynamic changes in those subclones during relapse and treatment is still understudied.

METHOD

Fifteen longitudinal bone marrow (BM) samples were collected from three relapsed and refractory (R/R) AML patients. Using droplet digital polymerase chain reaction (ddPCR), the frequencies of patient's leukemic variants were assessed in seven cell populations that were isolated from each BM sample based on cellular phenotypes. By quantifying mutant clones at the diagnosis, remission, and relapse stages, the distribution of AML subclones was sequentially monitored.

RESULTS

Minimal residual (MR) leukemic subclones exhibit heterogeneous distribution among BM cell populations, including mature leukocyte populations. During AML progression, these subclones undergo active phenotypic transitions and repopulate into distinct cell population regardless of normal hematopoiesis hierarchic order. Of these, MR subclones in progenitor populations of patient BM predominantly carry MR leukemic properties, leading to more robust expansion and stubborn persistence than those in mature populations. Moreover, a minor subset of MR leukemic subclones could be sustained at an extremely low frequency without clonal expansion during relapse.

CONCLUSIONS

In this study, we observed treatment persistent MR leukemic subclones and their phenotypic changes during the treatment process of R/R AML patients. This underscores the importance of preemptive inhibition of clonal promiscuity in R/R AML, proposing a practical method for monitoring AML MR subclones.

Selective Manipulation of the Mitochondria Oxidative Stress in Different Cells Using Intelligent Mesoporous Silica Nanoparticles to Activate On-Demand Immunotherapy for Cancer Treatment

Herein, the vitamin K2 (VK2)/maleimide (MA) coloaded mesoporous silica nanoparticles (MSNs), functional molecules including folic acid (FA)/triphenylphosphine (TPP)/tetrapotassium hexacyanoferrate trihydrate (THT), as well as CaCO3 are explored to fabricate a core-shell-corona nanoparticle (VMMFTTC) for on-demand anti-tumor immunotherapy. After application, the tumor-specific acidic environment first decomposed CaCO3 corona, which significantly levitates the pH value of tumor tissue to convert M2 type macrophage to the antitumor M1 type. The resulting VMMFTT would then internalize in both tumor cells and macrophages via FA-assisted endocytosis and free endocytosis, respectively. These distinct processes generate different amount of VMMFTT in above two cells followed by 1) TPP-induced accumulation in the mitochondria, 2) THT-mediated effective capture of various signal ions to cut off signal transmission and further inhibit glutathione (GSH) generation, 3) ions catalyzed reactive oxygen species (ROS) production through Fenton reaction, 4) sustained release of VK2 and MA to further enhance the ROS production and GSH depletion, which caused significant apoptosis of tumor cells and additional M2-to-M1 macrophage polarization via different processes of oxidative stress. Moreover, the primary tumor apoptosis further matures surrounding immature dendritic cells and activates T cells to continuously promote the antitumor immunotherapy.

Self-adjuvant Astragalus polysaccharide-based nanovaccines for enhanced tumor immunotherapy: a novel delivery system candidate for tumor vaccines

The study of tumor nanovaccines (NVs) has gained interest because they specifically recognize and eliminate tumor cells. However, the poor recognition and internalization by dendritic cells (DCs) and insufficient immunogenicity restricted the vaccine efficacy. Herein, we extracted two molecular-weight Astragalus polysaccharides (APS, 12.19 kD; APSHMw, 135.67 kD) from Radix Astragali and made them self-assemble with OVA257-264 directly forming OVA/APS integrated nanocomplexes through the microfluidic method. The nanocomplexes were wrapped with a sheddable calcium phosphate layer to improve stability. APS in the formed nanocomplexes served as drug carriers and immune adjuvants for potent tumor immunotherapy. The optimal APS-NVs were approximately 160 nm with uniform size distribution and could remain stable in physiological saline solution. The FITC-OVA in APS-NVs could be effectively taken up by DCs, and APS-NVs could stimulate the maturation of DCs, improving the antigen cross-presentation efficiency in vitro. The possible mechanism was that APS can induce DC activation via multiple receptors such as dectin-1 and Toll-like receptors 2 and 4. Enhanced accumulation of APS-NVs both in draining and distal lymph nodes were observed following s.c. injection. Smaller APS-NVs could easily access the lymph nodes. Furthermore, APS-NVs could markedly promote antigen delivery efficiency to DCs and activate cytotoxic T cells. In addition, APS-NVs achieve a better antitumor effect in established B16-OVA melanoma tumors compared with the OVA+Alum treatment group. The antitumor mechanism correlated with the increase in cytotoxic T cells in the tumor region. Subsequently, the poor tumor inhibitory effect of APS-NVs on the nude mouse model of melanoma also confirmed the participation of antitumor adaptive immune response induced by NVs. Therefore, this study developed a promising APS-based tumor NV that is an efficient tumor immunotherapy without systemic side effects.

In situ characterization of the tumor microenvironment

The development of new therapies for cancer is underpinned by an increasing need to comprehensively characterize the tumor microenvironment (TME). While traditional approaches have relied on bulk or single-cell approaches, these are limited in their ability to provide cellular context. Deconvolution of the complex TME is fundamental to understanding tumor dynamics and treatment resistance. Spatially resolved characterization of the TME is likely to provide greater insights into the cellular architecture, tumor-immune cell interactions, receptor-ligand interactions, and cell niches. In turn, these aid in dictating the optimal way in which to target each patient's individual cancer. In this review, we discuss a number of cutting-edge in situ spatial profiling methods giving us new insights into tumor biology.

Radiomics in cone-beam breast CT for the prediction of axillary lymph node metastasis in breast cancer: a multi-center multi-device study

OBJECTIVES

To develop a radiomics model in contrast-enhanced cone-beam breast CT (CE-CBBCT) for preoperative prediction of axillary lymph node (ALN) status and metastatic burden of breast cancer.

METHODS

Two hundred and seventy-four patients who underwent CE-CBBCT examination with two scanners between 2012 and 2021 from two institutions were enrolled. The primary tumor was annotated in each patient image, from which 1781 radiomics features were extracted with PyRadiomics. After feature selection, support vector machine models were developed to predict ALN status and metastatic burden. To avoid overfitting on a specific patient subset, 100 randomly stratified splits were made to assign the patients to either training/fine-tuning or test set. Area under the receiver operating characteristic curve (AUC) of these radiomics models was compared to those obtained when training the models only with clinical features and combined clinical-radiomics descriptors. Ground truth was established by histopathology.

RESULTS

One hundred and eighteen patients had ALN metastasis (N + (≥ 1)). Of these, 74 had low burden (N + (1~2)) and 44 high burden (N + (≥ 3)). The remaining 156 patients had none (N0). AUC values across the 100 test repeats in predicting ALN status (N0/N + (≥ 1)) were 0.75 ± 0.05 (0.67~0.93, radiomics model), 0.68 ± 0.07 (0.53~0.85, clinical model), and 0.74 ± 0.05 (0.67~0.88, combined model). For metastatic burden prediction (N + (1~2)/N + (≥ 3)), AUC values were 0.65 ± 0.10 (0.50~0.88, radiomics model), 0.55 ± 0.10 (0.40~0.80, clinical model), and 0.64 ± 0.09 (0.50~0.90, combined model), with all the ranges spanning 0.5. In both cases, the radiomics model was significantly better than the clinical model (both p < 0.01) and comparable with the combined model (p = 0.56 and 0.64).

CONCLUSIONS

Radiomics features of primary tumors could have potential in predicting ALN metastasis in CE-CBBCT imaging.

CLINICAL RELEVANCE STATEMENT

The findings support potential clinical use of radiomics for predicting axillary lymph node metastasis in breast cancer patients and addressing the limited axilla coverage of cone-beam breast CT.

KEY POINTS

• Contrast-enhanced cone-beam breast CT-based radiomics could have potential to predict N0 vs. N + (≥ 1) and, to a limited extent, N + (1~2) vs. N + (≥ 3) from primary tumor, and this could help address the limited axilla coverage, pending future verifications on larger cohorts. • The average AUC of radiomics and combined models was significantly higher than that of clinical models but showed no significant difference between themselves. • Radiomics features descriptive of tumor texture were found informative on axillary lymph node status, highlighting a higher heterogeneity for tumor with positive axillary lymph node.

Personalized anti-tumor drug efficacy prediction based on clinical data

Anti-tumor drug efficacy prediction poses an unprecedented challenge to realizing personalized medicine. This paper proposes to predict personalized anti-tumor drug efficacy based on clinical data. Specifically, we encode the clinical text as numeric vectors featured with hidden topics for patients using Latent Dirichlet Allocation model. Then, to classify patients into two classes, responsive or non-responsive to a drug, drug efficacy predictors are established by machine learning based on the Latent Dirichlet Allocation topic representation. To evaluate the proposed method, we collected and collated clinical records of lung and bowel cancer patients treated with platinum. Experimental results on the data sets show the efficacy and effectiveness of the proposed method, suggesting the potential value of clinical data in cancer precision medicine. We hope that it will promote the research of drug efficacy prediction based on clinical data.

CoT: a transformer-based method for inferring tumor clonal copy number substructure from scDNA-seq data

Single-cell DNA sequencing (scDNA-seq) has been an effective means to unscramble intra-tumor heterogeneity, while joint inference of tumor clones and their respective copy number profiles remains a challenging task due to the noisy nature of scDNA-seq data. We introduce a new bioinformatics method called CoT for deciphering clonal copy number substructure. The backbone of CoT is a Copy number Transformer autoencoder that leverages multi-head attention mechanism to explore correlations between different genomic regions, and thus capture global features to create latent embeddings for the cells. CoT makes it convenient to first infer cell subpopulations based on the learned embeddings, and then estimate single-cell copy numbers through joint analysis of read counts data for the cells belonging to the same cluster. This exploitation of clonal substructure information in copy number analysis helps to alleviate the effect of read counts non-uniformity, and yield robust estimations of the tumor copy numbers. Performance evaluation on synthetic and real datasets showcases that CoT outperforms the state of the arts, and is highly useful for deciphering clonal copy number substructure.

Clonal differences underlie variable responses to sequential and prolonged treatment

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

Impact of Resistance on Therapeutic Design: A Moran Model of Cancer Growth

Resistance of cancers to treatments, such as chemotherapy, largely arise due to cell mutations. These mutations allow cells to resist apoptosis and inevitably lead to recurrence and often progression to more aggressive cancer forms. Sustained-low dose therapies are being considered as an alternative over maximum tolerated dose treatments, whereby a smaller drug dosage is given over a longer period of time. However, understanding the impact that the presence of treatment-resistant clones may have on these new treatment modalities is crucial to validating them as a therapeutic avenue. In this study, a Moran process is used to capture stochastic mutations arising in cancer cells, inferring treatment resistance. The model is used to predict the probability of cancer recurrence given varying treatment modalities. The simulations predict that sustained-low dose therapies would be virtually ineffective for a cancer with a non-negligible probability of developing a sub-clone with resistance tendencies. Furthermore, calibrating the model to in vivo measurements for breast cancer treatment with Herceptin, the model suggests that standard treatment regimens are ineffective in this mouse model. Using a simple Moran model, it is possible to explore the likelihood of treatment success given a non-negligible probability of treatment resistant mutations and suggest more robust therapeutic schedules.

Canopy2: tumor phylogeny inference by bulk DNA and single-cell RNA sequencing

Tumors are comprised of a mixture of distinct cell populations that differ in terms of genetic makeup and function. Such heterogeneity plays a role in the development of drug resistance and the ineffectiveness of targeted cancer therapies. Insight into this complexity can be obtained through the construction of a phylogenetic tree, which illustrates the evolutionary lineage of tumor cells as they acquire mutations over time. We propose Canopy2, a Bayesian framework that uses single nucleotide variants derived from bulk DNA and single-cell RNA sequencing to infer tumor phylogeny and conduct mutational profiling of tumor subpopulations. Canopy2 uses Markov chain Monte Carlo methods to sample from a joint probability distribution involving a mixture of binomial and beta-binomial distributions, specifically chosen to account for the sparsity and stochasticity of the single-cell data. Canopy2 demystifies the sources of zeros in the single-cell data and separates zeros categorized as non-cancerous (cells without mutations), stochastic (mutations not expressed due to bursting), and technical (expressed mutations not picked up by sequencing). Simulations demonstrate that Canopy2 consistently outperforms competing methods and reconstructs the clonal tree with high fidelity, even in situations involving low sequencing depth, poor single-cell yield, and highly-advanced and polyclonal tumors. We further assess the performance of Canopy2 through application to breast cancer and glioblastoma data, benchmarking against existing methods. Canopy2 is an open-source R package available at https://github.com/annweideman/canopy2.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma

Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter- and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of antigens. In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface antigens that could be targeted by antibodies and chimeric antigen receptor-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas) and 9166 normal tissue samples (from the Genotype-Tissue Expression project), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN, which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative antigens. Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

Single-cell transcriptome analysis upon ECM-remodeling meningioma cells

Meningiomas are the most common tumours that primarily arise in the central nervous system, but their intratumoural heterogeneity has not yet been thoroughly studied. We aimed to investigate the transcriptome characteristics and biological properties of ECM-remodeling meningioma cells. Single-cell RNA sequencing (ScRNA-seq) data from meningioma samples were acquired and used for analyses. We conducted comprehensive bioinformatics analyses, including screening for differentially expressed genes (DEGs), Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway and Gene Ontology (GO) term enrichment analyses, Gene Set Enrichment Analysis (GSEA), protein-protein interaction (PPI) analysis, and copy number variation (CNV) analysis on single-cell sequencing data from meningiomas. Eighteen cell types, including six meningioma subtypes, were identified in the data. ECM-remodeling meningioma cells (MGCs) were mainly distributed in brain-tumour interface tissues. KEGG and GO enrichment analyses revealed that 908 DEGs were mainly related to cell adhesion, extracellular matrix organization, and ECM-receptor interaction. GSEA analysis demonstrated that homophilic cell adhesion via plasma membrane adhesion molecules was significantly enriched (NES = 2.375, P < 0.001). CNV analysis suggested that ECM-remodeling MGCs showed considerably lower average CNV scores. ECM-remodeling MGCs predominantly localized at the brain-tumour interface area and adhere stably to the basement membrane with a lower degree of malignancy. This study provides novel insights into the malignancy of meningiomas.

Bulk anda single-cell transcriptome profiling reveals the molecular characteristics of T cell-mediated tumor killing in pancreatic cancer

Background

Despite the potential of immune checkpoint blockade (ICB) as a promising treatment for Pancreatic adenocarcinoma (PAAD), there is still a need to identify specific subgroups of PAAD patients who may benefit more from ICB. T cell-mediated tumor killing (TTK) is the primary concept behind ICB. We explored subtypes according to genes correlated with the sensitivity to TKK and unraveled their underlying associations for PAAD immunotherapies.

Methods

Genes that control the responsiveness of T cell-induced tumor destruction (GSTTK) were examined in PAAD, focusing on their varying expression levels and association with survival results. Moreover, samples with PAAD were separated into two subsets using unsupervised clustering based on GSTTK. Variability was evident in the tumor immune microenvironment, genetic mutation, and response to immunotherapy among different groups. In the end, we developed TRGscore, an innovative scoring system, and investigated its clinical and predictive significance in determining sensitivity to immunotherapy.

Results

Patients with PAAD were categorized into 2 clusters based on the expression of 52 GSTTKs, which showed varying levels and prognostic relevance, revealing unique TTK patterns. Survival outcome, immune cell infiltration, immunotherapy responses, and functional enrichment are also distinguished among the two clusters. Moreover, we found the CATSPER1 gene promotes the progression of PAAD through experiments. In addition, the TRGscore effectively predicted the responses to chemotherapeutics or immunotherapy in patients with PAAD and overall survival.

Conclusions

TTK exerted a vital influence on the tumor immune environment in PAAD. A greater understanding of TIME characteristics was gained through the evaluation of the variations in TTK modes across different tumor types. It highlights variations in the performance of T cells in PAAD and provides direction for improved treatment approaches.

TOP1 and R-loops facilitate transcriptional DSBs at hypertranscribed cancer driver genes

DNA double-stranded breaks (DSBs) pose a significant threat to genomic integrity, and their generation during essential cellular processes like transcription remains poorly understood. In this study, we employ several techniques to map DSBs, R-loops, and topoisomerase 1 cleavage complex (TOP1cc) to comprehensively investigate the interplay between transcription, DSBs, topoisomerase 1 (TOP1), and R-loops. Our findings reveal the presence of DSBs at highly expressed genes enriched with TOP1 and R-loops. Remarkably, transcription-associated DSBs at these loci are significantly reduced upon depletion of R-loops and TOP1, uncovering the pivotal roles of TOP1 and R-loops in transcriptional DSB formation. By elucidating the intricate interplay between TOP1cc trapping, R-loops, and DSBs, our study provides insights into the mechanisms underlying transcription-associated genomic instability. Moreover, we establish a link between transcriptional DSBs and early molecular changes driving cancer development, highlighting the distinct etiology and molecular characteristics of driver mutations compared to passenger mutations.

A novel type-2 innate lymphoid cell-based immunotherapy for cancer

Cell-based cancer immunotherapy has achieved significant advancements, providing a source of hope for cancer patients. Notwithstanding the considerable progress in cell-based immunotherapy, the persistently low response rates and the exorbitant costs associated with their implementation still present a formidable challenge in clinical settings. In the landscape of cell-based cancer immunotherapies, an uncharted territory involves Type 2 innate lymphoid cells (ILC2s) and interleukin-33 (IL-33) which promotes ILC2 functionality, recognized for their inherent ability to enhance immune responses. Recent discoveries regarding their role in actuating cytolytic T lymphocyte responses, including curbing tumor growth rates and hindering metastasis, have added a new dimension to our understanding of the IL-33/ILC2 axis. These recent insights may hold significant promise for ILC2 cell-based immunotherapy. Nevertheless, the prospect of adoptively transferring ILC2s to confer immune protection against tumors has yet to be investigated. The present study addresses this hypothesis, revealing that ILC2s isolated from the lungs of tumor-bearing mice, and tumor infiltrating ILC2s when adoptively transferred after tumor establishment at a ratio of one ILC2 per sixty tumor cells, leads to an influx of tumor infiltrating CD4+ and CD8+ T lymphocytes as well as tumor infiltrating eosinophils resulting in a remarkable reduction in tumor growth. Moreover, we find that post-adoptive transfer of ILC2s, the number of tumor infiltrating ILC2s is inversely proportional to tumor size. Finally, we find corollaries of the IL-33/ILC2 axis enhancing the infiltration of eosinophils in human prostate carcinomas patients' expressing high levels of IL-33 versus those expressing low levels of IL-33. Our results underscore the heightened efficacy of adoptively transferred ILC2s compared to alternative approaches, revealing an approximately one hundred fifty-fold superiority on a cell-per-cell basis over CAR T-cells in the specific targeting and elimination of tumors within the same experimental model. Overall, this study demonstrates the functional significance of ILC2s in cancer immunosurveillance and provides the proof of concept of the potential utility of ILC2 cell-based cancer immunotherapies.

Circulating Oncometabolite 2-hydroxyglutarate as a Potential Biomarker for Isocitrate Dehydrogenase (IDH1/2) Mutant Cholangiocarcinoma

Isocitrate dehydrogenase (IDH) enzymes catalyze the decarboxylation of isocitrate to alpha-ketoglutarate (αKG). IDH1/2 mutations preferentially convert αKG to R-2-hydroxyglutarate (R2HG), resulting in R2HG accumulation in tumor tissues. We investigated circulating 2-hydroxyglutate (2HG) as potential biomarkers for patients with IDH-mutant (IDHmt) cholangiocarcinoma (CCA). R2HG and S-2-hydroxyglutarate (S2HG) levels in blood and tumor tissues were analyzed in a discovery cohort of patients with IDHmt glioma and CCA. Results were validated in cohorts of patients with CCA and clear-cell renal cell carcinoma. The R2HG/S2HG ratio (rRS) was significantly elevated in tumor tissues, but not in blood for patients with IDHmt glioma, while circulating rRS was elevated in patients with IDHmt CCA. There were overlap distributions of circulating R2HG and total 2HG in patients with both IDHmt and wild-type (IDHwt) CCA, while there was minimal overlap in rRS values between patients with IDHmt and IDHwt CCA. Using the rRS cut-off value of 1.5, the sensitivity of rRS was 90% and specificity was 96.8%. Circulating rRS is significantly increased in patients with IDHmt CCA compare with patients with IDHwt CCA. Circulating rRS is a sensitive and specific surrogate biomarker for IDH1/2 mutations in CCA. It can potentially be used as a tool for monitoring IDH-targeted therapy.

Domain generalization enables general cancer cell annotation in single-cell and spatial transcriptomics

Single-cell and spatial transcriptome sequencing, two recently optimized transcriptome sequencing methods, are increasingly used to study cancer and related diseases. Cell annotation, particularly for malignant cell annotation, is essential and crucial for in-depth analyses in these studies. However, current algorithms lack accuracy and generalization, making it difficult to consistently and rapidly infer malignant cells from pan-cancer data. To address this issue, we present Cancer-Finder, a domain generalization-based deep-learning algorithm that can rapidly identify malignant cells in single-cell data with an average accuracy of 95.16%. More importantly, by replacing the single-cell training data with spatial transcriptomic datasets, Cancer-Finder can accurately identify malignant spots on spatial slides. Applying Cancer-Finder to 5 clear cell renal cell carcinoma spatial transcriptomic samples, Cancer-Finder demonstrates a good ability to identify malignant spots and identifies a gene signature consisting of 10 genes that are significantly co-localized and enriched at the tumor-normal interface and have a strong correlation with the prognosis of clear cell renal cell carcinoma patients. In conclusion, Cancer-Finder is an efficient and extensible tool for malignant cell annotation.

Development and verification of a novel immunogenic cell death-related signature for predicting the prognosis and immune infiltration in triple-negative breast cancer

BACKGROUND

Insufficient understanding of the pathogenesis and tumor immunology of triple-negative breast cancer (TNBC) has limited the development of immunotherapy. The importance of tumor microenvironment (TME) in immunotyping, prognostic assessment and immunotherapy efficacy of cancer has been emphasized, however, potential immunogenic cell death (ICD) related genes function in TME of TNBC has been rarely investigated.

AIMS

To initially explore the role and related mechanisms of ICD in TNBC, especially the role played in the TME of TNBC, and to identify different relevant subtypes based on ICD, and then develop an ICD-related risk score to predict each TNBC patient TME status, prognosis and immunotherapy response.

METHODS AND RESULTS

In this study, we identified distinct ICD-related modification patterns based on 158 TNBC cases in the TCGA-TNBC cohort. We then investigated the possible correlation between ICD-related modification patterns and TME cell infiltration characteristics in TNBC. By using univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analysis, we created a risk scoring system (ICD score) to quantifiably evaluate the impact of ICD-related modification patterns in individual TNBC patient. Two different ICD-related modification patterns were found with significant differences in immune infiltration. Lower ICD score was correlated with higher immune infiltration, tumor mutational burden and significantly enriched in immune-related pathways, indicating a strong ability to activate immune response, which might account for relatively favorable prognosis of TNBC patients and could serve as a predictor to select suitable candidates for immunotherapy. We used two independent cohorts, GSE58812 cohort and Metabric cohort to validate prognosis and immunohistochemistry for preliminary in vitro validation.

CONCLUSION

This study evidenced that the ICD-related modification patterns might exert pivotal roles in the immune infiltration landscape of TNBC and ICD score might act as potential predictors of prognostic assessment and immunotherapy response. This research provides unique insights for individualize immune treatment strategies and promising immunotherapy candidates screening.

Advanced insights on tumor-associated macrophages revealed by single-cell RNA sequencing: The intratumor heterogeneity, functional phenotypes, and cellular interactions

Single-cell RNA sequencing (scRNA-seq) is an emerging technology used for cellular transcriptome analysis. The application of scRNA-seq has led to profoundly advanced oncology research, continuously optimizing novel therapeutic strategies. Intratumor heterogeneity extensively consists of all tumor components, contributing to different tumor behaviors and treatment responses. Tumor-associated macrophages (TAMs), the core immune cells linking innate and adaptive immunity, play significant roles in tumor progression and resistance to therapies. Moreover, dynamic changes occur in TAM phenotypes and functions subject to the regulation of the tumor microenvironment. The heterogeneity of TAMs corresponding to the state of the tumor microenvironment has been comprehensively recognized using scRNA-seq. Herein, we reviewed recent research and summarized variations in TAM phenotypes and functions from a developmental perspective to better understand the significance of TAMs in the tumor microenvironment.

Clonal origin and genomic diversity in Lynch syndrome-associated endometrial cancer with multiple synchronous tumors: Identification of the pathogenicity of MLH1 p.L582H

Lynch syndrome-associated endometrial cancer patients often present multiple synchronous tumors and this assessment can affect treatment strategies. We present a case of a 27-year-old woman with tumors in the uterine corpus, cervix, and ovaries who was diagnosed with endometrial cancer and exhibited cervical invasion and ovarian metastasis. Her family history suggested Lynch syndrome, and genetic testing identified a variant of uncertain significance, MLH1 p.L582H. We conducted immunohistochemical staining, microsatellite instability analysis, and Sanger sequencing for Lynch syndrome-associated cancers in three generations of the family and identified consistent MLH1 loss. Whole-exome sequencing for the corpus, cervical, and ovarian tumors of the proband identified a copy-neutral loss of heterozygosity (LOH) occurring at the MLH1 position in all tumors. This indicated that the germline variant and the copy-neutral LOH led to biallelic loss of MLH1 and was the cause of cancer initiation. All tumors shared a portion of somatic mutations with high mutant allele frequencies, suggesting a common clonal origin. There were no mutations shared only between the cervix and ovary samples. The profiles of mutant allele frequencies shared between the corpus and cervix or ovary indicated that two different subclones originating from the corpus independently metastasized to the cervix or ovary. Additionally, all tumors presented unique mutations in endometrial cancer-associated genes such as ARID1A and PIK3CA. In conclusion, we demonstrated clonal origin and genomic diversity in a Lynch syndrome-associated endometrial cancer, suggesting the importance of evaluating multiple sites in Lynch syndrome patients with synchronous tumors.

Chronic arsenic exposure induces malignant transformation of human HaCaT cells through both deterministic and stochastic changes in transcriptome expression

Biological processes are inherently stochastic, i.e., are partially driven by hard to predict random probabilistic processes. Carcinogenesis is driven both by stochastic and deterministic (predictable non-random) changes. However, very few studies systematically examine the contribution of stochastic events leading to cancer development. In differential gene expression studies, the established data analysis paradigms incentivize expression changes that are uniformly different across the experimental versus control groups, introducing preferential inclusion of deterministic changes at the expense of stochastic processes that might also play a crucial role in the process of carcinogenesis. In this study, we applied simple computational techniques to quantify: (i) The impact of chronic arsenic (iAs) exposure as well as passaging time on stochastic gene expression and (ii) Which genes were expressed deterministically and which were expressed stochastically at each of the three stages of cancer development. Using biological coefficient of variation as an empirical measure of stochasticity we demonstrate that chronic iAs exposure consistently suppressed passaging related stochastic gene expression at multiple time points tested, selecting for a homogenous cell population that undergo transformation. Employing multiple balanced removal of outlier data, we show that chronic iAs exposure induced deterministic and stochastic changes in the expression of unique set of genes, that populate largely unique biological pathways. Together, our data unequivocally demonstrate that both deterministic and stochastic changes in transcriptome-wide expression are critical in driving biological processes, pathways and networks towards clonal selection, carcinogenesis, and tumor heterogeneity.

Cancer cell plasticity, stem cell factors, and therapy resistance: how are they linked?

Cellular plasticity can occur naturally in an organism and is considered an adapting mechanism during the developmental stage. However, abnormal cellular plasticity is observed in different diseased conditions, including cancer. Cancer cell plasticity triggers the stimuli of epithelial-mesenchymal transition (EMT), abnormal epigenetic changes, expression of stem cell factors and implicated signaling pathways, etc., and helps in the maintenance of CSC phenotype. Conversely, CSC maintains the cancer cell plasticity, EMT, and epigenetic plasticity. EMT contributes to increased cell migration and greater diversity within tumors, while epigenetic changes, stem cell factors (OCT4, NANOG, and SOX2), and various signaling pathways allow cancer cells to maintain various phenotypes, giving rise to intra- and inter-tumoral heterogeneity. The intricate relationships between cancer cell plasticity and stem cell factors help the tumor cells adopt drug-tolerant states, evade senescence, and successfully acquire drug resistance with treatment dismissal. Inhibiting molecules/signaling pathways involved in promoting CSCs, cellular plasticity, EMT, and epigenetic plasticity might be helpful for successful cancer therapy management. This review discussed the role of cellular plasticity, EMT, and stem cell factors in tumor initiation, progression, reprogramming, and therapy resistance. Finally, we discussed how the intervention in this axis will help better manage cancers and improve patient survivability.

[Precision medicine in oncology]

Personalized oncology according to current practice is primarily based on tumor biology, which is translated into genomic biomarkers. Mutations in oncogenes and tumor suppressor genes are targeted by rationally designed drugs and, conversely, are used to inform tailored treatment strategies. Faster and cheaper technologies for DNA sequencing enable genomic medicine in a clinical routine setting. Genomic features, tumor biology and clinical implications are integrated into individual therapy recommendations by molecular tumor boards, which have been established at many cancer centers in Germany and worldwide throughout recent years. This article discusses the promises and limitations of genomics-centered precision oncology and highlights future avenues and alternative approaches to individualize cancer treatment.

Evolutionary trajectories of small cell lung cancer under therapy

The evolutionary processes that underlie the marked sensitivity of small cell lung cancer (SCLC) to chemotherapy and rapid relapse are unknown1-3. Here we determined tumour phylogenies at diagnosis and throughout chemotherapy and immunotherapy by multiregion sequencing of 160 tumours from 65 patients. Treatment-naive SCLC exhibited clonal homogeneity at distinct tumour sites, whereas first-line platinum-based chemotherapy led to a burst in genomic intratumour heterogeneity and spatial clonal diversity. We observed branched evolution and a shift to ancestral clones underlying tumour relapse. Effective radio- or immunotherapy induced a re-expansion of founder clones with acquired genomic damage from first-line chemotherapy. Whereas TP53 and RB1 alterations were exclusively part of the common ancestor, MYC family amplifications were frequently not constituents of the founder clone. At relapse, emerging subclonal mutations affected key genes associated with SCLC biology, and tumours harbouring clonal CREBBP/EP300 alterations underwent genome duplications. Gene-damaging TP53 alterations and co-alterations of TP53 missense mutations with TP73, CREBBP/EP300 or FMN2 were significantly associated with shorter disease relapse following chemotherapy. In summary, we uncover key processes of the genomic evolution of SCLC under therapy, identify the common ancestor as the source of clonal diversity at relapse and show central genomic patterns associated with sensitivity and resistance to chemotherapy.

Prostate cancer: Molecular aspects, consequences, and opportunities of the multifocal nature

Prostate cancer is unique compared to other major cancers due to the presence of multiple primary malignant foci in the majority of patients at the time of diagnosis. Each malignant focus has distinct somatic mutations and gene expression patterns, which represents a challenge for the development of prognostic tests for localized prostate cancer. Additionally, the molecular heterogeneity of advanced prostate cancer has important implications for management, particularly for patients with metastatic and locally recurrent cancer. Studies have shown that prostate cancers with mutations in DNA damage response genes are more sensitive to drugs inhibiting the poly ADP-ribose polymerase (PARP) enzyme. However, testing for such mutations should consider both spatial and temporal heterogeneity. Here, we summarize studies where multiregional genomics and transcriptomics analyses have been performed for primary prostate cancer. We further discuss the vast interfocal heterogeneity and how prognostic biomarkers and a molecular definition of the index tumor should be developed. The concept of focal treatments in prostate cancer has been evolving as a demand from patients and clinicians and is one example where there is a need for defining an index tumor. Here, biomarkers must have proven value for individual malignant foci. The potential discovery and implementation of biomarkers that are agnostic to heterogeneity are also explored as an alternative to multisample testing. Thus, deciding upon whole-organ treatment, such as radical prostatectomy, should depend on information from biomarkers which are informative for the whole organ.

Molecular heterogeneity in histomorphologic subtypes of lung adeno carcinoma represents a challenge for treatment decision

Lung cancer is the leading cause in cancer related death, with non-small cell lung cancer (NSCLC) being the most frequent subtype. The importance of NSCLC is reflected by the various targeted therapy options especially for NSCLC adenocarcinomas (lung adeno carcinoma (LUAD)) as well as a set of options for immune therapies. However, despite these therapy advances, the majority of patients do not show a long-term response to either targeted therapy or immune checkpoint inhibition. One reason for treatment failure appears to be the NSCLC tumor heterogeneity. NSCLC heterogeneity might lead to an insufficient molecular characterization of a given sample due to the limited tumor material used for pathological assessment as the majority of analyses is performed on small biopsies. To get a more detailed insight into the tumor heterogeneity of NSCLC LUAD, especially in the light of its different histomorphological growth patterns, we analysed isolated NSCLC growth pattern areas and the corresponding entire tumor samples of a cohort of 31 NSLCS LUAD patients and compared their mutational landscape and their expression profiles. While significant differences of complex biomarkers, like tumor mutational burden (TMB) or microsatellite instability (MSI), were not detected between the five growth patterns -lepidic, papillary, micropapillary, acinar, and solid- we observed various subclonal mutations and copy number variants. Moreover, RNASeq analysis revealed growth pattern specific expression profiles affecting cellular processes like apoptosis, metastasis and proliferation. Taken together, our data provide novel insights into the tumor heterogeneity of LUAD required to overcome tumor heterogeneity related therapy resistance.

Tumour-infiltrating lymphocyte therapy for patients with advanced-stage melanoma

Immunotherapy with immune-checkpoint inhibitors (ICIs) and targeted therapy with BRAF and MEK inhibitors have revolutionized the treatment of melanoma over the past decade. Despite these breakthroughs, the 5-year survival rate of patients with advanced-stage melanoma is at most 50%, emphasizing the need for additional therapeutic strategies. Adoptive cell therapy with tumour-infiltrating lymphocytes (TILs) is a therapeutic modality that has, in the past few years, demonstrated long-term clinical benefit in phase II/III trials involving patients with advanced-stage melanoma, including those with disease progression on ICIs and/or BRAF/MEK inhibitors. In this Review, we summarize the current status of TIL therapies for patients with advanced-stage melanoma, including potential upcoming marketing authorization, the characteristics of TIL therapy products, as well as future strategies that are expected to increase the efficacy of this promising cellular immunotherapy.

Engineering Logic DNA Nanoprobes on Live Cell Membranes for Simultaneously Monitoring Extracellular pH and Precise Drug Delivery

It remains a challenge to use a single probe to simultaneously detect extracellular pH fluctuations and specifically recognize cancer cells for precise drug delivery. Here, we engineered a tetrahedral framework nucleic acid-based logic nanoprobe (isgc8-tFNA) on live cell membranes for simultaneously monitoring extracellular pH and targeted drug delivery. Isgc8-tFNA was anchored stably on the cell surface through three cholesterol molecules inserting into the bilayer of the cell membrane. Once responding to the acidic tumor microenvironment, isgc8-tFNA formed an i-motif structure, leading to turn-on FRET signals for monitoring changes of extracellular pH. The nanoprobe exhibited a narrow pH-response window and excellent reversibility. Moreover, the nanoprobe could execute logic identification on the cell surface for precise drug delivery. Only if both in the acidic microenvironment and aptamer-targeting marker are present on the cell surface, the sgc8-ASO-chimera strand, carrying an antisense oligonucleotide drug, was released from the nanoprobe and entered into targeted cancer cells for gene silence. Additionally, the in situ drug release facilitated the uptake of drugs mediated by the interaction between sgc8 aptamer and membrane proteins, resulting in enhanced inhibition of cancer cell migration and proliferation. This logic nanoprobe will provide inspiration for designing smart devices for diagnosis of pH-related diseases and targeted drug delivery.

Pre-existing subclones determine radioresistance in rectal cancer organoids

More than half of all patients with cancer receive radiation therapy, but resistance is commonly observed. Currently, it is unknown whether resistance to radiation therapy is acquired or inherently present. Here, we employed organoids derived from rectal cancer and single-cell whole-genome sequencing to investigate the long-term evolution of subclones in response to radiation. Comparing single-cell whole-genome karyotypes between in-vitro-unirradiated and -irradiated organoids revealed three patterns of subclonal evolution: (1) subclonal persistence, (2) subclonal extinction, and (3) subclonal expansion. Organoids in which subclonal shifts occurred (i.e., expansion or extinction) became more resistant to radiation. Although radioresistant subclones did not share recurrent copy-number alterations that could explain their radioresistance, resistance was associated with reduced chromosomal instability, an association that was also observed in 529 human cancer cell lines. These data suggest that resistance to radiation is inherently present and associated with reduced chromosomal instability.

Circulating microRNA Panel for Prediction of Recurrence and Survival in Early-Stage Lung Adenocarcinoma

Early-stage lung adenocarcinoma (LUAD) patients remain at substantial risk for recurrence and disease-related death, highlighting the unmet need of biomarkers for the assessment and identification of those in an early stage who would likely benefit from adjuvant chemotherapy. To identify circulating miRNAs useful for predicting recurrence in early-stage LUAD, we performed miRNA microarray analysis with pools of pretreatment plasma samples from patients with stage I LUAD who developed recurrence or remained recurrence-free during the follow-up period. Subsequent validation in 85 patients with stage I LUAD resulted in the development of a circulating miRNA panel comprising miR-23a-3p, miR-320c, and miR-125b-5p and yielding an area under the curve (AUC) of 0.776 in predicting recurrence. Furthermore, the three-miRNA panel yielded an AUC of 0.804, with a sensitivity of 45.8% at 95% specificity in the independent test set of 57 stage I and II LUAD patients. The miRNA panel score was a significant and independent factor for predicting disease-free survival (p < 0.001, hazard ratio [HR] = 1.64, 95% confidence interval [CI] = 1.51-4.22) and overall survival (p = 0.001, HR = 1.51, 95% CI = 1.17-1.94). This circulating miRNA panel is a useful noninvasive tool to stratify early-stage LUAD patients and determine an appropriate treatment plan with maximal efficacy.

Management of Renal Cell Carcinoma: Promising Biomarkers and the Challenges to Reach the Clinic

The incidence of renal cell carcinoma (RCC) is increasing worldwide, yet research within this field is lagging behind other cancers. Despite increased detection of early disease as a consequence of the widespread use of diagnostic CT scans, 25% of patients have disseminated disease at diagnosis. Similarly, around 25% progress to metastatic disease following curatively intended surgery. Surgery is the cornerstone in the treatment of RCC; however, when the disease is disseminated, immunotherapy or immunotherapy in combination with a tyrosine kinase inhibitor is the patient's best option. Immunotherapy is a potent treatment, with durable treatment responses and potential to cure the patient, but only half of the patients benefit from the administered treatment, and there are currently no methods that can identify which patients will respond to immunotherapy. Moreover, there is a need to identify the patients in greatest risk of relapsing after surgery for localized disease and direct adjuvant treatment there. Even though several molecular biomarkers have been published to date, we are still lacking routinely used biomarkers to guide optimal clinical management. The purpose of this review is to highlight some of the most promising biomarkers, discuss the efforts made within this field to date, and describe the barriers needed to be overcome to have reliable and robust predictive and prognostic biomarkers in the clinic for renal cancer.

Metabolomic biomarkers in liquid biopsy: accurate cancer diagnosis and prognosis monitoring

Liquid biopsy, a novel detection method, has recently become an active research area in clinical cancer owing to its unique advantages. Studies on circulating free DNA, circulating tumor cells, and exosomes obtained by liquid biopsy have shown great advances and they have entered clinical practice as new cancer biomarkers. The metabolism of the body is dynamic as cancer originates and progresses. Metabolic abnormalities caused by cancer can be detected in the blood, sputum, urine, and other biological fluids via systemic or local circulation. A considerable number of recent studies have focused on the roles of metabolic molecules in cancer. The purpose of this review is to provide an overview of metabolic markers from various biological fluids in the latest clinical studies, which may contribute to cancer screening and diagnosis, differentiation of cancer typing, grading and staging, and prediction of therapeutic response and prognosis.