Rothman, Nathaniel

Abstract Background. Follicular lymphoma (FL) is an indolent B-cell malignancy with a variable clinical course that may transform to more aggressive forms of lymphoma. Previous modestly powered GWAS have reported multiple FL susceptibility loci in the HLA class I and II regions at 6p21.32-33, but novel non-HLA susceptibility loci and the genetic architecture of FL associations in the HLA region remain to be elucidated. Methods. We conducted the largest GWAS of FL to date consisting of 2,142 cases and 6,221 controls of European ancestry from 22 studies scanned on the Illumina OmniExpress. After imputation of common SNPs using IMPUTE2 and 1,000 Genomes Project v3 data, we conducted a meta-analysis of these data with two previous FL GWAS (n=586 cases, 1,537 controls) and replicated top hits in an additional 629 cases and 4,283 controls. Expression quantitative trait loci (eQTL) analyses were performed to evaluate the effects of associated variants on gene expression. Classical HLA allele imputations and stepwise regression analyses are underway to further characterize the HLA associations. Results. In the joint meta-analysis, the strongest association with FL risk was observed in the HLA region, where a large number of SNPs showed genome-wide significance (P<5x10-8 to P= 1.84x10-84), particularly in HLA class II at 6p21.32 where several HLA eQTLs were identified (P < 0.05). Outside HLA, we identified four novel loci associated with FL at 11q23.3 (CXCR5, rs4938573; P=3.14x10-15), 11q24.3 (ETS1, rs4937362; P=3.33x10-9), 3q13.33 (CD86, rs2681416; P=9.79x10-11) and 3q28 (LPP, rs6444305, P=2.55x10-8). These novel susceptibility loci were located in four biologically relevant genes involved in immune regulation: CXCR5 and its ligand, CXCL13, are involved in guiding B-cells into the B-cell zones of secondary lymphoid organs as well as T-cell migration. ETS1 is expressed in lymphoid cells and regulates immune cell function including differentiation, survival, and proliferation. CD86 and CD80 are classic members of the B7 costimulatory pathway that is important in maintaining immune function, suppression of autoimmunity and antitumor surveillance. LPP encodes a LIM domain-containing protein of the zyxin family and participates in cell adhesion, cell migration, proliferation and transcription dynamics. Conclusions. This large GWAS provides further support for the important role of common genetic variation in non-HLA genes and further evidence of the key role that HLA immune-regulatory genes play in the pathogenesis of FL. Citation Format: Christine F. Skibola, Sonja I. Berndt, James R. Cerhan, Zhaoming Wang, Joseph Vijai, Lucia Conde, Paul de Bakker, Sophia S. Wang, Claire M. Vajdic, Brenda M. Birmann, Susan L. Slager, James McKay, Paige M. Bracci, Alexandra Nieters, Qing Lan, Angela R. Brooks-Wilson, Martha S. Linet, Demetrius Albanes, John J. Spinelli, Roel C.H. Vermeulen, Mark P. Purdue, Meredith Yaeger, Lauren R. Teras, Silvia de Sanjose, Alain Monnereau, Simon Crouch, Jia Nee Foo, Henrik Hjalgrim, Gianluca Severi, Brian K. Link, Kimberly A. Bertrand, Yawei Zhang, Karin E. Smedby, Stephen J. Chanock, Nathaniel Rothman, NHL GWAS Consortium. Meta-analysis of genome-wide association studies identifies novel susceptibility loci for follicular lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5072. doi:10.1158/1538-7445.AM2014-5072

Abstract Marginal zone lymphomas (MZL) are a group of indolent non-Hodgkin's B-cell lymphomas (NHL) comprising extranodal MZL of mucosa-associated lymphoid tissue (MALT) type, splenic MZL and nodal MZL. MZL comprises 7-12% of adult NHLs. The genetic etiology of MZL is poorly understood. To better understand the genetic risk factors for MZL, we conducted the first genome-wide association study (GWAS) of this lymphoid malignancy. The GWAS included 824 MZL cases and 6,220 controls of European ancestry from 22 studies of NHL. Genotyping was carried out at the Cancer Genomics Research Laboratory (CGR) of NCI using Illumina Human OmniExpress arrays. Standard genotyping quality controls were performed. Imputation of common SNPs was performed using IMPUTE2 software and 1000 Genomes Project release version 3 data. Association testing was subsequently performed using SNPTEST, assuming a log-additive genetic model adjusting for age, gender and three significant eigenvectors. We identified one locus at 6p21.32 within the HLA class II region that attained genome-wide significance (P<5x10-8). When a small number of additional cases (n=58) and controls (n=66) genotyped using Taqman assays were included, rs9461741 (BTNL2) was associated with an increased risk of MZL (per-allele odds ratio = 2.38, 95% confidence interval: 1.74-3.24; P=4.32x10-8). Other SNPs at the 6p21.32 locus were also highly associated with MZL, suggesting that this is unlikely to be a chance finding. Promising signals yet to be confirmed were also observed on chromosomes 3 and 14. The BTNL2 protein has extracellular immunoglobulin domains and is expressed in multiple tissues including various immune organs and the gut and it inhibits T-cell activation. Genetic variability at BTNL2 has been associated with inflammatory and infectious phenotypes, highlighting the importance of this locus for immune diseases. In conclusion, findings from this GWAS identify the HLA class II region as an important susceptibility locus for MZL. Citation Format: Joseph Vijai, Zhaoming Wang, Sonja I. Berndt, Susan L. Slager, James R. Cerhan, Christine Skibola, Sophia S. Wang, Angela R. Brooks-Wilson, Silvia de Sanjose, Paige M. Bracci, Mads Melbye, Bengt Glimelius, Rebecca D. Jackson, Wendy Cozen, Nikolaus Becker, Lauren R. Teras, John J. Spinelli, Jenny Turner, Qing Lan, Mark P. Purdue, Kimberly A. Bertrand, Kenneth Offit, Paolo Vineis, Brian K. Link, Graham G. Giles, Anne Zelenuich-Jacquotte, Alain Monnereau, Maria Grazia Ennas, Demetrius Albanes, Laurie Burdett, Nathaniel Rothman, Stephen J. Chanock, Alexandra Nieters, On behalf of the NHL GWAS Consortium. A genome-wide association study suggests evidence of variants at 6p21.32 associated with marginal zone lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5071. doi:10.1158/1538-7445.AM2014-5071

AbstractObjectivesDiffuse large B‐cell lymphoma (DLBCL) is an aggressive subtype of non‐Hodgkin lymphoma (NHL) and is the most common NHL subtype diagnosed worldwide. The first large‐scale genome‐wide association study (GWAS) of DLBCL with over 4000 cases conducted among individuals of European ancestry recently identified five independent SNPs that achieved genome‐wide significance, and two SNPs that showed a suggestive association with DLBCL risk.MethodsTo evaluate whether Eastern Asians and individuals of European ancestry share similar genetic risk factors for this disease, we attempted to replicate these GWAS findings in a pooled series of 1124 DLBCL cases and 3596 controls from Hong Kong, South Korea, and Thailand.ResultsThree of the five genome‐wide significant SNPs from the DLBCL GWAS were significantly associated with DLBCL in our study population, including the top finding from the GWAS, EXOC2 rs116446171, which achieved genome‐wide significance in our data (per allele OR = 2.04, 95% CI = 1.63–2.56; ptrend = 3.9 × 10−10). Additionally, we observed a significant association with PVT1 rs13255292 (per allele OR = 1.34, 95% CI = 1.19–1.52; ptrend = 2.1 × 10−6), which was the second strongest finding in the GWAS, and with HLA‐B rs2523607 (per allele OR = 3.05, 95% CI = 1.32–7.05; ptrend = 0.009).ConclusionsOur study, which provides the first evaluation in Eastern Asians of SNPs definitively associated with DLBCL risk in individuals of European ancestry, indicates that at least some of the genetic factors associated with risk of DLBCL are similar between these populations.

Abstract Despite limited discovery stages (<1,125 cases), genome-wide association studies (GWAS) have successfully identified 13 loci associated with risk of chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL); however, these loci only explain a portion of the heritability. To identify additional genetic loci associated with CLL susceptibility, we conducted the largest meta-analysis of GWAS, to date, including 3,100 CLL cases and 7,667 controls of European ancestry from 25 independent studies genotyped on four different SNP microarrays. Association testing was conducted assuming a log-additive genetic model, adjusting for age, sex, and significant principal components. Meta-analyses were performed using the fixed effects inverse variance method based on the beta estimates and standard errors from each study. Promising loci were taken forward for replication in three additional studies totaling 788 cases and 4,872 controls. In the meta-analysis of both the discovery and replication studies, we conclusively identified ten independent SNPs in nine novel loci that reached genome-wide significance (P<5 x 10−8): 10q23.31 (FAS; P=1.22x10−14), 18q21.33 (BCL2; P=2.66x10−12 and P=7.76x10−11), 11p15.5 (C11orf21; P=2.15x10−10), 4q25 (LEF1; P=4.24x10−10), 2q33.1 (CASP8; P=2.50x10−9), 9p21.3 (CDKN2B-AS1; P=1.27x10−8), 18q21.32 (PMAIP1; P=2.51x10−8), 15q15.1 (BMF; P=2.71x10−10), and 2p22.2 (QPCT; P=1.68x10−8) as well as an independent signal at an established locus (2q13, ACOXL, P=2.08x10−18). We also found evidence for two additional promising loci that reached marginal genome-wide significance (P<2.0x10−7) at 8q22.3 (ODF1; P=5.40x10−8) and 5p15.33 (TERT; P=1.92x10−7). Although further studies are required, the proximity of five of the novel loci to genes involved in apoptosis suggests a possible unifying underlying biological mechanism. Together with the previously established loci, these genetic loci explain ~17% of the familial risk of CLL. Citation Format: Sonja I. Berndt, Christine F. Skibola, Vijai Joseph, Nicola J. Camp, Alexandra Nieters, Zhaoming Wang, Wendy Cozen, Alain Monnereau, Rachel S. Kelly, Silvia de Sanjose, Martha S. Linet, Elio Riboli, Paolo Vineis, Henrik Hjalgrim, James R. Cerhan, Stephen J. Chanock, Nathaniel Rothman, Susan L. Slager. Meta-analysis of genome-wide association studies identifies multiple loci associated with chronic lymphocytic leukemia. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr LB-23. doi:10.1158/1538-7445.AM2013-LB-23

Objectives Growing evidence suggests that gender-blind assessment of exposure may introduce exposure misclassification, but few studies have characterised gender differences across occupations and industries. We pooled control responses to job-specific, industry-specific and exposure-specific questionnaires (modules) that asked detailed questions about work activities from three US population-based case–control studies to examine gender differences in work tasks and their frequencies. Methods We calculated the ratio of female-to-male controls that completed each module. For four job modules (assembly worker, machinist, health professional, janitor/cleaner) and for subgroups of jobs that completed those modules, we evaluated gender differences in task prevalence and frequency using χ2 and Mann–Whitney U tests, respectively. Results The 1360 female and 2245 male controls reported 6033 and 12 083 jobs, respectively. Gender differences in female:male module completion ratios were observed for 39 of 45 modules completed by ≥20 controls. Gender differences in task prevalence varied in direction and magnitude. For example, female janitors were significantly more likely to polish furniture (79% vs 44%), while male janitors were more likely to strip floors (73% vs 50%). Women usually reported more time spent on tasks than men. For example, the median hours per week spent degreasing for production workers in product manufacturing industries was 6.3 for women and 3.0 for men. Conclusions Observed gender differences may reflect actual differences in tasks performed or differences in recall, reporting or perception, all of which contribute to exposure misclassification and impact relative risk estimates. Our findings reinforce the need to capture subject-specific information on work tasks.

AbstractMany common genetic variants have been associated with non-Hodgkin lymphoma (NHL), but individual study results are often conflicting. To confirm the role of putative risk alleles in B-cell NHL etiology, we performed a validation genotyping study of 67 candidate single nucleotide polymorphisms within InterLymph, a large international consortium of NHL case-control studies. A meta-analysis was performed on data from 5633 B-cell NHL cases and 7034 controls from 8 InterLymph studies. rs3789068 in the proapoptotic BCL2L11 gene was associated with an increased risk for B-cell NHL (odds ratio = 1.21, P random = 2.21 × 10−11), with similar risk estimates for common B-cell subtypes. PRRC2A rs3132453 in the HLA complex class III region conferred a reduced risk of B-cell NHL (odds ratio = 0.68, P random = 1.07 × 10−9) and was likewise evident for common B-cell subtypes. These results are consistent with the known biology of NHL and provide insights into shared pathogenic components, including apoptosis and immune regulation, for the major B-cell lymphoma subtypes.

Genetic variation in the 6p21 chromosomal region, including human leukocyte antigen (HLA) genes and tumor necrosis factor (TNF), has been linked to both etiology and clinical outcomes of lymphomas. We estimated the effects ofHLAclass I (A, B, and C), class IIDRB1alleles, and the ancestral haplotype (AH) 8.1 (HLAA*01-B*08-DRB1*03-TNF-308A) on overall survival (OS) among patients with diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL) in a population-based study of non-Hodgkin lymphoma. During a median followup of 89 months, 31% (52 of 166) DLBCL and 28% (46 of 165) FL patients died. Using multivariate Cox regression models, we observed statistically significant associations between genetic variants and survival:HLA-Cw*07:01was associated with poorer OS among DLBCL patients (Hazard ratio [HR] = 1.76, 95% confidence interval [CI] = 1.01–3.05);HLA-A*01:01was associated with poorer OS (HR = 2.23, 95% CI = 1.24–4.01), andHLA-DRB1*13(HR = 0.12, 95% CI = 0.02–0.90) andHLA-B Bw4(HR = 0.36, 95% CI = 0.20–0.63) with better OS among FL patients. These results support a role forHLAin the prognosis of DLBCL and FL and represent a promising class of prognostic factors that warrants further evaluation.

Abstract Background. Non-Hodgkin lymphoma (NHL) case-control studies consistently demonstrate the G-308A promoter polymorphism in the tumor necrosis factor gene on chromosome 6p21.3 to be associated with increased risk of the NHL subtype, diffuse large B-cell lymphoma (DLBCL). Initial results from a genome-wide association study have implicated a polymorphism located on 6p21.33 with increased risk of the follicular lymphoma subtype. In the present analysis, we evaluated 488 candidate gene regions (defined as 20kb 5’, 10kb 3’) and their relation to risk for NHL and the DLBCL and follicular lymphoma subtypes with a particular focus on genes within the 6p21 chromosomal region. Methods. We genotyped 6679 tag single nucleotide polymorphisms (SNPs) in 947 cases and 826 population-based controls from a multicenter U.S.-based case-control study. Genotyping was conducted on a custom-designed Illumina Infinium assay at the NCI Core Genotyping Facility. We obtained a gene region-level summary of association by computing the minimum p-value (“minP test”). We used logistic regression to evaluate the association between single nucleotide polymorphism (SNP) genotypes and haplotypes with NHL (adjusted for sex, race and age). For NHL subtypes, we conducted polytomous multivariate unconditional logistic regression. We evaluated tests for trend under the co-dominant model used a three-level ordinal variable for each SNP. Results. We identified 34 gene regions significantly associated with all NHL (p <0.05). The three most significant gene regions were all found on chromosome 6p21.3 (RING1/RXRB; AIF1/BAT2; BAT4, p=0.001). Within those gene regions, the most significant SNPs (p<0.0002) were rs2855429 and rs213213 (RING1/RXRB), rs2857597 and rs3130071 (AIF1/BAT2), and rs3115667 (BAT4). All five SNPs were statistically significantly associated with both DLBCL and follicular lymphoma. Additional gene regions on 6p21.3 associated with NHL include MSH5 (p=0.009), FAM46A (p=0.01) and TNXB/CREBL1 (p=0.04). Stratification by age and gender did not alter our results; consistent results were observed in analyses restricted to non-Hispanic Caucasians. Haplotype analyses yielded the same regions and SNPs of significance. Detailed analysis of the 6p21.3 region will be presented. Conclusions. Our results further implicate gene variants within the 6p21.3 chromosomal region in NHL etiology. Full characterization of the 6p21.3 region and identification of functional variants within this region are warranted for understanding lymphomagenesis. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 935.

Abstract Background: There is growing evidence linking genetic variations to non–Hodgkin lymphoma (NHL) etiology. To complement ongoing agnostic approaches for identifying susceptibility genes, we evaluated 488 candidate gene regions and their relation to risk for NHL and NHL subtypes. Methods: We genotyped 6,679 tag single nucleotide polymorphisms (SNPs) in 947 cases and 826 population-based controls from a multicenter U.S. case–control study. Gene-level summary of associations were obtained by computing the minimum P value (“minP test”) on the basis of 10,000 permutations. We used logistic regression to evaluate the association between genotypes and haplotypes with NHL. For NHL subtypes, we conducted polytomous multivariate unconditional logistic regression (adjusted for sex, race, age). We calculated P-trends under the codominant model for each SNP. Results: Fourteen gene regions were associated with NHL (P < 0.01). The most significant SNP associated with NHL maps to the SYK gene (rs2991216, P-trend = 0.00005). The three most significant gene regions were on chromosome 6p21.3 (RING1/RXRB; AIF1; BAT4). Accordingly, SNPs in RING1/RXRB (rs2855429), AIF1 (rs2857597), and BAT4 (rs3115667) were associated with NHL (P-trends ≤ 0.0002) and both diffuse large B-cell and follicular lymphomas (P-trends < 0.05). Conclusions: Our results suggest potential importance for SYK on chromosome 9 with NHL etiology. Our results further implicate 6p21.3 gene variants, supporting the need for full characterization of this chromosomal region in relation to lymphomagenesis. Impact: Gene variants on chromosome 9 may represent a new region of interesting for NHL etiology. The independence of the reported variants in 6p21.3 from implicated variants (TNF/HLA) supports the need to confirm causal variants in this region Cancer Epidemiol Biomarkers Prev; 20(1); 42–9. ©2011 AACR.