| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Epidemiology and Prevention |
1 Division of Cancer Epidemiology and Genetics and 2 Center for Cancer Research, National Cancer Institute, NIH, Department of Health and Human Services, Bethesda, Maryland; 3 School of Public Health, University of California, Berkeley, California; 4 Chinese Center for Disease Control and Prevention, Beijing, China; 5 School of Public Health, University of North Carolina, Chapel Hill, North Carolina; and 6 New York Blood Center, Clinical Services, White Plains, New York
Requests for reprints: Qing Lan, Occupational and Environmental Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, NIH, Department of Health and Human Services, MSC 7240, 6120 Executive Boulevard, EPS 8109, Bethesda, MD 20892-7240. Phone: 301-435-4706; Fax: 301-402-1819; E-mail: qingl{at}mail.nih.gov.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Interindividual variation to benzene-induced hematotoxicity has been observed in occupational studies, in which highly disparate responses have been noted among workers with similar benzene exposure (15, 16), suggesting a role for individual genetic susceptibility. The first reports to evaluate mechanisms of genetic susceptibility for benzene hematotoxicity identified genetic variations in enzymes responsible for benzene metabolism (e.g., CYP2E1, MPO, NQO1; refs. 3, 16). We hypothesized that another important source of susceptibility to benzene-induced hematotoxicity would be variation in selected cytokine, chemokine, and cellular adhesion genes because mounting evidence suggests that control of hematopoietic cell growth and differentiation, both during healthy conditions and in response to stress, is partly regulated by the cytokine, chemokine, and adhesion molecule networks (17, 18).
We recently reported the presence of hematotoxicity in Chinese workers exposed to low levels of benzene (<1 ppm, the U.S. occupational standard; ref. 3). Exposure to benzene and other solvents was extensively evaluated, and detailed information was collected on potential confounders including lifestyle, medication, and medical history. Here, we report the influence of single-nucleotide polymorphisms (SNPs) in hematopoiesis regulatory genes on benzene-induced hematotoxicity.
| Materials and Methods |
|---|
|
|
|---|
Interviews were carried out and blood samples were collected from 88 workers in June 2000 and from the remaining workers in May and June 2001 (28 subjects enrolled in the first year had repeat samples in the second year). A detailed questionnaire was administered requesting information on lifetime occupational and environmental exposures, lifestyle, medical history, and other potential confounding factors. Each subject also had a brief physical exam (3, 19) and provided a 29 mL blood sample. Complete blood cell counts and differentials were analyzed by a Beckman-Coulter T540 blood counter. The major lymphocyte subsets were measured by a Becton Dickinson FACSCalibur flow cytometer (software: SimulSET v3.1; refs. 3, 19). Colony-forming unit progenitor cells were cultured from peripheral blood from 29 workers selected for a wide range of measured benzene exposures (3).
Genotyping. We selected 40 SNPs in 20 genes coding for cytokines, chemokines, and adhesion molecules that play a role in hematopoiesis, based on the following criteria: minor allele frequencies of >5%, evidence of functionality in experimental or human studies, known association with disease, or to facilitate haplotype analysis (Table 1). Genotyping was done on an ABI 7900HT detection system using TaqMan end point reads as described on the website (20).7 Missing genotype data resulted from nonspecific failure of samples to amplify or to yield unambiguous genotype calls. Blind replicate samples were randomly interspersed throughout the study sample plates and showed intrasubject concordance rates of >99% for all assays.
|
2 test, with 1 degree of freedom. All genotypes were in Hardy-Weinberg equilibrium among the benzene-exposed workers except LTA (IVS1+90G>A; P = 0.016) and IL-16 (Ex22+889G>T; P = 0.001). Quality control data were rechecked to ensure the accuracy of these two assays. We tested the influence of each SNP on total WBC counts separately in workers exposed to benzene and among unexposed controls. Total WBC count was the key end point in this study because a decreased WBC count (under 4,000 cells/µL) measured repeatedly is the primary component of a diagnosis of benzene poisoning in China, which we have previously linked to risk for developing hematologic malignancies and related disorders among workers exposed to benzene (16).
For analysis of gene variants, the most prevalent homozygous genotype was considered the reference group. Data for minor homozygotes and heterozygotes were combined in the analysis if the former contained fewer than five subjects. Cell counts are presented in tables as unadjusted means (± SD). Linear regression was used to analyze the relationship between each genotype and the natural log (ln) of the total WBC count, adjusting for age (continuous variable), sex, current cigarette smoking status (yes/no), current alcohol consumption (yes/no), recent infections (yes/no), and body mass index (BMI), and in exposed workers and additionally for the ln mean air benzene exposure and ln mean air toluene exposure in the month before phlebotomy (3). The selection of potential confounding factors was based on the literature and associations observed in the current study population (3). Further adjustment for variants in MPO and NQO1, which we have previously reported influence WBC counts among exposed workers in this study (3), did not affect the associations reported here. Tests for gene-dosage effect (i.e., trend) were conducted by treating genotype as a continuous variable (0, 1, and 2 for the most prevalent homozygous, heterozygous, and variant homozygous genotypes, respectively). Gene-environment interactions were tested by adding a product term between the genotype (variant homozygous and heterozygous carriers combined versus most prevalent homozygous carriers) and benzene exposure (yes/no) into each model. Data from the 28 benzene-exposed workers studied in both years are treated independently by using generalized estimating equations to adjust for potential correlation between the repeated measurements (21). Results were negligibly affected when data from only the first or second year of study were used for these 28 subjects.
We controlled for multiple hypothesis testing by calculating the false discovery rate, using the Benjamini-Hochberg method (22), for results from the initial screening of the effect of each SNP on the total WBC count, the primary end point in our study. The false discovery rate is defined as the expected ratio of erroneous rejections of the null hypothesis to the total number of rejected hypotheses; we used a false discovery rate of 0.05 to report associations for which the proportion of expected false positives was <5%. False discovery rateadjusted P values were calculated for change in WBC counts for homozygous carriers of the rare versus common allele, as this provides the greatest potential contrast in effects across genotypes.
Those SNPs with false discovery rateadjusted P values of <0.05 that also showed significant gene-dosage effects (if testable) among workers exposed to benzene were further tested for their influence on specific WBC types, colony-forming unit granulocyte-macrophage (CFU-GM) progenitor cells, and colony-forming unit granulocyte-erythroid-macrophage-megakaryocyte (CFU-GEMM) progenitor cells. Negative binomial regression was used to evaluate each SNP association with CFU-GM colony counts, and logistic regression was used to test the association with CFU-GEMM colonies, as there were multiple subjects with no measurable colonies (3). The progenitor cell models were adjusted for only age, sex, and ln benzene to avoid overparameterizing a model with 29 observations, as other potential confounders were not significant. All P values presented are two sided and all analyses were carried out using SAS version 8.02 software (SAS Institute, Cary, NC).
When possible, a preliminary analysis of haplotype block structure was examined with HaploView8 using the four gamete rule with a minimum frequency of 0.005 for the fourth gamete. Haplotypes were estimated using the expectation-maximization algorithm (23). The association between blood cell counts and estimated haplotypes was assessed separately for controls and exposed subjects, adjusting for potential confounders using the HaploStats program (24) in R (Version 2.0.1; ref. 25).
| Results |
|---|
|
|
|---|
|
|
We conducted further analyses to evaluate the influence of the six SNPs shown in Table 3 on specific WBC populations in benzene-exposed workers (Table 4). The effect of SNPs in IL-1A, IL-10, and CSF3 was limited to granulocyte counts, with decreased cell counts observed for IL-1A and IL-10 and increased cell counts for CSF3. The IL-4 promoter SNP was significantly associated with decreased granulocyte and total lymphocyte counts; and the IL-12A variant was associated with decreased granulocytes, total lymphocyte count, and CD4+ and CD8+ T-cell subsets. Interestingly, the VCAM1 SNP displayed the broadest effect on WBC subtypes, with significant decreases for total lymphocyte count, CD4+ T cells, B cells, NK cells, and monocytes, and a borderline significant decrease for granulocytes (Table 4).
|
Because a number of SNPs in Table 4 influenced several mature WBC subtypes (e.g., granulocytes, CD4+ T cells, B cells, and monocytes), we carried out an exploratory analysis of the effect of each SNP on CFU-GM and CFU-GEMM progenitor cell counts, which were available from 29 workers exposed to benzene (the IL-1A variant could not be evaluated as there was only one carrier in this group). There were no significant effects observed for any SNP on CFU-GM progenitor cells. The VCAM1 variant was, however, significantly associated with a decrease in CFU-GEMM counts (Fig. 1A(i)) and the CSF3 SNP showed a borderline significant increase with CFU-GEMM counts (Fig. 1B(i)), consistent with the direction of their effect on total WBC counts (Fig. 1A(ii) and B(ii); Table 3).
|
| Discussion |
|---|
|
|
|---|
Benzene is metabolized to toxic compounds, which are known to cause bone marrow suppression, perhaps by complimentary mechanisms. Benzene can directly damage hematopoietic progenitor cells (58, 26), which could lead to apoptosis or decreased responsiveness to cytokines and cellular adhesion molecules. Alternatively, benzene toxicity to stromal cells or mature blood cells could disrupt the regulation of hematopoiesis, including hematopoietic commitment, maturation, or mobilization, through the network of cytokines, chemokines, and adhesion molecules (1114, 27). Accordingly, hematotoxic effects could be enhanced among individuals exposed to benzene who have genetic variants that alter key pathways that regulate hematopoiesis.
There is mounting evidence that benzene may alter the gene expression, production, or processing of several cytokines in vitro or in vivo. Both benzene and its metabolite hydroquinone were found to interact with signaling pathways regulated by granulocyte colony-stimulating factor (27). Benzoquinone and hydroquinone were shown to alter interleukin (IL)-1
processing, possibly by decreasing the content of the IL-1
processing enzyme calpain (11, 13, 28). Hydroquinone was reported to interrupt B-cell lymphopoiesis, possibly by inhibiting production of IL-4 by fibroblastic stromal cells (29), and was found to enhance IL-4 mRNA expression and IL-4 production in CD4+ T cells (30). Further, a recent report from our study population showed a highly significant increase in IL-4R mRNA expression in peripheral mononuclear cells from benzene-exposed workers (31). The influence of benzene on the production and function of IL-10, IL-12, or VCAM1, and how benzene might alter stem, progenitor, or stromal cell responses to these compounds, is not known.
Published data support the importance of these genes in hematopoiesis; for instance, IL-1, IL-4, IL-10, and IL-12 have each been implicated in myelopoiesis and lymphopoiesis (29, 3234). IL-12 and IL-4 play a role in the proliferation of T cells (34, 35) and VCAM1 has been shown to be critical for successful hematopoiesis in the bone marrow for cells of both myeloid and lymphoid lineage (17, 36, 37). It is also notable that several of our most significant findings were in genes that regulate the Th1-Th2 lymphocyte balance. In particular, we observed that variation in cytokines IL-4, IL-10, and IL-12A was each associated with decreased WBC. This raises the possibility that a shift in the balance of Th1-Th2 could have important consequences for the challenge of benzene to granulopoiesis and lymphopoiesis.
For three of the genes, IL-1A, IL-4, and IL-10, the variants studied have been reported to have functional consequences, either directly or as part of a haplotype. For instance, the IL-1A (889C>T) is in the 5' untranslated region, has been shown to produce a higher amount of IL-1A mRNA (38), and has been associated with autoimmune, degenerative, and infectious diseases (39, 40). The IL-4 promoter SNPs and haplotypes are particularly interesting because of emerging studies that indicate that variation in the expression of this gene could affect a range of diseases such as severity of infection with respiratory syncytial virus in young children, asthma, atopy, fungal infection with Candida albicans in leukemia patients, and inflammatory bowel disease (4144); moreover, this same promoter SNP has been shown to increase reporter gene expression in vitro (41, 42, 45). Lastly, proximal promoter IL-10 haplotypes have been shown to alter IL-10 secretion (4648) and a specific haplotype containing IL-10 (819T) has been associated with elevated IL-10 mRNA levels (47). IL-10 promoter haplotypes have been associated with progression of HIV infection, autoimmune diseases such as asthma, lupus erythematosus, graft-versus-host disease following marrow transplantation, and Alzheimer's disease (4952).
We detected associations between WBC counts and variants in three additional genes, VCAM1, CSF3, and IL-12A, but the specific function of these alleles is unknown. Recent reports found that SNPs in the VCAM1 gene, including the 1591C allele, were associated with increased risk of stroke (53, 54), suggesting that this variant, or a SNP in linkage disequilibrium, could be functional. The CSF3 (Ex4165C>T) and IL-12A (8685G>A) variants are in 3' untranslated regions, which could alter gene translation by several processes including ribosome binding, initiation, and elongation (55), but additional work is needed to determine their functionality.
In summary, this is the first report to show that common SNPs in cytokine genes and VCAM1 could be associated with hematotoxicity in workers exposed to benzene. Selected variants seemed to influence only granulocytes, whereas others altered cell types of both myeloid and lymphoid lineage, suggesting effects that could extend to earlier progenitor and possibly stem cells. Follow-up studies will require analysis with more comprehensive haplotypes containing the SNPs profiled in this study; in parallel, functional analysis should proceed in the laboratory to advance the identification of causal variants in benzene-induced hematotoxicity and to define specific mechanisms that explain their effects on peripheral WBC counts in the presence of benzene. At the same time, results of this study, although highly statistically significant, need to be confirmed in other populations exposed to benzene and extended to studies of benzene-related hematologic neoplasms and related disorders (2). Finally, our findings suggest that it may also be worthwhile to explore how these and related SNPs affect hematopoiesis and risk of hematologic disorders in people exposed to other environmental or pharmacologic hematotoxins and leukemogens.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Dr. Nilanjan Chatterjee for the false discovery rate calculation and the participants for taking part in this study.
| Footnotes |
|---|
7 http://snp500cancer.nci.nih.gov ![]()
8 http://www.broad.mit.edu/personal/jcbarret/haploview/ ![]()
Received 4/26/05. Revised 8/ 8/05. Accepted 8/ 9/05.
| References |
|---|
|
|
|---|
and 1ß to active cytokines by inhibition of the processing enzymes, calpain, and interleukin-1ß converting enzyme. Environ Health Perspect 1996;104 Suppl 6:12516.
processing enzyme calpain in murine bone marrow-derived macrophages after treatment with the benzene metabolite hydroquinone. Toxicol Lett 1994;74:17784.[CrossRef][Medline]
T mutation and rapid fractional excretion of chlorzoxazone. Cancer Res 1997;57:283942.
to mature cytokine in murine macrophages by hydroquinone and prevention of benzene-induced hematotoxicity in mice by IL-1
. Blood 1991;78:93844.
. Immunogenetics 2002;54:826.[CrossRef][Medline]
(889) promoter polymorphism is a risk factor for osteomyelitis. Am J Med Genet A 2003;119:1326.[Medline]
gene and Alzheimer's disease: a meta-analysis. Neurobiol Aging 2004;25:12938.[CrossRef][Medline]
polymorphisms. J Infect Dis 2003;187:211.[CrossRef][Medline]
chain genes confer susceptibility to asthma and atopy in a Caucasian population. Clin Exp Allergy 2003;33:11117.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
Q. Lan, L. Zhang, M. Shen, W. J. Jo, R. Vermeulen, G. Li, C. Vulpe, S. Lim, X. Ren, S. M. Rappaport, et al. Large-scale evaluation of candidate genes identifies associations between DNA repair and genomic maintenance and development of benzene hematotoxicity Carcinogenesis, January 1, 2009; 30(1): 50 - 58. [Abstract] [Full Text] [PDF] |
||||
![]() |
J-Q. He, K. Shumansky, J. E. Connett, N. R. Anthonisen, P. D. Pare, and A. J. Sandford Association of genetic variations in the CSF2 and CSF3 genes with lung function in smoking-induced COPD Eur. Respir. J., July 1, 2008; 32(1): 25 - 34. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Vineis and F. Perera Molecular Epidemiology and Biomarkers in Etiologic Cancer Research: The New in Light of the Old Am. Assoc. Cancer Res. Educ. Book, April 12, 2008; 2008(1): 547 - 567. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. S. Young and D. W. Kaufman The epidemiology of acquired aplastic anemia Haematologica, April 1, 2008; 93(4): 489 - 492. [Full Text] [PDF] |
||||
![]() |
P A Schulte The contributions of genetics and genomics to occupational safety and health Occup. Environ. Med., November 1, 2007; 64(11): 717 - 718. [Full Text] [PDF] |
||||
![]() |
P. Vineis and F. Perera Molecular Epidemiology and Biomarkers in Etiologic Cancer Research: The New in Light of the Old Cancer Epidemiol. Biomarkers Prev., October 1, 2007; 16(10): 1954 - 1965. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |