Visible and Invisible: Causal Variable Learning and its Application in a Cancer Study

Causal visual discovery is a fundamental yet challenging problem in many research fields. Given visual data and the outcome of interest, the goal is to infer the cause-effect relation. Aside from rich visual ('visible') variables, oftentimes, the outcome is also determined by 'invisible' variables, i.e. the variables from non-visual modalities that do not have visual counterparts. This combination is particularly common in the clinical domain. Built upon the promising invariant causal prediction (ICP) framework, we propose a novel -ICP algorithm to resolve the (visible, invisible) setting. To efficiently discover -plausible causal variables and to estimate the cause-effect relation, the -ICP is learned under a min-min optimisation scheme. Driven by the need for clinical reliability and interpretability, the -ICP is implemented with a typed neural-symbolic functional language. With the built-in program synthesis method, we can synthesize a type-safe program that is comprehensible to the clinical experts. For concept validation of the -ICP, we carefully design a series of synthetic experiments on the type of visual-perception tasks that are encountered in daily life. To further substantiate the proposed method, we demonstrate the application of -ICP on a real-world cancer study dataset, Swiss CRC. This population-based cancer study has spanned over two decades, including 25 fully annotated tissue micro-array (TMA) images with at least resolution and a broad spectrum of clinical meta data for 533 patients. Both the synthetic and clinical experiments demonstrate the advantages of -ICP over the state-of-the-art methods. Finally, we discuss the limitations and challenges to be addressed in the future. Code Of Ethics: I acknowledge that I and all co-authors of this work have read and commit to adhering to the ICLR Code of Ethic

Integrated functional and spatial profiling of tumour immune responses induced by immunotherapy: the iPROFILER platform

Background: Cancer immunotherapy elicits functional activation and changes in immune cell distribution in cancer. Tumour heterogeneity is a reason for treatment failure but is difficult to capture in experimental settings. This proof-of-principle study describes the integrated functional and digital spatial profiling platform iPROFILER to capture in-situ immune activation patterns with high precision. Materials and methods: iPROFILER combines an algorithm-based image analysis approach for spatial profiling with functional analyses of patient-derived tumour fragments (PDTFs). This study utilized a folate receptor 1 (FOLR1)xCD3 bispecific antibody in dual-affinity re-targeting (DART) format as a tool for inducing T-cell responses in patient tumour samples, and an in-depth investigation of the immune perturbations induced in the tumour microenvironment was performed. Results: Ex-vivo DART stimulation induces upregulation of multiple activation markers in CD4+ and CD8+ T-cell populations and secretion of pro-inflammatory cytokines in FOLR1-positive tumour specimens. This response was reduced or absent in tissue samples that did not express FOLR1. Immunological responses were driven by a strong induction of interferon gamma (IFNγ) and IFNγ-induced chemokines suggestive of activation of cytotoxic or Th1-like T cells. Ex-vivo DART treatment led to a numerical increase in effector T cells and an upregulation of immune activation markers in the tumour microenvironment as captured by digital image analysis. Analysis of immune activation in tumour and stromal regions further supported the potential of the platform to measure local differences in cell-type-specific activation patterns. Conclusions: iPROFILER effectively combines functional and spatial readouts to investigate immune responses ex vivo in human tumour samples. Keywords: Ex-vivo models; computer-assisted image processing; digital pathology; immunotherapy; tumour biomarkers; tumour microenvironment

Germline loss-of-function variants in the base-excision repair gene MBD4 cause a Mendelian recessive syndrome of adenomatous colorectal polyposis and acute myeloid leukaemia

Inherited defects in base-excision repair (BER) predispose to adenomatous polyposis and colorectal cancer (CRC), yet our understanding of this important DNA repair pathway remains incomplete. By combining detailed clinical, histological and molecular profiling, we reveal biallelic germline loss-of-function (LOF) variants in the BER gene MBD4 to predispose to adenomatous polyposis and -uniquely amongst CRC predisposition syndromes- to myeloid neoplasms. Neoplasms from MBD4-deficient patients almost exclusively accumulate somatic CpG>TpG mutations, resembling mutational signature SBS1. MBD4-deficient adenomas harbour mutations in known CRC driver genes, although AMER1 mutations were more common and KRAS mutations less frequent. We did not find an increased risk for colorectal tumours in individuals with a monoallelic MBD4 LOF variant. We suggest that this condition should be termed MBD4-associated neoplasia syndrome (MANS) and that MBD4 is included in testing for the genetic diagnosis of polyposis and/or early-onset AM

Molecular phenotyping of colorectal neoplasia shows dynamic and adaptive cancer stem cell population admixture

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), characterised by fetal transcriptional profiles. Neoplasia hijacks regenerative signalling, 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 signalling 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

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

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management. Keywords: 5′-methylcytosine deamination; colorectal cancer; mutational signature; mutator phenotype; polyposi

Tertiary lymphoid structures critical for prognosis in endometrial cancer patients

B-cells play a key role in cancer suppression, particularly when aggregated in tertiary lymphoid structures (TLS). Here, we investigate the role of B-cells and TLS in endometrial cancer (EC). Single cell RNA-sequencing of B-cells shows presence of naïve B-cells, cycling/germinal center B-cells and antibody-secreting cells. Differential gene expression analysis shows association of TLS with L1CAM overexpression. Immunohistochemistry and co-immunofluorescence show L1CAM expression in mature TLS, independent of L1CAM expression in the tumor. Using L1CAM as a marker, 378 of the 411 molecularly classified ECs from the PORTEC-3 biobank are evaluated, TLS are found in 19%. L1CAM expressing TLS are most common in mismatch-repair deficient (29/127, 23%) and polymerase-epsilon mutant EC (24/47, 51%). Multivariable Cox regression analysis shows strong favorable prognostic impact of TLS, independent of clinicopathological and molecular factors. Our data suggests a pivotal role of TLS in outcome of EC patients, and establishes L1CAM as a simple biomarker

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. Keywords: cell plasticity; colorectal cancer; colorectal neoplasia; intestinal polyps; intestinal stem cells; molecular phenotyping; stem cells

An R-based reproducible and user-friendly preprocessing pipeline for CyTOF data

Mass cytometry (CyTOF) has become a method of choice for in-depth characterization of tissue heterogeneity in health and disease, and is currently implemented in multiple clinical trials, where higher quality standards must be met. Currently, preprocessing of raw files is commonly performed in independent standalone tools, which makes it difficult to reproduce. Here, we present an R pipeline based on an updated version of CATALYST that covers all preprocessing steps required for downstream mass cytometry analysis in a fully reproducible way. This new version of CATALYST is based on Bioconductor’s SingleCellExperiment class and fully unit tested. The R-based pipeline includes file concatenation, bead-based normalization, single-cell deconvolution, spillover compensation and live cell gating after debris and doublet removal. Importantly, this pipeline also includes different quality checks to assess machine sensitivity and staining performance while allowing also for batch correction. This pipeline is based on open source R packages and can be easily be adapted to different study designs. It therefore has the potential to significantly facilitate the work of CyTOF users while increasing the quality and reproducibility of data generated with this technology

Loss of Cdx2 Expression in Primary Tumors and Lymph Node Metastases is Specific for Mismatch Repair-Deficiency in Colorectal Cancer

Background: Approximately 20% of all colorectal cancers are hypothesized to arise from the "serrated pathway" characterized by mutation in BRAF, high-level CpG Island Methylator Phenotype, and microsatellite instability/mismatch repair (MMR)-deficiency. MMR-deficient cancers show frequent losses of Cdx2, a homeodomain transcription factor. Here, we determine the predictive value of Cdx2 expression for MMR-deficiency and investigate changes in expression between primary cancers and matched lymph node metastases. Methods: Immunohistochemistry for Cdx2, Mlh1, Msh2, Msh6, and Pms2 was performed on whole tissue sections from 201 patients with primary colorectal cancer and 59 cases of matched lymph node metastases. Receiver operating characteristic curve analysis and Area under the Curve (AUC) were investigated; association of Cdx2 with clinicopathological features and patient survival was carried out. Results: Loss of Cdx2 expression was associated with higher tumor grade (p = 0.0002), advanced pT (p = 0.0166), and perineural invasion (p = 0.0228). Cdx2 loss was an unfavorable prognostic factor in univariate (p = 0.0145) and multivariate [p = 0.0427; HR (95% CI): 0.58 (0.34-0.98)] analysis. The accuracy (AUC) for discriminating MMR-proficient and - deficient cancers was 87% [OR (95% CI): 0.96 (0.95-0.98); p < 0.0001]. Specificity and negative predictive value for MMR-deficiency was 99.1 and 96.3%. One hundred and seventy-four patients had MMR-proficient cancers, of which 60 (34.5%) showed Cdx2 loss. Cdx2 loss in metastases was related to MMR-deficiency (p < 0.0001). There was no difference in expression between primary tumors and matched metastases. Conclusion: Loss of Cdx2 is a sensitive and specific predictor of MMR-deficiency, but is not limited to these tumors, suggesting that events "upstream" of the development of microsatellite instability may impact Cdx2 expression