Summary To characterize somatic alterations in colorectal carcinoma (CRC), we conducted genome-scale analysis of 276 samples, analyzing exome sequence, DNA copy number, promoter methylation, mRNA and microRNA expression. A subset (97) underwent low-depth-of-coverage whole-genome sequencing. 16% of CRC have hypermutation, three quarters of which have the expected high microsatellite instability (MSI), usually with hypermethylation and MLH1 silencing, but one quarter has somatic mismatch repair gene mutations. Excluding hypermutated cancers, colon and rectum cancers have remarkably similar patterns of genomic alteration. Twenty-four genes are significantly mutated. In addition to the expected APC, TP53, SMAD4, PIK3CA and KRAS mutations, we found frequent mutations in ARID1A, SOX9, and FAM123B/WTX. Recurrent copy number alterations include potentially drug-targetable amplifications of ERBB2 and newly discovered amplification of IGF2. Recurrent chromosomal translocations include fusion of NAV2 and WNT pathway member TCF7L1. Integrative analyses suggest new markers for aggressive CRC and important role for MYC-directed transcriptional activation and repression.

Major international projects are now underway aimed at creating a comprehensive catalog of all genes responsible for the initiation and progression of cancer. These studies involve sequencing of matched tumor–normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here, we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false positive findings that overshadow true driver events. Here, we show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumor-normal pairs and discover extraordinary variation in (i) mutation frequency and spectrum within cancer types, which shed light on mutational processes and disease etiology, and (ii) mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and allow true cancer genes to rise to attention.

Gastric cancer is a leading cause of cancer deaths, but analysis of its molecular and clinical characteristics has been complicated by histological and aetiological heterogeneity. Here we describe a comprehensive molecular evaluation of 295 primary gastric adenocarcinomas as part of The Cancer Genome Atlas (TCGA) project. We propose a molecular classification dividing gastric cancer into four subtypes: tumours positive for Epstein–Barr virus, which display recurrent PIK3CA mutations, extreme DNA hypermethylation, and amplification of JAK2, CD274 (also known as PD-L1) and PDCD1LG2 (also knownasPD-L2); microsatellite unstable tumours, which show elevated mutation rates, including mutations of genes encoding targetable oncogenic signalling proteins; genomically stable tumours, which are enriched for the diffuse histological variant and mutations of RHOA or fusions involving RHO-family GTPase-activating proteins; and tumours with chromosomal instability, which show marked aneuploidy and focal amplification of receptor tyrosine kinases. Identification of these subtypes provides a roadmap for patient stratification and trials of targeted therapies.

A powerful way to discover key genes playing causal roles in oncogenesis is to identify genomic regions that undergo frequent alteration in human cancers. Here, we report high-resolution analyses of somatic copy-number alterations (SCNAs) from 3131 cancer specimens, belonging largely to 26 histological types. We identify 158 regions of focal SCNA that are altered at significant frequency across multiple cancer types, of which 122 cannot be explained by the presence of a known cancer target gene located within these regions. Several gene families are enriched among these regions of focal SCNA, including the BCL2 family of apoptosis regulators and the NF-κB pathway. We show that cancer cells harboring amplifications surrounding the MCL1 and BCL2L1 anti-apoptotic genes depend upon expression of these genes for survival. Finally, we demonstrate that a large majority of SCNAs identified in individual cancer types are present in multiple cancer types.

The tumor microenvironment of colorectal carcinoma is a complex community of genomically altered cancer cells, nonneoplastic cells, and a diverse collection of microorganisms. Each of these components may contribute to carcinogenesis; however, the role of the microbiota is the least well understood. We have characterized the composition of the microbiota in colorectal carcinoma using whole genome sequences from nine tumor/normal pairs. Fusobacterium sequences were enriched in carcinomas, confirmed by quantitative PCR and 16S rDNA sequence analysis of 95 carcinoma/normal DNA pairs, while the Bacteroidetes and Firmicutes phyla were depleted in tumors. Fusobacteria were also visualized within colorectal tumors using FISH. These findings reveal alterations in the colorectal cancer microbiota; however, the precise role of Fusobacteria in colorectal carcinoma pathogenesis requires further investigation.

Most tumors display increased glucose metabolism compared to that of normal tissues. The preferential conversion of glucose to lactate in cancer cells (the Warburg Effect) has been emphasized1; however, the extent to which metabolic fluxes originating from glucose are utilized for alternative processes is poorly understood2,3. Here we used a combination of mass spectrometry and NMR with stable isotope labeling to investigate the alternate pathways derived from glucose metabolism in cancer cells. We found that in some cancer cells, a relatively large amount of glycolytic carbon is diverted into serine and glycine biosynthesis through phosphoglycerate dehydrogenase (PHGDH). A bioinformatics analysis of 3131 human cancers revealed that the gene PHGDH at 1p12 is recurrently amplified in a genomic region of focal copy number gain most commonly found in melanoma in which amplification was associated with increased protein expression. Decreased PHGDH expression by RNA interference impaired growth and flux into serine metabolism in PHGDH-amplified cell lines. Increased expression was also associated with breast cancer subtypes and ectopic expression of PHGDH in mammary epithelial cells (MCF-10a) disrupted acinar morphogenesis, induced loss of polarity, and preserved the viability of the extracellular matrix-deprived cells, each being phenotypic alterations that may predispose cells to transformation. Our findings demonstrate that altered metabolic flux from glucose into a specific alternate pathway can be selected during tumor development and may contribute to the pathogenesis of human cancer.

Clinical studies support the efficacy of programmed cell death 1 (PD-1) targeted therapy in a subset of patients with metastatic gastric cancer (mGC). With the goal of identifying determinants of response, we performed molecular characterization of tissues and circulating tumor DNA (ctDNA) from 61 patients with mGC who were treated with pembrolizumab as salvage treatment in a prospective phase 2 clinical trial. In patients with microsatellite instability-high and Epstein-Barr virus-positive tumors, which are mutually exclusive, dramatic responses to pembrolizumab were observed (overall response rate (ORR) 85.7% in microsatellite instability-high mGC and ORR 100% in Epstein-Barr virus-positive mGC). For the 55 patients for whom programmed death-ligand 1 (PD-L1) combined positive score positivity was available (combined positive score cut-off value ≥1%), ORR was significantly higher in PD-L1(+) gastric cancer when compared to PD-L1(-) tumors (50.0% versus 0.0%, P value <0.001). Changes in ctDNA levels at six weeks post-treatment predicted response and progression-free survival, and decreased ctDNA was associated with improved outcomes. Our findings provide insight into the molecular features associated with response to pembrolizumab in patients with mGC and provide biomarkers potentially relevant for the selection of patients who may derive greater benefit from PD-1 inhibition.

Lineage survival oncogenes are activated by somatic DNA alterations in cancers arising from the cell lineages in which these genes play a role in normal development.1,2 Here we show that a peak of genomic amplification on chromosome 3q26.33, found in squamous cell carcinomas (SCCs) of the lung and esophagus, contains the transcription factor gene SOX2—which is mutated in hereditary human esophageal malformations3 and necessary for normal esophageal squamous development4, promotes differentiation and proliferation of basal tracheal cells5 and co-operates in induction of pluripotent stem cells.6,7,8 SOX2 expression is required for proliferation and anchorage-independent growth of lung and esophageal cell lines, as shown by RNA interference experiments. Furthermore, ectopic expression of SOX2 cooperated with FOXE1 or FGFR2 to transform immortalized tracheobronchial epithelial cells. SOX2-driven tumors show expression of markers of both squamous differentiation and pluripotency. These observations identify SOX2 as a novel lineage survival oncogene in lung and esophageal SCC.