The Interaction between Pesticide Use and Genetic Variants Involved in Lipid Metabolism on Prostate Cancer Risk.

Background. Lipid metabolism processes have been implicated in prostate carcinogenesis. Since several pesticides are lipophilic or are metabolized via lipid-related mechanisms, they may interact with variants of genes in the lipid metabolism pathway. Methods. In a nested case-control study of 776 cases and 1444 controls from the Agricultural Health Study (AHS), a prospective cohort study of pesticide applicators, we examined the interactions between 39 pesticides (none, low, and high exposure) and 220 single nucleotide polymorphisms (SNPs) in 59 genes. The false discovery rate (FDR) was used to account for multiple comparisons. Results. We found 17 interactions that displayed a significant monotonic increase in prostate cancer risk with pesticide exposure in one genotype and no significant association in the other genotype. The most noteworthy association was for ALOXE3 rs3027208 and terbufos, such that men carrying the T allele who were low users had an OR of 1.86 (95% CI = 1.16-2.99) and high users an OR of 2.00 (95% CI = 1.28-3.15) compared to those with no use of terbufos, while men carrying the CC genotype did not exhibit a significant association. Conclusion. Genetic variation in lipid metabolism genes may modify pesticide associations with prostate cancer; however our results require replication.

1. Background

Previous studies of prostate cancer have shown elevated rates in agricultural and pesticide manufacturing populations [1, 2]. In the Agricultural Health Study (AHS), a significant excess of prostate cancer was observed among private and commercial pesticide applicators compared to the general population [3, 4]. Also, use of pesticides, such as phorate [5], fonofos [6], butylate [7], and coumaphos [8], has been linked with an increased risk of prostate cancer among AHS participants with a family history of prostate cancer. We conducted a prostate cancer nested case-control study within the AHS to examine interactions between prespecified genetic pathways and pesticide exposure. Recent findings from this study have identified significant pesticide interactions for several genetic variants in the 8q24 region [9], xenobiotic metabolism pathway [10], and DNA repair pathways [11, 12]. These studies help elucidate exposure-effect associations by identifying potentially susceptible subgroups. This allows us to better understand potential carcinogenic hazards and furthers public health research on the human health effects of pesticides.

In this analysis, we evaluated single nucleotide polymorphisms (SNPs) in genes related to lipid metabolism since there is evidence that lipogenic enzymes and serum lipids may play a role in prostate carcinogenesis via inflammatory, oxidative stress and insulin-resistance mechanisms [13, 14]. Cellular and molecular studies have reported an overexpression of lipogenic enzymes, such as fatty acid synthase and acetyl-CoA carboxylase, and increased lipogenesis in prostate cancer cells [15, 16]. Genetic variants involved in fatty acid or lipid metabolism, which have been linked to altered expression of lipogenic enzymes and serum lipid levels [16-18], have also been associated with an excess risk of prostate cancer [19-21]. Furthermore, several pesticides are lipophilic or are metabolized via a lipid-related mechanism [22, 23], suggesting that variants of genes involved in lipid metabolism may interact with certain pesticides to increase the risk of prostate cancer. We examined potential gene-environment interactions between SNPs in genes involved in lipid metabolism and pesticide exposure in a nested case-control study within the AHS.

2. Methods

2.1. Study Subjects

A detailed description of the AHS cohort has been published [24]. Briefly, 57, 310, or 82%, of pesticide applicators seeking pesticide licensing in Iowa and North Carolina were enrolled in the study between December 13, 1993 and December 31, 1997. Participants were licensed private pesticide applicators (mostly farmers) residing in Iowa and North Carolina and commercial applicators residing in Iowa. To be eligible for this nested case-control study, participants had to provide a mouth-wash rinse buccal cell sample, and could not have had a previous cancer diagnosis except for non-melanoma skin cancer. Cancer diagnoses were identified by linkage with population-based state cancer registries. Cases were males diagnosed with incident prostate cancer between 1993 and 2004, and controls were males alive at the time of case diagnosis with no cancer diagnosis, except nonmelanoma skin cancer at the time of selection. Eligible controls were randomly selected and frequency matched 2 : 1 to cases by date of birth (±1 year). All participants were self-identified as Caucasian and confirmed to be of European ancestry using population substructure evaluation and principal component analysis as previously described [9].

2.2. Pesticide Use

Information on pesticide use was obtained from two self-administered questionnaires completed at enrollment (http://aghealth.nci.nih.gov/questionnaires.html). These questionnaires collected comprehensive data on lifetime use of 50 pesticides. Participants were asked how many years (1 year or less, 2–5, 6–10, 11–20, 21–30, or more than 30 years) and how many days (less than 5, 5–9, 10–19, 20–39, 40–59, 60–150, or more than 150 days) in an average year they applied each chemical. Pesticide use was quantified as ever/never use, lifetime exposure days (years of use × days per year), and intensity weighted lifetime exposure days (lifetime exposure days × intensity score). The intensity score was computed from an algorithm that took into account exposure-modifying factors of how pesticides were used and applied [25].

2.3. SNP Selection and Genotyping

Germline DNA was extracted from buccal cells using Autopure (QIAGEN) at Sera Care, Frederick, MD. Genotyping was conducted with a custom Infinium BeadChip Assays (iSelect)from Illumina Inc., at the NCI Core Genotyping Facility (http://cgf.nci.nih.gov/operations/multiplex-genotyping.html). Common SNPs (minor allele frequency ≥ 5%) with an r2 threshold >0.80 were selected from the HapMap Project population using the Tagzilla software package that implements a modified tag SNP selection method described by Carlson et al. [26]. The Illumina chip included 26, 512 SNPs, of which 1,858 tag SNPs in 85 genes were related to the lipid metabolism pathway. Blinded duplicate samples (2%) were genotyped; the concordance of these samples ranged from 96–100%. The overall genotyping rate for lipid metabolism SNPs was between 97% and 100%. Of the samples genotyped, 108 were removed due to insufficient or poor DNA quality (14 cases, 6 controls) or <90% completion rate (47 cases, 41 controls), leaving 2,220 (776 cases, 1,444 controls) for statistical analyses.

2.4. Statistical Analysis

We compared the distributions of selected characteristics (age (<40, 40–49, 50–59, 60–69, ≥70), state of residence (Iowa, North Carolina), applicator type (private, commercial), and first-degree family history of prostate cancer (no, yes)) between cases and controls using a chi-square test.

Associations between the 1,895 lipid metabolism SNPs in 85 genes and prostate cancer risk were assessed using PLINK software [27] using a log-additive genetic model to calculate odd ratios (OR) and 95% confidence intervals (CI). Variants with a Hardy Weinberg equilibrium P value >0.001, a minor allele frequency ≥0.1, and a P trend for the additive effect of each allele <0.2 were carried forward for interaction analyses, which totaled 220 SNPs in 59 genes. Pairwise linkage disequilibrium (e.g., D′ and r2) between SNPs was evaluated using Haploview 4.2 [28].

The statistical interactions between pesticide use and genetic variants were assessed using the likelihood ratio test in unconditional logistic regression models adjusting for age and state. We examined 39 pesticides that were ever used by at least 10% of the study population. Exposure for each pesticide was categorized as nonexposed, low, and high exposure, with low- and high-exposure categories based on the median intensity-weighted lifetime days. Genetic variants were coded using a dominant genetic model due to the potential for small sample sizes of variant alleles. We calculated ORs and 95% CIs for prostate cancer risk and pesticide exposure stratified by genotype (homozygous wild-type, heterozygous and homozygous variant) using unconditional logistic regression models adjusting for age and state. In essence, we adopted a two-stage interaction approach [29], where we first set a cut point on the SNPs we included in the interaction analyses, and then cut point on the P interaction value to increase statistical power.

We used the false discovery rate (FDR) [30] to account for multiple comparisons for the main effect associations and the interactions. For the main effect associations, FDR was calculated considering all SNPs within each gene. For the interaction analyses, we calculated FDR in two ways: (1) accounting for the 220 SNPs we evaluated per pesticide, which is the most stringent approach, and (2) accounting for the number of SNPs within each gene per pesticide. Since our a priori hypothesis was an increased risk of cancer with exposure, we focused on statistically significant monotonically increasing ORs with increasing pesticide exposure in one genotype and no significant association in the other genotype (quantitative interaction), versus an interaction with increasing pesticide exposure in one genotype and a decreasing pesticide exposure in the other genotype (qualitative interaction), since the former are considered to be more biologically plausible and less due to chance [31]. Unless otherwise indicated, analyses were conducted using SAS, version 9.1 (SAS Institute, Cary, NC) and AHS data release version P1REL0712.04.

3. Results

Characteristics of the study population (see Supplementary Table 1 in Supplementary Material available online at doi:10.1155/2012/358076) and associations between intensity-weighted lifetime days of pesticide use and prostate cancer (Supplementary Table 2) were previously published [9-12]. Briefly, cases and controls were similar in terms of age, state of residence and applicator type, but cases were more apt to have a family history of prostate cancer.

Of the 1,895 lipid metabolism SNPs in 85 genes we examined, 20 SNPs in 8 genes (ABCG5, ALOXE3, AMACR, FOXA3, LRP2, PPARA, TBXAS1, UCP3) were significant using a gene-based FDR <0.05 (Table 1). The two most significant SNP-cancer associations were for two highly correlated variants in ABCG5, rs1835815, and rs2278356 (r2 = 0.89, FDR P main effect = 0.01).

Table 1

Prostate cancer risk in relation to lipid metabolism variants with FDR P trend <0.05 (sorted by gene and FDR P  trend).

Gene					Minor allele
	Chromosome	SNP	SNP region	Minor allele	Frequency1	OR2	(95% CI)2	P trend2	FDR  P trend2
ABCG5 (ATP-binding cassette subfamily G member 5)	2	rs1835815rs2278356rs4952683rs10205816rs4148214rs6709904	IVS10−1132C>TEx13−265T>GIVS12−44A>GIVS11−2527T>CIVS4+43T>CIVS6+317A>G	TCGGTG	0.280.310.190.260.460.10	1.221.221.211.160.871.23	(1.07–1.39) (1.07–1.40) (1.04–1.41) (1.01–1.33) (0.77–0.99) (1.02–1.49)	0.0040.0030.020.040.030.03	0.010.010.040.040.040.04
ALOXE3 (epidermis-type lipoxygenase 3)	17	rs6503086	−15319C>A	G	0.44	0.83	(0.73–0.94)	0.003	0.02
AMACR (alpha-methylacyl-CoA racemase)	5	rs299602rs6863657	IVS4−1798C>TIVS4+4012C>T	AA	0.310.11	1.170.77	(1.03–1.34) (0.62–0.96)	0.020.02	0.040.04
FOXA3 (hepatocyte nuclear factor 3-gamma)	19	rs11669442	IVS1−776T>C	C	0.31	1.18	(1.03–1.35)	0.01	0.04
LRP2 (low-density lipoprotein-related protein 2)	2	rs10490131rs2239595rs2247506rs2268365rs2268372rs2673164	IVS46−765T>CIVS30+70C>TIVS4+23A>GIVS27−372A>GIVS38−824A>TIVS4−2306G>A	GACCAC	0.130.110.190.110.180.28	1.241.331.211.271.220.83	(1.05–1.47) (1.11–1.61) (1.04–1.41) (1.05–1.53) (1.04–1.42) (0.72–0.95)	0.010.0020.010.010.010.01	0.040.040.040.040.040.04
PPARA (peroxisome proliferator-activated receptor alpha)	22	rs4253655rs6007662	IVS3+1780G>AIVS10+2685A>G	AG	0.150.25	1.241.20	(1.05–1.46) (1.04–1.38)	0.010.01	0.020.02
TBXAS1 (thromboxane A synthase 1)	7	rs2284212	IVS9−11575G>A	A	0.28	1.22	(1.07–1.40)	0.004	0.04
UCP3 (mitochondrial uncoupling protein 3)	11	rs1726764	−13576A>G	C	0.47	1.17	(1.04–1.33)	0.01	0.03

1Frequency of variant allele among controls.

2Risk estimate for effect per allele assuming log-additive model adjusting for age and state.

Among the 220 SNPs in 59 genes and 39 pesticides examined in interaction analyses, we found 16 interactions significant at FDR P interaction <0.2 accounting for all 220 SNPs per pesticide. All of these interactions had a significant association between pesticide use and prostate cancer in at least one genotype group (P trend < 0.05), 9 also had an increased risk of prostate cancer (OR trend > 1.0); and three were monotonic associations (i.e., increasing cancer risk with increasing pesticide use) and quantitative interactions (Table 2). These associations involved one pesticide, terbufos, and three SNPs in three genes (ADIPOR1, ALOXE3, SEC14L2).

Table 2

Associations between intensity-weighted lifetime days of pesticide use and prostate cancer risk by lipid metabolism variants where FDR P interaction <0.21.

Pesticide	Gene	SNP	SNP region	Genotype	Nonexposed		Low exposed				High exposed				Exposure trend			P interaction2	FDR P interaction1,2
					Case/control		Case/control	OR2	(95% CI)2	P 2	Case/control	OR2	(95% CI)2	P 2	OR2	(95% CI)2	P 2
	ALOXE3 (epidermis-type lipoxygenase 3)	rs3027208	IVS7+369G>A	CC CT+TT	307/576 99/227	REFREF	102/195 43/55	0.991.86	(0.74–1.31) (1.16–2.99)	0.930.01	81/188 50/60	0.812.00	(0.60–1.10) (1.28–3.15)	0.180.00	0.91 1.45	(0.79–1.05) (1.16–1.81)	0.220.001	0.001	0.08
Terbufos	SEC14L2 (SEC14-like protein 2)	rs2072159	IVS2−93G>A	GG AG+AA	201/457 204/342	REFREF	83/13961/108	1.380.97	(0.99–1.92) (0.67–1.40)	0.050.88	82/132 47/115	1.430.70	(1.03–1.99) (0.48–1.03)	0.030.07	1.21 0.86	(1.03–1.42) (0.71–1.03)	0.020.10	0.003	0.15
	ADIPOR1 (adiponectin receptor 1)	rsl2733285	IVS1−1748G>A	CC CT+TT	183/421 223/372	REFREF	74/123 70/126	1.390.96	(0.98–1.96) (0.68–1.35)	0.070.80	80/131 51/117	1.430.75	(1.02–1.99) (0.52–1.10)	0.040.14	1.21 0.88	(1.03–1.43) (0.74–1.05)	0.02 0.16	0.01	0.18

1Accounts for 220 SNPs included in interaction analyses.

2Risk estimates adjusted for age and state.

Using a less stringent method for accounting for multiple comparisons, we found 139 interactions significant using a gene-based FDR for each of the 39 pesticides examined (FDR P interaction < 0.05). Of these, 116 had significant associations between pesticide use and prostate cancer in at least one genotype group (P trend < 0.05); 31 resulted in increased prostate cancer risk (OR trend > 1.0) (Supplemental Table 4), and 17 were monotonic associations and quantitative interactions (Table 3). These 17 associations involved seven pesticides (atrazine, carbofuran, EPTC, fonofos, glyphosate, petroleum oil/distillate, and terbufos) and 15 SNPs in 11 genes (ABCG8, ACSL5, ADIPOR1, ALOX5, ALOX5AP, ALOX15, ALOXE3, A4GALT, LDLR, LRP1, SEC14L2). The findings from this method are inclusive of those using the more stringent method.

Table 3

Associations between intensity-weighted lifetime days of pesticide use and prostate cancer risk by lipid metabolism variants where FDR P interaction <0.051.

Pesticide	Gene	SNP	SNP region	Genotype	Non-exposed		Low exposed				High exposed				Exposure trend			P interaction2	FDR P interaction1,2
					Case/control		Case/control	OR2	(95% CI)2	P 2	Case/control	OR2	(95% CI)22	P 2	OR2	(95% CI)2	P 2
Insecticides
Terbufos	ALOXE3 (epidermis-type lip oxygenase 3)	rs3027208	IVS7+369G>A	CC CT+TT	307/576 99/227	REFREF	102/195 43/55	0.99 1.86	(0.74–1.31) (1.16–2.99)	0.93 0.01	81/188 50/60	0.812.00	(0.60–1.10) (1.28–3.15)	0.180.00	0.91 1.45	(0.79–1.05) (1.16–1.81)	0.22 0.001	0.001	0.01
	SEC14L2 (SEC14-like protein 2)	rs2072159	IVS2−93G>A	GG AG+AA	201/457 204/342	REFREF	83/139 61/108	1.38 0.97	(0.99–1.92) (0.67–1.40)	0.05 0.88	82/132 47/115	1.430.70	(1.03–1.99) (0.48–1.03)	0.03 0.07	1.21 0.86	(1.03–1.42) (0.71–1.03)	0.02 0.10	0.003	0.01
	ADIPOR1 (adiponectin receptor 1)	rs12733285	IVS1−1748G>A	CC CT+TT	183/421 223/372	REFREF	74/123 70/126	1.39 0.96	(0.98–1.96) (0.68–1.35)	0.07 0.80	80/131 51/117	1.430.75	(1.02–1.99) (0.52–1.10)	0.040.14	1.21 0.88	(1.03–1.43) (0.74–1.05)	0.02 0.16	0.01	0.02
Fonofos	ABCG8 (ATP-binding cassette subfamily G member 8)	rs4953028	IVS10+205G>A	GG AG+AA	141/345 368/644	REFREF	26/38 59/120	1.73 0.89	(0.99–3.00) (0.63–1.25)	0.05 0.50	36/46 56/106	1.940.95	(1.17–3.20) (0.66–1.36)	0.01 0.77	1.43 0.96	(1.12–1.82) (0.81–1.14)	0.004 0.63	0.005	0.02
Carbofuran	A4GALT (lactosylceramide 4-alpha-galactosyltransferase)	rs8136914	IVS1+4129T>C	AA AG+GG	210/477 222/378	REFREF	69/128 54/93	1.22 1.01	(0.87–1.71) (0.69–1.47)	0.25 0.96	72/115 48/107	1.430.76	(1.02–2.00) (0.52–1.12)	0.04 0.17	1.20 0.89	(1.02–1.41) (0.75–1.07)	0.03 0.22	0.02	0.03
	LDLR (low-density lipoprotein receptor)	rs8110695	IVS1−4369T>A	TT AT+AA	285/509 145/342	REFREF	76/132 47/92	1.04 1.17	(0.75–1.43) (0.78–1.76)	0.82 0.44	64/138 55/83	0.831.58	(0.60–1.16) (1.06–2.33)	0.27 0.02	0.93 1.24	(0.80–1.09) (1.03–1.51)	0.36 0.02	0.02	0.04
Herbicides
EPTC	ALOX15 (arachidonate 15-lipoxygenase)	rs916055	Ex14+241T>C	AA AG+GG	242/433 274/592	REFREF	34/45 45/70	1.38 1.47	(0.85–2.23) (0.98–2.22)	0.19 0.06	18/58 42/58	0.561.63	(0.32–0.99) (1.06–2.50)	0.05 0.02	0.85 1.31	(0.67–1.09) (1.07–1.60)	0.19 0.01	0.01	0.01
Petroleum iil/distillates	ALOX5 (arachidonate 5-lipoxygenase)	rs7099684	IVS12−61T>C	TT AT+AA	281/602 200/346	REFREF	33/63 19/36	1.19 0.91	(0.76–1.86) (0.51–1.63)	0.45 0.76	43/60 14/41	1.560.60	(1.03–2.38) (0.32–1.13)	0.04 0.12	1.24 0.80	(1.02–1.51) (0.60–1.07)	0.03 0.13	0.02	0.02
	LRP1 (low-density lipoprotein receptor-related protein 1)	rs1800159	IVS62+86A>G	GG AG+AA	229/477 258/484	REFREF	18/51 34/52	0.74 1.26	(0.42–1.31) (0.80–2.00)	0.31 0.32	20/58 41/45	0.721.76	(0.42–1.23) (1.12–2.77)	0.230.01	0.83 1.32	(0.65–1.07) (1.07–1.63)	0.15 0.01	0.01	0.03
		rs7956957	IVS78−51C>G	GG CG+CC	209/417 279/547	REFREF	14/44 38/59	0.66 1.32	(0.35–1.23) (0.85–2.03)	0.19 0.22	17/46 44/57	0.751.57	(0.42–1.34) (1.03–2.40)	0.33 0.04	0.83 1.26	(0.63–1.09) (1.04–1.54)	0.17 0.02	0.02	0.03
		rs6581128	IVS42−916C>A	CC AC+AA	227/452 261/509	REFREF	21/49 31/54	0.87 1.17	(0.51–1.49) (0.73–1.87)	0.61 0.52	18/51 43/50	0.701.73	(0.40–1.23) (1.12–2.68)	0.220.01	0.84 1.29	(0.65–1.09) (1.05–1.59)	0.20 0.01	0.02	0.03
		rs7978567	IVS67−220C>T	CC CT+TT	219/437 268/527	REFREF	16/48 36/55	0.69 1.34	(0.38–1.25) (0.86–2.10)	0.22 0.20	20/49 41/54	0.831.56	(0.48–1.43) (1.01–2.41)	0.49 0.05	0.87 1.26	(0.67–1.13) (1.03–1.55)	0.29 0.02	0.03	0.04
	ACSL5 (acyl-CoA synthetase long hain)	rs12098576	IVS6−530C>A	GG GT+TT	379/715 109/249	REFREF	35/71 17/32	0.96 1.20	(0.63–1.47) (0.64–2.26)	0.86 0.57	40/84 21/19	0.932.54	(0.62–1.39) (1.31–4.92)	0.72 0.01	0.96 1.51	(0.80–1.16) (1.11–2.05)	0.70 0.01	0.01	0.04
		rs876873	IVS6−530C>A	TT CT+CC	369/691 119/269	REFREF	32/69 20/34	0.90 1.33	(0.58–1.39) (0.73–2.41)	0.62 0.35	38/80 23/21	0.932.46	(0.61–1.39) (1.31–4.63)	0.71 0.01	0.95 1.52	(0.79–1.15) (1.14–2.04)	0.61 0.005	0.01	0.04
Glyphosate	ALOX5AP (arachidonate 5-lipoxygenase-activating protein)	rs9579645	IVS1+694A>C	AA AC+CC	151/241 31/91	REFREF	227/433 48/109	0.84 1.25	(0.64–1.09) (0.73–2.16)	0.18 0.42	200/410 71/120	0.781.74	(0.59–1.02) (1.03–2.93)	0.07 0.04	0.89 1.33	(0.77–1.02) (1.03–1.72)	0.08 0.03	0.004	0.02
Atrazine	LDLR (low-density lipoprotein receptor)	rs8110695	IVS1−4369T>A	TT AT+AA	128/223 60/150	REFREF	181/307 92/207	1.04 1.11	(0.78–1.40) (0.75–1.66)	0.79 0.59	154/315 117/199	0.861.49	(0.64–1.16) (1.01–2.19)	0.33 0.04	0.92 1.23	(0.79–1.07) (1.02–1.49)	0.27 0.03	0.02	0.03
	LRP1 (low-density lipoprotein receptor-related protein 1)	rs1800159	IVS62+86A>G	GG AG+AA	91/173 98/202	REFREF	117/261 156/254	0.85 1.32	(0.60–1.19) (0.95–1.82)	0.34 0.10	108/273 165/241	0.741.46	(0.52–1.05) (1.06–2.02)	0.090.02	0.86 1.20	(0.73–1.02) (1.02–1.40)	0.09 0.02	0.01	0.04

1Accounts for the number of SNPs per gene.

2Risk estimates adjusted for age and state.

We observed significant interactions with three insecticides, terbufos, fonofos, and carbofuran. Overall, the most noteworthy association (i.e., smallest FDR P-interaction and pesticide P-trend) was for terbufos and ALOXE3 rs3027208 (Table 3). Men carrying the variant T allele who were low users of terbufos had an OR of 1.86 (95%  CI = 1.16–2.99) and high users an OR of 2.00 (95%  CI = 1.28–3.15) compared to no use of terbufos (P trend = 0.001), while men carrying the CC genotype did not exhibit a significant association (FDR P interaction = 0.01). This association persisted after adjusting for ALOXE3 rs6503086, which had a significant main effect association with prostate cancer (FDR P trend = 0.02) but was not correlated with rs3027208 (r2 = 0.02). Terbufos also interacted significantly with SEC14L2 rs2072159 (FDR P interaction = 0.01) and ADIPOR1 rs12733285 (FDR P-interaction = 0.02). Fonofos interacted significantly with ABCG8 rs4953028, such that men carrying the GG genotype who were low users of fonofos had an OR of 1.73 (95% CI = 0.99–3.00) and high users an OR of 1.94 (95% CI = 1.17–3.20) compared to no use of fonofos (P trend = 0.004), while men carrying the variant A allele did not exhibit a significant association (FDR P interaction = 0.02). Carbofuran interacted significantly with A4GALT rs8136914 (FDR P interaction = 0.03) and with LDLR rs8110695 (FDR P interaction = 0.04).

We observed significant interactions with four herbicides, EPTC, petroleum oil/distillates, glyphosate, and atrazine. Among these, the most noteworthy association (i.e., smallest FDR P interaction and pesticide P trend) was for EPTC and ALOX15 rs916055. Men carrying the G allele who were low users of EPTC had an OR of 1.47 (95% CI = 0.98–2.22) and high users an OR of 1.63 (95% CI = 1.06–2.50) compared to no use of EPTC (P trend = 0.01), while men carrying the A allele did not exhibit a significant association (FDR P interaction = 0.01). Petroleum oil interacted with variants in three genes (ALOX5, LRP1, and ACSL5). The most significant interaction for petroleum oil was with ALOX5 rs7099684, such that men carrying the TT genotype who were low users of petroleum oil had an OR of 1.19 (95% CI = 0.76–1.86) and high users an OR of 1.56 (95% CI = 1.03–2.38) compared to no use of petroleum oil (P trend = 0.03), while men carrying the A allele did not exhibit a significant association (FDR P-interaction = 0.02). Petroleum oil interacted with four LRP1 SNPs (rs1800159, rs6581128, rs795957, rs7978567) that are moderately to highly correlated (r2 = 0.56–0.88). The most significant association (i.e., smallest FDR P-interaction and P-trend) was for LRP1 rs1800159, such that men carrying the A allele who were low users of petroleum oil had an OR of 1.26 (95% CI = 0.80–2.00) and high users an OR of 1.76 (95% CI = 1.12–2.77) compared to no use of petroleum oil (P trend = 0.01), while men carrying the G allele did not exhibit a significant association (FDR P interaction = 0.03). Petroleum oil also interacted significantly with two ACSL5 SNPs (rs876873, rs120985760) that are highly correlated (r2 = 0.89). Other noteworthy associations included the interactions between glyphosate and ALOX5AP rs9579645 (FDR P interaction = 0.02), and between atrazine and LDLR rs8110695 (FDR P interaction = 0.03) as well as LRP1 rs1800159 (FDR P-interaction = 0.04), which, as mentioned above, also interacted with petroleum oil.

4. Discussion

In this nested case-control study, we evaluated the interaction between pesticide use and SNPs in genes related to lipid metabolism. We found 17 interactions that were robust to multiple comparison adjustment and displayed a significant monotonic increase in prostate cancer risk with increasing pesticide use in one genotype group and no significant association in the other genotype group. With these criteria, our analyses yielded seven pesticides (atrazine, carbofuran, EPTC, fonofos, glyphosate, petroleum oil/distillate, and terbufos) and 15 SNPs in 11 genes (ABCG8, ACSL5, ADIPOR1, ALOX5, ALOX5AP, ALOX15, ALOXE3, A4GALT, LDLR, LRP1, SEC14L2) of primary interest.

Increased prostate cancer risk has been associated with the insecticides (terbufos, fonofos, carbofuran) that we found to interact with lipid metabolism genes in previous studies in the AHS. Terbufos was suggestively linked to prostate cancer [32], but significantly with aggressive prostate cancer [4], and showed evidence of interaction with variants in the 8q24 region [9] and xenobiotic metabolism genes [10]. Fonofos has been linked to prostate cancer [3, 6] and aggressive prostate cancer [4], and interacted significantly with variants in the 8q24 region, and the nucleotide and base excision repair genes [11, 12]. Similarly, carbofuran has been linked to prostate cancer [3], and interacted with variants in the nucleotide excision repair pathway [12]. In contrast to the insecticides, none of the four herbicides we identified have been associated with prostate cancer risk in main effect studies in the AHS (petroleum oil [3, 4], EPTC [3, 4, 33], atrazine [34, 35], glyphosate [36]). However, in previous analyses in the prostate nested case-control study, petroleum oil interacted with xenobiotic metabolism genes [10].

Our findings support the hypothesis that genetic variation in lipid metabolism genes might modify the associations between pesticides and prostate cancer; however, the biologic mechanisms are unclear. Organochlorines are lipophilic and accumulate in animal and human adipose tissue [37, 38]; however, none of the six organochlorines we evaluated interacted significantly with the SNPs we examined. This may be due to organochlorines being of the market for approximately 30 years and having a very small number of participants who were exposed to these pesticides. Some pesticides in other classes have lipophilic properties or are metabolized via lipid-related mechanisms [37, 38]. For example, organosphosphates have been shown to inhibit triglycerides and lipoprotein lipases in laboratory mice [22, 39, 40]. Two (terbufos and fonofos) of the nine organophophates we examined interacted significantly with lipid metabolism gene variants and were associated with an increased risk of prostate cancer. Terbufos interacted with a variant in ALOXE3 (arachidonate lipoxygenase 3), which is involved in a lipoxygenase reaction of a fatty acid that is converted to a molecule involved in the proliferation and specialization of the lipid membrane of the epidermis [41]. Terbufos also interacted with ADIPOR1 (adiponectin receptor 1) that encodes the adiponectin receptor and mediates ligand activities, fatty acid oxidation, and glucose uptake by adiponectin [42]; and SEC14L2 (SEC14-like protein 2) that is involved in lipid-binding and cholesterol biosynthesis [43]. Fonofos interacted with a variant in ABCG8 (ATP-binding cassette, sub-family G member 8), which encodes a protein that is a member of the ATP-binding cassette transporters that promotes biliary excretion of sterols [44]. Carbofuran belongs to another class of insecticide, carbamates, that inhibits cholinesterase and has been shown to generate reactive oxygen species and induce lipid peroxidation [45]. Carbofuran interacted with A4GALT (alpha 1,4-galactosyltransferase) which is involved in the formation of a glycolipid that provides cellular energy [46].

Among the herbicides, EPTC, classified as a thiocarbamate, has been reported to disrupt the biosynthesis of lipid metabolism in plants [47]. EPTC interacted with ALOX15, another gene in the lipoxygenase family that is involved in the metabolism of arachidonic and linoleic acid. Petroleum oil, which contains a mixture of hydrocarbons and is also lipophilic, is unique in that it is used as a herbicide, but also as an additive to other pesticides, and so it may have a wide variability in use [48]. It interacted with variants in three genes: ALOX5 (arachidonate 5-lipoxygenase), which transforms essential fatty acids into leukotrienes [49]; LRP1 (low-density lipoprotein receptor related protein 1), which regulates a large endocytic receptor that functions in lipoprotein transport [50]; ACSL5 (acyl-CoA synthetase long-chain 5) which converts free long-chain fatty acid to fatty acyl-CoA esters, and thereby plays a key role in lipid biosynthesis and fatty acid degradation [51]. Atrazine, a triazine herbicide, reported to be an endocrine disruptor in animals, was also shown to affect the lipid profiles of liver and muscle in fish [52]. It interacted with two genes belonging to the low-density lipoprotein receptor family (LDLR and LRP1). The low-density lipoprotein receptor (LDLR) is involved in the regulation of cholesterol homeostasis by receptor-mediated endocytosis [53].

While there is biologic plausibility to the interactions we observed, the specific biological mechanisms underlying these interactions are unclear. None of the SNPs we identified in our analysis are known to be functional or have been identified in published GWAS findings for prostate cancer (http://www.genome.gov/gwastudies/) [54]; however, they may be in linkage disequilibrium with other functional variants. Research conducted on the role of lipid metabolism in prostate carcinogenesis suggests that the mechanisms are complex, involving an array of intermediate compounds that can be modified by other pathways as well as various genetic, anthropometric, dietary, and lifestyle factors [55]. It is also possible that mechanisms other than lipid metabolism might play a role in the interactions we observed since most of the genes we identified have multiple functions. For example, genes within the arachidonic acid family play an important role in immunity and inflammation [56]. Furthermore, considering the findings from previous analyses from the nested case-control study of prostate cancer within AHS, we are seeing some of the same pesticides, such as terbufos, fonofos, carbofuran, EPTC, and petroleum oil, interacting with SNPs in multiple genetic pathways, suggesting that the relationship between pesticides and prostate cancer may involve multiple biologic processes.

In this analysis, we examined SNPs with a minor allele frequency (MAF) >10% because of limited power below this MAF threshold and restricted analyses to SNPs with a main effect association with prostate cancer at P < 0.2 to focus on SNPs potentially more associated with prostate cancer. However, by doing this, we may have excluded important SNPs that modify risk. Interaction analyses were conducted using the dominant genetic model since the number of participants in the homozygous variant group was often small; however, this may have resulted in loss of statistical power if another genetic model was more appropriate. Although some of our findings may be false-positives, we adopted a two-stage interaction approach and utilized FDR to limit this possibility. In addition, we focused our interpretations on interactions with positive monotonic associations between pesticide use and prostate cancer in one genotype and no significant association in the other genotype, which have greater biological plausibility than qualitative interactions, but could have resulted in exclusions of potentially relevant interactions. Due to small sample size we were unable to evaluate interactions by prostate cancer stage or grade.

In conclusion, we observed several positive interactions between pesticide use and variants in lipid metabolism genes. These findings suggest that lipid metabolism may be involved in the associations between pesticides and prostate cancer; however, little is known about the biologic mechanisms underlying these interactions. Furthermore, since this is the first study to explore the interaction between specific pesticide use and variants of genes in the lipid metabolism pathway, our results require replication.

Supplemental Table 1: shows the distribution of the study population by selected characteristics along with the chi-square p-value.

Supplemental Table 2: shows the association between intensity-weighted lifetime days of pesticide use and prostate cancer.

Supplemental Table 3: shows prostate cancer risk in relation to the 220 variants in 59 genes included in our interaction analysis. These variants had a Hardy Weinberg Equilibrium p-value>0.001, a minor allele frequency >=0.1, and a p-trend for the additive effect of each allele<0.2.

Supplemental Table 4: shows the associations between intensity-weighted lifetime days of pesticide use and prostate cancer risk stratified by SNP genotype. These associations showed significant FDR interactions between pesticide use and SNPs, and increased prostate cancer risk in at least one genotype

Click here for additional data file.