Ensemble analyses improve signatures of tumour hypoxia and reveal inter-platform differences

BACKGROUND: The reproducibility of transcriptomic biomarkers across datasets remains poor, limiting clinical application. We and others have suggested that this is in-part caused by differential error-structure between datasets, and their incomplete removal by pre-processing algorithms. METHODS: To test this hypothesis, we systematically assessed the effects of pre-processing on biomarker classification using 24 different pre-processing methods and 15 distinct signatures of tumour hypoxia in 10 datasets (2,143 patients). RESULTS: We confirm strong pre-processing effects for all datasets and signatures, and find that these differ between microarray versions. Importantly, exploiting different pre-processing techniques in an ensemble technique improved classification for a majority of signatures. CONCLUSIONS: Assessing biomarkers using an ensemble of pre-processing techniques shows clear value across multiple diseases, datasets and biomarkers. Importantly, ensemble classification improves biomarkers with initially good results but does not result in spuriously improved performance for poor biomarkers. While further research is required, this approach has the potential to become a standard for transcriptomic biomarkers

Multi-scale characterisation of homologous recombination deficiency in breast cancer

BACKGROUND: Homologous recombination is a robust, broadly error-free mechanism of double-strand break repair, and deficiencies lead to PARP inhibitor sensitivity. Patients displaying homologous recombination deficiency can be identified using 'mutational signatures'. However, these patterns are difficult to reliably infer from exome sequencing. Additionally, as mutational signatures are a historical record of mutagenic processes, this limits their utility in describing the current status of a tumour. METHODS: We apply two methods for characterising homologous recombination deficiency in breast cancer to explore the features and heterogeneity associated with this phenotype. We develop a likelihood-based method which leverages small insertions and deletions for high-confidence classification of homologous recombination deficiency for exome-sequenced breast cancers. We then use multinomial elastic net regression modelling to develop a transcriptional signature of heterogeneous homologous recombination deficiency. This signature is then applied to single-cell RNA-sequenced breast cancer cohorts enabling analysis of homologous recombination deficiency heterogeneity and differential patterns of tumour microenvironment interactivity. RESULTS: We demonstrate that the inclusion of indel events, even at low levels, improves homologous recombination deficiency classification. Whilst BRCA-positive homologous recombination deficient samples display strong similarities to those harbouring BRCA1/2 defects, they appear to deviate in microenvironmental features such as hypoxic signalling. We then present a 228-gene transcriptional signature which simultaneously characterises homologous recombination deficiency and BRCA1/2-defect status, and is associated with PARP inhibitor response. Finally, we show that this signature is applicable to single-cell transcriptomics data and predict that these cells present a distinct milieu of interactions with their microenvironment compared to their homologous recombination proficient counterparts, typified by a decreased cancer cell response to TNFα signalling. CONCLUSIONS: We apply multi-scale approaches to characterise homologous recombination deficiency in breast cancer through the development of mutational and transcriptional signatures. We demonstrate how indels can improve homologous recombination deficiency classification in exome-sequenced breast cancers. Additionally, we demonstrate the heterogeneity of homologous recombination deficiency, especially in relation to BRCA1/2-defect status, and show that indications of this feature can be captured at a single-cell level, enabling further investigations into interactions between DNA repair deficient cells and their tumour microenvironment

Evaluation of dormancy states in cancer and associated therapeutic opportunities

Tumour mass dormancy and cancer cell quiescence represent the two facets of cancer dormancy and play key roles in cancer development and progression. Quiescence describes the reversible, proliferative arrest of individual cancer cells that has been observed as a contributing factor of resistance to chemotherapy and other treatments targeting cycling cells. In contrast, tumour mass dormancy describes the state of no net tumour growth, which can arise due to inadequate tumour vascularisation or anti- tumour immune response, during which tumours can acquire additional mutations and establish a microenvironment permissive for growth. Currently, both dormancy states remain poorly characterised. This thesis presents computational frameworks for evaluating the two states and comprehensively profiles their abundance and associated genomic and cellular features across 31 solid cancers from the Cancer Genome Atlas. Using machine learning approaches, I demonstrate that cancer cell quiescence preferentially arises in less mutated tumours with intact TP53 and DNA damage repair pathways. I also highlight novel genomic dependencies, such as CEP89 amplification, which drive an impairment of quiescence. Similarly, mutations within CASP8 and HRAS oncogenes are shown to be enriched and positively selected in samples with tumour mass dormancy. I also highlight an association between APOBEC mutagenesis and both dormancy states. Moreover, tumour mass dormancy is shown to be associated with infiltration with macrophages and cytotoxic and regulatory T cells but a decreased infiltration with Th17 cells. Lastly, using single-cell data, I demonstrate that quiescence underlies resistance to a wide range of therapies, including treatments targeting cell cycle regulation, proliferative kinase signalling and epigenetic regulation. Ultimately, this analysis sheds light on the underlying biology of cancer dormancy states, potentially highlighting vulnerabilities that can be targeted in the clinic. It also provides a transcriptional signature of therapy-tolerant quiescent cells that could be explored further in the clinic to monitor patient therapy response

Leveraging big data resources and data integration in biology: applying computational systems analyses and machine learning to gain insights into the biology of cancers

Recently, many "molecular profiling" projects have yielded vast amounts of genetic, epigenetic, transcription, protein expression, metabolic and drug response data for cancerous tumours, healthy tissues, and cell lines. We aim to facilitate a multi-scale understanding of these high-dimensional biological data and the complexity of the relationships between the different data types taken from human tumours. Further, we intend to identify molecular disease subtypes of various cancers, uncover the subtype-specific drug targets and identify sets of therapeutic molecules that could potentially be used to inhibit these targets. We collected data from over 20 publicly available resources. We then leverage integrative computational systems analyses, network analyses and machine learning, to gain insights into the pathophysiology of pancreatic cancer and 32 other human cancer types. Here, we uncover aberrations in multiple cell signalling and metabolic pathways that implicate regulatory kinases and the Warburg effect as the likely drivers of the distinct molecular signatures of three established pancreatic cancer subtypes. Then, we apply an integrative clustering method to four different types of molecular data to reveal that pancreatic tumours can be segregated into two distinct subtypes. We define sets of proteins, mRNAs, miRNAs and DNA methylation patterns that could serve as biomarkers to accurately differentiate between the two pancreatic cancer subtypes. Then we confirm the biological relevance of the identified biomarkers by showing that these can be used together with pattern-recognition algorithms to infer the drug sensitivity of pancreatic cancer cell lines accurately. Further, we evaluate the alterations of metabolic pathway genes across 32 human cancers. We find that while alterations of metabolic genes are pervasive across all human cancers, the extent of these gene alterations varies between them. Based on these gene alterations, we define two distinct cancer supertypes that tend to be associated with different clinical outcomes and show that these supertypes are likely to respond differently to anticancer drugs. Overall, we show that the time has already arrived where we can leverage available data resources to potentially elicit more precise and personalised cancer therapies that would yield better clinical outcomes at a much lower cost than is currently being achieved

Characterising Heterogeneity of Glioblastoma using Multi-parametric Magnetic Resonance Imaging

A better understanding of tumour heterogeneity is central for accurate diagnosis, targeted therapy and personalised treatment of glioblastoma patients. This thesis aims to investigate whether pre-operative multi-parametric magnetic resonance imaging (MRI) can provide a useful tool for evaluating inter-tumoural and intra-tumoural heterogeneity of glioblastoma. For this purpose, we explored: 1) the utilities of habitat imaging in combining multi-parametric MRI for identifying invasive sub-regions (I & II); 2) the significance of integrating multi-parametric MRI, and extracting modality inter-dependence for patient stratification (III & IV); 3) the value of advanced physiological MRI and radiomics approach in predicting epigenetic phenotypes (V). The following observations were made: I. Using a joint histogram analysis method, habitats with different diffusivity patterns were identified. A non-enhancing sub-region with decreased isotropic diffusion and increased anisotropic diffusion was associated with progression-free survival (PFS, hazard ratio [HR] = 1.08, P < 0.001) and overall survival (OS, HR = 1.36, P < 0.001) in multivariate models. II. Using a thresholding method, two low perfusion compartments were identified, which displayed hypoxic and pro-inflammatory microenvironment. Higher lactate in the low perfusion compartment with restricted diffusion was associated with a worse survival (PFS: HR = 2.995, P = 0.047; OS: HR = 4.974, P = 0.005). III. Using an unsupervised multi-view feature selection and late integration method, two patient subgroups were identified, which demonstrated distinct OS (P = 0.007) and PFS (P < 0.001). Features selected by this approach showed significantly incremental prognostic value for 12-month OS (P = 0.049) and PFS (P = 0.022) than clinical factors. IV. Using a method of unsupervised clustering via copula transform and discrete feature extraction, three patient subgroups were identified. The subtype demonstrating high inter-dependency of diffusion and perfusion displayed higher lactate than the other two subtypes (P = 0.016 and P = 0.044, respectively). Both subtypes of low and high inter-dependency showed worse PFS compared to the intermediate subtype (P = 0.046 and P = 0.009, respectively). V. Using a radiomics approach, advanced physiological images showed better performance than structural images for predicting O6-methylguanine-DNA methyltransferase (MGMT) methylation status. For predicting 12-month PFS, the model of radiomic features and clinical factors outperformed the model of MGMT methylation and clinical factors (P = 0.010). In summary, pre-operative multi-parametric MRI shows potential for the non-invasive evaluation of glioblastoma heterogeneity, which could provide crucial information for patient care.The Cambridge Trust and China Scholarship Council ; Clare College; the British Neuro-Oncology Society; the EG Fearnsides Trust; the International Society for Magnetic Resonance in Medicin

In Vitro and In Vivo Models of Colorectal Cancer for Clinical Application

The Special Issue "In Vitro and In Vivo Models of Colorectal Cancer for Clinical Application", edited by Marta Baiocchi and Ann Zeuner for Cancers, collects original research papers and reviews, depicting the current state and the perspectives of CRC models for preclinical and translational research. Original research papers published in this issue focus on some of the hottest topics in CRC research, such as circulating tumor cells, epigenetic regulation of stemness states, new therapeutic targets, molecular CRC classification and experimental CRC models such as organoids and PDXs. Additionally, four reviews on CRC stem cells, immunotherapy and drug discovery provide an updated viewpoint on key topics linking benchtop to bedside research in CRC