Extended correlation functions for spatial analysis of multiplex imaging data.

Imaging platforms for generating highly multiplexed histological images are being continually developed and improved. Significant improvements have also been made in the accuracy of methods for automated cell segmentation and classification. However, less attention has focused on the quantification and analysis of the resulting point clouds, which describe the spatial coordinates of individual cells. We focus here on a particular spatial statistical method, the cross-pair correlation function (cross-PCF), which can identify positive and negative spatial correlation between cells across a range of length scales. However, limitations of the cross-PCF hinder its widespread application to multiplexed histology. For example, it can only consider relations between pairs of cells, and cells must be classified using discrete categorical labels (rather than labeling continuous labels such as stain intensity). In this paper, we present three extensions to the cross-PCF which address these limitations and permit more detailed analysis of multiplex images: topographical correlation maps can visualize local clustering and exclusion between cells; neighbourhood correlation functions can identify colocalization of two or more cell types; and weighted-PCFs describe spatial correlation between points with continuous (rather than discrete) labels. We apply the extended PCFs to synthetic and biological datasets in order to demonstrate the insight that they can generate.

Pathway level subtyping identifies a slow-cycling biological phenotype associated with poor clinical outcomes in colorectal cancer.

Molecular stratification using gene-level transcriptional data has identified subtypes with distinctive genotypic and phenotypic traits, as exemplified by the consensus molecular subtypes (CMS) in colorectal cancer (CRC). Here, rather than gene-level data, we make use of gene ontology and biological activation state information for initial molecular class discovery. In doing so, we defined three pathway-derived subtypes (PDS) in CRC: PDS1 tumors, which are canonical/LGR5+ stem-rich, highly proliferative and display good prognosis; PDS2 tumors, which are regenerative/ANXA1+ stem-rich, with elevated stromal and immune tumor microenvironmental lineages; and PDS3 tumors, which represent a previously overlooked slow-cycling subset of tumors within CMS2 with reduced stem populations and increased differentiated lineages, particularly enterocytes and enteroendocrine cells, yet display the worst prognosis in locally advanced disease. These PDS3 phenotypic traits are evident across numerous bulk and single-cell datasets, and demark a series of subtle biological states that are currently under-represented in pre-clinical models and are not identified using existing subtyping classifiers.

Poor Diagnostic Reproducibility in the Identification of Nonconventional Dysplasia in Colitis Impacts the Application of Histologic Stratification Tools.

Due to their increased cancer risk, patients with longstanding inflammatory bowel disease are offered endoscopic surveillance with concomitant histopathologic assessments, aimed at identifying dysplasia as a precursor lesion of colitis-associated colorectal cancer. However, this strategy is beset with difficulties and limitations. Recently, a novel classification criterion for colitis-associated low-grade dysplasia has been proposed, and an association between nonconventional dysplasia and progression was reported, suggesting the possibility of histology-based stratification of patients with colitis-associated lesions. Here, a cohort of colitis-associated lesions was assessed by a panel of 6 experienced pathologists to test the applicability of the published classification criteria and try and validate the association between nonconventional dysplasia and progression. While confirming the presence of different morphologic patterns of colitis-associated dysplasia, the study demonstrated difficulties concerning diagnostic reproducibility between pathologists and was unable to validate the association of nonconventional dysplasia with cancer progression. Our study highlights the overall difficulty of using histologic assessment of precursor lesions for cancer risk prediction in inflammatory bowel disease patients and suggests the need for a different diagnostic strategy that can objectively identify high-risk phenotypes.

Improving radiotherapy in immunosuppressive microenvironments by targeting complement receptor C5aR1.

An immunosuppressive microenvironment causes poor tumor T cell infiltration and is associated with reduced patient overall survival in colorectal cancer. How to improve treatment responses in these tumors is still a challenge. Using an integrated screening approach to identify cancer-specific vulnerabilities, we identified complement receptor C5aR1 as a druggable target, which when inhibited improved radiotherapy, even in tumors displaying immunosuppressive features and poor CD8+ T cell infiltration. While C5aR1 is well-known for its role in the immune compartment, we found that C5aR1 is also robustly expressed on malignant epithelial cells, highlighting potential tumor cell-specific functions. C5aR1 targeting resulted in increased NF-κB-dependent apoptosis specifically in tumors and not normal tissues, indicating that, in malignant cells, C5aR1 primarily regulated cell fate. Collectively, these data revealed that increased complement gene expression is part of the stress response mounted by irradiated tumors and that targeting C5aR1 could improve radiotherapy, even in tumors displaying immunosuppressive features.

ß-Catenin Drives Butyrophilin-like Molecule Loss and γδ T-cell Exclusion in Colon Cancer.

Intraepithelial lymphocytes (IEL) expressing γδ T-cell receptors (γδTCR) play key roles in elimination of colon cancer. However, the precise mechanisms by which progressing cancer cells evade immunosurveillance by these innate T cells are unknown. Here, we investigated how loss of the Apc tumor suppressor in gut tissue could enable nascent cancer cells to escape immunosurveillance by cytotoxic γδIELs. In contrast with healthy intestinal or colonic tissue, we found that γδIELs were largely absent from the microenvironment of both mouse and human tumors, and that butyrophilin-like (BTNL) molecules, which can critically regulate γδIEL through direct γδTCR interactions, were also downregulated in tumors. We then demonstrated that ß-catenin activation through loss of Apc rapidly suppressed expression of the mRNA encoding the HNF4A and HNF4G transcription factors, preventing their binding to promoter regions of Btnl genes. Reexpression of BTNL1 and BTNL6 in cancer cells increased γδIEL survival and activation in coculture assays but failed to augment their cancer-killing ability in vitro or their recruitment to orthotopic tumors. However, inhibition of ß-catenin signaling via genetic deletion of Bcl9/Bcl9L in either Apc-deficient or mutant ß-catenin mouse models restored Hnf4a, Hnf4g, and Btnl gene expression and γδ T-cell infiltration into tumors. These observations highlight an immune-evasion mechanism specific to WNT-driven colon cancer cells that disrupts γδIEL immunosurveillance and furthers cancer progression.

Epithelial TGFß engages growth-factor signalling to circumvent apoptosis and drive intestinal tumourigenesis with aggressive features.

The pro-tumourigenic role of epithelial TGFß signalling in colorectal cancer (CRC) is controversial. Here, we identify a cohort of born to be bad early-stage (T1) colorectal tumours, with aggressive features and a propensity to disseminate early, that are characterised by high epithelial cell-intrinsic TGFß signalling. In the presence of concurrent Apc and Kras mutations, activation of epithelial TGFß signalling rampantly accelerates tumourigenesis and share transcriptional signatures with those of the born to be bad T1 human tumours and predicts recurrence in stage II CRC. Mechanistically, epithelial TGFß signalling induces a growth-promoting EGFR-signalling module that synergises with mutant APC and KRAS to drive MAPK signalling that re-sensitise tumour cells to MEK and/or EGFR inhibitors. Together, we identify epithelial TGFß signalling both as a determinant of early dissemination and a potential therapeutic vulnerability of CRC's with born to be bad traits.

Dynamic and adaptive cancer stem cell population admixture in colorectal neoplasia.

Intestinal homeostasis is underpinned by LGR5+ve crypt-base columnar stem cells (CBCs), but following injury, dedifferentiation results in the emergence of LGR5-ve regenerative stem cell populations (RSCs), characterized by fetal transcriptional profiles. Neoplasia hijacks regenerative signaling, so we assessed the distribution of CBCs and RSCs in mouse and human intestinal tumors. Using combined molecular-morphological analysis, we demonstrate variable expression of stem cell markers across a range of lesions. The degree of CBC-RSC admixture was associated with both epithelial mutation and microenvironmental signaling disruption and could be mapped across disease molecular subtypes. The CBC-RSC equilibrium was adaptive, with a dynamic response to acute selective pressure, and adaptability was associated with chemoresistance. We propose a fitness landscape model where individual tumors have equilibrated stem cell population distributions along a CBC-RSC phenotypic axis. Cellular plasticity is represented by position shift along this axis and is influenced by cell-intrinsic, extrinsic, and therapeutic selective pressures.

Clinicopathological and molecular characteristics of RSPO fusion-positive colorectal cancer.

BACKGROUND: RSPO fusions that lead to WNT pathway activation are potential therapeutic targets in colorectal cancer (CRC), but their clinicopathological significance remains unclear. METHODS: We screened 1019 CRCs for RSPO fusions using multiplex reverse transcription-PCR. The RSPO fusion-positive tumours were subjected to whole-exome sequencing (WES). RESULTS: Our analysis identified 29 CRCs with RSPO fusions (2.8%), consisting of five with an EIF3E-RSPO2 fusion and 24 with PTPRK-RSPO3 fusions. The patients were 17 women and 12 men. Thirteen tumours (45%) were right-sided. Histologically, approximately half of the tumours (13/29, 45%) had a focal or extensive mucinous component that was significantly more frequent than the RSPO fusion-negative tumours (13%; P = 8.1 × 10-7). Four tumours (14%) were mismatch repair-deficient. WES identified KRAS, BRAF, and NRAS mutations in a total of 27 tumours (93%). In contrast, pathogenic mutations in major WNT pathway genes, such as APC, CTNNB1 and RNF43, were absent. RSPO fusion status did not have a statistically significant influence on the overall or recurrence-free survival. These clinicopathological and genetic features were also confirmed in a pooled analysis of previous studies. CONCLUSION: RSPO fusion-positive CRCs constitute a rare subgroup of CRCs with several characteristic clinicopathological and genetic features.

Activation of innate-adaptive immune machinery by poly(I:C) exposes a therapeutic vulnerability to prevent relapse in stroma-rich colon cancer.

OBJECTIVE: Stroma-rich tumours represent a poor prognostic subtype in stage II/III colon cancer (CC), with high relapse rates and limited response to standard adjuvant chemotherapy. DESIGN: To address the lack of efficacious therapeutic options for patients with stroma-rich CC, we stratified our human tumour cohorts according to stromal content, enabling identification of the biology underpinning relapse and potential therapeutic vulnerabilities specifically within stroma-rich tumours that could be exploited clinically. Following human tumour-based discovery and independent clinical validation, we use a series of in vitro and stroma-rich in vivo models to test and validate the therapeutic potential of elevating the biology associated with reduced relapse in human tumours. RESULTS: By performing our analyses specifically within the stroma-rich/high-fibroblast (HiFi) subtype of CC, we identify and validate the clinical value of a HiFi-specific prognostic signature (HPS), which stratifies tumours based on STAT1-related signalling (High-HPS v Low-HPS=HR 0.093, CI 0.019 to 0.466). Using in silico, in vitro and in vivo models, we demonstrate that the HPS is associated with antigen processing and presentation within discrete immune lineages in stroma-rich CC, downstream of double-stranded RNA and viral response signalling. Treatment with the TLR3 agonist poly(I:C) elevated the HPS signalling and antigen processing phenotype across in vitro and in vivo models. In an in vivo model of stroma-rich CC, poly(I:C) treatment significantly increased systemic cytotoxic T cell activity (p<0.05) and reduced liver metastases (p<0.0002). CONCLUSION: This study reveals new biological insight that offers a novel therapeutic option to reduce relapse rates in patients with the worst prognosis CC.

Germline variant testing in serrated polyposis syndrome.

BACKGROUND AND AIM: Serrated polyposis syndrome (SPS) is now known to be the commonest polyposis syndrome. Previous analyses for germline variants have shown no consistent positive findings. To exclude other polyposis syndromes, 2019 British Society of Gastroenterology (BSG) guidelines advise gene panel testing if the patient is under 50 years, there are multiple affected individuals within a family, or there is dysplasia within any of the polyps. METHODS: A database of SPS patients was established at the Oxford University Hospitals NHS Foundation Trust. Patients were referred for genetic assessment based on personal and family history and patient preference. The majority were tested for a hereditary colorectal cancer panel including MUTYH, APC, PTEN, SMAD4, BMPR1A, STK11, NTLH1, POLD1, POLE, GREM1 (40-kb duplication), PMS2, and Lynch syndrome mismatch repair genes. RESULTS: One hundred and seventy-three patients were diagnosed with SPS based on World Health Organization 2019 criteria between February 2010 and December 2020. The mean age of diagnosis was 54.2 ± 16.8 years. Seventy-three patients underwent genetic testing and 15/73 (20.5%) were found to have germline variants, of which 7/73 (9.6%) had a pathogenic variant (MUTYH n = 2, SMAD4 n = 1, CHEK2 n = 2, POLD1 n = 1, and RNF43 n = 1). Only 60% (9/15) of these patients would have been recommended for gene panel testing according to current BSG guidelines. CONCLUSIONS: A total of 20.5% of SPS patients tested were affected by heterozygous germline variants, including previously unreported associations with CHEK2 and POLD1. This led to a change in management in seven patients (9.6%). Current recommendations may miss SPS associated with germline variants, which is more common than previously anticipated.

Fluctuating methylation clocks for cell lineage tracing at high temporal resolution in human tissues.

Molecular clocks that record cell ancestry mutate too slowly to measure the short-timescale dynamics of cell renewal in adult tissues. Here, we show that fluctuating DNA methylation marks can be used as clocks in cells where ongoing methylation and demethylation cause repeated 'flip-flops' between methylated and unmethylated states. We identify endogenous fluctuating CpG (fCpG) sites using standard methylation arrays and develop a mathematical model to quantitatively measure human adult stem cell dynamics from these data. Small intestinal crypts were inferred to contain slightly more stem cells than the colon, with slower stem cell replacement in the small intestine. Germline APC mutation increased the number of replacements per crypt. In blood, we measured rapid expansion of acute leukemia and slower growth of chronic disease. Thus, the patterns of human somatic cell birth and death are measurable with fluctuating methylation clocks (FMCs).

Morphogen regulation of stem cell plasticity in intestinal regeneration and carcinogenesis.

The intestinal epithelium is a tissue with high cell turnover, supported by adult intestinal stem cells. Intestinal homeostasis is underpinned by crypt basal columnar stem cells, marked by expression of the LGR5 gene. However, recent research has demonstrated considerable stem cell plasticity following injury, with dedifferentiation of a range of other intestinal cell populations, induced by a permissive microenvironment in the regenerating mucosa. The regulation of this profound adaptive cell reprogramming response is the subject of current research. There is a demonstrable contribution from disruption of key homeostatic signaling pathways such as wingless-related integration site and bone morphogenetic protein, and an emerging signaling hub role for the mechanoreceptor transducers Yes-associated protein 1/transcriptional coactivator with PDZ-binding motif, negatively regulated by the Hippo pathway. However, a number of outstanding questions remain, including a need to understand how tissues sense damage, and how pathways intersect to mediate dynamic changes in the stem cell population. Better understanding of these pathways, associated functional redundancies, and how they may be both enhanced for recovery of inflammatory diseases, and co-opted in neoplasia development, may have significant clinical implications, and could lead to development of more targeted molecular therapies which target individual stem or stem-like cell populations.

The Origin and Contribution of Cancer-Associated Fibroblasts in Colorectal Carcinogenesis.

BACKGROUND & AIMS: Cancer-associated fibroblasts (CAFs) play an important role in colorectal cancer (CRC) progression and predict poor prognosis in CRC patients. However, the cellular origins of CAFs remain unknown, making it challenging to therapeutically target these cells. Here, we aimed to identify the origins and contribution of colorectal CAFs associated with poor prognosis. METHODS: To elucidate CAF origins, we used a colitis-associated CRC mouse model in 5 different fate-mapping mouse lines with 5-bromodeoxyuridine dosing. RNA sequencing of fluorescence-activated cell sorting-purified CRC CAFs was performed to identify a potential therapeutic target in CAFs. To examine the prognostic significance of the stromal target, CRC patient RNA sequencing data and tissue microarray were used. CRC organoids were injected into the colons of knockout mice to assess the mechanism by which the stromal gene contributes to colorectal tumorigenesis. RESULTS: Our lineage-tracing studies revealed that in CRC, many ACTA2+ CAFs emerge through proliferation from intestinal pericryptal leptin receptor (Lepr)+ cells. These Lepr-lineage CAFs, in turn, express melanoma cell adhesion molecule (MCAM), a CRC stroma-specific marker that we identified with the use of RNA sequencing. High MCAM expression induced by transforming growth factor ß was inversely associated with patient survival in human CRC. In mice, stromal Mcam knockout attenuated orthotopically injected colorectal tumoroid growth and improved survival through decreased tumor-associated macrophage recruitment. Mechanistically, fibroblast MCAM interacted with interleukin-1 receptor 1 to augment nuclear factor κB-IL34/CCL8 signaling that promotes macrophage chemotaxis. CONCLUSIONS: In colorectal carcinogenesis, pericryptal Lepr-lineage cells proliferate to generate MCAM+ CAFs that shape the tumor-promoting immune microenvironment. Preventing the expansion/differentiation of Lepr-lineage CAFs or inhibiting MCAM activity could be effective therapeutic approaches for CRC.

A Pilot Study on Automatic Three-Dimensional Quantification of Barrett's Esophagus for Risk Stratification and Therapy Monitoring.

BACKGROUND & AIMS: Barrett's epithelium measurement using widely accepted Prague C&M classification is highly operator dependent. We propose a novel methodology for measuring this risk score automatically. The method also enables quantification of the area of Barrett's epithelium (BEA) and islands, which was not possible before. Furthermore, it allows 3-dimensional (3D) reconstruction of the esophageal surface, enabling interactive 3D visualization. We aimed to assess the accuracy of the proposed artificial intelligence system on both phantom and endoscopic patient data. METHODS: Using advanced deep learning, a depth estimator network is used to predict endoscope camera distance from the gastric folds. By segmenting BEA and gastroesophageal junction and projecting them to the estimated mm distances, we measure C&M scores including the BEA. The derived endoscopy artificial intelligence system was tested on a purpose-built 3D printed esophagus phantom with varying BEAs and on 194 high-definition videos from 131 patients with C&M values scored by expert endoscopists. RESULTS: Endoscopic phantom video data demonstrated a 97.2% accuracy with a marginal ± 0.9 mm average deviation for C&M and island measurements, while for BEA we achieved 98.4% accuracy with only ±0.4 cm2 average deviation compared with ground-truth. On patient data, the C&M measurements provided by our system concurred with expert scores with marginal overall relative error (mean difference) of 8% (3.6 mm) and 7% (2.8 mm) for C and M scores, respectively. CONCLUSIONS: The proposed methodology automatically extracts Prague C&M scores with high accuracy. Quantification and 3D reconstruction of the entire Barrett's area provides new opportunities for risk stratification and assessment of therapy response.

Oncogenic BRAF, unrestrained by TGFß-receptor signalling, drives right-sided colonic tumorigenesis.

Right-sided (proximal) colorectal cancer (CRC) has a poor prognosis and a distinct mutational profile, characterized by oncogenic BRAF mutations and aberrations in mismatch repair and TGFß signalling. Here, we describe a mouse model of right-sided colon cancer driven by oncogenic BRAF and loss of epithelial TGFß-receptor signalling. The proximal colonic tumours that develop in this model exhibit a foetal-like progenitor phenotype (Ly6a/Sca1+) and, importantly, lack expression of Lgr5 and its associated intestinal stem cell signature. These features are recapitulated in human BRAF-mutant, right-sided CRCs and represent fundamental differences between left- and right-sided disease. Microbial-driven inflammation supports the initiation and progression of these tumours with foetal-like characteristics, consistent with their predilection for the microbe-rich right colon and their antibiotic sensitivity. While MAPK-pathway activating mutations drive this foetal-like signature via ERK-dependent activation of the transcriptional coactivator YAP, the same foetal-like transcriptional programs are also initiated by inflammation in a MAPK-independent manner. Importantly, in both contexts, epithelial TGFß-receptor signalling is instrumental in suppressing the tumorigenic potential of these foetal-like progenitor cells.

Stromal composition predicts recurrence of early rectal cancer after local excision.

AIMS: After local excision of early rectal cancer, definitive lymph node status is not available. An alternative means for accurate assessment of recurrence risk is required to determine the most appropriate subsequent management. Currently used measures are suboptimal. We assess three measures of tumour stromal content to determine their predictive value after local excision in a well-characterised cohort of rectal cancer patients without prior radiotherapy. METHODS AND RESULTS: A total of 143 patients were included. Haematoxylin and eosin (H&E) sections were scanned for (i) deep neural network (DNN, a machine-learning algorithm) tumour segmentation into compartments including desmoplastic stroma and inflamed stroma; and (ii) digital assessment of tumour stromal fraction (TSR) and optical DNA ploidy analysis. 3' mRNA sequencing was performed to obtain gene expression data from which stromal and immune scores were calculated using the ESTIMATE method. Full results were available for 139 samples and compared with disease-free survival. All three methods were prognostic. Most strongly predictive was a DNN-determined ratio of desmoplastic to inflamed stroma >5.41 (P < 0.0001). A ratio of ESTIMATE stromal to immune score <1.19 was also predictive of disease-free survival (P = 0.00051), as was stromal fraction >36.5% (P = 0.037). CONCLUSIONS: The DNN-determined ratio of desmoplastic to inflamed ratio is a novel and powerful predictor of disease recurrence in locally excised early rectal cancer. It can be assessed on a single H&E section, so could be applied in routine clinical practice to improve the prognostic information available to patients and clinicians to inform the decision concerning further management.

NOTUM from Apc-mutant cells biases clonal competition to initiate cancer.

The tumour suppressor APC is the most commonly mutated gene in colorectal cancer. Loss of Apc in intestinal stem cells drives the formation of adenomas in mice via increased WNT signalling1, but reduced secretion of WNT ligands increases the ability of Apc-mutant intestinal stem cells to colonize a crypt (known as fixation)2. Here we investigated how Apc-mutant cells gain a clonal advantage over wild-type counterparts to achieve fixation. We found that Apc-mutant cells are enriched for transcripts that encode several secreted WNT antagonists, with Notum being the most highly expressed. Conditioned medium from Apc-mutant cells suppressed the growth of wild-type organoids in a NOTUM-dependent manner. Furthermore, NOTUM-secreting Apc-mutant clones actively inhibited the proliferation of surrounding wild-type crypt cells and drove their differentiation, thereby outcompeting crypt cells from the niche. Genetic or pharmacological inhibition of NOTUM abrogated the ability of Apc-mutant cells to expand and form intestinal adenomas. We identify NOTUM as a key mediator during the early stages of mutation fixation that can be targeted to restore wild-type cell competitiveness and provide preventative strategies for people at a high risk of developing colorectal cancer.

Bone Morphogenetic Protein Pathway Antagonism by Grem1 Regulates Epithelial Cell Fate in Intestinal Regeneration.

BACKGROUND & AIMS: In homeostasis, intestinal cell fate is controlled by balanced gradients of morphogen signaling. The bone morphogenetic protein (BMP) pathway has a physiological, prodifferentiation role, predominantly inferred through previous experimental pathway inactivation. Intestinal regeneration is underpinned by dedifferentiation and cell plasticity, but the signaling pathways that regulate this adaptive reprogramming are not well understood. We assessed the BMP signaling landscape and investigated the impact and therapeutic potential of pathway manipulation in homeostasis and regeneration. METHODS: A novel mouse model was generated to assess the effect of the autocrine Bmp4 ligand on individual secretory cell fate. We spatiotemporally mapped BMP signaling in mouse and human regenerating intestine. Transgenic models were used to explore the functional impact of pathway manipulation on stem cell fate and intestinal regeneration. RESULTS: In homeostasis, ligand exposure reduced proliferation, expedited terminal differentiation, abrogated secretory cell survival, and prevented dedifferentiation. After ulceration, physiological attenuation of BMP signaling arose through upregulation of the secreted antagonist Grem1 from topographically distinct populations of fibroblasts. Concomitant expression supported functional compensation after Grem1 deletion from tissue-resident cells. BMP pathway manipulation showed that antagonist-mediated BMP attenuation was obligatory but functionally submaximal, because regeneration was impaired or enhanced by epithelial overexpression of Bmp4 or Grem1, respectively. Mechanistically, Bmp4 abrogated regenerative stem cell reprogramming despite a convergent impact of YAP/TAZ on cell fate in remodeled wounds. CONCLUSIONS: BMP signaling prevents epithelial dedifferentiation, and pathway attenuation through stromal Grem1 upregulation was required for adaptive reprogramming in intestinal regeneration. This intercompartmental antagonism was functionally submaximal, raising the possibility of therapeutic pathway manipulation in inflammatory bowel disease.