Genomic Evolution and Transcriptional Changes in the Evolution of Prostate Cancer into Neuroendocrine and Ductal Carcinoma Types

Prostate cancer is typically of acinar adenocarcinoma type but can occasionally present as neuroendocrine and/or ductal type carcinoma. These are associated with clinically aggressive disease, and the former often arises on a background of androgen deprivation therapy, although it can also arise de novo. Two prostate cancer cases were sequenced by exome capture from archival tissue. Case 1 was de novo small cell neuroendocrine carcinoma and ductal adenocarcinoma with three longitudinal samples over 5 years. Case 2 was a single time point after the development of treatment-related neuroendocrine prostate carcinoma. Case 1 showed whole genome doubling in all samples and focal amplification of AR in all samples except the first time point. Phylogenetic analysis revealed a common ancestry for ductal and small cell carcinoma. Case 2 showed 13q loss (involving RB1) in both adenocarcinoma and small cell carcinoma regions, and 3p gain, 4p loss, and 17p loss (involving TP53) in the latter. By using highly curated samples, we demonstrate for the first time that small-cell neuroendocrine and ductal prostatic carcinoma can have a common ancestry. We highlight whole genome doubling in a patient with prostate cancer relapse, reinforcing its poor prognostic nature.

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Phenome-wide association study (PheWAS) of colorectal cancer risk SNP effects on health outcomes in UK Biobank

BACKGROUND

Associations between colorectal cancer (CRC) and other health outcomes have been reported, but these may be subject to biases, or due to limitations of observational studies.

METHODS

We set out to determine whether genetic predisposition to CRC is also associated with the risk of other phenotypes. Under the phenome-wide association study (PheWAS) and tree-structured phenotypic model (TreeWAS), we studied 334,385 unrelated White British individuals (excluding CRC patients) from the UK Biobank cohort. We generated a polygenic risk score (PRS) from CRC genome-wide association studies as a measure of CRC risk. We performed sensitivity analyses to test the robustness of the results and searched the Danish Disease Trajectory Browser (DTB) to replicate the observed associations.

RESULTS

Eight PheWAS phenotypes and 21 TreeWAS nodes were associated with CRC genetic predisposition by PheWAS and TreeWAS, respectively. The PheWAS detected associations were from neoplasms and digestive system disease group (e.g. benign neoplasm of colon, anal and rectal polyp and diverticular disease). The results from the TreeWAS corroborated the results from the PheWAS. These results were replicated in the observational data within the DTB.

CONCLUSIONS

We show that benign colorectal neoplasms share genetic aetiology with CRC using PheWAS and TreeWAS methods. Additionally, CRC genetic predisposition is associated with diverticular disease.

Genetically predicted physical activity levels are associated with lower colorectal cancer risk: a Mendelian randomisation study

BACKGROUND

We conducted a Mendelian randomisation (MR) study to investigate whether physical activity (PA) causes a reduction of colorectal cancer risk and to understand the contributions of effects mediated through changes in body fat.

METHODS

Common genetic variants associated with self-reported moderate-to-vigorous PA (MVPA), acceleration vector magnitude PA (AMPA) and sedentary time were used as instrumental variables. To control for confounding effects of obesity, we included instrumental variables for body mass index (BMI), body fat percentage, waist circumference and arm, trunk and leg fat ratios. We analysed the effect of these instrumental variables in a colorectal cancer genome-wide association study comprising 31,197 cases and 61,770 controls of European ancestry by applying two-sample and multivariable MR study designs.

RESULTS

We found decreased colorectal cancer risk for genetically represented measures of MVPA and AMPA that were additional to effects mediated through genetic measures of obesity. Odds ratio and 95% confidence interval (CI) per standard deviation increase in MVPA and AMPA was 0.56 (0.31, 1.01) and 0.60 (0.41, 0.88), respectively. No association has been found between sedentary time and colorectal cancer risk. The proportion of effect mediated through BMI was 2% (95% CI: 0, 14) and 32% (95% CI: 12, 46) for MVPA and AMPA, respectively.

CONCLUSION

These findings provide strong evidence to reinforce public health measures on preventing colorectal cancer that promote PA at a population level regardless of body fatness.

Prediction of colorectal cancer risk based on profiling with common genetic variants

Increasing numbers of common genetic variants associated with colorectal cancer (CRC) have been identified. Our study aimed to determine whether risk prediction based on common genetic variants might enable stratification for CRC risk. Meta-analysis of 11 genome-wide association studies comprising 16 871 cases and 26 328 controls was performed to capture CRC susceptibility variants. Genetic prediction models with several candidate polygenic risk scores (PRSs) were generated from Scottish CRC case-control studies (6478 cases and 11 043 controls) and the score with the best performance was then tested in UK Biobank (UKBB) (4800 cases and 20 287 controls). A weighted PRS of 116 CRC single nucleotide polymorphisms (wPRS116 ) was found with the best predictive performance, reporting a c-statistics of 0.60 and an odds ratio (OR) of 1.46 (95% confidence interval [CI] = 1.41-1.50, per SD increase) in Scottish data set. The predictive performance of this wPRS116 was consistently validated in UKBB data set with c-statistics of 0.61 and an OR of 1.49 (95% CI = 1.44-1.54, per SD increase). Modeling the levels of PRS with age and sex in the general UK population shows that employing genetic risk profiling can achieve a moderate degree of risk discrimination that could be helpful to identify a subpopulation with higher CRC risk due to genetic susceptibility.

A sensitive and affordable multiplex RT-qPCR assay for SARS-CoV-2 detection

With the ongoing COVID-19 (Coronavirus Disease 2019) pandemic, caused by the novel coronavirus SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2), there is a need for sensitive, specific, and affordable diagnostic tests to identify infected individuals, not all of whom are symptomatic. The most sensitive test involves the detection of viral RNA using RT-qPCR (quantitative reverse transcription PCR), with many commercial kits now available for this purpose. However, these are expensive, and supply of such kits in sufficient numbers cannot always be guaranteed. We therefore developed a multiplex assay using well-established SARS-CoV-2 targets alongside a human cellular control (RPP30) and a viral spike-in control (Phocine Herpes Virus 1 [PhHV-1]), which monitor sample quality and nucleic acid extraction efficiency, respectively. Here, we establish that this test performs as well as widely used commercial assays, but at substantially reduced cost. Furthermore, we demonstrate >1,000-fold variability in material routinely collected by combined nose and throat swabbing and establish a statistically significant correlation between the detected level of human and SARS-CoV-2 nucleic acids. The inclusion of the human control probe in our assay therefore provides a quantitative measure of sample quality that could help reduce false-negative rates. We demonstrate the feasibility of establishing a robust RT-qPCR assay at approximately 10% of the cost of equivalent commercial assays, which could benefit low-resource environments and make high-volume testing affordable.

Pervasive chromosomal instability and karyotype order in tumour evolution

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Why is cancer not more common? A changing microenvironment may help to explain why, and suggests strategies for anti-cancer therapy

One of the great unsolved puzzles in cancer biology is not why cancers occur, but rather explaining why so few cancers occur compared with the theoretical number that could occur, given the number of progenitor cells in the body and the normal mutation rate. We hypothesized that a contributory explanation is that the tumour microenvironment (TME) is not fixed due to factors such as immune cell infiltration, and that this could impair the ability of neoplastic cells to retain a high enough fitness to become a cancer. The TME has implicitly been assumed to be static in most cancer evolution models, and we therefore developed a mathematical model of spatial cancer evolution assuming that the TME, and thus the optimum cancer phenotype, changes over time. Based on simulations, we show how cancer cell populations adapt to diverse changing TME conditions and fitness landscapes. Compared with static TMEs, which generate neutral dynamics, changing TMEs lead to complex adaptations with characteristic spatio-temporal heterogeneity involving variable fitness effects of driver mutations, subclonal mixing, subclonal competition and phylogeny patterns. In many cases, cancer cell populations fail to grow or undergo spontaneous regression, and even extinction. Our analyses predict that cancer evolution in a changing TME is challenging, and can help to explain why cancer is neither inevitable nor as common as expected. Should cancer driver mutations with effects dependent of the TME exist, they are likely to be selected. Anti-cancer prevention and treatment strategies based on changing the TME are feasible and potentially effective.

Measuring single cell divisions in human tissues from multi-region sequencing data

Both normal tissue development and cancer growth are driven by a branching process of cell division and mutation accumulation that leads to intra-tissue genetic heterogeneity. However, quantifying somatic evolution in humans remains challenging. Here, we show that multi-sample genomic data from a single time point of normal and cancer tissues contains information on single-cell divisions. We present a new theoretical framework that, applied to whole-genome sequencing data of healthy tissue and cancer, allows inferring the mutation rate and the cell survival/death rate per division. On average, we found that cells accumulate 1.14 mutations per cell division in healthy haematopoiesis and 1.37 mutations per division in brain development. In both tissues, cell survival was maximal during early development. Analysis of 131 biopsies from 16 tumours showed 4 to 100 times increased mutation rates compared to healthy development and substantial inter-patient variation of cell survival/death rates.

Mendelian randomisation: A powerful and inexpensive method for identifying and excluding non-genetic risk factors for colorectal cancer

Colorectal cancer (CRC) is the third most common cancer in economically developed countries and a major cause of cancer-related mortality. The importance of lifestyle and diet as major determinants of CRC risk is suggested by differences in CRC incidence between countries and in migration studies. Previous observational epidemiological studies have identified associations between modifiable environmental risk factors and CRC, but these studies can be susceptible to reverse causation and confounding, and their results can therefore conflict. Mendelian randomisation (MR) analysis represents an approach complementary to conventional observational studies examining associations between exposures and disease. The MR strategy employs allelic variants as instrumental variables (IVs), which act as proxies for non-genetic exposures. These allelic variants are randomly assigned during meiosis and can therefore inform on life-long exposure, whilst not being subject to reverse causation. In previous studies MR frameworks have associated several modifiable factors with CRC risk, including adiposity, hyperlipidaemia, fatty acid profile and alcohol consumption. In this review we detail the use of MR to investigate and discover CRC risk factors, and its future applications.

Evolutionary history of human colitis-associated colorectal cancer

OBJECTIVE

IBD confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multiregion sequencing.

DESIGN

Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 patients with CA-CRC (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using single nucleotide polymorphism arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD; n=30), high-grade dysplasia (HGD; n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC.

RESULTS

10/12 tumours were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNA) per Mb, ~20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated ageing-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighbouring CA-CRC, indicating a precancer mutational field. CA-CRCs were often near triploid (40%) or near tetraploid (20%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated 'catastrophic' CNA increase.

CONCLUSIONS

Evolutionary genomic analysis revealed precancer clones bearing extensive SNAs and CNAs, with progression to cancer involving a dramatic accrual of CNAs at HGD. Detection of the cancerised field is an encouraging prospect for surveillance, but punctuated evolution may limit the window for early detection.

The evolutionary landscape of colorectal tumorigenesis

The evolutionary events that cause colorectal adenomas (benign) to progress to carcinomas (malignant) remain largely undetermined. Using multi-region genome and exome sequencing of 24 benign and malignant colorectal tumours, we investigate the evolutionary fitness landscape occupied by these neoplasms. Unlike carcinomas, advanced adenomas frequently harbour sub-clonal driver mutations-considered to be functionally important in the carcinogenic process-that have not swept to fixation, and have relatively high genetic heterogeneity. Carcinomas are distinguished from adenomas by widespread aneusomies that are usually clonal and often accrue in a 'punctuated' fashion. We conclude that adenomas evolve across an undulating fitness landscape, whereas carcinomas occupy a sharper fitness peak, probably owing to stabilizing selection.

Adjuvant Treatment for POLE Proofreading Domain-Mutant Cancers: Sensitivity to Radiotherapy, Chemotherapy, and Nucleoside Analogues

Purpose: Pathogenic POLE proofreading domain mutations are found in many malignancies where they are associated with ultramutation and favorable prognosis. The extent to which this prognosis depends on their sensitivity to adjuvant treatment is unknown, as is the optimal therapy for advanced-staged or recurrent POLE-mutant cancers.Experimental Design: We examined the recurrence-free survival of women with POLE-mutant and POLE-wild-type endometrial cancers (EC) in the observation arm of the randomized PORTEC-1 endometrial cancer trial (N = 245 patients with stage I endometrial cancer for analysis). Sensitivity to radiotherapy and selected chemotherapeutics was compared between Pole-mutant mouse-derived embryonic stem (mES) cells, generated using CRISPR-Cas9 (Pole mutations D275A/E275A, and cancer-associated P286R, S297F, V411L) and isogenic wild-type cell lines.Results: In the observation arm of the PORTEC-1 trial (N = 245), women with POLE-mutant endometrial cancers (N = 16) had an improved recurrence-free survival (10-year recurrence-free survival 100% vs. 80.1% for POLE-wild-type; HR, 0.143; 95% confidence interval, 0.001-0.996; P = 0.049). Pole mutations did not increase sensitivity to radiotherapy nor to chemotherapeutics in mES cells. In contrast, Pole-mutant cells displayed significantly increased sensitivity to cytarabine and fludarabine (IC50Pole P286R-mutant vs. wild-type: 0.05 vs. 0.17 μmol/L for cytarabine, 4.62 vs. 11.1 μmol/L for fludarabine; P < 0.001 for both comparisons).Conclusions: The favorable prognosis of POLE-mutant cancers cannot be explained by increased sensitivity to currently used adjuvant treatments. These results support studies exploring minimization of adjuvant therapy for early-stage POLE-mutant cancers, including endometrial and colorectal cancers. Conversely, POLE mutations result in hypersensitivity to nucleoside analogues, suggesting the use of these compounds as a potentially effective targeted treatment for advanced-stage POLE-mutant cancers. Clin Cancer Res; 24(13); 3197-203. ©2018 AACR.

The effects of mutational processes and selection on driver mutations across cancer types

Epidemiological evidence has long associated environmental mutagens with increased cancer risk. However, links between specific mutation-causing processes and the acquisition of individual driver mutations have remained obscure. Here we have used public cancer sequencing data from 11,336 cancers of various types to infer the independent effects of mutation and selection on the set of driver mutations in a cancer type. First, we detect associations between a range of mutational processes, including those linked to smoking, ageing, APOBEC and DNA mismatch repair (MMR) and the presence of key driver mutations across cancer types. Second, we quantify differential selection between well-known alternative driver mutations, including differences in selection between distinct mutant residues in the same gene. These results show that while mutational processes have a large role in determining which driver mutations are present in a cancer, the role of selection frequently dominates.

Recurrent Coding Sequence Variation Explains Only A Small Fraction of the Genetic Architecture of Colorectal Cancer

Whilst common genetic variation in many non-coding genomic regulatory regions are known to impart risk of colorectal cancer (CRC), much of the heritability of CRC remains unexplained. To examine the role of recurrent coding sequence variation in CRC aetiology, we genotyped 12,638 CRCs cases and 29,045 controls from six European populations. Single-variant analysis identified a coding variant (rs3184504) in SH2B3 (12q24) associated with CRC risk (OR = 1.08, P = 3.9 × 10(-7)), and novel damaging coding variants in 3 genes previously tagged by GWAS efforts; rs16888728 (8q24) in UTP23 (OR = 1.15, P = 1.4 × 10(-7)); rs6580742 and rs12303082 (12q13) in FAM186A (OR = 1.11, P = 1.2 × 10(-7) and OR = 1.09, P = 7.4 × 10(-8)); rs1129406 (12q13) in ATF1 (OR = 1.11, P = 8.3 × 10(-9)), all reaching exome-wide significance levels. Gene based tests identified associations between CRC and PCDHGA genes (P < 2.90 × 10(-6)). We found an excess of rare, damaging variants in base-excision (P = 2.4 × 10(-4)) and DNA mismatch repair genes (P = 6.1 × 10(-4)) consistent with a recessive mode of inheritance. This study comprehensively explores the contribution of coding sequence variation to CRC risk, identifying associations with coding variation in 4 genes and PCDHG gene cluster and several candidate recessive alleles. However, these findings suggest that recurrent, low-frequency coding variants account for a minority of the unexplained heritability of CRC.

Fine-scale mapping of the 4q24 locus identifies two independent loci associated with breast cancer risk

BACKGROUND

A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored.

METHODS

We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium.

RESULTS

Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10(-4); OR, 1.04; 95% confidence interval (CI), 1.02-1.07] and rs77928427 (P = 1.86 × 10(-4); OR, 1.04; 95% CI, 1.02-1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r(2) ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor-binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue.

CONCLUSION

Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2.

IMPACT

Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk.

Methylation changes in the TFAP2E promoter region are associated with BRAF mutation and poorer overall & disease free survival in colorectal cancer

INTRODUCTION

BRAF mutant colorectal cancer carries a poor prognosis which is thought to be related to poor response to conventional chemotherapy. BRAF mutation is associated with the serrated tumour phenotype. We hypothesised that one of the mechanisms by which BRAF mutant colorectal cancer demonstrate poor outcomes with chemotherapy is abnormal gene methylation.

METHODS

The Cancer Genome Atlas (TCGA) methylation data was analysed using a linear regression model with BRAF mutation as an independent variable. Expression datasets were also obtained to correlate functional changes. Top differentially methylated probes were taken forward for validation by methylation pyrosequencing. These probes were analysed on a cohort of patients enriched for BRAF mutations taken from the VICTOR and QUASAR2 studies.

RESULTS

In an analysis of 91 tumours (9 BRAF mutant, 82 wild type), the Illumina probe cg11835197 was the probe identified as top differentially methylated (p = 2.56×10-7, Bayes Factor (BF) =6.54). This probe covered a region -413bp from the promoter region of TFAP2E. We found a complex pattern of CpG specific methylation of this region which was associated with both overall (p=0.044) and disease free (p=0.046) survival.

DISCUSSION

BRAF mutant tumours may attain part of their chemoresistance from abnormal TFAP2E methylation, which has not previously been described.

Prognostic significance of POLE proofreading mutations in endometrial cancer

BACKGROUND

Current risk stratification in endometrial cancer (EC) results in frequent over- and underuse of adjuvant therapy, and may be improved by novel biomarkers. We examined whether POLE proofreading mutations, recently reported in about 7% of ECs, predict prognosis.

METHODS

We performed targeted POLE sequencing in ECs from the PORTEC-1 and -2 trials (n = 788), and analyzed clinical outcome according to POLE status. We combined these results with those from three additional series (n = 628) by meta-analysis to generate multivariable-adjusted, pooled hazard ratios (HRs) for recurrence-free survival (RFS) and cancer-specific survival (CSS) of POLE-mutant ECs. All statistical tests were two-sided.

RESULTS

POLE mutations were detected in 48 of 788 (6.1%) ECs from PORTEC-1 and-2 and were associated with high tumor grade (P < .001). Women with POLE-mutant ECs had fewer recurrences (6.2% vs 14.1%) and EC deaths (2.3% vs 9.7%), though, in the total PORTEC cohort, differences in RFS and CSS were not statistically significant (multivariable-adjusted HR = 0.43, 95% CI = 0.13 to 1.37, P = .15; HR = 0.19, 95% CI = 0.03 to 1.44, P = .11 respectively). However, of 109 grade 3 tumors, 0 of 15 POLE-mutant ECs recurred, compared with 29 of 94 (30.9%) POLE wild-type cancers; reflected in statistically significantly greater RFS (multivariable-adjusted HR = 0.11, 95% CI = 0.001 to 0.84, P = .03). In the additional series, there were no EC-related events in any of 33 POLE-mutant ECs, resulting in a multivariable-adjusted, pooled HR of 0.33 for RFS (95% CI = 0.12 to 0.91, P = .03) and 0.26 for CSS (95% CI = 0.06 to 1.08, P = .06).

CONCLUSION

POLE proofreading mutations predict favorable EC prognosis, independently of other clinicopathological variables, with the greatest effect seen in high-grade tumors. This novel biomarker may help to reduce overtreatment in EC.

Mutations in Exon 1 highlight the role of MED12 in uterine leiomyomas

Mediator regulates transcription by connecting gene-specific transcription factors to the RNA polymerase II initiation complex. We recently discovered by exome sequencing that specific exon 2 mutations in mediator complex subunit 12 (MED12) are extremely common in uterine leiomyomas. Subsequent screening studies have focused on this mutational hot spot, and mutations have been detected in uterine leiomyosarcomas, extrauterine leiomyomas and leiomyosarcomas, endometrial polyps, and colorectal cancers. All mutations have been missense changes or in-frame insertions/deletions. Here, we have analyzed 611 samples representing all above-mentioned tumor types for possible exon 1 mutations. Five mutations were observed, all of which were in-frame insertion/deletions in uterine leiomyomas. Transcriptome-wide expression data revealed that MED12 exon 1 and exon 2 mutations lead to the same unique global gene expression pattern with RAD51B being the most upregulated gene. Immunoprecipitation and kinase activity assays showed that both exon 1 and exon 2 mutations disrupt the interaction between MED12 and Cyclin C and CDK8/19 and abolish the mediator-associated CDK kinase activity. These results further emphasize the role of MED12 in uterine leiomyomas, show that exon 1 and exon 2 exert their tumorigenic effect in similar manner, and stress that exon 1 should be included in subsequent MED12 screenings.

Colorectal cancer cell lines are representative models of the main molecular subtypes of primary cancer

Human colorectal cancer cell lines are used widely to investigate tumor biology, experimental therapy, and biomarkers. However, to what extent these established cell lines represent and maintain the genetic diversity of primary cancers is uncertain. In this study, we profiled 70 colorectal cancer cell lines for mutations and DNA copy number by whole-exome sequencing and SNP microarray analyses, respectively. Gene expression was defined using RNA-Seq. Cell line data were compared with those published for primary colorectal cancers in The Cancer Genome Atlas. Notably, we found that exome mutation and DNA copy-number spectra in colorectal cancer cell lines closely resembled those seen in primary colorectal tumors. Similarities included the presence of two hypermutation phenotypes, as defined by signatures for defective DNA mismatch repair and DNA polymerase ε proofreading deficiency, along with concordant mutation profiles in the broadly altered WNT, MAPK, PI3K, TGFβ, and p53 pathways. Furthermore, we documented mutations enriched in genes involved in chromatin remodeling (ARID1A, CHD6, and SRCAP) and histone methylation or acetylation (ASH1L, EP300, EP400, MLL2, MLL3, PRDM2, and TRRAP). Chromosomal instability was prevalent in nonhypermutated cases, with similar patterns of chromosomal gains and losses. Although paired cell lines derived from the same tumor exhibited considerable mutation and DNA copy-number differences, in silico simulations suggest that these differences mainly reflected a preexisting heterogeneity in the tumor cells. In conclusion, our results establish that human colorectal cancer lines are representative of the main subtypes of primary tumors at the genomic level, further validating their utility as tools to investigate colorectal cancer biology and drug responses.

Evaluation of PIK3CA mutation as a predictor of benefit from nonsteroidal anti-inflammatory drug therapy in colorectal cancer

PURPOSE

Aspirin and other nonsteroidal anti-inflammatory drugs (NSAIDs) protect against colorectal cancer (CRC) and are associated with reduced disease recurrence and improved outcome after primary treatment. However, toxicities of NSAIDs have limited their use as antineoplastic therapy. Recent data have suggested that the benefit of aspirin after CRC diagnosis is limited to patients with PIK3CA-mutant cancers. We sought to determine the predictive utility of PIK3CA mutation for benefit from both cyclooxygenase-2 inhibition and aspirin.

METHODS

We performed molecular analysis of tumors from 896 participants in the Vioxx in Colorectal Cancer Therapy: Definition of Optimal Regime (VICTOR) trial, a large randomized trial comparing rofecoxib with placebo after primary CRC resection. We compared relapse-free survival and overall survival between rofecoxib therapy and placebo and between the use and nonuse of low-dose aspirin, according to tumor PIK3CA mutation status.

RESULTS

We found no evidence of a greater benefit from rofecoxib treatment compared with placebo in patients whose tumors had PIK3CA mutations (multivariate adjusted hazard ratio [HR], 1.2; 95% CI, 0.53 to 2.72; P = .66; (P)INTERACTION = .47) compared with patients with PIK3CA wild-type cancers (HR, 0.87; 95% CI, 0.64 to 1.16; P = .34). In contrast, regular aspirin use after CRC diagnosis was associated with a reduced rate of CRC recurrence in patients with PIK3CA-mutant cancers (HR, 0.11; 95% CI, 0.001 to 0.832; P = .027; (P)INTERACTION = .024) but not in patients lacking tumor PIK3CA mutation (HR, 0.92; 95% CI, 0.60 to 1.42; P = .71).

CONCLUSION

Although tumor PIK3CA mutation does not predict benefit from rofecoxib treatment, it merits further evaluation as a predictive biomarker for aspirin therapy. Our findings are concordant with recent data and support the prospective investigation of adjuvant aspirin in PIK3CA-mutant CRC.

A study of genomic instability in early preneoplastic colonic lesions

It is difficult to explain the differential rates of progression of premalignant colonic lesions and differences in behaviour of morphologically similar lesions. Heterogeneity for microsatellite instability (MSI) and promoter methylation in driving these phenomena forward may explain this; however, no previous analysis has examined this in detail at the gland level, the smallest unit of colorectal premalignant lesions. We aimed to carry out an analysis of gland level genomic instability for MSI and promoter methylation. MSI occurred significantly more frequently (20%) in colonic glands than has previously been observed in whole colorectal polyps. Significant promoter methylation was seen in MLH1, PMS2, MLH3 and MSH3 as well as significant heterogeneity for both MSI and promoter methylation. Methylation and MSI may have a significant role in driving forward colorectal carcinogenesis, although in the case of MSI, this association is less clear as it occurs significantly more frequently than previously thought, and may simply be a passenger in the adenoma-carcinoma sequence. Promoter methylation in MLH1, MLH3, MSH3 and PMS2 was also found to be significantly associated with MSI and should be investigated further. A total of 273 colorectal glands (126 hyperplastic, 147 adenomatous) were isolated via laser capture microdissection (targeted at regions of MLH1 loss) from 93 colonic polyps and tested for MSI, and promoter methylation of the DNA mismatch repair genes MLH1, MSH2, MLH3, MSH6, PMS2, MGMT and MLH3 via methylation specific multiplex ligation-dependent probe amplification. Logistic regression modelling was then used to identify significant associations between promoter methylation and gland histological type and MSI status.

Survival in stage II/III colorectal cancer is independently predicted by chromosomal and microsatellite instability, but not by specific driver mutations

OBJECTIVES

Microsatellite instability (MSI) is an established marker of good prognosis in colorectal cancer (CRC). Chromosomal instability (CIN) is strongly negatively associated with MSI and has been shown to be a marker of poor prognosis in a small number of studies. However, a substantial group of "double-negative" (MSI-/CIN-) CRCs exists. The prognosis of these patients is unclear. Furthermore, MSI and CIN are each associated with specific molecular changes, such as mutations in KRAS and BRAF, that have been associated with prognosis. It is not known which of MSI, CIN, and the specific gene mutations are primary predictors of survival.

METHODS

We evaluated the prognostic value (disease-free survival, DFS) of CIN, MSI, mutations in KRAS, NRAS, BRAF, PIK3CA, FBXW7, and TP53, and chromosome 18q loss-of-heterozygosity (LOH) in 822 patients from the VICTOR trial of stage II/III CRC. We followed up promising associations in an Australian community-based cohort (N=375).

RESULTS

In the VICTOR patients, no specific mutation was associated with DFS, but individually MSI and CIN showed significant associations after adjusting for stage, age, gender, tumor location, and therapy. A combined analysis of the VICTOR and community-based cohorts showed that MSI and CIN were independent predictors of DFS (for MSI, hazard ratio (HR)=0.58, 95% confidence interval (CI) 0.36-0.93, and P=0.021; for CIN, HR=1.54, 95% CI 1.14-2.08, and P=0.005), and joint CIN/MSI testing significantly improved the prognostic prediction of MSI alone (P=0.028). Higher levels of CIN were monotonically associated with progressively poorer DFS, and a semi-quantitative measure of CIN was a better predictor of outcome than a simple CIN+/- variable. All measures of CIN predicted DFS better than the recently described Watanabe LOH ratio.

CONCLUSIONS

MSI and CIN are independent predictors of DFS for stage II/III CRC. Prognostic molecular tests for CRC relapse should currently use MSI and a quantitative measure of CIN rather than specific gene mutations.

Different APC genotypes in proximal and distal sporadic colorectal cancers suggest distinct WNT/β-catenin signalling thresholds for tumourigenesis

Biallelic protein-truncating mutations in the adenomatous polyposis coli (APC) gene are prevalent in sporadic colorectal cancer (CRC). Mutations may not be fully inactivating, instead producing WNT/β-catenin signalling levels 'just-right' for tumourigenesis. However, the spectrum of optimal APC genotypes accounting for both hits, and the influence of clinicopathological features on genotype selection remain undefined. We analysed 630 sporadic CRCs for APC mutations and loss of heterozygosity (LOH) using sequencing and single-nucleotide polymorphism microarrays, respectively. Truncating APC mutations and/or LOH were detected in 75% of CRCs. Most truncating mutations occurred within a mutation cluster region (MCR; codons 1282-1581) leaving 1-3 intact 20 amino-acid repeats (20AARs) and abolishing all Ser-Ala-Met-Pro (SAMP) repeats. Cancers commonly had one MCR mutation plus either LOH or another mutation 5' to the MCR. LOH was associated with mutations leaving 1 intact 20AAR. MCR mutations leaving 1 vs 2-3 intact 20AARs were associated with 5' mutations disrupting or leaving intact the armadillo-repeat domain, respectively. Cancers with three hits had an over-representation of mutations upstream of codon 184, in the alternatively spliced region of exon 9, and 3' to the MCR. Microsatellite unstable cancers showed hyper-mutation at MCR mono- and di-nucleotide repeats, leaving 2-3 intact 20AARs. Proximal and distal cancers exhibited different preferred APC genotypes, leaving a total of 2 or 3 and 0 to 2 intact 20AARs, respectively. In conclusion, APC genotypes in sporadic CRCs demonstrate 'fine-tuned' interdependence of hits by type and location, consistent with selection for particular residual levels of WNT/β-catenin signalling, with different 'optimal' thresholds for proximal and distal cancers.

DNA polymerase ε and δ exonuclease domain mutations in endometrial cancer

Accurate duplication of DNA prior to cell division is essential to suppress mutagenesis and tumour development. The high fidelity of eukaryotic DNA replication is due to a combination of accurate incorporation of nucleotides into the nascent DNA strand by DNA polymerases, the recognition and removal of mispaired nucleotides (proofreading) by the exonuclease activity of DNA polymerases δ and ε, and post-replication surveillance and repair of newly synthesized DNA by the mismatch repair (MMR) apparatus. While the contribution of defective MMR to neoplasia is well recognized, evidence that faulty DNA polymerase activity is important in cancer development has been limited. We have recently shown that germline POLE and POLD1 exonuclease domain mutations (EDMs) predispose to colorectal cancer (CRC) and, in the latter case, to endometrial cancer (EC). Somatic POLE mutations also occur in 5-10% of sporadic CRCs and underlie a hypermutator, microsatellite-stable molecular phenotype. We hypothesized that sporadic ECs might also acquire somatic POLE and/or POLD1 mutations. Here, we have found that missense POLE EDMs with good evidence of pathogenic effects are present in 7% of a set of 173 endometrial cancers, although POLD1 EDMs are uncommon. The POLE mutations localized to highly conserved residues and were strongly predicted to affect proofreading. Consistent with this, POLE-mutant tumours were hypermutated, with a high frequency of base substitutions, and an especially large relative excess of G:C>T:A transversions. All POLE EDM tumours were microsatellite stable, suggesting that defects in either DNA proofreading or MMR provide alternative mechanisms to achieve genomic instability and tumourigenesis.

What can be learnt about disease progression in breast cancer dormancy from relapse data?

Breast cancer patients have an anomalously high rate of relapse many years–up to 25 years–after apparently curative surgery removed the primary tumour. Disease progression during the intervening years between resection and relapse is poorly understood. There is evidence that the disease persists as dangerous, tiny metastases that remain at a growth restricted, clinically undetectable size until a transforming event restarts growth. This is the starting point for our study, where patients who have metastases that are all tiny and growth-restricted are said to have cancer dormancy. Can long-term follow-up relapse data from breast cancer patients be used to extract knowledge about the progression of the undetected disease? Here, we evaluate whether this is the case by introducing and analysing four simple mathematical models of cancer dormancy. These models extend the common assumption that a random transforming event, such as a mutation, can restart growth of a tiny, growth-restricted metastasis; thereafter, cancer dormancy progresses to detectable metastasis. We find that physiopathological details, such as the number of random transforming events that metastases must undergo to escape from growth restriction, cannot be extracted from relapse data. This result is unsurprising. However, the same analysis suggested a natural question that does have a surprising answer: why are interesting trends in long-term relapse data not more commonly observed? Further, our models indicate that (a) therapies which induce growth restriction among metastases but do not prevent increases in metastases' tumourigenicity may introduce a time post-surgery when more patients are prone to relapse; and (b), if a number of facts about disease progression are first established, how relapse data might be used to estimate clinically relevant variables, such as the likely numbers of undetected growth-restricted metastases. This work is a necessary, early step in building a quantitative mechanistic understanding of cancer dormancy.

Loss of expression and promoter methylation of SLIT2 are associated with sessile serrated adenoma formation

Serrated adenomas form a distinct subtype of colorectal pre-malignant lesions that may progress to malignancy along a different molecular pathway than the conventional adenoma-carcinoma pathway. Previous studies have hypothesised that BRAF mutation and promoter hypermethylation plays a role, but the evidence for this is not robust. We aimed to carry out a whole-genome loss of heterozygosity analysis, followed by targeted promoter methylation and expression analysis to identify potential pathways in serrated adenomas. An initial panel of 9 sessile serrated adenomas (SSA) and one TSA were analysed using Illumina Goldengate HumanLinkage panel arrays to ascertain regions of loss of heterozygosity. This was verified via molecular inversion probe analysis and microsatellite analysis of a further 32 samples. Methylation analysis of genes of interest was carried out using methylation specific PCR (verified by pyrosequencing) and immunohistochemistry used to correlate loss of expression of genes of interest. All experiments used adenoma samples and normal tissue samples as control. SSA samples were found on whole-genome analysis to have consistent loss of heterozygosity at 4p15.1-4p15.31, which was not found in the sole TSA, adenomas, or normal tissues. Genes of interest in this region were PDCH7 and SLIT2, and combined MSP/IHC analysis of these genes revealed significant loss of SLIT2 expression associated with promoter methylation of SLIT2. Loss of expression of SLIT2 by promoter hypermethylation and loss of heterozygosity events is significantly associated with serrated adenoma development, and SLIT2 may represent a epimutated tumour suppressor gene according to the Knudson "two hit" hypothesis.

Whole-genome methylation analysis of benign and malignant colorectal tumours

Changes in DNA methylation, whether hypo- or hypermethylation, have been shown to be associated with the progression of colorectal cancer. Methylation changes substantially in the progression from normal mucosa to adenoma and to carcinoma. This phenomenon has not been studied extensively and studies have been restricted to individual CpG islands, rather than taking a whole-genome approach. We aimed to study genome-wide methylation changes in colorectal cancer. We obtained 10 fresh-frozen normal tissue-cancer sample pairs, and five fresh-frozen adenoma samples. These were run on the lllumina HumanMethylation27 whole-genome methylation analysis system. Differential methylation between normal tissue, adenoma and carcinoma was analysed using Bayesian regression modelling, gene set enrichment analysis (GSEA) and hierarchical clustering (HC). The highest-rated individual gene for differential methylation in carcinomas versus normal tissue and adenomas versus normal tissue was GRASP (padjusted  = 1.59 × 10(-5) , BF = 12.62, padjusted  = 1.68 × 10(-6) , BF = 14.53). The highest-rated gene when comparing carcinomas versus adenomas was ATM (padjusted  = 2.0 × 10(-4) , BF = 10.17). Hierarchical clustering demonstrated poor clustering by the CIMP criteria for methylation. GSEA demonstrated methylation changes in the Netrin-DCC and SLIT-ROBO pathways. Widespread changes in DNA methylation are seen in the transition from adenoma to carcinoma. The finding that GRASP, which encodes the general receptor for phosphoinositide 1-associated scaffold protein, was differentially methylated in colorectal cancer is interesting. This may be a potential biomarker for colorectal cancer.

A basal gradient of Wnt and stem-cell number influences regional tumour distribution in human and mouse intestinal tracts

OBJECTIVE

Wnt signalling is critical for normal intestinal development and homeostasis. Wnt dysregulation occurs in almost all human and murine intestinal tumours and an optimal but not excessive level of Wnt activation is considered favourable for tumourigenesis. The authors assessed effects of pan-intestinal Wnt activation on tissue homeostasis, taking into account underlying physiological Wnt activity and stem-cell number in each region of the bowel.

DESIGN

The authors generated mice that expressed temporally controlled, stabilised β-catenin along the crypt-villus axis throughout the intestines. Physiological Wnt target gene activity was assessed in different regions of normal mouse and human tissue. Human intestinal tumour mutation spectra were analysed.

RESULTS

In the mouse, β-catenin stabilisation resulted in a graduated neoplastic response, ranging from dysplastic transformation of the entire epithelium in the proximal small bowel to slightly enlarged crypts of non-dysplastic morphology in the colorectum. In contrast, stem and proliferating cell numbers were increased in all intestinal regions. In the normal mouse and human intestines, stem-cell and Wnt gradients were non-identical, but higher in the small bowel than large bowel in both species. There was also variation in the expression of some Wnt modulators. Human tumour analysis confirmed that different APC mutation spectra are selected in different regions of the bowel.

CONCLUSIONS

There are variable gradients in stem-cell number, physiological Wnt activity and response to pathologically increased Wnt signalling along the crypt-villus axis and throughout the length of the intestinal tract. The authors propose that this variation influences regional mutation spectra, tumour susceptibility and lesion distribution in mice and humans.

MED12 exon 2 mutations are common in uterine leiomyomas from South African patients

Uterine leiomyomas, or fibroids, are extremely common tumors. Regardless of their benign nature, fibroids can cause considerable morbidity. Women with African ancestry have a threefold increased risk of developing uterine leiomyomas with a greater symptom severity when compared to white women. Recently, we demonstrated that exon 2 of the MED12 gene is somatically altered in up to 70 per cent of uterine leiomyomas in a series of Finnish (Caucasian) patients. To validate these results in other populations, we sequenced a set of 28 uterine leiomyomas for MED12 exon 2 mutations from 18 different Black African or Coloured South African patients. We observed 14 mutation positive lesions (50%). When corrected by tumor size, these results are very similar to those derived in the Finnish material. This study confirms a major role of MED12 in the genesis of leiomyomas, regardless of ethnicity.

Investigation of the effects of DNA repair gene polymorphisms on the risk of colorectal cancer

Despite their prime candidate status, polymorphisms near genes involved in DNA repair or in other functions related to genome stability have been conspicuously under-represented in the significant associations reported from genome-wide association studies (GWAS) of cancer susceptibility. In this study, we assessed a set of single-nucleotide polymorphisms (SNPs) near 157 DNA repair genes in three colorectal cancer (CRC) GWAS. Although no individual SNP showed evidence of association, the set of SNPs as a whole was associated with colorectal cancer risk. When candidate SNPs were examined, our data did not support most of the previously reported associations with CRC susceptibility, an exception being an effect of the MLH1 promoter SNP -93G>A (rs1800734). Rare variants in CHEK2 (I157T and possibly del1100C) also appear to be associated with CRC risk. Overall, the absence to date of disease-associated DNA repair SNPs in cancer GWAS may be explained by a combination of the following: (i) many loci with individually very small effects on risk; (ii) rare alleles of moderate effect and (iii) subgroups of CRC, such as those with microsatellite instability, associated with specific variants. It will be particularly intriguing to determine whether any GWAS across cancer types identify DNA variants that predispose to cancers of more than one site.

Common genetic variation contributes significantly to the risk of childhood B-cell precursor acute lymphoblastic leukemia

Recent genome-wide association studies (GWAS) have provided the first unambiguous evidence that common genetic variation influences the risk of childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL), identifying risk single-nucleotide polymorphisms (SNPs) localizing to 7p12.2, 9p21.3, 10q21.2 and 14q11.2. The testing of SNPs individually for an association in GWA studies necessitates the imposition of a very stringent P-value to address the issue of multiple testing. While this reduces false positives, real associations may be missed and therefore any estimate of the total heritability will be negatively biased. Using GWAS data on 823 BCP-ALL cases by considering all typed SNPs simultaneously, we have calculated that 24% of the total variation in BCP-ALL risk is accounted for common genetic variation (95% confidence interval 6-42%). Our findings provide support for a polygenic basis for susceptibility to BCP-ALL and have wider implications for future searches for novel disease-causing risk variants.

TERC polymorphisms are associated both with susceptibility to colorectal cancer and with longer telomeres

BACKGROUND AND AIMS

Shorter telomeres have been associated with increased risk of malignancy, including colorectal cancer (CRC). Telomere length is heritable and may be an intermediate phenotype linked to genetic susceptibility to CRC.

METHODS

In a large sample, the study investigated whether candidate single nucleotide polymorphisms (SNP) in 'telomere biology' genes were associated with telomere length in leucocytes. SNP associated with an increased risk of CRC were searched for separately.

RESULTS

Carriers of the common allele at SNP rs10936599, near the telomerase RNA component (TERC) locus, had significantly longer telomeres. It was independently found that the same rs10936599 allele was associated with increased risk of both CRC and colorectal adenomas. Neither telomere length nor CRC risk was associated with variation near telomerase reverse transcriptase or other telomere biology genes. In silico analysis showed that SNP rs2293607 was strongly correlated with rs10936599, mapped within TERC transcripts, had a predicted effect on messenger RNA folding and lay at a reported transcription factor binding site. TERC mRNA were expressed, differing only at the alleles of rs2293607, in CRC cell line HCT116. The long-telomere/CRC-risk allele was associated with higher levels of TERC mRNA and the formation of longer telomeres.

CONCLUSIONS

Common genetic variation at TERC is associated with both longer telomeres and an increased risk of CRC, a potential mechanism being reduced levels of cell senescence or death. This finding is somewhat paradoxical, given retrospective studies reporting that CRC cases have shorter telomeres than controls. One possibility is that that association actually results from poorer survival in patients with longer telomeres.

Refinement of the associations between risk of colorectal cancer and polymorphisms on chromosomes 1q41 and 12q13.13

In genome-wide association studies (GWASs) of colorectal cancer, we have identified two genomic regions in which pairs of tagging-single nucleotide polymorphisms (tagSNPs) are associated with disease; these comprise chromosomes 1q41 (rs6691170, rs6687758) and 12q13.13 (rs7163702, rs11169552). We investigated these regions further, aiming to determine whether they contain more than one independent association signal and/or to identify the SNPs most strongly associated with disease. Genotyping of additional sample sets at the original tagSNPs showed that, for both regions, the two tagSNPs were unlikely to identify a single haplotype on which the functional variation lay. Conversely, one of the pair of SNPs did not fully capture the association signal in each region. We therefore undertook more detailed analyses, using imputation, logistic regression, genealogical analysis using the GENECLUSTER program and haplotype analysis. In the 1q41 region, the SNP rs11118883 emerged as a strong candidate based on all these analyses, sufficient to account for the signals at both rs6691170 and rs6687758. rs11118883 lies within a region with strong evidence of transcriptional regulatory activity and has been associated with expression of PDGFRB mRNA. For 12q13.13, a complex situation was found: SNP rs7972465 showed stronger association than either rs11169552 or rs7136702, and GENECLUSTER found no good evidence for a two-SNP model. However, logistic regression and haplotype analyses supported a two-SNP model, in which a signal at the SNP rs706793 was added to that at rs11169552. Post-GWAS fine-mapping studies are challenging, but the use of multiple tools can assist in identifying candidate functional variants in at least some cases.

Haem oxygenase is synthetically lethal with the tumour suppressor fumarate hydratase

Fumarate hydratase (FH) is an enzyme of the tricarboxylic acid cycle (TCA cycle) that catalyses the hydration of fumarate into malate. Germline mutations of FH are responsible for hereditary leiomyomatosis and renal-cell cancer (HLRCC). It has previously been demonstrated that the absence of FH leads to the accumulation of fumarate, which activates hypoxia-inducible factors (HIFs) at normal oxygen tensions. However, so far no mechanism that explains the ability of cells to survive without a functional TCA cycle has been provided. Here we use newly characterized genetically modified kidney mouse cells in which Fh1 has been deleted, and apply a newly developed computer model of the metabolism of these cells to predict and experimentally validate a linear metabolic pathway beginning with glutamine uptake and ending with bilirubin excretion from Fh1-deficient cells. This pathway, which involves the biosynthesis and degradation of haem, enables Fh1-deficient cells to use the accumulated TCA cycle metabolites and permits partial mitochondrial NADH production. We predicted and confirmed that targeting this pathway would render Fh1-deficient cells non-viable, while sparing wild-type Fh1-containing cells. This work goes beyond identifying a metabolic pathway that is induced in Fh1-deficient cells to demonstrate that inhibition of haem oxygenation is synthetically lethal when combined with Fh1 deficiency, providing a new potential target for treating HLRCC patients.

Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.

Chromosomal instability confers intrinsic multidrug resistance

Aneuploidy is associated with poor prognosis in solid tumors. Spontaneous chromosome missegregation events in aneuploid cells promote chromosomal instability (CIN) that may contribute to the acquisition of multidrug resistance in vitro and heighten risk for tumor relapse in animal models. Identification of distinct therapeutic agents that target tumor karyotypic complexity has important clinical implications. To identify distinct therapeutic approaches to specifically limit the growth of CIN tumors, we focused on a panel of colorectal cancer (CRC) cell lines, previously classified as either chromosomally unstable (CIN(+)) or diploid/near-diploid (CIN(-)), and treated them individually with a library of kinase inhibitors targeting components of signal transduction, cell cycle, and transmembrane receptor signaling pathways. CIN(+) cell lines displayed significant intrinsic multidrug resistance compared with CIN(-) cancer cell lines, and this seemed to be independent of somatic mutation status and proliferation rate. Confirming the association of CIN rather than ploidy status with multidrug resistance, tetraploid isogenic cells that had arisen from diploid cell lines displayed lower drug sensitivity than their diploid parental cells only with increasing chromosomal heterogeneity and isogenic cell line models of CIN(+) displayed multidrug resistance relative to their CIN(-) parental cancer cell line derivatives. In a meta-analysis of CRC outcome following cytotoxic treatment, CIN(+) predicted worse progression-free or disease-free survival relative to patients with CIN(-) disease. Our results suggest that stratifying tumor responses according to CIN status should be considered within the context of clinical trials to minimize the confounding effects of tumor CIN status on drug sensitivity.

Systematic search for enhancer elements and somatic allelic imbalance at seven low-penetrance colorectal cancer predisposition loci

BACKGROUND

Common single-nucleotide polymorphisms (SNPs) in ten chromosomal loci have been shown to predispose to colorectal cancer (CRC) in genome-wide association studies. A plausible biological mechanism of CRC susceptibility associated with genetic variation has so far only been proposed for three loci, each pointing to variants that affect gene expression through distant regulatory elements. In this study, we aimed to gain insight into the molecular basis of seven low-penetrance CRC loci tagged by rs4779584 at 15q13, rs10795668 at 10p14, rs3802842 at 11q23, rs4444235 at 14q22, rs9929218 at 16q22, rs10411210 at 19q13, and rs961253 at 20p12.

METHODS

Possible somatic gain of the risk allele or loss of the protective allele was studied by analyzing allelic imbalance in tumour and corresponding normal tissue samples of heterozygous patients. Functional variants were searched from in silico predicted enhancer elements locating inside the CRC-associating linkage-disequilibrium regions.

RESULTS

No allelic imbalance targeting the SNPs was observed at any of the seven loci. Altogether, 12 SNPs that were predicted to disrupt potential transcription factor binding sequences were genotyped in the same population-based case-control series as the seven tagging SNPs originally. None showed association with CRC.

CONCLUSIONS

The results of the allelic imbalance analysis suggest that the seven CRC risk variants are not somatically selected for in the neoplastic progression. The bioinformatic approach was unable to pinpoint cancer-causing variants at any of the seven loci. While it is possible that many of the predisposition loci for CRC are involved in control of gene expression by targeting transcription factor binding sites, also other possibilities, such as regulatory RNAs, should be considered.

KRAS mutation is associated with lung metastasis in patients with curatively resected colorectal cancer

PURPOSE

Oncogene mutations contribute to colorectal cancer development. We searched for differences in oncogene mutation profiles between colorectal cancer metastases from different sites and evaluated these as markers for site of relapse.

EXPERIMENTAL DESIGN

One hundred colorectal cancer metastases were screened for mutations in 19 oncogenes, and further 61 metastases and 87 matched primary cancers were analyzed for genes with identified mutations. Mutation prevalence was compared between (a) metastases from liver (n = 65), lung (n = 50), and brain (n = 46), (b) metastases and matched primary cancers, and (c) metastases and an independent cohort of primary cancers (n = 604). Mutations differing between metastasis sites were evaluated as markers for site of relapse in 859 patients from the VICTOR trial.

RESULTS

In colorectal cancer metastases, mutations were detected in 4 of 19 oncogenes: BRAF (3.1%), KRAS (48.4%), NRAS (6.2%), and PIK3CA (16.1%). KRAS mutation prevalence was significantly higher in lung (62.0%) and brain (56.5%) than in liver metastases (32.3%; P = 0.003). Mutation status was highly concordant between primary cancer and metastasis from the same individual. Compared with independent primary cancers, KRAS mutations were more common in lung and brain metastases (P < 0.005), but similar in liver metastases. Correspondingly, KRAS mutation was associated with lung relapse (HR = 2.1; 95% CI, 1.2 to 3.5, P = 0.007) but not liver relapse in patients from the VICTOR trial.

CONCLUSIONS

KRAS mutation seems to be associated with metastasis in specific sites, lung and brain, in colorectal cancer patients. Our data highlight the potential of somatic mutations for informing surveillance strategies.

Evaluation of association of HNF1B variants with diverse cancers: collaborative analysis of data from 19 genome-wide association studies

BACKGROUND

Genome-wide association studies have found type 2 diabetes-associated variants in the HNF1B gene to exhibit reciprocal associations with prostate cancer risk. We aimed to identify whether these variants may have an effect on cancer risk in general versus a specific effect on prostate cancer only.

METHODOLOGY/PRINCIPAL FINDINGS

In a collaborative analysis, we collected data from GWAS of cancer phenotypes for the frequently reported variants of HNF1B, rs4430796 and rs7501939, which are in linkage disequilibrium (r(2) = 0.76, HapMap CEU). Overall, the analysis included 16 datasets on rs4430796 with 19,640 cancer cases and 21,929 controls; and 21 datasets on rs7501939 with 26,923 cases and 49,085 controls. Malignancies other than prostate cancer included colorectal, breast, lung and pancreatic cancers, and melanoma. Meta-analysis showed large between-dataset heterogeneity that was driven by different effects in prostate cancer and other cancers. The per-T2D-risk-allele odds ratios (95% confidence intervals) for rs4430796 were 0.79 (0.76, 0.83)] per G allele for prostate cancer (p<10(-15) for both); and 1.03 (0.99, 1.07) for all other cancers. Similarly for rs7501939 the per-T2D-risk-allele odds ratios (95% confidence intervals) were 0.80 (0.77, 0.83) per T allele for prostate cancer (p<10(-15) for both); and 1.00 (0.97, 1.04) for all other cancers. No malignancy other than prostate cancer had a nominally statistically significant association.

CONCLUSIONS/SIGNIFICANCE

The examined HNF1B variants have a highly specific effect on prostate cancer risk with no apparent association with any of the other studied cancer types.

Breast cancer dormancy can be maintained by small numbers of micrometastases

Late relapse of breast cancer can occur more than 25 years after primary diagnosis. During the intervening years between initial treatment and relapse, occult cancers are maintained in an apparent state of dormancy that is poorly understood. In this study, we applied a probabilistic mathematical model to long-term follow-up studies of postresection patients to investigate the factors involved in mediating breast cancer dormancy. Our results suggest that long-term dormancy is maintained most often by just one growth-restricted dangerous micrometastasis. Analysis of the empirical data by Approximate Bayesian Computation indicated that patients in dormancy have between 1 and 5 micrometastases at 10 years postresection, when they escape growth restriction with a half-life of <69 years and are >0.4 mm in diameter. Before resection, primary tumors seed at most an average of 6 dangerous micrometastases that escape from growth restriction with a half-life of at least 12 years. Our findings suggest that effective preventive treatments will need to eliminate these small numbers of micrometastases, which may be preangiogenic and nonvascularized until they switch to growth due to one oncogenic mutation or tumor suppressor gene inactivation. In summary, breast cancer dormancy seems to be maintained by small numbers of sizeable micrometastases that escape from growth restriction with a half-life exceeding 12 years.

Low-penetrance susceptibility variants in familial colorectal cancer

BACKGROUND

Genomewide association studies have identified 10 low-penetrance loci that confer modestly increased risk for colorectal cancer (CRC). Although they underlie a significant proportion of CRC in the general population, their impact on the familial risk for CRC has yet to be formally enumerated. The aim of this study was to examine the combined contribution of the 10 variants, rs6983267, rs4779584, rs4939827, rs16892766, rs10795668, rs3802842, rs4444235, rs9929218, rs10411210, and rs961253, on familial CRC.

METHODS

The population-based series of CRC samples included in this study consisted of 97 familial cases and 691 sporadic cases. Genotypes in the 10 loci and clinical data, including family history of cancer verified from the Finnish Cancer Registry, were available. The overall number of risk alleles (0-20) was determined, and its association with familial CRC was analyzed. Excess familial risk was estimated using cancer incidence data from the first-degree relatives of the cases.

RESULTS

A linear association between the number of risk alleles and familial CRC was observed (P = 0.006). With each risk-allele addition, the odds of having an affected first-degree relative increased by 1.16 (95% confidence interval, 1.04-1.30). The 10 low-penetrance loci collectively explain approximately 9% of the variance in familial risk for CRC.

CONCLUSIONS

This study provides evidence to support the previous indirect estimations that these low-penetrance variants account for a relatively small proportion of the familial aggregation of CRC.

IMPACT

Our results emphasize the need to characterize the remaining molecular basis of familial CRC, which should eventually yield in individualized targeting of preventive interventions.

The stem cell marker CD133 associates with enhanced colony formation and cell motility in colorectal cancer

CD133 is a membrane molecule that has been, controversially, reported as a CSC marker in colorectal cancer (CRC). In this study, we sought to clarify the expression and role of CD133 in CRC. Initially the size of the CD133-expressing (CD133+) population in eight well-described CRC cell lines was measured by flow cytometry and was found to range from 0% to >95%. The cell line HT29 has a CD133+ population of >95% and was chosen for functional evaluation of CD133 after gene knockdown by RNA interference. A time course assay showed that CD133 inhibition had no significant effect on cell proliferation or apoptosis. However, CD133 knockdown did result in greater susceptibility to staurosporine-induced apoptosis (p = 0.01) and reduction in cell motility (p<0.04). Since gene knockdown may cause "off-target" effects, the cell line SW480 (which has a CD133+ population of 40%) was sorted into pure CD133+ and CD133- populations to allow functional comparison of isogenic populations separated only by CD133 expression. In concordance with the knockdown experiments, a time course assay showed no significant proliferative differences between the CD133+/CD133- populations. Also greater resistance to staurosporine-induced apoptosis (p = 0.008), greater cell motility (p = 0.03) and greater colony forming efficiency was seen in the CD133+ population than the CD133- population in both 2D and 3D culture (p<0.0001 and p<0.003 respectively). Finally, the plasticity of CD133 expression in tumour cells was tested. Quantitative PCR analysis showed there was transcriptional repression in the CD133- population of SW480. Prolonged culture of a pure CD133- population resulted in re-emergence of CD133+ cells. We conclude that CD133 expression in CRCs is associated with some features attributable to stemness and that there is plasticity of CD133 expression. Further studies are necessary to delineate the mechanistic basis of these features.

CD24 is upregulated in inflammatory bowel disease and stimulates cell motility and colony formation

BACKGROUND

We investigated whether CD24 (reportedly a stem cell marker and adhesion molecule) was expressed in regenerative mucosa in inflammatory bowel disease (IBD) and whether it could be functionally relevant.

METHODS

CD24 expression was examined in 10 cases of IBD and the relationship of CD24 with Wnt signaling was tested using dominant negative (DN)-TCF4 expression. For functional evaluation, CD24 was 1) cloned and forcibly expressed in HCT116 (which expresses very low levels of CD24) and 2) knocked-down by RNA interference in HT29 (which expresses high levels of CD24). The effect of altered CD24 expression on proliferation/apoptosis, staurosporine-induced apoptosis, colony formation in soft agar, migration, and invasion was examined.

RESULTS

CD24 was not expressed in normal tissue, while 10/10 cases of IBD showed CD24 upregulation. Inhibition of Wnt signaling with DN-TCF4 caused CD24 downregulation. Forced expression of CD24 did not influence cell proliferation, apoptosis, or staurosporine-induced apoptosis but it did significantly enhance colony forming efficiency (P < 0.01). Furthermore, there was increased transwell migration (P < 0.001) and invasion (P < 0.03) and there was increased cell migration in wounding assays. Conversely, knockdown of CD24 reduced transwell migration (P < 0.01) and invasion (P < 0.01) and reduced cell motility in wounding assays. CD24 knockdown did not influence proliferation, apoptosis resistance, or staurosporine-induced apoptosis.

CONCLUSIONS

This is the first study to report upregulation of CD24 in regenerating tissue in IBD. This may be regulated by Wnt signaling and can confer enhanced colony forming ability and enhanced cell motility-features that may be important in tissue healing in the colon.

Clonality assessment and clonal ordering of individual neoplastic crypts shows polyclonality of colorectal adenomas

BACKGROUND & AIMS

According to the somatic mutation theory, monoclonal colorectal lesions arise from sequential mutations in the progeny of a single stem cell. However, studies in a sex chromosome mixoploid mosaic (XO/XY) patient indicated that colorectal adenomas were polyclonal. We assessed adenoma clonality on an individual crypt basis and completed a genetic dependency analysis in carcinomas-in-adenomas to assess mutation order and timing.

METHODS

Polyp samples were analyzed from the XO/XY individual, patients with familial adenomatous polyposis and attenuated familial adenomatous polyposis, patients with small sporadic adenomas, and patients with sporadic carcinoma-in-adenomas. Clonality was analyzed using X/Y chromosome fluorescence in situ hybridization, analysis of 5q loss of heterozygosity in XO/XY tissue, and sequencing of adenomatous polyposis coli. Individual crypts and different phenotypic areas of carcinoma-in-adenoma lesions were analyzed for mutations in adenomatous polyposis coli, p53, and K-RAS; loss of heterozygosity at 5q, 17p, and 18q; and aneuploidy. Phylogenetic trees were constructed.

RESULTS

All familial adenomatous polyposis-associated adenomas and some sporadic lesions had polyclonal genetic defects. Some independent clones appeared to be maintained in advanced adenomas. No clear obligate order of genetic events was established. Top-down growth of dysplastic tissue into neighboring crypts was a possible mechanism of clonal competition.

CONCLUSIONS

Human colorectal microadenomas are polyclonal and may arise from a combination of host genetic features, mucosal exposures, and active crypt interactions. Analyses of tumor phylogenies show that most lesions undergo intermittent genetic homogenization, but heterotypic mutation patterns indicate that independent clonal evolution can occur throughout adenoma development. Based on observations of clonal ordering the requirement and timing of genetic events during neoplastic progression may be more variable than previously thought.

COGENT (COlorectal cancer GENeTics): an international consortium to study the role of polymorphic variation on the risk of colorectal cancer

It is now recognised that a part of the inherited risk of colorectal cancer (CRC) can be explained by the co-inheritance of low-penetrance genetic variants. The accumulated experience to date in identifying these variants has served to highlight difficulties in conducting statistically and methodologically rigorous studies and follow-up analyses. The COGENT (COlorectal cancer GENeTics) consortium includes 20 research groups in Europe, Australia, the Americas, China and Japan. The overarching goal of COGENT is to identify and characterise low-penetrance susceptibility variants for CRC through association-based analyses. In this study, we review the rationale for identifying low-penetrance variants for CRC and our proposed strategy for establishing COGENT.

DNA copy-number alterations underlie gene expression differences between microsatellite stable and unstable colorectal cancers

PURPOSE

About 15% of colorectal cancers harbor microsatellite instability (MSI). MSI-associated gene expression changes have been identified in colorectal cancers, but little overlap exists between signatures hindering an assessment of overall consistency. Little is known about the causes and downstream effects of differential gene expression.

EXPERIMENTAL DESIGN

DNA microarray data on 89 MSI and 140 microsatellite-stable (MSS) colorectal cancers from this study and 58 MSI and 77 MSS cases from three published reports were randomly divided into test and training sets. MSI-associated gene expression changes were assessed for cross-study consistency using training samples and validated as MSI classifier using test samples. Differences in biological pathways were identified by functional category analysis. Causation of differential gene expression was investigated by comparison to DNA copy-number data.

RESULTS

MSI-associated gene expression changes in colorectal cancers were found to be highly consistent across multiple studies of primary tumors and cancer cell lines from patients of different ethnicities (P < 0.001). Clustering based on consistent changes separated additional test cases by MSI status, and classification of individual samples predicted MSI status with a sensitivity of 96% and specificity of 85%. Genes associated with immune response were up-regulated in MSI cancers, whereas genes associated with cell-cell adhesion, ion binding, and regulation of metabolism were down-regulated. Differential gene expression was shown to reflect systematic differences in DNA copy-number aberrations between MSI and MSS tumors (P < 0.001).

CONCLUSIONS

Our results show cross-study consistency of MSI-associated gene expression changes in colorectal cancers. DNA copy-number alterations partly cause the differences in gene expression between MSI and MSS cancers.