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.

Virtual alignment of pathology image series for multi-gigapixel whole slide images.

Interest in spatial omics is on the rise, but generation of highly multiplexed images remains challenging, due to cost, expertise, methodical constraints, and access to technology. An alternative approach is to register collections of whole slide images (WSI), generating spatially aligned datasets. WSI registration is a two-part problem, the first being the alignment itself and the second the application of transformations to huge multi-gigapixel images. To address both challenges, we developed Virtual Alignment of pathoLogy Image Series (VALIS), software which enables generation of highly multiplexed images by aligning any number of brightfield and/or immunofluorescent WSI, the results of which can be saved in the ome.tiff format. Benchmarking using publicly available datasets indicates VALIS provides state-of-the-art accuracy in WSI registration and 3D reconstruction. Leveraging existing open-source software tools, VALIS is written in Python, providing a free, fast, scalable, robust, and easy-to-use pipeline for registering multi-gigapixel WSI, facilitating downstream spatial analyses.

ß-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.

Metastatic recurrence in colorectal cancer arises from residual EMP1+ cells.

Around 30-40% of patients with colorectal cancer (CRC) undergoing curative resection of the primary tumour will develop metastases in the subsequent years1. Therapies to prevent disease relapse remain an unmet medical need. Here we uncover the identity and features of the residual tumour cells responsible for CRC relapse. An analysis of single-cell transcriptomes of samples from patients with CRC revealed that the majority of genes associated with a poor prognosis are expressed by a unique tumour cell population that we named high-relapse cells (HRCs). We established a human-like mouse model of microsatellite-stable CRC that undergoes metastatic relapse after surgical resection of the primary tumour. Residual HRCs occult in mouse livers after primary CRC surgery gave rise to multiple cell types over time, including LGR5+ stem-like tumour cells2-4, and caused overt metastatic disease. Using Emp1 (encoding epithelial membrane protein 1) as a marker gene for HRCs, we tracked and selectively eliminated this cell population. Genetic ablation of EMP1high cells prevented metastatic recurrence and mice remained disease-free after surgery. We also found that HRC-rich micrometastases were infiltrated with T cells, yet became progressively immune-excluded during outgrowth. Treatment with neoadjuvant immunotherapy eliminated residual metastatic cells and prevented mice from relapsing after surgery. Together, our findings reveal the cell-state dynamics of residual disease in CRC and anticipate that therapies targeting HRCs may help to avoid metastatic relapse.

Homeostasis limits keratinocyte evolution.

Recent studies have revealed that normal human tissues accumulate many somatic mutations. In particular, human skin is riddled with mutations, with multiple subclones of variable sizes. Driver mutations are frequent and tend to have larger subclone sizes, suggesting selection. To begin to understand the histories encoded by these complex somatic mutations, we incorporated genomes into a simple agent-based skin-cell model whose prime directive is homeostasis. In this model, stem-cell survival is random and dependent on proximity to the basement membrane. This simple homeostatic skin model recapitulates the observed log-linear distributions of somatic mutations, where most mutations are found in increasingly smaller subclones that are typically lost with time. Hence, neutral mutations are "passengers" whose fates depend on the random survival of their stem cells, where a rarer larger subclone reflects the survival and spread of mutations acquired earlier in life. The model can also maintain homeostasis and accumulate more frequent and larger driver subclones if these mutations (NOTCH1 and TP53) confer relatively higher persistence in normal skin or during tissue damage (sunlight). Therefore, a relatively simple model of epithelial turnover indicates how observed passenger and driver somatic mutations could accumulate without violating the prime directive of homeostasis in normal human tissues.

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.

The emerging era of personalized medicine in advanced colorectal cancer.

Colorectal cancer (CRC) is a genetically heterogeneous disease with its pathogenesis often driven by varying genetic or epigenetic alterations. This has led to a substantial number of patients developing chemoresistance and treatment failure, resulting in a high mortality rate for advanced disease. Deep molecular analysis has allowed for the discovery of key intestinal signaling pathways which impacts colonic epithelial cell fate, and the integral role of the tumor microenvironment on cancer growth and dissemination. Through transitioning pre-clinical knowledge in research into clinical practice, many potential druggable targets within these pathways have been discovered in the hopes of overcoming the roadblocks encountered by conventional therapies. A personalized approach tailoring treatment according to the histopathological and molecular features of individual tumors can hopefully translate to better patient outcomes, and reduce the rate of recurrence in patients with advanced CRC. Herein, the latest understanding on the molecular science behind CRC tumorigenesis, and the potential treatment targets currently at the forefront of research are summarized.

Biological Misinterpretation of Transcriptional Signatures in Tumor Samples Can Unknowingly Undermine Mechanistic Understanding and Faithful Alignment with Preclinical Data.

PURPOSE: Precise mechanism-based gene expression signatures (GES) have been developed in appropriate in vitro and in vivo model systems, to identify important cancer-related signaling processes. However, some GESs originally developed to represent specific disease processes, primarily with an epithelial cell focus, are being applied to heterogeneous tumor samples where the expression of the genes in the signature may no longer be epithelial-specific. Therefore, unknowingly, even small changes in tumor stroma percentage can directly influence GESs, undermining the intended mechanistic signaling. EXPERIMENTAL DESIGN: Using colorectal cancer as an exemplar, we deployed numerous orthogonal profiling methodologies, including laser capture microdissection, flow cytometry, bulk and multiregional biopsy clinical samples, single-cell RNA sequencing and finally spatial transcriptomics, to perform a comprehensive assessment of the potential for the most widely used GESs to be influenced, or confounded, by stromal content in tumor tissue. To complement this work, we generated a freely-available resource, ConfoundR; https://confoundr.qub.ac.uk/, that enables users to test the extent of stromal influence on an unlimited number of the genes/signatures simultaneously across colorectal, breast, pancreatic, ovarian and prostate cancer datasets. RESULTS: Findings presented here demonstrate the clear potential for misinterpretation of the meaning of GESs, due to widespread stromal influences, which in-turn can undermine faithful alignment between clinical samples and preclinical data/models, particularly cell lines and organoids, or tumor models not fully recapitulating the stromal and immune microenvironment. CONCLUSIONS: Efforts to faithfully align preclinical models of disease using phenotypically-designed GESs must ensure that the signatures themselves remain representative of the same biology when applied to clinical samples.

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.

Immunosuppressive niche engineering at the onset of human colorectal cancer.

The evolutionary dynamics of tumor initiation remain undetermined, and the interplay between neoplastic cells and the immune system is hypothesized to be critical in transformation. Colorectal cancer (CRC) presents a unique opportunity to study the transition to malignancy as pre-cancers (adenomas) and early-stage cancers are frequently resected. Here, we examine tumor-immune eco-evolutionary dynamics from pre-cancer to carcinoma using a computational model, ecological analysis of digital pathology data, and neoantigen prediction in 62 patient samples. Modeling predicted recruitment of immunosuppressive cells would be the most common driver of transformation. As predicted, ecological analysis reveals that progressed adenomas co-localized with immunosuppressive cells and cytokines, while benign adenomas co-localized with a mixed immune response. Carcinomas converge to a common immune "cold" ecology, relaxing selection against immunogenicity and high neoantigen burdens, with little evidence for PD-L1 overexpression driving tumor initiation. These findings suggest re-engineering the immunosuppressive niche may prove an effective immunotherapy in CRC.

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.

Gattaca: Base-Pair Resolution Mutation Tracking for Somatic Evolution Studies using Agent-based Models.

Research over the past two decades has made substantial inroads into our understanding of somatic mutations. Recently, these studies have focused on understanding their presence in homeostatic tissue. In parallel, agent-based mechanistic models have emerged as an important tool for understanding somatic mutation in tissue; yet no common methodology currently exists to provide base-pair resolution data for these models. Here, we present Gattaca as the first method for introducing and tracking somatic mutations at the base-pair resolution within agent-based models that typically lack nuclei. With nuclei that incorporate human reference genomes, mutational context, and sequence coverage/error information, Gattaca is able to realistically evolve sequence data, facilitating comparisons between in silico cell tissue modeling with experimental human somatic mutation data. This user-friendly method, incorporated into each in silico cell, allows us to fully capture somatic mutation spectra and evolution.

Molecular Subtyping Resource: a user-friendly tool for rapid biological discovery from transcriptional data.

Generation of transcriptional data has dramatically increased in the past decade, driving the development of analytical algorithms that enable interrogation of the biology underpinning the profiled samples. However, these resources require users to have expertise in data wrangling and analytics, reducing opportunities for biological discovery by 'wet-lab' users with a limited programming skillset. Although commercial solutions exist, costs for software access can be prohibitive for academic research groups. To address these challenges, we have developed an open source and user-friendly data analysis platform for on-the-fly bioinformatic interrogation of transcriptional data derived from human or mouse tissue, called Molecular Subtyping Resource (MouSR). This internet-accessible analytical tool, https://mousr.qub.ac.uk/, enables users to easily interrogate their data using an intuitive 'point-and-click' interface, which includes a suite of molecular characterisation options including quality control, differential gene expression, gene set enrichment and microenvironmental cell population analyses from RNA sequencing. The MouSR online tool provides a unique freely available option for users to perform rapid transcriptomic analyses and comprehensive interrogation of the signalling underpinning transcriptional datasets, which alleviates a major bottleneck for biological discovery. This article has an associated First Person interview with the first author of the paper.

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.