Cozen, Wendy

Abstract Background: Persons of African ancestry (AA) experience a 1.5-2-fold risk of multiple myeloma (MM) compared to persons of European ancestry (EA). We assembled a set of MM patients with self-reported AA in order to evaluate the contribution of genetics to etiology in this high-risk group. Methods: Here we present the results of a meta-analysis of two GWAS in 1,813 cases and 8,871 controls of AA. We also conducted an admixture mapping scan to identify risk alleles associated with local ancestry, fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA, and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. Finally, we conducted an eQTL analysis measuring gene expression in those genes harboring a risk variant in malignant plasma cells from 292 of the patients from a single site. Results: In GWAS analysis, we identified two suggestive novel loci located at 9p24.3 and 9p13.1 at P<1 × 10-6, but no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. 20 of the 23 known EA risk variants showed directional consistency and 9 replicated at P<0.05 in AA individuals. In eight regions, we identified markers that better capture MM risk in persons of AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95%CI: 1.56, 2.11) increased MM risk compared to those with average risk (25-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P= 5.1 × 10–12). Conclusion: Our study shows that common genetic variation contributes to MM risk individuals of AA. This abstract is also being presented as Poster C040. Citation Format: Zhaohui Du, Niels Weinhold, Gregory Chi Song, Kristen A. Rand, David J. Van Den Berg, Amie E. Hwang, Xin Sheng, Victor Hom, Sikander Ailawadhi, Ajay K. Nooka, Seema Singhal, Karen Pawlish, Edward Peters, Cathryn Bock, Ann Mohrbacher, Alexander Stram, Sonja I. Berndt, William J. Blot, Graham Casey, Victoria L. Stevens, Rick Kittles, Phyllis J. Goodman, W. Ryan Diver, Anselm Hennis, Barbara Nemesure, Eric A. Klein, Benjamin A. Rybicki, Janet L. Stanford, John S. Witte, Lisa Signorello, Esther M. John, Leslie Bernstein, Antoinette Stroup, Owen W. Stephens, Maurizio Zangari, Frits Van Rhee, Andrew Olshan, Wei Zheng, Jennifer J. Hu, Regina Ziegler, Sarah J. Nyante, Sue Ann Ingles, Michael Press, John David Carpten, Stephen Chanock, Jayesh Mehta, Graham A Colditz, Jeffrey Wolf, Thomas G. Martin, Michael Tomasson, Mark A. Fiala, Howard Terebelo, Nalini Janakiraman, Laurence Kolonel, Kenneth C. Anderson, Loic Le Marchand, Daniel Auclair, Brian C.-H. Chiu, Elad Ziv, Daniel Stram, Ravi Vij, Leon Bernal-Mizrachi, Gareth J. Morgan, Jeffrey A. Zonder, Carol Ann Huff, Sagar Lonial, Robert Z. Orlowski, David V. Conti, Christopher A. Haiman, Wendy Cozen. A meta-analysis of genome-wide association study and eQTL analysis of multiple myeloma among African Americans [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr PR05.

Abstract So far, 23 germline susceptibility loci have been associated with multiple myeloma (MM) risk. It is unclear whether the genetic variation associated with MM susceptibility also predisposes to its precursor, monoclonal gammopathy of undetermined significance (MGUS). Leveraging 2434 MM cases, 754 MGUS cases, and 2 independent sets of controls (2567/879), we investigated potential shared genetic susceptibility of MM and MGUS by (1) performing MM and MGUS genome-wide association studies (GWAS); (2) validating the association of a polygenic risk score (PRS) based on 23 established MM loci (MM-PRS) with risk of MM, and for the first time with MGUS; and (3) examining genetic correlation of MM and MGUS. Heritability and genetic estimates yielded 17% (standard error [SE] ±0.04) and 15% (SE ±0.11) for MM and MGUS risk, respectively, and a 55% (SE ±0.30) genetic correlation. The MM-PRS was associated with risk of MM when assessed continuously (odds ratio [OR], 1.17 per SD; 95% confidence interval [CI], 1.13-1.21) or categorically (OR, 1.70; 95% CI, 1.38-2.09 for highest; OR, 0.71; 95% CI, 0.55-0.90 for lowest compared with middle quintile). The MM-PRS was similarly associated with MGUS (OR, 1.19 per SD; 95% CI, 1.14-1.26 as a continuous measure, OR, 1.77, 95%CI: 1.29-2.43 for highest and OR, 0.70, 95%CI: 0.50-0.98 for lowest compared with middle quintile). MM and MGUS associations did not differ by age, sex, or MM immunoglobulin isotype. We validated a 23-SNP MM-PRS in an independent series of MM cases and provide evidence for its association with MGUS. Our results suggest shared common genetic susceptibility to MM and MGUS.

Purpose: Multiple myeloma (MM) treatment has changed tremendously, with significant improvement in patient out-comes. One group with a suboptimal benefit is patients with high-risk cytogenetics, as tested by conventional karyotyping or fluorescence in situ hybridization (FISH). Methodology for these tests has been published, but not necessarily standardized. Methods: We address variability in the testing and reporting methodology for MM cytogenetics in the United States using the ongoing African American Multiple Myeloma Study (AAMMS). We evaluated clinical and cytogenetic data from 1,221 patients (1,161 with conventional karyotyping and 976 with FISH) tested between 1998 and 2016 across 58 laboratories nationwide. Results: Interlab and intralab variability was noted for the number of cells analyzed for karyotyping, with a significantly higher number of cells analyzed in patients in whom cytogenetics were normal (P 5.0025). For FISH testing, CD138-positive cell enrichment was used in 29.7% of patients and no enrichment in 50% of patients, whereas the remainder had unknown status. A significantly smaller number of cells was analyzed for patients in which CD138 cell enrichment was used compared with those without such enrichment (median, 50 v 200; P, .0001). A median of 7 loci probes (range, 1-16) were used for FISH testing across all laboratories, with variability in the loci probed even within a given laboratory. Chromosome 13–related abnormalities were the most frequently tested abnormality (n5956; 97.9%), and t(14;16) was the least frequently tested abnormality (n 5 119; 12.2%). Conclusions: We report significant variability in cytogenetic testing across the United States for MM, potentially leading to variability in risk stratification, with possible clinical implications and personalized treatment approaches.

Abstract Background. The strongest established prognostic factors for follicular lymphoma are age and clinical factors. Biologic studies suggest a major role for the host immunologic environment in follicular lymphomagenesis, which is in part determined by host genetic background. Cytokines are key regulators of immune function and regulation, are highly polymorphic, and have been implicated in lymphoma etiology and prognosis. We evaluated the hypothesis that inherited variability in cytokine and related immune genes impact overall survival in follicular lymphoma. Methods. We genotyped 73 single nucleotide polymorphisms (SNPs) from 44 candidate cytokine and immune genes in 278 follicular lymphoma patients who participated in a population-based case-control study conducted from 1998–2000 in the Detroit, Iowa, Los Angeles and Seattle SEER (Surveillance, Epidemiology and End Results) cancer registries. Baseline clinical data and survival through early 2005 were obtained from cancer registry files. Cox proportional hazards analysis was used to estimate hazard ratios (HR) and corresponding 95% confidence intervals for the association between individual SNPs and overall survival, adjusting for age, demographic and clinical factors. We also used parallel modeling strategies to identify the best summary multi-SNP risk score to predict survival, and applied these results in a time-dependent receiver-operator characteristics (ROC) analysis. Results. The median age at diagnosis was 57 years (range, 25–74), and 59 (21%) of the patients died during follow-up, with a median follow-up of 59 months (range, 27 – 78 months) for surviving patients. In multivariate modeling, SNPs in IL8 (rs4073; HRTT=2.14, 1.26–3.63), IL2 (rs2069762; HRGT/TT=1.80, 1.06–3.05), IL12B (rs3212227; HRAA/CC=1.83, 1.06–3.06), and IL1RN (rs454078; HRAA=1.93, 1.11–3.34) were the strongest and most robust predictors of survival. A summary score of the number of deleterious genotypes from these four genes was strongly associated with survival (p=1.4 × 10−5) after accounting for demographic and clinical variables (HR=2.06 per deleterious genotype, 1.52–2.79). A combination of the four SNPs and the demographic and clinical risk scores was strongly associated with survival (p=1.8 × 10−11); the 5-year Kaplan Meier survival estimates were 96% (93%–100%), 72% (62%–83%) and 58% (48%–72%) for low, intermediate, and high risk groups respectively. In a time-dependent ROC analysis the three risk groups had an area under the curve (AUC) of 0.83 at 72 months of follow-up (95% CI 0.72–0.93). Conclusion: Host genetic variability in immune genes, particularly IL8, IL2, IL12B, and IL1RN, individually and particularly in combination, was associated with overall survival in follicular lymphoma after accounting for demographic and clinical factors. Host immunogenetics is a promising class of prognostic factors in follicular lymphoma that warrants further evaluation.

Abstract Subtle differences in immune response may play a role in non-Hodgkin lymphoma (NHL) etiology. Because adult immune response may be influenced by early childhood exposures, we examined the role of childhood crowding, history of atopic disease, and other childhood immune-related exposures on the risk of non-Hodgkin lymphoma in a multi-center case-control study. Interviews were completed with 1,321 cases ascertained from population-based cancer registries in Seattle, Detroit, Los Angeles and Iowa, and with 1,057 frequency-matched controls, selected by random-digit dialing and from the Health Care Financing Administration (HCFA) database. The association between NHL risk in relation to atopy and other exposures was assessed using multivariable logistic regression methods. Most types of allergy were associated with protection from NHL, with hay fever especially protective against all NHL combined (Odds Ratio [OR] = 0.71, 95% confidence interval [CI]= 0.54–0.94), diffuse large B-cell lymphoma [DLBCL] (OR=0.61, 95% CI=0.41–0.91), and follicular lymphoma (OR=0.70, 95% CI=0.45–1.09). A history of eczema increased risk of follicular lymphoma (OR=1.92, 95% CI= 1.08–3.41) but not DLBCL (OR=1.06, 95% CI= 0.55.2.04). Asthma in childhood was not associated with risk of NHL. Risk of DLBCL (OR =1.72, 95% CI=1.17–2.52), but not follicular lymphoma (OR=1.15, 95% CI=0.75–1.76) was elevated for the youngest compared to the oldest of siblings. Neither number of siblings nor years between births of siblings were significantly associated with risk. These results suggest that some immune-related exposures may affect NHL risk.

Abstract Background: Several studies suggest that smoking and obesity may increase the risk of developing NHL, but the impact of these factors on survival is relatively unexplored. One recent population-based Italian study found that smoking may negatively impact overall survival. Methods: We evaluated the association of body mass index (BMI) and cigarette smoking on overall survival in 1,286 patients who participated in a population-based case-control study conducted through the Surveillance, Epidemiology, and End Results (SEER) cancer registries in Detroit, Iowa, Seattle and Los Angeles. BMI (on 1189 participants) and a smoking history (on 471 participants) were obtained through in-person interviews. Histology, stage, presence of B symptoms, first course of therapy, date of last follow-up, and vital status were derived from linkage to registry databases at each study site in the spring of 2005. Hazard Ratios (HRs) and 95% confidence intervals were estimated using Cox proportional hazards models, while adjusting for demographic and clinical factors. Results: The median age at diagnosis was 58 years (range: 16 – 74), and the most common histologies were DLBCL (33%), follicular (25%), SLL/CLL (12%), marginal zone (9%), mantle cell (4%) and Peripheral T-cell (3%). The median BMI was 26.6 (range: 16.2 – 47.3), and 26% of patients were classified as obese (i.e. BMI >30). Of patients with a smoking history, 47% had never smoked, 19% were active smokers at diagnosis, and 34% were former smokers. No significant differences were observed between those with and without tobacco use data in terms of primary lymphoma characteristics. Overall survival (OS) was defined as time from diagnosis to death due to any cause. Median follow-up on all living patients was 58 months (range: 4 – 78), and 333 patients had died. When BMI was assessed as a continuous measure, it was borderline significantly associated with OS (p=0.050). An indicator of whether or not BMI > 30 was also significantly associated with OS (HR=1.30; 95% CI: 1.01 – 1.68; p=0.042). Ever vs. never smokers also exhibited a survival disadvantage (HR=1.66; 95% CI: 1.12 – 2.44; p=0.011). When broken down further, former smokers who had smoked for >10 years (HR=1.61; 95% CI: 1.02 – 2.52; p=0.039) and current smokers (HR=1.79; 95% CI: 1.07 – 3.02; p=0.028) were at a survival disadvantage over those who had never smoked. When both smoking and BMI status were included in the model, only the indicator of being a current smoker remained statistically significant (p=0.031), with a trend toward decreased survival for former smokers who had smoked >10 years (p=0.061) and those with BMI >30 (p=0.058). Conclusions: NHL patients who smoke or are obese had a poorer overall survival. These observations warrant further investigation, including whether interventions for smoking cessation or weight loss can impact NHL survival.

Abstract Background: Mantle cell lymphoma (MCL) is an incurable non-Hodgkin lymphoma. MCL is characterized by the chromosomal translocation t(11;14) leading to overexpression of cyclin D1 (CCND1), a regulator of the cell cycle. NF-κB plays a central role in key cellular processes, including cell cycle regulation and apoptosis, and it has been shown to be constitutively activated in MCL cell lines and in patient biopsy samples and to play a key role in the growth and survival of MCL cells. In an exploratory study, we evaluated the hypothesis that germline variation in candidate genes from the cell cycle and NF-κB pathways predict overall survival in MCL. Methods: We genotyped 235 SNPs from 29 cell cycle genes and 447 SNPs from 55 NF-κB genes in 39 MCL patients aged 38–64 years who participated in a population-based study conducted from 1998–2000 through the Surveillance, Epidemiology, and End Results (SEER) cancer registries in Detroit, Seattle, Iowa, and Los Angeles county. Tagging single nucleotide polymorphisms (SNPs) were selected from HapMap. Stage, presence of B-symptoms, first course of therapy, date of last follow-up, and vital status were derived from linkage to registry databases at each study site. To assess the statistical significance of a gene, we first used principal components analysis to capture the genetic variation of the SNPs within a gene. Cox proportional hazards analysis was then used to assess the association between the top principal components from a gene with overall survival. Results: Through early 2005, there were 22 deaths in 39 patients (56%). The median follow-up of the 17 surviving patients was 56 months (range 36–70). Eleven of the 55 genes (TNFSF14 p=0.0007, NF-κBIA p=0.003, IL1R2 p=0.01, TNFRSF25 p=0.02, TNFSF13B p=0.02, FAS p=0.03, TLR2 p=0.04, TNFSF10D p=0.04, TNFSF10 p=0.06, MAP3K2 p=0.06) in the NF-κB pathway were associated with overall survival. In the cell cycle pathway, only 3 of the 29 genes (CASP5 p=0.06, CASP9 p=0.09, BCL2 p=0.09) were associated with overall survival, and a previously reported SNP in CCND1 (rs603965) showed no association (p=0.5). Permutation p-values for observing at least this number of associated tests in each pathway due to chance (from multiple testing) was p=0.8 for the cell cycle pathway and p=0.03 for NF-κB pathway. The small sample size precluded multi-gene modeling. Conclusion: In a population series of MCL patients, we found suggestive evidence that host genetic variation in candidate NF-κB pathway genes affects overall survival. These results complement tumor expression data implicating this pathway in the pathogenesis of MCL. In addition, such findings further support the value of developing therapeutic targets from this pathway.

Abstract Background. Intracellular one-carbon transfer reactions are essential for nucleotide synthesis and methylation of biologic compounds including DNA. Previous studies have linked genetic variants in one-carbon metabolism genes with risk of developing NHL, but little is known regarding the impact of these variants on disease outcome. We evaluated the hypothesis that inherited genetic variation in genes involved in one-carbon metabolism is associated with overall survival in DLBCL. Methods. We genotyped 30 single nucleotide polymorphisms (SNPs) from 18 candidate one-carbon metabolism genes in 215 DLBCL patients who participated in a population-based case-control study conducted from 1998–2000 using the SEER (Surveillance, Epidemiology and End Results) cancer registries in the Detroit, Iowa, Los Angeles and Seattle. Stage, B-symptoms, first course of therapy, date of last follow-up and vital status through early 2005 were obtained from cancer registry files. Cox proportional hazards analysis was used to estimate hazard ratios (HR) and corresponding 95% confidence intervals for the association between individual SNPs and overall survival, adjusting for age, demographic and clinical factors. We also used parallel modeling strategies to identify the best summary multi-SNP risk score to predict survival. Results. The median age at diagnosis was 57 years (range, 20–74), and 50 (23%) of the patients died during follow-up, with a median follow-up of 57 months (range, 31–78 months) for surviving patients. After adjusting for demographic and clinical variables, SNPs in SHMT1 (rs1979276; HRCT/TT=2.47, 1.31–4.67), BHMT (rs585800; HRAT/TT=2.02, 1.16–3.54), and TCN1 (rs526934; HRTT=1.86, 1.04–3.33) were the strongest and most robust predictors of survival. A summary score of the number of deleterious genotypes (0–3) from these three genes was strongly associated with survival (p=7.9 x 10−5) after accounting for demographic and clinical variables (HR=2.58 per deleterious genotype, 95% CI 1.75–3.80). A risk score combining the three SNPs with clinical and demographic variables (score of 0 to 5) was even more strongly associated with survival (p=1.4 x 10−13); the Kaplan Meier survival curves are shown in the Figure. In a time-dependent ROC analysis, the combined risk score had a concordance index of 0.75 at 5 years of follow-up (95% CI 0.69–0.81). Conclusion: Host genetic variation in the one-carbon metabolism genes SHMT1, BHMT, and TCN1, individually and particularly in combination, was associated with overall DLBCL survival after accounting for clinical and demographic factors, supporting a role for this pathway in disease progression. Future work should evaluate interactions of genes from this pathway with dietary nutrients and therapeutic agents in DLBCL prognosis. Figure Figure

Abstract Persons of African ancestry (AA) have a twofold higher risk for multiple myeloma (MM) compared with persons of European ancestry (EA). Genome-wide association studies (GWASs) support a genetic contribution to MM etiology in individuals of EA. Little is known about genetic risk factors for MM in individuals of AA. We performed a meta-analysis of 2 GWASs of MM in 1813 cases and 8871 controls and conducted an admixture mapping scan to identify risk alleles. We fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. In GWAS meta-analysis, we identified 2 suggestive novel loci located at 9p24.3 and 9p13.1 at P < 1 × 10−6; however, no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. Of the 23 known EA risk variants, 20 showed directional consistency, and 9 replicated at P < .05 in AA individuals. In 8 regions, we identified markers that better capture MM risk in persons with AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95% confidence interval, 1.56-2.11) increased MM risk compared with those with average risk (25%-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P = 5.1 × 10−12). Our study shows that common genetic variation contributes to MM risk in individuals with AA.

Abstract Background Recent genome wide association studies (GWAS) of Hodgkin lymphoma (HL) have identified several associations at both HLA and non-HLA loci. However, much of HL heritability remains unexplained. Methods To identify novel risk loci, we performed a meta-analysis of 3 HL GWAS including a total of 1,810 cases and 7,879 controls. Results were replicated in an independent set of 1,163 cases and 2,580 controls, for a total of 3,097 and 11,097 cases and controls combined, respectively. participants in discovery and replication stages were of European descent. quality control and imputation we conducted a meta-analysis addressing 1,004,829 variants (λ= 1.10, λ1000= 1.03). Associations between SNP genotypes and HL risk were evaluated under a log-additive model of inheritance adjusting for sex, study center and significant principal components to control for population stratification. We performed an analysis with all HL cases and then conducted stratified analyses by histological subtype (classical, nodular sclerosis and mixed cellularity), age at diagnosis (nodular sclerosis among those diagnosed at 15- 35 years in all studies, and those diagnosed at 35 and older years in the European Study only) and EBV tumor status (negative and positive). We then used a bioinformatic approach (FunciSNP) to identify potential functional variants associated with HD risk correlated with risk loci of interest. We extracted publically available ENCODE data on biofeatures to identify potential functional motifs associated with the index SNP or correlated SNPs. Finally, we measured expression levels of the two alternative mRNA transcripts in lymphoblastoid cell lines (LCLs) from 49 post-therapy HL patients and 25 unaffected controls. RT-PCR was carried out in triplicate. Relative expression levels were calculated relative to TBP as housekeeping gene. Linear models were used to assess correlation between genotype and TCF3expression levels. Results The meta-analysis identified a novel susceptibility variant at chromosome 19p13.3 (rs1860661) associated with HL risk (Odds Ratio [OR]= 0.78, P=2.0*10-8, I2=0%). variant is located in intron 2 of TCF3 (also known as E2A), a regulator of B- and T-cell lineage commitment. was also significantly associated with HL (OR= 0.85, P=0.002) in the replication series of 1,281 cases and 3,218 controls. the combined analysis consisting of the discovery and replication sets, rs1860661was strongly associated with HL (OR=0.81,=3.5*10-10), with no evidence of heterogeneity between contributing studies (Phom=0.41, I2=0%). The number of G alleles defined by rs1860661 was significantly associated with a reduced risk of each HL subtype except EBV positive HL. rs1860661 and two correlated SNPs, rs10413888 (r2=0.90) and rs8103453 (r2=0.89) identified by FunciSNP analysis map in or near marks of open chromatin and in DNAse hypersensitivity sites in TCF3 in CD20+ B cell lines., the protective minor alleles of these SNPs as defined by the G-G-G haplotype map to the binding sites for ZBTB7a (rs10413888 and rs1860661) and (rs8103453) transcription factors, likely improving the binding efficiency to the sites which may result in increased transcription rates of TCF3. TCF3 is encoded by two alternative transcripts (E12 and E47). Higher expression levels of TCF3-E47, whose transcription start site is located close to rs1860661, was associated with the rs1860661-G allele in controls (P=0.02), but not in HL patients (P=0.22). Conclusion/Discussion TCF3 is essential for the commitment of lymphoid progenitors to both B-cell and T-cell lineage development. A molecular and phenotypic hallmark of classical HL is the loss of the B-cell phenotype in HRS cells, including lack of demonstrable B-cell receptor and most B-cell specific markers such as CD19 or CD20. HRS cells have a low level of TCF3, particularly homodimers of the isoform E47, due to expression of the ABF-1 and ID2 inhibitors that bind to TCF3. Thus, higher TCF3 levels in HRS precursor cells may lead to enhanced retention of the B cell phenotype, thereby conferring a protective effect. These data suggest a link between the 19p13.3 locus including TCF3 and HL risk, indicating that TCF3 could be relevant to HL etiology and pathogenesis. Disclosures: Link: Genentech: Consultancy; Millenium: Consultancy; Pharmacyclics: Consultancy; Spectrum: Consultancy.