Organochlorine exposure was linked to non-Hodgkin lymphoma (NHL) risk. To determine whether this relation is modified by immune gene variation, we genotyped 61 polymorphisms in 36 immune genes in 1172 NHL cases and 982 controls from the National Cancer Institute–Surveillance, Epidemiology, and End Results (NCI-SEER) study. We examined 3 exposures with elevated risk in this study: PCB180 (plasma, dust measurements), the toxic equivalency quotient (an integrated functional measure of several organochlorines) in plasma, and α-chlordane (dust measurements, self-reported termiticide use). Plasma (100 cases, 100 controls) and dust (682 cases, 513 controls) levels were treated as natural log-transformed continuous variables. Unconditional logistic regression was used to calculate β coefficients and odds ratios, stratified by genotype. Associations between all 3 exposures and NHL risk were limited to the same genotypes for IFNG (C−1615T) TT and IL4 (5′-UTR, Ex1-168C>T) CC. Associations between PCB180 in plasma and dust and NHL risk were limited to the same genotypes for IL16 (3′-UTR, Ex22+871A>G) AA, IL8 (T−251A) TT, and IL10 (A−1082G) AG/GG. This shows that the relation between organochlorine exposure and NHL risk may be modified by particular variants in immune genes and provides one of the first examples of a potential gene-environment interaction for NHL.
Introduction
Organochlorine compounds such as polychlorinated biphenyls (PCBs), dioxins, furans, and certain pesticides are ubiquitous environmental pollutants that accumulate in adipose tissue. Several epidemiologic studies have suggested that organochlorine compounds may play a role in the development of non-Hodgkin lymphoma (NHL).1,,,,,,,,,–11 The National Cancer Institute–Surveillance, Epidemiology, and End Results (NCI-SEER) population-based case-control study of NHL, the parent study for the current analysis, examined these compounds with the use of data from questionnaires about termiticide use and analysis of plasma and carpet dust samples for organochlorine compounds. PCB180 levels in plasma4 and carpet dust2 were positively associated with NHL risk, as were plasma levels of PCB156 and PCB194, the coplanar PCB169, total furans and several individual furan congeners, and the toxic equivalency quotient (TEQ, a toxicity-weighted summary metric indicating biologic activity through the arylhydrocarbon receptor pathway).4 NCI-SEER study participants were also at increased risk of developing NHL if levels of the termiticide chlordane were elevated in their carpet dust, or if their homes had been treated for termites before the 1988 ban on chlordane in the United States.3 Although the association between particular organochlorines, or organochlorines as a class, and risk of NHL has not been established, the epidemiologic data are highly suggestive. There is a substantial amount of data, predominantly experimental, showing that these compounds can influence immune function.12,,,,–17
Common variants in immune and inflammatory response genes have been associated with NHL risk in several studies.18,,,,–23 Consistent main effect associations with NHL have been reported for one or more polymorphisms in tumor necrosis factor (TNF)20,21,23 and interleukin-10 (IL10),19,–21,23 with growing evidence for interleukin-4 (IL4).19,20
The identification of variants in one or more immune genes influencing NHL risk, and the epidemiologic and experimental evidence suggesting an association between organochlorine exposures and NHL, led us to evaluate the relation between organochlorine exposures and NHL risk according to the presence or absence of polymorphisms in a wide spectrum of immune and inflammatory response genes. We previously reported the effects of organochlorine exposure limited to genotypes of IL10 and TNF24 ; here, we evaluate more broadly the relation between these exposures and genetic variants in 36 genes.
Methods
Study population
The study population was described previously.25 Cases included 1321 patients newly diagnosed with NHL at ages 20 to 74 between July 1, 1998, and June 30, 2000, at 4 SEER registries (the state of Iowa, Los Angeles County, and the metropolitan areas of Detroit and Seattle). Population controls (n = 1057), frequency-matched to cases on age, sex, race, and SEER registry were identified by random digit dialing (younger than age 65) or Medicare files (age 65 and older). Overall participation rates were 76% in cases and 52% in controls; overall response rates were 59% and 44%, respectively. The study was approved by human subjects review boards at all participating institutions, and written informed consent was obtained from each participant before interview in accordance with the Declaration of Helsinki. All cases were histologically confirmed and coded according to the International Classification of Diseases for Oncology, 2nd edition.26
Exposure assessment
Participants were mailed a residential and occupational history calendar. During the subsequent home visit, the interviewer administered a computer-assisted personal interview that included a detailed history of pest treatment in each residence occupied for 2 years or more since 1970. Among the questions asked was whether any of the homes had been treated for termites and, if so, whether the treatment occurred before or after 1988. Used vacuum cleaner bags were collected from eligible participants and sent to Southwest Research Institute (San Antonio, TX), where dust samples (682 cases and 513 controls meeting eligibility criteria) were analyzed for insecticides and PCBs.2,3,27 A multiple imputation procedure was used to assign values to missing data2,28 ; in the present analysis, we use values from one of those imputations. We include in the current analysis the 2 analytes whose levels in carpet dust were significantly associated with NHL risk: PCB1802 and α-chlordane.3
Interviewed participants were asked to provide a venous blood or mouthwash buccal cell sample. From a subset of 100 untreated cases and 100 controls with blood samples, we measured persistent organochlorine chemical levels in plasma4 (approximately one-half of these participants also had dust measurements for organochlorines). A multiple imputation procedure was used to assign values to missing data determined to be below the detection limits; as with dust, we used one of those imputations for the current analysis. We include PCB180 in the current analysis because its levels in both plasma and carpet dust were significantly associated with NHL risk. We also include the TEQ, which was calculated as the sum of the lipid-adjusted levels of dioxin, furan, coplanar PCB, and mono-ortho-chlorinated PCB congeners, each weighted by its congener-specific toxicity equivalency factor (TEF, congener potency relative to TCDD). The TEQ calculation is based on 11 congeners for which World Health Organization TEFs have been assigned (PCB180 is not among them) and that were detected in at least 30% of the plasma samples.4,29,30 We chose to include the TEQ rather than its individual components because the latter are highly correlated.
We genotyped blood samples from 773 cases and 668 controls and buccal cell samples from 399 cases and 314 controls, as previously described.23 We identified for analysis 61 single nucleotide polymorphisms (SNPs) in 36 proinflammatory and other immunoregulatory genes (Table 1). As previously reported,23 these variants were selected according to a priori laboratory evidence that suggested functional consequences for an allele or reported associations with NHL or related risk factors (eg, autoimmune disorders, infectious diseases). Genotyping was conducted at the NCI Core Genotyping Facility (Gaithersburg, MD), as previously described.23 Agreement for 40 blinded quality control replicates and 100 blinded duplicates was at least 99% for all assays. Successful genotyping was achieved for 96% to 100% of DNA samples for all SNPs.
Statistical analysis
Measurements in carpet dust (PCB180 and α-chlordane) and plasma (PCB180 and TEQ) were log-normally distributed and were treated as natural log-transformed continuous variables in the analysis. Unconditional logistic regression models were used to calculate the β coefficient and P value for each risk factor, stratified by dichotomized genotype. We report the results as the percentage increase in NHL risk for each 10% increase in the exposure level, which was calculated as exp(β × ln[1.1]). We calculated the P for interaction based on the log-transformed continuous variable for each risk factor and the scored variable for dichotomized genotype.
Termite treatment before 1988 was analyzed as a categorical variable. The referent group comprised persons whose homes had never been treated for termites or had been treated only after 1988, and the exposed group was composed of persons with at least one termite treatment before 1988. We excluded from the analysis people who responded “don't know” for all of their homes, people who could not recall whether one or more of their homes had been treated for termites (but knew that the others had not), and people who could not remember when the treatments occurred. We used unconditional logistic regression models to calculate ORs and 95% confidence intervals (CIs) for the joint effect of termite treatment and dichotomized genes. We calculated P for interaction based on the scored variable for termite treatment and for dichotomized genotype.
All risk estimates were adjusted for the study design variables sex, age (<45, 45-64, ≥65 years), race (white, other/unknown), education (<12, 12-15, >15 years), and study center. Logistic regression models were conducted using SAS 9.13 (SAS Institute, Cary, NC).
We used the following criteria to identify findings of possible significance. For PCB180, we first identified SNPs for which (1) PCB180 levels in both plasma and carpet dust were significantly associated at P less than .05 with NHL risk in the same genotype, or (2) PCB 180 levels in either plasma or carpet dust were highly significantly associated at P values less than .01 with NHL risk for one of the genotypes. Similarly, for α-chlordane, which was also assessed 2 ways (carpet dust measurements and termiticide use), we identified SNPs for which both measures were significantly associated at P values less than .05 with NHL risk in the same genotype, or either measure was highly significantly associated at P values less than .01 with NHL risk in one of the genotypes. We further imposed the criterion that the genotype conferring increased NHL risk must be consistent for both measures of the same exposure. For plasma TEQ, we selected SNPs for which one genotype was significantly associated at P values less than .01 with NHL risk.
Sixteen SNPs met at least one of these criteria (see Tables S1,S2, available on the Blood website; see the Supplemental Tables link at the top of the online article); results for the remaining 45 SNPs are shown in Tables S3 and S4. To assist in interpretation of these findings, we summarize the main effects of the exposures among cases and controls with genotype data in Table 2.
Results
Among the 16 SNPs meeting the stated criteria, IFNG was noteworthy for its consistent effect on NHL risk across all of the exposures (Figure 1, Table 3). Specifically, we observed significant elevations in risk from PCB180 exposure among IFNG (C−1615T) TT homozygotes for measurements in both plasma (16.9% increase in risk for each 10% increase in pg/g lipid) and carpet dust (1.2% increase in risk for each 10% increase in ng/g), from plasma TEQ (19.2% increase in risk), and from α-chlordane in carpet dust (1.1% increase in risk); for termite treatment there was a 50% elevation in risk that was of borderline statistical significance. In contrast, no significant elevations in risk for any of the exposures evaluated were observed among carriers of the variant allele. No statistically significant interactions were observed between genotype and these exposures for this SNP.
IL4 also showed evidence of a broad effect across exposures. The IL4 5′-untranslated region (5′-UTR), Ex1−168C>T polymorphism interacted significantly with plasma TEQ, with significantly elevated risk observed only among persons with the more common genotype (CC; 12.5% increase in risk for each 10% increase in pg/g lipid of TEQ). CC homozygotes, but not those with the CT or TT genotype, had significantly increased NHL risk from PCB180 exposure, regardless of whether exposure was evaluated in plasma or carpet dust. The associations between NHL risk and both measures of chlordane exposure were also higher among CC homozygotes; the association for termite treatment was of borderline statistical significance.
Associations specific to PCB180, as measured in both plasma and carpet dust, and risk of NHL were limited to the same genotype for 3 additional genes: IL16 (3′-UTR, Ex22+871A>G), IL8 (T−251A), and IL10 (A−1082G; Figure 2). The IL16 3′-UTR, Ex22+871A>G SNP significantly modified the association between plasma PCB180 levels and NHL risk, with elevated risk restricted to the more common genotype (AA; 15.0% increase in risk for each 10% increase in pg/g lipid of PCB180); PCB180 levels in carpet dust were also significantly associated with NHL risk only for this genotype (1.1% increase), although the interaction was not significant. We also observed significant elevations in risk from PCB180 exposure only among IL8 (T−251A) TT homozygotes for measurements in both plasma (28.9% increase in risk for each 10% increase in pg/g lipid) and carpet dust (1.4% increase in risk for each 10% increase in ng/g). Finally, only those persons carrying a variant allele for the IL10 (A−1082G) polymorphism had significantly elevated risk from PCB180 exposure, regardless of whether exposure was measured in plasma or carpet dust.
For each polymorphism listed in Table 3, we also examined potential interactions between genotype and plasma levels of PCB156 and PCB194, which were significantly associated with NHL risk in the main study4 but were not measured in carpet dust. For each SNP, the genotype that had a significant association for PCB180 levels in plasma and carpet dust also had a significant association for PCB194, with significant interactions for IL8 (T−251A; P = .025) and IL16 (3′-UTR, Ex22+871A>G; P = .001; data not shown). For all of the SNPs except IL10 (A−1082G), significant associations were also observed in the same genotype for plasma PCB156, with significant interactions for IFNG (C−1615T; P = .007) and IL16 (3′-UTR, Ex22+871A>G; P = .030; data not shown).
Discussion
We examined potential interactions between exposure to organochlorine compounds and SNPs in genes involved in immune and inflammatory response. Our analysis suggests that IFNG and IL4 broadly affect the relation between organochlorine exposure and NHL risk, whereas the effects of IL16, IL8, and IL10 appear to be limited to PCB180. Overall, our findings suggest that organochlorines may affect NHL risk at least in part through biologic mechanisms involving immune effects.
Exposure to PCBs, dioxins, and furans can lead to adverse effects on immune, reproductive, neurobehavioral, and endocrine functions.17,31,32 Their effects on the immune system are of particular interest because immune alteration is a known risk factor for NHL. Experiments with laboratory animals and nonhuman primates have shown that dioxin exposure can lead to cellular and humoral immune suppression, increased susceptibility to infectious disease, thymus atrophy, and depressed antibody and lymphoproliferative responses.15 Exposure to PCBs and dioxins has been found to diminish host resistance to infection, resulting in increased disease incidence and severity.15 Prenatal exposure to these compounds was associated with changes in the T-cell lymphocyte population in Dutch infants.16 In Dutch preschool children, PCB body burden was associated with recurrent middle-ear infections and chickenpox.17 In Flemish adolescents, biomarkers of internal exposure to dioxin-like compounds were related to biomarkers of immune status, and exposure to these compounds was associated with a lower prevalence of allergic diseases15 ; interestingly, a history of allergies was associated with a reduced risk of NHL in the NCI-SEER study.33
Data also indicate that chlordane and its constituents have effects on the immune system. In an experimental mouse model, Tryphonas et al14 found that cis-nonachlor increased IgG levels, although chlordane itself had no effect. Both cis-nonachlor and trans-nonachlor led to an increased susceptibility for mice to bacterial infection, suggesting that the chemicals are immunosuppressing. In humans, exposure to chlordane was associated with decreased lymphocyte responsiveness and the frequent presence of autoantibodies.12
IL4 and IL10 were previously reported to be associated with risk of NHL or its subtypes.19,–21 Both play roles in one or more processes associated with lymphocyte development and immune function. IL-4 regulates B-cell proliferation and immunoglobulin (Ig) class switching,34 and it is a key component in the induction of the Th2 lymphocyte phenotype and the down-regulation of the Th1 lymphocyte phenotype.35 There is compelling evidence that variants in IL10 play an important role in the cause of NHL.19,–21 The IL10 A-1082G polymorphism is believed to result in higher production of IL-10,36,37 and it has been postulated that IL-10, as a B-cell stimulatory cytokine, may promote lymphomagenesis.38 In support of this hypothesis is evidence that a putatively high IL-10 expressing genotype (−592 CC) is associated with increased risk of AIDS-associated NHL.20,38 In addition, elevated levels of IL-10 in the vitreous have been correlated with primary intraocular lymphoma.19,39 We are unaware of previous reports of associations between NHL risk and IFNG, IL16, or IL8 polymorphisms.
Interaction between organochlorine exposure and immune gene variation is plausible, given the observed effects of organochlorines on the immune system, the known relation between immune alteration and NHL risk, and the observed associations between common variants in immune and inflammatory response genes and NHL risk. We believe that our findings are notable but require replication in a large independent effort before further speculation about the potential mechanisms involved in NHL risk.
A strength of our study is the systematic approach taken to evaluate the joint effects of exposure to organochlorines and a wide variety of genetic polymorphisms. Case ascertainment using SEER tumor registries ensured that few were missed, and the use of population controls minimized biases typically associated with hospital-based case-control studies. Our multifaceted approach toward exposure assessment (ie, questionnaire, blood samples, carpet dust samples) allowed us to evaluate the consistency of our results across various exposure assessment methods.
A limitation of the study is the low statistical power for estimating the combined effects of genetic polymorphisms and some of the less common risk factors. Another shortcoming is that genetic variability of inflammatory/immune pathways is not completely explained by the genes and SNPs studied here. Other limitations include low response rates (typical of current population-based studies), possible misclassification in the exposure data, possible confounding, and the possibility that multiple comparisons testing produced false-positive findings, although this is unlikely for the SNPs we have highlighted because of their consistency across multiple exposures. Finally, it is possible that the SNPs producing the effects observed in the study are in linkage disequilibrium with other variants that could confer the observed biologic effect.
In summary, our exploratory results require replication in additional large studies and in pooled analyses. However, they provide evidence that organochlorine exposure may pose a differential risk of NHL to persons according to variability in genes involved in immune and inflammatory response.
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
We thank our study participants for their involvement. We thank Carol Haines (Westat) for development of study materials and procedures, for selection of younger controls, and for study coordination; Lonn Irish (IMS) for computer programming support; Carla Chorley (BBI Biotech Research Laboratories) for specimen handling; and Geoffrey Tobias for research assistance.
This work was supported by the Intramural Research Program of the National Institutes of Health at the NCI and by contracts with the NCI (NCI-SEER contracts N01-CN-67008, N01-CN-67010, N01-PC-67009, and N01-PC-65064; NCI-Southwest Research Institute contract N02-CP-19 114; and NCI-CDC interagency agreement Y1-CP-8028-02).
National Institutes of Health
Authorship
Contribution: P.H. was the principal investigator of the NCI-SEER NHL study; R.K.S., J.R.C., W.C., and S.D. obtained questionnaire and cancer outcome data from the 4 NCI-SEER registries; J.S.C. obtained dust sampling data; N.R. and A.J.D.R. obtained plasma organochlorine data; S.C. supervised genotyping at the NCI Core Genotyping Facility; J.S.C. performed the statistical analysis with input from S.S.W., N.R., P.H., L.M.M, and N.C.; and J.S.C., S.S.W., and N.R. drafted and revised the paper. All authors reviewed and approved the final manuscript.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Joanne S. Colt, National Cancer Institute, Occupational and Environmental Epidemiology Branch, 6120 Executive Boulevard, Room 8112, Bethesda, MD 20892-7240; e-mail: coltj@mail.nih.gov.