Spatially organized multicellular immune hubs in human colorectal cancer

SPACEL: deep learning-based characterization of spatial transcriptome architectures

Spatial transcriptomics (ST) technologies detect mRNA expression in single cells/spots while preserving their two-dimensional (2D) spatial coordinates, allowing researchers to study the spatial distribution of the transcriptome in tissues; however, joint analysis of multiple ST slices and aligning them to construct a three-dimensional (3D) stack of the tissue still remain a challenge. Here, we introduce spatial architecture characterization by deep learning (SPACEL) for ST data analysis. SPACEL comprises three modules: Spoint embeds a multiple-layer perceptron with a probabilistic model to deconvolute cell type composition for each spot in a single ST slice; Splane employs a graph convolutional network approach and an adversarial learning algorithm to identify spatial domains that are transcriptomically and spatially coherent across multiple ST slices; and Scube automatically transforms the spatial coordinate systems of consecutive slices and stacks them together to construct a 3D architecture of the tissue. Comparisons against 19 state-of-the-art methods using both simulated and real ST datasets from various tissues and ST technologies demonstrate that SPACEL outperforms the others for cell type deconvolution, for spatial domain identification, and for 3D alignment, thus showcasing SPACEL as a valuable integrated toolkit for ST data processing and analysis.

Single cell dynamics of tumor specificity vs bystander activity in CD8+ T cells define the diverse immune landscapes in colorectal cancer

CD8+ T cell activation via immune checkpoint blockade (ICB) is successful in microsatellite instable (MSI) colorectal cancer (CRC) patients. By comparison, the success of immunotherapy against microsatellite stable (MSS) CRC is limited. Little is known about the most critical features of CRC CD8+ T cells that together determine the diverse immune landscapes and contrasting ICB responses. Hence, we pursued a deep single cell mapping of CRC CD8+ T cells on transcriptomic and T cell receptor (TCR) repertoire levels in a diverse patient cohort, with additional surface proteome validation. This revealed that CRC CD8+ T cell dynamics are underscored by complex interactions between interferon-γ signaling, tumor reactivity, TCR repertoire, (predicted) TCR antigen-specificities, and environmental cues like gut microbiome or colon tissue-specific 'self-like' features. MSI CRC CD8+ T cells showed tumor-specific activation reminiscent of canonical 'T cell hot' tumors, whereas the MSS CRC CD8+ T cells exhibited tumor unspecific or bystander-like features. This was accompanied by inflammation reminiscent of 'pseudo-T cell hot' tumors. Consequently, MSI and MSS CRC CD8+ T cells showed overlapping phenotypic features that differed dramatically in their TCR antigen-specificities. Given their high discriminating potential for CD8+ T cell features/specificities, we used the single cell tumor-reactive signaling modules in CD8+ T cells to build a bulk tumor transcriptome classification for CRC patients. This "Immune Subtype Classification" (ISC) successfully distinguished various tumoral immune landscapes that showed prognostic value and predicted immunotherapy responses in CRC patients. Thus, we deliver a unique map of CRC CD8+ T cells that drives a novel tumor immune landscape classification, with relevance for immunotherapy decision-making.

Intercellular Molecular Crosstalk Networks within Invasive and Immunosuppressive Tumor Microenvironment Subtypes Associated with Clinical Outcomes in Four Cancer Types

Heterogeneity is a critical basis for understanding how the tumor microenvironment (TME) contributes to tumor progression. However, an understanding of the specific characteristics and functions of TME subtypes (subTMEs) in the progression of cancer is required for further investigations into single-cell resolutions. Here, we analyzed single-cell RNA sequencing data of 250 clinical samples with more than 200,000 cells analyzed in each cancer datum. Based on the construction of an intercellular infiltration model and unsupervised clustering analysis, four, three, three, and four subTMEs were revealed in breast, colorectal, esophageal, and pancreatic cancer, respectively. Among the subTMEs, the immune-suppressive subTME (subTME-IS) and matrix remodeling with malignant cells subTME (subTME-MRM) were highly enriched in tumors, whereas the immune cell infiltration subTME (subTME-ICI) and precancerous state of epithelial cells subTME (subTME-PSE) were less in tumors, compared with paracancerous tissues. We detected and compared genes encoding cytokines, chemokines, cytotoxic mediators, PD1, and PD-L1. The results showed that these genes were specifically overexpressed in different cell types, and, compared with normal tissues, they were upregulated in tumor-derived cells. In addition, compared with other subTMEs, the expression levels of PDCD1 and TGFB1 were higher in subTME-IS. The Cox proportional risk regression model was further constructed to identify possible prognostic markers in each subTME across four cancer types. Cell-cell interaction analysis revealed the distinguishing features in molecular pairs among different subTMEs. Notably, ligand-receptor gene pairs, including COL1A1-SDC1, COL6A2-SDC1, COL6A3-SDC1, and COL4A1-ITGA2 between stromal and tumor cells, associated with tumor invasion phenotypes, poor patient prognoses, and tumor advanced progression, were revealed in subTME-MRM. C5AR1-RPS19, LGALS9-HAVCR2, and SPP1-PTGER4 between macrophages and CD8+ T cells, associated with CD8+ T-cell dysfunction, immunosuppressive status, and tumor advanced progression, were revealed in subTME-IS. The spatial co-location information of cellular and molecular interactions was further verified by spatial transcriptome data from colorectal cancer clinical samples. Overall, our study revealed the heterogeneity within the TME, highlighting the potential pro-invasion and pro-immunosuppressive functions and cellular infiltration characteristics of specific subTMEs, and also identified the key cellular and molecular interactions that might be associated with the survival, invasion, immune escape, and classification of cancer patients across four cancer types.

Molecular portraits of colorectal cancer morphological regions

Heterogeneity of colorectal carcinoma (CRC) represents a major hurdle towards personalized medicine. Efforts based on whole tumor profiling demonstrated that the CRC molecular subtypes were associated with specific tumor morphological patterns representing tumor subregions. We hypothesize that whole-tumor molecular descriptors depend on the morphological heterogeneity with significant impact on current molecular predictors. We investigated intra-tumor heterogeneity by morphology-guided transcriptomics to better understand the links between gene expression and tumor morphology represented by six morphological patterns (morphotypes): complex tubular, desmoplastic, mucinous, papillary, serrated, and solid/trabecular. Whole-transcriptome profiling by microarrays of 202 tumor regions (morphotypes, tumor-adjacent normal tissue, supportive stroma, and matched whole tumors) from 111 stage II-IV CRCs identified morphotype-specific gene expression profiles and molecular programs and differences in their cellular buildup. The proportion of cell types (fibroblasts, epithelial and immune cells) and differentiation of epithelial cells were the main drivers of the observed disparities with activation of EMT and TNF-α signaling in contrast to MYC and E2F targets signaling, defining major gradients of changes at molecular level. Several gene expression-based (including single-cell) classifiers, prognostic and predictive signatures were examined to study their behavior across morphotypes. Most exhibited important morphotype-dependent variability within same tumor sections, with regional predictions often contradicting the whole-tumor classification. The results show that morphotype-based tumor sampling allows the detection of molecular features that would otherwise be distilled in whole tumor profile, while maintaining histopathology context for their interpretation. This represents a practical approach at improving the reproducibility of expression profiling and, by consequence, of gene-based classifiers.

Concepts of multi-level dynamical modelling: understanding mechanisms of squamous cell carcinoma development in Fanconi anemia

Fanconi anemia (FA) is a rare disease (incidence of 1:300,000) primarily based on the inheritance of pathogenic variants in genes of the FA/BRCA (breast cancer) pathway. These variants ultimately reduce the functionality of different proteins involved in the repair of DNA interstrand crosslinks and DNA double-strand breaks. At birth, individuals with FA might present with typical malformations, particularly radial axis and renal malformations, as well as other physical abnormalities like skin pigmentation anomalies. During the first decade of life, FA mostly causes bone marrow failure due to reduced capacity and loss of the hematopoietic stem and progenitor cells. This often makes hematopoietic stem cell transplantation necessary, but this therapy increases the already intrinsic risk of developing squamous cell carcinoma (SCC) in early adult age. Due to the underlying genetic defect in FA, classical chemo-radiation-based treatment protocols cannot be applied. Therefore, detecting and treating the multi-step tumorigenesis process of SCC in an early stage, or even its progenitors, is the best option for prolonging the life of adult FA individuals. However, the small number of FA individuals makes classical evidence-based medicine approaches based on results from randomized clinical trials impossible. As an alternative, we introduce here the concept of multi-level dynamical modelling using large, longitudinally collected genome, proteome- and transcriptome-wide data sets from a small number of FA individuals. This mechanistic modelling approach is based on the "hallmarks of cancer in FA", which we derive from our unique database of the clinical history of over 750 FA individuals. Multi-omic data from healthy and diseased tissue samples of FA individuals are to be used for training constituent models of a multi-level tumorigenesis model, which will then be used to make experimentally testable predictions. In this way, mechanistic models facilitate not only a descriptive but also a functional understanding of SCC in FA. This approach will provide the basis for detecting signatures of SCCs at early stages and their precursors so they can be efficiently treated or even prevented, leading to a better prognosis and quality of life for the FA individual.

Toward the clinical development of synthetic immunity to cancer

Synthetic biology (synbio) tools, such as chimeric antigen receptors (CARs), have been designed to target, activate, and improve immune cell responses to tumors. These therapies have demonstrated an ability to cure patients with blood cancers. However, there are significant challenges to designing, testing, and efficiently translating these complex cell therapies for patients who do not respond or have immune refractory solid tumors. The rapid progress of synbio tools for cell therapy, particularly for cancer immunotherapy, is encouraging but our development process should be tailored to increase translational success. Particularly, next-generation cell therapies should be rooted in basic immunology, tested in more predictive preclinical models, engineered for potency with the right balance of safety, educated by clinical findings, and multi-faceted to combat a range of suppressive mechanisms. Here, we lay out five principles for engineering future cell therapies to increase the probability of clinical impact, and in the context of these principles, we provide an overview of the current state of synbio cell therapy design for cancer. Although these principles are anchored in engineering immune cells for cancer therapy, we posit that they can help guide translational synbio research for broad impact in other disease indications with high unmet need.

Single-Cell Profiling Reveals Sustained Immune Infiltration, Surveillance, and Tumor Heterogeneity in Infiltrative Basal Cell Carcinoma

Infiltrative basal cell carcinoma (iBCC) is a particularly aggressive subtype of basal cell carcinoma that tends to progress and recur after surgery, and its malignancy is closely related to the tumor microenvironment. In this study, we performed a comprehensive single-cell RNA analysis to profile 29,334 cells from iBCC and adjacent normal skin. We found active immune collaborations enriched in iBCC. Specifically, SPP1+CXCL9/10high macrophage 1 had strong BAFF signaling with plasma cells, and T follicular helper-like cells highly expressed the B-cell chemokine CXCL13. Heterogeneous proinflammatory SPP1+CXCL9/10high macrophage 1 and angiogenesis-related SPP1+CCL2high macrophage 1 were identified within the tumor microenvironment. Interestingly, we found an upregulation of major histocompatibility complex I molecules in fibroblasts in iBCC compared with those in adjacent normal skin. Moreover, MDK signals derived from malignant basal cells were markedly increased, and their expression was an independent factor in predicting the infiltration depth of iBCC, emphasizing its role in driving malignancy and remodeling the tumor microenvironment. In addition, we identified differentiation-associated SOSTDC1+IGFBP5+CTSV+ malignant basal subtype 1 and epithelial-mesenchymal transition-associated TNC+SFRP1+CHGA+ malignant basal subtype 2 cells. The high expression of malignant basal 2 cell markers was associated with the invasion and recurrence of iBCC. Altogether, our study helps to elucidate the cellular heterogeneity in iBCC and provides potential therapeutic targets for clinical research.

Depletion of slow-cycling PDGFRα+ADAM12+ mesenchymal cells promotes antitumor immunity by restricting macrophage efferocytosis

The capacity to survive and thrive in conditions of limited resources and high inflammation is a major driver of tumor malignancy. Here we identified slow-cycling ADAM12+PDGFRα+ mesenchymal stromal cells (MSCs) induced at the tumor margins in mouse models of melanoma, pancreatic cancer and prostate cancer. Using inducible lineage tracing and transcriptomics, we demonstrated that metabolically altered ADAM12+ MSCs induced pathological angiogenesis and immunosuppression by promoting macrophage efferocytosis and polarization through overexpression of genes such as Gas6, Lgals3 and Csf1. Genetic depletion of ADAM12+ cells restored a functional tumor vasculature, reduced hypoxia and acidosis and normalized CAFs, inducing infiltration of effector T cells and growth inhibition of melanomas and pancreatic neuroendocrine cancer, in a process dependent on TGF-β. In human cancer, ADAM12 stratifies patients with high levels of hypoxia and innate resistance mechanisms, as well as factors associated with a poor prognosis and drug resistance such as AXL. Altogether, our data show that depletion of tumor-induced slow-cycling PDGFRα+ MSCs through ADAM12 restores antitumor immunity.

DNA mismatch repair deficient cancer - Emerging biomarkers of resistance to immune checkpoint inhibition

The DNA mismatch repair pathway is involved in the identification, excision, and repair of base-base mismatches and indel loops in the genome. Mismatch repair deficiency occurs in approximately 20% of all cancers and results in a type of DNA damage called microsatellite instability. In 2017, the immune checkpoint inhibitor, Pembrolizumab, an anti-PD-1 therapy, was approved for use in all unresectable or metastatic tumours that were mismatch repair deficient or had high microsatellite instability regardless of tissue origin. This landmark approval was the first time a drug had been approved in a site agnostic way, but accumulating data has revealed that up to 50% of mismatch repair deficient tumours are refractory to treatment and there is a huge amount of variability in the therapeutic benefit amongst responders. Several mechanisms of resistance to immune checkpoint blockade for mismatch repair deficient cancers have been identified but our understanding of what is driving resistance in a proportion of patients remains lacking. In this review article, we discuss the emerging mechanisms of resistance which may enable optimal stratification of patients for treatment with immune checkpoint inhibitors in the future.

Bridging live-cell imaging and next-generation cancer treatment

By providing spatial, molecular and morphological data over time, live-cell imaging can provide a deeper understanding of the cellular and signalling events that determine cancer response to treatment. Understanding this dynamic response has the potential to enhance clinical outcome by identifying biomarkers or actionable targets to improve therapeutic efficacy. Here, we review recent applications of live-cell imaging for uncovering both tumour heterogeneity in treatment response and the mode of action of cancer-targeting drugs. Given the increasing uses of T cell therapies, we discuss the unique opportunity of time-lapse imaging for capturing the interactivity and motility of immunotherapies. Although traditionally limited in the number of molecular features captured, novel developments in multidimensional imaging and multi-omics data integration offer strategies to connect single-cell dynamics to molecular phenotypes. We review the effect of these recent technological advances on our understanding of the cellular dynamics of tumour targeting and discuss their implication for next-generation precision medicine.

Microdissection of cancer-associated fibroblast infiltration subtypes unveils the secreted SERPINE2 contributing to immunosuppressive microenvironment and immuotherapeutic resistance in gastric cancer: A large-scale study integrating bulk and single-cell transcriptome profiling

In the era of immunotherapy, the suboptimal response rate and the development of acquired resistance among the initial beneficiaries continue to present significant challenges across multiple malignancies, including gastric cancer (GC). Considering that the interactions of tumor stroma, especially the cancer-associated fibroblasts (CAFs), with immune and tumor cells, play indispensable roles in tumor progression, tumor microenvironment remodeling and therapeutic responsiveness, in-depth exploration on the roles of CAFs and pivotal mediators of their functions may provide novel clues to increase the effectiveness of current immunotherapeutic drugs and further achieve synergistic antitumor response. Herein, through the consensus clustering of canonical biomarkers, three GC subclasses with different abundance of CAFs were virtually microdissected in four integrated bulk cohorts encompassing 2148 GC patients from 11 independent datasets. An extensive immunogenomic analysis revealed that tumors with high CAFs infiltration were characterized with unfavorable outcomes, aggressive phenotypes, decreased tumor immunogenicity, high risk of immune evasion and thus immunotherapeutic resistance. By leveraging large-scale single-cell transcriptomic profiling, a series of CAF-secreted proteins were identified, among which the SERPINE2 was confirmed to be restrictively enriched in stromal fibroblasts of GC tissues and contribute to promoting a protumor milieu and fostering an immunosuppressive microenvironment via bioinformatics computations and tissue microarray analysis. Moreover, pan-cancer investigations generalized the immunological roles of SERPINE2, especially in pan-gastrointestinal malignancies, with multiple real-world immunotherapy cohorts further confirming its implications on predicting immunotherapeutic efficacy. In conclusion, these findings suggest that the CAF-derived SERPINE2 is a promising immune-oncology target with therapeutic implications to further synergize the immunotherapeutic combinations.

Crosstalk between gut microbiota and metastasis in colorectal cancer: implication of neutrophil extracellular traps

Primary colorectal cancer (CRC) often leads to liver metastasis, possibly due to the formation of pre-metastatic niche (PMN) in liver. Thus, unravelling the key modulator in metastasis is important for the development of clinical therapies. Gut microbiota dysregulation is a key event during CRC progression and metastasis. Numerous studies have elucidated the correlation between specific gut bacteria strains (e.g., pks + E. coli and Bacteroides fragilis) and CRC initiation, and gut bacteria translocation is commonly witnessed during CRC progression. Gut microbiota shapes tumor microenvironment (TME) through direct contact with immune cells or through its functional metabolites. However, how gut microbiota facilitates CRC metastasis remains controversial. Meanwhile, recent studies identify the dissemination of bacteria from gut lumen to liver, suggesting the role of gut microbiota in shaping tumor PMN. A pro-tumoral PMN is characterized by the infiltration of immunosuppressive cells and increased pro-inflammatory immune responses. Notably, neutrophils form web-like structures known as neutrophil extracellular traps (NETs) both in primary TME and metastatic sites, NETs are involved in cancer progression and metastasis. In this review, we focus on the role of gut microbiota in CRC progression and metastasis, highlight the multiple functions of different immune cell types in TME, especially neutrophils and NETs, discuss the possible mechanisms of gut microbiota in shaping PMN formation, and provide therapeutical indications in clinic.

Cancer-associated fibroblasts undergoing neoadjuvant chemotherapy suppress rectal cancer revealed by single-cell and spatial transcriptomics

Neoadjuvant chemotherapy (NAC) for rectal cancer (RC) shows promising clinical response. The modulation of the tumor microenvironment (TME) by NAC and its association with therapeutic response remain unclear. Here, we use single-cell RNA sequencing and spatial transcriptome sequencing to examine the cell dynamics in 29 patients with RC, who are sampled pairwise before and after treatment. We construct a high-resolution cellular dynamic landscape remodeled by NAC and their associations with therapeutic response. NAC markedly reshapes the populations of cancer-associated fibroblasts (CAFs), which is strongly associated with therapeutic response. The remodeled CAF subsets regulate the TME through spatial recruitment and crosstalk to activate immunity and suppress tumor progression through multiple cytokines, including CXCL12, SLIT2, and DCN. In contrast, the epithelial-mesenchymal transition of malignant cells is upregulated by CAF_FAP through MIR4435-2HG induction, resulting in worse outcomes. Our study demonstrates that NAC inhibits tumor progression and modulates the TME by remodeling CAFs.

Suppression of epithelial proliferation and tumorigenesis by immunoglobulin A

Immunoglobulin A (IgA) is the most abundant antibody isotype produced across mammals and plays a specialized role in mucosal homeostasis 1 . Constantly secreted into the lumen of the intestine, IgA binds commensal microbiota to regulate their colonization and function 2,3 , with unclear implications for health. IgA deficiency is common in humans but is difficult to study due to its complex etiology and comorbidities 4-8 . Using genetically and environmentally controlled mice, here we show that IgA-deficient animals have a baseline alteration in the colon epithelium that increases susceptibility to multiple models of colorectal cancer. Transcriptome, imaging, and flow cytometry-based analyses revealed that, in the absence of IgA, colonic epithelial cells induce antibacterial factors and accelerate cell cycling in response to the microbiota. Oral treatment with IgA was sufficient to suppress aberrant epithelial proliferation independently of bacterial binding, suggesting that IgA provides a feedback signal to epithelial cells in parallel with its known roles in microbiome shaping. In a primary colonic organoid culture system, IgA directly suppresses epithelial growth. Conversely, the susceptibility of IgA-deficient mice to colorectal cancer was reversed by Notch inhibition to suppress the absorptive colonocyte developmental program, or by inhibition of the cytokine MIF, the receptor for which was upregulated in stem cells of IgA-deficient animals. These studies demonstrate a homeostatic function for IgA in tempering physiological epithelial responses to microbiota to maintain mucosal health.

CD8+ T cells in the cancer-immunity cycle

CD8+ T cells are end effectors of cancer immunity. Most forms of effective cancer immunotherapy involve CD8+ T cell effector function. Here, we review the current understanding of T cell function in cancer, focusing on key CD8+ T cell subtypes and states. We discuss factors that influence CD8+ T cell differentiation and function in cancer through a framework that incorporates the classic three-signal model and a fourth signal-metabolism-and also consider the impact of the tumor microenvironment from a T cell perspective. We argue for the notion of immunotherapies as "pro-drugs" that act to augment or modulate T cells, which ultimately serve as the drug in vivo, and for the importance of overall immune health in cancer treatment and prevention. The progress in understanding T cell function in cancer has and will continue to improve harnessing of the immune system across broader tumor types to benefit more patients.

Cancer cell-intrinsic mechanisms driving acquired immune tolerance

Immune evasion is a hallmark of cancer, enabling tumors to survive contact with the host immune system and evade the cycle of immune recognition and destruction. Here, we review the current understanding of the cancer cell-intrinsic factors driving immune evasion. We focus on T cells as key effectors of anti-cancer immunity and argue that cancer cells evade immune destruction by gaining control over pathways that usually serve to maintain physiological tolerance to self. Using this framework, we place recent mechanistic advances in the understanding of cancer immune evasion into broad categories of control over T cell localization, antigen recognition, and acquisition of optimal effector function. We discuss the redundancy in the pathways involved and identify knowledge gaps that must be overcome to better target immune evasion, including the need for better, routinely available tools that incorporate the growing understanding of evasion mechanisms to stratify patients for therapy and trials.

Tertiary lymphoid structures and B cells: An intratumoral immunity cycle

The generation of anti-tumor immunity in the draining lymph nodes is known as the cancer immunity cycle. Accumulating evidence supports the occurrence of such a cycle at tumor sites in the context of chronic inflammation. Here, we review the role of tertiary lymphoid structures (TLS) in the generation of T and B cell immunities, focusing on the impact of B cells that undergo full maturation, resulting in the generation of plasma cells (PCs) producing high-affinity IgG and IgA antibodies. In this context, we propose that antibodies binding to tumor cells induce macrophage or natural killer (NK)-cell-dependent apoptosis. Subsequently, released antigen-antibody complexes are internalized and processed by dendritic cells (DCs), amplifying antigen presentation to T cells. Immune complexes may also be fixed by follicular DCs (FDCs) in TLS, thereby increasing memory B cell responses. This amplification loop creates an intra-tumoral immunity cycle, capable of increasing sensitivity of tumors to immunotherapy even in cancers with low mutational burden.

TCAF2 in Pericytes Promotes Colorectal Cancer Liver Metastasis via Inhibiting Cold-Sensing TRPM8 Channel

Hematogenous metastasis is the main approach for colorectal cancer liver metastasis (CRCLM). However, as the gatekeepers in the tumor vessels, the role of TPCs in hematogenous metastasis remains largely unknown, which may be attributed to the lack of specific biomarkers for TPC isolation. Here, microdissection combined with a pericyte medium-based approach is developed to obtain TPCs from CRC patients. Proteomic analysis reveals that TRP channel-associated factor 2 (TCAF2), a partner protein of the transient receptor potential cation channel subfamily M member 8 (TRPM8), is overexpressed in TPCs from patients with CRCLM. TCAF2 in TPCs is correlated with liver metastasis, short overall survival, and disease-free survival in CRC patients. Gain- and loss-of-function experiments validate that TCAF2 in TPCs promotes tumor cell motility, epithelial-mesenchymal transition (EMT), and CRCLM, which is attenuated in pericyte-conditional Tcaf2-knockout mice. Mechanistically, TCAF2 inhibits the expression and activity of TRPM8, leading to Wnt5a secretion in TPCs, which facilitates EMT via the activation of the STAT3 signaling pathway in tumor cells. Menthol, a TRPM8 agonist, significantly suppresses Wnt5a secretion in TPCs and CRCLM. This study reveals the previously unidentified pro-metastatic effects of TPCs from the perspective of cold-sensory receptors, providing a promising diagnostic biomarker and therapeutic target for CRCLM.

The diverse nature of intestinal fibroblasts in development, homeostasis, and disease

Only in recent years have we begun to appreciate the involvement of fibroblasts in intestinal development, tissue homeostasis, and disease. These insights followed the advent of single-cell transcriptomics that allowed researchers to explore the heterogeneity of intestinal fibroblasts in unprecedented detail. Since researchers often defined cell types and their associated function based on the biological process they studied, there are a plethora of partially overlapping markers for different intestinal fibroblast populations. This ambiguity complicates putting different research findings into context. Here, we provide a census on the function and identity of intestinal fibroblasts in mouse and human. We propose a simplified framework consisting of three colonic and four small intestinal fibroblast populations to aid navigating the diversity of intestinal fibroblasts.

ATR inhibition induces synthetic lethality in mismatch repair-deficient cells and augments immunotherapy

The mismatch repair (MMR) deficiency of cancer cells drives mutagenesis and offers a useful biomarker for immunotherapy. However, many MMR-deficient (MMR-d) tumors do not respond to immunotherapy, highlighting the need for alternative approaches to target MMR-d cancer cells. Here, we show that inhibition of the ATR kinase preferentially kills MMR-d cancer cells. Mechanistically, ATR inhibitor (ATRi) imposes synthetic lethality on MMR-d cells by inducing DNA damage in a replication- and MUS81 nuclease-dependent manner. The DNA damage induced by ATRi is colocalized with both MSH2 and PCNA, suggesting that it arises from DNA structures recognized by MMR proteins during replication. In syngeneic mouse models, ATRi effectively reduces the growth of MMR-d tumors. Interestingly, the antitumor effects of ATRi are partially due to CD8+ T cells. In MMR-d cells, ATRi stimulates the accumulation of nascent DNA fragments in the cytoplasm, activating the cGAS-mediated interferon response. The combination of ATRi and anti-PD-1 antibody reduces the growth of MMR-d tumors more efficiently than ATRi or anti-PD-1 alone, showing the ability of ATRi to augment the immunotherapy of MMR-d tumors. Thus, ATRi selectively targets MMR-d tumor cells by inducing synthetic lethality and enhancing antitumor immunity, providing a promising strategy to complement and augment MMR deficiency-guided immunotherapy.

Spatially Resolved Multi-Omics Single-Cell Analyses Inform Mechanisms of Immune Dysfunction in Pancreatic Cancer

BACKGROUND & AIMS

As pancreatic ductal adenocarcinoma (PDAC) continues to be recalcitrant to therapeutic interventions, including poor response to immunotherapy, albeit effective in other solid malignancies, a more nuanced understanding of the immune microenvironment in PDAC is urgently needed. We aimed to unveil a detailed view of the immune micromilieu in PDAC using a spatially resolved multimodal single-cell approach.

METHODS

We applied single-cell RNA sequencing, spatial transcriptomics, multiplex immunohistochemistry, and mass cytometry to profile the immune compartment in treatment-naïve PDAC tumors and matched adjacent normal pancreatic tissue, as well as in the systemic circulation. We determined prognostic associations of immune signatures and performed a meta-analysis of the immune microenvironment in PDAC and lung adenocarcinoma on single-cell level.

RESULTS

We provided a spatially resolved fine map of the immune landscape in PDAC. We substantiated the exhausted phenotype of CD8 T cells and immunosuppressive features of myeloid cells, and highlighted immune subsets with potentially underappreciated roles in PDAC that diverged from immune populations within adjacent normal areas, particularly CD4 T cell subsets and natural killer T cells that are terminally exhausted and acquire a regulatory phenotype. Differential analysis of immune phenotypes in PDAC and lung adenocarcinoma revealed the presence of extraordinarily immunosuppressive subtypes in PDAC, along with a distinctive immune checkpoint composition.

CONCLUSIONS

Our study sheds light on the multilayered immune dysfunction in PDAC and presents a holistic view of the immune landscape in PDAC and lung adenocarcinoma, providing a comprehensive resource for functional studies and the exploration of therapeutically actionable targets in PDAC.

Spatiotemporal and cell-state control of antigen presentation during tolerance and immunity

Effective adaptive immunity is rendered possible by highly organized tissue architecture and coordinated cellular crosstalk. While detailed spatiotemporal analyses of antigen presentation and adaptive immune activation in secondary lymphoid tissues have been a major focus of study, it is clear that antigen presentation in other tissues also plays a critical role in shaping the immune response. In this article, we concentrate on two opposing aspects of adaptive immunity: tolerance and antitumor immunity, to illustrate how a complex set of antigen presentation mechanisms contributes to maintaining a delicate balance between robust immunity and avoidance of autoimmune pathology. We emphasize the importance of how immune cell identity, state, and location collectively determine the nature of adaptive immune responses.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

Intertumoral lineage diversity and immunosuppressive transcriptional programs in well-differentiated gastroenteropancreatic neuroendocrine tumors

Neuroendocrine tumors (NETs) are rare cancers that most often arise in the gastrointestinal tract and pancreas. The fundamental mechanisms driving gastroenteropancreatic (GEP)-NET growth remain incompletely elucidated; however, the heterogeneous clinical behavior of GEP-NETs suggests that both cellular lineage dynamics and tumor microenvironment influence tumor pathophysiology. Here, we investigated the single-cell transcriptomes of tumor and immune cells from patients with gastroenteropancreatic NETs. Malignant GEP-NET cells expressed genes and regulons associated with normal, gastrointestinal endocrine cell differentiation, and fate determination stages. Tumor and lymphoid compartments sparsely expressed immunosuppressive targets commonly investigated in clinical trials, such as the programmed cell death protein-1/programmed death ligand-1 axis. However, infiltrating myeloid cell types within both primary and metastatic GEP-NETs were enriched for genes encoding other immune checkpoints, including VSIR (VISTA), HAVCR2 (TIM3), LGALS9 (Gal-9), and SIGLEC10. Our findings highlight the transcriptomic heterogeneity that distinguishes the cellular landscapes of GEP-NET anatomic subtypes and reveal potential avenues for future precision medicine therapeutics.

An immune cell atlas reveals the dynamics of human macrophage specification during prenatal development

Macrophages are heterogeneous and play critical roles in development and disease, but their diversity, function, and specification remain inadequately understood during human development. We generated a single-cell RNA sequencing map of the dynamics of human macrophage specification from PCW 4-26 across 19 tissues. We identified a microglia-like population and a proangiogenic population in 15 macrophage subtypes. Microglia-like cells, molecularly and morphologically similar to microglia in the CNS, are present in the fetal epidermis, testicle, and heart. They are the major immune population in the early epidermis, exhibit a polarized distribution along the dorsal-lateral-ventral axis, and interact with neural crest cells, modulating their differentiation along the melanocyte lineage. Through spatial and differentiation trajectory analysis, we also showed that proangiogenic macrophages are perivascular across fetal organs and likely yolk-sac-derived as microglia. Our study provides a comprehensive map of the heterogeneity and developmental dynamics of human macrophages and unravels their diverse functions during development.

Cracking the pattern of tumor evolution based on single-cell copy number alterations

Copy number alterations (CNAs) are a key characteristic of tumor development and progression. The accumulation of various CNAs during tumor development plays a critical role in driving tumor evolution. Heterogeneous clones driven by distinct CNAs have different selective advantages, leading to differential patterns of tumor evolution that are essential for developing effective cancer therapies. Recent advances in single-cell sequencing technology have enabled genome-wide copy number profiling of tumor cell populations at single-cell resolution. This has made it possible to explore the evolutionary patterns of CNAs and accurately discover the mechanisms of intra-tumor heterogeneity. Here, we propose a two-step statistical approach that distinguishes neutral, linear, branching and punctuated evolutionary patterns for a tumor cell population based on single-cell copy number profiles. We assessed our approach using a variety of simulated and real single-cell genomic and transcriptomic datasets, demonstrating its high accuracy and robustness in predicting tumor evolutionary patterns. We applied our approach to single-cell DNA sequencing data from 20 breast cancer patients and observed that punctuated evolution is the dominant evolutionary pattern in breast cancer. Similar conclusions were drawn when applying the approach to single-cell RNA sequencing data obtained from 132 various cancer patients. Moreover, we found that differential immune cell infiltration is associated with specific evolutionary patterns. The source code of our study is available at https://github.com/FangWang-SYSU/PTEM.

Transcriptomic Signatures of MSI-High Metastatic Colorectal Cancer Predict Efficacy of Immune Checkpoint Inhibitors

PURPOSE

Microsatellite instability (MSI) is currently the only predictive biomarker of efficacy of immune checkpoint inhibitors (ICI) in metastatic colorectal cancers (mCRC). However, 10% to 40% of patients with MSI mCRC will experience a primary resistance to ICI.

EXPERIMENTAL DESIGN

In two cohorts of patients with MSI mCRC treated with ICI (exploratory, N = 103; validation, N = 35), 3' RNA sequencing was performed from primary tumors. Previously described single-cell transcriptomic signatures of tumor microenvironment (TME) were analyzed.

RESULTS

In the exploratory cohort, the unsupervised clustering allowed the identification of three clusters of tumors with distinct transcriptional profiles: cluster A ("stromalHIGH-proliferationLOW"), cluster B ("stromalHIGH-proliferationMED"), and cluster C ("stromalLOW-proliferationHIGH"), with an enrichment of patients progressing at first disease assessment under ICI in cluster A (30% vs. 12% in cluster B and 8.1% in cluster C; P = 0.074). Progression-free survival (PFS) was also significantly shorter in patients belonging to cluster A, compared with clusters B or C (P < 0.001) with 2-year PFS rates of 33.5%, 80.5%, and 78.3%, respectively. In multivariate analysis, PFS was still significantly longer in patients belonging to cluster B [HR, 0.19; 95% confidence interval (CI), 0.08-0.45; P < 0.001] and cluster C (HR, 0.25; 95% CI, 0.10-0.59; P = 0.02), compared with patients belonging to cluster A. The association of this clustering with PFS under ICI was confirmed in the validation cohort. PFS related to non-ICI-based regimens was not significantly different according to cluster.

CONCLUSIONS

This unsupervised transcriptomic classification identified three groups of MSI mCRCs with different compositions of TME cells and proliferative capacities of TME/tumor cells. The "stromalHIGH-proliferationLOW" cluster is associated with a poorer prognosis with ICI treatment.

Global proteomic identifies multiple cancer-related signaling pathways altered by a gut pathobiont associated with colorectal cancer

In this work, we investigated the oncogenic role of Streptococcus gallolyticus subsp. gallolyticus (SGG), a gut bacterium associated with colorectal cancer (CRC). We showed that SGG UCN34 accelerates colon tumor development in a chemically induced CRC murine model. Full proteome and phosphoproteome analysis of murine colons chronically colonized by SGG UCN34 revealed that 164 proteins and 725 phosphorylation sites were differentially regulated. Ingenuity Pathway Analysis (IPA) indicates a pro-tumoral shift specifically induced by SGG UCN34, as ~ 90% of proteins and phosphoproteins identified were associated with digestive cancer. Comprehensive analysis of the altered phosphoproteins using ROMA software revealed up-regulation of several cancer hallmark pathways such as MAPK, mTOR and integrin/ILK/actin, affecting epithelial and stromal colonic cells. Importantly, an independent analysis of protein arrays of human colon tumors colonized with SGG showed up-regulation of PI3K/Akt/mTOR and MAPK pathways, providing clinical relevance to our findings. To test SGG's capacity to induce pre-cancerous transformation of the murine colonic epithelium, we grew ex vivo organoids which revealed unusual structures with compact morphology. Taken together, our results demonstrate the oncogenic role of SGG UCN34 in a murine model of CRC associated with activation of multiple cancer-related signaling pathways.

IL-2Rα-biased agonist enhances antitumor immunity by invigorating tumor-infiltrating CD25+CD8+ T cells

To avoid regulatory T cell promotion and vascular toxicity, the interleukin-2 receptor-β/interleukin-2 receptor-γ (IL-2Rβγ)-biased approach is used by most IL-2 analogs in immuno-oncology. However, recent clinical disappointments in these IL-2 agonists have questioned this strategy. Here we show that both wild-type (IL-2wt) and IL-2Rβγ-attenuated (IL-2α-bias) agonists that preserve IL-2Rα (CD25) activities can effectively expand tumor-specific CD8+ T cells (TSTs) and exhibit better antitumor efficacy and safety than the 'non-α' counterpart (IL-2nα). Mechanistically, TSTs coexpress elevated CD25 and PD-1 and are more susceptible to stimulation by IL-2Rα-proficient agonists. Moreover, the antitumor efficacy of anti-PD-1 depends on activation of PD-1+CD25+ TSTs through autocrine IL-2-CD25 signaling. In individuals with cancer, a low IL-2 signature correlates with non-responsiveness to anti-PD-1 treatment. In mouse models, IL-2α-bias, but not IL-2nα, restores the IL-2 signature and synergizes with anti-PD-1 to eradicate large established tumors. These findings underscore the indispensable function of CD25 in IL-2-based immunotherapy and provide rationales for evaluating IL-2Rα-biased agonists in individuals with cancer.

Deep learning exploration of single-cell and spatially resolved cancer transcriptomics to unravel tumour heterogeneity

Tumour heterogeneity is one of the critical confounding aspects in decoding tumour growth. Malignant cells display variations in their gene transcription profiles and mutation spectra even when originating from a single progenitor cell. Single-cell and spatial transcriptomics sequencing have recently emerged as key technologies for unravelling tumour heterogeneity. Single-cell sequencing promotes individual cell-type identification through transcriptome-wide gene expression measurements of each cell. Spatial transcriptomics facilitates identification of cell-cell interactions and the structural organization of heterogeneous cells within a tumour tissue through associating spatial RNA abundance of cells at distinct spots in the tissue section. However, extracting features and analyzing single-cell and spatial transcriptomics data poses challenges. Single-cell transcriptome data is extremely noisy and its sparse nature and dropouts can lead to misinterpretation of gene expression and the misclassification of cell types. Deep learning predictive power can overcome data challenges, provide high-resolution analysis and enhance precision oncology applications that involve early cancer prognosis, diagnosis, patient survival estimation and anti-cancer therapy planning. In this paper, we provide a background to and review of the recent progress of deep learning frameworks to investigate tumour heterogeneity using both single-cell and spatial transcriptomics data types.

A Roadmap for the Human Gut Cell Atlas

The number of studies investigating the human gastrointestinal tract using various single-cell profiling methods has increased substantially in the past few years. Although this increase provides a unique opportunity for the generation of the first comprehensive Human Gut Cell Atlas (HGCA), there remains a range of major challenges ahead. Above all, the ultimate success will largely depend on a structured and coordinated approach that aligns global efforts undertaken by a large number of research groups. In this Roadmap, we discuss a comprehensive forward-thinking direction for the generation of the HGCA on behalf of the Gut Biological Network of the Human Cell Atlas. Based on the consensus opinion of experts from across the globe, we outline the main requirements for the first complete HGCA by summarizing existing data sets and highlighting anatomical regions and/or tissues with limited coverage. We provide recommendations for future studies and discuss key methodologies and the importance of integrating the healthy gut atlas with related diseases and gut organoids. Importantly, we critically overview the computational tools available and provide recommendations to overcome key challenges.

T follicular helper cells in cancer, tertiary lymphoid structures, and beyond

With the emergence and success of checkpoint blockade immunotherapy, immuno-oncology has primarily focused on CD8 T cells, whose cytotoxic programs directly target tumor cells. However, the limited response rate of current immunotherapy regimens has prompted investigation into other types of tumor-infiltrating immune cells, such as CD4 T cells and B cells, and how they interact with CD8 T cells in a coordinated network. Recent studies have demonstrated the potential therapeutic benefits of CD4 T follicular helper (TFH) cells and B cells in cancer, highlighting the important role of their crosstalk and interactions with other immune cell components in the tumor microenvironment. These interactions also occur in tumor-associated tertiary lymphoid structures (TLS), which resemble secondary lymphoid organs (SLOs) with orchestrated vascular, chemokine, and cellular infrastructures that support the developmental pathways of functional immune cells. In this review, we discuss recent breakthroughs on TFH biology and T cell-B cell interactions in tumor immunology, and their potential as novel therapeutic targets to advance cancer treatment.

Bacteria in cancer initiation, promotion and progression

Cancer cells originate from a series of acquired genetic mutations that can drive their uncontrolled cell proliferation and immune evasion. Environmental factors, including the microorganisms that colonize the human body, can shift the metabolism, growth pattern and function of neoplastic cells and shape the tumour microenvironment. Dysbiosis of the gut microbiome is now recognized as a hallmark of cancer by the scientific community. However, only a few microorganisms have been identified that directly initiate tumorigenesis or skew the immune system to generate a tumour-permissive milieu. Over the past two decades, research on the human microbiome and its functionalities within and across individuals has revealed microbiota-focused strategies for health and disease. Here, we review the evolving understanding of the mechanisms by which the microbiota acts in cancer initiation, promotion and progression. We explore the roles of bacteria in gastrointestinal tract malignancies and cancers of the lung, breast and prostate. Finally, we discuss the promises and limitations of targeting or harnessing bacteria in personalized cancer prevention, diagnostics and treatment.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Intercellular cross-talk through lineage-specific gap junction of cancer-associated fibroblasts related to stromal fibrosis and prognosis

Stromal fibrosis in cancer is usually associated with poor prognosis and chemotherapy resistance. It is thought to be caused by fibroblasts; however, the exact mechanism is not yet well understood. The study aimed to identify lineage-specific cancer-associated fibroblast (CAF) subgroup and their associations with extracellular matrix remodeling and clinical significances in various tumor types using single-cell and bulk RNA sequencing data. Through unsupervised clustering, six subclusters of CAFs were identified, including a cluster with exclusively high gap junction protein beta-2 (GJB2) expression. This cluster was named GJB2-positive CAF. It was found to be a unique subgroup of terminally differentiated CAFs associated with collagen gene expression and extracellular matrix remodeling. GJB2-positive CAFs showed higher communication frequency with vascular endothelial cells and cancer cells than GJB2-negative CAFs. Moreover, GJB2 was poorly expressed in normal tissues, indicating that its expression is dependent on interaction with other cells, including vascular endothelial cells and cancer cells. Finally, the study investigated the clinical significance of GJB2 signature score for GJB2-positive CAFs in cancer and found a correlation with poor prognosis. These results suggest that GJB2-positive CAF is a unique fibroblast subtype involved in extracellular matrix remodeling, with significant clinical implications in cancer.

scGEM: Unveiling the Nested Tree-Structured Gene Co-Expressing Modules in Single Cell Transcriptome Data

BACKGROUND

Single-cell transcriptome analysis has fundamentally changed biological research by allowing higher-resolution computational analysis of individual cells and subsets of cell types. However, few methods have met the need to recognize and quantify the underlying cellular programs that determine the specialization and differentiation of the cell types.

METHODS

In this study, we present scGEM, a nested tree-structured nonparametric Bayesian model, to reveal the gene co-expression modules (GEMs) reflecting transcriptome processes in single cells.

RESULTS

We show that scGEM can discover shared and specialized transcriptome signals across different cell types using peripheral blood mononuclear single cells and early brain development single cells. scGEM outperformed other methods in perplexity and topic coherence (p < 0.001) on our simulation data. Larger datasets, deeper trees and pre-trained models are shown to be positively associated with better scGEM performance. The GEMs obtained from triple-negative breast cancer single cells exhibited better correlations with lymphocyte infiltration (p = 0.009) and the cell cycle (p < 0.001) than other methods in additional validation on the bulk RNAseq dataset.

CONCLUSIONS

Altogether, we demonstrate that scGEM can be used to model the hidden cellular functions of single cells, thereby unveiling the specialization and generalization of transcriptomic programs across different types of cells.

Progressive plasticity during colorectal cancer metastasis

Metastasis is the principal cause of cancer death, yet we lack an understanding of metastatic cell states, their relationship to primary tumor states, and the mechanisms by which they transition. In a cohort of biospecimen trios from same-patient normal colon, primary and metastatic colorectal cancer, we show that while primary tumors largely adopt LGR5 + intestinal stem-like states, metastases display progressive plasticity. Loss of intestinal cell states is accompanied by reprogramming into a highly conserved fetal progenitor state, followed by non-canonical differentiation into divergent squamous and neuroendocrine-like states, which is exacerbated by chemotherapy and associated with poor patient survival. Using matched patient-derived organoids, we demonstrate that metastatic cancer cells exhibit greater cell-autonomous multilineage differentiation potential in response to microenvironment cues than their intestinal lineage-restricted primary tumor counterparts. We identify PROX1 as a stabilizer of intestinal lineage in the fetal progenitor state, whose downregulation licenses non-canonical reprogramming.

Evolution of immune and stromal cell states and ecotypes during gastric adenocarcinoma progression

Understanding tumor microenvironment (TME) reprogramming in gastric adenocarcinoma (GAC) progression may uncover novel therapeutic targets. Here, we performed single-cell profiling of precancerous lesions, localized and metastatic GACs, identifying alterations in TME cell states and compositions as GAC progresses. Abundant IgA+ plasma cells exist in the premalignant microenvironment, whereas immunosuppressive myeloid and stromal subsets dominate late-stage GACs. We identified six TME ecotypes (EC1-6). EC1 is exclusive to blood, while EC4, EC5, and EC2 are highly enriched in uninvolved tissues, premalignant lesions, and metastases, respectively. EC3 and EC6, two distinct ecotypes in primary GACs, associate with histopathological and genomic characteristics, and survival outcomes. Extensive stromal remodeling occurs in GAC progression. High SDC2 expression in cancer-associated fibroblasts (CAFs) is linked to aggressive phenotypes and poor survival, and SDC2 overexpression in CAFs contributes to tumor growth. Our study provides a high-resolution GAC TME atlas and underscores potential targets for further investigation.

Single-Cell RNA Sequencing for Studying Human Cancers

Since the first publication a decade ago describing the use of single-cell RNA sequencing (scRNA-seq) in the context of cancer, over 200 datasets and thousands of scRNA-seq studies have been published in cancer biology. scRNA-seq technologies have been applied across dozens of cancer types and a diverse array of study designs to improve our understanding of tumor biology, the tumor microenvironment, and therapeutic responses, and scRNA-seq is on the verge of being used to improve decision-making in the clinic. Computational methodologies and analytical pipelines are key in facilitating scRNA-seq research. Numerous computational methods utilizing the most advanced tools in data science have been developed to extract meaningful insights. Here, we review the advancements in cancer biology gained by scRNA-seq and discuss the computational challenges of the technology that are specific to cancer research.

CXCL9:SPP1 macrophage polarity identifies a network of cellular programs that control human cancers

Tumor microenvironments (TMEs) influence cancer progression but are complex and often differ between patients. Considering that microenvironment variations may reveal rules governing intratumoral cellular programs and disease outcome, we focused on tumor-to-tumor variation to examine 52 head and neck squamous cell carcinomas. We found that macrophage polarity-defined by CXCL9 and SPP1 (CS) expression but not by conventional M1 and M2 markers-had a noticeably strong prognostic association. CS macrophage polarity also identified a highly coordinated network of either pro- or antitumor variables, which involved each tumor-associated cell type and was spatially organized. We extended these findings to other cancer indications. Overall, these results suggest that, despite their complexity, TMEs coordinate coherent responses that control human cancers and for which CS macrophage polarity is a relevant yet simple variable.

Applications of multi-omics analysis in human diseases

Multi-omics usually refers to the crossover application of multiple high-throughput screening technologies represented by genomics, transcriptomics, single-cell transcriptomics, proteomics and metabolomics, spatial transcriptomics, and so on, which play a great role in promoting the study of human diseases. Most of the current reviews focus on describing the development of multi-omics technologies, data integration, and application to a particular disease; however, few of them provide a comprehensive and systematic introduction of multi-omics. This review outlines the existing technical categories of multi-omics, cautions for experimental design, focuses on the integrated analysis methods of multi-omics, especially the approach of machine learning and deep learning in multi-omics data integration and the corresponding tools, and the application of multi-omics in medical researches (e.g., cancer, neurodegenerative diseases, aging, and drug target discovery) as well as the corresponding open-source analysis tools and databases, and finally, discusses the challenges and future directions of multi-omics integration and application in precision medicine. With the development of high-throughput technologies and data integration algorithms, as important directions of multi-omics for future disease research, single-cell multi-omics and spatial multi-omics also provided a detailed introduction. This review will provide important guidance for researchers, especially who are just entering into multi-omics medical research.

Cuproptosis-related miRNAs signature and immune infiltration characteristics in colorectal cancer

BACKGROUND

A novel form of cell death termed cuproptosis was proposed recently. miRNAs play important roles in colorectal cancer (CRC). However, their relationships have not been reported.

METHODS

miRNAs that negatively regulate 16 cuproptosis regulators were predicted using Targetscan database. The univariate Cox, LASSO, and multivariate Cox regression analyses were performed to select cuproptosis-related miRNAs. GSEA and ssGSEA analysis was carried out for functional enrichment analysis. The immune cell proportion score (IPS) and the efficiencies of multiple chemotherapy drugs were compared between different risk groups. The CCK8, cell colony, edu, and flow cytometry assays were performed to validate the roles of miRNA. Luciferase reporter assay confirmed the regulatory mechanism of miRNA on cuproptosis.

RESULTS

Six cuproptosis-related miRNAs (hsa-miR-653, hsa-miR-216a, hsa-miR-3684, hsa-miR-4437, hsa-miR-641, and hsa-miR-552) were screened out for model construction. The risk score could act as an independent prognostic indicator in CRC (p < 0.001, 95% HR = 1.243 (1.129-1.369)). The nomogram could efficiently predict the overall survival rate (AUC = 0.836). Then, the level of immunosuppressive pathways, immunosuppressive cells, stromal-activated genes, and stromal score was higher in the high-risk group. The IPS analysis showed a better response to immunotherapy in the low-risk group. Also, the risk score was closely correlated with efficiencies of multiple chemotherapy drugs. Furthermore, miR-653 was highly expressed in CRC tissues (p < 0.001), closely correlated with T stage (p < 0.001), metastasis (p < 0.001), and tumor stage (p < 0.001). High expression of miR-653 predicted a shorter overall survival (p = 0.0282) and disease-free survival (p = 0.0056). In addition, miR-653 promoted cell proliferation, inhibited apoptosis, and negatively regulated the expression of DLD through directly binding to the 3'-UTR of DLD mRNA.

CONCLUSION

We constructed a cuproptosis-related miRNA signature for the prediction of CRC patient survival and immunotherapy sensitivity. miR-653 was highly expressed in CRC tissues, promoted cell proliferation, and inhibited apoptosis by negatively regulating the expression of DLD.

γδ T cell dichotomy with opposing cytotoxic and wound healing functions in human solid tumors

γδ T cells are important tissue-resident, innate T cells that are critical for tissue homeostasis. γδ cells are associated with positive prognosis in most tumors; however, little is known about their heterogeneity in human cancers. Here, we phenotyped innate and adaptive cells in human colorectal (CRC) and endometrial cancer. We found striking differences in γδ subsets and function in tumors compared to normal tissue, and in the γδ subsets present in tumor types. In CRC, an amphiregulin (AREG)-producing subset emerges, while endometrial cancer is infiltrated by cytotoxic cells. In humanized CRC models, tumors induced this AREG phenotype in Vδ1 cells after adoptive transfer. To exploit the beneficial roles of γδ cells for cell therapy, we developed an expansion method that enhanced cytotoxic function and boosted metabolic flexibility, while eliminating AREG production, achieving greater tumor infiltration and tumor clearance. This method has broad applications in cellular therapy as an 'off-the-shelf' treatment option.

Metformin Reprograms Tryptophan Metabolism to Stimulate CD8+ T-cell Function in Colorectal Cancer

Colorectal carcinogenesis coincides with immune cell dysfunction. Metformin has been reported to play a role in stimulating antitumor immunity, suggesting it could be used to overcome immunosuppression in colorectal cancer. Herein, using single-cell RNA sequencing (scRNA-seq), we showed that metformin remodels the immune landscape of colorectal cancer. In particular, metformin treatment expanded the proportion of CD8+ T cells and potentiated their function. Analysis of the metabolic activities of cells in the colorectal cancer tumor microenvironment (TME) at a single-cell resolution demonstrated that metformin reprogrammed tryptophan metabolism, which was reduced in colorectal cancer cells and increased in CD8+ T cells. Untreated colorectal cancer cells outcompeted CD8+ T cells for tryptophan, leading to impaired CD8+ T-cell function. Metformin in turn reduced tryptophan uptake by colorectal cancer cells, thereby restoring tryptophan availability for CD8+ T cells and increasing their cytotoxicity. Metformin inhibited tryptophan uptake in colorectal cancer cells by downregulating MYC, which led to a reduction in the tryptophan transporter SLC7A5. This work highlights metformin as an essential regulator of T-cell antitumor immunity by reprogramming tryptophan metabolism, suggesting it could be a potential immunotherapeutic strategy for treating colorectal cancer.

SIGNIFICANCE

Analysis of the impact of metformin on the colorectal cancer immunometabolic landscape at a single-cell resolution shows that metformin alters cancer cell tryptophan metabolism to stimulate CD8+ T-cell antitumor activity.

Immune escape and resistance to immunotherapy in mismatch repair deficient tumors

Up to 30% of colorectal, endometrial and gastric cancers have a deficiency in mismatch repair (MMR) protein expression due to either germline or epigenetic inactivation. Patients with Lynch Syndrome who inherit an inactive MMR allele have an up to 80% risk for developing a mismatch repair deficient (MMRd) cancer. Due to an inability to repair DNA, MMRd tumors present with genomic instability in microsatellite regions (MS). Tumors with high MS instability (MSI-H) are characterized by an increased frequency of insertion/deletions (indels) that can encode novel neoantigens if they occur in coding regions. The high tumor antigen burden for MMRd cancers is accompanied by an inflamed tumor microenvironment (TME) that contributes to the clinical effectiveness of anti-PD-1 therapy in this patient population. However, between 40 and 70% of MMRd cancer patients do not respond to treatment with PD-1 blockade, suggesting that tumor-intrinsic and -extrinsic resistance mechanisms may affect the success of checkpoint blockade. Immune evasion mechanisms that occur during early tumorigenesis and persist through cancer development may provide a window into resistance pathways that limit the effectiveness of anti-PD-1 therapy. Here, we review the mechanisms of immune escape in MMRd tumors during development and checkpoint blockade treatment, including T cell dysregulation and myeloid cell-mediated immunosuppression in the TME. Finally, we discuss the development of new therapeutic approaches to tackle resistance in MMRd tumors, including cancer vaccines, therapies targeting immunosuppressive myeloid programs, and immune checkpoint combination strategies.

The tumor-stroma ratio and the immune microenvironment improve the prognostic prediction of pancreatic ductal adenocarcinoma

Tumor-infiltrating immune cells and fibroblasts are significant components of the tumor microenvironment (TME) of pancreatic ductal adenocarcinoma (PDAC), and they participate in tumor progression as closely as tumor cells. However, the relationship between the features of the TME and patient outcomes and the interactions among TME components are still unclear. In this study, we evaluated the PDAC TME in terms of the quantity and location of cluster of differentiation (CD)4⁺ T cells, CD8⁺ T cells, macrophages, stromal maturity, and tumor-stroma ratio (TSR), as evaluated by immunohistochemical staining of serial whole-tissue sections from 116 patients with PDAC. The density of T cells and macrophages (mainly activated macrophages) was significantly higher at the invasive margins (IMs) than at the tumor center (TC). CD4⁺ T cells were significantly association with all the other tumor-associated immune cells (TAIs) including CD8, CD68 and CD206 positive cells. Tumors of the non-mature (intermediate and immature) stroma type harbored significantly more CD8⁺ T cells at the IMs and more CD68⁺ macrophages at the IMs and the TC. The density of CD4⁺, CD8⁺, and CD206⁺ cells at the TC; CD206⁺ cells at the IMs; and tumor-node-metastasis (TNM) staging were independent risk factors for patient outcomes, and the c-index of the risk nomogram for predicting the survival probability based on the TME features and TNM staging was 0.772 (95% confidence interval: 0.713-0.832). PDAC harbored a significantly immunosuppressive TME, of which the IMs were the hot zones for TAIs, while cells at the TC were more predictive of prognosis. Our results indicated that the model based on the features of the TME and TNM staging could predict patient outcomes.

SELENOP modifies sporadic colorectal carcinogenesis and WNT signaling activity through LRP5/6 interactions

Although selenium deficiency correlates with colorectal cancer (CRC) risk, the roles of the selenium-rich antioxidant selenoprotein P (SELENOP) in CRC remain unclear. In this study, we defined SELENOP's contributions to sporadic CRC. In human single-cell cRNA-Seq (scRNA-Seq) data sets, we discovered that SELENOP expression rose as normal colon stem cells transformed into adenomas that progressed into carcinomas. We next examined the effects of Selenop KO in a mouse adenoma model that involved conditional, intestinal epithelium-specific deletion of the tumor suppressor adenomatous polyposis coli (Apc) and found that Selenop KO decreased colon tumor incidence and size. We mechanistically interrogated SELENOP-driven phenotypes in tumor organoids as well as in CRC and noncancer cell lines. Selenop-KO tumor organoids demonstrated defects in organoid formation and decreases in WNT target gene expression, which could be reversed by SELENOP restoration. Moreover, SELENOP increased canonical WNT signaling activity in noncancer and CRC cell lines. In defining the mechanism of action of SELENOP, we mapped protein-protein interactions between SELENOP and the WNT coreceptors low-density lipoprotein receptor-related proteins 5 and 6 (LRP5/6). Last, we confirmed that SELENOP-LRP5/6 interactions contributed to the effects of SELENOP on WNT activity. Overall, our results position SELENOP as a modulator of the WNT signaling pathway in sporadic CRC.

Spatial transcriptome unveils a discontinuous inflammatory pattern in proficient mismatch repair colorectal adenocarcinoma

The preexistence of immune cells in the tumor microenvironment substantiates the efficacy of immunotherapy in cancer patients. Although the complex intratumoral immune heterogeneity has been extensively studied in single cell resolution, hi-res spatial investigations are limited. In this study, we performed a spatial transcriptome analysis of 4 colorectal adenocarcinoma specimens and 2 paired distant normal specimens to identify the molecular pattern involved in a discontinuous inflammatory response in pathologically annotated cancer regions. Based on the location of spatially varied gene expression, we unmasked the spatially-varied immune ecosystem and identified the locoregional "warmed-up" immune response in predefined "cold" tumor with substantial infiltration of immune components. This "warmed-up" immune profile was found to be associated with the in-situ copy number variance and the tissue remodeling process. Further, "warmed-up" signature genes indicated improved overall survival in CRC patients obtained from TCGA database.

Drivers of heterogeneity in synovial fibroblasts in rheumatoid arthritis

Inflammation of non-barrier immunologically quiescent tissues is associated with a massive influx of blood-borne innate and adaptive immune cells. Cues from the latter are likely to alter and expand activated states of the resident cells. However, local communications between immigrant and resident cell types in human inflammatory disease remain poorly understood. Here, we explored drivers of fibroblast-like synoviocyte (FLS) heterogeneity in inflamed joints of patients with rheumatoid arthritis using paired single-cell RNA and ATAC sequencing, multiplexed imaging and spatial transcriptomics along with in vitro modeling of cell-extrinsic factor signaling. These analyses suggest that local exposures to myeloid and T cell-derived cytokines, TNF, IFN-γ, IL-1β or lack thereof, drive four distinct FLS states some of which closely resemble fibroblast states in other disease-affected tissues including skin and colon. Our results highlight a role for concurrent, spatially distributed cytokine signaling within the inflamed synovium.

B-cell-specific checkpoint molecules that regulate anti-tumour immunity

The role of B cells in anti-tumour immunity is still debated and, accordingly, immunotherapies have focused on targeting T and natural killer cells to inhibit tumour growth1,2. Here, using high-throughput flow cytometry as well as bulk and single-cell RNA-sequencing and B-cell-receptor-sequencing analysis of B cells temporally during B16F10 melanoma growth, we identified a subset of B cells that expands specifically in the draining lymph node over time in tumour-bearing mice. The expanding B cell subset expresses the cell surface molecule T cell immunoglobulin and mucin domain 1 (TIM-1, encoded by Havcr1) and a unique transcriptional signature, including multiple co-inhibitory molecules such as PD-1, TIM-3, TIGIT and LAG-3. Although conditional deletion of these co-inhibitory molecules on B cells had little or no effect on tumour burden, selective deletion of Havcr1 in B cells both substantially inhibited tumour growth and enhanced effector T cell responses. Loss of TIM-1 enhanced the type 1 interferon response in B cells, which augmented B cell activation and increased antigen presentation and co-stimulation, resulting in increased expansion of tumour-specific effector T cells. Our results demonstrate that manipulation of TIM-1-expressing B cells enables engagement of the second arm of adaptive immunity to promote anti-tumour immunity and inhibit tumour growth.