Cancer drivers and clonal dynamics in acute lymphoblastic leukaemia subtypes.

To obtain a comprehensive picture of composite genetic driver events and clonal dynamics in subtypes of paediatric acute lymphoblastic leukaemia (ALL) we analysed tumour-normal whole genome sequencing and expression data from 361 newly diagnosed patients. We report the identification of both structural drivers, as well as recurrent non-coding variation in promoters. Additionally we found the transcriptional profile of histone gene cluster 1 and CTCF altered tumours shared hallmarks of hyperdiploid ALL suggesting a 'hyperdiploid like' subtype. ALL subtypes are driven by distinct mutational processes with AID mutagenesis being confined to ETV6-RUNX1 tumours. Subclonality is a ubiquitous feature of ALL, consistent with Darwinian evolution driving selection and expansion of tumours. Driver mutations in B-cell developmental genes (IKZF1, PAX5, ZEB2) tend to be clonal and RAS/RTK mutations subclonal. In addition to identifying new avenues for therapeutic exploitation, this analysis highlights that targeted therapies should take into account composite mutational profile and clonality.

Identification of four novel associations for B-cell acute lymphoblastic leukaemia risk.

There is increasing evidence for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children. To identify new risk variants for B-cell ALL (B-ALL) we conducted a meta-analysis with four GWAS (genome-wide association studies), totalling 5321 cases and 16,666 controls of European descent. We herein describe novel risk loci for B-ALL at 9q21.31 (rs76925697, P = 2.11 × 10-8), for high-hyperdiploid ALL at 5q31.1 (rs886285, P = 1.56 × 10-8) and 6p21.31 (rs210143 in BAK1, P = 2.21 × 10-8), and ETV6-RUNX1 ALL at 17q21.32 (rs10853104 in IGF2BP1, P = 1.82 × 10-8). Particularly notable are the pleiotropic effects of the BAK1 variant on multiple haematological malignancies and specific effects of IGF2BP1 on ETV6-RUNX1 ALL evidenced by both germline and somatic genomic analyses. Integration of GWAS signals with transcriptomic/epigenomic profiling and 3D chromatin interaction data for these leukaemia risk loci suggests deregulation of B-cell development and the cell cycle as central mechanisms governing genetic susceptibility to ALL.

Genetic predisposition to B-cell acute lymphoblastic leukemia at 14q11.2 is mediated by a CEBPE promoter polymorphism.

Acute lymphoblastic leukaemia (ALL) is the most common paediatric malignancy. Genome-wide association studies have shown variation at 14q11.2 influences ALL risk. We sought to decipher causal variant(s) at 14q11.2 and the mechanism of tumorigenesis. We show rs2239630 G>A resides in the promoter of the CCAT enhancer-binding protein epsilon (CEBPE) gene. The rs2239630-A risk allele is associated with increased promotor activity and CEBPE expression. Depletion of CEBPE in ALL cells reduces cell growth, correspondingly CEBPE binds to the promoters of electron transport and energy generation genes. RNA-seq in CEBPE depleted cells demonstrates CEBPE regulates the expression of genes involved in B-cell development (IL7R), apoptosis (BCL2), and methotrexate resistance (RASS4L). CEBPE regulated genes significantly overlapped in CEBPE depleted cells, ALL blasts and IGH-CEBPE translocated ALL. This suggests CEBPE regulates a similar set of genes in each, consistent with a common biological mechanism of leukemogenesis for rs2239630 associated and CEBPE translocated ALL. Finally, we map IGH-CEBPE translocation breakpoints in two cases, implicating RAG recombinase activity in their formation.

Author Correction: Genetic and regulatory mechanism of susceptibility to high-hyperdiploid acute lymphoblastic leukaemia at 10q21.2.

This corrects the article DOI: 10.1038/ncomms14616.

Genetic and regulatory mechanism of susceptibility to high-hyperdiploid acute lymphoblastic leukaemia at 10p21.2.

Despite high-hyperdiploid acute lymphoblastic leukaemia (HD-ALL) being the most common subgroup of paediatric ALL, its aetiology remains unknown. Genome-wide association studies have demonstrated association at 10q21.2. Here, we sought to determine how this region influences HD-ALL risk. We impute genotypes across the locus, finding the single nucleotide polymorphism rs7090445 highly associated with HD-ALL (P=1.54 × 10-38), and residing in a predicted enhancer element. We show this region physically interacts with the transcription start site of ARID5B, that alleles of rs7090445 have differential enhancer activity and influence RUNX3 binding. RUNX3 knock-down reduces ARID5B expression and rs7090445 enhancer activity. Individuals carrying the rs7090445-C risk allele also have reduced ARID5B expression. Finally, the rs7090445-C risk allele is preferentially retained in HD-ALL blasts consistent with inherited genetic variation contributing to arrest of normal lymphocyte development, facilitating leukaemic clonal expansion. These data provide evidence for a biological mechanism underlying hereditary risk of HD-ALL at 10q21.2.

The non-coding variant rs1800734 enhances DCLK3 expression through long-range interaction and promotes colorectal cancer progression.

Genome-wide association studies have identified a great number of non-coding risk variants for colorectal cancer (CRC). To date, the majority of these variants have not been functionally studied. Identification of allele-specific transcription factor (TF) binding is of great importance to understand regulatory consequences of such variants. A recently developed proteome-wide analysis of disease-associated SNPs (PWAS) enables identification of TF-DNA interactions in an unbiased manner. Here we perform a large-scale PWAS study to comprehensively characterize TF-binding landscape that is associated with CRC, which identifies 731 allele-specific TF binding at 116 CRC risk loci. This screen identifies the A-allele of rs1800734 within the promoter region of MLH1 as perturbing the binding of TFAP4 and consequently increasing DCLK3 expression through a long-range interaction, which promotes cancer malignancy through enhancing expression of the genes related to epithelial-to-mesenchymal transition.

Multiple myeloma risk variant at 7p15.3 creates an IRF4-binding site and interferes with CDCA7L expression.

Genome-wide association studies have identified several risk loci for multiple myeloma (MM); however, the mechanisms by which they influence MM are unknown. Here by using genetic association data and functional characterization, we demonstrate that rs4487645 G>T, the most highly associated variant (P = 5.30 × 10-25), resides in an enhancer element 47 kb upstream of the transcription start site of c-Myc-interacting CDCA7L. The G-risk allele, associated with increased CDCA7L expression (P=1.95 × 10-36), increases IRF4 binding and the enhancer interacts with the CDCA7L promoter. We show that suppression of CDCA7L limits MM proliferation through apoptosis, and increased CDCA7L expression is associated with adverse patient survival. These findings implicate IRF4-mediated CDCA7L expression in MM biology and indicate how germline variation might confer susceptibility to MM.

Genetic Predisposition to Chronic Lymphocytic Leukemia Is Mediated by a BMF Super-Enhancer Polymorphism.

Chronic lymphocytic leukemia (CLL) is an adult B cell malignancy. Genome-wide association studies show that variation at 15q15.1 influences CLL risk. We deciphered the causal variant at 15q15.1 and the mechanism by which it influences tumorigenesis. We imputed all possible genotypes across the locus and then mapped highly associated SNPs to areas of chromatin accessibility, evolutionary conservation, and transcription factor binding. SNP rs539846 C>A, the most highly associated variant (p = 1.42 × 10(-13), odds ratio = 1.35), localizes to a super-enhancer defined by extensive histone H3 lysine 27 acetylation in intron 3 of B cell lymphoma 2 (BCL2)-modifying factor (BMF). The rs539846-A risk allele alters a conserved RELA-binding motif, disrupts RELA binding, and is associated with decreased BMF expression in CLL. These findings are consistent with rs539846 influencing CLL susceptibility through differential RELA binding, with direct modulation of BMF expression impacting on anti-apoptotic BCL2, a hallmark of oncogenic dependency in CLL.

The role of interleukin-8 (IL-8) and IL-8 receptors in platinum response in high grade serous ovarian carcinoma.

Platinum based drugs are the cornerstone of chemotherapy for ovarian cancer, however the development of chemoresistance hinders its success. IL-8 is involved in regulating several pro-survival pathways in cancer. We studied the expression of IL-8 and IL-8 receptors in platinum sensitive and resistant cell lines. Using qRT-PCR and immunohistochemistry, both platinum sensitive (PEA1, PEO14) and resistant (PEA2, PEO23) show increased expression of IL-8 and IL-8 receptors. IL-8RA shows nuclear and cytoplasmic expression, whilst IL-8RB is present solely in the cytoplasm. Knockdown of IL-8 increased sensitivity to cisplatin in platinum sensitive and reversed platinum resistance in resistant cell lines, decreased the expression of anti-apoptotic Bcl-2 and decreased inhibitory phosphorylation of pro-apoptotic Bad. IL-8 receptor antagonist treatment also enhanced platinum sensitivity. Nuclear localisation of IL-8RA was only detected in platinum resistant tumours. Inhibition of IL-8 signalling can enhance response in platinum sensitive and resistant disease. Nuclear IL-8RA may have potential as a biomarker of resistant disease.

HDAC4-regulated STAT1 activation mediates platinum resistance in ovarian cancer.

Ovarian cancer frequently acquires resistance to platinum chemotherapy, representing a major challenge for improving patient survival. Recent work suggests that resistant clones exist within a larger drug-sensitive cell population prior to chemotherapy, implying that resistance is selected for rather than generated by treatment. We sought to compare clinically derived, intrapatient paired models of initial platinum response and subsequent resistant relapse to define molecular determinants of evolved resistance. Transcriptional analysis of a matched cell line series from three patients with high-grade serous ovarian cancer before and after development of clinical platinum resistance (PEO1/PEO4/PEO6, PEA1/PEA2, PEO14/PEO23) identified 91 up- and 126 downregulated genes common to acquired resistance. Significantly enhanced apoptotic response to platinum treatment in resistant cells was observed following knockdown of histone deacetylase (HDAC) 4, FOLR2, PIK3R1, or STAT1 (P < 0.05). Interestingly, HDAC4 and STAT1 were found to physically interact. Acetyl-STAT1 was detected in platinum-sensitive cells but not in HDAC4 overexpressing platinum-resistant cells from the same patient. In resistant cells, STAT1 phosphorylation/nuclear translocation was seen following platinum exposure, whereas silencing of HDAC4 increased acetyl-STAT1 levels, prevented platinum-induced STAT1 activation, and restored cisplatin sensitivity. Conversely, matched sensitive cells were refractory to STAT1 phosphorylation on platinum treatment. Analysis of 16 paired tumor biopsies taken before and after development of clinical platinum resistance showed significantly increased HDAC4 expression in resistant tumors [n = 7 of 16 (44%); P = 0.04]. Therefore, clinical selection of HDAC4-overexpressing tumor cells upon exposure to chemotherapy promotes STAT1 deacetylation and cancer cell survival. Together, our findings identify HDAC4 as a novel, therapeutically tractable target to counter platinum resistance in ovarian cancer.