- Research article
- Open Access
A comprehensive evaluation of the role of genetic variation in follicular lymphoma survival
BMC Medical Genetics volume 15, Article number: 113 (2014)
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.
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.
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-8). The strongest association was observed for tightly linked SNPs on 17q24 near the ABCA10 and ABCA6 genes (rs10491178 HRrandom = 3.17, 95% CI 2.09-4.79, prandom = 5.24 ×10-8). 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.
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.
For most patients with follicular lymphoma (FL) the disease is incurable and treatment aims to relieve symptoms and prolong survival . In the clinical setting, established prognostic markers include tumor grade and the Follicular Lymphoma International Prognostic Index (FLIPI) . However, marked variation in outcome remains within each risk group , and new prognostic markers are therefore desirable for more accurate prognostication and personalized treatment. In the search for such markers, several studies have assessed the impact of inherited genetic variation in FL survival -.
There is growing evidence that the tumor microenvironment and the host immune response, partly determined by host genetics, are important for the clinical course of FL . Germline genetic variation may further affect the metabolism  or receptor affinity of anti-cancer drugs  with possible implications for prognosis. To date, at least 30 specific germline single nucleotide polymorphisms (SNPs), mostly targeting cytokine and immune function genes (e.g., C9, CD46, CD55, CFH, IL2, IL4R, IL8 and FCGR2A), have been associated with FL outcome in previous studies ,,-,-. However, most studies were small in size (<150 FL cases ,,,,-), and the majority of the reported findings are unconfirmed. In recent genome-wide association studies (GWAS), an important role for variation in the human leukocyte antigen (HLA) region has been revealed for risk of FL, and one variant in the HLA class I region has been associated both with FL risk and survival -.
To gain additional knowledge of the role of genetic variation in FL outcome, we explored the association of ~300.000 SNPs with lymphoma-specific death in two independent GWAS in Sweden/Denmark  and the USA , followed by a full meta-analysis encompassing a total of 586 FL cases. We also investigated genetic variants previously reported to be associated with FL survival as well as recently established FL risk-associated SNPs for the association with lymphoma-specific death and progression in the Swedish/Danish (N = 373) and Swedish datasets (N = 231), respectively.
This study was approved by the Ethical Review Board in each country (Sweden: The Regional Ethical Review Board in Stockholm; Denmark: Scientific Ethics Committee for the Capital Region of Copenhagen; California, USA: UCSF Human Research Protection Program, Committee on Human Subjects Research) and all study participants gave informed consent.
SCALE study participants
We included cases with FL in the Scandinavian Lymphoma Etiology (SCALE) study, described in detail elsewhere . Briefly, SCALE is a population-based case-control study of the etiology of malignant lymphomas conducted in Sweden and Denmark 1999 to 2002. Overall, 3,740 incident malignant lymphoma patients 18-74 years of age were included. Through a rapid case ascertainment network, the patients were identified shortly after diagnosis. Lymphoma subtypes were reviewed and reclassified according to the WHO classification . Study participants were restricted to individuals with sufficient knowledge of the Danish or Swedish language to answer questions in a telephone interview and without a history of organ transplantation, HIV infection, or other hematopoietic malignancy. In the study overall, 85% of eligible NHL patients were included. In Sweden, a detailed registration of subtypes of all eligible patients permitted evaluation of participation rates by NHL subtype, which was 94% for FL. Early death was an uncommon reason for non-participation among the FL patients (1%). Four-hundred ninety-nine (85%) of all interviewed FL patients also gave blood. The median time from diagnosis to venipuncture was three months (interquartile range: one to five months). In SCALE overall, patients who gave blood were similar to all participating patients with regard to age, sex and educational level.
Clinical data and follow-up
In Sweden, information about established prognostic factors (Ann Arbor stage, performance status, number of involved nodal areas, hemoglobin and lactate dehydrogenase levels), treatments and signs of progression, was collected from medical records. Lymphoma progression was defined as the time to progression, date of start of second-line treatment or lymphoma-specific death, if no second-line treatment was given (data only available in Sweden). Dates and causes of death were obtained through the Cause of Death registry , with complete follow-up through February 2, 2012. In Denmark, clinical parameters (Ann Arbor stage, performance status, treatment) and complete follow-up of survival through December 17, 2009, were obtained from the Danish national lymphoma database (LYFO) . Lymphoma-specific death was defined as having lymphoma as the main underlying cause of death. Based on the number of FLIPI risk factors, the Swedish cases were classified as having low, intermediate or high risk of death (0-1, 2 or ≥3 risk factors, respectively). The Danish patients were classified according to a modified FLIPI based on age (≥60) and stage (I-II vs. III-IV) only, as low risk if they had none of these factors, intermediate if they had one, and high risk if they had two risk factors.
In 400 FL cases for whom a sufficient amount of DNA was available, genotyping of 317,503 SNPs was conducted at the Genome Institute of Singapore within a GWAS using the Illumina HumanHap 300 (version 1.0) array . The resulting dataset was filtered on the basis of genotyping call rates (≥95%), sample completion rate (≥90%), minor allele frequency (≥0.03) and non-deviation from the Hardy-Weinberg equilibrium (p < 10-6). SNPs on the X or Y chromosomes or with cluster plot problems were excluded. Study participants with gender discrepancies and/or labeling errors were removed, as were samples with evidence of cryptic family relationship (identified using the genome command in PLINK ) and population outliers on the basis of their values of the first three principal components (identified using EIGENSTRAT software and unlinked SNPs only ). Complete genotyping was obtained for 298,702 SNPs in 373 (94%) of the genotyped cases.
The University of California San Francisco (UCSF) study
The UCSF study population consisted of 213 FL cases included in a population-based case-control study of NHL (>2,055 cases, 2,081 controls) conducted in the San Francisco Bay area of incident NHL at ages 20-84 years 2001 to 2006. The cases were identified through rapid case ascertainment methods by the Cancer Prevention Institute of California with case reporting supplemented by Surveillance, Epidemiology and End Results data . Overall, 69% of eligible cases participated in the study. Diagnostic material was re-reviewed and classified according to the WHO lymphoma classification by the study pathologist . Blood and/or buccal specimens were collected from 87% of the cases within median 27 days from diagnosis. Genotyping of 329,294 SNPs was performed using the Illumina HumanCNV370-Duo BeadChip (Illumina, Inc., San Diego, CA). SNPs were excluded for low call rates (<90%), low minor allele frequency (<0.03), location on the sex chromosomes or cluster plot problems. Non-European samples according to MDS plot inspection were removed. Vital status and cause of death were updated through August 2012 through record linkage with the Greater Bay Area Cancer registry. Patients were censored at date of death or date of last known contact.
Selection of SNPs associated with FL outcome and established FL risk SNPs
Published studies of germline polymorphisms and FL outcome were found by searching PubMed and Web of Science (through June 30, 2013). We selected SNPs statistically significantly (p < 0.05) associated with overall or lymphoma-specific survival or with time to progression/treatment in at least one study and investigated the association of these with lymphoma-specific death and progression in SCALE and the Swedish cohort, respectively. Fifteen studies and 30 SNPs representing 29 separate loci were identified (Table 1). Genotype data on 22 of these markers (21 loci) were available through genotyping (N = 9) or imputation (N = 13). We used 1000 genomes multi-ethnic reference panel and Impute 2 for imputation (http://mathgen.stats.ox.ac.uk/impute/impute_v2.html#download_reference_data). A strict threshold was set to genotypes with probabilities >0.9, SNPs with information scores >0.8 and call rates >0.9. Six SNPs strongly associated with FL risk in recent GWAS -, of which one (rs6457327 in the HLA class II region) also has been associated with FL prognosis in two studies ,, were also selected for survival analysis.
Genome-wide association studies
In SCALE and UCSF separately, we estimated the association of SNPs with lymphoma-specific death using the Cox proportional hazards model. Time was measured from diagnosis to lymphoma-specific death or end of follow-up. Deaths from other causes than lymphoma were censored. Patients without an event were censored at the time of last known follow-up. The SNPs were coded as 0, 1 or 2 based on the number of minor alleles, and treated as continuous variables in the model. Adjustment was made for age at diagnosis, sex and population stratification (using the first three dimensions calculated by using EIGENSTRAT  in SCALE and multidimensional scaling  in UCSF; in both cohorts unlinked SNPs only were used for these calculations). Data on FLIPI and first-line rituximab was available in SCALE and the Swedish cohort, respectively, but not in UCSF, and hence these covariates were not included in the model. A full meta-analysis of the GWAS results was performed using the DerSimonian-Laird random effects method . PLINK v 1.07 (http://pngu.mgh.harvard.edu/purcell/plink/)  and R  were used for these analyses. In complementary analyses, top variants in the meta-analysis were also adjusted for FLIPI risk groups (in SCALE) and first-line rituximab (in the Swedish cohort). The top locus in the meta-analysis was further explored in a regional association plot, combining the p-values of association in the meta-analysis, LD data from 1000 Genomes European population, gene information from the UCSC browser, and estimated recombination rates (http://csg.sph.umich.edu/locuszoom/) . The potential functionality of the top SNPs and SNPs in strong LD (r2 ≥ 0.8) with these in the 1000 Genomes European population were explored using RegulomeDB, which integrates multiple types of functional data generated by ENCODE and other sources (http://RegulomeDB.org/) .
Investigation of SNPs with reported prognostic impact in FL and established FL risk SNPs
We investigated the association of the selected candidate SNPs and established FL risk SNPs with lymphoma-specific death in SCALE and lymphoma progression in SCALE Sweden using the same Cox model as in the GWAS. The proportional hazard assumption was tested by plotting Schoenfeld residuals against follow-up time; no violations were found. In complementary analyses, associated variants were also adjusted for FLIPI risk groups and first-line rituximab. Also, potential immortal time bias was assessed in the Swedish cohort by starting follow-up at date of venipuncture instead of diagnosis. We further applied the multi-SNP modeling of risk alleles performed in three previous studies (rs5361, rs3799488, rs1799864 and rs1800796 ; rs4073, rs2069762, rs3212227 and rs454078 ; rs1801131, rs1127717 and rs719235 ) on the SCALE dataset. The SNPs were then coded as 0 or 1, based on the presence of the risk (protective/deleterious) allele ,,. The impact of having one, two or three or more as compared to no risk alleles on lymphoma-specific death was assessed in a Cox model, adjusting for age and sex. Multiple testing was adjusted for with the Bonferroni method . The analyses were done in SAS 9.2. Power calculations were performed in R (using survSNP) .
The Danish-Swedish FL cohort was followed for a median of 8.9 years from diagnosis (range 4.3 months to 12.1 years). Median age at diagnosis was 57 years and 65% had Ann Arbor stage III or IV disease (Table 2). One-hundred thirty-seven patients (37%) died of which 88 deaths (64%) were classified as lymphoma-specific. The UCSF cohort was followed for a median of 7.5 years (range 6.9 months to 10.3 years). Twenty-six percent (N = 56) died, of which 43% (N = 24) were classified as lymphoma-specific deaths. In Swedish patients, progression occurred in 67% (N = 155), and rituximab was used in first-line treatment in 10% and overall in 47%. Overall survival differed as expected between risk groups defined by FLIPI (low, intermediate, high; plog-rank < 0.05) in the total SCALE study population and in the Swedish and Danish populations separately (data not shown).
Genome-wide association studies and meta-analysis
In SCALE and UCSF separately, no variants were associated with lymphoma-specific death at a genome-wide significant level (p ≤ 5.0 × 10-8). Thirteen variants in SCALE (on chromosomes 1, 4, 8 and 17; Additional file 1: Figure S1a) and five variants in UCSF (on chromosomes 1, 2, 6, 9 and 20; Additional file 1: Figure S1b) were associated with lymphoma-specific death with a p ≤ 10-6, which was more than would be expected by chance (Additional file 1: Figure S2a and S2b).
In the meta-analysis of the 298,702 variants in SCALE and 319,693 in UCSF, 295,134 variants were present in both cohorts. No variant reached genome-wide significant level of association (p < 5.0 × 10-8). The strongest signals were observed for SNPs on chromosome 1, 4, 17 and 19 (Additional fle 1: Figure S3), for which the p-values were smaller than would be expected by chance (Additional file 1: Figure S4). The strongest association with lymphoma-specific death was observed at rs10491178 (HRrandom = 3.17, 95% CI 2.09-4.79, prandom = 5.24 × 10-8), located on 17q24 in the ATP binding cassette A10 gene (ABCA10) (Table 3). An additional five SNPs on 17q24, close to ABCA6, showed similar associations and were in strong linkage disequilibrium (LD) with rs10491178 (r2 ≤ 0.96 in the SCALE cohort; Additional file 1: Figure S5). One SNP on chromosome 1 showed a suggestive association of prandom = 10-7 (rs3131729 HRrandom = 2.45, 95% CI 1.75-3.44, prandom = 2.22 × 10-7; Table 3). Another SNP on chromosome 1 showed similar association and was linked with rs3131729 (r2 = 0.7) (Additional file 1: Table S1). One intergenic variant on chromosome 4 (rs11932201 HRrandom = 2.16, 95% CI 1.57-2.97, prandom = 2.24 × 10-6) and one variant located on 19q13 (rs2250066 HRrandom = 2.11, 95% CI 1.52-2.94, prandom = 9.38 ×10-6) also showed suggestive associations at prandom = 10-6 (Table 3). For these ten top SNPs, there was no indication of heterogeneity between the studies (pheterogeneity > 0.05, I2 = 0; Additional file 1: Table S1).
Among the top 48 SNPs in the UCSF GWAS with p ≤ 10-5 seven were not genotyped in SCALE. We were able to impute six of these with high confidence (information score >0.8). The meta-analysis of these SNPs gave no additional suggestive associations (data not shown).
There was no correlation between rs10491178 and performance status, lactate dehydrogenase level (Swedish cohort, N = 231) or FLIPI risk groups (SCALE cohort, N = 373; Fisher's Exact test p ≤ 0.15; data not shown). For rs10491178, rs3131729, rs11932201 and rs2250066, additional adjustment for FLIPI and first-line rituximab in SCALE and the Swedish cohort, respectively, did not alter the results meaningfully (data not shown). The association of the top SNP at 17q24 with overall survival (SCALE) was weaker compared with lymphoma-specific death (rs10491178 HR = 2.00, 95% CI 1.35-2.97, p = 5.68 × 10-4), and there was no association with lymphoma progression (Swedish cohort; Additional file 1: Table S2).
Combining data in the RegulomeDB suggested a putative functional role for six of the ten top SNPs, of which rs10491178 (chr 17), rs11932201 (chr 4) and rs2250066 (chr 19) had RegulomeDB scores of 5, indicating that they may affect transcription factor binding or are located in a DNAse hypersensitivity site. There was no data for rs3131729 (chr 1). The lowest RegulomeDB score (i.e. strongest evidence of being in a regulatory site) among our top SNPs and highly linked variants (r2 ≤ 0.8) was found at the SNP rs113464685 in ABCA10 (score 3b), in LD with rs10491178 (data not shown). After imputing the genotypes for rs113464685 in SCALE (as described in Methods for the candidate gene SNPs), we ran the Cox regression for this SNP and found estimates for lymphoma-specific death virtually identical to those of the other SNPs on 17q24 (HR = 3.10, 95% CI 1.98-4.89, p = 1.12 × 10-6).
Investigation of SNPs with reported prognostic impact in FL and established FL risk SNPs
For lymphoma-specific death (SCALE cohort, N = 373), we had ≤80% power to detect HR ≤ 2.0 (or ≤0.5) for rare variants (minor allele frequency (MAF) =0.09 ) and HR ≤ 1.5 (or ≤0.67) for more common variants (MAF = 0.49) (Additionl file 1: Figure S6). For lymphoma progression (Swedish cohort, N = 231), we had ≤80% power to detect an HR ≤ 1.75 (or ≤0.57) for variants with MAF ≤0.09. Among the 22 candidate SNPs evaluated, two imputed variants showed nominally significant association (p < 0.05) with lymphoma-specific death: rs1801131 in MTHFR (HR = 0.69, 95% CI 0.49-0.97, p = 0.03) and rs2069762 in IL2 (HR = 0.63, 95% CI 0.43-0.92, p = 0.02; Table 4). rs1801131 was also associated with lymphoma progression (HR = 0.59, 95% CI 0.45-0.77, p = 0.0001). The direction of association for rs1801131 and rs2069762 with the two lymphoma-specific outcomes was opposite of what was previously reported for overall survival (Table 1) ,. In line with a previous study investigating event-free survival in FL , we found that the C allele of rs2466571 in the CD46 gene was associated with shorter time to lymphoma progression (HR = 1.37, 95% CI 1.10-1.72, p = 0.006; Table 4). Likewise, the estimates for two imputed tightly linked SNPs (r2 = 0.99) in the IL8 gene, rs4073 and rs2227307 (HRrs4073 = 0.78, 95% CI 0.62-0.97, p = 0.02; HRrs2227307 = 0.75, 95% CI 0.60-0.94, p = 0.01), and lymphoma progression were in the same direction as in two previous studies investigating overall survival in FL (Table 4) ,. Additional adjustment with FLIPI categories or first-line rituximab only marginally altered these results (data not shown). Starting follow-up at date of venipuncture instead of diagnosis also did not change the results (data not shown). After correcting for performing 54 tests (27 SNPs and two outcomes), the association of rs1801131 (MTHFR) with lymphoma progression remained significant (Bonferroni corrected p = 0.01).
We finally stratified patients by the number of deleterious/protective alleles, as defined in three previous studies ,,, but found no significant associations with the outcomes in the present study (data not shown). None of the six SNPs previously reported to be associated with risk of FL - were significantly associated with lymphoma-specific death or progression (Table 4 (rs6457327) and Additional file 1: Table S3).
In this GWAS meta-analysis of common genetic variation in FL prognosis in a total of 586 patients, we did not identify any variants associated with lymphoma-specific death at a genome-wide significant level (p < 5.0 ×10-8). The strongest association was observed for tightly linked SNPs located on 17q24 near the ABCA10 and ABCA6 genes (rs10491178 HRrandom = 3.17, 95% CI 2.09-4.79, prandom = 5.24 ×10-8). Investigating previously reported candidate genes, we found further support of a role for rs2466571 in CD46 in FL progression. Two SNPs in IL8 (rs4073 and rs2227307) previously linked with overall survival in FL, where here associated with FL progression. We also observed an association with MTHFR and FL progression but the direction was opposite of a previous report on overall survival in FL.
Genome-wide association study and meta-analysis
The increased rate of lymphoma-specific death for carriers of the A allele of rs10491178 in ABCA10 on 17p24 was consistent in the SCALE and UCSF cohorts, strengthening the notion of a possible causal association with this region. The association was independent of established prognostic risk factors (lactate dehydrogenase levels, performance status and FLIPI). ABCA10 is clustered among four other members of the ATP binding cassette (ABC) 1 family on 17q24, including ABCA6 harboring another five linked associated SNP markers. The ABC transporters - a family of proteins responsible for the movement of a wide variety of xenobiotics, lipids and metabolic products across the cell membranes - are implicated in multidrug resistance . They are overexpressed in several tumor types , including FL . The exact function of ABCA10 and ABCA6 remains to be elucidated , but the ABCA3 transporter, similar in terms of amino acid sequences and structural organization to ABCA10 and ABCA6 , has been shown to impede the efficacy of rituximab in aggressive B-cell lymphomas, and vincristine, anthracyclines and etopside in acute myeloid leukemia when expressed in high levels . Interestingly, rs10491178 (C > T) encodes a premature stop codon (Arg (CGA) -- > Stop (TGA)) in the ABCA10 transcript (Refseq NM_080282; Ensembl ABCA10-001) .
We also looked up SNPs in high LD (r2 > 0.8) with rs10491178 for potential regulatory functions in RegulomeDB, and identified another SNP rs113464685 in ABCA10, in strong LD (r2 = 1.00) with the top SNPs at 17q24, that lies in the binding site of the transcription factor paired box 5 (PAX5). PAX5 has a pivotal regulatory function in B-cell development and its aberrant expression is correlated with aggressive subsets of B-cell NHL . This SNP is thus an alternative candidate to explain the observed association with the region.
The suggestive associations (p ≤ 10-6) of rs3131729 in the Dab, reelin signal transducer, homolog 1 (DAB1) gene on 1p32-p31, the intergenic SNP rs11932201 on chromosome 4, and rs2250066 in the kallikrein (KLK) 11 gene on 19q13, were also consistent in the SCALE and UCSF studies. The role of the DAB1 protein in the developing nervous system is well studied but for other tissues knowledge is scarce . rs11932201 is flanked by two genes (long intergenic non-protein coding RNA 1085 (LINC01085) and CPEB2 antisense RNA 1 (CPEB2-AS1)) of which little is known. The KLK locus on 19q13 harbors genes encoding a family of serine proteases associated with cancer risk and prognosis in several studies .
Investigation of SNPs with reported prognostic impact in FL and established FL risk SNPs
Upon investigating SNPs linked to FL survival in previous studies, we observed a shorter time to lymphoma progression for C compared to A allele carriers at rs2466571 in CD46 on 1q32 (HR = 1.37, ptrend = 0.006), at a nominally significant level. This is consistent with a recent study of 107 FL patients investigating event-free survival in FL (HR = 1.49, ptrend < 0.05) . CD46 is a complement inhibitor, protecting the cell against complement-mediated lysis . High expression levels of membrane-bound complement regulatory proteins such as CD46 in tumors has been shown to suppress anti-tumor T-cell responses, and may inhibit anti-tumor therapeutic activity of monoclonal antibodies, including rituximab in B-cell NHL .
The minor alleles of two SNPs in IL8, rs4073 and rs2227307, on 4q12-q21 were associated with longer time to progression (rs4073 HRA = 0.78, p = 0.02; rs2227307 HRG = 0.75, p = 0.01). Two previous studies reported an association for the same two SNPs with overall survival in FL; Cerhan et al. observed a statistically significant longer overall survival with homozygosity of the same alleles (rs4073 HRAA = 0.47, 95% CI 0.27-0.79, rs2227307 HRGG = 0.53, 95% CI 0.31-0.91) , whereas Aschebrook-Kilfoy et al. observed that the opposite alleles were associated with shorter overall survival (rs4073 HRTT = 2.60, 95% CI 1.10-6.15, rs2227307 HRCC = 2.57, 95% CI 1.07-6.17)  (Table 1). Time to FL progression was not evaluated in these former investigations. The A allele of rs4073 has been associated with increased production of IL8 .
In our study, rs1801131 in MTHFR was associated with time to progression and lymphoma-specific death, and rs2069762 in the promoter of IL2 with lymphoma-specific death. However, these results are conflicting with previous reports of overall survival in FL and are therefore difficult to interpret ,. For rs1801274, located in FCGR2A on 1q23, one of the two most studied variants -,,,,,, we observed no associations with lymphoma progression or lymphoma-specific death. This is in line with the results of most previous studies, irrespective of how outcome was defined (progression-free ,,,,, event-free , or overall survival ,).
Several SNPs, including rs6457327, rs10484561 and rs2647012 in the HLA class I and II regions on 6p21, are convincingly associated with FL risk ,,, suggesting an important role of this genetic region in the development of FL, and thereby possibly also in FL progression. Indeed, rs6457327 has been associated with FL transformation , and overall survival  in two previous studies (Table 1). However, we did not observe any clear associations for rs6457327, rs10484561 and rs2647012 with FL progression or lymphoma-specific death in our cohort.
Strengths and weaknesses
The strengths of the present analysis include the population-based design of the two cohorts, the large number of patients compared to previous studies, and the comprehensive evaluation of the impact of genetic variation in FL outcome. Population stratification was adjusted for in the GWAS and we observed little or no evidence of inflation of the statistics (lambda GCSCALE = lambda GCUCSF = 1.00) . In the analysis of candidate gene variants and lymphoma progression or death, adjustment for established risk factors including age, sex, FLIPI categories and first-line rituximab resulted in marginal changes of the point estimates only.
Although candidate gene studies have obvious weaknesses compared with genome-wide studies, a strength include a higher prior probability of true findings, since they are based on a biological understanding of cancer survival pathways . Still, many candidate gene findings cannot be replicated, suggesting some false-positive results in the existing literature . Most of the previous studies were small (N < 150) and the number of SNPs tested were sometimes large (up to 6679 variants), resulting in a considerable risk of chance findings. However, the present study was limited to detect moderate to strong effects (HR >1.5-2.0) for the investigated SNPs with acceptable probability (≤0.80). Hence, chance findings cannot be excluded.
In the study of genetic variation and FL prognosis, different outcome measures have been used, including overall and lymphoma-specific survival and event-free survival. Progression-free survival, defined as the time from study entry to lymphoma progression or all-cause death, is generally recognized as the most valid measure in intervention studies . We argue that lymphoma-specific survival/death may be superior to overall survival (but not to progression-free survival) if the study aims to test an association with lymphoma progression leading to death. In view of the high median age at diagnosis of FL and the often indolent clinical course , a relatively large proportion of the patients are expected to die of non-lymphoma related causes, which could dilute associations with lymphoma progression and death if overall survival is used instead. On the other hand, with the use of lymphoma-specific survival some deaths due in part to progression may have been missed potentially leading to a lower specificity. We did not observe similar results for three out of four top genetic variants and lymphoma progression, as with lymphoma-specific death, although this analysis could only be performed in a subset of the patients. Explanations for such a difference could lie in the definitions of these outcomes. While time to progression and progression-free survival primarily reflect the first part of the follow-up period, including mainly time to first-line or second-line treatment, time to lymphoma-specific death reflects a larger part of the follow-up and thus, to a higher degree, treatment response or lack of such. In the future, attempts to validate the current findings, and when comparing results between investigations, the use of similar definitions of outcome will be important for the interpretation.
In the present GWAS meta-analysis, there was suggestive evidence that inherited polymorphisms in the ABCA10 or ABCA6 genes may be associated with risk of lymphoma-specific death. In candidate-gene analysis, an association with lymphoma progression was observed for SNPs in CD46 and IL8, previously linked with lymphoma progression (CD46) and overall survival (IL8). These findings need further confirmation in other independent populations or in a larger multicenter GWAS.
KES was the principal investigator and takes primary responsibility for the paper; HOA, BG and MM recruited the patients; HH and MM provided data for the Danish cohort; JNF and JL performed and directed the laboratory work for SCALE; FB collected clinical data from medical records in Sweden; FB, JNF, HD, RK and KH participated in the statistical analysis; FB and KES coordinated the research; PB, LC and CS provided data for the UCSF validation cohort and performed the survival analysis in this cohort; FB and KES drafted the paper; JNF, HH, HOA, BG, MM, HD, RR, KH, JL, PB, LC and CS revised the manuscript and approved of its final version before submission. All authors read and approved the final manuscript.
Freedman A: Follicular lymphoma: 2014 update on diagnosis and management. Am J Hematol. 2014, 89 (4): 429-436. 10.1002/ajh.23674.
Relander T, Johnson NA, Farinha P, Connors JM, Sehn LH, Gascoyne RD: Prognostic factors in follicular lymphoma. J Clin Oncol. 2010, 28 (17): 2902-2913. 10.1200/JCO.2009.26.1693.
Aschebrook-Kilfoy B, Zheng T, Foss F, Ma S, Han X, Lan Q, Holford T, Chen Y, Leaderer B, Rothman N, Zhang Y: Polymorphisms in immune function genes and non-Hodgkin lymphoma survival. J Cancer Surviv. 2012, 6 (1): 102-114. 10.1007/s11764-010-0164-4.
Berglund M, Enblad G, Thunberg U: SNP rs6457327 is a predictor for overall survival in follicular lymphoma as well as survival after transformation. Blood. 2011, 118 (16): 4489-10.1182/blood-2011-06-362921.
Carlotti E, Palumbo GA, Oldani E, Tibullo D, Salmoiraghi S, Rossi A, Golay J, Pulsoni A, Foa R, Rambaldi A: FcgammaRIIIA and FcgammaRIIA polymorphisms do not predict clinical outcome of follicular non-Hodgkin's lymphoma patients treated with sequential CHOP and rituximab. Haematologica. 2007, 92 (8): 1127-1130. 10.3324/haematol.11288.
Cartron G, Dacheux L, Salles G, Solal-Celigny P, Bardos P, Colombat P, Watier H: Therapeutic activity of humanized anti-CD20 monoclonal antibody and polymorphism in IgG Fc receptor FcgammaRIIIa gene. Blood. 2002, 99 (3): 754-758. 10.1182/blood.V99.3.754.
Cerhan JR, Wang S, Maurer MJ, Ansell SM, Geyer SM, Cozen W, Morton LM, Davis S, Severson RK, Rothman N, Lynch CF, Wacholder S, Chanock SJ, Habermann TM, Hartge P: Prognostic significance of host immune gene polymorphisms in follicular lymphoma survival. Blood. 2007, 109 (12): 5439-5446. 10.1182/blood-2006-11-058040.
Charbonneau B, Maurer MJ, Fredericksen ZS, Zent CS, Link BK, Novak AJ, Ansell SM, Weiner GJ, Wang AH, Witzig TE, Dogan A, Slager SL, Habermann TM, Cerhan JR: Germline variation in complement genes and event-free survival in follicular and diffuse large B-cell lymphoma. Am J Hematol. 2012, 87 (9): 880-885. 10.1002/ajh.23273.
Charbonneau B, Wang AH, Maurer MJ, Asmann YW, Zent CS, Link BK, Ansell SM, Weiner GJ, Ozsan N, Feldman AL, Witzig TE, Cunningham JM, Dogan A, Habermann TM, Slager SL, Novak AJ, Cerhan JR: CXCR5 polymorphisms in non-Hodgkin lymphoma risk and prognosis. Cancer Immunol Immunother. 2013, 62 (9): 1475-1484. 10.1007/s00262-013-1452-4.
Ghesquieres H, Cartron G, Seymour JF, Delfau-Larue MH, Offner F, Soubeyran P, Perrot A, Brice P, Bouabdallah R, Sonet A, Dupuis J, Casasnovas O, Catalano JV, Delmer A, Jardin F, Verney A, Dartigues P, Salles G: Clinical outcome of patients with follicular lymphoma receiving chemoimmunotherapy in the PRIMA study is not affected by FCGR3A and FCGR2A polymorphisms. Blood. 2012, 120 (13): 2650-2657. 10.1182/blood-2012-05-431825.
Ghielmini M, Rufibach K, Salles G, Leoncini-Franscini L, Leger-Falandry C, Cogliatti S, Fey M, Martinelli G, Stahel R, Lohri A, Ketterer N, Wernli M, Cerny T, Schmitz SF: Single agent rituximab in patients with follicular or mantle cell lymphoma: clinical and biological factors that are predictive of response and event-free survival as well as the effect of rituximab on the immune system: a study of the Swiss Group for Clinical Cancer Research (SAKK). Ann Oncol. 2005, 16 (10): 1675-1682. 10.1093/annonc/mdi320.
Gibson TM, Wang SS, Cerhan JR, Maurer MJ, Hartge P, Habermann TM, Davis S, Cozen W, Lynch CF, Severson RK, Rothman N, Chanock SJ, Morton LM: Inherited genetic variation and overall survival following follicular lymphoma. Am J Hematol. 2012, 87 (7): 724-726. 10.1002/ajh.23184.
Persky DO, Dornan D, Goldman BH, Braziel RM, Fisher RI, Leblanc M, Maloney DG, Press OW, Miller TP, Rimsza LM: Fc gamma receptor 3a genotype predicts overall survival in follicular lymphoma patients treated on SWOG trials with combined monoclonal antibody plus chemotherapy but not chemotherapy alone. Haematologica. 2012, 97 (6): 937-942. 10.3324/haematol.2011.050419.
Racila E, Link BK, Weng WK, Witzig TE, Ansell S, Maurer MJ, Huang J, Dahle C, Halwani A, Levy R, Weiner GJ: A polymorphism in the complement component C1qA correlates with prolonged response following rituximab therapy of follicular lymphoma. Clin Cancer Res. 2008, 14 (20): 6697-6703. 10.1158/1078-0432.CCR-08-0745.
Wang SS, Maurer MJ, Morton LM, Habermann TM, Davis S, Cozen W, Lynch CF, Severson RK, Rothman N, Chanock SJ, Hartge P, Cerhan JR: Polymorphisms in DNA repair and one-carbon metabolism genes and overall survival in diffuse large B-cell lymphoma and follicular lymphoma. Leukemia. 2009, 23 (3): 596-602. 10.1038/leu.2008.240.
Weng WK, Czerwinski D, Timmerman J, Hsu FJ, Levy R: Clinical outcome of lymphoma patients after idiotype vaccination is correlated with humoral immune response and immunoglobulin G Fc receptor genotype. J Clin Oncol. 2004, 22 (23): 4717-4724. 10.1200/JCO.2004.06.003.
Weng WK, Levy R: Two immunoglobulin G fragment C receptor polymorphisms independently predict response to rituximab in patients with follicular lymphoma. J Clin Oncol. 2003, 21 (21): 3940-3947. 10.1200/JCO.2003.05.013.
Wrench D, Leighton P, Skibola CF, Conde L, Cazier JB, Matthews J, Iqbal S, Carlotti E, Bodor C, Montoto S, Calaminici M, Gribben JG, Lister TA, Fitzgibbon J: SNP rs6457327 in the HLA region on chromosome 6p is predictive of the transformation of follicular lymphoma. Blood. 2011, 117 (11): 3147-3150. 10.1182/blood-2010-10-315382.
Gribben JG: Implications of the tumor microenvironment on survival and disease response in follicular lymphoma. Curr Opin Oncol. 2010, 22 (5): 424-430. 10.1097/CCO.0b013e32833d5938.
Phan VH, Tan C, Rittau A, Xu H, McLachlan AJ, Clarke SJ: An update on ethnic differences in drug metabolism and toxicity from anti-cancer drugs. Expert Opin Drug Metab Toxicol. 2011, 7 (11): 1395-1410. 10.1517/17425255.2011.624513.
Bruhns P, Iannascoli B, England P, Mancardi DA, Fernandez N, Jorieux S, Daeron M: Specificity and affinity of human Fcgamma receptors and their polymorphic variants for human IgG subclasses. Blood. 2009, 113 (16): 3716-3725. 10.1182/blood-2008-09-179754.
Conde L, Halperin E, Akers NK, Brown KM, Smedby KE, Rothman N, Nieters A, Slager SL, Brooks-Wilson A, Agana L, Riby J, Liu J, Adami HO, Darabi H, Hjalgrim H, Low HQ, Humphreys K, Melbye M, Chang ET, Glimelius B, Cozen W, Davis S, Hartge P, Morton LM, Schenk M, Wang SS, Armstrong B, Kricker A, Milliken S, Purdue MP, et al: Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32. Nat Genet. 2010, 42 (8): 661-664. 10.1038/ng.626.
Morton LM, Purdue MP, Zheng T, Wang SS, Armstrong B, Zhang Y, Menashe I, Chatterjee N, Davis S, Lan Q, Vajdic CM, Severson RK, Holford TR, Kricker A, Cerhan JR, Leaderer B, Grulich A, Yeager M, Cozen W, Hoar Zahm S, Chanock SJ, Rothman N, Hartge P: Risk of non-Hodgkin lymphoma associated with germline variation in genes that regulate the cell cycle, apoptosis, and lymphocyte development. Cancer Epidemiol Biomarkers Prev. 2009, 18 (4): 1259-1270. 10.1158/1055-9965.EPI-08-1037.
Skibola CF, Bracci PM, Halperin E, Conde L, Craig DW, Agana L, Iyadurai K, Becker N, Brooks-Wilson A, Curry JD, Spinelli JJ, Holly EA, Riby J, Zhang L, Nieters A, Smith MT, Brown B: Genetic variants at 6p21.33 are associated with susceptibility to follicular lymphoma. Nat Genet. 2009, 41 (8): 873-875. 10.1038/ng.419.
Smedby KE, Foo JN, Skibola CF, Darabi H, Conde L, Hjalgrim H, Kumar V, Chang ET, Rothman N, Cerhan JR, Brooks-Wilson AR, Rehnberg E, Irwan ID, Ryder LP, Brown PN, Bracci PM, Agana L, Riby J, Cozen W, Davis S, Hartge P, Morton LM, Severson RK, Wang SS, Slager SL, Fredericksen ZS, Novak AJ, Kay NE, Habermann TM, Armstrong B, et al: GWAS of follicular lymphoma reveals allelic heterogeneity at 6p21.32 and suggests shared genetic susceptibility with diffuse large B-cell lymphoma. PLoS Genet. 2011, 7 (4): e1001378-10.1371/journal.pgen.1001378.
Wang SS, Purdue MP, Cerhan JR, Zheng T, Menashe I, Armstrong BK, Lan Q, Hartge P, Kricker A, Zhang Y, Morton LM, Vajdic CM, Holford TR, Severson RK, Grulich A, Leaderer BP, Davis S, Cozen W, Yeager M, Chanock SJ, Chatterjee N, Rothman N: Common gene variants in the tumor necrosis factor (TNF) and TNF receptor superfamilies and NF-kB transcription factors and non-Hodgkin lymphoma risk. PLoS One. 2009, 4 (4): e5360-10.1371/journal.pone.0005360.
Smedby KE, Hjalgrim H, Melbye M, Torrang A, Rostgaard K, Munksgaard L, Adami J, Hansen M, Porwit-MacDonald A, Jensen BA, Roos G, Pedersen BB, Sundstrom C, Glimelius B, Adami HO: Ultraviolet radiation exposure and risk of malignant lymphomas. J Natl Cancer Inst. 2005, 97 (3): 199-209. 10.1093/jnci/dji022.
Jaffe E: Pathology and genetics of tumours of hematopoietic and lymphoid tissues. 2001, IARC Press, International Agency for Research on Cancer, Lyon, France
The Swedish Cause of Death Register, The National Board of Health and Welfare, Sweden.Available at:. Accessed May 2, 2012. In., [http://www.socialstyrelsen.se/register/dodsorsaksregistret]
D'Amore F, Christensen BE, Brincker H, Pedersen NT, Thorling K, Hastrup J, Pedersen M, Jensen MK, Johansen P, Andersen E, Bach B, Sørensen E: Clinicopathological features and prognostic factors in extranodal non-Hodgkin lymphomas. Danish LYFO Study Group. Eur J Cancer. 1991, 27 (10): 1201-1208. 10.1016/0277-5379(91)90081-N.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38 (8): 904-909. 10.1038/ng1847.
Skibola CF, Bracci PM, Halperin E, Nieters A, Hubbard A, Paynter RA, Skibola DR, Agana L, Becker N, Tressler P, Forrest MS, Sankararaman S, Conde L, Holly EA, Smith MT: Polymorphisms in the estrogen receptor 1 and vitamin C and matrix metalloproteinase gene families are associated with susceptibility to lymphoma. PLoS One. 2008, 3 (7): e2816-10.1371/journal.pone.0002816.
DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7 (3): 177-188. 10.1016/0197-2456(86)90046-2.
R: A Language and Environment for Statistical Computing. 2010, R Foundation for Statistical Computing, Vienna, Austria
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ: LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010, 26 (18): 2336-2337. 10.1093/bioinformatics/btq419.
Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, Kasowski M, Karczewski KJ, Park J, Hitz BC, Weng S, Cherry JM, Snyder M: Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 2012, 22 (9): 1790-1797. 10.1101/gr.137323.112.
Bland JM, Altman DG: Multiple significance tests: the Bonferroni method. BMJ. 1995, 310 (6973): 170-10.1136/bmj.310.6973.170.
Falasca M, Linton KJ: Investigational ABC transporter inhibitors. Expert Opin Investig Drugs. 2012, 21 (5): 657-666. 10.1517/13543784.2012.679339.
Yagi K, Yamamoto K, Umeda S, Abe S, Suzuki S, Onishi I, Kirimura S, Fukayama M, Arai A, Kitagawa M, Kurata M: Expression of multidrug resistance 1 gene in B-cell lymphomas: association with follicular dendritic cells. Histopathology. 2013, 62 (3): 414-420. 10.1111/his.12035.
Coleman JA, Quazi F, Molday RS: Mammalian P4-ATPases and ABC transporters and their role in phospholipid transport. Biochim Biophys Acta. 2013, 1831 (3): 555-574. 10.1016/j.bbalip.2012.10.006.
Aung T, Chapuy B, Vogel D, Wenzel D, Oppermann M, Lahmann M, Weinhage T, Menck K, Hupfeld T, Koch R, Trumper L, Wulf GG: Exosomal evasion of humoral immunotherapy in aggressive B-cell lymphoma modulated by ATP-binding cassette transporter A3. Proc Natl Acad Sci U S A. 2011, 108 (37): 15336-15341. 10.1073/pnas.1102855108.
Flicek P, Amode MR, Barrell D, Beal K, Billis K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fitzgerald S, Gil L, Giron CG, Gordon L, Hourlier T, Hunt S, Johnson N, Juettemann T, Kahari AK, Keenan S, Kulesha E, Martin FJ, Maurel T, McLaren WM, Murphy DN, Nag R, Overduin B, Pignatelli M, Pritchard B, Pritchard E, Riat HS, et al: Ensembl 2014. Nucleic Acids Res. 2014, 42 (Database issue): D749-755. 10.1093/nar/gkt1196.
O'Brien P, Morin P, Ouellette RJ, Robichaud GA: The Pax-5 gene: a pluripotent regulator of B-cell differentiation and cancer disease. Cancer Res. 2011, 71 (24): 7345-7350. 10.1158/0008-5472.CAN-11-1874.
Vazquez-Carretero MD, Garcia-Miranda P, Calonge ML, Peral MJ, Ilundain AA: Dab1 and reelin participate in a common signal pathway that controls intestinal crypt/villus unit dynamics. Biol Cell. 2014, 106 (3): 83-96. 10.1111/boc.201300078.
Batra J, O'Mara T, Patnala R, Lose F, Clements JA: Genetic polymorphisms in the human tissue kallikrein (KLK) locus and their implication in various malignant and non-malignant diseases. Biol Chem. 2012, 393 (12): 1365-1390. 10.1515/hsz-2012-0211.
Yan J, Allendorf DJ, Li B, Yan R, Hansen R, Donev R: The role of membrane complement regulatory proteins in cancer immunotherapy. Adv Exp Med Biol. 2008, 632: 159-174.
Wang N, Zhou R, Wang C, Guo X, Chen Z, Yang S, Li Y: -251 T/A polymorphism of the interleukin-8 gene and cancer risk: a HuGE review and meta-analysis based on 42 case-control studies. Mol Biol Rep. 2012, 39 (3): 2831-2841. 10.1007/s11033-011-1042-5.
Wang D, Sun Y, Stang P, Berlin JA, Wilcox MA, Li Q: Comparison of methods for correcting population stratification in a genome-wide association study of rheumatoid arthritis: principal-component analysis versus multidimensional scaling. BMC Proc. 2009, 3 (Suppl 7): S109-10.1186/1753-6561-3-s7-s109.
Pasche B, Yi N: Candidate gene association studies: successes and failures. Curr Opin Genet Dev. 2010, 20 (3): 257-261. 10.1016/j.gde.2010.03.006.
Cheson BD, Pfistner B, Juweid ME, Gascoyne RD, Specht L, Horning SJ, Coiffier B, Fisher RI, Hagenbeek A, Zucca E, Rosen ST, Stroobants S, Lister TA, Hoppe RT, Dreyling M, Tobinai K, Vose JM, Connors JM, Federico M, Diehl V: Revised response criteria for malignant lymphoma. J Clin Oncol. 2007, 25 (5): 579-586. 10.1200/JCO.2006.09.2403.
The study was financed by: Swedish Cancer Society 2009/1084, 2012/774 (KES); National Institutes of Health 5R01 CA69669-02 (MM, HOA) and R01CA87014, R03CA143947, R03CA150037, R01CA104682, RO1CA122663, RO1CA154643 (CFS, PMB); Lundbeck Foundation Grant R19-A2364 (HH, MM); Danish Cancer Society Grant DP 08-155 (HH, MM). FB was supported by a doctoral student grant and KES by the Strategic Research Program in Epidemiology at Karolinska Institutet. The financiers had no role in the study design, collection, analysis and interpretation of data, in the writing of the manuscript, or in the decision to submit the manuscript for publication.
The authors declare that they have no competing interest.
Electronic supplementary material
Additional file 1: Table S1.: Top SNPs associated with lymphoma-specific death (prandom<5.0 x10-5) in the meta-analysis of FL patients in SCALE and UCSF (N=586). The minor allele (A1) was investigated for association. A2=major allele, MAF= minor allele frequency, NA= not applicable. Relative risk of lymphoma-specific and all-cause death in SCALE (N=373), and lymphoma progression in SCALE Sweden (N=231) for the top SNPs in the pooled analysis of SCALE and UCSF. The minor allele was investigated for association. Table S3. Relative risk of lymphoma-specific death and lymphoma progression for SNPs previously associated with follicular lymphoma risk. Figure S1. Manhattan plot of the -log10(p-values) for the association with lymphoma-specific death in SCALE (A) and UCSF (B), by chromosome and within chromosome location. The top ten SNPs in the meta-analysis are highlighted (light green). The red line marks -log10(5.0 x10-8). Figure S2. QQ plot of the observed versus the expected -log(p-values) for the association with lymphoma-specific death in SCALE (A) and UCSF (B). Figure S3. Manhattan plot of the -log10(p-values) for the association with lymphoma-specific death in the meta-analysis of SCALE and UCSF. The top ten SNPs are highlighted (light green). The red line marks -log10(5.0 x10-8). Figure S4. QQ plot of the observed versus the expected -log(p-values) for the association with lymphoma-specific death in the meta-analysis of SCALE and UCSF. Figure S5. Regional association plot of the top locus in the combined analysis on 17q, displaying the p-values of association in SCALE, extent of LD with rs10491178 as reference, relation of SNPs and genes, and estimated recombination rates at each position. Figure S6. Estimations of the power to detect associations between selected candidate SNPs and lymphoma-specific death and lymphoma progression, respectively, for different effect sizes and minor allele frequencies (MAF). (PDF 729 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Baecklund, F., Foo, JN., Bracci, P. et al. A comprehensive evaluation of the role of genetic variation in follicular lymphoma survival. BMC Med Genet 15, 113 (2014). https://doi.org/10.1186/s12881-014-0113-6
- Follicular lymphoma
- Single nucleotide polymorphism
- Genome-wide association study
- Candidate gene study