Phenome-wide Mendelian randomisation analysis of 378,142 cases reveals risk factors for eight common cancers.

For many cancers there are only a few well-established risk factors. Here, we use summary data from genome-wide association studies (GWAS) in a Mendelian randomisation (MR) phenome-wide association study (PheWAS) to identify potentially causal relationships for over 3,000 traits. Our outcome datasets comprise 378,142 cases across breast, prostate, colorectal, lung, endometrial, oesophageal, renal, and ovarian cancers, as well as 485,715 controls. We complement this analysis by systematically mining the literature space for supporting evidence. In addition to providing supporting evidence for well-established risk factors (smoking, alcohol, obesity, lack of physical activity), we also find sex steroid hormones, plasma lipids, and telomere length as determinants of cancer risk. A number of the molecular factors we identify may prove to be potential biomarkers. Our analysis, which highlights aetiological similarities and differences in common cancers, should aid public health prevention strategies to reduce cancer burden. We provide a R/Shiny app to visualise findings.

Risk factors for eight common cancers revealed from a phenome-wide Mendelian randomisation analysis of 378,142 cases and 485,715 controls.

For many cancers there are few well-established risk factors. Summary data from genome-wide association studies (GWAS) can be used in a Mendelian randomisation (MR) phenome-wide association study (PheWAS) to identify causal relationships. We performed a MR-PheWAS of breast, prostate, colorectal, lung, endometrial, oesophageal, renal, and ovarian cancers, comprising 378,142 cases and 485,715 controls. To derive a more comprehensive insight into disease aetiology we systematically mined the literature space for supporting evidence. We evaluated causal relationships for over 3,000 potential risk factors. In addition to identifying well-established risk factors (smoking, alcohol, obesity, lack of physical activity), we provide evidence for specific factors, including dietary intake, sex steroid hormones, plasma lipids and telomere length as determinants of cancer risk. We also implicate molecular factors including plasma levels of IL-18, LAG-3, IGF-1, CT-1, and PRDX1 as risk factors. Our analyses highlight the importance of risk factors that are common to many cancer types but also reveal aetiological differences. A number of the molecular factors we identify have the potential to be biomarkers. Our findings should aid public health prevention strategies to reduce cancer burden. We provide a R/Shiny app (https://mrcancer.shinyapps.io/mrcan/) to visualise findings.

Risk factors for eight common cancers revealed from a phenome-wide Mendelian randomisation analysis of 378,142 cases and 485,715 controls.

For many cancers there are few well-established risk factors. Summary data from genome-wide association studies (GWAS) can be used in a Mendelian randomisation (MR) phenome-wide association study (PheWAS) to identify causal relationships. We performed a MR-PheWAS of breast, prostate, colorectal, lung, endometrial, oesophageal, renal, and ovarian cancers, comprising 378,142 cases and 485,715 controls. To derive a more comprehensive insight into disease aetiology we systematically mined the literature space for supporting evidence. We evaluated causal relationships for over 3,000 potential risk factors. In addition to identifying well-established risk factors (smoking, alcohol, obesity, lack of physical activity), we provide evidence for specific factors, including dietary intake, sex steroid hormones, plasma lipids and telomere length as determinants of cancer risk. We also implicate molecular factors including plasma levels of IL-18, LAG-3, IGF-1, CT-1, and PRDX1 as risk factors. Our analyses highlight the importance of risk factors that are common to many cancer types but also reveal aetiological differences. A number of the molecular factors we identify have the potential to be biomarkers. Our findings should aid public health prevention strategies to reduce cancer burden. We provide a R/Shiny app (https://mrcancer.shinyapps.io/mrcan/) to visualise findings.

Chemotherapy-related hyperbilirubinemia in pediatric acute lymphoblastic leukemia: a genome-wide association study from the AIEOP-BFM ALL study group.

BACKGROUND: Characterization of clinical phenotypes in context with tumor and host genomic information can aid in the development of more effective and less toxic risk-adapted and targeted treatment strategies. To analyze the impact of therapy-related hyperbilirubinemia on treatment outcome and to identify contributing genetic risk factors of this well-recognized adverse effect we evaluated serum bilirubin levels in 1547 pediatric patients with acute lymphoblastic leukemia (ALL) and conducted a genome-wide association study (GWAS). PATIENTS AND METHODS: Patients were treated in multicenter trial AIEOP-BFM ALL 2000 for pediatric ALL. Bilirubin toxicity was graded 0 to 4 according to the Common Toxicity Criteria (CTC) of the National Cancer Institute. In the GWAS discovery cohort, including 650 of the 1547 individuals, genotype frequencies of 745,895 single nucleotide variants were compared between 435 patients with hyperbilirubinemia (CTC grades 1-4) during induction/consolidation treatment and 215 patients without it (grade 0). Replication analyses included 224 patients from the same trial. RESULTS: Compared to patients with no (grade 0) or moderate hyperbilirubinemia (grades 1-2) during induction/consolidation, patients with grades 3-4 had a poorer 5-year event free survival (76.6 ± 3% versus 87.7 ± 1% for grades 1-2, P = 0.003; 85.2 ± 2% for grade 0, P < 0.001) and a higher cumulative incidence of relapse (15.6 ± 3% versus 9.0 ± 1% for grades 1-2, P = 0.08; 11.1 ± 1% for grade 0, P = 0.007). GWAS identified a strong association of the rs6744284 variant T allele in the UGT1A gene cluster with risk of hyperbilirubinemia (allelic odds ratio (OR) = 2.1, P = 7 × 10- 8). TT-homozygotes had a 6.5-fold increased risk of hyperbilirubinemia (grades 1-4; 95% confidence interval (CI) = 2.9-14.6, P = 7 × 10- 6) and a 16.4-fold higher risk of grade 3-4 hyperbilirubinemia (95% CI 6.1-43.8, P = 2 × 10- 8). Replication analyses confirmed these associations with joint analysis yielding genome-wide significance (allelic OR = 2.1, P = 6 × 10- 11; 95% CI 1.7-2.7). Moreover, rs6744284 genotypes were strongly linked to the Gilbert's syndrome-associated UGT1A1*28/*37 allele (r2 = 0.70), providing functional support for study findings. Of clinical importance, the rs6744284 TT genotype counterbalanced the adverse prognostic impact of high hyperbilirubinemia on therapy outcome. CONCLUSIONS: Chemotherapy-related hyperbilirubinemia is a prognostic factor for treatment outcome in pediatric ALL and genetic variation in UGT1A aids in predicting the clinical impact of hyperbilirubinemia. TRIAL REGISTRATION: http://www. CLINICALTRIALS: gov ; #NCT00430118.

Leveraging Genome and Phenome-Wide Association Studies to Investigate Genetic Risk of Acute Lymphoblastic Leukemia.

BACKGROUND: Genome-wide association studies (GWAS) of childhood cancers remain limited, highlighting the need for novel analytic strategies. We describe a hybrid GWAS and phenome-wide association study (PheWAS) approach to uncover genotype-phenotype relationships and candidate risk loci, applying it to acute lymphoblastic leukemia (ALL). METHODS: PheWAS was performed for 12 ALL SNPs identified by prior GWAS and two control SNP-sets using UK Biobank data. PheWAS-traits significantly associated with ALL SNPs compared with control SNPs were assessed for association with ALL risk (959 cases, 2,624 controls) using polygenic score and Mendelian randomization analyses. Trait-associated SNPs were tested for association with ALL risk in single-SNP analyses, with replication in an independent case-control dataset (1,618 cases, 9,409 controls). RESULTS: Platelet count was the trait most enriched for association with known ALL risk loci. A polygenic score for platelet count (223 SNPs) was not associated with ALL risk (P = 0.82) and Mendelian randomization did not suggest a causal relationship. However, twelve platelet count-associated SNPs were nominally associated with ALL risk in COG data and three were replicated in UK data (rs10058074, rs210142, rs2836441). CONCLUSIONS: In our hybrid GWAS-PheWAS approach, we identify pleiotropic genetic variation contributing to ALL risk and platelet count. Three SNPs known to influence platelet count were reproducibly associated with ALL risk, implicating genomic regions containing IRF1, proapoptotic protein BAK1, and ERG in platelet production and leukemogenesis. IMPACT: Incorporating PheWAS data into association studies can leverage genetic pleiotropy to identify cancer risk loci, highlighting the utility of our novel approach.

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.

Author Correction: Genome-wide association study identifies susceptibility loci for B-cell childhood acute lymphoblastic leukemia.

The original version of this Article contained an error in the spelling of a member of the PRACTICAL Consortium, Manuela Gago-Dominguez, which was incorrectly given as Manuela Gago Dominguez. This has now been corrected in both the PDF and HTML versions of the Article. Furthermore, in the original HTML version of this Article, the order of authors within the author list was incorrect. The PRACTICAL consortium was incorrectly listed after Richard S. Houlston and should have been listed after Nora Pashayan. This error has been corrected in the HTML version of the Article; the PDF version was correct at the time of publication.

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.

Genome-wide association study identifies susceptibility loci for B-cell childhood acute lymphoblastic leukemia.

Genome-wide association studies (GWAS) have advanced our understanding of susceptibility to B-cell precursor acute lymphoblastic leukemia (BCP-ALL); however, much of the heritable risk remains unidentified. Here, we perform a GWAS and conduct a meta-analysis with two existing GWAS, totaling 2442 cases and 14,609 controls. We identify risk loci for BCP-ALL at 8q24.21 (rs28665337, P = 3.86 × 10-9, odds ratio (OR) = 1.34) and for ETV6-RUNX1 fusion-positive BCP-ALL at 2q22.3 (rs17481869, P = 3.20 × 10-8, OR = 2.14). Our findings provide further insights into genetic susceptibility to ALL and its biology.

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.

Genome-wide association study identifies multiple susceptibility loci for glioma.

Previous genome-wide association studies (GWASs) have shown that common genetic variation contributes to the heritable risk of glioma. To identify new glioma susceptibility loci, we conducted a meta-analysis of four GWAS (totalling 4,147 cases and 7,435 controls), with imputation using 1000 Genomes and UK10K Project data as reference. After genotyping an additional 1,490 cases and 1,723 controls we identify new risk loci for glioblastoma (GBM) at 12q23.33 (rs3851634, near POLR3B, P=3.02 × 10(-9)) and non-GBM at 10q25.2 (rs11196067, near VTI1A, P=4.32 × 10(-8)), 11q23.2 (rs648044, near ZBTB16, P=6.26 × 10(-11)), 12q21.2 (rs12230172, P=7.53 × 10(-11)) and 15q24.2 (rs1801591, near ETFA, P=5.71 × 10(-9)). Our findings provide further insights into the genetic basis of the different glioma subtypes.

The 9p21.3 risk of childhood acute lymphoblastic leukaemia is explained by a rare high-impact variant in CDKN2A.

Genome-wide association studies (GWAS) have provided strong evidence for inherited predisposition to childhood acute lymphoblastic leukaemia (ALL) identifying a number of risk loci. We have previously shown common SNPs at 9p21.3 influence ALL risk. These SNP associations are generally not themselves candidates for causality, but simply act as markers for functional variants. By means of imputation of GWAS data and subsequent validation SNP genotyping totalling 2,177 ALL cases and 8,240 controls, we have shown that the 9p21.3 association can be ascribed to the rare high-impact CDKN2A p.Ala148Thr variant (rs3731249; Odds ratio = 2.42, P = 3.45 × 10(-19)). The association between rs3731249 genotype and risk was not specific to particular subtype of B-cell ALL. The rs3731249 variant is associated with predominant nuclear localisation of the CDKN2A transcript suggesting the functional effect of p.Ala148Thr on ALL risk may be through compromised ability to inhibit cyclin D within the cytoplasm.

Inherited genetic susceptibility to monoclonal gammopathy of unknown significance.

Monoclonal gammopathy of undetermined significance (MGUS) is present in ∼2% of individuals age >50 years. The increased risk of multiple myeloma (MM) in relatives of individuals with MGUS is consistent with MGUS being a marker of inherited genetic susceptibility to MM. Common single-nucleotide polymorphisms (SNPs) at 2p23.3 (rs6746082), 3p22.1 (rs1052501), 3q26.2 (rs10936599), 6p21.33 (rs2285803), 7p15.3 (rs4487645), 17p11.2 (rs4273077), and 22q13.1 (rs877529) have recently been shown to influence MM risk. To examine the impact of these 7 SNPs on MGUS, we analyzed two case-control series totaling 492 cases and 7306 controls. Each SNP independently influenced MGUS risk with statistically significant associations (P < .02) for rs1052501, rs2285803, rs4487645, and rs4273077. SNP associations were independent, with risk increasing with a larger number of risk alleles carried (per allele odds ratio, 1.18; P < 10(-7)). Collectively these data are consistent with a polygenic model of disease susceptibility to MGUS.

Variation at 3p24.1 and 6q23.3 influences the risk of Hodgkin's lymphoma.

In addition to HLA, recent genome-wide association studies (GWASs) of Hodgkin's lymphoma (HL) have identified susceptibility loci for HL at 2p16.1, 8q24.21 and 10p14. In this study, we perform a GWAS meta-analysis with published GWAS (totalling 1,465 cases and 6,417 controls of European background), and follow-up the most significant association signals in 2,024 cases and 1,853 controls. A combined analysis identifies new HL susceptibility loci mapping to 3p24.1 (rs3806624; P=1.14 × 10(-12), odds ratio (OR)=1.26) and 6q23.3 (rs7745098; P=3.42 × 10(-9), OR=1.21). rs3806624 localizes 5' to the EOMES (eomesodermin) gene within a p53 response element affecting p53 binding. rs7745098 maps intergenic to HBS1L and MYB, a region previously associated with haematopoiesis. These findings provide further insight into the genetic and biological basis of inherited susceptibility to HL.

Variation at 10p12.2 and 10p14 influences risk of childhood B-cell acute lymphoblastic leukemia and phenotype.

Acute lymphoblastic leukemia (ALL) is the major pediatric cancer diagnosed in economically developed countries with B-cell precursor (BCP)-ALL, accounting for approximately 70% of ALL. Recent genome-wide association studies (GWAS) have provided the first unambiguous evidence for common inherited susceptibility to BCP-ALL, identifying susceptibility loci at 7p12.2, 9p21.3, 10q21.2, and 14q11.2. To identify additional BCP-ALL susceptibility loci, we conducted a GWAS and performed a meta-analysis with a published GWAS totaling 1658 cases and 4723 controls, with validation in 1449 cases and 1488 controls. Combined analysis identified novel loci mapping to 10p12.2 (rs10828317, odds ratio [OR] = 1.23; P = 2.30 × 10(-9)) and 10p14 marked by rs3824662 (OR = 1.31; P = 8.62 × 10(-12)). The single nucleotide polymorphism rs10828317 is responsible for the N215S polymorphism in exon 7 of PIP4K2A, and rs3824662 localizes to intron 3 of the transcription factor and putative tumor suppressor gene GATA3. The rs10828317 association was shown to be specifically associated with hyperdiploid ALL, whereas the rs3824662-associated risk was confined to nonhyperdiploid non-TEL-AML1 + ALL. The risk allele of rs3824662 was correlated with older age at diagnosis (P < .001) and significantly worse event-free survivorship (P < .0001). These findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to BCP-ALL and the influence of constitutional genotype on disease development.

Common variation at 3q26.2, 6p21.33, 17p11.2 and 22q13.1 influences multiple myeloma risk.

To identify variants for multiple myeloma risk, we conducted a genome-wide association study with validation in additional series totaling 4,692 individuals with multiple myeloma (cases) and 10,990 controls. We identified four risk loci at 3q26.2 (rs10936599, P = 8.70 × 10(-14)), 6p21.33 (rs2285803, PSORS1C2, P = 9.67 × 10(-11)), 17p11.2 (rs4273077, TNFRSF13B, P = 7.67 × 10(-9)) and 22q13.1 (rs877529, CBX7, P = 7.63 × 10(-16)). These data provide further evidence for genetic susceptibility to this B-cell hematological malignancy, as well as insight into the biological basis of predisposition.

Common variation near CDKN1A, POLD3 and SHROOM2 influences colorectal cancer risk.

We performed a meta-analysis of five genome-wide association studies to identify common variants influencing colorectal cancer (CRC) risk comprising 8,682 cases and 9,649 controls. Replication analysis was performed in case-control sets totaling 21,096 cases and 19,555 controls. We identified three new CRC risk loci at 6p21 (rs1321311, near CDKN1A; P = 1.14 × 10(-10)), 11q13.4 (rs3824999, intronic to POLD3; P = 3.65 × 10(-10)) and Xp22.2 (rs5934683, near SHROOM2; P = 7.30 × 10(-10)) This brings the number of independent loci associated with CRC risk to 20 and provides further insight into the genetic architecture of inherited susceptibility to CRC.

Variation in TP63 is associated with lung adenocarcinoma in the UK population.

BACKGROUND: Variation at TP63 has recently been shown to be associated with lung adenocarcinoma in the Asian population. METHODS: To investigate how this finding translates to the European population we compared the genotypes of SNPs annotating the TP63 locus at 3q28 in 4,462 lung cancer patients, including 911 with adenocarcinoma, and 8,235 controls from the United Kingdom. RESULTS: A statistically significant association between adenocarcinoma risk and SNP genotype was shown: rs10937405, OR = 1.21, P = 1.82 × 10(-4); rs17429138, OR = 1.23, P = 7.49 × 10(-5); and rs4396880, OR = 1.21, P = 2.03 × 10(-4). Haplotype analysis was consistent with a single TP63 risk locus defined by SNPs rs10937405, rs17429138, and rs4396880. While no association between SNPs and small cell lung cancer was shown, the rs10937405 and rs439680 associations were significant for squamous cancer (respective P-values, 0.0022 and 0.02). CONCLUSIONS: These findings show TP63 variation is a risk factor for the development of lung adenocarcinoma in the UK population. Furthermore, they provide additional insight into the subtype-specificity of the 3q28 lung cancer association. IMPACT: Our data confirm the association of 3q28 with lung adenocarcinoma and that this association is not confined to the Asian population. Elucidating the functional basis of this association will be contingent on future fine mapping of the TP63 loci.

Rationale for an international consortium to study inherited genetic susceptibility to childhood acute lymphoblastic leukemia.

Acute lymphoblastic leukemia is the major pediatric cancer in developed countries. To date most association studies of acute lymphoblastic leukemia have been based on the candidate gene approach and have evaluated a restricted number of polymorphisms. Such studies have served to highlight difficulties in conducting statistically and methodologically rigorous investigations into acute lymphoblastic leukemia risk. Recent genome-wide association studies of childhood acute lymphoblastic leukemia have provided robust evidence that common variation at four genetic loci confers a modest increase in risk. The accumulated experience to date and relative lack of success of initial efforts to identify novel acute lymphoblastic leukemia predisposition loci emphasize the need for alternative study designs and methods. The International Childhood Acute Lymphoblastic Leukaemia Genetics Consortium includes 12 research groups in Europe, Asia, the Middle East and the Americas engaged in studying the genetics of acute lymphoblastic leukemia. The initial goal of this consortium is to identify and characterize low-penetrance susceptibility variants for acute lymphoblastic leukemia through association-based analyses. Efforts to develop genome-wide association studies of acute lymphoblastic leukemia, in terms of both sample size and single nucleotide polymorphism coverage, and to increase the number of single nucleotide polymorphisms taken forward to large-scale replication should lead to the identification of additional novel risk variants for acute lymphoblastic leukemia. Ethnic differences in the risk of acute lymphoblastic leukemia are well recognized and thus in assessing the interplay between inherited and non-genetic risk factors, analyses using different population cohorts with different incidence rates are likely to be highly informative. Given that the frequency of many acute lymphoblastic leukemia subgroups is small, identifying differential effects will realistically only be possible through multi-center pooled analyses. Here, we review the rationale for identifying genetic risk variants for acute lymphoblastic leukemia and our proposed strategy for establishing the International Childhood Acute Lymphoblastic Leukaemia Genetics Consortium.