SummaryLarge-scale genomic characterization of tumors from prospective cohort studies may yield new insights into cancer pathogenesis. We performed whole-exome sequencing of 619 incident colorectal cancers (CRCs) and integrated the results with tumor immunity, pathology, and survival data. We identified recurrently mutated genes in CRC, such as BCL9L, RBM10, CTCF, and KLF5, that were not previously appreciated in this disease. Furthermore, we investigated the genomic correlates of immune-cell infiltration and found that higher neoantigen load was positively associated with overall lymphocytic infiltration, tumor-infiltrating lymphocytes (TILs), memory T cells, and CRC-specific survival. The association with TILs was evident even within microsatellite-stable tumors. We also found positive selection of mutations in HLA genes and other components of the antigen-processing machinery in TIL-rich tumors. These results may inform immunotherapeutic approaches in CRC. More generally, this study demonstrates a framework for future integrative molecular epidemiology research in colorectal and other malignancies.

Table S3 has been updated. The updated version includes all somatic mutations used in our analyses.

Cisplatin-based chemotherapy is the standard of care for patients with muscle invasive urothelial carcinoma. Pathologic downstaging to pT0/pTis after neoadjuvant cisplatin-based chemotherapy is associated with improved survival, although molecular determinants of cisplatin response are incompletely understood. We performed whole exome sequencing on pre-treatment tumor and germline DNA from 50 patients with muscle invasive urothelial carcinoma who received neoadjuvant cisplatin-based chemotherapy followed by cystectomy (25 pT0/pTis “responders”, 25 pT2+ “non-responders”) to identify somatic mutations that occurred preferentially in responders. ERCC2, a nucleotide excision repair gene, was the only significantly mutated gene enriched in the cisplatin responders compared with non-responders (q < 0.01). Expression of representative ERCC2 mutations in an ERCC2-deficient cell line failed to rescue cisplatin and UV sensitivity compared to wild-type ERCC2. Lack of normal ERCC2 function may contribute to cisplatin sensitivity in urothelial cancer and somatic ERCC2 mutation status may inform cisplatin-containing regimen usage in muscle invasive urothelial carcinoma.

Patients with chromosome 13q deletion or normal cytogenetics represent the majority of chronic lymphocytic leukaemia (CLL) cases, yet have relatively few driver mutations. To better understand their genomic landscape, here we perform whole-genome sequencing on a cohort of patients enriched with these cytogenetic characteristics. Mutations in known CLL drivers are seen in only 33% of this cohort, and associated with normal cytogenetics and unmutated IGHV. The most commonly mutated gene in our cohort, IGLL5, shows a mutational pattern suggestive of activation-induced cytidine deaminase (AID) activity. Unsupervised analysis of mutational signatures demonstrates the activities of canonical AID (c-AID), leading to clustered mutations near active transcriptional start sites; non-canonical AID (nc-AID), leading to genome-wide non-clustered mutations, and an ageing signature responsible for most mutations. Using mutation clonality to infer time of onset, we find that while ageing and c-AID activities are ongoing, nc-AID-associated mutations likely occur earlier in tumour evolution.

BRAF mutations occur in ~10% of colorectal cancer (CRC). While RAF inhibitor monotherapy is highly effective in BRAF-mutant melanoma, response rates in BRAF-mutant CRC are poor. Recent clinical trials of combined RAF/EGFR or RAF/MEK inhibition have produced improved efficacy, but patients ultimately develop resistance. To identify molecular alterations driving clinical acquired resistance, we performed whole-exome sequencing on paired pre-treatment and post-progression tumor biopsies from BRAF-mutant CRC patients treated with RAF inhibitor combinations. We identified alterations in MAPK pathway genes in resistant tumors not present in matched pre-treatment tumors, including KRAS amplification, BRAF amplification, and a MEK1 mutation. These alterations conferred resistance to RAF/EGFR or RAF/MEK combinations through sustained MAPK pathway activity, but an ERK inhibitor could suppress MAPK activity and overcome resistance. Identification of MAPK pathway reactivating alterations upon clinical acquired resistance underscores the MAPK pathway as a critical target in BRAF-mutant CRC and suggests therapeutic options to overcome resistance.

African-American men have the highest incidence and mortality from prostate cancer. Whether a biological basis exists for this disparity remains unclear. Exome sequencing (n=102) and targeted validation (n = 90) of localized primary hormone-naïve prostate cancer in African-American men identified several gene mutations not previously observed in this context, including recurrent loss-of-function mutations in ERF, an ETS transcriptional repressor, in 5% of cases. Analysis of existing prostate cancer cohorts revealed ERF deletions in 3% of primary prostate cancers and mutations or deletions in ERF in 3–5% of lethal castration-resistant prostate cancers. Knockdown of ERF confers increased anchorage-independent growth and generates a gene expression signature associated with oncogenic ETS activation and androgen signaling. Together, these results suggest that ERF is a prostate cancer tumor suppressor gene. More generally, our findings support the application of systematic cancer genomic characterization in settings of broader ancestral diversity to enhance discovery and, eventually, therapeutic applications.

Analyses of frequency profiles of markers on disease or drug-response related genes in diverse populations are important for the dissection of common diseases. We report the results of analyses of data on 405 SNPs from 75 such genes and a 5.2 Mb chromosome, 22 genomic region in 1871 individuals from diverse 55 endogamous Indian populations. These include 32 large (>10 million individuals) and 23 isolated populations, representing a large fraction of the people of India. We observe high levels of genetic divergence between groups of populations that cluster largely on the basis of ethnicity and language. Indian populations not only overlap with the diversity of HapMap populations, but also contain population groups that are genetically distinct. These data and results are useful for addressing stratification and study design issues in complex traits especially for heterogeneous populations.

SummarySpinocerebellar ataxia type 12 (SCA12) is an autosomal dominant cerebellar ataxia associated with the expansion of an unstable CAG repeat in the 5 region of the PPP2R2B gene on chromosome 5q31-5q32. We found that it accounts for ∼ 16% (20/124) of all the autosomal dominant ataxia cases diagnosed in AIIMS, a major tertiary referral centre in North India. The length of the expanded allele in this population ranges from 51-69 CAG triplets. Interestingly, all the affected families belong to an endogamous population, which originated in the state of Haryana, India. We identified four novel SNPs and a dinucleotide marker spanning ∼ 137 kb downstream of CAG repeat in the PPP2R2B gene. Analysis of 20 Indian SCA12 families and ethnically matched normal unrelated individuals revealed one haplotype to be significantly associated with the affected alleles (P = 0.000), clearly indicating the presence of a common founder for SCA12 in the Indian population. This haplotype was not shared by the American pedigree with SCA12. Therefore, the SCA12 expansion appears to have originated at least twice.