Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather ‘plastic’. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

CRISPR guide RNA libraries have been iteratively improved to provide increasingly efficient reagents, although their large size is a barrier for many applications. We design an optimised minimal genome-wide human CRISPR-Cas9 library (MinLibCas9) by mining existing large-scale gene loss-of-function datasets, resulting in a greater than 42% reduction in size compared to other CRISPR-Cas9 libraries while preserving assay sensitivity and specificity. MinLibCas9 provides backward compatibility with existing datasets, increases the dynamic range of CRISPR-Cas9 screens and extends their application to complex models and assays.

Background CRISPR-Cas9 genome-wide screens are being increasingly performed, allowing systematic explorations of cancer dependencies at unprecedented accuracy and scale. One of the major computational challenges when analysing data derived from such screens is to identify genes that are essential for cell survival invariantly across tissues, conditions, and genomic-contexts (core-fitness genes), and to distinguish them from context-specific essential genes. This is of paramount importance to assess the safety profile of candidate therapeutic targets and for elucidating mechanisms involved in tissue-specific genetic diseases. Results We have developed CoRe: an R package implementing existing and novel methods for the identification of core-fitness genes (at two different level of stringency) from joint analyses of multiple CRISPR-Cas9 screens. We demonstrate, through a fully reproducible benchmarking pipeline, that CoRe outperforms state-of-the-art tools, yielding more reliable and biologically relevant sets of core-fitness genes. Conclusions CoRe offers a flexible pipeline, compatible with many pre-processing methods for the analysis of CRISPR data, which can be tailored onto different use-cases. The CoRe package can be used for the identification of high-confidence novel core-fitness genes, as well as a means to filter out potentially cytotoxic hits while analysing cancer dependency datasets for identifying and prioritising novel selective therapeutic targets.

CRISPR-Cas9 recessive genome-wide pooled screens have allowed systematic explorations of weaknesses and vulnerabilities existing in cancer cells, across different tissue lineages at unprecedented accuracy and scale. The identification of novel genes essential for selective cancer cell survival is currently one of the main applications of this technology. Towards this aim, distinguishing genes that are constitutively essential (invariantly across tissues and genomic contexts, i.e. core-fitness genes) from those whose essentiality is associated with molecular features peculiar to certain cancers is of paramount importance for identifying new oncology therapeutic targets. This is crucial to assess the risk of a candidate target's suppression impacting critical cellular processes that are unspecific to cancer. On the other hand, identifying new human core-fitness genes might also elucidate new mechanisms involved in tissue-specific genetic diseases. We present CoRe: an open-source R package implementing established and novel methods for the identification of core-fitness genes based on joint analyses of data from multiple CRISPR-Cas9 screens. In addition, we present results from a fully reproducible benchmarking pipeline demonstrating that CoRe outperforms other state-of-the-art methods, and it yields more reliable sets of core-fitness and common-essential genes with respect to existing reference sets and methods.

Cancer genomic medicine relies on targeting driver genes. However, current catalogues of cancer drivers are mostly based on indirect measurements of mutation frequencies, positions or types, rather than their effect on clonal expansions in vivo. Moreover, nongenetic drivers are largely unknown, as are the epigenetic and transcriptomic effects of genetic drivers. Here we perform spatial computational inference on multiomic data with matched whole-genome sequencing, ATAC-seq and RNA-seq. Using 436 samples, we directly quantify the contribution, or lack thereof, of putative driver genes to subclonal expansions in vivo in 30 colorectal carcinomas (4-33 samples per patient, median=15). Although subclonal neutral evolution was widespread (13/26 cases with sufficient data), there were cases with clear evidence of subclonal selection (6/26) in which we measured epigenetic and transcriptomic differences between subclones in vivo. In 7/26 cases we could not distinguish between neutral or selective evolution with the available data. We identified expanding subclones that were not driven by known genetic alterations, and propose candidate epigenetic drivers. We identified the distinguishing patterns of genomic heterogeneity produced in fast, exponentially growing tumours (7/26) versus neoplasms growing only at the periphery (19/26), as well as identifying clonally intermixed (16/28 cases with sufficient data) versus segregated malignancies (10/28). Our model-based approach measures genetic and non-genetic subclonal selection, or lack thereof, in space and time and allows in vivo comparisons of the emergent phenotypic properties of subclones within human tumours.