Yan, Li

Abstract Genome-wide association studies (GWAS) of esophageal adenocarcinoma (EAC) and its precursor, Barrett’s esophagus (BE), have uncovered significant genetic components of risk, but most heritability remains unexplained. Targeted assessment of genetic variation in biologically relevant pathways using novel analytical approaches may identify missed susceptibility signals. Central obesity, a key BE/EAC risk factor, is linked to systemic inflammation, altered hormonal signaling and insulin-like growth factor (IGF) axis dysfunction. Here, we assessed IGF-related genetic variation and risk of BE and EAC. Principal component analysis was employed to evaluate pathway-level and gene-level associations with BE/EAC, using genotypes for 270 single-nucleotide polymorphisms (SNPs) in or near 12 IGF-related genes, ascertained from 3295 BE cases, 2515 EAC cases and 3207 controls in the Barrett’s and Esophageal Adenocarcinoma Consortium (BEACON) GWAS. Gene-level signals were assessed using Multi-marker Analysis of GenoMic Annotation (MAGMA) and SNP summary statistics from BEACON and an expanded GWAS meta-analysis (6167 BE cases, 4112 EAC cases, 17 159 controls). Global variation in the IGF pathway was associated with risk of BE (P = 0.0015). Gene-level associations with BE were observed for GHR (growth hormone receptor; P = 0.00046, false discovery rate q = 0.0056) and IGF1R (IGF1 receptor; P = 0.0090, q = 0.0542). These gene-level signals remained significant at q < 0.1 when assessed using data from the largest available BE/EAC GWAS meta-analysis. No significant associations were observed for EAC. This study represents the most comprehensive evaluation to date of inherited genetic variation in the IGF pathway and BE/EAC risk, providing novel evidence that variation in two genes encoding cell-surface receptors, GHR and IGF1R, may influence risk of BE.

Abstract The first two authors and last two authors contributed equally. Genome-wide association studies (GWAS) have identified risk loci for Acute Lymphoblastic Leukemia (ALL), Chronic Lymphoblastic Leukemia (CLL) and Non-Hodgkin Lymphoma, however an Acute Myeloid Leukemia (AML) GWAS has not been published to date. We performed a GWAS to identify AML and Myelodysplastic Syndrome (MDS) risk loci using a nested case-control study design in the DISCOVeRY-BMT cohorts which includes almost 2000 AML and MDS patients as cases and 2813 unrelated donors as controls. Genotyping was performed using the Illumina Human OmniExpress BeadChip and imputed using the Haplotype Reference Consortium, yielding > 8 million high-quality variants for statistical analysis. Logistic regression models with AML (de novo AML with normal cytogenetics, de novo AML with abnormal cytogenetics and therapy-related AML) and MDS cases and European American healthy donor controls adjusted for age and sex were used to test the association of each SNP with disease status. To identify the strongest association signal with disease we conducted a summary statistic SNP-based association analysis (ASSET) using non-overlapping AML and MDS cases implemented in R statistical software. ASSET uses an exhaustive search for SNPs with small but common pleiotropic effects across groups of traits while accounting for the multiple tests required by the subset search, as well as any shared controls between groups. This approach allowed us to further investigate the heterogeneity within AML subtypes and to gain increased power by pooling subtypes that show pleiotropic effects. ASSET genome wide (GW) significance is defined as P<5.0x10-8, however, we present results of each subset analysis. ASSET analyses identified an association of T allele at rs12203592, Interferon Regulatory Factor 4 (IRF4), with increased risk of de novo AML and MDS (Figure 1) (OR=1.37; 95% CI, 1.25-1.51, Pmeta=1.59x10-10). The variant has been reported significant in GWAS of skin pigmentation, hair color, non-melanoma skin cancer, squamous cell carcinoma, actinic keratosis, and childhood ALL; GW case-control studies of neuroblastoma and breast cancer show suggestive association signals (P<5 x 10-6). IRF4 belongs to the IRF family of transcription factors and is a key regulator of differentiation stages in hematopoiesis. The T allele at rs12203592 is associated with significantly increased expression of IRF4 in whole blood and lung tissue and in vitro mouse studies have shown the T allele to be associated with higher levels of IRF4 transcription. rs12203592 is <80bp from an IRF4 transcription start site and in an important position within NF-κB motifs in multiple blood and immune cell lines, supporting the hypothesis that this SNP modulates NF-κB repression of IRF4 expression. Another variant in the IRF4 regulatory region, rs62389423 (A allele), showed a putative association with subsets of de novo AML and MDS (OR=1.36; 95% CI,1.21-1.52, Pmeta=1.2x10-7). Although not in linkage disequilibrium with rs12203592, the A allele has also been previously associated with susceptibility to skin cancer and melanoma in multiple GWAS. Analysis of 596 UKbiobank CLL cases and >300,000 controls, shows the A allele correlates with a 50% increased risk of CLL (95% CI, 1.32-1.73, P= 2.4x10-9).Several other GWAS also show an association between additional IRF4 variants and CLL. Analyses in individual subtypes also revealed an intronic variant, rs10098598, in RAD21 to be associated with de novo AML (OR=2.18; 95% CI, 1.60-2.97; Pmeta=8.9x10-7). RAD21, on8q24, a region known to be associated with multiple cancer types, is involved in DNA double-stranded break repair and in chromatid cohesion in mitotic cells. Loss of function would theoretically lead to chromosomal instability and tumorigenesis. Components of the cohesion complex (including RAD21) are somatically mutated in ~12% of de novo AML and MDS patients, where it is frequently an early event. We provide the first GW evidence of association between a common variant and AML susceptibility. This SNP has been shown to be associated with multiple phenotypes, suggesting that there are pleiotropic effects at work. Our RAD21 finding is consistent with the role of cohesion in leukemogenesis and provides some evidence of its role in de novo AML susceptibility. Replication and further subset analyses with genome-wide data in >2000 AML cases are ongoing. Disclosures Griffiths: Alexion Inc.: Honoraria, Research Funding; Astex/Otsuka Pharmaceuticals: Honoraria, Research Funding; Celgene, Inc: Honoraria, Research Funding; Novartis, Inc.: Research Funding; Pfizer, Inc.: Research Funding. McCarthy:Bristol Myers Squibb: Consultancy, Honoraria; Celgene: Consultancy, Honoraria, Research Funding; Karyopharm: Consultancy, Honoraria; Takeda: Consultancy, Honoraria; Janssen: Consultancy, Honoraria.

Abstract Both genome wide association studies (GWAS) of common variation and exome wide association studies (EXWAS) of rare variation have successfully identified disease susceptibility variants for a variety of diseases. One GWAS of inherited susceptibility to Acute Myeloid Leukemia (AML) has been conducted, but no EXWAS have been performed to measure risk of AML attributable to low-frequency constitutional genetic variation. We performed the first EXWAS of risk of AML as a nested case-control study in the DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to one-Year mortality after BMT) cohorts. The DISCOVeRY-BMT parent study examined transplant-related mortality in leukemia patients undergoing unrelated donor allogeneic BMT. To identify low frequency variants and genes contributing to increased susceptibility to AML we used genotype data from the Illumina HumanExome BeadChip typed in the DISCOVeRY-BMT cohorts; the HumanExome BeadChip contains 242,901 variants, which are mainly protein-coding variants. The optimal sequence kernel association test (SKAT-O) was used to analyze gene-level associations with risk of AML. These gene-based tests evaluate the cumulative effects of multiple single gene variants on risk of AML. Analyses were performed in all European American AML cases and two subtypes: 1) de novo AML, 2) de novo AML with normal cytogenetics. Models were adjusted for age at transplant and principal components to control for population stratification. For gene-based tests at least 2 variants with minor allele frequency (MAF) ≤ 5%, were required to be present in the gene. This yielded a total of 13,687 genes tested, and a Bonferroni corrected significance level of P<3.65 x 10-6. Association tests were performed in 1,189 AML cases reported to CIBMTR 2000-08 (Cohort 1) and 327 AML cases reported to CIBMTR from 2009-11 (Cohort 2). Controls in Cohorts 1 (n=1,986) and 2 (n= 515) were 10/10 HLA-matched unrelated donors who passed a comprehensive medical exam and deemed healthy. We used metaSKAT to combine Cohorts 1 and 2 and obtain p-values of association with AML. We present the results of gene-level tests significant in both cohorts. The likely pathogenicity of these variants was determined in silico using SIFT, PolyPhen and MutationTaster. Patient characteristics are in Table 1. DNMT3A, on chromosome 2, was associated in the gene-based test with risk of AML (Pmeta=1.70x10-9, Table 2). Three missense variants at MAF <1% comprise both overall AML and de novo AML gene-based association: exm177559 (Asn->Ser), exm177507 (Arg->His), and exm177543 (Arg->Trp). Normal cytogenetics de novo AML gene-based assocations consisted of only 2 of these variants: exm177559 and exm177507 (Table 2). While prevalence of exm177507 is <1% for all AML cases, in de novo AML with normal cytogenetics the MAF was higher at 3%. The other 2 variants had a MAF<1% irrespective of subtype. Somatically, DNMT3A is most frequently mutated in hematologic malignancies, with >30% of de novo AML cases with a normal karyotype and >10% of MDS patients having DNMT3A mutations. Although these are germline gene associations all three of the variants found have been reported somatically in hematologic malignancies. In 200 AML cases from The Cancer Genome Atlas (TCGA) p.R882H (represented as exm177507 on the exome chip) was a frequent somatic mutation (25%). Exm177543 (p.R635W) and exm177559 (p.N501S) are reported in the Catalogue of Somatic Mutations in Cancer (COSMIC) as somatic mutations involved in hematopoietic and lymphoid tissue in both cell lines and humans. Exm177507 and exm177543 show evidence of pathogenicity in all three in silico tools, while exm177559 was reported as deleterious and disease causing by Sift and MutationTaster, respectively. Our results show that multiple potentially pathogenic missense germline variants in DNMT3A comprise the gene-based association with AML, specifically de novo AML with normal cytogenetics. Given the functional nature of these variants it is possible germline risk stratification could be informative in determining AML risk, and subsequently development of AML harboring DNMT3A mutations. Confirmation of these findings in additional cohorts could have implications for individualized risk screening, prediction and prognosis. Additional cytogenetic subgroup analyses, including treatment-related AML, are underway. Disclosures Hahn: Novartis: Equity Ownership; NIH: Research Funding. McCarthy:Sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; Onyx: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bristol Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees; The Binding Site: Honoraria, Membership on an entity's Board of Directors or advisory committees; Karyopharm: Honoraria, Membership on an entity's Board of Directors or advisory committees; Gamida Cell: Honoraria, Membership on an entity's Board of Directors or advisory committees. Sucheston-Campbell:NIH/NCI: Research Funding.

Abstract Over the past decade, there have been numerous published candidate gene association studies of non-HLA variants, genes and/or pathways in relation to BMT survival outcomes. Most of these studies investigated a small number of candidate single nucleotide polymorphisms (SNPs), chosen based on known gene function, in a limited number of cases without replication. Here we test the association of previously reported SNPs with survival outcomes (overall survival (OS), progression-free survival (PFS), death due to disease (DD) or transplant-related mortality (TRM)) after URD BMT using data from the DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to one-Year mortality after BMT) study which included 2,052 European descent recipient-unrelated donor pairs (Cohort 1) and 763 recipient-unrelated donor pairs (Cohort 1) reported to the Center for International Blood and Marrow Transplant Research (CIBMTR) between 2000-2009 and 2009-2011, respectively. DISCOVeRY-BMT is a Genome Wide Association Study (GWAS) that includes over 9 million typed or imputed SNPs. This represents the largest study to date and is well powered to test both SNP and gene-based association hypotheses. We conducted an extensive literature review of PubMed using MESH terms related to transplant and hematologic disease which yielded 76 candidate gene association studies that investigated a total of 46 genes with any survival outcome. Of the 46 genes, 11 genes (CTLA4, CYP2B6, CYP2C19, GSTM1, GSTT1, IL23R, IL6, MTHFR, NOD2/CARD15, TNFSRF1B, and VDR) were studied in at least 2 publications. A total of 16 SNPs in 8 of these genes reported a significant (P< 0.05) association with at least one survival outcome (OS, PFS, DD or TRM). We tested these 16 SNPs for an association with the survival outcome originally reported in each article (Table 1). All association models included age at BMT, diagnosis (AML, ALL, MDS), disease status at BMT, cell source (blood, bone marrow) and year of BMT. P-values from Cohorts 1 and 2 were combined using a fixed-effects model in METAL software. Since the candidate gene studies have a gene-based hypothesis, we also performed gene-level testing using the VErsatile Gene-Based Association Study (VEGAS2) software to measure the aggregate effect of all available typed and imputed SNPs from our GWAS data in these 8 genes. The candidate gene association studies included heterogeneous patient populations of autologous and/or allogeneic (related or unrelated donors) BMT for a variety of hematologic diseases. Most candidate gene studies performed genotyping on both donors and recipients, and SNPs were tested either independently for donor or recipient genomes, or considering the total number of variant alleles the recipient-donor pair had at a given locus (0-4) (Table 1). The 16 SNPs from the previous candidate gene association studies were not significantly associated at a nominal P<0.05 to either the survival outcome reported in the original published paper, or to any survival outcome (Table 1). Likewise, the VEGAS2 gene-based analysis did not reveal any statistically significant associations for any of the 8 genes at P<0.05. Our study also corroborated prior reports of SNPs and genes that were not associated with survival outcomes. Candidate gene association studies use a priori biological knowledge to select genes that are likely associated with outcome and subsequently SNPs within the gene thought likely to affect function based on knowledge at the time of the investigation. Our findings suggest that neither the genetic variation tested in these candidate genes, nor other individual SNPs in the gene independently or collectively are not associated with survival outcomes in recipients of URD BMT. Our failure to replicate the prior candidate gene studies may be due to these SNPs being dependent on the underlying disease (ie. AML, ALL, etc), the donor relation (related, unrelated, autologous), or that there is no true association of these SNPs with the outcomes studied. #These authors contributed equally Table 1. Table 1. Disclosures McCarthy: Gamida Cell: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Karyopharm: Honoraria, Membership on an entity's Board of Directors or advisory committees; Onyx: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bristol Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; The Binding Site: Honoraria, Membership on an entity's Board of Directors or advisory committees; Sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees. Hahn:Novartis: Equity Ownership; NIH: Research Funding. Sucheston-Campbell:NIH/NCI: Research Funding.

Key PointsIKZF1 associations with high-risk B-ALL may differ by age and sex. A novel variant on chromosome 14, rs189434316, is associated with over a 3.5-fold risk of normal cytogenetic B-ALL.

Key Points We conducted the first exome-wide association study between germ line variant genotype and survival outcomes after unrelated-donor BMT. A number of novel genes were found to significantly affect survival outcomes.

Key PointsCandidate SNP associations with survival outcomes after URD transplant are most likely false-positive findings. Over 85% of candidate SNPs are not linked to a biochemical function; of those that are, about half are not linked to the candidate gene.