
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
1 School of Public Health, University of California, Berkeley, California and 2 Department of Epidemiology and Biostatistics, University of California, San Francisco, California
Requests for reprints: Martyn T. Smith, Division of Environmental Health Sciences, School of Public Health, University of California, 140 Earl Warren Hall #7360, Berkeley, CA 94720-7360. Phone: 510-642-8770; Fax: 510-642-0427. E-mail: martynts{at}berkeley.edu
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
2-fold increased risk estimate for NHL (7). We hypothesized that the chronic states of inflammation and altered immune responses associated with obesity that influence B- and T-cell function may contribute to development of NHL. We examined polymorphisms in the neuropeptide Y (NPY) and ghrelin (GHRL) genes to explore the association further between NHL and polymorphisms in genes related to energy homeostasis, obesity, and immune function in a population-based case-control study in the San Francisco Bay Area. NPY acts on the central nervous system as a potent appetite stimulator controlled by the feedback action of both leptin from adipose tissue and ghrelin from the stomach (8). Recent studies show that these opposing hormonal signals regulate the secretion of NPY in the hypothalamus and that the NPY network is the primary pathway involved in appetite stimulation (9). Ghrelin stimulates appetite during fasting directly through increased gut motility and through the release of NPY (ref. 9; Fig. 1), leading to increased adiposity and energy intake (10, 11). Conversely, the anorexigenic properties of leptin, a hormone released after eating, results in inhibition of NPY secretion (ref. 8; Fig. 1). NPY also regulates immune function via its expression in primary and secondary lymphoid organs and expression of its receptors on macrophages, natural killer cells, peripheral blood mononuclear cells, and B and T lymphocytes (12). A polymorphism in the NPY gene (1128T>C: rs16139) that results in an amino acid change (Leu > Pro) at codon 7 is associated with elevated levels of NPY, exercise-induced growth hormone secretion (13), and high serum cholesterol and low-density lipoprotein cholesterol levels (14). We examined this SNP in the present study along with three other SNPs distributed across the NPY gene (Fig. 2). Additionally, we examined four SNPs across the GHRL gene (Fig. 2), including a nonsynonymous SNP (GHRL 408 C>A: rs696217) or (Leu72Met) that is associated with lower BMI, fat mass, and abdominal visceral fat (10) and higher levels of circulating lipoprotein A associated with the wild-type genotype in diabetics (15). Ghrelin is the endogenous ligand for the growth hormone secretagogue receptor and potently stimulates the release of growth hormone (ref. 10; Fig. 1). Ghrelin levels are inversely correlated with BMI and therefore are decreased in obese individuals and increased in individuals with anorexia nervosa (11). This down-regulation of ghrelin levels in human obesity may be the result of elevated leptin and insulin levels (16).
|
|
| Materials and Methods |
|---|
|
|
|---|
65 years of age and that include
98% of U.S. residents ages
65 years. Controls were frequency matched to cases by 5-year age group, sex, and county. There were 1,593 eligible patients (72%) and 2,515 eligible controls (78%) who completed in-person interviews. Approximately 65% of patients who died before we could contact them were HIV-positive cases (23) and 12% (including eligibles and ineligibles) were HIV-negative cases. Therefore, those who died before contact were unlikely to have biased the results presented here that focus on HIV-negative participants only.
Biological Samples
Details of blood collection and specimen DNA isolation have been published (7) and will be presented briefly here. Blood specimens were obtained from 63% of eligible patients and 66% of controls for viral testing (individuals who had no chemotherapy within 3 months and no contraindications to venipuncture). However, because some specimen samples were depleted after testing that was part of the original laboratory analyses, we had fewer specimens with DNA available than originally were collected in the main study. The blood was processed using Ficoll-Paque separation and the lymphocytes were cryopreserved in liquid nitrogen. Coded specimens (458 cases, 812 controls), with all personal identifiers removed, were sent to the University of California at Berkeley laboratory for DNA isolation. DNA was isolated using a modified QIAamp DNA Blood Maxi Kit protocol (QIAgen, Inc., Santa Clarita, CA), and was quantified using PicoGreen dsDNA Quantitation kits (Molecular Probes, Eugene, OR) according to the manufacturers' specifications. All protocols and procedures were approved by the University of California at San Francisco Committee on Human Research and by the University of California at Berkeley Committee for the Protection of Human Subjects. All participants provided written informed consent before interview and a separate consent before venipuncture.
Histopathology
NHL histologic subtype and grade were re-reviewed by an expert pathologist for 97% of all NHL study patients and classified using the Working Formulation (NHL Classification Project). To approximate the Revised European-American Classification of Lymphoid Neoplasms in these analyses, Working Formulation diffuse large cell, and immunoblastic lymphoma were combined for the DLCL subtype and Working Formulation follicular small, mixed, and large cell lymphomas were combined for the follicular lymphoma subtype.
SNP Selection
SNPs were selected for the GHRL and NPY genes using the SNPper (http://snpper.chip.org/) and SNP500Cancer (http://snp500cancer.nci.nih.gov/) web sites and are described in Table 1. We attempted to cover the haplotype of each gene by choosing tag SNPs using Haploview (http://www.broad.mit.edu/personal/jcbarret/haploview) for haplotype analysis. All available Taqman assays designed by Applied Biosystems (Foster City, CA) were identified at http://www.appliedbiosystems.com. SNPs were chosen based on a minor allele frequency exceeding 5%, and preference was given to SNPs residing in the untranslated region or exonic coding regions. When no acceptable coding SNPs were found, intronic SNPs were used to ensure gene coverage.
|
|
2 statistics and for low-frequency alleles, an exact test using a Markov-chain method (Genepop v3.4). Odds ratios (OR) and 95% confidence intervals (95% CI) as estimates of the relative risk (hereafter referred to as risk) were computed from unconditional logistic regression. In genotype analyses, the wild-type category (chosen either as the most common homozygotic genotype or arbitrarily if the same) was the reference group and models were adjusted for age and sex. BMI was categorized as <25 kg/m2 (lean to reference range) versus
25 kg/m2 (overweight to obese) or <30 kg/m2 (nonobese) versus
30 kg/m2 (obese) to evaluate gene-environment interactions. A likelihood ratio test (comparing nested models with and without the relevant interaction terms) was used to test for statistical evidence of gene-environment interactions between BMI and NPY or GHRL SNP genotypes. However, the power to detect interactions for all NHL and by NHL subtypes is low after adjusting for multiple testing. For each set of analyses (e.g., all tests of interaction between all pairs of SNPs, including LEP 2548G>A, LEP 19A>G, and LEPR Q223R from our previous analyses, ref. 7, and their association with NHL), the false discovery rate was estimated to adjust for multiple comparisons (26). The false discovery rate is the expected proportion of false positives among a set of tests where the hypotheses were rejected. Haplotype frequencies for case and control groups were estimated and the ORs for common haplotypes were computed with the reference group defined as the haplotype carrying the most frequent alleles at each loci. The probability of each haplotype, for each participant, was determined using the E-M algorithm as implemented in the gap package (gc.em function) available as an add-on in the statistical package, R (27). HIV-positive participants were excluded and all analyses were restricted to White non-Hispanic participants to avoid potential bias due to population stratification (28). Results from all statistical analyses were considered significant for two-sided Ps
0.05 and were considered borderline associated for 0.05 < P
0.10. | Results |
|---|
|
|
|---|
Effect of NPY and GHRL Genotypes on NHL Risk
Figure 2 displays the genomic locations of the individual polymorphisms genotyped in this study. Aside from NPY 1128T>C (borderline result from exact test P = 0.05), all control genotype distributions were in Hardy-Weinberg equilibrium.
The NPY 485T>C, 1258G>A, and 5671C>T SNPs were found to be in total linkage disequilibrium (D' = 0.99), and there also was evidence of linkage disequilibrium between these SNPs and the low frequency nonsynonymous NPY 1128T>C (Leu > Pro) polymorphism (P < 0.01). In univariate analyses, homozygous variant genotypes for NPY 485T>C, 1258G>A, and 5671C>T were overrepresented in the NHL cases resulting in increased risk for NHL in NPY 485T>C (OR, 1.7; 95% CI, 1.1-2.5), 1258G>A (OR, 1.7; 95% CI, 1.1-2.5), and 5671C>T (OR, 1.9; 95% CI, 1.3-2.8; Table 2). In analyses stratified by NHL subtype, results were even stronger for follicular lymphoma in homozygous variants for NPY 485T>C (OR, 2.4; 95% CI, 1.3-4.4), NPY 1258G>A (OR, 2.2; 95% CI, 1.1-4.1), and NPY 5671C>T (OR, 2.4; 95% CI, 1.3-4.3). The variant NPY 1128C allele also was associated with an increased risk for follicular lymphoma (OR, 2.3; 95% CI, 1.1-4.9; Table 2). No associations between NPY SNPs and DLCL were observed in these analyses (Table 2).
Strong linkage disequilibrium was found between the nonsynonymous GHRL 408C>A SNP and both the 4427A>G (D' = 0.84) and 5179A>G SNPs (D' = 1.0). The GHRL 4427GG genotype was associated with a lower risk of all NHL, but confidence limits overlapped unity (OR, 0.77; 95% CI, 0.51-1.2). However, in analyses stratified by NHL subtype, the homozygous variant GHRL 4427AA genotype was associated with a reduced risk for DLCL (OR, 0.31; 95% CI, 0.13-0.74; Table 2), with heterozygotes having an intermediary reduced level of risk (OR, 0.71). Few participants were homozygous variant for GHRL 5179A>G, but a borderline reduced risk for DLCL was associated with this genotype. No associations between GHRL SNPs and follicular lymphoma were observed in these analyses.
We found no interaction between BMI and variant allele frequencies. We did find significant two-way SNP-SNP interactions between NPY 1128T>C and GHRL 9344G>A and between GHRL 5179G>A and LEP 2548G>A; however, none of these results were significant after adjusting for multiple comparisons using the false discovery rate (the false discovery rate in no set of tests of interaction was <0.45) indicating relatively weak evidence of statistical interaction in this data set.
Haplotype Analyses
Common haplotypes (>5%) frequency was estimated for GHRL and NPY using tagSNP implementation of the EM algorithm (Table 3). Five common haplotypes were predicted for GHRL. Haplotypes containing variant alleles for any of the GHRL SNPs were associated with reduced risk for NHL and DLCL, although all estimates were imprecise. However, using the model where HapA was compared with all other haplotypes, HapA was associated with NHL (OR, 1.8; 95% CI, 1.2-2.8). For NPY, due to the low frequency of the NPY 1128T>C SNP and the tight linkage disequilibrium across the gene, the haplotype analysis did not reveal any further relevant information beyond that found in the univariate analysis.
|
| Discussion |
|---|
|
|
|---|
65% decreased risk was observed for homozygous variant genotypes. These results may suggest alternative mechanisms in the etiology of follicular lymphoma and DLCL related to the actions of NPY and GHRL in immunoregulation. However, the associations in NPY and GHRL by NHL subtype should be interpreted with caution given the wide confidence intervals and number of participants in some genotype categories when stratified by disease entity. These results warrant exploration in larger data sets and in multiple studies.
Whereas there has been major interest in the actions of NPY as a neurotransmitter, there is increasing evidence of the role of NPY in immune modulation. In immune cells such as B cells, monocytes, and macrophages, NPY levels are up-regulated upon cell activation (29). In macrophages, NPY increases adhesion, chemotaxis, phagocytosis, and superoxide anion production. However, NPY also suppresses innate immunity by inhibiting natural killer cell activity (30). Studies show NPY-induced reduction in natural killer activity in tumor cells of young animals (31), and that elevated NPY levels associated with stress and depression lead to reduced natural killer cytotoxicity (32). Furthermore, NPY stimulates lymphocyte proliferation and the release of interleukin 4, interleukin 6, and tumor necrosis factor-
cytokines (12), whereas it inhibits IFN-
production (33). Previous studies have shown that the nonsynonymous NPY Leu7Pro polymorphism (1128T>C) is associated with increased NPY secretion, elevated cholesterol levels, enhanced angiogenesis, and lymphocyte proliferation (14). Thus, in the present study, the associations found in the four NPY SNPs may be driven by the functional SNP 1128T>C, given that the four SNPs are in close linkage disequilibrium. Whether elevated cholesterol is a risk factor for NHL remains to be explored further, but this would be consistent with our previous results (21) and those of others (34) that show an inverse relationship between cholesterol-lowering drugs and NHL.
Because NPY suppresses natural killer cell activity, the increased NPY secretion associated with the 1128T>C variant could hinder normal actions of the innate immune system, which is the first line of defense against viral and bacterial infections. NPY also stimulates estradiol release by human granulosa cells in a dose-dependent fashion (35). Previous studies have shown that estradiol promotes B-cell proliferation and T-cell suppression and inhibits apoptosis (36). These actions could enhance the growth of transformed B-cells that may favor lymphomagenesis.
Ghrelin is a peptide hormone that is recognized as an important regulator of growth hormone release and energy homeostasis. As the only known circulating orexigen, ghrelin exerts antagonistic effects on the leptin-induced decrease in food intake through activation of the NPY pathway. Interestingly, as modulators in immune function, a mutually antagonistic relationship exists between ghrelin and leptin. Ghrelin inhibits the proliferation of inflammatory cytokines such as interleukin 6, interleukin 1ß, and tumor necrosis factor-
, whereas leptin promotes proinflammatory cytokine release (37). In addition, ghrelin levels are decreased in obese individuals, whereas leptin resistance results in elevated leptin levels found among obese individuals. Hypothetically, low levels of ghrelin and high levels of leptin present in obese individuals may induce a chronic proinflammatory state, potentially increasing NHL risk. Previous studies have shown that the homozygous variant of GHRL 408C>A (Leu72Met) is associated with lower BMI, fat mass, and abdominal visceral fat (10), but we did not detect a similar association between GHRL SNPs and BMI in the present study. The small number of carriers of the variant 408C>A allele in our study population may have limited our ability to detect an association if one exists.
Strengths of this study included the population-based design that used cancer registry data to identify incident NHL patients shortly after diagnosis and centralized expert re-review of diagnostic pathology materials for 97% of cases. These data were supplemented by receipt of the Surveillance, Epidemiology, End Results abstracts to identify additional cases that may have been missed by rapid case ascertainment. Controls from the same population base as the cases were identified through random digit dial supplemented by random sampling of Healthcare Financing Administration lists that include
98% of U.S. residents
65 years. No proxy interviews were conducted. Bias associated with exclusion of potential cases due to death is a concern in studies of cancer. However,
65% of patients who died before we could contact them were HIV-positive cases (23) and therefore are unlikely to have biased these results that focused on HIV-negative participants only. There is a potential for bias associated with misclassification of BMI given that self-reported usual adult weight and height used to compute BMI are subject to misclassification. The National Health and Nutrition Examination Survey data showed that the prevalence of overweight U.S. adults is likely to be underestimated because both men and women tend to overestimate their height whereas men overestimate and women underestimate their weight, especially among those who are overweight or obese and among those >60 years old resulting in underestimates of overweight (38). If this misclassification was nondifferential then the estimates for BMI may be biased toward the null. Analyses were restricted to white non-Hispanic participants to control for the potential effects of population stratification (28). Due to few participants with homozygous variant genotypes, estimates associated with this genotype are unstable and may be due to chance, especially for histologic subtypes. In addition, there is some uncertainty in haplotypes from association studies because they must be inferred from unphased genotype data. However, based on a recent comparison of standard haplotype analysis techniques that showed that the E-M algorithm method used in our analyses performs well for uncomplicated haplotypes and haplotype risk estimates of
2 (39), the potential bias for these haplotype estimates is likely to be small. Despite concerns regarding small sample size in some analyses and multiple testing issues, results tended to be consistent and warrant replication and confirmation in larger studies.
In conclusion, this study provides further evidence of a relationship between NHL and genes associated with energy regulation and homeostasis and suggests that elevated NPY levels may promote immune dysfunction that could enhance the growth of transformed B-cells that may progress to lymphoma. Larger studies are needed to replicate these findings, to further explore potential interactions among NPY, GHRL, and LEP polymorphisms, and to clarify the interplay between these genes and environmental factors. Animal studies and functional studies of specific SNPs also are needed to identify the mechanistic pathways to explain these observations.
| Footnotes |
|---|
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.
Received 12/ 6/04; revised 2/15/05; accepted 3/ 3/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. F. Skibola, J. D. Curry, and A. Nieters Genetic susceptibility to lymphoma Haematologica, July 1, 2007; 92(7): 960 - 969. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Wagner, K. Hemminki, E. Grzybowska, R. Klaes, B. Burwinkel, P. Bugert, R. K. Schmutzler, B. Wappenschmidt, D. Butkiewicz, J. Pamula, et al. Polymorphisms in genes involved in GH1 release and their association with breast cancer risk Carcinogenesis, September 1, 2006; 27(9): 1867 - 1875. [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 | Cell Growth & Differentiation |