Clay‐Gilmour, Alyssa

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.

Abstract Genome-wide association studies (GWAS) conducted among populations of European ancestry (EA) have identified 23 common single nucleotide polymorphisms (SNPs) associated with multiple myeloma (MM) risk. We hypothesize that the combination of these SNPs in a polygenic risk score (PRS) is likely to be a strong risk factor for MM. However, it is unclear whether the genetic variation associated with MM susceptibility also predisposes to monoclonal gammopathy of undetermined significance (MGUS). Thus, we calculated a PRS and evaluated the association with risk of MM and its precursor, MGUS. We pooled genotype data for 2434 MM and 3446 controls from ten MM GWAS of individuals of EA within the Interlymph Consortium, for 23 MM risk SNPs identified by prior GWAS. An additional 754 MGUS cases were ascertained from Mayo Clinic and MD Anderson clinical practices. To calculate the PRS, we used the risk estimates corresponding to the 23 SNP associations from the largest published MM GWAS. The log of the odds ratio (OR) for each SNP was multiplied by the respective number of risk alleles and summed to generate a PRS for each individual. The PRS was examined continuously, per one standard deviation (SD), and as quintiles, based on the PRS distribution in the controls. Associations of PRS with MM and MGUS risk were examined separately, using multivariable logistic regression assuming an additive model to assess ORs and 95% confidence intervals adjusted for age, sex, and site. We also evaluated age and sex stratified models. The distribution of sex within MM cases, MGUS cases and controls were each ~60% male and ~40% female. The median age was 61, 66, and 66 years for MM cases, MGUS cases and controls, respectively. PRS ranged from 1.52-4.91, with a median PRS of 3.21 for MM cases, 3.19 for MGUS cases, and 3.05 for controls. PRS was significantly associated with MM risk when assessed continuously (OR=1.19 per SD, p=2.2x10-16) and categorically; compared with the middle quintile (Q3), individuals in the highest quintile (Q5) had a 66% increased MM risk (OR=1.66, p=2.3x10-9) and those in the lowest quintile (Q1) had a 38% decreased MM risk (OR=0.62, p=1.3x10-6). PRS was also significantly associated with MGUS risk (OR=1.19 per SD, p=1.7x10-11); individuals with the highest PRS (Q5) had a 77% increased risk (OR=1.77, p=4.0x10-4) and those with lowest PRS (Q1) had 30% decreased risk (OR=0.70, p=0.04), compared with Q3. When stratified by age and sex, similar associations and trends were found. Using an independent sample of MM / MGUS cases and controls, we showed that a PRS constructed from 23 common genetic variants for MM risk is associated with risk of both MM and MGUS, regardless of age or sex. A future direction of this work is testing associations with PRS and clinical characteristics of the MM cases, as well as differences between MGUS cases that progress and those that do not. Our results suggest that common genetic variation may predispose to MGUS as the precursor to MM. Citation Format: Alyssa I. Clay-Gilmour, Michelle A. Hildebrandt, Nicola J. Camp, Elad Ziv, Elizabeth E. Brown, Jonathan N. Hofmann, John J. Spinelli, Graham G. Giles, Parveen Bhatti, Wendy Cozen, Xifeng Wu, Dennis P. Robinson, Aaron D. Norman, Jason P. Sinnwell, Shaji K. Kumar, S Vincent Rajkumar, Susan L. Slager, Celine M. Vachon. Associations between a polygenic risk score and risk of multiple myeloma and its precursor [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2686.

Abstract CLL is a highly heritable cancer with first degree relatives of CLL cases having a 7.5-fold increased CLL risk. Genome-wide association studies (GWAS) and linkage studies have been performed to study inherited predisposition; however a larger proportion of heritability to CLL remains unexplained. Rare coding variants might account for the missing heritability information. Inherited loss of function variants in shelterin complex genes (POT1, ACD, TERF1, TINF2, TERF2, TERF2IP- involved in telomere regulation), CDK1 (critical for cell division) and ATM (tumor suppressor gene) have been found to co-segregate in CLL families and be enriched in CLL cases using exome-wide sequencing data. Our study evaluates rare germline variants from these suspect genes segregating in CLL families who are followed at the Mayo Clinic. Using whole exome sequencing (WES), we sequenced 93 CLL families with at least 2 reported CLL cases consisting of 443 individuals: 160 with CLL, 73 with monoclonal B-cell lymphocytosis (MBL), and 210 relatives. DNA was extracted from buccal cells, coding exons were selectively captured using Agilent 50Mb and SureSelect Human All Exon V4 capture kits; sequencing was performed using Illumina HiSeq 2000. Mayo Clinic's DNASeq pipeline uses Novoalign (initial read alignment), Picard (marking duplicate reads), and the Genome Analysis Toolkit (GATK) for local realignment, recalibration, and variant calling. The variant discovery step leverages GATK's HaplotypeCaller in per sample mode and all of the samples across the cohort are jointly genotyped together. All called variants are evaluated with GATK's Variant Quality Score Recalibration tool and annotated for biological relevance (BIOR). Quality control included removing variants that had &lt;75% call rate across the two capture kits, &lt;8x coverage, or phred score&lt;10, resulting in 317,666 remaining variants. Of these, over 80% of the coding sequence had a median read depth of 23 reads. In our pedigrees, we searched for rare variants within the genes described above. We identified suspect variants with the following criteria: 1) enriched in CLL and MBL samples compared to unaffected samples; 2) multiple affected members with the variant within a family; 3) variants present in all sequenced affecteds within the family; 4) rarely seen in an in-house database of non-cancer controls or 1K Genomes; and 5) predicted to have a functional damaging effect (using SIFT). We identified three novel rare missense variants, defined as functionally deleterious, which each co-segregated within a CLL family. Specifically, these variants from shelterin complex genes; POT1 (rs116916706), TERF2IP (rs138458227), and TERF2 (rs749171225), met the criteria. This study further highlights telomere dysregulation as a key process in CLL development. Investigating rare variants within CLL pedigrees with WES can help identify germline variants impacting predisposition to familial CLL. Citation Format: Alyssa I. Clay-Gilmour, Daniel R. O'Brien, Sara J. Achenbach, Celine M. Vachon, Kari G. Chaffee, Timothy G. Call, Jose F. Leis, Aaron D. Norman, Brian F. Kabat, Sameer A. Parikh, Neil E. Kay, Esteban Braggio, James R. Cerhan, Susan L. Slager. Rare germline variants segregating in chronic lymphocytic leukemia (CLL) families [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1226.

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.

Abstract Background: The precise mechanism by which the immune system is adversely affected in cancer patients remains poorly understood, but the accumulation of immunosuppressive/protumorigenic myeloid-derived suppressor cells (MDSCs) is thought to be a prominent mechanism contributing to immunologic tolerance of malignant cells in epithelial ovarian cancer (EOC). To this end, we hypothesized genetic variation in MDSC pathway genes would be associated with survival after EOC diagnoses. Methods: We measured the hazard of death due to EOC within 10 years of diagnosis, overall and by invasive subtype, attributable to SNPs in 24 genes relevant in the MDSC pathway in 10,751 women diagnosed with invasive EOC. Versatile Gene-based Association Study and the admixture likelihood method were used to test gene and pathway associations with survival. Results: We did not identify individual SNPs that were significantly associated with survival after correction for multiple testing (P &lt; 3.5 × 10−5), nor did we identify significant associations between the MDSC pathway overall, or the 24 individual genes and EOC survival. Conclusions: In this well-powered analysis, we observed no evidence that inherited variations in MDSC-associated SNPs, individual genes, or the collective genetic pathway contributed to EOC survival outcomes. Impact: Common inherited variation in genes relevant to MDSCs was not associated with survival in women diagnosed with invasive EOC. Cancer Epidemiol Biomarkers Prev; 26(3); 420–4. ©2016 AACR.

Abstract Introduction: The use of multi-drug regimens, including taxanes (T), has increased survival for breast cancer patients. However, T-induced neuropathy, one of the most common side effects, cannot be predicted, prevented or effectively treated. Using genome wide analyses (GWA), we sought to identify common genetic variants that infer susceptibility to clinical and self-reported neuropathy. Methods: Women with high risk breast cancer were enrolled in SWOG 0221, a trial of dosing schedules for AC and T. Genotyping was performed using the Illumina 1M chip and SNPs examined in relation to CTCAE grade 3 and 4 neurotoxicity and self-reported symptoms (FACT-TAXANE) among women providing blood samples. Significant associations were tested for replication in CALGB40101. Results: Among 1269 European Americans (EA) and 139 African Americans (AA) in S0221, there were 147 and 31 grade 3+ neuropathies, respectively. In CALGB40101, grade 3 toxicity was less frequent (59/855 EA and 15/117 AA). In S0221 EA patients, loci on chromosomes 7, 10, 16 and 17 showed associations with grades 3 and 4 neuropathy at p&lt;5×10e-06, with the most significant genetic association on chromosome 17, rs910920, p = 6.8×10e-08, odds ratio (OR) = 0.44 (95% confidence interval (CI), 0.32-0.61); this SNP is a cis eQTL, thought likely to affect transcription factor binding and linked to expression of VPS53 in lymphoblastoid cell lines. While the rs910920 association was specific to EA, the A allele in rs1858826 on chromosome 7 was associated with reduced odds of neuropathy in both EA (p = 7×10e-07) and AA (p = .02) groups. The S0221 associations at p&lt;5×10e-06 did not replicate in CALGB40101. A subset of patients (n = 677 EA and n = 42 AA) completed the FACT-TAXANE. The most significant associations with this neuropathy phenotype (p&lt;5×10e-07) were on chromosomes 2 (rs4267555 and rs4443044) and 3 (rs2421468). Of the 168 SNPs tested for association with neuropathy in previous studies, rs1695 in GSTP1 (p = 0.0009) and rs228591 in ATM (p = 0.0006), showed some evidence of association with FACT-TAXANE scores; the latter SNP was significant at p = .005 in ordinal regression analysis of neuropathy grades 0-4 in CALGB40101. None of these SNPs were significantly associated with grade 3 CTCAE neuropathy at p&lt;.05, nor were the top CTCAE associations in relation to FACT-TAXANE scores. Replication of the most significant S0221 FACT-TAXANE findings in CALGB40101 is ongoing. Conclusions: Failure to replicate S0221 in CALGB40101, and the lack of coherence between self-report and clinical grade in S0221 could be attributable to differences in trial designs as well as differences in defining a complex phenotype. Larger homogeneous trials with similar dosing and criteria for defining neuropathy are needed to assess the relationship of genetics and the neuropathy symptom spectrum. Citation Format: Lara Sucheston-Campbell, Alyssa Clay, William E. Barlow, G Thomas Budd, Dan Stram, Chris Haiman, Li Yan, Gary Zirpoli, Song Yao, Dawn Hershman, Kathy S. Albain, Daniel F. Hayes, Halle Moore, Timothy J. Hobday, James A. Stewart, Claudine Isaacs, Muhammad Salim, Julie R. Gralow, Gabriel N. Hortobagyi, Robert B. Livingston, Xin Sheng, Deanna L. Kroetz, Christine B. Ambrosone. A genome-wide association study identifies novel loci associated with taxane-related sensory neuropathy in breast cancer patients enrolled in a cooperative group clinical trial (SWOG S0221). [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 5489. doi:10.1158/1538-7445.AM2015-5489