Spatially organized multicellular immune hubs in human colorectal cancer

Single-Cell Profiling Reveals the Impact of Genetic Alterations on the Differentiation of Inflammation-Induced Murine Colon Tumors

Genetic mutations and chronic inflammation of the colon contribute to the development of colorectal cancer (CRC). Using a murine model of inflammation-induced colon tumorigenesis, we determined how genetic mutations alter colon tumor cell differentiation. Inflammation induced by enterotoxigenic Bacteroides fragilis (ETBF) colonization of multiple intestinal neoplasia (MinApcΔ716/+) mice triggers loss of heterozygosity of Apc causing colon tumor formation. Here, we report that the addition of BRAFV600E mutation (BRAFF-V600ELgr5tm1(Cre/ERT2)CleMinApcΔ716/+, BLM) or knocking out Msh2 (Msh2LoxP/LoxPVil1-creMinApcΔ716/+, MSH2KO) in the Min model altered colon tumor differentiation. Using single-cell RNA sequencing, we uncovered the differences between BLM, Min, and MSH2KO tumors at a single-cell resolution. BLM tumors showed an increase in differentiated tumor epithelial cell lineages and a reduction in the tumor stem cell population. Interestingly, the tumor stem cell population of BLM tumors had revival colon stem cell characteristics with low WNT signaling and an increase in RevCSC marker gene expression. In contrast, MSH2KO tumors were characterized by an increased tumor stem cell population that had higher WNT signaling activity compared to Min tumors. Furthermore, overall BLM tumors had higher expression of transcription factors that drive differentiation, such as Cdx2, than Min tumors. Using RNA velocity, we identified additional potential regulators of BLM tumor differentiation such as NDRG1. The role of CDX2 and NDRG1 as putative regulators for BLM tumor cell differentiation was verified using organoids derived from BLM tumors. Our results demonstrate the critical connections between genetic mutations and cell differentiation in inflammation-induced colon tumorigenesis. Understanding such roles will deepen our understanding of inflammation-associated colon cancer.

Resistance to immune checkpoint inhibitors in colorectal cancer with deficient mismatch repair/microsatellite instability: misdiagnosis, pseudoprogression and/or tumor heterogeneity?

Colorectal carcinoma (CRC) with deficiency of the deficient mismatch repair (dMMR) pathway/microsatellite instability (MSI) is characterized by a high mutation load and infiltration of immune cells in the tumor microenvironment. In agreement with these findings, clinical trials have demonstrated a significant activity of immune checkpoint inhibitors (ICIs) in dMMR/MSI metastatic CRC (mCRC) patients and, more recently, in CRC patients with early disease undergoing neoadjuvant therapy. However, despite high response rates and durable clinical benefits, a fraction of mCRC patients, up to 30%, showed progressive disease when treated with single agent anti-programmed cell death 1 (PD-1) antibody. This article discusses the three main causes that have been associated with early progression of dMMR/MSI mCRC patients while on treatment with ICIs, i.e., misdiagnosis, pseudoprogression and tumor heterogeneity. While pseudoprogression probably does not play a relevant role, data from clinical studies demonstrate that some dMMR/MSI CRC cases with rapid progression on ICIs may be misdiagnosed, underlining the importance of correct diagnostics. More importantly, evidence suggests that dMMR/MSI mCRC is a heterogeneous group of tumors with different sensitivity to ICIs. Therefore, we propose novel diagnostic and therapeutic strategies to improve the outcome of dMMR/MSI CRC patients.

IGSF8 is an innate immune checkpoint and cancer immunotherapy target

Antigen presentation defects in tumors are prevalent mechanisms of adaptive immune evasion and resistance to cancer immunotherapy, whereas how tumors evade innate immunity is less clear. Using CRISPR screens, we discovered that IGSF8 expressed on tumors suppresses NK cell function by interacting with human KIR3DL2 and mouse Klra9 receptors on NK cells. IGSF8 is normally expressed in neuronal tissues and is not required for cell survival in vitro or in vivo. It is overexpressed and associated with low antigen presentation, low immune infiltration, and worse clinical outcomes in many tumors. An antibody that blocks IGSF8-NK receptor interaction enhances NK cell killing of malignant cells in vitro and upregulates antigen presentation, NK cell-mediated cytotoxicity, and T cell signaling in vivo. In syngeneic tumor models, anti-IGSF8 alone, or in combination with anti-PD1, inhibits tumor growth. Our results indicate that IGSF8 is an innate immune checkpoint that could be exploited as a therapeutic target.

New Insight into Intestinal Mast Cells Revealed by Single-Cell RNA Sequencing

Mast cells (MCs) are tissue-resident immune cells distributed in all tissues and strategically located close to blood and lymphatic vessels and nerves. Thanks to the expression of a wide array of receptors, MCs act as tissue sentinels, able to detect the presence of bacteria and parasites and to respond to different environmental stimuli. MCs originate from bone marrow (BM) progenitors that enter the circulation and mature in peripheral organs under the influence of microenvironment factors, thus differentiating into heterogeneous tissue-specific subsets. Even though MC activation has been traditionally linked to IgE-mediated allergic reactions, a role for these cells in other pathological conditions including tumor progression has recently emerged. However, several aspects of MC biology remain to be clarified. The advent of single-cell RNA sequencing platforms has provided the opportunity to understand MCs' origin and differentiation as well as their phenotype and functions within different tissues, including the gut. This review recapitulates how single-cell transcriptomic studies provided insight into MC development as well as into the functional role of intestinal MC subsets in health and disease.

Multiregional transcriptomics identifies congruent consensus subtypes with prognostic value beyond tumor heterogeneity of colorectal cancer

Intra-tumor heterogeneity compromises the clinical value of transcriptomic classifications of colorectal cancer. We investigated the prognostic effect of transcriptomic heterogeneity and the potential for classifications less vulnerable to heterogeneity in a single-hospital series of 1093 tumor samples from 692 patients, including multiregional samples from 98 primary tumors and 35 primary-metastasis sets. We show that intra-tumor heterogeneity of the consensus molecular subtypes (CMS) is frequent and has poor-prognostic associations independently of tumor microenvironment markers. Multiregional transcriptomics uncover cancer cell-intrinsic and low-heterogeneity signals that recapitulate the intrinsic CMSs proposed by single-cell sequencing. Further subclassification identifies congruent CMSs that explain a larger proportion of variation in patient survival than intra-tumor heterogeneity. Plasticity is indicated by discordant intrinsic phenotypes of matched primary and metastatic tumors. We conclude that multiregional sampling reconciles the prognostic power of tumor classifications from single-cell and bulk transcriptomics in the context of intra-tumor heterogeneity, and phenotypic plasticity challenges the reconciliation of primary and metastatic subtypes.

Immune Cell Migration to Cancer

Immune cell migration is required for the development of an effective and robust immune response. This elegant process is regulated by both cellular and environmental factors, with variables such as immune cell state, anatomical location, and disease state that govern differences in migration patterns. In all cases, a major factor is the expression of cell surface receptors and their cognate ligands. Rapid adaptation to environmental conditions partly depends on intrinsic cellular immune factors that affect a cell's ability to adjust to new environment. In this review, we discuss both myeloid and lymphoid cells and outline key determinants that govern immune cell migration, including molecules required for immune cell adhesion, modes of migration, chemotaxis, and specific chemokine signaling. Furthermore, we summarize tumor-specific elements that contribute to immune cell trafficking to cancer, while also exploring microenvironment factors that can alter these cellular dynamics within the tumor in both a pro and antitumor fashion. Specifically, we highlight the importance of the secretome in these later aspects. This review considers a myriad of factors that impact immune cell trajectory in cancer. We aim to highlight the immunotherapeutic targets that can be harnessed to achieve controlled immune trafficking to and within tumors.

A blueprint for tumor-infiltrating B cells across human cancers

B lymphocytes are essential mediators of humoral immunity and play multiple roles in human cancer. To decode the functions of tumor-infiltrating B cells, we generated a B cell blueprint encompassing single-cell transcriptome, B cell-receptor repertoire, and chromatin accessibility data across 20 different cancer types (477 samples, 269 patients). B cells harbored extraordinary heterogeneity and comprised 15 subsets, which could be grouped into two independent developmental paths (extrafollicular versus germinal center). Tumor types grouped into the extrafollicular pathway were linked with worse clinical outcomes and resistance to immunotherapy. The dysfunctional extrafollicular program was associated with glutamine-derived metabolites through epigenetic-metabolic cross-talk, which promoted a T cell-driven immunosuppressive program. These data suggest an intratumor B cell balance between extrafollicular and germinal-center responses and suggest that humoral immunity could possibly be harnessed for B cell-targeting immunotherapy.

ILT2 and ILT4 Drive Myeloid Suppression via Both Overlapping and Distinct Mechanisms

Solid tumors are dense three-dimensional (3D) multicellular structures that enable efficient receptor-ligand trans interactions via close cell-cell contact. Immunoglobulin-like transcript (ILT)2 and ILT4 are related immune-suppressive receptors that play a role in the inhibition of myeloid cells within the tumor microenvironment. The relative contribution of ILT2 and ILT4 to immune inhibition in the context of solid tumor tissue has not been fully explored. We present evidence that both ILT2 and ILT4 contribute to myeloid inhibition. We found that although ILT2 inhibits myeloid cell activation in the context of trans-engagement by MHC-I, ILT4 efficiently inhibits myeloid cells in the presence of either cis- or trans-engagement. In a 3D spheroid tumor model, dual ILT2/ILT4 blockade was required for the optimal activation of myeloid cells, including the secretion of CXCL9 and CCL5, upregulation of CD86 on dendritic cells, and downregulation of CD163 on macrophages. Humanized mouse tumor models showed increased immune activation and cytolytic T-cell activity with combined ILT2 and ILT4 blockade, including evidence of the generation of immune niches, which have been shown to correlate with clinical response to immune-checkpoint blockade. In a human tumor explant histoculture system, dual ILT2/ILT4 blockade increased CXCL9 secretion, downregulated CD163 expression, and increased the expression of M1 macrophage, IFNγ, and cytolytic T-cell gene signatures. Thus, we have revealed distinct contributions of ILT2 and ILT4 to myeloid cell biology and provide proof-of-concept data supporting the combined blockade of ILT2 and ILT4 to therapeutically induce optimal myeloid cell reprogramming in the tumor microenvironment.

Microenvironmental changes in familial adenomatous polyposis during colorectal cancer carcinogenesis

Familial adenomatous polyposis (FAP) is a heritable disease that increases the risk of colorectal cancer (CRC) development because of heterozygous mutations in APC. Little is known about the microenvironment of FAP. Here, single-cell RNA sequencing was performed on matched normal tissues, adenomas, and carcinomas from four patients with FAP. We analyzed the transcriptomes of 56,225 unsorted single cells, revealing the heterogeneity of each cell type, and compared gene expression among tissues. Then we compared the gene expression with that of sporadic CRC. Furthermore, we analyzed specimens of 26 FAP patients and 40 sporadic CRC patients by immunohistochemistry. Immunosuppressiveness of myeloid cells, fibroblasts, and regulatory T cells was upregulated even in the early stages of carcinogenesis. CD8+ T cells became exhausted only in carcinoma, although the cytotoxicity of CD8+ T cells was gradually increased according to the carcinogenic step. When compared with those in the sporadic CRC microenvironment, the composition and function of each cell type in the FAP-derived CRC microenvironment had differences. Our findings indicate that an immunosuppressive microenvironment is constructed from a precancerous stage in FAP.

A spatiotemporal atlas of mouse liver homeostasis and regeneration

The mechanism by which mammalian liver cell responses are coordinated during tissue homeostasis and perturbation is poorly understood, representing a major obstacle in our understanding of many diseases. This knowledge gap is caused by the difficulty involved with studying multiple cell types in different states and locations, particularly when these are transient. We have combined Stereo-seq (spatiotemporal enhanced resolution omics-sequencing) with single-cell transcriptomic profiling of 473,290 cells to generate a high-definition spatiotemporal atlas of mouse liver homeostasis and regeneration at the whole-lobe scale. Our integrative study dissects in detail the molecular gradients controlling liver cell function, systematically defining how gene networks are dynamically modulated through intercellular communication to promote regeneration. Among other important regulators, we identified the transcriptional cofactor TBL1XR1 as a rheostat linking inflammation to Wnt/β-catenin signaling for facilitating hepatocyte proliferation. Our data and analytical pipelines lay the foundation for future high-definition tissue-scale atlases of organ physiology and malfunction.

Tissue-of-origin for cancers determines HIF-1 activation induced phenotypic heterogeneity

Hypoxia-inducible factor-1 (HIF-1) is the master regulator of cellular response to hypoxia, and is activated in many cancers contributing to many steps in the metastatic cascade by acting as a key transcription co-regulator for a large number of downstream genes. Presence of hypoxia within a tumor is spatially nonuniform, and can also by dynamic. Further, although HIF-1 is primarily stabilized and activated by lack of molecular O2, its stability is also affected by other factors present in the tumor microenvironment. HIF-1 also crosstalks with other transcription factors in co-regulating gene expression. Consequently, it is nontrivial to predict the gene expression patterns in cells in response to hypoxia, or HIF-1 activation. Additionally, cancers originating from tissue origins with different basal level of partial oxygen tension may activate HIF-1 at different threshold of hypoxia. We analyzed large published single cell RNAseq data for colorectal, lung, and pancreatic cancers to investigate the phenotypic outcome of HIF-1 activation in cancer cells. We found that cancers from tissues with different partial O2 tension levels exhibit HIF-1 activation at different stages of metastasis, and phenotypically respond differently to HIF-1 activation, likely by contextual co-option of different transcription factors. We experimentally confirmed these predictions by using cell lines representative of colorectal, lung, and pancreatic cancers, finding that while hypoxia enhances growth of colorectal cancer, it induces increased invasion of lung, and pancreatic cancers. Our analysis suggest that HIF-1 activation may act as a rheostat regulating downstream gene expression towards phenotypic outcomes differently in various cancers.

Early Immune Remodeling Steers Clinical Response to First-Line Chemoimmunotherapy in Advanced Gastric Cancer

Adding anti-programmed cell death protein 1 (anti-PD-1) to 5-fluorouracil (5-FU)/platinum improves survival in some advanced gastroesophageal adenocarcinomas (GEA). To understand the effects of chemotherapy and immunotherapy, we conducted a phase II first-line trial (n = 47) sequentially adding pembrolizumab to 5-FU/platinum in advanced GEA. Using serial biopsy of the primary tumor at baseline, after one cycle of 5-FU/platinum, and after the addition of pembrolizumab, we transcriptionally profiled 358,067 single cells to identify evolving multicellular tumor microenvironment (TME) networks. Chemotherapy induced early on-treatment multicellular hubs with tumor-reactive T-cell and M1-like macrophage interactions in slow progressors. Faster progression featured increased MUC5A and MSLN containing treatment resistance programs in tumor cells and M2-like macrophages with immunosuppressive stromal interactions. After pembrolizumab, we observed increased CD8 T-cell infiltration and development of an immunity hub involving tumor-reactive CXCL13 T-cell program and epithelial interferon-stimulated gene programs. Strategies to drive increases in antitumor immune hub formation could expand the portion of patients benefiting from anti-PD-1 approaches.

SIGNIFICANCE

The benefit of 5-FU/platinum with anti-PD-1 in first-line advanced gastric cancer is limited to patient subgroups. Using a trial with sequential anti-PD-1, we show coordinated induction of multicellular TME hubs informs the ability of anti-PD-1 to potentiate T cell-driven responses. Differential TME hub development highlights features that underlie clinical outcomes. This article is featured in Selected Articles from This Issue, p. 695.

De novo detection of somatic mutations in high-throughput single-cell profiling data sets

Characterization of somatic mutations at single-cell resolution is essential to study cancer evolution, clonal mosaicism and cell plasticity. Here, we describe SComatic, an algorithm designed for the detection of somatic mutations in single-cell transcriptomic and ATAC-seq (assay for transposase-accessible chromatin sequence) data sets directly without requiring matched bulk or single-cell DNA sequencing data. SComatic distinguishes somatic mutations from polymorphisms, RNA-editing events and artefacts using filters and statistical tests parameterized on non-neoplastic samples. Using >2.6 million single cells from 688 single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) data sets spanning cancer and non-neoplastic samples, we show that SComatic detects mutations in single cells accurately, even in differentiated cells from polyclonal tissues that are not amenable to mutation detection using existing methods. Validated against matched genome sequencing and scRNA-seq data, SComatic achieves F1 scores between 0.6 and 0.7 across diverse data sets, in comparison to 0.2-0.4 for the second-best performing method. In summary, SComatic permits de novo mutational signature analysis, and the study of clonal heterogeneity and mutational burdens at single-cell resolution.

Single-cell transcriptomic analyses reveal distinct immune cell contributions to epithelial barrier dysfunction in checkpoint inhibitor colitis

Immune checkpoint inhibitor (ICI) therapy has revolutionized oncology, but treatments are limited by immune-related adverse events, including checkpoint inhibitor colitis (irColitis). Little is understood about the pathogenic mechanisms driving irColitis, which does not readily occur in model organisms, such as mice. To define molecular drivers of irColitis, we used single-cell multi-omics to profile approximately 300,000 cells from the colon mucosa and blood of 13 patients with cancer who developed irColitis (nine on anti-PD-1 or anti-CTLA-4 monotherapy and four on dual ICI therapy; most patients had skin or lung cancer), eight controls on ICI therapy and eight healthy controls. Patients with irColitis showed expanded mucosal Tregs, ITGAEHi CD8 tissue-resident memory T cells expressing CXCL13 and Th17 gene programs and recirculating ITGB2Hi CD8 T cells. Cytotoxic GNLYHi CD4 T cells, recirculating ITGB2Hi CD8 T cells and endothelial cells expressing hypoxia gene programs were further expanded in colitis associated with anti-PD-1/CTLA-4 therapy compared to anti-PD-1 therapy. Luminal epithelial cells in patients with irColitis expressed PCSK9, PD-L1 and interferon-induced signatures associated with apoptosis, increased cell turnover and malabsorption. Together, these data suggest roles for circulating T cells and epithelial-immune crosstalk critical to PD-1/CTLA-4-dependent tolerance and barrier function and identify potential therapeutic targets for irColitis.

Identification of the molecular subtypes and signatures to predict the prognosis, biological functions, and therapeutic response based on the anoikis-related genes in colorectal cancer

BACKGROUND

Tumors that resist anoikis, a programmed cell death triggered by detachment from the extracellular matrix, promote metastasis; however, the role of anoikis-related genes (ARGs) in colorectal cancer (CRC) stratification, prognosis, and biological functions remains unclear.

METHODS

We obtained transcriptomic profiles of CRC and 27 ARGs from The Cancer Genome Atlas, the Gene Expression Omnibus, and MSigDB databases, respectively. CRC tissue samples were classified into two clusters based on the expression pattern of ARGs, and their functional differences were explored. Hub genes were screened using weighted gene co-expression network analysis, univariate analysis, and least absolute selection and shrinkage operator analysis, and validated in cell lines, tissues, or the Human Protein Atlas database. We constructed an ARG-risk model and nomogram to predict prognosis in patients with CRC, which was validated using an external cohort. Multifaceted landscapes, including stemness, tumor microenvironment (TME), immune landscape, and drug sensitivity, between high- and low-risk groups were examined.

RESULTS

Patients with CRC were divided into C1 and C2 clusters. Cluster C1 exhibited higher TME scores, whereas cluster C2 had favorable outcomes and a higher stemness index. Eight upregulated hub ARGs (TIMP1, P3H1, SPP1, HAMP, IFI30, ADAM8, ITGAX, and APOC1) were utilized to construct the risk model. The qRT-PCR, Western blotting, and immunohistochemistry results were consistent with those of the bioinformatics analysis. Patients with high risk exhibited worse overall survival (p < 0.01), increased stemness, TME, immune checkpoint expression, immune infiltration, tumor mutation burden, and drug susceptibility compared with the patients with low risk.

CONCLUSION

Our results offer a novel CRC stratification based on ARGs and a risk-scoring system that could predict the prognosis, stemness, TME, immunophenotypes, and drug susceptibility of patients with CRC, thereby improving their prognosis. This stratification may facilitate personalized therapies.

Human NLRC4 expression promotes cancer survival and associates with type I interferon signaling and immune infiltration

The immune system can control cancer progression. However, even though some innate immune sensors of cellular stress are expressed intrinsically in epithelial cells, their potential role in cancer aggressiveness and subsequent overall survival in humans is mainly unknown. Here, we show that nucleotide-binding oligomerization domain-like receptor (NLR) family CARD domain-containing 4 (NLRC4) is downregulated in epithelial tumor cells of patients with colorectal cancer (CRC) by using spatial tissue imaging. Strikingly, only the loss of tumor NLRC4, but not stromal NLRC4, was associated with poor immune infiltration (mainly DCs and CD4+ and CD8+ T cells) and accurately predicted progression to metastatic stage IV and decrease in overall survival. By combining multiomics approaches, we show that restoring NLRC4 expression in human CRC cells triggered a broad inflammasome-independent immune reprogramming consisting of type I interferon (IFN) signaling genes and the release of chemokines and myeloid growth factors involved in the tumor infiltration and activation of DCs and T cells. Consistently, such reprogramming in cancer cells was sufficient to directly induce maturation of human DCs toward a Th1 antitumor immune response through IL-12 production in vitro. In multiple human carcinomas (colorectal, lung, and skin), we confirmed that NLRC4 expression in patient tumors was strongly associated with type I IFN genes, immune infiltrates, and high microsatellite instability. Thus, we shed light on the epithelial innate immune sensor NLRC4 as a therapeutic target to promote an efficient antitumor immune response against the aggressiveness of various carcinomas.

Spatiotemporal architecture of immune cells and cancer-associated fibroblasts in high-grade serous ovarian carcinoma

High-grade serous ovarian carcinoma (HGSOC), the deadliest form of ovarian cancer, is typically diagnosed after it has metastasized and often relapses after standard-of-care platinum-based chemotherapy, likely due to advanced tumor stage, heterogeneity, and immune evasion and tumor-promoting signaling from the tumor microenvironment. To understand how spatial heterogeneity contributes to HGSOC progression and early relapse, we profiled an HGSOC tissue microarray of patient-matched longitudinal samples from 42 patients. We found spatial patterns associated with early relapse, including changes in T cell localization, malformed tertiary lymphoid structure (TLS)-like aggregates, and increased podoplanin-positive cancer-associated fibroblasts (CAFs). Using spatial features to compartmentalize the tissue, we found that plasma cells distribute in two different compartments associated with TLS-like aggregates and CAFs, and these distinct microenvironments may account for the conflicting reports about the role of plasma cells in HGSOC prognosis.

Molecular Landscape and Therapeutic Strategies against Colorectal Cancer

Colorectal cancer (CRC) is the second leading cause of cancer deaths worldwide. Although the overall incidence of CRC is decreasing, the incidence of young-onset CRC, characterized by a diagnosis of CRC before age 50, is increasing. Outcomes for CRC patients are improving, partly due to comprehensive molecular characterization of tumors and novel therapeutic strategies. Advances in genomic and transcriptomic analyses using blood- and tumor-tissue-based sequencing have facilitated identification of distinct tumor subtypes harboring unique biological characteristics and therapeutic vulnerabilities. These insights have led to the development and incorporation of targeted therapies and immunotherapy in CRC treatment. In this review, we discuss the molecular landscape and key oncogenes/tumor suppressors contributing to CRC tumorigenesis, metastasis, and therapeutic resistance. We also discuss personalized therapeutic strategies for subsets of CRC patients and provide an overview of evolving novel treatments being evaluated in clinical trials.

Single-cell transcriptome analysis profiling lymphatic invasion-related TME in colorectal cancer

Lymphatic invasion (LI) is extremely aggressive and induces worse prognosis among patients with colorectal cancer (CRC). Thus, it is critical to characterize the cellular and molecular mechanisms underlying LI in order to establish novel and efficacious therapeutic targets that enhance the prognosis of CRC patients. RNA-seq data, clinical and survival information of colon adenocarcinoma (COAD) patients were obtained from the TCGA database. In addition, three scRNA-seq datasets of CRC patients were acquired from the GEO database. Data analyses were conducted with the R packages. We assessed the tumor microenvironment (TME) differences between LI+ and LI- based scRNA-seq data, LI+ cells exhibited augmented abundance of immunosuppression and invasive subset. Marked extracellular matrix network activation was also observed in LI+ cells within SPP1+ macrophages. We revealed that an immunosuppressive and pro-angiogenic TME strongly enhanced LI, as was evidenced by the CD4+ Tregs, CD8+ GZMK+, SPP1+ macrophages, e-myCAFs, and w-myCAFs subcluster infiltrations. Furthermore, we identified potential LI targets that influenced tumor development, metastasis, and immunotherapeutic response. Finally, a novel LIRS model was established based on the expression of 14 LI-related signatures, and in the two testing cohorts, LIRS was also proved to have accurate prognostic predictive ability. In this report, we provided a valuable resource and extensive insights into the LI of CRC. Our conclusions can potentially benefit the establishment of highly efficacious therapeutic targets as well as diagnostic biomarkers that improve patient outcomes.

Th17 Cells Secrete TWEAK to Trigger Epithelial-Mesenchymal Transition and Promote Colorectal Cancer Liver Metastasis

Liver metastasis is the leading cause of mortality in patients with colorectal cancer. Given the significance of both epithelial-mesenchymal transition (EMT) of tumor cells and the immune microenvironment in colorectal cancer liver metastasis (CRLM), the interplay between them could hold the key for developing improved treatment options. We employed multiomics analysis of 130 samples from 18 patients with synchronous CRLM integrated with external datasets to comprehensively evaluate the interaction between immune cells and EMT of tumor cells in liver metastasis. Single-cell RNA sequencing analysis revealed distinct distributions of nonmalignant cells between primary tumors from patients with metastatic colorectal cancer (mCRC) and non-metastatic colorectal cancer, showing that Th17 cells were predominantly enriched in the primary lesion of mCRC. TWEAK, a cytokine secreted by Th17 cells, promoted EMT by binding to receptor Fn14 on tumor cells, and the TWEAK-Fn14 interaction enhanced tumor migration and invasion. In mouse models, targeting Fn14 using CRISPR-induced knockout or lipid nanoparticle-encapsulated siRNA alleviated metastasis and prolonged survival. Mice lacking Il17a or Tnfsf12 (encoding TWEAK) exhibited fewer metastases compared with wild-type mice, while cotransfer of Th17 with tumor cells promoted liver metastasis. Higher TWEAK expression was associated with a worse prognosis in patients with colorectal cancer. In addition, CD163L1+ macrophages interacted with Th17 cells, recruiting Th17 via the CCL4-CCR5 axis. Collectively, this study unveils the role of immune cells in the EMT process and identifies TWEAK secreted by Th17 as a driver of CRLM.

SIGNIFICANCE

TWEAK secreted by Th17 cells promotes EMT by binding to Fn14 on colorectal cancer cells, suggesting that blocking the TWEAK-Fn14 interaction may be a promising therapeutic approach to inhibit liver metastasis.

A specific enterotype derived from gut microbiome of older individuals enables favorable responses to immune checkpoint blockade therapy

Immunotherapy has revolutionized cancer treatment, but inconsistent responses persist. Our study delves into the intriguing phenomenon of enhanced immunotherapy sensitivity in older individuals with cancers. Through a meta-analysis encompassing 25 small-to-mid-sized trials of immune checkpoint blockade (ICB), we demonstrate that older individuals exhibit heightened responsiveness to ICB therapy. To understand the underlying mechanism, we reanalyze single-cell RNA sequencing (scRNA-seq) data from multiple studies and unveil distinct upregulation of exhausted and cytotoxic T cell markers within the tumor microenvironment (TME) of older patients. Recognizing the potential role of gut microbiota in modulating the efficacy of immunotherapy, we identify an aging-enriched enterotype linked to improved immunotherapy outcomes in older patients. Fecal microbiota transplantation experiments in mice confirm the therapeutic potential of the aging-enriched enterotype, enhancing treatment sensitivity and reshaping the TME. Our discoveries confront the prevailing paradox and provide encouraging paths for tailoring cancer immunotherapy strategies according to an individual's gut microbiome profile.

CD177 drives the transendothelial migration of Treg cells enriched in human colorectal cancer

Objectives

Regulatory T (Treg) cells regulate immunity in autoimmune diseases and cancers. However, immunotherapies that target tumor-infiltrating Treg cells often induce unwanted immune responses and tissue inflammation. Our research focussed on exploring the expression pattern of CD177 in tumor-infiltrating Treg cells with the aim of identifying a potential target that can enhance immunotherapy effectiveness.

Methods

Single-cell RNA sequencing (scRNA-seq) data and survival data were obtained from public databases. Twenty-one colorectal cancer patient samples, including fresh tumor tissues, peritumoral tissues and peripheral blood mononuclear cells (PBMCs), were analysed using flow cytometry. The transendothelial activity of CD177+ Treg cells was substantiated using in vitro experiments.

Results

ScRNA-seq and flow cytometry results indicated that CD177 was exclusively expressed in intratumoral Treg cells. CD177+ Treg cells exhibited greater activation status and expressed elevated Treg cell canonical markers and immune checkpoint molecules than CD177- Treg cells. We further discovered that both intratumoral CD177+ Treg cells and CD177-overexpressing induced Treg (iTreg) cells had lower levels of PD-1 than their CD177- counterparts. Moreover, CD177 overexpression significantly enhanced the transendothelial migration of Treg cells in vitro.

Conclusions

These results demonstrated that Treg cells with higher CD177 levels exhibited an enhanced activation status and transendothelial migration capacity. Our findings suggest that CD177 may serve as an immunotherapeutic target and that overexpression of CD177 may improve the efficacy of chimeric antigen receptor T (CAR-T) cell therapy.

scRNAseq and High-Throughput Spatial Analysis of Tumor and Normal Microenvironment in Solid Tumors Reveal a Possible Origin of Circulating Tumor Hybrid Cells

Metastatic cancer is a leading cause of death in cancer patients worldwide. While circulating hybrid cells (CHCs) are implicated in metastatic spread, studies documenting their tissue origin remain sparse, with limited candidate approaches using one-two markers. Utilizing high-throughput single-cell and spatial transcriptomics, we identified tumor hybrid cells (THCs) co-expressing epithelial and macrophage markers and expressing a distinct transcriptome. Rarely, normal tissue showed these cells (NHCs), but their transcriptome was easily distinguishable from THCs. THCs with unique transcriptomes were observed in breast and colon cancers, suggesting this to be a generalizable phenomenon across cancer types. This study establishes a framework for HC identification in large datasets, providing compelling evidence for their tissue residence and offering comprehensive transcriptomic characterization. Furthermore, it sheds light on their differential function and identifies pathways that could explain their newly acquired invasive capabilities. THCs should be considered as potential therapeutic targets.

Molecular insights into clinical trials for immune checkpoint inhibitors in colorectal cancer: Unravelling challenges and future directions

Colorectal cancer (CRC) is a complex disease with diverse etiologies and clinical outcomes. Despite considerable progress in development of CRC therapeutics, challenges remain regarding the diagnosis and management of advanced stage metastatic CRC (mCRC). In particular, the five-year survival rate is very low since mCRC is currently rarely curable. Over the past decade, cancer treatment has significantly improved with the introduction of cancer immunotherapies, specifically immune checkpoint inhibitors. Therapies aimed at blocking immune checkpoints such as PD-1, PD-L1, and CTLA-4 target inhibitory pathways of the immune system, and thereby enhance anti-tumor immunity. These therapies thus have shown promising results in many clinical trials alone or in combination. The efficacy and safety of immunotherapy, either alone or in combination with CRC, have been investigated in several clinical trials. Clinical trials, including KEYNOTE-164 and CheckMate 142, have led to Food and Drug Administration approval of the PD-1 inhibitors pembrolizumab and nivolumab, respectively, for the treatment of patients with unresectable or metastatic microsatellite instability-high or deficient mismatch repair CRC. Unfortunately, these drugs benefit only a small percentage of patients, with the benefits of immunotherapy remaining elusive for the vast majority of CRC patients. To this end, primary and secondary resistance to immunotherapy remains a significant issue, and further research is necessary to optimize the use of immunotherapy in CRC and identify biomarkers to predict the response. This review provides a comprehensive overview of the clinical trials involving immune checkpoint inhibitors in CRC. The underlying rationale, challenges faced, and potential future steps to improve the prognosis and enhance the likelihood of successful trials in this field are discussed.

Identification of Colorectal Cancer Cell Stemness from Single-Cell RNA Sequencing

Cancer stem cells (CSC) play a critical role in metastasis, relapse, and therapy resistance in colorectal cancer. While characterization of the normal lineage of cell development in the intestine has led to the identification of many genes involved in the induction and maintenance of pluripotency, recent studies suggest significant heterogeneity in CSC populations. Moreover, while many canonical colorectal cancer CSC marker genes have been identified, the ability to use these classical markers to annotate stemness at the single-cell level is limited. In this study, we performed single-cell RNA sequencing on a cohort of 6 primary colon, 9 liver metastatic tumors, and 11 normal (nontumor) controls to identify colorectal CSCs at the single-cell level. Finding poor alignment of the 11 genes most used to identify colorectal CSC, we instead extracted a single-cell stemness signature (SCS_sig) that robustly identified "gold-standard" colorectal CSCs that expressed all marker genes. Using this SCS_sig to quantify stemness, we found that while normal epithelial cells show a bimodal distribution, indicating distinct stem and differentiated states, in tumor epithelial cells stemness is a continuum, suggesting greater plasticity in these cells. The SCS_sig score was quite variable between different tumors, reflective of the known transcriptomic heterogeneity of CRC. Notably, patients with higher SCS_sig scores had significantly shorter disease-free survival time after curative intent surgical resection, suggesting stemness is associated with relapse.

IMPLICATIONS

This study reveals significant heterogeneity of expression of genes commonly used to identify colorectal CSCs, and identifies a novel stemness signature to identify these cells from scRNA-seq data.

q-Diffusion leverages the full dimensionality of gene coexpression in single-cell transcriptomics

Unlocking the full dimensionality of single-cell RNA sequencing data (scRNAseq) is the next frontier to a richer, fuller understanding of cell biology. We introduce q-diffusion, a framework for capturing the coexpression structure of an entire library of genes, improving on state-of-the-art analysis tools. The method is demonstrated via three case studies. In the first, q-diffusion helps gain statistical significance for differential effects on patient outcomes when analyzing the CALGB/SWOG 80405 randomized phase III clinical trial, suggesting precision guidance for the treatment of metastatic colorectal cancer. Secondly, q-diffusion is benchmarked against existing scRNAseq classification methods using an in vitro PBMC dataset, in which the proposed method discriminates IFN-γ stimulation more accurately. The same case study demonstrates improvements in unsupervised cell clustering with the recent Tabula Sapiens human atlas. Finally, a local distributional segmentation approach for spatial scRNAseq, driven by q-diffusion, yields interpretable structures of human cortical tissue.

Tissue-specific reprogramming leads to angiogenic neutrophil specialization and tumor vascularization in colorectal cancer

Neutrophil (PMN) tissue accumulation is an established feature of ulcerative colitis (UC) lesions and colorectal cancer (CRC). To assess the PMN phenotypic and functional diversification during the transition from inflammatory ulceration to CRC we analyzed the transcriptomic landscape of blood and tissue PMNs. Transcriptional programs effectively separated PMNs based on their proximity to peripheral blood, inflamed colon, and tumors. In silico pathway overrepresentation analysis, protein-network mapping, gene signature identification, and gene-ontology scoring revealed unique enrichment of angiogenic and vasculature development pathways in tumor-associated neutrophils (TANs). Functional studies utilizing ex vivo cultures, colitis-induced murine CRC, and patient-derived xenograft models demonstrated a critical role for TANs in promoting tumor vascularization. Spp1 (OPN) and Mmp14 (MT1-MMP) were identified by unbiased -omics and mechanistic studies to be highly induced in TANs, acting to critically regulate endothelial cell chemotaxis and branching. TCGA data set and clinical specimens confirmed enrichment of SPP1 and MMP14 in high-grade CRC but not in patients with UC. Pharmacological inhibition of TAN trafficking or MMP14 activity effectively reduced tumor vascular density, leading to CRC regression. Our findings demonstrate a niche-directed PMN functional specialization and identify TAN contributions to tumor vascularization, delineating what we believe to be a new therapeutic framework for CRC treatment focused on TAN angiogenic properties.

Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Dysregulation of Wnt/β-catenin signaling contributes to intestinal inflammation through regulation of group 3 innate lymphoid cells

RORγt+ group 3 innate lymphoid cells (ILC3s) are essential for intestinal homeostasis. Dysregulation of ILC3s has been found in the gut of patients with inflammatory bowel disease and colorectal cancer, yet the specific mechanisms still require more investigation. Here we observe increased β-catenin in intestinal ILC3s from inflammatory bowel disease and colon cancer patients compared with healthy donors. In contrast to promoting RORγt expression in T cells, activation of Wnt/β-catenin signaling in ILC3s suppresses RORγt expression, inhibits its proliferation and function, and leads to a deficiency of ILC3s and subsequent intestinal inflammation in mice. Activated β-catenin and its interacting transcription factor, TCF-1, cannot directly suppress RORγt expression, but rather alters global chromatin accessibility and inhibits JunB expression, which is essential for RORγt expression in ILC3s. Together, our findings suggest that dysregulated Wnt/β-catenin signaling impairs intestinal ILC3s through TCF-1/JunB/RORγt regulation, further disrupting intestinal homeostasis, and promoting inflammation and cancer.

Epithelial regulation of microbiota-immune cell dynamics

The mammalian gastrointestinal tract hosts a diverse community of trillions of microorganisms, collectively termed the microbiota, which play a fundamental role in regulating tissue physiology and immunity. Recent studies have sought to dissect the cellular and molecular mechanisms mediating communication between the microbiota and host immune system. Epithelial cells line the intestine and form an initial barrier separating the microbiota from underlying immune cells, and disruption of epithelial function has been associated with various conditions ranging from infection to inflammatory bowel diseases and cancer. From several studies, it is now clear that epithelial cells integrate signals from commensal microbes. Importantly, these non-hematopoietic cells also direct regulatory mechanisms that instruct the recruitment and function of microbiota-sensitive immune cells. In this review, we discuss the central role that has emerged for epithelial cells in orchestrating intestinal immunity and highlight epithelial pathways through which the microbiota can calibrate tissue-intrinsic immune responses.

Colorectal Cancer-Associated Myofibroblasts Exhibit Enhanced Angiogenin Expression and Signaling via the PLXNB2 Receptor

INTRODUCTION

Dynamic cell-cell interactions shape the tumor microenvironment to regulate tumor growth and invasiveness. Myofibroblasts are gastrointestinal stromal cells that are upregulated in the setting of colorectal cancer (CRC) and may play an important role in tumor-stromal cell communication. Angiogenin is a 14-kDa ribonuclease that regulates myofibroblast function and has been implicated in myofibroblast-CRC cell communication in mouse models. However, its role in human patients has not been well established.

METHODS

Open access, annotated single-cell RNA sequencing data of paired normal human colon and CRC tissue were available in the National Center for Biotechnology Information Gene Expression Omnibus Database. We supplemented and verified these data by analyzing scRNA-seq data from an independent set of paired normal human colon and CRC tissue. CellChat was used to quantitatively infer biologically meaningful cell-cell communication networks from scRNA-seq data. PLXNB2 and α-2 actin (ACTA2) are cell surface angiogenin receptors that regulate angiogenin signaling. Ligand-receptor interactions involving angiogenin, PLXNB2, and ACTA2 were analyzed between cell populations in each sample.

RESULTS

We found no difference in overall angiogenin expression comparing normal colon and CRC tissue. In normal colon tissue, myofibroblasts do not express angiogenin or the PLXNB2 receptor. In the presence of CRC, there was a striking increase in the number of myofibroblast cells within the surrounding stroma. CRC-associated myofibroblasts were characterized by a significant upregulation of both angiogenin and PLXNB2 receptor expression (P < 0.05), while no difference was seen in ACTA2. CRC cells not only use angiogenin for autocrine signaling but also communicate with myofibroblasts via the PLXNB2 receptor.

CONCLUSIONS

Compared to normal human colon tissue, CRC tissue is associated with an enrichment of myofibroblasts that exhibit upregulated expression of angiogenin and the angiogenin receptor PLXNB2. CRC cells engage in autocrine signaling via angiogenin and paracrine signaling with myofibroblasts via PLXNB2. Angiogenin appears to be directly involved in tumor-stromal cell communication in human CRC tissue and may play an important role in disease progression.

Quantitative Multiplexed Analysis of Indoleamine 2,3-Dioxygenase (IDO) and Arginase-1 (ARG1) Expression and Myeloid Cell Infiltration in Colorectal Cancer

Indoleamine 2,3-dioxygenase (IDO) and arginase-1 (ARG1) are amino acid-metabolizing enzymes, frequently highly expressed in cancer. Their expression may deplete essential amino acids, lead to immunosuppression, and promote cancer growth. Still, their expression patterns, prognostic significance, and spatial localization in the colorectal cancer microenvironment are incompletely understood. Using a custom 10-plex immunohistochemistry assay and supervised machine learning-based digital image analysis, we characterized IDO and ARG1 expression in monocytic cells, granulocytes, mast cells, and tumor cells in 833 colorectal cancer patients. We evaluated the prognostic value and spatial arrangement of IDO- and ARG1-expressing myeloid and tumor cells. IDO was mainly expressed not only by monocytic cells but also by some tumor cells, whereas ARG1 was predominantly expressed by granulocytes. Higher density of IDO+ monocytic cells was an independent prognostic factor for improved cancer-specific survival both in the tumor center (Ptrend = .0002; hazard ratio [HR] for the highest ordinal category Q4 [vs Q1], 0.51; 95% CI, 0.33-0.79) and the invasive margin (Ptrend = .0015). Higher density of granulocytes was associated with prolonged cancer-specific survival in univariable models, and higher FCGR3+ARG1+ neutrophil density in the tumor center also in multivariable analysis (Ptrend = .0020). Granulocytes were, on average, located closer to tumor cells than monocytic cells. Furthermore, IDO+ monocytic cells and ARG1- granulocytes were closer than IDO- monocytic cells and ARG1+ granulocytes, respectively. The mRNA expression of the IDO1 gene was assessed in myeloid and tumor cells using publicly available single-cell RNA sequencing data for 62 colorectal cancers. IDO1 was mainly expressed in monocytes and dendritic cells, and high IDO1 activity in monocytes was associated with enriched immunostimulatory pathways. Our findings provided in-depth information about the infiltration patterns and prognostic value of cells expressing IDO and/or ARG1 in the colorectal cancer microenvironment, highlighting the significance of host immune response in tumor progression.

Cellular senescence is associated with the spatial evolution toward a higher metastatic phenotype in colorectal cancer

In this study, we explore the dynamic process of colorectal cancer progression, emphasizing the evolution toward a more metastatic phenotype. The term "evolution" as used in this study specifically denotes the phenotypic transition toward a higher metastatic potency from well-formed glandular structures to collective invasion, ultimately resulting in the development of cancer cell buddings at the invasive front. Our findings highlight the spatial correlation of this evolution with tumor cell senescence, revealing distinct types of senescent tumor cells (types I and II) that play different roles in the overall cancer progression. Type I senescent tumor cells (p16INK4A+/CXCL12+/LAMC2-/MMP7-) are identified in the collective invasion region, whereas type II senescent tumor cells (p16INK4A+/CXCL12+/LAMC2+/MMP7+), representing the final evolved form, are prominently located in the partial-EMT region. Importantly, type II senescent tumor cells associate with local invasion and lymph node metastasis in colorectal cancer, potentially affecting patient prognosis.

The CLCF1-CNTFR axis drives an immunosuppressive tumor microenvironment and blockade enhances the effects of established cancer therapies

Tumors comprise a complex ecosystem consisting of many cell types that communicate through secreted factors. Targeting these intercellular signaling networks remains an important challenge in cancer research. Cardiotrophin-like cytokine factor 1 (CLCF1) is an interleukin-6 (IL-6) family member secreted by cancer-associated fibroblasts (CAFs) that binds to ciliary neurotrophic factor receptor (CNTFR), promoting tumor growth in lung and liver cancer1,2. A high-affinity soluble receptor (eCNTFR-Fc) that sequesters CLCF1 has anti-oncogenic effects3. However, the role of CLCF1 in mediating cell-cell interactions in cancer has remained unclear. We demonstrate that eCNTFR-Fc has widespread effects on both tumor cells and the tumor microenvironment and can sensitize cancer cells to KRAS inhibitors or immune checkpoint blockade. After three weeks of treatment with eCNTFR-Fc, there is a shift from an immunosuppressive to an immunostimulatory macrophage phenotype as well as an increase in activated T, NKT, and NK cells. Combination of eCNTFR-Fc and αPD1 was significantly more effective than single-agent therapy in a syngeneic allograft model, and eCNTFR-Fc sensitizes tumor cells to αPD1 in a non-responsive GEM model of lung adenocarcinoma. These data suggest that combining eCNTFR-Fc with KRAS inhibition or with αPD1 is a novel therapeutic strategy for lung cancer and potentially other cancers in which these therapies have been used but to date with only modest effect. Overall, we demonstrate the potential of cancer therapies that target cytokines to alter the immune microenvironment.

Unveiling spatial complexity in solid tumor immune microenvironments through multiplexed imaging

Deciphering cellular components and the spatial interaction network of the tumor immune microenvironment (TIME) of solid tumors is pivotal for understanding biologically relevant cross-talks and, ultimately, advancing therapies. Multiplexed tissue imaging provides a powerful tool to elucidate spatial complexity in a holistic manner. We established and cross-validated a comprehensive immunophenotyping panel comprising over 121 markers for multiplexed tissue imaging using MACSima™ imaging cyclic staining (MICS) alongside an end-to-end analysis workflow. Applying this panel and workflow to primary cancer tissues, we characterized tumor heterogeneity, investigated potential therapeutical targets, conducted in-depth profiling of cell types and states, sub-phenotyped T cells within the TIME, and scrutinized cellular neighborhoods of diverse T cell subsets. Our findings highlight the advantage of spatial profiling, revealing immunosuppressive molecular signatures of tumor-associated myeloid cells interacting with neighboring exhausted, PD1high T cells in the TIME of hepatocellular carcinoma (HCC). This study establishes a robust framework for spatial exploration of TIMEs in solid tumors and underscores the potency of multiplexed tissue imaging and ultra-deep cell phenotyping in unraveling clinically relevant tumor components.

An immune cell map of human lung adenocarcinoma development reveals an anti-tumoral role of the Tfh-dependent tertiary lymphoid structure

The immune responses during the initiation and invasion stages of human lung adenocarcinoma (LUAD) development are largely unknown. Here, we generated a single-cell RNA sequencing map to decipher the immune dynamics during human LUAD development. We found that T follicular helper (Tfh)-like cells, germinal center B cells, and dysfunctional CD8+ T cells increase during tumor initiation/invasion and form a tertiary lymphoid structure (TLS) inside the tumor. This TLS starts with an aggregation of CD4+ T cells and the generation of CXCL13-expressing Tfh-like cells, followed by an accumulation of B cells, and then forms a CD4+ T and B cell aggregate. TLS and its associated cells are correlated with better patient survival. Inhibiting TLS formation by Tfh or B cell depletion promotes tumor growth in mouse models. The anti-tumoral effect of the Tfh-dependent TLS is mediated through interleukin-21 (IL-21)-IL-21 receptor signaling. Our study establishes an anti-tumoral role of the Tfh-dependent TLS in the development of LUAD.

Highly sensitive spatial transcriptomics using FISHnCHIPs of multiple co-expressed genes

High-dimensional, spatially resolved analysis of intact tissue samples promises to transform biomedical research and diagnostics, but existing spatial omics technologies are costly and labor-intensive. We present Fluorescence In Situ Hybridization of Cellular HeterogeneIty and gene expression Programs (FISHnCHIPs) for highly sensitive in situ profiling of cell types and gene expression programs. FISHnCHIPs achieves this by simultaneously imaging ~2-35 co-expressed genes (clustered into modules) that are spatially co-localized in tissues, resulting in similar spatial information as single-gene Fluorescence In Situ Hybridization (FISH), but with ~2-20-fold higher sensitivity. Using FISHnCHIPs, we image up to 53 modules from the mouse kidney and mouse brain, and demonstrate high-speed, large field-of-view profiling of a whole tissue section. FISHnCHIPs also reveals spatially restricted localizations of cancer-associated fibroblasts in a human colorectal cancer biopsy. Overall, FISHnCHIPs enables fast, robust, and scalable cell typing of tissues with normal physiology or undergoing pathogenesis.

Comparative analysis of single-cell transcriptome reveals heterogeneity and commonality in the immune microenvironment of colorectal cancer and inflammatory bowel disease

Background

During aging, chronic inflammation can promote tumor development and metastasis. Patients with chronic inflammatory bowel diseases (IBD) are at an increased risk of developing colorectal cancer (CRC). However, the molecular mechanism underlying is still unclear.

Methods

We conducted a large-scale single-cell sequencing analysis comprising 432,314 single cells from 92 CRC and 24 IBD patients. The analysis focused on the heterogeneity and commonality of CRC and IBD with respect to immune cell landscape, cellular communication, aging and inflammatory response, and Meta programs.

Results

The CRC and IBD had significantly different propensities in terms of cell proportions, differential genes and their functions, and cellular communication. The progression of CRC was mainly associated with epithelial cells, fibroblasts, and monocyte-macrophages, which displayed pronounced metabolic functions. In particular, monocyte-macrophages were enriched for the aging and inflammation-associated NF-κB pathway. And IBD was enriched in immune-related functions with B cells and T cells. Cellular communication analysis in CRC samples displayed an increase in MIF signaling from epithelial cells to T cells, and an increase in the efferent signal of senescence-associated SPP1 signaling from monocyte-macrophages. Notably, we also found some commonalities between CRC and IBD. The efferent and afferent signals showed that the pro-inflammatory cytokine played an important role. And the activity of aging and inflammatory response with AUCell analysis also showed a high degree of commonality. Furthermore, using the Meta programs (MPs) with the NMF algorithm, we found that the CRC non-malignant cells shared a substantial proportion of the MP genes with CRC malignant cells (68% overlap) and IBD epithelial cells (52% overlap), respectively. And it was extensively involved in functions of cell cycle and immune response, revealing its dual properties of inflammation and cancer. In addition, CRC malignant and non-malignant cells were enriched for the senescence-related cell cycle G2M phase transition and the p53 signaling pathway.

Conclusion

Our study highlights the characteristics of aging, inflammation and tumor in CRC and IBD at the single-cell level, and the dual property of inflammation-cancer in CRC non-malignant cells may provide a more up-to-date understanding of disease transformation.

CDSKNNXMBD: a novel clustering framework for large-scale single-cell data based on a stable graph structure

BACKGROUND

Accurate and efficient cell grouping is essential for analyzing single-cell transcriptome sequencing (scRNA-seq) data. However, the existing clustering techniques often struggle to provide timely and accurate cell type groupings when dealing with datasets with large-scale or imbalanced cell types. Therefore, there is a need for improved methods that can handle the increasing size of scRNA-seq datasets while maintaining high accuracy and efficiency.

METHODS

We propose CDSKNNXMBD (Community Detection based on a Stable K-Nearest Neighbor Graph Structure), a novel single-cell clustering framework integrating partition clustering algorithm and community detection algorithm, which achieves accurate and fast cell type grouping by finding a stable graph structure.

RESULTS

We evaluated the effectiveness of our approach by analyzing 15 tissues from the human fetal atlas. Compared to existing methods, CDSKNN effectively counteracts the high imbalance in single-cell data, enabling effective clustering. Furthermore, we conducted comparisons across multiple single-cell datasets from different studies and sequencing techniques. CDSKNN is of high applicability and robustness, and capable of balancing the complexities of across diverse types of data. Most importantly, CDSKNN exhibits higher operational efficiency on datasets at the million-cell scale, requiring an average of only 6.33 min for clustering 1.46 million single cells, saving 33.3% to 99% of running time compared to those of existing methods.

CONCLUSIONS

The CDSKNN is a flexible, resilient, and promising clustering tool that is particularly suitable for clustering imbalanced data and demonstrates high efficiency on large-scale scRNA-seq datasets.

SOX17 enables immune evasion of early colorectal adenomas and cancers

A hallmark of cancer is the avoidance of immune destruction. This process has been primarily investigated in locally advanced or metastatic cancer1-3; however, much less is known about how pre-malignant or early invasive tumours evade immune detection. Here, to understand this process in early colorectal cancers (CRCs), we investigated how naive colon cancer organoids that were engineered in vitro to harbour Apc-null, KrasG12D and Trp53-null (AKP) mutations adapted to the in vivo native colonic environment. Comprehensive transcriptomic and chromatin analyses revealed that the endoderm-specifying transcription factor SOX17 became strongly upregulated in vivo. Notably, whereas SOX17 loss did not affect AKP organoid propagation in vitro, its loss markedly reduced the ability of AKP tumours to persist in vivo. The small fraction of SOX17-null tumours that grew displayed notable interferon-γ (IFNγ)-producing effector-like CD8+ T cell infiltrates in contrast to the immune-suppressive microenvironment in wild-type counterparts. Mechanistically, in both endogenous Apc-null pre-malignant adenomas and transplanted organoid-derived AKP CRCs, SOX17 suppresses the ability of tumour cells to sense and respond to IFNγ, preventing anti-tumour T cell responses. Finally, SOX17 engages a fetal intestinal programme that drives differentiation away from LGR5+ tumour cells to produce immune-evasive LGR5- tumour cells with lower expression of major histocompatibility complex class I (MHC-I). We propose that SOX17 is a transcription factor that is engaged during the early steps of colon cancer to orchestrate an immune-evasive programme that permits CRC initiation and progression.

Sirpα on tumor-associated myeloid cells restrains antitumor immunity in colorectal cancer independent of its interaction with CD47

Immunosuppressive myeloid cells hinder immunotherapeutic efficacy in tumors, but the precise mechanisms remain undefined. Here, by performing single-cell RNA sequencing in colorectal cancer tissues, we found tumor-associated macrophages and granulocytic myeloid-derived suppressor cells increased most compared to their counterparts in normal tissue and displayed the highest immune-inhibitory signatures among all immunocytes. These cells exhibited significantly increased expression of immunoreceptor tyrosine-based inhibitory motif-bearing receptors, including SIRPA. Notably, Sirpa-/- mice were more resistant to tumor progression than wild-type mice. Moreover, Sirpα deficiency reprogramed the tumor microenvironment through expansion of TAM_Ccl8hi and gMDSC_H2-Q10hi subsets showing strong antitumor activity. Sirpa-/- macrophages presented strong phagocytosis and antigen presentation to enhance T cell activation and proliferation. Furthermore, Sirpa-/- macrophages facilitated T cell recruitment via Syk/Btk-dependent Ccl8 secretion. Therefore, Sirpα deficiency enhances innate and adaptive immune activation independent of expression of CD47 and Sirpα blockade could be a promising strategy to improve cancer immunotherapy efficacy.

BANKSY unifies cell typing and tissue domain segmentation for scalable spatial omics data analysis

Spatial omics data are clustered to define both cell types and tissue domains. We present Building Aggregates with a Neighborhood Kernel and Spatial Yardstick (BANKSY), an algorithm that unifies these two spatial clustering problems by embedding cells in a product space of their own and the local neighborhood transcriptome, representing cell state and microenvironment, respectively. BANKSY's spatial feature augmentation strategy improved performance on both tasks when tested on diverse RNA (imaging, sequencing) and protein (imaging) datasets. BANKSY revealed unexpected niche-dependent cell states in the mouse brain and outperformed competing methods on domain segmentation and cell typing benchmarks. BANKSY can also be used for quality control of spatial transcriptomics data and for spatially aware batch effect correction. Importantly, it is substantially faster and more scalable than existing methods, enabling the processing of millions of cell datasets. In summary, BANKSY provides an accurate, biologically motivated, scalable and versatile framework for analyzing spatially resolved omics data.

Targeted deletion of CD244 on monocytes promotes differentiation into anti-tumorigenic macrophages and potentiates PD-L1 blockade in melanoma

BACKGROUND

In the myeloid compartment of the tumor microenvironment, CD244 signaling has been implicated in immunosuppressive phenotype of monocytes. However, the precise molecular mechanism and contribution of CD244 to tumor immunity in monocytes/macrophages remains elusive due to the co-existing lymphoid cells expressing CD244.

METHODS

To directly assess the role of CD244 in tumor-associated macrophages, monocyte-lineage-specific CD244-deficient mice were generated using cre-lox recombination and challenged with B16F10 melanoma. The phenotype and function of tumor-infiltrating macrophages along with antigen-specific CD8 T cells were analyzed by flow cytometry and single cell RNA sequencing data analysis, and the molecular mechanism underlying anti-tumorigenic macrophage differentiation, antigen presentation, phagocytosis was investigated ex vivo. Finally, the clinical feasibility of CD244-negative monocytes as a therapeutic modality in melanoma was confirmed by adoptive transfer experiments.

RESULTS

CD244fl/flLysMcre mice demonstrated a significant reduction in tumor volume (61% relative to that of the CD244fl/fl control group) 14 days after tumor implantation. Within tumor mass, CD244fl/flLysMcre mice also showed higher percentages of Ly6Clow macrophages, along with elevated gp100+IFN-γ+ CD8 T cells. Flow cytometry and RNA sequencing data demonstrated that ER stress resulted in increased CD244 expression on monocytes. This, in turn, impeded the generation of anti-tumorigenic Ly6Clow macrophages, phagocytosis and MHC-I antigen presentation by suppressing autophagy pathways. Combining anti-PD-L1 antibody with CD244-/- bone marrow-derived macrophages markedly improved tumor rejection compared to the anti-PD-L1 antibody alone or in combination with wild-type macrophages. Consistent with the murine data, transcriptome analysis of human melanoma tissue single-cell RNA-sequencing dataset revealed close association between CD244 and the inhibition of macrophage maturation and function. Furthermore, the presence of CD244-negative monocytes/macrophages significantly increased patient survival in primary and metastatic tumors.

CONCLUSION

Our study highlights the novel role of CD244 on monocytes/macrophages in restraining anti-tumorigenic macrophage generation and tumor antigen-specific T cell response in melanoma. Importantly, our findings suggest that CD244-deficient macrophages could potentially be used as a therapeutic agent in combination with immune checkpoint inhibitors. Furthermore, CD244 expression in monocyte-lineage cells serve as a prognostic marker in cancer patients.

Targeting the activated microenvironment with endosialin (CD248)-directed CAR-T cells ablates perivascular cells to impair tumor growth and metastasis

BACKGROUND

Targeting of solid cancers with chimeric antigen receptor (CAR)-T cells is limited by the lack of suitable tumor-specific antigens and the immunosuppressive, desmoplastic tumor microenvironment that impedes CAR-T cell infiltration, activity and persistence. We hypothesized that targeting the endosialin (CD248) receptor, strongly expressed by tumor-associated pericytes and perivascular cancer-associated fibroblasts, would circumvent these challenges and offer an exciting antigen for CAR-T cell therapy due to the close proximity of target cells to the tumor vasculature, the limited endosialin expression in normal tissues and the lack of phenotype observed in endosialin knockout mice.

METHODS

We generated endosialin-directed E3K CAR-T cells from three immunocompetent mouse strains, BALB/c, FVB/N and C57BL/6. E3K CAR-T cell composition (CD4+/CD8+ ratio), activity in vitro against endosialin+ and endosialin- cells, and expansion and activity in vivo in syngeneic tumor models as well as in tumor-naive healthy and wounded mice and tumor-bearing endosialin knockout mice was assessed.

RESULTS

E3K CAR-T cells were active in vitro against both mouse and human endosialin+, but not endosialin-, cells. Adoptively transferred E3K CAR-T cells exhibited no activity in endosialin knockout mice, tumor-naive endosialin wildtype mice or in wound healing models, demonstrating an absence of off-target and on-target/off-tumor activity. By contrast, adoptive transfer of E3K CAR-T cells into BALB/c, FVB/N or C57BL/6 mice bearing syngeneic breast or lung cancer lines depleted target cells in the tumor stroma resulting in increased tumor necrosis, reduced tumor growth and a substantial impairment in metastatic outgrowth.

CONCLUSIONS

Together these data highlight endosialin as a viable antigen for CAR-T cell therapy and that targeting stromal cells closely associated with the tumor vasculature avoids CAR-T cells having to navigate the harsh immunosuppressive tumor microenvironment. Further, the ability of E3K CAR-T cells to recognize and target both mouse and human endosialin+ cells makes a humanized and optimized E3K CAR a promising candidate for clinical development applicable to a broad range of solid tumor types.

Identification and validation of a pyroptosis-related prognostic model for colorectal cancer based on bulk and single-cell RNA sequencing data

BACKGROUND

Pyroptosis impacts the development of malignant tumors, yet its role in colorectal cancer (CRC) prognosis remains uncertain.

AIM

To assess the prognostic significance of pyroptosis-related genes and their association with CRC immune infiltration.

METHODS

Gene expression data were obtained from The Cancer Genome Atlas (TCGA) and single-cell RNA sequencing dataset GSE178341 from the Gene Expression Omnibus (GEO). Pyroptosis-related gene expression in cell clusters was analyzed, and enrichment analysis was conducted. A pyroptosis-related risk model was developed using the LASSO regression algorithm, with prediction accuracy assessed through K-M and receiver operating characteristic analyses. A nomogram predicting survival was created, and the correlation between the risk model and immune infiltration was analyzed using CIBERSORTx calculations. Finally, the differential expression of the 8 prognostic genes between CRC and normal samples was verified by analyzing TCGA-COADREAD data from the UCSC database.

RESULTS

An effective pyroptosis-related risk model was constructed using 8 genes-CHMP2B, SDHB, BST2, UBE2D2, GJA1, AIM2, PDCD6IP, and SEZ6L2 (P < 0.05). Seven of these genes exhibited differential expression between CRC and normal samples based on TCGA database analysis (P < 0.05). Patients with higher risk scores demonstrated increased death risk and reduced overall survival (P < 0.05). Significant differences in immune infiltration were observed between low- and high-risk groups, correlating with pyroptosis-related gene expression.

CONCLUSION

We developed a pyroptosis-related prognostic model for CRC, affirming its correlation with immune infiltration. This model may prove useful for CRC prognostic evaluation.

Identifying CTH and MAP1LC3B as ferroptosis biomarkers for prognostic indication in gastric cancer decoding

Gastric cancer (GC), known for its high incidence and poor prognosis, urgently necessitates the identification of reliable prognostic biomarkers to enhance patient outcomes. We scrutinized data from 375 GC patients alongside 32 non-cancer controls, sourced from the TCGA database. A univariate Cox Proportional Hazards Model (COX) regression was employed to evaluate expressions of ferroptosis-related genes. This was followed by the application of Least Absolute Shrinkage and Selection Operator (LASSO) and multivariate COX regression for the development of prognostic models. The composition of immune cell subtypes was quantified utilizing CIBERSORT, with their distribution in GC versus control samples being comparatively analyzed. Furthermore, the correlation between the expressions of Cystathionine Gamma-Lyase (CTH) and Microtubule Associated Protein 1 Light Chain 3 Beta (MAP1LC3B) and the abundance of immune cell subtypes was explored. Our bioinformatics findings underwent validation through immunohistochemical analysis. Our prognostic models integrated CTH and MAP1LC3B. Survival analysis indicated that patients categorized as high-risk, as defined by the model, exhibited significantly lower survival rates compared to their low-risk counterparts. Notably, CTH expression inversely correlated with monocyte levels, while MAP1LC3B expression showed an inverse relationship with the abundance of M2 macrophages. Immunohistochemical validation corroborated lower expressions of CTH and MAP1LC3B in GC tissues relative to control samples, in concordance with our bioinformatics predictions. Our study suggests that the dysregulation of CTH, MAP1LC3B, and the accompanying monocyte-macrophage dynamics could be pivotal in the prognosis of GC. These elements present potential targets for prognostic assessment and therapeutic intervention.

Mapping the core senescence phenotype of primary human colon fibroblasts

Advanced age is the largest risk factor for many diseases and several types of cancer, including colorectal cancer (CRC). Senescent cells are known to accumulate with age in various tissues, where they can modulate the surrounding tissue microenvironment through their senescence associated secretory phenotype (SASP). Recently, we showed that there is an increased number of senescent cells in the colons of CRC patients and demonstrated that senescent fibroblasts and their SASP create microniches in the colon that are conducive to CRC onset and progression. However, the composition of the SASP is heterogenous and cell-specific, and the precise senescence profile of colon fibroblasts has not been well-defined. To generate a SASP atlas of human colon fibroblasts, we induced senescence in primary human colon fibroblasts using various in vitro methods and assessed the resulting transcriptome. Using RNASequencing and further validation by quantitative RT-PCR and Luminex assays, we define and validate a 'core senescent profile' that might play a significant role in shaping the colon microenvironment. We also performed KEGG analysis and GO analyses to identify key pathways and biological processes that are differentially regulated in colon fibroblast senescence. These studies provide insights into potential driver proteins involved in senescence-associated diseases, like CRC, which may lead to therapies to improve overall health in the elderly and to prevent CRC.

Clinical and biomarker results from a phase II trial of combined cabozantinib and durvalumab in patients with chemotherapy-refractory colorectal cancer (CRC): CAMILLA CRC cohort

CAMILLA is a basket trial (NCT03539822) evaluating cabozantinib plus the ICI durvalumab in chemorefractory gastrointestinal cancer. Herein, are the phase II colorectal cohort results. 29 patients were evaluable. 100% had confirmed pMMR/MSS tumors. Primary endpoint was met with ORR of 27.6% (95% CI 12.7-47.2%). Secondary endpoints of 4-month PFS rate was 44.83% (95% CI 26.5-64.3%); and median OS was 9.1 months (95% CI 5.8-20.2). Grade≥3 TRAE occurred in 39%. In post-hoc analysis of patients with RAS wild type tumors, ORR was 50% and median PFS and OS were 6.3 and 21.5 months respectively. Exploratory spatial transcriptomic profiling of pretreatment tumors showed upregulation of VEGF and MET signaling, increased extracellular matrix activity and preexisting anti-tumor immune responses coexisting with immune suppressive features like T cell migration barriers in responders versus non-responders. Cabozantinib plus durvalumab demonstrated anti-tumor activity, manageable toxicity, and have led to the activation of the phase III STELLAR-303 trial.

A multiomics analysis-assisted deep learning model identifies a macrophage-oriented module as a potential therapeutic target in colorectal cancer

Colorectal cancer (CRC) is a common malignancy involving multiple cellular components. The CRC tumor microenvironment (TME) has been characterized well at single-cell resolution. However, a spatial interaction map of the CRC TME is still elusive. Here, we integrate multiomics analyses and establish a spatial interaction map to improve the prognosis, prediction, and therapeutic development for CRC. We construct a CRC immune module (CCIM) that comprises FOLR2+ macrophages, exhausted CD8+ T cells, tolerant CD8+ T cells, exhausted CD4+ T cells, and regulatory T cells. Multiplex immunohistochemistry is performed to depict the CCIM. Based on this, we utilize advanced deep learning technology to establish a spatial interaction map and predict chemotherapy response. CCIM-Net is constructed, which demonstrates good predictive performance for chemotherapy response in both the training and testing cohorts. Lastly, targeting FOLR2+ macrophage therapeutics is used to disrupt the immunosuppressive CCIM and enhance the chemotherapy response in vivo.