The genetic basis for cancer treatment decisions.

Personalized cancer medicine is based on increased knowledge of the cancer mutation repertoire and availability of agents that target altered genes or pathways. Given advances in cancer genetics, technology, and therapeutics development, the timing is right to develop a clinical trial and research framework to move future clinical decisions from heuristic to evidence-based decisions. Although the challenges of integrating genomic testing into cancer treatment decision making are wide-ranging and complex, there is a scientific and ethical imperative to realize the benefits of personalized cancer medicine, given the overwhelming burden of cancer and the unprecedented opportunities for advancements in outcomes for patients.

Tracing the origins of relapse in acute myeloid leukaemia to stem cells.

In acute myeloid leukaemia, long-term survival is poor as most patients relapse despite achieving remission. Historically, the failure of therapy has been thought to be due to mutations that produce drug resistance, possibly arising as a consequence of the mutagenic properties of chemotherapy drugs. However, other lines of evidence have pointed to the pre-existence of drug-resistant cells. For example, deep sequencing of paired diagnosis and relapse acute myeloid leukaemia samples has provided direct evidence that relapse in some cases is generated from minor genetic subclones present at diagnosis that survive chemotherapy, suggesting that resistant cells are generated by evolutionary processes before treatment and are selected by therapy. Nevertheless, the mechanisms of therapy failure and capacity for leukaemic regeneration remain obscure, as sequence analysis alone does not provide insight into the cell types that are fated to drive relapse. Although leukaemia stem cells have been linked to relapse owing to their dormancy and self-renewal properties, and leukaemia stem cell gene expression signatures are highly predictive of therapy failure, experimental studies have been primarily correlative and a role for leukaemia stem cells in acute myeloid leukaemia relapse has not been directly proved. Here, through combined genetic and functional analysis of purified subpopulations and xenografts from paired diagnosis/relapse samples, we identify therapy-resistant cells already present at diagnosis and two major patterns of relapse. In some cases, relapse originated from rare leukaemia stem cells with a haematopoietic stem/progenitor cell phenotype, while in other instances relapse developed from larger subclones of immunophenotypically committed leukaemia cells that retained strong stemness transcriptional signatures. The identification of distinct patterns of relapse should lead to improved methods for disease management and monitoring in acute myeloid leukaemia. Moreover, the shared functional and transcriptional stemness properties that underlie both cellular origins of relapse emphasize the importance of developing new therapeutic approaches that target stemness to prevent relapse.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns.

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

A balance score between immune stimulatory and suppressive microenvironments identifies mediators of tumour immunity and predicts pan-cancer survival.

BACKGROUND: The balance between immune-stimulatory and immune-suppressive mechanisms in the tumour microenvironment is associated with tumour rejection and can predict the efficacy of immune checkpoint-inhibition therapies. METHODS: We consider the observed differences between the transcriptional programmes associated with cancer types where the levels of immune infiltration predict a favourable prognosis versus those in which the immune infiltration predicts an unfavourable prognosis and defined a score named Mediators of Immune Response Against Cancer in soLid microEnvironments (MIRACLE). MIRACLE deconvolves T cell infiltration, from inhibitory mechanisms, such as TGFß, EMT and PI3Kγ signatures. RESULTS: Our score outperforms current state-of-the-art immune signatures as a predictive marker of survival in TCGA (n = 9305, HR: 0.043, p value: 6.7 × 10-36). In a validation cohort (n = 7623), MIRACLE predicts better survival compared to other immune metrics (HR: 0.1985, p value: 2.73 × 10-38). MIRACLE also predicts response to checkpoint-inhibitor therapies (n = 333). The tumour-intrinsic factors inversely associated with the reported score such as EGFR, PRKAR1A and MAP3K1 are frequently associated with immune-suppressive phenotypes. CONCLUSIONS: The association of cancer outcome with the level of infiltrating immune cells is mediated by the balance of activatory and suppressive factors. MIRACLE accounts for this balance and predicts favourable cancer outcomes.

Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia.

Targeted sequencing of 103 leukemia-associated genes in leukemia cells from 841 treatment-naive patients with chronic lymphocytic leukemia (CLL) identified 89 (11%) patients as having CLL cells with mutations in genes encoding proteins that putatively are involved in hedgehog (Hh) signaling. Consistent with this finding, there was a significant association between the presence of these mutations and the expression of GLI1 (χ2 test, P < .0001), reflecting activation of the Hh pathway. However, we discovered that 38% of cases without identified mutations also were GLI1+ Patients with GLI1+ CLL cells had a shorter median treatment-free survival than patients with CLL cells lacking expression of GLI1 independent of IGHV mutation status. We found that GANT61, a small molecule that can inhibit GLI1, was highly cytotoxic for GLI1+ CLL cells relative to that of CLL cells without GLI1. Collectively, this study shows that a large proportion of patients have CLL cells with activated Hh signaling, which is associated with early disease progression and enhanced sensitivity to inhibition of GLI1.

An epigenome-wide association study of total serum immunoglobulin E concentration.

Immunoglobulin E (IgE) is a central mediator of allergic (atopic) inflammation. Therapies directed against IgE can alleviate hay fever and allergic asthma. Genetic association studies have not yet identified novel therapeutic targets or pathways underlying IgE regulation. We therefore surveyed epigenetic associations between serum IgE concentrations and methylation at loci concentrated in CpG islands genome wide in 95 nuclear pedigrees, using DNA from peripheral blood leukocytes. We validated positive results in additional families and in subjects from the general population. Here we show replicated associations--with a meta-analysis false discovery rate less than 10(-4)--between IgE and low methylation at 36 loci. Genes annotated to these loci encode known eosinophil products, and also implicate phospholipid inflammatory mediators, specific transcription factors and mitochondrial proteins. We confirmed that methylation at these loci differed significantly in isolated eosinophils from subjects with and without asthma and high IgE levels. The top three loci accounted for 13% of IgE variation in the primary subject panel, explaining the tenfold higher variance found compared with that derived from large single-nucleotide polymorphism genome-wide association studies. This study identifies novel therapeutic targets and biomarkers for patient stratification for allergic diseases.

Identification of pre-leukaemic haematopoietic stem cells in acute leukaemia.

In acute myeloid leukaemia (AML), the cell of origin, nature and biological consequences of initiating lesions, and order of subsequent mutations remain poorly understood, as AML is typically diagnosed without observation of a pre-leukaemic phase. Here, highly purified haematopoietic stem cells (HSCs), progenitor and mature cell fractions from the blood of AML patients were found to contain recurrent DNMT3A mutations (DNMT3A(mut)) at high allele frequency, but without coincident NPM1 mutations (NPM1c) present in AML blasts. DNMT3A(mut)-bearing HSCs showed a multilineage repopulation advantage over non-mutated HSCs in xenografts, establishing their identity as pre-leukaemic HSCs. Pre-leukaemic HSCs were found in remission samples, indicating that they survive chemotherapy. Therefore DNMT3A(mut) arises early in AML evolution, probably in HSCs, leading to a clonally expanded pool of pre-leukaemic HSCs from which AML evolves. Our findings provide a paradigm for the detection and treatment of pre-leukaemic clones before the acquisition of additional genetic lesions engenders greater therapeutic resistance.

Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer.

Genome-wide association studies have reported 56 independently associated colorectal cancer (CRC) risk variants, most of which are non-coding and believed to exert their effects by modulating gene expression. The computational method PrediXcan uses cis-regulatory variant predictors to impute expression and perform gene-level association tests in GWAS without directly measured transcriptomes. In this study, we used reference datasets from colon (n = 169) and whole blood (n = 922) transcriptomes to test CRC association with genetically determined expression levels in a genome-wide analysis of 12,186 cases and 14,718 controls. Three novel associations were discovered from colon transverse models at FDR ≤ 0.2 and further evaluated in an independent replication including 32,825 cases and 39,933 controls. After adjusting for multiple comparisons, we found statistically significant associations using colon transcriptome models with TRIM4 (discovery P = 2.2 × 10- 4, replication P = 0.01), and PYGL (discovery P = 2.3 × 10- 4, replication P = 6.7 × 10- 4). Interestingly, both genes encode proteins that influence redox homeostasis and are related to cellular metabolic reprogramming in tumors, implicating a novel CRC pathway linked to cell growth and proliferation. Defining CRC risk regions as one megabase up- and downstream of one of the 56 independent risk variants, we defined 44 non-overlapping CRC-risk regions. Among these risk regions, we identified genes associated with CRC (P < 0.05) in 34/44 CRC-risk regions. Importantly, CRC association was found for two genes in the previously reported 2q25 locus, CXCR1 and CXCR2, which are potential cancer therapeutic targets. These findings provide strong candidate genes to prioritize for subsequent laboratory follow-up of GWAS loci. This study is the first to implement PrediXcan in a large colorectal cancer study and findings highlight the utility of integrating transcriptome data in GWAS for discovery of, and biological insight into, risk loci.

Correction to: Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer.

Every author has erroneously been assigned to the affiliation "62". The affiliation 62 belongs to the author Graham Casey.

Genome-Wide Interaction Analyses between Genetic Variants and Alcohol Consumption and Smoking for Risk of Colorectal Cancer.

Genome-wide association studies (GWAS) have identified many genetic susceptibility loci for colorectal cancer (CRC). However, variants in these loci explain only a small proportion of familial aggregation, and there are likely additional variants that are associated with CRC susceptibility. Genome-wide studies of gene-environment interactions may identify variants that are not detected in GWAS of marginal gene effects. To study this, we conducted a genome-wide analysis for interaction between genetic variants and alcohol consumption and cigarette smoking using data from the Colon Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Interactions were tested using logistic regression. We identified interaction between CRC risk and alcohol consumption and variants in the 9q22.32/HIATL1 (Pinteraction = 1.76×10-8; permuted p-value 3.51x10-8) region. Compared to non-/occasional drinking light to moderate alcohol consumption was associated with a lower risk of colorectal cancer among individuals with rs9409565 CT genotype (OR, 0.82 [95% CI, 0.74-0.91]; P = 2.1×10-4) and TT genotypes (OR,0.62 [95% CI, 0.51-0.75]; P = 1.3×10-6) but not associated among those with the CC genotype (p = 0.059). No genome-wide statistically significant interactions were observed for smoking. If replicated our suggestive finding of a genome-wide significant interaction between genetic variants and alcohol consumption might contribute to understanding colorectal cancer etiology and identifying subpopulations with differential susceptibility to the effect of alcohol on CRC risk.

Assessment of TP53 Polymorphisms and MDM2 SNP309 in Premenopausal Breast Cancer Risk.

Germline polymorphic variants in cancer predisposition genes such as TP53 have been shown to impact the risk of premenopausal cancer. Accordingly, the aim of this study was to assess the spectrum of polymorphisms in TP53 and its negative regulatory gene, MDM2 (SNP309:T>G) in patients with premenopausal breast cancer. Our findings in a cohort of 40 female patients demonstrate no significant correlation between the studied polymorphisms and risk of premenopausal breast cancer. Although one polymorphism is found in high frequency in this cohort (rs1800372:A>G, 9.0%), it was not associated with the risk of developing cancer before the age of 35 years in an extended cohort of 1,420 breast cancer cases. Functional studies of the rs1800372:A>G polymorphic allele reveal that it does not affect p53 transactivation function. Further study of variants or mutations in other cancer susceptibility genes is warranted to refine our understanding of the germline contribution to premenopausal breast cancer susceptibility.

Genetic determinants of telomere length and risk of common cancers: a Mendelian randomization study.

Epidemiological studies have reported inconsistent associations between telomere length (TL) and risk for various cancers. These inconsistencies are likely attributable, in part, to biases that arise due to post-diagnostic and post-treatment TL measurement. To avoid such biases, we used a Mendelian randomization approach and estimated associations between nine TL-associated SNPs and risk for five common cancer types (breast, lung, colorectal, ovarian and prostate cancer, including subtypes) using data on 51 725 cases and 62 035 controls. We then used an inverse-variance weighted average of the SNP-specific associations to estimate the association between a genetic score representing long TL and cancer risk. The long TL genetic score was significantly associated with increased risk of lung adenocarcinoma (P = 6.3 × 10(-15)), even after exclusion of a SNP residing in a known lung cancer susceptibility region (TERT-CLPTM1L) P = 6.6 × 10(-6)). Under Mendelian randomization assumptions, the association estimate [odds ratio (OR) = 2.78] is interpreted as the OR for lung adenocarcinoma corresponding to a 1000 bp increase in TL. The weighted TL SNP score was not associated with other cancer types or subtypes. Our finding that genetic determinants of long TL increase lung adenocarcinoma risk avoids issues with reverse causality and residual confounding that arise in observational studies of TL and disease risk. Under Mendelian randomization assumptions, our finding suggests that longer TL increases lung adenocarcinoma risk. However, caution regarding this causal interpretation is warranted in light of the potential issue of pleiotropy, and a more general interpretation is that SNPs influencing telomere biology are also implicated in lung adenocarcinoma risk.

Functional organization of an Mbp enhancer exposes striking transcriptional regulatory diversity within myelinating glia.

In mammals, large caliber axons are ensheathed by myelin, a glial specialization supporting axon integrity and conferring accelerated and energy-efficient action potential conduction. Myelin basic protein (MBP) is required for normal myelin elaboration with maximal mbp transcription in oligodendrocytes requiring the upstream M3 enhancer. To further characterize the mechanism regulating mbp transcription, we defined M3 structure/function relationships by evaluating its evolutionary conservation, DNA footprints and the developmental programing conferred in mice by M3 derivatives. Multiple M3 regulatory element combinations were found to drive expression in oligodendrocytes and Schwann cells with a minimal 129 bp sequence conferring expression in oligodendrocytes throughout myelin elaboration, maintenance and repair. Unexpectedly, M3 derivatives conferred markedly different spatial and temporal expression programs thus illuminating striking transcriptional heterogeneity within post-mitotic oligodendrocytes. Finally, one M3 derivative engaged only during primary myelination, not during adult remyelination, demonstrating that transcriptional regulation in the two states is not equivalent.

CYP24A1 variant modifies the association between use of oestrogen plus progestogen therapy and colorectal cancer risk.

BACKGROUND: Menopausal hormone therapy (MHT) use has been consistently associated with a decreased risk of colorectal cancer (CRC) in women. Our aim was to use a genome-wide gene-environment interaction analysis to identify genetic modifiers of CRC risk associated with use of MHT. METHODS: We included 10 835 postmenopausal women (5419 cases and 5416 controls) from 10 studies. We evaluated use of any MHT, oestrogen-only (E-only) and combined oestrogen-progestogen (E+P) hormone preparations. To test for multiplicative interactions, we applied the empirical Bayes (EB) test as well as the Wald test in conventional case-control logistic regression as primary tests. The Cocktail test was used as secondary test. RESULTS: The EB test identified a significant interaction between rs964293 at 20q13.2/CYP24A1 and E+P (interaction OR (95% CIs)=0.61 (0.52-0.72), P=4.8 × 10(-9)). The secondary analysis also identified this interaction (Cocktail test OR=0.64 (0.52-0.78), P=1.2 × 10(-5) (alpha threshold=3.1 × 10(-4)). The ORs for association between E+P and CRC risk by rs964293 genotype were as follows: C/C, 0.96 (0.61-1.50); A/C, 0.61 (0.39-0.95) and A/A, 0.40 (0.22-0.73), respectively. CONCLUSIONS: Our results indicate that rs964293 modifies the association between E+P and CRC risk. The variant is located near CYP24A1, which encodes an enzyme involved in vitamin D metabolism. This novel finding offers additional insight into downstream pathways of CRC etiopathogenesis.

International network of cancer genome projects.

The International Cancer Genome Consortium (ICGC) was launched to coordinate large-scale cancer genome studies in tumours from 50 different cancer types and/or subtypes that are of clinical and societal importance across the globe. Systematic studies of more than 25,000 cancer genomes at the genomic, epigenomic and transcriptomic levels will reveal the repertoire of oncogenic mutations, uncover traces of the mutagenic influences, define clinically relevant subtypes for prognosis and therapeutic management, and enable the development of new cancer therapies.

The effect of 5-fluorouracil/leucovorin chemotherapy on CpG methylation, or the confounding role of leukocyte heterogeneity: An illustration.

Blood-based epigenome-wide association studies that aim at comparing CpG methylation between colorectal cancer (CRC) patients and controls can lead to the discovery of diagnostic or prognostic biomarkers. Numerous confounders can lead to spurious associations. We aimed to see if 5-fluorouracil (5-FU)/leucovorin chemotherapy administered to cases prior to the collection of their blood has an effect on methylation. 304 patients who received treatment and 273 who did not were profiled on the HumanMethylation450 array. Association tests were adjusted for confounders, including proxies for leukocyte cell counts. There were substantial methylation differences between these two groups that vanished once the leukocyte heterogeneity was accounted for. We observed a significant decrease of T cells in the treatment group (CD4+: p=10(-6); CD8+: p=0.036) and significant increase of NK cells (p=0.05) and monocytes (p=0.0006). 5-FU/leucovorin has no effect on global and local blood-based methylation profiles, other than through differences in the leukocyte compositions that the treatment induced.

Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24.

Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage 1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 x 10(-8)) and 1.14 (trend; P = 1.32 x 10(-5)), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 x 10(-4)). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 x 10(-11)). This locus has also been implicated in prostate cancer.

A genome-wide association study for colorectal cancer identifies a risk locus in 14q23.1.

Over 50 loci associated with colorectal cancer (CRC) have been uncovered by genome-wide association studies (GWAS). Identifying additional loci has the potential to help elucidate aspects of the underlying biological processes leading to better understanding of the pathogenesis of the disease. We re-evaluated a GWAS by excluding controls that have family history of CRC or personal history of colorectal polyps, as we hypothesized that their inclusion reduces power to detect associations. This is supported empirically and through simulations. Two-phase GWAS analysis was performed in a total of 16,517 cases and 14,487 controls. We identified rs17094983, a SNP associated with risk of CRC [p = 2.5 × 10(-10); odds ratio estimated by re-including all controls (OR) = 0.87, 95% confidence interval (CI) 0.83-0.91; minor allele frequency (MAF) = 13%]. Results were replicated in samples of African descent (1894 cases and 4703 controls; p = 0.01; OR = 0.86, 95% CI 0.77-0.97; MAF = 16 %). Gene expression data in 195 colon adenocarcinomas and 59 normal colon tissues from two different studies revealed that this locus has genotypes that are associated with RTN1 (Reticulon 1) expression (p = 0.001), a protein-coding gene involved in survival and proliferation of cancer cells which is highly expressed in normal colon tissues but has significantly reduced expression in tumor cells (p = 1.3 × 10(-8)).