Visualizing Human Colorectal Cancer Intratumor Heterogeneity with Phylogeography

Establishment of patient-derived organoids for guiding personalized therapies in breast cancer patients

Breast cancer has become the most commonly diagnosed cancer. The intra- and interpatient heterogeneity induced a considerable variation in treatment efficacy. There is an urgent requirement for preclinical models to anticipate the effectiveness of individualized drug responses. Patient-derived organoids (PDOs) can accurately recapitulate the architecture and biological characteristics of the origin tumor, making them a promising model that can overtake many limitations of cell lines and PDXs. However, it is still unclear whether PDOs-based drug testing can benefit breast cancer patients, particularly those with tumor recurrence or treatment resistance. Fresh tumor samples were surgically resected for organoid culture. Primary tumor samples and PDOs were subsequently subjected to H&E staining, immunohistochemical (IHC) analysis, and whole-exome sequencing (WES) to make comparisons. Drug sensitivity tests were performed to evaluate the feasibility of this model for predicting patient drug response in clinical practice. We established 75 patient-derived breast cancer organoid models. The results of H&E staining, IHC, and WES revealed that PDOs inherited the histologic and genetic characteristics of their parental tumor tissues. The PDOs successfully predicted the patient's drug response, and most cases exhibited consistency between PDOs' drug susceptibility test results and the clinical response of the matched patient. We conclude that the breast cancer organoids platform can be a potential preclinical tool used for the selection of effective drugs and guided personalized therapies for patients with advanced breast cancer.

Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors

Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments.

An Intratumor Heterogeneity-Related Signature for Predicting Prognosis, Immune Landscape, and Chemotherapy Response in Colon Adenocarcinoma

Background

Colon adenocarcinoma (COAD) is a frequent malignancy of the digestive system with a poor prognosis and high mortality rate worldwide. Intratumor heterogeneity (ITH) is associated with tumor progression, poor prognosis, immunosuppression, and therapy resistance. However, the relationship between ITH and prognosis, the immune microenvironment, and the chemotherapy response in COAD patients remains unknown, and this knowledge is urgently needed.

Methods

We obtained clinical information and gene expression data for COAD patients from The Cancer Genome Atlas (TCGA) database. The DEPTH2 algorithm was utilized to evaluate the ITH score. X-tile software was used to determine the optimal cutoff value of the ITH score. The COAD patients were divided into high- and low-ITH groups based on the cutoff value. We analyzed prognosis, tumor mutation burden (TMB), gene mutations, and immune checkpoint expression between the high- and low-ITH groups. Differentially expressed genes (DEGs) in the high- and low-ITH groups were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. We performed univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses to screen the prognosis-related genes for the construction of an ITH-related prognostic signature. The nomogram was used to predict the overall survival (OS) of COAD patients. The protein–protein interaction (PPI) network was constructed by using the GeneMANIA database. Principal component analysis (PCA) and single-sample gene set enrichment analysis (ssGSEA) were employed to explore the differences in biological pathway activation status between the high- and low-risk groups. The proportion and type of tumor-infiltrating immune cells were evaluated by the CIBERSORT and ESTIMATE algorithms. Additionally, we assessed the chemotherapy response and predicted small-molecule drugs for treatment. Finally, the expression of the prognosis-related genes was validated by using the UALCAN database and Human Protein Atlas (HPA) database.

Results

The OS of the high-ITH group was worse than that of the low-ITH group. A positive correlation between ITH and TMB was identified. In subgroups stratified by age, gender, and tumor stage, the OS of the low-ITH group remained better than that of the high-ITH group. There were dramatic differences in the mutated genes, single nucleotide variant classes, variant types, immune checkpoints and cooccurring and mutually exclusive mutations of the DEGs between the high- and low-ITH groups. Based on the DEGs between the high- and low-ITH groups, we constructed a five-gene signature consisting of CEACAM5, ENO2, GABBR1, MC1R, and SLC44A4. The COAD patients were divided into high- and low-risk groups according to the median risk score. The OS of the high-risk group was worse than that of the low-risk group. The nomogram was used to accurately predict the 1-, 3- and 5-year OS of COAD patients and showed good calibration and moderate discrimination ability. The stromal score, immune score, and ESTIMATE score of the high-risk group were significantly higher than those of the low-risk group, whereas tumor purity showed the opposite trend. The patients classified by the risk score had distinguishable sensitivity to chemotherapeutic drugs. Finally, two public databases confirmed that CEACAM5 and SLC44A4 were upregulated in normal tissues compared with COAD tissues, and ENO2, GABBR1, and MC1R were upregulated in COAD tissues compared with normal tissues.

Conclusion

Overall, we identified an ITH-related prognostic signature for COAD that was closely related to the tumor microenvironment and chemotherapy response. This signature may help clinicians make more personalized and precise treatment decisions for COAD patients.

Pancreatic cancer pathology viewed in the light of evolution

One way to understand ductal adenocarcinoma of the pancreas (pancreatic cancer) is to view it as unimaginably large numbers of evolving living organisms interacting with their environment. This “evolutionary view” creates both expected and surprising perspectives in all stages of neoplastic progression. Advances in the field will require greater attention to this critical evolutionary prospective.

Genetic Alterations of Metastatic Colorectal Cancer

Genome sequencing studies have characterized the genetic alterations of different tumor types, highlighting the diversity of the molecular processes driving tumor development. Comprehensive sequencing studies have defined molecular subtypes of colorectal cancers (CRCs) through the identification of genetic events associated with microsatellite stability (MSS), microsatellite-instability-high (MSI-H), and hypermutation. Most of these studies characterized primary tumors. Only recent studies have addressed the characterization of the genetic and clinical heterogeneity of metastatic CRC. Metastatic CRC genomes were found to be not fundamentally different from primary CRCs in terms of the mutational landscape or of genes that drive tumorigenesis, and a genomic heterogeneity associated with tumor location of primary tumors helps to define different clinical behaviors of metastatic CRCs. Although CRC metastatic spreading was traditionally seen as a late-occurring event, growing evidence suggests that this process can begin early during tumor development and the clonal architecture of these tumors is consistently influenced by cancer treatment. Although the survival rate of patients with metastatic CRC patients improved in the last years, the response to current treatments and prognosis of many of these patients remain still poor, indicating the need to discover new improvements for therapeutic vulnerabilities and to formulate a rational prospective of personalized therapies.