Marc J Williams, Benjamin Werner, Christina Curtis, Chris P Barnes, Andrea Sottoriva, Trevor A. Graham

Recent studies have identified prevalent subclonal architectures within many cancer types. However, the temporal evolutionary dynamics that produce these subclonal architectures remain unknown. Here we measure evolutionary dynamics in primary human cancers using computational modelling of clonal selection applied to high throughput sequencing data. Our approach simultaneously determines the subclonal architecture of a tumour sample, and measures the mutation rate, the selective advantage, and the time of appearance of subclones. Simulations demonstrate the accuracy of the method, and revealed the degree to which evolutionary dynamics are recorded in the genome. Application of our method to high-depth sequencing data from gastric and lung cancers revealed that detectable subclones consistently emerged early during tumour growth and had considerably large fitness advantages (>20% growth advantage). Our quantitative platform provides new insight into the evolutionary history of cancers by facilitating the measurement of fundamental evolutionary parameters in individual patients.

Fredrik Salmen, Sanja Vickovic, Ludvig Larsson, Linnea Stenbeck, Johan Vallon-Christersson, Anna Ehinger, Jari Hakkinen, Åke Borg, Jonas Frisen, Patrik Stahl, Joakim Lundeberg

The comprehensive analysis of tumor tissue heterogeneity is crucial for determining specific disease states and establishing suitable treatment regimes. Here, we analyze tumor tissue sections from ten patients diagnosed with HER2+ breast cancer. We obtain and analyze multidimensional, genome-wide transcriptomics data to resolve spatial immune cell distribution and identity within the tissue sections. Furthermore, we determine the extent of immune cell infiltration in different regions of the tumor tissue, including invasive cancer regions. We combine cross-sectioning and computational alignment to build three-dimensional images of the transcriptional landscape of the tumor and its microenvironment. The three-dimensional data clearly demonstrates the heterogeneous nature of tumor-immune interactions and reveal interpatient differences in immune cell infiltration patterns. Our study shows the potential for an improved stratification and description of the tumor-immune interplay, which is likely to be essential in treatment decisions.

Michele Olivieri, Tiffany Cho, Alejandro Álvarez-Quilón, Kejiao Li, Matthew J Schellenberg, Michal Zimmermann, Nicole Hustedt, Silvia Emma Rossi, Salomé Adam, Henrique Melo, Anne Margriet Heijink, Guillermo Sastre-Moreno, Nathalie Moatti, Rachel Szilard, Andrea McEwan, Alexanda K Ling, Almudena Serrano-Benitez, Tajinder Ubhi, Irene Delgado-Sainz, Michael W Ferguson, Grant W Brown, Felipe Cortés-Ledesma, R. Scott Williams, Alberto Martin, Dongyi Xu, Daniel Durocher

The response to DNA damage is critical for cellular homeostasis, tumor suppression, immunity and gametogenesis. In order to provide an unbiased and global view of the DNA damage response in human cells, we undertook 28 CRISPR/Cas9 screens against 25 genotoxic agents in the retinal pigment epithelium-1 (RPE1) cell line. These screens identified 840 genes whose loss causes either sensitivity or resistance to DNA damaging agents. Mining this dataset, we uncovered that ERCC6L2, which is mutated in a bone-marrow failure syndrome, codes for a canonical non-homologous end-joining pathway factor; that the RNA polymerase II component ELOF1 modulates the response to transcription-blocking agents and that the cytotoxicity of the G-quadruplex ligand pyridostatin involves trapping topoisomerase II on DNA. This map of the DNA damage response provides a rich resource to study this fundamental cellular system and has implications for the development and use of genotoxic agents in cancer therapy.

Yotam E Bar-Ephraim, Kai Kretzschmar, Priyanca Asra, Evelien de Jongh, Kim E Boonekamp, Jarno Drost, Joost van Gorp, Apollo Pronk, Niels Smakman, Inez J. Gan, Zsolt Sebestyen, Jurgen HE Kuball, Robert GJ Vries, Hans Clevers

Immune escape has been recognised as one of the hallmarks of cancer. Overcoming this immunomodulatory process by tumour cells has become a major therapeutic target. Here we utilize organoid technology to study immune-cancer interactions and assess immunomodulation by colorectal cancer (CRC). Transcriptional profiling and flow cytometry revealed that organoids maintain differential expression of immunomodulatory molecules present in primary tumours. Finally, we established a method to model antigen-specific epithelial cell killing and cancer immunomodulation in vitro using CRC organoids co-cultured with cytotoxic T cells. Our method may serve as a first step to rebuilding the tumor microenvironment in vitro.

Zachary Steinhart, Traver Hart, Megha Chandrashekhar, Zvezdan Pavlovic, Mélanie Robitaille, Xiaowei Wang, Jarrett Adams, James Pan, Sachdev Sidhu, Jason Moffat, Stephane Angers

CRISPR-Cas9 genome editing enables high-resolution detection of genetic vulnerabilities of cancer cells. We conducted a genome-wide CRISPR-Cas9 screen in RNF43 mutant pancreatic ductal adenocarcinoma (PDAC) cells, which rely on Wnt signaling for proliferation, and discovered a unique requirement for a WNT7B-FZD5 signaling circuit. Our results highlight an underappreciated level of functional specificity at the ligand-receptor level. We derived a panel of recombinant antibodies that reports the expression of nine out of ten human Frizzled receptors and confirm that WNT7B-FZD5 functional specificity cannot be explained by protein expression patterns. We developed two human antibodies that target FZD5 and robustly inhibited the growth of RNF43 mutant PDAC cells grown in vitro and as xenografts, providing strong orthogonal support for the functional specificity observed genetically. Proliferation of a patient-derived PDAC cell line harboring a RNF43 variant previously associated with PDAC was also selectively inhibited by the FZD5 antibodies, further demonstrating their use as a potential targeted therapy.

Ron S Gejman, Heather F Jones, Martin G Klatt, Aaron Y Chang, Claire Y Oh, Smita S Chandran, Tatiana Korontsvit, Viktoriya Zakahleva, Tao Dao, Christopher A Klebanoff, David A Scheinberg

T cell receptor (TCR)-based therapeutic cells and agents have emerged as a new class of effective cancer therapeutics. These therapies work on cells that express intracellular cancer-associated proteins by targeting peptides displayed on major histocompatibility complex receptors. However, cross-reactivities of these agents to off-target cells and tissues have resulted in serious, sometimes fatal, adverse events. We have developed a high throughput genetic platform (termed "PresentER") that encodes MHC-I peptide minigenes for functional immunological assays as well as for determining the reactivities of TCR-like therapeutic agents against large libraries of MHC-I ligands. In this report, we demonstrate that PresentER can be used to identify the on-and-off targets of T cells and TCR mimic antibodies using in vitro co-culture assays or binding assays. We find dozens of MHC-I ligands that are cross-reactive with two TCR mimic antibodies and two native TCRs and that are not easily predictable by other methods.

Alexandre Coudray, Anna M Battenhouse, Philipp Bucher, Vishwanath R Iyer

To detect functional somatic mutations in tumor samples, whole-exome sequencing (WES) is often used for its reliability and relative low cost. RNA-seq, while generally used to measure gene expression, can potentially also be used for identification of somatic mutations. However there has been little systematic evaluation of the utility of RNA-seq for identifying somatic mutations. Here, we develop and evaluate a pipeline for processing RNA-seq data from glioblastoma multiforme (GBM) tumors in order to identify somatic mutations. The pipeline entails the use of the STAR aligner 2-pass procedure jointly with MuTect2 from GATK to detect somatic variants. Variants identified from RNA-seq data were evaluated by comparison against the COSMIC and dbSNP databases, and also compared to somatic variants identified by exome sequencing. We also estimated the putative functional impact of coding variants in the most frequently mutated genes in GBM. Interestingly, variants identified by RNA-seq alone showed better representation of GBM-related mutations cataloged by COSMIC. RNA-seq-only data substantially outperformed the ability of WES to reveal potentially new somatic mutations in known GBM-related pathways, and allowed us to build a high-quality set of somatic mutations common to exome and RNA-seq calls. Using RNA-seq data in parallel with WES data to detect somatic mutations in cancer genomes can thus broaden the scope of discoveries and lend additional support to somatic variants identified by exome sequencing alone.

Peter Ulz, Samantha Perakis, Qing Zhou, Tina Moser, Jelena Belic, Isaac Lazzeri, Albert Wölfler, Armin Zebisch, Armin Gerger, Gunda Pristauz, Edgar Petru, Brandon White, Charles E.S. Roberts, John St. John, Michael G Schimek, Jochen B Geigl, Thomas Bauernhofer, Heinz Sill, Christoph Bock, Ellen Heitzer, Michael R Speicher

Deregulation of transcription factors (TFs) is an important driver of tumorigenesis. We developed and validated a minimally invasive method for assessing TF activity based on cell-free DNA sequencing and nucleosome footprint analysis. We analyzed whole genome sequencing data for >1,000 cell-free DNA samples from cancer patients and healthy controls using a newly developed bioinformatics pipeline that infers accessibility of TF binding sites from cell-free DNA fragmentation patterns. We observed patient-specific as well as tumor-specific patterns, including accurate prediction of tumor subtypes in prostate cancer, with important clinical implications for the management of patients. Furthermore, we show that cell-free DNA TF profiling is capable of early detection of colorectal carcinomas. Our approach for mapping tumor-specific transcription factor binding in vivo based on blood samples makes a key part of the noncoding genome amenable to clinical analysis.

Charles P. Couturier, Shamini Ayyadhury, Phuong U. Le, Jean Monlong, Gabriele Riva, Redouane Allache, Salma Baig, Xiaohua Yan, Mathieu Bourgey, Changseok Lee, Yu Chang David Wang, V. Wee Yong, Marie-Christine Guiot, Bratislav Misic, Jack Antel, Guillaume Bourque, Jiannis Ragoussis, Kevin Petrecca

Cancer stem cells are critical for cancer initiation, development, and resistance to treatments. Our understanding of these processes, and how they relate to glioblastoma heterogeneity, is limited. To overcome these limitations, we performed single-cell RNA-sequencing on 38 296 glioblastoma cells and 22 637 normal human fetal brain cells. Using an unbiased approach, we mapped the lineage hierarchy of the developing human brain and compared the transcriptome of each cancer cell to this roadmap. We discovered a conserved neural trilineage cancer hierarchy with glial progenitor-like cells at the apex. We also found that this progenitor population contains the majority of cancer's cycling cells and is the origin of heterogeneity. Finally, we show that this hierarchal map can be used to identify therapeutic targets specific to progenitor cancer stem cells. Our analyses show that normal brain development reconciles glioblastoma development, unravels the origin of glioblastoma heterogeneity, and helps to identify cancer stem cell-specific targets.

Nathan Wan, David Weinberg, Tzu-Yu Liu, Katherine Niehaus, Daniel Delubac, Ajay Kannan, Brandon White, Eric A. Ariazi, Mitch Bailey, Marvin Bertin, Nathan Boley, Derek Bowen, James Cregg, Adam M. Drake, Riley Ennis, Signe Fransen, Erik Gafni, Loren Hansen, Yaping Liu, Gabriel L Otte, Jennifer Pecson, Brandon Rice, Gabriel E Sanderson, Aarushi Sharma, John St. John, Catherina Tang, Abraham Tzou, Leilani Young, Girish Putcha, Imran S. Haque

Background: Blood-based methods using cell-free DNA (cfDNA) are under development as an alternative to existing screening tests. However, early-stage detection of cancer using tumor-derived cfDNA has proven challenging because of the small proportion of cfDNA derived from tumor tissue in early-stage disease. A machine learning approach to discover signatures in cfDNA, potentially reflective of both tumor and non-tumor contributions, may represent a promising direction for the early detection of cancer. Methods: Whole-genome sequencing was performed on cfDNA extracted from plasma samples (N=546 colorectal cancer and 271 non-cancer controls). Reads aligning to protein-coding gene bodies were extracted, and read counts were normalized. cfDNA tumor fraction was estimated using IchorCNA. Machine learning models were trained using k-fold cross-validation and confounder-based cross-validation to assess generalization performance. Results: In a colorectal cancer cohort heavily weighted towards early-stage cancer (80% stage I/II), we achieved a mean AUC of 0.92 (95% CI 0.91-0.93) with a mean sensitivity of 85% (95% CI 83-86%) at 85% specificity. Sensitivity generally increased with tumor stage and increasing tumor fraction. Stratification by age, sequencing batch, and institution demonstrated the impact of these confounders and provided a more accurate assessment of generalization performance. Conclusions: A machine learning approach using cfDNA achieved high sensitivity and specificity in a large, predominantly early-stage, colorectal cancer cohort. The possibility of systematic technical and institution-specific biases warrants similar confounder analyses in other studies. Prospective validation of this machine learning method and evaluation of a multi-analyte approach are underway.

Luiza Moore, Daniel Leongamornlert, Tim H. H. Coorens, Mathijs A Sanders, Peter Ellis, Kevin Dawson, Francesco Maura, Jyoti Nangalia, Patrick S Tarpey, Simon F Brunner, Henry Lee-Six, Raheleh Rahbari, Sarah Moody, Yvette Hooks, Krishnaa Mahbubani, Mercedes Jimenez-Linan, Jan J Brosens, Christine A. Iacobuzio-Donahue, Inigo Martincorena, Kourosh Saeb-Parsy, Peter J. Campbell, Michael R. Stratton

All normal somatic cells are thought to acquire mutations. However, characterisation of the patterns and consequences of somatic mutation in normal tissues is limited. Uterine endometrium is a dynamic tissue that undergoes cyclical shedding and reconstitution and is lined by a gland-forming epithelium. Whole genome sequencing of normal endometrial glands showed that most are clonal cell populations derived from a recent common ancestor with mutation burdens differing from other normal cell types and manyfold lower than endometrial cancers. Mutational signatures found ubiquitously account for most mutations. Many, in some women potentially all, endometrial glands are colonised by cell clones carrying driver mutations in cancer genes, often with multiple drivers. Total and driver mutation burdens increase with age but are also influenced by other factors including body mass index and parity. Clones with drivers often originate during early decades of life. The somatic mutational landscapes of normal cells differ between cell types and are revealing the procession of neoplastic change leading to cancer.

Steven M. Corsello, Rohith T Nagari, Ryan D Spangler, Jordan Rossen, Mustafa Kocak, Jordan G Bryan, Ranad Humeidi, David Peck, Xiaoyun Wu, Andrew A Tang, Vickie M Wang, Samantha A Bender, Evan Lemire, Rajiv Narayan, Philip Montgomery, Uri Ben-David, Yejia Chen, Matthew G Rees, Nicholas J. Lyons, James M McFarland, Bang T Wong, Li Wang, Nancy Dumont, Patrick J. O’Hearn, Eric Stefan, John G. Doench, Heidi Greulich, Matthew Meyerson, Francisca Vazquez, Aravind Subramanian, Jennifer A Roth, Joshua A. Bittker, Jesse S Boehm, Christopher C Mader, Aviad Tsherniak, Todd R. Golub

Anti-cancer uses of non-oncology drugs have been found on occasion, but such discoveries have been serendipitous and rare. We sought to create a public resource containing the growth inhibitory activity of 4,518 drugs tested across 578 human cancer cell lines. To accomplish this, we used PRISM, which involves drug treatment of molecularly barcoded cell lines in pools. Relative barcode abundance following treatment thus reflects viability of each cell line. We found that an unexpectedly large number of non-oncology drugs selectively inhibited subsets of cancer cell lines. Moreover, the killing activity of the majority of these drugs was predictable based on the molecular features of the cell lines. Follow-up of several of these compounds revealed novel mechanisms. For example, compounds that kill by inducing PDE3A-SLFN12 complex formation; vanadium-containing compounds whose killing is dependent on the sulfate transporter SLC26A2; the alcohol dependence drug disulfiram, which kills cells with low expression of metallothioneins; and the anti-inflammatory drug tepoxalin, whose killing is dependent on high expression of the multi-drug resistance gene ABCB1. These results illustrate the potential of the PRISM drug repurposing resource as a starting point for new oncology therapeutic development. The resource is available at https://depmap.org.

Michael K Dame, Durga Attili, Shannon D McClintock, Priya H Dedhia, Peter Ouilette, Olaf Hardt, Alana M Chin, Xiang Xue, Julie Laliberte, Erica L Katz, Gina M Newsome, David R. Hill, Alyssa J. Miller, Yu-Hwai Tsai, David Agorku, Christopher H Altheim, Andreas Bosio, Becky Simon, Linda C Samuelson, Jay A Stoerker, Henry D Appelman, James Varani, Max S Wicha, Dean E Brenner, Yatrik M Shah, Jason R Spence, Justin A Colacino

The intestine is maintained by stem cells located at the base of crypts and distinguished by the expression of LGR5. Genetically engineered mouse models have provided a wealth of information about intestinal stem cells, while less is known about human intestinal stem cells due to difficulty detecting and isolating these cells. We established an organoid repository from patient-derived adenomas, adenocarcinomas, and normal colon, which we analyzed for variants in 71 colorectal cancer (CRC) associated genes. Normal and neoplastic colon tissue organoids were analyzed by immunohistochemistry and fluorescent-activated cell sorting for LGR5. LGR5-positive cells were isolated from 4 adenoma organoid lines and were subjected to RNA-sequencing. We found that LGR5 expression in the epithelium and stroma was associated with tumor stage, and by integrating functional experiments with LGR5-sorted cell RNA-seq data from adenoma and normal organoids, we found correlations between LGR5 and CRC-specific genes, including DKK4 (dickkopf WNT signaling pathway inhibitor 4) and SMOC2 (SPARC related modular calcium binding 2). Collectively, this work provides resources, methods and new markers to isolate and study stem cells in human tissue homeostasis and carcinogenesis.

Yu Kong, Chris Rose, Ashley A Cass, Martine Darwish, Steve Lianoglou, Pete M Haverty, Ann-Jay Tong, Craig Blanchette, Ira Mellman, Richard Bourgon, John Greally, Suchit Jhunjhunwala, Matthew L. Albert, Haiyin Chen-Harris

Profound loss of DNA methylation is a well-recognized hallmark of cancer. Given its role in silencing transposable elements (TEs), we hypothesized that extensive TE expression occurs in tumors with highly demethylated DNA. We developed REdiscoverTE, a computational method for quantifying genome-wide TE expression in RNA sequencing data. Using The Cancer Genome Atlas database, we observed increased expression of over 400 TE subfamilies, of which 262 appeared to result from a proximal loss of DNA methylation. The most recurrent TEs were among the evolutionarily youngest in the genome, predominantly expressed from intergenic loci, and associated with antiviral or DNA damage responses. Treatment of glioblastoma cells with a demethylation agent resulted in both increased TE expression and de novo presentation of TE-derived peptides on MHC class I molecules. Therapeutic reactivation of tumor-specific TEs may synergize with immunotherapy by inducing both inflammation and the display of potentially immunogenic neoantigens.

Alicia C Smart, Claire A Margolis, Harold Pimentel, Meng Xiao He, Diana Miao, Dennis Adeegbe, Tim Fugmann, Kwok-Kin Wong, Eliezer M Van Allen

Personalized cancer vaccine strategies directed at tumor neoantigens derived from somatic mutations in the DNA are currently under prospective evaluation. Alterations in tumor RNA, rather than DNA, may also represent a previously-unexplored source of neoantigens. Here, we show that intron retention, a widespread feature of cancer transcriptomes, represents a novel source of tumor neoantigens. We developed an in silico approach to identify retained intron neoantigens from RNA sequencing data and applied this methodology to tumor samples from patients with melanoma treated with immune checkpoint blockade, discovering that the retained intron neoantigen burden in these samples augments the DNA-derived, somatic neoantigen burden. We validated the existence of retained intron derived neoantigens by implementing this technique on cancer cell lines with mass spectrometry-derived immunopeptidome data, revealing that retained intron neoantigens were complexed with MHC I experimentally. Unexpectedly, we observed a trend toward lack of clinical benefit from immune checkpoint blockade in high retained intron load-tumors, which harbored transcriptional signatures consistent with cell cycle dysregulation and DNA damage repair. Our results demonstrate the contribution of transcriptional dysregulation to the overall burden of tumor neoantigens, provide a foundation for augmenting personalized cancer vaccine development with a new class of tumor neoantigens, and demonstrate how global transcriptional dysregulation may impact selective response to immune checkpoint blockade.

Lina Wadi, Liis Uusküla-Reimand, Keren Isaev, Shimin Shuai, Vincent Huang, Minggao Liang, J. Drew Thompson, Yao Li, Luyao Ruan, Marta Paczkowska, Michal Krassowski, Irakli Dzneladze, Ken Kron, Alexander Murison, Parisa Mazrooei, Robert G Bristow, Jared T Simpson, Mathieu Lupien, Michael D Wilson, Lincoln D. Stein, Paul C. Boutros, Jüri Reimand

A comprehensive catalogue of the mutations that drive tumorigenesis and progression is essential to understanding tumor biology and developing therapies. Protein-coding driver mutations have been well-characterized by large exome-sequencing studies, however many tumors have no mutations in protein-coding driver genes. Non-coding mutations are thought to explain many of these cases, however few non-coding drivers besides TERT promoter are known. To fill this gap, we analyzed 150,000 cis-regulatory regions in 1,844 whole cancer genomes from the ICGC-TCGA PCAWG project. Using our new method, ActiveDriverWGS, we found 41 frequently mutated regulatory elements (FMREs) enriched in non-coding SNVs and indels (FDR<0.05) characterized by aging-associated mutation signatures and frequent structural variants. Most FMREs are distal from genes, reported here for the first time and also recovered by additional driver discovery methods. FMREs were enriched in super-enhancers, H3K27ac enhancer marks of primary tumors and long-range chromatin interactions, suggesting that the mutations drive cancer by distally controlling gene expression through three-dimensional genome organization. In support of this hypothesis, the chromatin interaction network of FMREs and target genes revealed associations of mutations and differential gene expression of known and novel cancer genes (e.g., CNNB1IP1, RCC1), activation of immune response pathways and altered enhancer marks. Thus distal genomic regions may include additional, infrequently mutated drivers that act on target genes via chromatin loops. Our study is an important step towards finding such regulatory regions and deciphering the somatic mutation landscape of the non-coding genome.

Gabriela S Kinker, Alissa C Greenwald, Rotem Tal, Zhanna Orlova, Michael S Cuoco, James M McFarland, Allison Warren, Christopher Rodman, Jennifer A Roth, Samantha A Bender, Bhavna Kumar, James W. Rocco, Pedro ACM Fernandes, Christopher C Mader, Hadas Keren-Shaul, Alexander Plotnikov, Haim Barr, Aviad Tsherniak, Orit Rozenblatt-Rosen, Valery Krizhanovsky, Sidharth V Puram, Aviv Regev, Itay Tirosh

Cultured cell lines are the workhorse of cancer research, but it is unclear to what extent they recapitulate the cellular heterogeneity observed among malignant cells in tumors, given the absence of a native tumor microenvironment. Here, we used multiplexed single cell RNA-Seq to profile ~200 cancer cell lines. We uncovered expression programs that are recurrently heterogeneous within many cancer cell lines and are largely independent of observed genetic diversity. These programs of heterogeneity are associated with diverse biological processes, including cell cycle, senescence, stress and interferon responses, epithelial-to-mesenchymal transition (EMT), and protein maturation and degradation. Notably, some of these recurrent programs recapitulate those seen in human tumors, suggesting a prominent role of intrinsic plasticity in generating intra-tumoral heterogeneity. Moreover, the data allowed us to prioritize specific cell lines as model systems of cellular plasticity. We used two such models to demonstrate the dynamics, regulation and vulnerabilities associated with a cancer senescence program observed both in cell lines and in human tumors. Our work describes the landscape of cellular heterogeneity in diverse cancer cell lines, and identifies recurrent patterns of expression heterogeneity that are shared between tumors and specific cell lines and can thus be further explored in follow up studies.

Heiko Horn, Michael S Lawrence, Jessica Xin Hu, Elizabeth Worstell, Nina Ilic, Yashaswi Shrestha, Eejung Kim, Atanas Kamburov, Alireza Kashani, William C Hahn, Jesse S Boehm, Gad Getz, Kasper Lage

Heterogeneity across cancer makes it difficult to find driver genes with intermediate (2-20%) and low frequency (<2%) mutations, and we are potentially missing entire classes of networks (or pathways) of biological and therapeutic value. Here, we quantify the extent to which cancer genes across 21 tumor types have an increased burden of mutations in their immediate gene network derived from functional genomics data. We formalize a classifier that accurately calculates the significance level of a gene’s network mutation burden (NMB) and show it can accurately predict known cancer genes and recently proposed driver genes in the majority of tested tumours. Our approach predicts 62 putative cancer genes, including 35 with clear connection to cancer and 27 genes, which point to new cancer biology. NMB identifies proportionally more (4x) low-frequency mutated genes as putative cancer genes than gene-based tests, and provides molecular clues in patients without established driver mutations. Our quantitative and comparative analysis of pan-cancer networks across 21 tumour types gives new insights into the biological and genetic architecture of cancers and enables additional discovery from existing cancer genomes. The framework we present here should become increasingly useful with more sequencing data in the future.

Jeremiah Wala, Ofer Shapira, Marcin Imielinski, David Craft, Steven E Schumacher

Cancer cells can acquire profound alterations to the structure of their genomes, including rearrangements that fuse distant DNA breakpoints. We analyze the distribution of somatic rearrangements across the cancer genome, using whole-genome sequencing data from 2,693 tumor-normal pairs. We observe substantial variation in the density of rearrangement breakpoints, with enrichment in open chromatin and sites with high densities of repetitive elements. After accounting for these patterns, we identify significantly recurrent breakpoints (SRBs) at 52 loci, including novel SRBs near BRD4 and AKR1C3. Taking into account both loci fused by a rearrangement, we observe different signatures resembling either single breaks followed by strand invasion or two separate breaks that become joined. Accounting for these signatures, we identify 90 pairs of loci that are significantly recurrently juxtaposed (SRJs). SRJs are primarily tumor-type specific and tend to involve genes with tissue-specific expression. SRJs were frequently associated with disruption of topology-associated domains, juxtaposition of enhancer elements, and increased expression of neighboring genes. Lastly, we find that the power to detect SRJs decreases for short rearrangements, and that reliable detection of all driver SRJs will require whole-genome sequencing data from an order of magnitude more cancer samples than currently available.

Julien Racle, Kaat de Jonge, Petra Baumgaertner, Daniel E. Speiser, David Gfeller

Immune cells infiltrating tumors can have important impact on tumor progression and response to therapy. We present an efficient algorithm to simultaneously estimate the fraction of cancer and immune cell types from bulk tumor gene expression data. Our method integrates novel gene expression profiles from each major non-malignant cell type found in tumors, renormalization based on cell-type specific mRNA content, and the ability to consider uncharacterized and possibly highly variable cell types. Feasibility is demonstrated by validation with flow cytometry, immunohistochemistry and single-cell RNA-Seq analyses of human melanoma and colorectal tumor specimens. Altogether, our work not only improves accuracy but also broadens the scope of absolute cell fraction predictions from tumor gene expression data, and provides a unique novel experimental benchmark for immunogenomics analyses in cancer research.