Wang, Sophia

Abstract A growing number of loci within the human leukocyte antigen (HLA) region have been implicated in non-Hodgkin lymphoma (NHL) etiology. Here, we test a complementary hypothesis of “heterozygote advantage” regarding the role of HLA and NHL, whereby HLA diversity is beneficial and homozygous HLA loci are associated with increased disease risk. HLA alleles at class I and II loci were imputed from genome-wide association studies (GWAS) using SNP2HLA for 3,617 diffuse large B-cell lymphomas (DLBCL), 2,686 follicular lymphomas (FL), 2,878 chronic lymphocytic leukemia/small lymphocytic lymphomas (CLL/SLL), 741 marginal zone lymphomas (MZL), and 8,753 controls of European descent. Both DLBCL and MZL risk were elevated with homozygosity at class I HLA-B and -C loci (OR DLBCL = 1.31, 95% CI = 1.06–1.60; OR MZL = 1.45, 95% CI = 1.12–1.89) and class II HLA-DRB1 locus (OR DLBCL = 2.10, 95% CI = 1.24–3.55; OR MZL = 2.10, 95% CI = 0.99–4.45). Increased FL risk was observed with the overall increase in number of homozygous HLA class II loci (P trend < 0.0001, FDR = 0.0005). These results support a role for HLA zygosity in NHL etiology and suggests that distinct immune pathways may underly the etiology of the different NHL subtypes. Significance: HLA gene diversity reduces risk for non-Hodgkin lymphoma. Cancer Res; 78(14); 4086–96. ©2018 AACR.

Objective Determinants of the increased risk of diffuse large B-cell lymphoma (DLBCL) in SLE are unclear. Using data from a recent lymphoma genome-wide association study (GWAS), we assessed whether certain lupus-related single nucleotide polymorphisms (SNPs) were also associated with DLBCL. Methods GWAS data on European Caucasians from the International Lymphoma Epidemiology Consortium (InterLymph) provided a total of 3857 DLBCL cases and 7666 general-population controls. Data were pooled in a random-effects meta-analysis. Results Among the 28 SLE-related SNPs investigated, the two most convincingly associated with risk of DLBCL included the CD40 SLE risk allele rs4810485 on chromosome 20q13 (OR per risk allele=1.09, 95% CI 1.02 to 1.16, p=0.0134), and the HLA SLE risk allele rs1270942 on chromosome 6p21.33 (OR per risk allele=1.17, 95% CI 1.01 to 1.36, p=0.0362). Of additional possible interest were rs2205960 and rs12537284. The rs2205960 SNP, related to a cytokine of the tumour necrosis factor superfamily TNFSF4, was associated with an OR per risk allele of 1.07, 95% CI 1.00 to 1.16, p=0.0549. The OR for the rs12537284 (chromosome 7q32, IRF5 gene) risk allele was 1.08, 95% CI 0.99 to 1.18, p=0.0765. Conclusions These data suggest several plausible genetic links between DLBCL and SLE.

Abstract Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk. Cancer Epidemiol Biomarkers Prev; 24(11); 1680–91. ©2015 AACR.

Abstract Background: Excess adiposity has been associated with lymphomagenesis, possibly mediated by increased cytokine production causing a chronic inflammatory state. The relationship between obesity, cytokine polymorphisms, and selected mature B-cell neoplasms is reported. Method: Data on 4,979 cases and 4,752 controls from nine American/European studies from the InterLymph consortium (1988–2008) were pooled. For diffuse large B-cell lymphoma (DLBCL), follicular lymphoma (FL), and chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL), joint associations of body mass index (from self-reported height and weight) and 12 polymorphisms in cytokines IL1A (rs1800587), IL1B (rs16944, rs1143627), IL1RN (rs454078), IL2 (rs2069762), IL6 (rs1800795, rs1800797), IL10 (rs1800890, rs1800896), TNF (rs1800629), LTA (rs909253), and CARD15 (rs2066847) were investigated using unconditional logistic regression. BMI-polymorphism interaction effects were estimated using the relative excess risk due to interaction (RERI). Results: Obesity (BMI ≥ 30 kg/m2) was associated with DLBCL risk [OR = 1.33; 95% confidence interval (CI), 1.02–1.73], as was TNF-308GA+AA (OR = 1.24; 95% CI, 1.07–1.44). Together, being obese and TNF-308GA+AA increased DLBCL risk almost 2-fold relative to those of normal weight and TNF-308GG (OR = 1.93; 95% CI, 1.27–2.94), with a RERI of 0.41 (95% CI, −0.05–0.84; Pinteraction = 0.13). For FL and CLL/SLL, no associations with obesity or TNF-308GA+AA, either singly or jointly, were observed. No evidence of interactions between obesity and the other polymorphisms were detected. Conclusions: Our results suggest that cytokine polymorphisms do not generally interact with BMI to increase lymphoma risk but obesity and TNF-308GA+AA may interact to increase DLBCL risk. Impact: Studies using better measures of adiposity are needed to further investigate the interactions between obesity and TNF-308G>A in the pathogenesis of lymphoma. Cancer Epidemiol Biomarkers Prev; 24(7); 1061–70. ©2015 AACR.

Abstract Background: Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. Relatives of DLBCL patients are at elevated risk for DLBCL. In small genome-wide association studies (GWAS) in which all NHL subtypes were combined, no conclusive loci for DLBCL have been previously identified in individuals of European background, whereas a recent study conducted in East Asia identified a locus at 3q27. Methods: To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of three new GWAS and one prior scan, totaling 3,857 DLBCL cases and 7,666 controls of European ancestry. Imputation was used to combine all the data using the 1000 Genomes Project release v3 and IMPUTE2. Based on these findings, we selected 6 promising SNPs for genotyping in an additional 202 cases and 3,431 controls, which were then included in a second meta-analysis. SNPTEST was used to estimate odds ratios (ORs), using the additive model and adjusting for age, gender and significant eigenvectors. Results: In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; OR=2.33, P=1.60x10-20), rs13255292 and rs4733601 at 8q24.21 (PVT1; OR=1.23, P=4.86x10-12 and OR=1.18, P=1.63x10-8, respectively), rs2523607 at 6p21.33 (HLA-B; OR=1.34, 3.35x10-9) and rs2681416 at 3q13.33 (CD86; OR=1.20, P=6.04x10-9). The two 8q24.21 SNPs displayed minimal linkage disequilibrium (r2=0.03) and conditional analysis supported their independence. We did not observe a notable signal for a locus previously reported for DLBCL on 3q27 in East Asia, rs6773854 (OR=1.06, P=0.81). We estimated that common SNPs, including but not limited to the loci discovered in this study, explain approximately 16% of the variance for DLBCL risk overall. Discussion: The susceptibility locus at 6p25.3 maps near a plausible DLBCL candidate gene, EXOC2 (exocyst complex component 2), which interacts with Ral proteins, and disruption of the Ral-exocyst regulatory node impacts cellular and developmental processes associated with malignant transformation and progression. The two 8q24.21 variants are approximately 1Mb telomeric to the 8q24 region linked with multiple cancers including CLL. Both variants are positioned in close proximity to PVT1, which is a non-coding RNA implicated in the MYC activation, and MYC is known to be deregulated in some DLBCLs. The locus at 3q13.33 maps to CD86 (B7-2), a B-7 family member that along with CD80 (B7-1) is essential in cognate interactions between T- and B-lymphocytes; both play key roles in immunodeficiency, autoimmune diseases, and in anti-tumor immunity. Finally, we had multiple signals in the HLA region, the strongest of which was at 6p21.33 near HLA-B. HLA-B encodes the HLA class I heavy chain paralogue, which heterodimerizes with a light chain (β2 microglobulin) and plays a central role in presenting intracellularly processed self or foreign antigens to CD8+ cytotoxic T lymphocytes. Class I molecules have been extensively linked to the risk of a variety of immune-mediated diseases and cancers including DLBCL. Conclusions: These data provide substantial new evidence for genetic susceptibility to this B-cell malignancy, and point towards pathways involved in immune recognition and immune function in the pathogenesis of DLBCL. Citation Format: James R. Cerhan, Sonja I. Berndt, Joseph Vijai, Hervé Ghesquières, James McKay, Sophia S. Wang, Zhaoming Wang, Meredith Yeager, Alexandra Nieters, David Cox, Alain Monnereau, Christopher Flowers, Anneclaire J. De Roos, Angela R. Brooks-Wilson, Qing Lan, Gianluca Severi, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark Purdue, Claire Vajdic, Demetrius Albanes, Kimberly A. Bertrand, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Paolo Vineis, Susan L. Slager, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, on behalf of the NHL GWAS Project. Genome-wide association study identifies multiple susceptibility loci for diffuse large B-cell 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 LB-272. doi:10.1158/1538-7445.AM2014-LB-272

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

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.