GWAS of follicular lymphoma reveals allelic heterogeneity at 6p21.32 and suggests shared genetic susceptibility with diffuse large B-cell lymphoma

HLA variations in patients with diffuse large B-cell lymphoma and association with disease risk and prognosis: a case-control study

Introduction

Human leukocyte antigen (HLA) polymorphisms have been associated with the development of various autoimmune diseases, as well as malignant neoplasms. Non-Hodgkin lymphomas (NHLs) are a heterogenous group of lymphoid malignancies in which a genetic substrate has been established and is deemed to play a crucial role in disease pathogenesis. This study aimed to identify whether variations in the HLA gene region were associated with diffuse large B-cell lymphoma (DLBCL) risk and prognosis.

Methods

We defined HLA class I (HLA-A, HLA-B, HLA-C) and class II (HLA-DRB1, HLA-DQB1) alleles in 60 patients with DLBCL and compared the results to those found by 236 healthy adult donors from the bone marrow bank of Northern Greece. HLA typing was performed by two molecular methods, Sequence - Specific Oligonucleotide HLA typing (SSO) and Sequence - Specific Primer HLA typing (SSP), from white blood cells recovered from peripheral blood. The phenotypic frequencies of HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 between patients and controls were compared with the 2-sided Fisher's exact test. Results with p-value <0.05 were considered statistically significant. Odds Ratios with 95% Confidence Intervals were calculated to further strengthen the results. The 2-sided Fisher's exact test was also applied to alleles found only in one of the two groups, while the odds ratios together with the confidence intervals were corrected with Haldane-Anscombe method.

Results

Among the studied HLA polymorphisms, the frequency HLA-C*12 allele was significantly lower in patients with DLBCL compared with control subjects (6.7% vs. 34.7%, OR = 0.16, 95% CI: 0.04-0.44). Frequency of HLA-B*39 was significantly lower in patients with DLBCL compared with controls, but due to the low frequency of this polymorphism in the studied population and small sample size, determinations regarding the significance of this findings were limited. Survival analysis revealed that the presence of HLA-C*12 was not associated with improved or worsened overall and progression-free survival. No statistically significant associations were observed in the phenotypic frequencies of HLA-A, HLA-DQB1, HLA-DRB1 and the rest of HLA-B alleles between the control and DLBCL groups.

Discussion

Collectively, our results provide valuable insight regarding the role of HLA variations on DLBCL risk. Further studies are required to consolidate our findings and ascertain the clinical implications of these genetic variations on DLBCL management and prognosis.

Epidemiology and etiology of diffuse large B-cell lymphoma

As the most common non-Hodgkin lymphoma subtype, diffuse large B-cell lymphoma (DLBCL) incidence patterns generally parallel that for NHL overall. Globally, DLBCL accounts for a third of all NHLs, ranging between 20% and 50% by country. Based on United States (U.S.) cancer registry data, age-standardized incidence rate for DLBCL was 7.2 per 100,000. DLBCL incidence rises with age and is generally higher in males than females; in the U.S., incidence is highest among non-Hispanic whites (9.2/100,000). Like NHL incidence, DLBCL incidence rose in the first half of the 20th century but has largely plateaued. However, there is some evidence that incidence rates are rising in areas of historically low rates, such as Asia; there are also estimates for rising DLBCL incidence in the near future due to the changing demographics in developed countries whose aging population is growing. Established risk factors for DLBCL include those that result in severe immune deficiency such as HIV/AIDS, inherited immunodeficiency syndromes, and organ transplant recipients. Factors that lead to chronic immune dysregulations are also established risk factors, and include a number of autoimmune conditions (eg, Sjögren syndrome, systemic lupus erythematosus, rheumatoid arthritis), viral infections (eg, HIV, KSHV/HHV8, HCV, EBV), and obesity. Family history of NHL/DLBCL, personal history of cancer, and multiple genetic susceptibility loci are also well-established risk factors for DLBCL. There is strong evidence for multiple environmental exposures in DLBCL etiology, including exposure to trichloroethylene, benzene, and pesticides and herbicides, with recent associations noted with glyphosate. There is also strong evidence for associations with other viruses, such as HBV. Recent estimates suggest that obesity accounts for nearly a quarter of DLBCLs that develop, but despite recent gains in the understanding of DLBCL etiology, the majority of disease remain unexplained. An understanding of the host and environmental contributions to disease etiology, and concerted efforts to expand our understanding to multiple race/ethnic groups, will be essential for constructing clinically relevant risk prediction models and develop effective strategies for disease prevention.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

Assessment of genetic susceptibility to multiple primary cancers through whole-exome sequencing in two large multi-ancestry studies

BACKGROUND

Up to one of every six individuals diagnosed with one cancer will be diagnosed with a second primary cancer in their lifetime. Genetic factors contributing to the development of multiple primary cancers, beyond known cancer syndromes, have been underexplored.

METHODS

To characterize genetic susceptibility to multiple cancers, we conducted a pan-cancer, whole-exome sequencing study of individuals drawn from two large multi-ancestry populations (6429 cases, 165,853 controls). We created two groupings of individuals diagnosed with multiple primary cancers: (1) an overall combined set with at least two cancers across any of 36 organ sites and (2) cancer-specific sets defined by an index cancer at one of 16 organ sites with at least 50 cases from each study population. We then investigated whether variants identified from exome sequencing were associated with these sets of multiple cancer cases in comparison to individuals with one and, separately, no cancers.

RESULTS

We identified 22 variant-phenotype associations, 10 of which have not been previously discovered and were significantly overrepresented among individuals with multiple cancers, compared to those with a single cancer.

CONCLUSIONS

Overall, we describe variants and genes that may play a fundamental role in the development of multiple primary cancers and improve our understanding of shared mechanisms underlying carcinogenesis.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

Background

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

Case presentation

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

Conclusions

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

Genetic predisposition to lymphomas: Overview of rare syndromes and inherited familial variants

Approximately 10 % of malignancies occur in carriers of germline mutations predisposing to cancer. A high risk of developing lymphomas has been noted in many primary immunodeficiencies, including DNA repair disorders. Moreover, implementation of next-generation sequencing has recently enabled to uncover rare genetic variants predisposing patients to lymphoid neoplasms. Some patients harboring inherited predisposition to lymphomas require dedicated clinical management, which will contribute to effective cancer treatment and to the prevention of potential severe toxicities and secondary malignancies. In line with that, our review summarizes the natural history of lymphoid tumors developing on different germline genetic backgrounds and discusses the progress that has been made toward successfully treating these malignancies.

Genome-wide homozygosity and risk of four non-Hodgkin lymphoma subtypes

AIM

Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.

METHODS

We used data from eight genome-wide association studies of four common NHL subtypes: 3061 chronic lymphocytic leukemia (CLL), 3814 diffuse large B-cell lymphoma (DLBCL), 2784 follicular lymphoma (FL), and 808 marginal zone lymphoma (MZL) cases, as well as 9374 controls. We examined the effect of homozygous variation on risk by: (1) estimating the fraction of the autosome containing runs of homozygosity (FROH); (2) calculating an inbreeding coefficient derived from the correlation among uniting gametes (F3); and (3) examining specific autosomal regions containing ROH. For each, we calculated beta coefficients and standard errors using logistic regression and combined estimates across studies using random-effects meta-analysis.

RESULTS

We discovered positive associations between FROH and CLL (β = 21.1, SE = 4.41, P = 1.6 × 10-6) and FL (β = 11.4, SE = 5.82, P = 0.02) but not DLBCL (P = 1.0) or MZL (P = 0.91). For F3, we observed an association with CLL (β = 27.5, SE = 6.51, P = 2.4 × 10-5). We did not find evidence of associations with specific ROH, suggesting that the associations observed with FROH and F3 for CLL and FL risk were not driven by a single region of homozygosity.

CONCLUSION

Our findings support the role of recessive genetic variation in the etiology of CLL and FL; additional research is needed to identify the specific loci associated with NHL risk.

Sequencing at lymphoid neoplasm susceptibility loci maps six myeloma risk genes

Inherited genetic risk factors play a role in multiple myeloma (MM), yet considerable missing heritability exists. Rare risk variants at genome-wide association study (GWAS) loci are a new avenue to explore. Pleiotropy between lymphoid neoplasms (LNs) has been suggested in family history and genetic studies, but no studies have interrogated sequencing for pleiotropic genes or rare risk variants. Sequencing genetically enriched cases can help discover rarer variants. We analyzed exome sequencing in familial or early-onset MM cases to identify rare, functionally relevant variants near GWAS loci for a range of LNs. A total of 149 distinct and significant LN GWAS loci have been published. We identified six recurrent, rare, potentially deleterious variants within 5 kb of significant GWAS single nucleotide polymorphisms in 75 MM cases. Mutations were observed in BTNL2, EOMES, TNFRSF13B, IRF8, ACOXL and TSPAN32. All six genes replicated in an independent set of 255 early-onset MM or familial MM or precursor cases. Expansion of our analyses to the full length of these six genes resulted in a list of 39 rare and deleterious variants, seven of which segregated in MM families. Three genes also had significant rare variant burden in 733 sporadic MM cases compared with 935 control individuals: IRF8 (P = 1.0 × 10-6), EOMES (P = 6.0 × 10-6) and BTNL2 (P = 2.1 × 10-3). Together, our results implicate six genes in MM risk, provide support for genetic pleiotropy between LN subtypes and demonstrate the utility of sequencing genetically enriched cases to identify functionally relevant variants near GWAS loci.

Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts

Even distinct cancer types share biological hallmarks. Here, we investigate polygenic risk score (PRS)-specific pleiotropy across 16 cancers in European ancestry individuals from the Genetic Epidemiology Research on Adult Health and Aging cohort (16,012 cases, 50,552 controls) and UK Biobank (48,969 cases, 359,802 controls). Within cohorts, each PRS is evaluated in multivariable logistic regression models against all other cancer types. Results are then meta-analyzed across cohorts. Ten positive and one inverse cross-cancer associations are found after multiple testing correction. Two pairs show bidirectional associations; the melanoma PRS is positively associated with oral cavity/pharyngeal cancer and vice versa, whereas the lung cancer PRS is positively associated with oral cavity/pharyngeal cancer, and the oral cavity/pharyngeal cancer PRS is inversely associated with lung cancer. Overall, we validate known, and uncover previously unreported, patterns of pleiotropy that have the potential to inform investigations of risk prediction, shared etiology, and precision cancer prevention strategies.

Tracking the Genetic Susceptibility Background of B-Cell Non-Hodgkin's Lymphomas from Genome-Wide Association Studies

B-cell non-Hodgkin’s lymphoma (NHL) risk associations had been mainly attributed to family history of the disease, inflammation, and immune components including human leukocyte antigen (HLA) genetic variations. Nevertheless, a broad range of genome-wide association studies (GWAS) have shed light into the identification of several genetic variants presumptively associated with B-cell NHL etiologies, survival or shared genetic risk with other diseases. The present review aims to overview HLA structure and diversity and summarize the evidence of genetic variations, by GWAS, on five NHL subtypes (diffuse large B-cell lymphoma DLBCL, follicular lymphoma FL, chronic lymphocytic leukemia CLL, marginal zone lymphoma MZL, and primary central nervous system lymphoma PCNSL). Evidence indicates that the HLA zygosity status in B-cell NHL might promote immune escape and that genome-wide significance variants can give biological insight but also potential therapeutic markers such as WEE1 in DLBCL. However, additional studies are needed, especially for non-DLBCL, to replicate the associations found to date.

Pan-cancer study detects genetic risk variants and shared genetic basis in two large cohorts

Deciphering the shared genetic basis of distinct cancers has the potential to elucidate carcinogenic mechanisms and inform broadly applicable risk assessment efforts. Here, we undertake genome-wide association studies (GWAS) and comprehensive evaluations of heritability and pleiotropy across 18 cancer types in two large, population-based cohorts: the UK Biobank (408,786 European ancestry individuals; 48,961 cancer cases) and the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging cohorts (66,526 European ancestry individuals; 16,001 cancer cases). The GWAS detect 21 genome-wide significant associations independent of previously reported results. Investigations of pleiotropy identify 12 cancer pairs exhibiting either positive or negative genetic correlations; 25 pleiotropic loci; and 100 independent pleiotropic variants, many of which are regulatory elements and/or influence cross-tissue gene expression. Our findings demonstrate widespread pleiotropy and offer further insight into the complex genetic architecture of cross-cancer susceptibility.

Infectious mononucleosis, immune genotypes, and non-Hodgkin lymphoma (NHL): an InterLymph Consortium study

PURPOSE

We explored the interaction between non-Hodgkin lymphoma (NHL), infectious mononucleosis (IM) history, and immune-related genotypes in a pooled case-control analysis.

METHODS

A total of 7,926 NHL patients and 10,018 controls from 12 case-control studies were included. Studies were conducted during various time periods between 1988 and 2008, and participants were 17-96 years of age at the time of ascertainment/recruitment. Self-reported IM history and immune response genotypes were provided by the InterLymph Data Coordinating Center at Mayo Clinic. Odds ratios (OR) were estimated using multivariate logistic regression, and interactions were estimated using the empirical Bayes method. PACT was used to account for multiple comparisons.

RESULTS

There was evidence of an interaction effect between IM history and two variants on T-cell lymphoma (TCL) risk: rs1143627 in interleukin-1B (IL1B) (pinteraction = 0.04, ORinteraction = 0.09, 95% confidence interval [CI] 0.01, 0.87) and rs1800797 in interleukin-6 (IL6) (pinteraction = 0.03, ORinteraction = 0.08, 95% CI 0.01, 0.80). Neither interaction effect withstood adjustment for multiple comparisons. There were no statistically significant interactions between immune response genotypes and IM on other NHL subtypes.

CONCLUSIONS

Genetic risk variants in IL1B and IL6 may affect the association between IM and TCL, possibly by influencing T-cell activation, growth, and differentiation in the presence of IM, thereby decreasing risk of immune cell proliferation.

The role of HLA variation in lymphoma aetiology and survival

Epidemiologic and laboratory evidence has consistently supported a strong inflammatory and immune component for lymphoma aetiology. These studies have consistently implicated variation in the immune gene, human leucocyte antigen (HLA), to be associated with lymphoma risk. In this review, we summarize the historical and recent evidence of HLA in both lymphoma aetiology and survival. The recent momentum in uncovering HLA associations has been propelled by the conduct of genome-wide association studies (GWAS), which has permitted the evaluation of imputed HLA alleles in much larger sample sizes than historically feasible with allelotyping studies. Based on the culmination of smaller HLA typing studies and larger GWAS, we now recognize several HLA associations with Hodgkin (HL) and non-Hodgkin lymphomas (NHLs) and their subtypes. Although other genetic variants have also been implicated with lymphoma risk, it is notable that HLA associations have been reported in every NHL and HL subtype evaluated to date. Both HLA class I and class II alleles have been linked with NHL and HL risk. It is notable that the associations identified are largely specific to each lymphoma subtype. However, pleiotropic HLA associations have also been observed. For example, rs10484561, which is in linkage disequilibrium with HLA-DRB1*01:01˜DQA1*01:01˜DQB1*05:01, has been implicated in increased FL and DLBCL risk. Opposing HLA associations across subtypes have also been reported, such as for HLA-A*01:01 which is associated with increased risk of EBV-positive cHL but decreased risk of EBV-negative cHL and chronic lymphocytic leukaemia/small cell lymphoma. Due to extensive linkage disequilibrium and allele/haplotypic variation across race/ethnicities, identification of causal alleles/haplotypes remains challenging. Follow-up functional studies are needed to identify the specific immunological pathways responsible in the multifactorial aetiology of HL and NHL. Correlative studies linking HLA alleles with known molecular subtypes and HLA expression in the tumours are also needed. Finally, additional association studies investigating HLA diversity and lymphoma survival are also required to replicate initial associations reported to date.

The association between HLA and non-Hodgkin lymphoma subtypes, among a transplant-indicated population

Several studies have implicated HLA in non-Hodgkin lymphoma (NHL) subtype etiology. However, NHL patients indicated for stem cell transplants are underrepresented in these reports. We therefore evaluated the association between HLA and NHL subtypes among a transplant-indicated population. One thousand three hundred and sixty-six NHL patients HLA-typed and indicated for transplant at the City of Hope National Medical Center (Duarte, CA) were compared to 10,271 prospective donors. Odds ratios and 95% confidence intervals were calculated for HLA haplotype and alleles, adjusted for sex and age. The HLA-A*0201∼C*0602∼B*1302∼DRB1*0701∼DQB1*0201 haplotype was significantly associated with follicular lymphoma (FL) risk among Caucasians. Several haplotypes were associated with diffuse large B-cell lymphoma (DLBCL) risk among Caucasians, including the previously implicated DLBCL risk loci, HLA-B*0801. The HLA-A*0101 allele was also observed to be associated with mantle cell lymphoma (MCL) risk. Our results support the association between previously reported susceptibility loci and FL and suggest potentially new DLBCL and MCL risk loci.

Genome-wide association studies identify susceptibility loci for epithelial ovarian cancer in east Asian women

OBJECTIVE

Genome-wide association studies (GWASs) for epithelial ovarian cancer (EOC) have focused largely on populations of European ancestry. We aimed to identify common germline variants associated with EOC risk in Asian women.

METHODS

Genotyping was performed as part of the OncoArray project. Samples with >60% Asian ancestry were included in the analysis. Genotyping was performed on 533,631 SNPs in 3238 Asian subjects diagnosed with invasive or borderline EOC and 4083 unaffected controls. After imputation, genotypes were available for 11,595,112 SNPs to identify associations.

RESULTS

At chromosome 6p25.2, SNP rs7748275 was associated with risk of serous EOC (odds ratio [OR] = 1.34, P = 8.7 × 10-9) and high-grade serous EOC (HGSOC) (OR = 1.34, P = 4.3 × 10-9). SNP rs6902488 at 6p25.2 (r2 = 0.97 with rs7748275) lies in an active enhancer and is predicted to impact binding of STAT3, P300 and ELF1. We identified additional risk loci with low Bayesian false discovery probability (BFDP) scores, indicating they are likely to be true risk associations (BFDP <10%). At chromosome 20q11.22, rs74272064 was associated with HGSOC risk (OR = 1.27, P = 9.0 × 10-8). Overall EOC risk was associated with rs10260419 at chromosome 7p21.3 (OR = 1.33, P = 1.2 × 10-7) and rs74917072 at chromosome 2q37.3 (OR = 1.25, P = 4.7 × 10-7). At 2q37.3, expression quantitative trait locus analysis in 404 HGSOC tissues identified ESPNL as a putative candidate susceptibility gene (P = 1.2 × 10-7).

CONCLUSION

While some risk loci were shared between East Asian and European populations, others were population-specific, indicating that the landscape of EOC risk in Asian women has both shared and unique features compared to women of European ancestry.

Genome-wide association studies and CRISPR/Cas9-mediated gene editing identify regulatory variants influencing eyebrow thickness in humans

Hair plays an important role in primates and is clearly subject to adaptive selection. While humans have lost most facial hair, eyebrows are a notable exception. Eyebrow thickness is heritable and widely believed to be subject to sexual selection. Nevertheless, few genomic studies have explored its genetic basis. Here, we performed a genome-wide scan for eyebrow thickness in 2961 Han Chinese. We identified two new loci of genome-wide significance, at 3q26.33 near SOX2 (rs1345417: P = 6.51×10(-10)) and at 5q13.2 near FOXD1 (rs12651896: P = 1.73×10(-8)). We further replicated our findings in the Uyghurs, a population from China characterized by East Asian-European admixture (N = 721), the CANDELA cohort from five Latin American countries (N = 2301), and the Rotterdam Study cohort of Dutch Europeans (N = 4411). A meta-analysis combining the full GWAS results from the three cohorts of full or partial Asian descent (Han Chinese, Uyghur and Latin Americans, N = 5983) highlighted a third signal of genome-wide significance at 2q12.3 (rs1866188: P = 5.81×10(-11)) near EDAR. We performed fine-mapping and prioritized four variants for further experimental verification. CRISPR/Cas9-mediated gene editing provided evidence that rs1345417 and rs12651896 affect the transcriptional activity of the nearby SOX2 and FOXD1 genes, which are both involved in hair development. Finally, suitable statistical analyses revealed that none of the associated variants showed clear signals of selection in any of the populations tested. Contrary to popular speculation, we found no evidence that eyebrow thickness is subject to strong selective pressure.

HLA Class I and II Diversity Contributes to the Etiologic Heterogeneity of Non-Hodgkin Lymphoma Subtypes

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.

Associations between genetic variants in immunoregulatory genes and risk of non-Hodgkin lymphoma in a Chinese population

We undertook a hospital-based case-control study to examine the associations between single nucleotide polymorphisms (SNPs) in selected immunoregulatory genes and non-Hodgkin lymphoma (NHL) risk in a Chinese population. One hundred and sixty-nine NHL patients diagnosed according to the World Health Organization (WHO) 2001 standard and 421 controls were recruited. Nine SNPs in three genes (IL-10, IL-1RN, and TNF-α) were selected based on predicted functions and previous study findings. Genetic association analysis was performed using the Cochran-Armitage trend test and multiple logistic regression. Four SNPs were associated with an increased risk of overall NHL: odds ratio per minor allele [ORper-minor-allele] and 95% confidence interval [CI] were 2.64 (1.75-3.98) for IL-10 rs1800893, 2.67 (1.72-4.16) for IL-1RN rs4251961, 1.80 (1.24-2.63) for TNF- α rs1800630, and 1.55 (1.02-2.37) for TNF- α rs2229094. These SNPs were also associated with an increased risk of diffuse large B-cell lymphoma (DLBCL). In addition, another SNP (TNF- α rs1041981) was associated with an increased risk of DLBCL (ORper-minor-allele=1.73, 95% CI 1.14-2.61). The findings provide evidence on the role of these immunoregulatory gene variants in NHL etiology.

Germline variation in the 3'-untranslated region of the POU2AF1 gene is associated with susceptibility to lymphoma

Genetic variations in certain genes may alter the susceptibility to lymphoma. We searched electronic databases and selected candidate single-nucleotide polymorphisms (SNPs) located within 3'-untranslated regions (3'-UTRs) that might affect miRNA-binding ability in the 50 most dysregulated genes in lymphoma for further study. We found that rs1042752-located in the 3'-UTR of POU2AF1, which plays a vital role in lymphomagenesis-was significantly associated with lymphoma risk in a case-control study with 793 patients and 939 controls. Compared with individuals with the rs1042752TT genotype, those with the rs1042752CC genotype had a higher risk of lymphoma (OR = 2.14, 95% CI: 1.55-2.95, P < 0.001), even in stratified analysis for non-Hodgkin lymphoma (OR = 4.58, 95% CI: 2.38-8.81, P < 0.001), B-cell lymphoma (OR = 4.89, 95% CI: 2.46-9.73, P < 0.001), T-cell lymphoma (OR = 4.20, 95% CI: 1.76-10.05, P = 0.001), and Hodgkin lymphoma (OR = 3.62, 95% CI: 1.25-10.46, P = 0.018). Similar results were also observed in a recessive genetic model. Database findings suggested that rs1042752 might affect the interaction of POU2AF1 mRNA with hsa-miR-633. Functional assays confirmed that rs1042752C altered the binding site of hsa-miR-633 and increased POU2AF1 expression in Ramos, HuT 102, and Jurkat E6-1 cell lines. These findings demonstrate for the first time that functional polymorphism in the 3'-UTR of POU2AF1 is associated with susceptibility, and that SNP interaction with hsa-miR-633 affects gene expression and increases the risk of lymphoma.

Possible Interaction Between Cigarette Smoking and HLA-DRB1 Variation in the Risk of Follicular Lymphoma

Follicular lymphoma (FL) risk is strongly associated with germline genetic variation in human leukocyte antigen (HLA) class II. Cigarette smoking has been suggested to increase FL risk, primarily among women. We hypothesized that amino acids in HLA-antigen D-related β1 subunit (DRB1) interact with smoking in FL risk, as shown for rheumatoid arthritis. We analyzed 373 patients with FL and 818 controls from 2 population-based case-control studies in Sweden and Denmark (1999-2003). Haplotypes in HLA-DRB1 were imputed at amino acid positions 11, 13, 28, 30, and 70-74 (shared epitope). We estimated the relative risk of FL as odds ratios with 95% confidence intervals for different smoking status/haplotype combinations. Interaction was defined as departure from additivity of effects and quantified by the attributable proportion (AP). Relative to never-smokers carrying no shared epitope alleles, smoking was associated with the risk of FL among all subjects (for former smokers, odds ratio (OR) = 2.20, 95% confidence interval (CI): 1.10, 4.41; ORcurrent = 3.56, 95% CI: 1.60, 7.92) and women (ORformer = 2.95, 95% CI: 1.18, 7.37; ORcurrent = 5.63, 95% CI: 2.07, 15.3) carrying 2 shared epitope alleles but not among those carrying zero or 1 shared epitope allele. Smoking and shared epitope status interacted significantly as measured by AP (overall, AP = 0.6, 95% CI: 0.15, 1.0; for women, AP = 0.5, 95% CI: 0.005, 1.0). These results suggest a possible interaction between smoking and HLA-DRB1-associated antigen presentation in FL risk and provide a model to further unravel FL etiology.

SNP variants associated with non-Hodgkin lymphoma (NHL) correlate with human leukocyte antigen (HLA) class II expression

Large consortia efforts and genome-wide association studies (GWASs) have linked a number of genetic variants within the 6p21 chromosomal region to non-Hodgkin lymphoma (NHL). Complementing these efforts, we genotyped previously reported SNPs in the human leukocyte antigen (HLA) class I (rs6457327) and class II (rs9271100, rs2647012 and rs10484561) regions in a total of 1,145 subjects (567 NHL cases and 578 healthy controls) from two major ethnic groups in Malaysia, the Malays and the Chinese. We identified a NHL-associated (P_{NHL_add} = 0.0008; OR_{NHL_add} = 0.54; 95% CI = 0.37–0.77) and B-cell associated (P_{Bcell_add} = 0.0007; OR_{Bcell_add} = 0.51; 95% CI = 0.35–0.76) SNP rs2647012 in the Malaysian Malays. In silico cis-eQTL analysis of rs2647012 suggests potential regulatory function of nearby HLA class II molecules. Minor allele rs2647012-T is linked to higher expression of HLA-DQB1, rendering a protective effect to NHL risk. Our findings suggest that the HLA class II region plays an important role in NHL etiology.

A Genome-Wide SNP Linkage Analysis Suggests a Susceptibility Locus on 6p21 for Ankylosing Spondylitis and Inflammatory Back Pain Trait

Objectives

To screen susceptibility loci for ankylosing spondylitis (AS) using an affected-only linkage analysis based on high-density single nucleotide polymorphisms (SNPs) in a genome-wide manner.

Patients and Methods

AS patients from ten families with Cantonese origin of China were enrolled in the study. Blood samples were genotyped using genomic DNA derived from peripheral blood leukocytes by Illumina HumanHap 610-Quad SNP Chip. Genotype data were generated using the Illumina BeadStudio 3.2 software. PLINK package was used to remove non-autosomal SNPs and to further eliminate markers of typing errors. An affected-only linkage analysis was carried out using both non-parametric and parametric linkage analyses, as implemented in MERLIN.

Result

Seventy-eight AS patients (48 males and 30 females, mean age: 39±16 years) were enrolled in the study. The mean age of onset was 23±10 years and mean duration of disease was 16.7±12.2 years. Iritis (2/76, 2.86%), dactylitis (5/78, 6.41%), hip joint involvement (9/78, 11.54%), peripheral arthritis (22/78, 28.21%), inflammatory back pain (IBP) (69/78, 88.46%) and HLA-B27 positivity (70/78, 89.74%) were observed in these patients. Using non-parameter linkage analysis, we found one susceptibility locus for AS, IBP and HLA-B27 in 6p21 respectively, spanning about 13.5Mb, 20.9Mb and 21.2Mb, respectively No significant results were found in the other clinical trait groups including dactylitis, hip involved and arthritis. The identical susceptibility locus region spanning above 9.44Mb was detected in AS IBP and HLA-B27 by the parametric linkage analysis.

Conclusion

Our genome-wide SNP linkage analysis in ten families with ankylosing spondylitis suggests a susceptibility locus on 6p21 in AS, which is a risk locus for IBP in AS patients.

HLA and Risk of Diffuse Large B cell Lymphoma After Solid Organ Transplantation

BACKGROUND

Solid organ transplant recipients have heightened risk for diffuse large B cell lymphoma (DLBCL). The role of donor-recipient HLA mismatch and recipient HLA type on DLBCL risk are not well established.

METHODS

We examined 172 231 kidney, heart, pancreas, and lung recipients transplanted in the United States between 1987 and 2010, including 902 with DLBCL. Incidence rate ratios (IRRs) were calculated using Poisson regression for DLBCL risk in relation to HLA mismatch, types, and zygosity, adjusting for sex, age, race/ethnicity, year, organ, and transplant number.

RESULTS

Compared with recipients who had 2 HLA-DR mismatches, those with zero or 1 mismatch had reduced DLBCL risk, (zero: IRR, 0.76, 95% confidence interval [95% CI], 0.61-0.95; one: IRR, 0.83; 95% CI, 0.69-1.00). In stratified analyses, recipients matched at either HLA-A, -B, or -DR had a significantly reduced risk of late-onset (>2 years after transplantation), but not early-onset DLBCL, and there was a trend for decreasing risk with decreasing mismatch across all 3 loci (P = 0.0003). Several individual recipient HLA-A, -B, -C, -DR, and -DQ antigens were also associated with DLBCL risk, including DR13 (IRR, 0.74; 95% CI, 0.57-0.93) and B38 (IRR, 1.48; 95% CI, 1.10-1.93), confirming prior findings that these 2 antigens are associated with risk of infection-associated cancers.

CONCLUSIONS

In conclusion, variation in HLA is related to susceptibility to DLBCL, perhaps reflecting intensity of immunosuppression, control of Epstein-Barr virus infection among transplant recipients or chronic immune stimulation.

Comprehensive genomic profiling of orbital and ocular adnexal lymphomas identifies frequent alterations in MYD88 and chromatin modifiers: new routes to targeted therapies

Non-Hodgkin lymphoma of the orbit and ocular adnexa is the most common primary orbital malignancy. Treatments for low- (extra-nodal marginal zone and follicular lymphomas) and high-grade (diffuse large B-cell lymphoma) are associated with local and vision-threatening toxicities. High-grade lymphomas relapse frequently and exhibit poor survival rates. Despite advances in genomic profiling and precision medicine, orbital and ocular adnexal lymphomas remain poorly characterized molecularly. We performed targeted next-generation sequencing (NGS) profiling of 38 formalin-fixed, paraffin-embedded orbital and ocular adnexal lymphomas obtained from a single-center using a panel targeting near-term, clinically relevant genes. Potentially actionable mutations and copy number alterations were prioritized based on gain- and loss-of-function analyses, and catalogued, approved, and investigational therapies. Of 36 informative samples, including marginal zone lymphomas (n=20), follicular lymphomas (n=9), and diffuse large B-cell lymphomas (n=7), 53% harbored a prioritized alteration (median=1, range 0-5/sample). MYD88 was the most frequently altered gene in our cohort, with potentially clinically relevant hotspot gain-of-function mutations identified in 71% of diffuse large B-cell lymphomas and 25% of marginal zone lymphomas. Prioritized alterations in epigenetic modulators were common and included gain-of-function EZH2 and loss-of-function ARID1A mutations (14% of diffuse large B-cell lymphomas and 22% of follicular lymphomas contained alterations in each of these two genes). Single prioritized alterations were also identified in the histone methyltransferases KMT2B (follicular lymphoma) and KMT3B (diffuse large B-cell lymphoma). Loss-of-function mutations and copy number alterations in the tumor suppressors TP53 (diffuse large B-cell and follicular lymphoma), CDKN2A (diffuse large B-cell and marginal zone lymphoma), PTEN (diffuse large B-cell lymphoma), ATM (diffuse large B-cell lymphoma), and NF1 (diffuse large B-cell lymphoma), and gain-of-function mutations in the oncogenes HRAS (follicular lymphoma) and NRAS (diffuse large B-cell lymphoma) were also observed. Together, our study demonstrates that NGS can be used to profile routine formalin-fixed, paraffin-embedded orbital and ocular adnexal lymphomas for identification of somatic-driving alterations and nomination of potential therapeutic strategies.

Emerging therapies provide new opportunities to reshape the multifaceted interactions between the immune system and lymphoma cells

The acquisition of a complete neoplastic phenotype requires cancer cells to develop escape mechanisms from the host immune system. This phenomenon, commonly referred to as 'immune evasion,' represents a hallmark of cancers and results from a Darwinian selection of the fittest tumor clones. First reported in solid tumors, cancer immunoescape characterizes several hematological malignancies. The biological bases of cancer immunoescape have recently been disclosed and include: (i) impaired human leukocyte antigen-mediated cancer cell recognition (B2M, CD58, CTIIA, CD80/CD86, CD28 and CTLA-4 mutations); (ii) deranged apoptotic mechanisms (reduced pro-apoptotic signals and/or increased expression of anti-apoptotic molecules); and (iii) changes in the tumor microenvironment involving regulatory T cells and tumor-associated macrophages. These immune-escape mechanisms characterize both Hodgkin and non-Hodgkin (B and T cell) lymphomas and represent a promising target for new anti-tumor therapies. In the present review, the principles of cancer immunoescape and their role in human lymphomagenesis are illustrated. Current therapies targeting these pathways and possible applications for lymphoma treatment are also addressed.

Human leukocyte antigens and genetic susceptibility to lymphoma

Familial aggregation, coupled with ethnic variation in incidence, suggests that inherited susceptibility plays a role in the development of lymphoma, and the search for genetic risk factors has highlighted the contribution of the human leukocyte antigen (HLA) complex. In a landmark study published almost 50 years ago, Hodgkin lymphoma (HL) was the first disease to be associated with HLA variation. It is now clear that Epstein-Barr virus (EBV)-positive and -negative HL are strongly associated with specific HLA polymorphisms but these differ by EBV status of the tumours. HLA class I alleles are consistently associated with EBV-positive HL while a polymorphism in HLA class II is the strongest predictor of risk of EBV-negative HL. Recent investigations, particularly genome-wide association studies (GWAS), have also revealed associations between HLA and common types of non-Hodgkin lymphoma (NHL). Follicular lymphoma is strongly associated with two distinct haplotypes in HLA class II whereas diffuse large B-cell lymphoma is most strongly associated with HLA-B*08. Although chronic lymphocytic leukaemia is associated with variation in HLA class II, the strongest signals in GWAS are from non-HLA polymorphisms, suggesting that inherited susceptibility is explained by co-inheritance of multiple low risk variants. Associations between B-cell derived lymphoma and HLA variation suggest that antigen presentation, or lack of, plays an important role in disease pathogenesis but the precise mechanisms have yet to be elucidated.

CD19 chimeric antigen receptor (CD19 CAR)-redirected adoptive T-cell immunotherapy for the treatment of relapsed or refractory B-cell Non-Hodgkin's Lymphomas

Recovery rates for B-cell Non-Hodgkin's Lymphoma (NHL) are up to 70% with current standard-of-care treatments including rituximab (chimeric anti-CD20 monoclonal antibody) in combination with chemotherapy (R-CHOP). However, patients who do not respond to first-line treatment or develop resistance have a very poor prognosis. This signifies the need for the development of an optimal treatment approach for relapsed/refractory B-NHL. Novel CD19- chimeric antigen receptor (CAR) T-cell redirected immunotherapy is an attractive option for this subset of patients. Anti-CD19 CAR T-cell therapy has already had remarkable efficacy in various leukemias as well as encouraging outcomes in phase I clinical trials of relapsed/refractory NHL. In going forward with additional clinical trials, complementary treatments that may circumvent potential resistance mechanisms should be used alongside anti-CD19 T-cells in order to prevent relapse with resistant strains of disease. Some such supplementary tactics include conditioning with lymphodepletion agents, sensitizing with kinase inhibitors and Bcl-2 inhibitors, enhancing function with multispecific CAR T-cells and CD40 ligand-expressing CAR T-cells, and safeguarding with lymphoma stem cell-targeted treatments. A therapy regimen involving anti-CD19 CAR T-cells and one or more auxiliary treatments could dramatically improve prognoses for patients with relapsed/refractory B-cell NHL. This approach has the potential to revolutionize B-NHL salvage therapy in much the same way rituximab did for first-line treatments.

Inherited Inflammatory Response Genes Are Associated with B-Cell Non-Hodgkin's Lymphoma Risk and Survival

BACKGROUND

Malignant B-cell clones are affected by both acquired genetic alterations and by inherited genetic variations changing the inflammatory tumour microenvironment.

METHODS

We investigated 50 inflammatory response gene polymorphisms in 355 B-cell non-Hodgkin's lymphoma (B-NHL) samples encompassing 216 diffuse large B cell lymphoma (DLBCL) and 139 follicular lymphoma (FL) and 307 controls. The effect of single genes and haplotypes were investigated and gene-expression analysis was applied for selected genes. Since interaction between risk genes can have a large impact on phenotype, two-way gene-gene interaction analysis was included.

RESULTS

We found inherited SNPs in genes critical for inflammatory pathways; TLR9, IL4, TAP2, IL2RA, FCGR2A, TNFA, IL10RB, GALNT12, IL12A and IL1B were significantly associated with disease risk and SELE, IL1RN, TNFA, TAP2, MBL2, IL5, CX3CR1, CHI3L1 and IL12A were, associated with overall survival (OS) in specific diagnostic entities of B-NHL. We discovered noteworthy interactions between DLBCL risk alleles on IL10 and IL4RA and FL risk alleles on IL4RA and IL4. In relation to OS, a highly significant interaction was observed in DLBCL for IL4RA (rs1805010) * IL10 (rs1800890) (HR = 0.11 (0.02-0.50)). Finally, we explored the expression of risk genes from the gene-gene interaction analysis in normal B-cell subtypes showing a different expression of IL4RA, IL10, IL10RB genes supporting a pathogenetic effect of these interactions in the germinal center.

CONCLUSIONS

The present findings support the importance of inflammatory genes in B-cell lymphomas. We found association between polymorphic sites in inflammatory response genes and risk as well as outcome in B-NHL and suggest an effect of gene-gene interactions during the stepwise oncogenesis.

Tissue-Specific Enrichment of Lymphoma Risk Loci in Regulatory Elements

Though numerous polymorphisms have been associated with risk of developing lymphoma, how these variants function to promote tumorigenesis is poorly understood. Here, we report that lymphoma risk SNPs, especially in the non-Hodgkin's lymphoma subtype chronic lymphocytic leukemia, are significantly enriched for co-localization with epigenetic marks of active gene regulation. These enrichments were seen in a lymphoid-specific manner for numerous ENCODE datasets, including DNase-hypersensitivity as well as multiple segmentation-defined enhancer regions. Furthermore, we identify putatively functional SNPs that are both in regulatory elements in lymphocytes and are associated with gene expression changes in blood. We developed an algorithm, UES, that uses a Monte Carlo simulation approach to calculate the enrichment of previously identified risk SNPs in various functional elements. This multiscale approach integrating multiple datasets helps disentangle the underlying biology of lymphoma, and more broadly, is generally applicable to GWAS results from other diseases as well.

Familial predisposition and genetic risk factors for lymphoma

Our understanding of familial predisposition to lymphoma (collectively defined as non-Hodgkin lymphoma [NHL], Hodgkin lymphoma [HL], and chronic lymphocytic leukemia [CLL]) outside of rare hereditary syndromes has progressed rapidly during the last decade. First-degree relatives of NHL, HL, and CLL patients have an ∼1.7-fold, 3.1-fold, and 8.5-fold elevated risk of developing NHL, HL, and CLL, respectively. These familial risks are elevated for multiple lymphoma subtypes and do not appear to be confounded by nongenetic risk factors, suggesting at least some shared genetic etiology across the lymphoma subtypes. However, a family history of a specific subtype is most strongly associated with risk for that subtype, supporting subtype-specific genetic factors. Although candidate gene studies have had limited success in identifying susceptibility loci, genome-wide association studies (GWAS) have successfully identified 67 single nucleotide polymorphisms from 41 loci, predominately associated with specific subtypes. In general, these GWAS-discovered loci are common (minor allele frequency >5%), have small effect sizes (odds ratios, 0.60-2.0), and are of largely unknown function. The relatively low incidence of lymphoma, modest familial risk, and the lack of a screening test and associated intervention, all argue against active clinical surveillance for lymphoma in affected families at this time.

To Each Its Own: Linking the Biology and Epidemiology of NHL Subtypes

Non-Hodgkin lymphoma (NHL) constitutes a diverse group of more than 40 subtypes, each characterized by distinct biologic and clinical features. Until recently, pinpointing genetic and epidemiologic risk factors for individual subtypes has been limited by the relative rarity of each. However, several large pooled case-control studies have provided sufficient statistical power for detecting etiologic differences and commonalities between subtypes and thus yield new insight into their unique epidemiologic backgrounds. Here, we review the subtype-specific medical, lifestyle, and biologic components identified in these studies, which suggest that a complex interplay between host genetics, autoimmune disorders, modifiable risk factors, and occupation contributes to lymphomagenesis.

Inherited variation in immune response genes in follicular lymphoma and diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL) both depend on immune-mediated survival and proliferation signals from the tumor microenvironment. Inherited genetic variation influences this complex interaction. A total of 89 studies investigating immune-response genes in DLBCL and FL were critically reviewed. Relatively consistent association exists for variation in the tumor necrosis factor alpha (TNFA) and interleukin-10 loci and DLBCL risk; for DLBCL outcome association with the TNFA locus exists. Variations at chromosome 6p31-32 were associated with FL risk. Importantly, individual risk alleles have been shown to interact with each other. We suggest that the pathogenetic impact of polymorphic genes should include gene-gene interaction analysis and should be validated in preclinical model systems of normal B lymphopoiesis and B-cell malignancies. In the future, large cohort studies of interactions and genome-wide association studies are needed to extend the present findings and explore new risk alleles to be studied in preclinical models.

Rationale and Design of the International Lymphoma Epidemiology Consortium (InterLymph) Non-Hodgkin Lymphoma Subtypes Project

BACKGROUND

Non-Hodgkin lymphoma (NHL), the most common hematologic malignancy, consists of numerous subtypes. The etiology of NHL is incompletely understood, and increasing evidence suggests that risk factors may vary by NHL subtype. However, small numbers of cases have made investigation of subtype-specific risks challenging. The International Lymphoma Epidemiology Consortium therefore undertook the NHL Subtypes Project, an international collaborative effort to investigate the etiologies of NHL subtypes. This article describes in detail the project rationale and design.

METHODS

We pooled individual-level data from 20 case-control studies (17471 NHL cases, 23096 controls) from North America, Europe, and Australia. Centralized data harmonization and analysis ensured standardized definitions and approaches, with rigorous quality control.

RESULTS

The pooled study population included 11 specified NHL subtypes with more than 100 cases: diffuse large B-cell lymphoma (N = 4667), follicular lymphoma (N = 3530), chronic lymphocytic leukemia/small lymphocytic lymphoma (N = 2440), marginal zone lymphoma (N = 1052), peripheral T-cell lymphoma (N = 584), mantle cell lymphoma (N = 557), lymphoplasmacytic lymphoma/Waldenström macroglobulinemia (N = 374), mycosis fungoides/Sézary syndrome (N = 324), Burkitt/Burkitt-like lymphoma/leukemia (N = 295), hairy cell leukemia (N = 154), and acute lymphoblastic leukemia/lymphoma (N = 152). Associations with medical history, family history, lifestyle factors, and occupation for each of these 11 subtypes are presented in separate articles in this issue, with a final article quantitatively comparing risk factor patterns among subtypes.

CONCLUSIONS

The International Lymphoma Epidemiology Consortium NHL Subtypes Project provides the largest and most comprehensive investigation of potential risk factors for a broad range of common and rare NHL subtypes to date. The analyses contribute to our understanding of the multifactorial nature of NHL subtype etiologies, motivate hypothesis-driven prospective investigations, provide clues for prevention, and exemplify the benefits of international consortial collaboration in cancer epidemiology.

Associations of non-Hodgkin Lymphoma (NHL) risk with autoimmune conditions according to putative NHL loci

Autoimmune conditions and immune system-related genetic variations are associated with risk of non-Hodgkin lymphoma (NHL). In a pooled analysis of 8,692 NHL cases and 9,260 controls from 14 studies (1988-2007) within the International Lymphoma Epidemiology Consortium, we evaluated the interaction between immune system genetic variants and autoimmune conditions in NHL risk. We evaluated the immunity-related single nucleotide polymorphisms rs1800629 (tumor necrosis factor gene (TNF) G308A), rs1800890 (interleukin-10 gene (IL10) T3575A), rs6457327 (human leukocyte antigen gene (HLA) class I), rs10484561 (HLA class II), and rs2647012 (HLA class II)) and categorized autoimmune conditions as primarily mediated by B-cell or T-cell responses. We constructed unconditional logistic regression models to measure associations between autoimmune conditions and NHL with stratification by genotype. Autoimmune conditions mediated by B-cell responses were associated with increased NHL risk, specifically diffuse large B-cell lymphoma (odds ratio (OR) = 3.11, 95% confidence interval (CI): 2.25, 4.30) and marginal zone lymphoma (OR = 5.80, 95% CI: 3.82, 8.80); those mediated by T-cell responses were associated with peripheral T-cell lymphoma (OR = 2.14, 95% CI: 1.35, 3.38). In the presence of the rs1800629 AG/AA genotype, B-cell-mediated autoimmune conditions increased NHL risk (OR = 3.27, 95% CI: 2.07, 5.16; P-interaction = 0.03) in comparison with the GG genotype (OR = 1.82, 95% CI: 1.31, 2.53). This interaction was consistent across major B-cell NHL subtypes, including marginal zone lymphoma (P-interaction = 0.02) and follicular lymphoma (P-interaction = 0.04).

Genome-wide association study identifies shared risk loci common to two malignancies in golden retrievers

Dogs, with their breed-determined limited genetic background, are great models of human disease including cancer. Canine B-cell lymphoma and hemangiosarcoma are both malignancies of the hematologic system that are clinically and histologically similar to human B-cell non-Hodgkin lymphoma and angiosarcoma, respectively. Golden retrievers in the US show significantly elevated lifetime risk for both B-cell lymphoma (6%) and hemangiosarcoma (20%). We conducted genome-wide association studies for hemangiosarcoma and B-cell lymphoma, identifying two shared predisposing loci. The two associated loci are located on chromosome 5, and together contribute ~20% of the risk of developing these cancers. Genome-wide p-values for the top SNP of each locus are 4.6×10^-7 and 2.7×10^-6, respectively. Whole genome resequencing of nine cases and controls followed by genotyping and detailed analysis identified three shared and one B-cell lymphoma specific risk haplotypes within the two loci, but no coding changes were associated with the risk haplotypes. Gene expression analysis of B-cell lymphoma tumors revealed that carrying the risk haplotypes at the first locus is associated with down-regulation of several nearby genes including the proximal gene TRPC6, a transient receptor Ca²⁺-channel involved in T-cell activation, among other functions. The shared risk haplotype in the second locus overlaps the vesicle transport and release gene STX8. Carrying the shared risk haplotype is associated with gene expression changes of 100 genes enriched for pathways involved in immune cell activation. Thus, the predisposing germ-line mutations in B-cell lymphoma and hemangiosarcoma appear to be regulatory, and affect pathways involved in T-cell mediated immune response in the tumor. This suggests that the interaction between the immune system and malignant cells plays a common role in the tumorigenesis of these relatively different cancers.

To shed light on the genetic predisposition to cancers of the hematologic system, we performed genome-wide association analysis of affected and non-affected pet dogs. Dogs naturally develop the same diseases as humans, including cancer, and the relatively limited genetic diversity within different breeds makes genetic studies easier compared to in humans. By doing genome-wide association, we identified loci predisposing to hemangiosarcoma and B-cell lymphoma. To our surprise, we found two shared loci predisposing to both diseases. Within these two regions we identified several partially overlapping haplotypes, predisposing somewhat differently to the two cancers. We found no coding mutations that followed the risk or non-risk haplotypes suggesting that regulatory mutations exert the effect on disease. We also looked at gene expression in B-cell lymphomas, comparing samples from individuals with risk or non-risk haplotypes. This analysis showed differential expression associated with the haplotypes at both loci, suggesting the risk haplotypes are associated with an effect on T-cell response.

A genome-wide association study of marginal zone lymphoma shows association to the HLA region

Marginal zone lymphoma (MZL) is the third most common subtype of B-cell non-Hodgkin lymphoma. Here we perform a two-stage GWAS of 1,281 MZL cases and 7,127 controls of European ancestry and identify two independent loci near BTNL2 (rs9461741, P=3.95 × 10(-15)) and HLA-B (rs2922994, P=2.43 × 10(-9)) in the HLA region significantly associated with MZL risk. This is the first evidence that genetic variation in the major histocompatibility complex influences MZL susceptibility.

Epidemiology and etiology of non-hodgkin lymphoma

Non-Hodgkin lymphoma (NHL) consists of many histologically and biologically distinct lymphoid malignancies with poorly understood, but possibly distinct, etiologies. The patterns of incidence and time trend vary not only by age, sex, and race/ethnicity in the USA, but also show significant geographic differences, suggesting the potential role of infectious agents, environmental factors, and lifestyle factors in addition to host genetic status in the development of NHL. Important pathogenetic mechanisms include immune modulation and chronic antigen stimulation. Epidemiologic studies in the past two decades have provided intriguing new insights on the possible causes of lymphoma and support the idea that there is some mechanistic commonality of lymphomagenesis, but significant etiologic heterogeneity clearly exists. This review presents a summary of the current understanding of the descriptive epidemiology and etiology of NHL and suggests areas of focus for future epidemiologic research.

The association between the polymorphisms of TNF-α and non-Hodgkin lymphoma: a meta-analysis

Many genetic variations in the promoter region of tumor necrosis factor alpha (TNF-α) may confer host susceptibility to cancer by influencing TNF-α expression. Nevertheless, the results remain inconclusive. The current meta-analysis was performed to investigate the association between three common TNF-α promoter polymorphisms and the risk of non-Hodgkin lymphoma (NHL). A literature search was conducted mainly from PubMed for all eligible studies. The pooled odds ratios (ORs) and corresponding 95 % confidence intervals (CIs) were used to assess the association of TNF-α polymorphisms with the risk of NHL. TNF-α -308 A allele showed a statistically significant increased risk for NHL under the homozygous (AA vs. GG, OR = 1.51, 95 % CI = 1.26-1.80) and recessive (OR = 1.47, 95 % CI = 1.23-1.75) models, respectively. The stratified analyses showed an increased risk of NHL with the presence of TNF-α -308 A allele among Africans and Caucasians, but a decreased risk among Asians. No association was observed between -238 G/A polymorphism and NHL risk either in the overall analysis or in the stratified analysis. Similarly, pooled analysis did not reveal an altered risk of NHL with -857 C/T polymorphism. Nonetheless, a statistically significant association was observed among Asians when stratified by ethnicity. Among the three genetic variations of interest, TNF-α -308 G/A polymorphism was significantly associated with the risk of NHL; neither -238 G/A nor -857 C/T polymorphism was shown to alter the overall NHL risk; however, stratified analysis by ethnicity observed a statistically significant association between -857 C/T polymorphism and the risk of NHL among Asians.

Genome-wide association study identifies multiple susceptibility loci for diffuse large B cell lymphoma

Diffuse large B cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of 3 new genome-wide association studies (GWAS) and 1 previous scan, totaling 3,857 cases and 7,666 controls of European ancestry, with additional genotyping of 9 promising SNPs in 1,359 cases and 4,557 controls. In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; P = 2.33 × 10(-21)), rs2523607 at 6p21.33 (HLA-B; P = 2.40 × 10(-10)), rs79480871 at 2p23.3 (NCOA1; P = 4.23 × 10(-8)) and two independent SNPs, rs13255292 and rs4733601, at 8q24.21 (PVT1; P = 9.98 × 10(-13) and 3.63 × 10(-11), respectively). These data provide substantial new evidence for genetic susceptibility to this B cell malignancy and point to pathways involved in immune recognition and immune function in the pathogenesis of DLBCL.

Discovery of new risk loci for IgA nephropathy implicates genes involved in immunity against intestinal pathogens

We performed a genome-wide association study (GWAS) of IgA nephropathy (IgAN), the most common form of glomerulonephritis, with discovery and follow-up in 20,612 individuals of European and East Asian ancestry. We identified six novel genome-wide significant associations, four in ITGAM-ITGAX, VAV3 and CARD9 and two new independent signals at HLA-DQB1 and DEFA. We replicated the nine previously reported signals, including known SNPs in the HLA-DQB1 and DEFA loci. The cumulative burden of risk alleles is strongly associated with age at disease onset. Most loci are either directly associated with risk of inflammatory bowel disease (IBD) or maintenance of the intestinal epithelial barrier and response to mucosal pathogens. The geo-spatial distribution of risk alleles is highly suggestive of multi-locus adaptation and the genetic risk correlates strongly with variation in local pathogens, particularly helminth diversity, suggesting a possible role for host-intestinal pathogen interactions in shaping the genetic landscape of IgAN.

A comprehensive evaluation of the role of genetic variation in follicular lymphoma survival

Background

Survival in follicular lymphoma (FL) is highly variable, even within prognostic groups defined by tumor grade and the Follicular Lymphoma International Prognostic Index. Studies suggest that germline single nucleotide polymorphisms (SNPs) may hold prognostic information but further investigation is needed.

Methods

We explored the association between SNPs and FL outcome using two approaches: 1) Two independent genome-wide association studies (GWAS) of ~300.000 SNPs followed by a meta-analysis encompassing 586 FL patients diagnosed in Denmark/Sweden 1999–2002 and in the United States 2001–2006; and 2) Investigation of 22 candidate-gene variants previously associated with FL outcome in the Danish/Swedish cohort (N = 373). We estimated time to lymphoma-specific death (approach 1 and 2) and lymphoma progression (approach 2) with hazard ratios (HR) and 95% confidence intervals (CI) in a multivariable Cox regression model.

Results

In the GWAS meta-analysis, using a random effects model, no variants were associated with lymphoma-specific death at a genome-wide significant level (p < 5.0 ×10⁻⁸). The strongest association was observed for tightly linked SNPs on 17q24 near the ABCA10 and ABCA6 genes (rs10491178 HR_random = 3.17, 95% CI 2.09-4.79, p_random = 5.24 ×10⁻⁸). The ABCA10 and ABCA6 genes belong to a family of genes encoding for ABC transporter proteins, implicated in multidrug resistance. In line with a previous study, rs2466571 in CD46 (HR = 0.73, 95% CI 0.58-0.91, p = 0.006) showed nominal association with lymphoma progression, as did two highly linked SNPs in IL8 (rs4073 HR = 0.78, 95% CI 0.62-0.97, p = 0.02; rs2227307 HR = 0.75, 95% CI 0.60-0.94, p = 0.01) previously associated with overall survival.

Conclusions

The results suggest a possible role for multidrug resistance in FL survival and add to the evidence that genetic variation in CD46 and IL8 may have prognostic implications in FL. Our findings need further confirmation in other independent populations or in a larger multicenter GWAS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-014-0113-6) contains supplementary material, which is available to authorized users.

Genome-wide association study identifies five susceptibility loci for follicular lymphoma outside the HLA region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10(-20)) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10(-11)) near ETS1; 3q28 (rs6444305, p = 1.10 × 10(-10)) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10(-10)) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10(-8)) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (pomnibus = 4.20 × 10(-67) to 2.67 × 10(-70)). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR(per-allele)] = 1.44; p = 4.59 × 10(-16)) and rs3130437 in HLA class I (OR(per-allele) = 1.23; p = 8.23 × 10(-9)). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.

Copy number variation analysis on a non-Hodgkin lymphoma case-control study identifies an 11q25 duplication associated with diffuse large B-cell lymphoma

Recent GWAS have identified several susceptibility loci for NHL. Despite these successes, much of the heritable variation in NHL risk remains to be explained. Common copy-number variants are important genomic sources of variability, and hence a potential source to explain part of this missing heritability. In this study, we carried out a CNV analysis using GWAS data from 681 NHL cases and 749 controls to explore the relationship between common structural variation and lymphoma susceptibility. Here we found a novel association with diffuse large B-cell lymphoma (DLBCL) risk involving a partial duplication of the C-terminus region of the LOC283177 long non-coding RNA that was further confirmed by quantitative PCR. For chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL), known somatic deletions were identified on chromosomes 13q14, 11q22-23, 14q32 and 22q11.22. Our study shows that GWAS data can be used to identify germline CNVs associated with disease risk for DLBCL and somatic CNVs for CLL/SLL.

Genetic variation in DNA repair pathways and risk of non-Hodgkin's lymphoma

Molecular and genetic evidence suggests that DNA repair pathways may contribute to lymphoma susceptibility. Several studies have examined the association of DNA repair genes with lymphoma risk, but the findings from these reports have been inconsistent. Here we provide the results of a focused analysis of genetic variation in DNA repair genes and their association with the risk of non-Hodgkin's lymphoma (NHL). With a population of 1,297 NHL cases and 1,946 controls, we have performed a two-stage case/control association analysis of 446 single nucleotide polymorphisms (SNPs) tagging the genetic variation in 81 DNA repair genes. We found the most significant association with NHL risk in the ATM locus for rs227060 (OR = 1.27, 95% CI: 1.13-1.43, p = 6.77×10(-5)), which remained significant after adjustment for multiple testing. In a subtype-specific analysis, associations were also observed for the ATM locus among both diffuse large B-cell lymphomas (DLBCL) and small lymphocytic lymphomas (SLL), however there was no association observed among follicular lymphomas (FL). In addition, our study provides suggestive evidence of an interaction between SNPs in MRE11A and NBS1 associated with NHL risk (OR = 0.51, 95% CI: 0.34-0.77, p = 0.0002). Finally, an imputation analysis using the 1,000 Genomes Project data combined with a functional prediction analysis revealed the presence of biologically relevant variants that correlate with the observed association signals. While the findings generated here warrant independent validation, the results of our large study suggest that ATM may be a novel locus associated with the risk of multiple subtypes of NHL.