Skip to main content

Effect modification of air pollution on Urinary 8-Hydroxy-2'-Deoxyguanosine by genotypes: an application of the multiple testing procedure to identify significant SNP interactions



Air pollution is associated with adverse human health, but mechanisms through which pollution exerts effects remain to be clarified. One suggested pathway is that pollution causes oxidative stress. If so, oxidative stress-related genotypes may modify the oxidative response defenses to pollution exposure.


We explored the potential pathway by examining whether an array of oxidative stress-related genes (twenty single nucleotide polymorphisms, SNPs in nine genes) modified associations of pollutants (organic carbon (OC), ozone and sulfate) with urinary 8-hydroxy-2-deoxygunosine (8-OHdG), a biomarker of oxidative stress among the 320 aging men. We used a Multiple Testing Procedure in R modified by our team to identify the significance of the candidate genes adjusting for a priori covariates.


We found that glutathione S-tranferase P1 (GSTP1, rs1799811), M1 and catalase (rs2284367) and group-specific component (GC, rs2282679, rs1155563) significantly or marginally significantly modified effects of OC and/or sulfate with larger effects among those carrying the wild type of GSTP1 , catalase, non-wild type of GC and the non-null of GSTM1.


Polymorphisms of oxidative stress-related genes modified effects of OC and/or sulfate on 8-OHdG, suggesting that effects of OC or sulfate on 8-OHdG and other endpoints may be through the oxidative stress pathway.

Peer Review reports


Many studies have shown that ambient pollution is consistently associated with adverse health outcomes [16], but mechanisms accountable for these associations have not been fully elucidated. Suggested biological mechanisms linking air pollution and cardiovascular diseases include direct effect on the myocardium, disturbance of the cardiac autonomic nervous system, pulmonary and systematic oxidative stress and inflammatory response that triggers endothelial dysfunction, atherosclerosis and coagulation/thrombosis [7]. Understanding relative roles of such potential is a priority of recent air pollution epidemiology.

Some studies have demonstrated that exposures to particulate matter (aerodynamic diameter ≤2.5 μm, PM2.5) and ozone are associated with global oxidative stress [711]. Others reported that the exposures were associated with heart rate variability (HRV), plasma homocysteine and C-reactive protein and such effects were modified by genetic polymorphisms related to oxidative defenses [1216]. In living cells, reactive oxygen species (ROS) are continuously generated as a consequence of metabolic reactions, which may cause oxidative damage to nucleic acids. DNA damage may be repaired by the base excision repair pathway. The resulting repair product, 8-Hydroxy-2'-deoxyguanosine (8-OHdG), is the most common DNA lesion [17] and is not affected directly by either diet or cell turnover [18]. Therefore, 8-OHdG is a good biomarker for ROS or oxidative stress.

A limited number of epidemiological studies reported that 8-OHdG was associated with exposures to indoor and ambient pollution or smoking, but they were conducted among a small number of children or occupationally exposed employees [9, 10, 19]. Oxidative stress caused by air pollution may be implicated in the development of respiratory disease, cardiovascular disease, lung cancer and other diseases [2022]. Our recent study found that the elevated urinary 8-OHdG was associated with pollutants often thought of as secondary or formed through photochemical reactions after emission (PM2.5, nitrogen dioxide, NO2, maximal one-hour ozone, O3, sulfate, SO4 2-or organic carbon, OC), but not with directly emitted primary pollutants (black carbon, BC, carbon monoxide, CO or elemental carbon, EC), suggesting that secondary pollution plays a stronger role in oxidative stress [23].

Several studies have demonstrated that certain genetic polymorphisms related to oxidative stress modified effects of PM on cardiovascular responses [6, 13, 14], but a set of examined single nucleotide polymorphisms (SNPs) was very limited. Further, these studies only indirectly implicated oxidative stress as none of these outcomes was a direct measure of oxidative stress. For example, some studies reported that associations between exposure to PM2.5 and heart rate variability (HRV) were modified by polymorphisms of the glutathione-S-transferase M1 (GSTM1) gene [14] or heme oxygenase-1 (HMOX) [15], enzymes that reduce impacts of ROS. Our previous studies examined a set of genotypes related to oxidative stress and found that polymorphisms of hemochromatosis (HFE) and glutathione S-transferase T1 (GSTT1) significantly modified associations of PM2.5 with plasma homocysteine [12]. Anh et al. [24] reported that vitamin D-related genes (group-specific component, GC) were significantly associated with the serum D-vitamin concentrations that related to prostate cancer.

However, the selection of certain genes is somewhat arbitrary and the use of an array of genes is vulnerable to false positives from multiple comparisons, a major issue in genetic association studies. In this study, we aimed to examine whether daily ambient OC, SO4 2- and maximal one-hour O3 were associated with urinary 8-OHdG based on our previous findings [23] and such associations were modified by genotypes related to oxidative stress in the Normative Aging Study population (NAS). Because of multiple comparisons, we used the Multiple Testing Procedures (MTP) modified by our team, multtest in the R project to identify significant SNPs from a set of candidate genes [2528].


Study population

Data were obtained from a longitudinal NAS [29]. Briefly, the NAS is a longitudinal aging population initiated by the Veterans Administration (VA) in 1963. A total of 2,280 men from the greater Boston area free of known chronic medical conditions were enrolled. Subjects were asked to return for examinations every three to five years in the study center, including routine physical examinations, laboratory tests, collection of medical history, social status information, and administration of questionnaires on smoking history, food intake and other factors that may influence health. All participants provided written informed consents and the study protocol was approved by the institutions. By 2006, only did a small proportion of participants remain in the cohort, as many participants had died or were lost to follow up. A total of 320 participants, who still remained in this cohort, were included in our analyses, visiting the clinic between January 2006 and December 2008 for measurement of urinary 8-OHdG and other covariates (no repeated measurements).

8-hydroxy-2'-deoxyguanosine and plasma analysis of B vitamins

Urinary 8-OHdG analysis was conducted by Genox Corp (Baltimore, MD). A competitive enzyme-linked immunosorbent assay was used to analyze urinary 8-OHdG [30, 31]. The measurement methods have been described elsewhere [23]. Folate, vitamin B6 and B12 in fasting plasma were analyzed at the USDA Human Nutrition Research Center on Aging at Tufts University. Folate and vitamin B12 were examined by radioassay using a commercially available kit from Bio-Rad (Hercules, CA); vitamin B6 (as pyridoxal-5-phosphate) by an enzymatic method using tyrosine decarboxylase. Further details are described elsewhere [32, 33]. Plasma creatinine was measured with urine 8-OHdG using spectrophotometric assay. The method has been described elsewhere in details [34].

Air pollution and Weather Data

Averages of daily OC, SO4 2- and maximal one-hour O3 were used in this study. O3 and OC were provided by the Massachusetts Department of Environmental Protection and SO4 2- was measured at Harvard School Public Health monitoring station. For each day, SO4 2-, OC and O3 values were averaged for periods for up to four weeks before the visit as these averaging periods were shown to be most relevant in our previous analyses. Findings from our previous study showed that 8-OHdG were only associated with the secondary pollutants [23]. To adjust for weather condition, we used apparent temperature as an index, defined as a person's perceived air temperature, given the humidity [35].


In order to avoid the arbitrary selection of genes, we selected all 20 oxidative stress-related SNPs available in the NAS database. We examined effect modification using the array of candidate SNPs, including catalase (CAT, rs480575, rs1001179, rs2284367 and rs2300181), HFE H63 D (rs1799945), HFE C282Y (rs1800562), GSTM1, GSTT1, GSTP1 I105V (rs1695), GSTP1 A114V (rs1799811), HMOX (rs2071746, rs2071747, rs2071749, rs5995098), HMOX-1 VNTR, GC (rs2282679, rs1155563), glutamate cysteine ligase catalytic subunit (GCLC, rs17883901) and glutamate cysteine ligase modifier (GCLM, rs2301022 and rs3170633). HFE is related to cellular uptake of metals that are related to ROS generation and inflammation [8, 36]. Glutathione pathways play a vital role in cellular defenses against ROS [14, 3739]. Similarly, GC, GCLC and GCLM are related to glutathione-related metabolism [40, 41]. CAT helps catalyze hydrogen peroxide, a powerful ROS into water and molecular oxygen to maintain oxidative balance [39, 42]. HMOX-1 was categorized into two levels (any short and both long) based on repeated number of microsatellite (GTn) because previous studies have shown that a high GT repeats at 5'-flanking region may reduce HMOX-1 inducibility by ROS and has been associated with increased risk of cardiovascular diseases [15, 43, 44]. Previous studies have shown that variations of HFE C282Y, HFE H63 D, HMOX-1, GSTs genes modify associations of PM2.5 or BC with HRV or homocysteine [1215].

Multiplex polymerase chain reaction assays were designed using Sequenom SpectroDESIGNER software (Sequenom Inc, San Diego, Calif) by inputting sequence containing the SNP site and 100 bp of flanking sequence on either side of the SNP. Assays were genotyped using the Sequenom MassArray MALDI-TOF mass spectrometer (Sequonom, CA, USA) with semiautomated primer design (SpectroDESIGNER, Sequenom) and implementation of the very short extension method [45]. Assays failing to multiplex were genotyped using the TaqMan 5' exonuclease [Applied Biosystems (ABI), Foster City, CA, USA] with primers from ABI using radioactive labeled probes detected with ABI PRISM 7900 Sequence Detector System [46].

Statistical analyses

Statistical analyses were performed with R version 2.9.1. First, we fitted linear regression models to separately examine the association of a single pollutant with urinary 8-OHdG at different day moving averages up to four weeks during the study period to decide which day moving averages for each pollutant were strongly associated with 8-OHdG for effect modification assessment. We used the log-transformation of 8-OHdG to minimize residuals and to stabilize the variance. We identified a priori the following variables as important potential confounders based on our previous NAS studies and other studies [9, 12, 14]: age, body mass index (BMI), alcohol consumption (≥2 drinks/day; yes/no), smoking status (never, former, current), pack-years of cigarettes smoked, plasma folate, vitamin B6, B12, use of statin medication (yes/no) and season and chronic disease status (cardiovascular disease, diabetes and chronic cough). We controlled plasma folate, vitamin B6, B12, age, BMI and pack-years of cigarettes smoked as continuous variables and adjusted for alcohol consumption, smoking status, use of statin medication and season as categorical variables. We adjusted for temperature using three-day moving average of apparent temperature with linear and quadratic terms due to the potential nonlinear relationship between temperature and 8-OHdG. In addition, we adjusted for creatinine clearance rate using the Cockcroft-Gault formula ([140 - age(year)]*weight(kg)]/[72* serum creatinine(mg/dL)]) [47]. We also adjusted for chronic disease status (cardiovascular disease or chronic respiratory diseases) as a dummy variable [23].

We examined effect modification by each of candidate SNP via adding an interaction term of the SNP and the pollutant simultaneously with both the main effect terms adjusting for the same covariates as the above [12, 23]. Because two dozens of candidate SNPs were involved in the analyses, results were vulnerable to the multiple comparison problem. To decrease type I errors, we used MTP model to identify the significance of interaction terms of individual SNP and pollutant. The current version of MTP allows one to identify the significance of a group of candidate variables to reduce the false discovery rate meanwhile adjusting for a group of fixed covariates. We used MTP to identify the significance of the group of interaction terms. Because the current version of MTP in R can only include one term that varied across models, our team modified it to include two terms, i.e., the main effect term of genes and the interaction term of one pollutant and genes.

We used the family-wise error rate (fwer) for type I error adjustment, step-down max T (sd.maxT) for method and default values for others in MTP. We briefly described the rationale here. More details about the rationale are described elsewhere [2527]. MTP is based on Bootstrap estimation of the null distribution samples and the data generating distribution P. Samples are drawn at random with replacement from the observed data. The program generates B bootstrap samples from hypotheses M and obtains M × B samples or M × B matrix of test statistics. Then, based on the M × B matrix of test statistics, the bootstrap estimates or test statistics are induced. There are several methods to define type I error and calculate adjusted p-values in MTP. We selected family-wise error rate and step-down maxT methods in this study. In step-down procedures, the hypotheses corresponding to the most significant test statistics are considered successively, with further tests depending on the outcomes of earlier ones. Therefore, it is more powerful than a single-step. The adjusted p-values for the step-down maxT procedures are given by [26]

p ˜ r j = max h = 1 , j { Pr ( max l { r k , , r m } | T l | | t r k | | H 0 C ) }

where Pr refers to p-value, H denotes hypothesis, and T means test statistic.

MTP directly reported adjusted p-values. An advantage of this method as opposed to only rejection or not of hypotheses, is that it is not needed to determine the level of the test in advance. This study reported adjusted p-values. Then, we quantitatively estimated associations between the pollutants and 8-OHdG across those carrying variants of the significant genes identified by MTP with significant interactions using the bootstrap method with the combination of coefficients of the main effect and the interaction [6].


Table 1 shows the descriptive statistics of the demographic characteristics, health and environmental variables among the NAS population during 2006-2008 at visit (n = 320). There were no repeated measurements in this study. Table 2 shows distributions of polymorphisms of candidate genes. Among 320 participants, wild types were dominant for CATs, HFEs, GSTP1 (rs1799811), HMOX (rs2071749) and GCLC, but the situation varied for other candidate genes. There were no obvious differences for the distributions of wild and heterozygous types in GCLM, GC and GSTP1 (rs1695). Heterozygous types for HMOX (rs2071746 and rs2071749) consisted of large components. 80.9% and 48.8% of subjects were classified as non-deletions for GSTT1 and GSTM1, respectively. Mean of the HMOX-1 GC repeated number was 25.8 (SD: 3.3) with median 24.

Table 1 Descriptive statistics of the demographic characteristics, health and environmental variables among the male Normative Study Aging population at their visits during 2006-2008 at visit (n = 320)
Table 2 Genotype distribution of participants (N = 320)*

We first fit the linear regression model to estimate associations of OC, SO4 2- and maximal one-hour O3 with 8-OHdG using moving averages of pollutants up to four weeks. Results show that main effects varied across different day moving averages and 24-, 20- and 18-day moving averages were strongest associated with SO4 2-, OC and maximal one-hour O3, respectively, which were used to assess effect modifications. The detailed information has been reported elsewhere [23]. For an IQR increases in 24-, 20- and 18-day moving averages of daily SO4 2-, OC and maximal one-hour O3, urinary 8-OHdG increased by 29.0% (95% CI: 5.9%, 52.1%), 27.6% (95% CI: 3.6%, 51.6%) and 54.3% (95% CI: 7.6%, 100.9%), respectively.

Before examining effect modification, we categorized each candidate gene into a dummy variable so that the gene and the pollutant of interest only have one interaction term. We combined the homozygous and heterozygous types for appropriate genes known as the non-wild type (dominant model) due to small number of the homozygous type. We also combined the homozygous and heterozygous short repeat for HMOX-1, referred to as any short (Table 2). Then, we identified candidate genes that executed significant effect modification as aforementioned. Adjusted p-values in MTP model show that GSTP1 A114V (rs1799811) marginally significantly modified the effect of SO4 2- on 8-OHdG (adjusted p = 0.091). CAT (rs2286367) (adjusted p = 0.037), GSTM1 (adjusted p = 0.037), GC (rs2282679) (adjusted p = 0.025) and GC (rs1155563) (adjusted p = 0.027) significantly modified effects of OC on 8-OHdG. There was no significant effect modification for O3 (Table 3). As sensitive analyses, we used different options in MTP for typeone (type I error) (tail probabilities for error rate, TPPER; false discovery rate, FDR) and methods (single-step maximum T, ss.maxT; single-step minimum P ss.minP; step-down minimum P, ss.minP). Similar trends were found in spite of some variations. We also categorized pack-years of cigarettes smoked using tertiles as cut-off and re-ran MTP model. Results were similar to those using continuous variable for pack-years of cigarettes smoked. Figure 1 shows the estimated effects of OC or SO4 2- on 8-OHdG across subpopulations carrying different genotypes, for those SNPs where an interaction with p < 0.10 was found.

Table 3 Statistical p-values for the interaction between pollutants and SNPs from MTP model using family-wise error rate and step-down max T method *
Figure 1
figure 1

Estimated percent changes in 8-OHdG (log) (95% confident interval) associated with a unit increase of 17- and 20-day moving averages of organic carbon and sulfate, respectively by gene polymorphisms. Adjusting for apparent temperature, age, body mass index, smoking status, pack-years of cigarettes smoked, alcohol consumption, use of statin medication, plasma folate, vitamin B6 and B12, season, chronic disease and creatinine clearance rate. Wild^: non-wild; Delet: deletion, delet^: non-deletion.


We found that associations of the secondary pollutants, specifically OC and SO4 2-, with 8-OHdG, a direct oxidative stress-related biomarker, were modified by polymorphisms in genes related to oxidative defenses. This is significant for several reasons. First, the finding that genetic polymorphisms in the oxidative defense pathway modified the association suggests that it is not due to chance or confounding, since neither should be associated with the genotypes of the individuals. Second, while considerable focus has been placed recently on freshly generated traffic particles, such as BC or ultrafine particle number, this study confirms that particles, including particles from coal burning power plants, play a role in increasing systemic oxidative stress.

The specific polymorphisms that modified the associations were GSTP1 (rs1799811), GSTM1, CAT (rs1799811) and GC (rs22826799, rs1155563). We found 8-OHdG was more strongly associated with SO4 2- among those carrying the wild type of the GSPT1, and more strongly associated with OC among those carrying the wild type of CAT (rs2284367), the non-deletion of GSTM1 and the non-wild type of the GCs (rs2282679 and rs1155563) comparing with other types of the corresponding genes (Figure 1). Based on our knowledge, it is the first time that MTP has been used to identify significant gene-environment interactions. MTP has advantages over some other approaches to controlling for false discovery rates in which a group of fixed covariates are adjusted for while a set of variables were compared.

Several studies have examined effect modification and found that people carrying variants of oxidative stress-related genes are differentially susceptible to air [1214, 16, 48]. Human GSTs are subdivided into several classes, among which GSTT1, GSTM1 and GSTP1 have been extensively investigated [12, 14, 49, 50]. GSTM1 or GSTT1 catalyzes the conjugation of glutathione to numerous potentially genotoxic compounds [50]. Individuals with the deletion of GSTM1 or GSTT1 have been shown to reduce GST activity and thus may be unable to eliminate toxins as efficiently when they expose to oxidative pollutants [50]. Schwartz et al. [14] found that PM2.5 was significantly associated with high frequency of HRV among those without the GSTM1 allele, but not for those with the allele. Gilliland et al. [48] reported that exposure to in utero maternal smoking was associated with increased prevalence of early onset asthma among those without GSTM1 allele, but not for those with GTSM1 allele. Similarly, Romieu et al. [51] found that GSTM1 null children were more sensitive to ozone exposure. However, all the aforementioned studies did not report whether there were significant effect modifications. Differential results from these stratification analyses might also be attributed to statistical powers across subpopulations or differential distributions of other controlled or uncontrolled covariates across subpopulations. This study observed that GSTM1 significantly modified associations of OC with 8-OHdG, but paradoxically that the GSTM1 null allele provided protection against exposure. Our recent study examined whether variations of a set of genes altered effects of black carbon and PM2.5 on plasma homocysteine in this population and found that GSTT1 (but not GSTM1) significantly modified associations between pollutants and homocysteine. PM2.5 and black carbon were more strongly associated with homocysteine among those carrying GSTM1 allele comparing those without the allele although no significant interactive effects were found [12]. Different findings of effect modification by GSTM1 variation across studies may reflect differences of exposure, outcome and population, measurement errors in exposure or phenotype, and by chance. Similar situations also appeared in other studies [52, 53]. Therefore, statistical effect modification may be inconsistent with biological interaction. Further research or meta-analysis is needed for GSTM1.

In contrast, few studies have examined the function of GSTP1 A114V (rs1799811) on diseases with inconsistent results [5457]. None of these studies found the GSTP1 is significantly associated with the outcomes of interest although some studies found positive trends. Therefore, the functions of the polymorphisms have not been determined. Several studies examined effect modifications of GSTT1 on various endpoints but no significant effect modification was found [5860]. For example, Melén et al. [59] examined whether GST modified traffic-related pollution effect on childhood allergic disease and found that carriers with variants of GSTP1 (rs1799811) were higher susceptible to NOx. Our study found the variation of GSTP1 showed a protective effect of SO4 2- on 8-OHdG. However, other two studies did not find any evidence that the GSTP1 modified effects of black carbon or smoking on blood pressure or Parkinson's disease occurrence [58, 60]. Inconsistent observed findings may be attributable to various sources as aforementioned. In this study, it may also related to the small number of variants in this population, which probably lead to unstable estimates. Therefore, its functions remain to be clarified by others (Table 2).

GC, vitamin D-related genes, is related to the vitamin D metabolism [61]. Vitamin D is activated to form 1, 25-dihydroxyvitamin D in the liver and kidney and then transported in serum to different tissues by the vitamin D-binding protein, which is encoded by GC [61]. Studies show that polymorphisms of vitamin D-related genes are associated with various cancers, cardiovascular diseases and respiratory diseases [6264]. Ahn et al. [61] examined variations of 212 SNPs related to vitamin D metabolism and found that all four SNPs of GC (rs1212631, rs2282679, rs7041, rs1155563) are significantly associated with the concentration of serum vitamin D. When these four SNPs were simultaneously included in the multivariate model, only two SNPs (rs22679, rs1155563) were significantly associated with vitamin D. In this study, we found that the two SNPs of GC (rs22679, rs1155563) were associated with 8-OHdG in this study. The mechanisms remain to be clarified yet.

Catalase is a protein of 526 amino acids, encoded by the catalase gene with 34 kb pairs of nuclear acids [65]. Catalase is the main regulator of hydrogen peroxide metabolism [66]. Catalase enzyme mutations may reduce its activity and probably results in the increase of the hydrogen peroxide concentrations in the tissues [62]. Inherited catalase deficiency results in acatalasemia (homozygous state) and hypocatalasemia (heterozygous) and is related to increased plasma homocysteine concentrations [42, 67, 68]. Our previous study reported that the variation of CAT modified associations between particle matter and plasma homocysteine concentrations [12].

Experimental toxicology studies have shown that air pollutants act via the oxidative stress pathway [8, 36, 69]. Ghio et al. [36] found that homozygous Belgrade rats functionally deficient in divalent metal transporter-1 display decreased metal transport from the lower respiratory tract and have stronger lung injury than control littermates, when exposed to oil fly ash containing iron. Belgrade rats cannot transport iron and other divalent metals across membranes via HFE gene regulated processes. They also reported that healthy volunteers exposed to concentrated ambient air particles had increased concentrations of blood fibrinogen and induced mild pulmonary inflammation [8]. Tamagawa et al. [69] reported that five-day and four-week exposures to PM10 caused acute and chronic lung and systematic inflammation of New Zealand rabbits.

There are several strengths in this study. First, we used MTP model to identify the significance of a group of candidate genes while we examined effect modification by genes on air pollution effects. This method overcame some problems in this kind of studies, such as arbitrary selection of a few significant genes or high false discovery rate when individually examining a set of genes. Secondly, this study was conducted in a relatively large population. Information of participants was well measured and collected. However, several limitations also exist with this study. First, we used air pollution data collected from a single monitoring site for personal pollution exposure and therefore, some extent misclassification might happen. A recent study compared ambient concentrations with personal exposures with monitoring measurement and results show that ambient measures were good surrogates for PM2.5 and SO4 2- in both winter and summer, but O3 was only good in summer, not well in winter [70]. Nevertheless, with non-differential misclassification, any potential bias would be expected toward the null. Second, MTP has several options to select type I error and several methods to calculate adjusted p-values. Using bootstrap re-sampling methods will result in different estimates when a MTP model is rerun. These will introduce the uncertainties in model selections [2528]. In addition, the NAS consists of an aged population and non-Hispanic white men were dominant. Thus, the findings are not well generalizable to other populations.


This study found that variations of oxidative stress-related genes modified effects of OC or SO4 2- on 8-OHdG. This suggests that effects of OC or SO4 2- on 8-OHdG and other endpoints may be through the oxidative stress pathway.



black carbon


organic carbon


element of carbon


single nucleotide polymorphism


nitrogen dioxide


carbon monoxide





PM2.5 :

particulate matter ≤2.5 μm in aerodynamic diameter


glutathione S-tranferase




group-specific component




heme oxygenase-1


glutamate cysteine ligase catalytic subunit


glutamate cysteine ligase modifier


  1. Schwartz J: The effects of particulate air pollution on daily deaths: a multi-city case-crossover analysis. Occup Environ Med. 2004, 61: 956-961. 10.1136/oem.2003.008250.

    Article  CAS  Google Scholar 

  2. Bell ML, McDermott A, Zeger SL, Samet JM, Dominici F: Ozone and short-term mortality in 95 US urban communities, 1987-2000. JAMA. 2004, 292: 2372-2378. 10.1001/jama.292.19.2372.

    Article  CAS  Google Scholar 

  3. Dominici F, Peng RD, Bell ML, Pham L, McDermott A, Zeger SL, Samet JM: Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA. 2006, 295: 1127-1134. 10.1001/jama.295.10.1127.

    Article  CAS  Google Scholar 

  4. Zanobetti A, Schwartz J: Particulate air pollution, progression, survival after myocardial infarction. Environ Health Perspect. 2007, 115: 769-775. 10.1289/ehp.9201.

    Article  CAS  Google Scholar 

  5. Ren C, Willims GM, Morawska L, Mengersen K, Tong S: Ozone modifies associations between temperature and cardiovascular mortality: analysis the NMMAPS data. Occup Environ Med. 2008, 65: 255-260. 10.1136/oem.2007.033878.

    Article  CAS  Google Scholar 

  6. Ren C, Baccarelli A, Wilker E, Suh H, Sparrow D, Vokonas P, Wright R, Schwartz J: Lipid and endothelial related genes, ambient particulate matter, and heart rate variability --the VA Normative Aging Study. J Epidemiol Community Health. 2010, 64: 49-56. 10.1136/jech.2008.083295.

    Article  CAS  Google Scholar 

  7. Brook RD: Cardiovascular effects of air pollution. Clin Sci. 2008, 115: 175-187. 10.1042/CS20070444.

    Article  CAS  Google Scholar 

  8. Ghio AJ, Kim C, Devlin RB: Concentrated ambient air particles induce mild pulmonary inflammation in healthy human volunteers. Am J Respir Crit Care Med. 2000, 162: 981-988.

    Article  CAS  Google Scholar 

  9. Kim JY, Mukherjee S, Ngo L, Christiani DC: Urinary 8-hydroxy-2'-deoxyguanosine as a biomarker of oxidative DNA damage in workers exposure to fine particles. Environ Health Perspect. 2004, 112: 666-671. 10.1289/ehp.6827.

    Article  CAS  Google Scholar 

  10. Gurgueira SA, Lawrence J, Coull B, Murthy GK, González-Flecha B: Rapid increases in the steady-state concentration of reactive oxygen species in the lungs and heart after particulate air pollution inhalation. Environ Health Perspect. 2002, 110: 749-755. 10.1289/ehp.02110749.

    Article  CAS  Google Scholar 

  11. Vinzents PS, Møller P, Sørensen M, Knudsen LE, Hertel O, Jensen FP, Schibye B, Loft S: Personal exposure to ultrafine particles and oxidative DNA damage. Environ Health Perspect. 2005, 113: 1485-1490. 10.1289/ehp.7562.

    Article  CAS  Google Scholar 

  12. Ren C, Park SK, Vokonas PS, Sparrow D, Wilker E, Baccarelli A, Suh H, Schwartz J: Air pollution and homocysteine: more evidence that oxidative stress-related genes modify effects of particulate air pollution. Epidemiology. 2010, 21: 198-206. 10.1097/EDE.0b013e3181cc8bfc.

    Article  Google Scholar 

  13. Park SK, O'Neill MS, Wright RO, Hu H, Vokonas PS, Sparrow D, Suh H, Schwartz J: HFE genotype, particulate air pollution, and heart rate variability--a gene-environment interaction. Circulation. 2006, 114: 2798-2805. 10.1161/CIRCULATIONAHA.106.643197.

    Article  CAS  Google Scholar 

  14. Schwartz J, Park SK, O'Neill MS, Vokonas PS, Sparrow D, Welss S, Kelsey K: Glutathione-S-transferase M1, obesity, statins, and autonomic effects of particles: gene-by-drug-by-environment interaction. Am J Respir Crit Care Med. 2005, 172: 1529-1533. 10.1164/rccm.200412-1698OC.

    Article  Google Scholar 

  15. Chahine T, Baccarelli A, Litonjua A, Write RO, Suh H, Gold DR, Sparrow D, Vokonas P, Schwartz J: Particulate air pollution, oxidative stress genes, and heart rate variability in an elderly cohort. Environ Health Perspect. 2007, 115: 1617-1622. 10.1289/ehp.10318.

    Article  CAS  Google Scholar 

  16. Zeka A, Sullivan JR, Vokonas PS, Sparrow D, Schwartz J: Inflammatory markers and particulate air pollution: characterizing the pathway to disease. Int J Epidemiol. 2006, 35: 1347-1354. 10.1093/ije/dyl132.

    Article  Google Scholar 

  17. Kasai H, Crain PF, Kuchino Y, Nishimura S, Ootsuyama A, Tanooka H: Formation of 8-hydroxygunine moiety in cellular DNA by agents producing oxygen radicals and evidence for its repair. Carcinogenesis. 1986, 7: 1849-1851. 10.1093/carcin/7.11.1849.

    Article  CAS  Google Scholar 

  18. Cooke MS, Evans MD, Dove R, Rozalski R, Gackowski D, Siomek A, Lunec J, Olinski R: DNA repair is responsible for the presence of oxidative damaged DNA lesions in urine. Mutat Res. 2005, 574 (1-2): 58-66.

    Article  CAS  Google Scholar 

  19. Lu CY, Ma YC, Lin JM, Chuang CY, Sung FC: Oxidative DNA damage estimated by urinary 8-hydroxydeoxyguanosine and indoor air pollution among non-smoking office employees. Environ Res. 2007, 103: 331-337. 10.1016/j.envres.2006.08.009.

    Article  CAS  Google Scholar 

  20. Chuang KJ, Chang CC, Su TC, Lee CT, Tang CS: The effect of urban air pollution on inflammation, oxidative stress, coagulation, and autonomic dysfunction in young adults. Am J Respir Crit Care Med. 2007, 176: 370-376. 10.1164/rccm.200611-1627OC.

    Article  CAS  Google Scholar 

  21. Wiseman H, Halliwell B: Damage to DNA by reactive oxygen and nitrogen species: role in inflammatory disease and progression to cancer. Biochem J. 1996, 313: 17-29.

    Article  CAS  Google Scholar 

  22. Higashi Y, Noma K, Yoshizumi M, Kihara Y: Endothelial function and oxidative stress in cardiovascular diseases. Circ J. 2009, 73: 411-418. 10.1253/circj.CJ-08-1102.

    Article  CAS  Google Scholar 

  23. Ren C, Fang S, Wright RO, Suh H, Schwartz J: Urinary 8-hydroxy-2'-deoxyguanosine as a biomarker of oxidative DNA damage induced by ambient pollution in the Normative Aging Study. Occup Environ Med. 2010, [Online Oct 27, 2010]

    Google Scholar 

  24. Anh J, Albanes D, Berndt SI, Peters U, Chatterjee N, Freedman ND, Abnet CC, Huang WY, Kibel AS, Crawford DE, Weinstein SJ, Chanock SJ, Schatzki A, Hayes RB: Vitamin D-related genes, serum vitamin D concentrations and prostate cancer risk. Carcinogenesis. 2009, 30 (5): 769-776. 10.1093/carcin/bgp055.

    Article  Google Scholar 

  25. Pollard KS, Dudiot S, van der Lann MJ: Multiple testing procedures: R multtest package and application to genetics. 2005, []

    Google Scholar 

  26. Dudoit S, Shaffer JP, Boldrick JC: Multiple hypothesis testing in microarray experiments. Stat Sci. 2003, 18: 71-103. 10.1214/ss/1056397487.

    Article  Google Scholar 

  27. Dudoit S, van der Laan M, Pollard KS: Multiple testing part I: single-step procedures for control of general type I error rates. Stat Appl Genet Mol Biol. 2004, 3: 13-

    Google Scholar 

  28. van der Laan M, Dudoit S, Pollard KS: Multiple testing part II: step down procedures for control of family-wise error rate. Stat Appl Genet Mol Biol. 2004, 3: 14-

    Google Scholar 

  29. Bell B, Rose C, Damon A: The veterans Administration longitudinal study of healthy aging. Gerontologist. 1966, 6: 179-184.

    Article  CAS  Google Scholar 

  30. Erhola M, Toyokuni S, Okada K, Tanaka T, Hiai H, Ochi H, Uchida K, Osawa T, Nieminen MM, Alho H, Kellokumpu-Lehtinen Pl: Biomarker evidence of DNA oxidation in lung cancer patients: association of urinary 8-hydroxy-2'-deoxyguanosine excretion with radiotherapy, chemotherapy, and response to treatment. FEBS lett. 1997, 409: 287-291. 10.1016/S0014-5793(97)00523-1.

    Article  CAS  Google Scholar 

  31. Leinonen J, Lehtimaki T, Toyokuni S, Okada K, Tanaka T, Hiai H, Ochi H, Laippala P, Rantalaiho V, Virta O, Pasternack A, Alho H: New biomarker evidence of oxidative DNA damage in patients with non-insulin-dependent diabetes mellitus. FEBS Lett. 1997, 417: 150-152. 10.1016/S0014-5793(97)01273-8.

    Article  CAS  Google Scholar 

  32. Park SK, O'Neill MS, Vokonas PS, Sparrow D, Spiro A III, Tucker KL, Suh H, Hu H, Schwartz J: Traffic-related particles are associated with elevated homocysteine - the VA Normative Aging Study. Am J Respir Crit Care Med. 2008, 178: 283-289. 10.1164/rccm.200708-1286OC.

    Article  Google Scholar 

  33. Tucker KL, Qiao N, Scott T, Rosenberg I, Spiro A: High homocysteine and low B vitamins predict cognitive decline in aging men: the Veterans Affairs Normative Aging Study. Am J Clin Nutr. 2005, 82: 627-635.

    CAS  Google Scholar 

  34. Bowers L, Wong E: Kinetic serum creatinine assays. II. A critcal evaluation and review. Clin Chem. 1980, 26: 555-

    CAS  Google Scholar 

  35. Kalkstein L, Valamont K: An evaluation of summer discomfort in the United States using a relative climatologic index. Bull Am Meteorol Soc. 1986, 67: 842-848. 10.1175/1520-0477(1986)067<0842:AEOSDI>2.0.CO;2.

    Article  Google Scholar 

  36. Ghio AJ, Piantadosi CA, Wang X, et al.: Divalent metal transporter-1 decreases metal-related injury in the lung. Am J Physiol Lung Cell Mol Physiol. 2005, 289: 460-467. 10.1152/ajplung.00154.2005.

    Article  Google Scholar 

  37. Hayes JD, McLellan LI: Glutathione and glutathione dependent enzymes respresent a co-ordinately regulated defense against oxidative stress. Free Radic Res. 1999, 31: 273-300. 10.1080/10715769900300851.

    Article  CAS  Google Scholar 

  38. Gilliland FD, LI YF, Saxon A, Diaz-Sanchez D: Effect of glutathione-S-transferase M1 and P1 genotypes on xenobiotic enhancement of allergic responses: randomized, placebo-controlled crossover study. Lancet. 2004, 363: 119-125. 10.1016/S0140-6736(03)15262-2.

    Article  CAS  Google Scholar 

  39. Forsberg L, Lyrenäs L, de Faire U, Morgenstern R: A common functional C-T substitution polymorphisms in the promoter region of the human catalase gene influences transcription factor binding, reporter gene transcription and is correlated to blood catalase levels. Free Radic Biol Med. 2001, 30: 500-505. 10.1016/S0891-5849(00)00487-1.

    Article  CAS  Google Scholar 

  40. Engstöm KS, Strömberg U, Lundh T, Johansson I, Vessby B, Hallmans G, Skerfving S, Broberg K: Genetic variation in glutathione-related genes and body burden of methylmercury. Environ Health Perspect. 2008, 116: 734-739.

    Article  Google Scholar 

  41. Siedlinski M, Postma DS, van Diemen CC, Blokstra A, Smit HA, Boezen HM: Lung function loss, smoking, vitamin C intake, and polymorphisms of the glutamate-cysteine ligase genes. Am J Respir Crit Care Med. 2008, 178: 13-19. 10.1164/rccm.200711-1749OC.

    Article  CAS  Google Scholar 

  42. Góth L, Vitai M: The effects of hydrogen peroxide promoted by homocysteine and inherited catalase deficiency on human hypocatalasemic patients. Free Radic Biol Med. 2003, 35: 882-888.

    Article  Google Scholar 

  43. Chen YH, Lin SJ, Lin MW, Tsai HL, Kuo SS, chen JW, Charng MJ, Wu TC, Chen LC, Ding PYA, Pan WH, Jou YS, Chau LY: Microsatellite polymorphism in promoter of heme oxygenase-1 gene is associated with susceptibility to coronary artery disease in type 2 diabetes patients. Hum Genet. 2002, 111: 1-8. 10.1007/s00439-002-0769-4.

    Article  CAS  Google Scholar 

  44. Kaneda H, Ohno M, Taguchi J, Togo M, Hashimoto H, Ogasawara K, Aizawa T, Ishizaka N, Nagai R: Heme oxygenase-1 gene promoter polymorphism is associated with coronary artery disease in Japanese patients with coronary risk factors. Arterioscler Thromb Vasc Biol. 2002, 22: 1680-1685. 10.1161/01.ATV.0000033515.96747.6F.

    Article  CAS  Google Scholar 

  45. Sun X, Ding H, Hung K, Guo B: A new MALDI-TOF based mini-sequencing assay for genotyping of SNPs. Nucleic Acids Res. 2000, 28: e68-10.1093/nar/28.12.e68.

    Article  CAS  Google Scholar 

  46. Lee LG, Connell CR, Bloch W: Allelic discrimination by nick-translation PCR with fluorogenic probes. Nucleic Acids Res. 1993, 21: 3761-3766. 10.1093/nar/21.16.3761.

    Article  CAS  Google Scholar 

  47. Cockcroft DW, Gault MH: Prediction of creatinine clearance from serum creatinine. Nephron. 1976, 16: 31-41. 10.1159/000180580.

    Article  CAS  Google Scholar 

  48. Gilliland FD, Li Y, Dubeau L, Berhane K, Avol E, Gauderman WJ, Peters JM: Effects of glutathione-S-transferase M1, maternal smoking during pregnancy, and environmental tobacco smoke on asthma and wheezing in children. Am J Respir Crit Care Med. 2002, 457-463. 10.1164/rccm.2112064. 166

  49. Bergamaschi E, De Palma G, Mozzoni P, Vanni S, Vettori MV, Broeckaert F, Bernard A, Mutti A: Polymorphism of quinone-metabolizing enzymes and susceptibility to ozone-induced acute effects. Am J Respir Crit Care Med. 2001, 163: 1426-1431.

    Article  CAS  Google Scholar 

  50. Couphlin SS, Hall IJ: Glutahione S-transferase polymorphisms and risk of ovarian cancer: a HuGE review. Genet Med. 2002, 4: 250-257. 10.1097/00125817-200207000-00003.

    Article  Google Scholar 

  51. Romieu I, Sienra-Monge JJ, Ramírez-Aguilar M, Moreno-Macías H, Reyes-Ruiz NI, Estela del Rio-Navarro B, Hernández-Avila M, London SJ: Genetic polymorphism of GSTM1 and antioxidante supplementation influence lung function in relation to ozone exposure in asthmatic children in Mexico City. Thorax. 2004, 59: 8-10.

    CAS  Google Scholar 

  52. Gilliland FD, Rappaport EB, Berhane K, Islam T, Dubeau L, Gauderman WJ, McConnell R: Effects of glutathione S-transferase P1, M1, and T1 on acute respiratory illness in school children. Am J Respir Crit Care Med. 2002, 346-351. 10.1164/rccm.2111048. 166

  53. Gilliland FD, Gauderman WJ, Vora H, Rappaport E, Dubeau L: Effects of glutathione-S-transferase M1, T, and P1 on childhood lung function growth. Am J Respir Crit Care Med. 2002, 166: 710-716. 10.1164/rccm.2112065.

    Article  Google Scholar 

  54. Al-Dyyel F, Al-Rasheed M, Ibrahim M, Bu R, Bavi P, Abubaker J, Al-Jomah N, Mohamed GH, Moorji A, Uddin S, Siral AK, Al-Kuraya K: Polymorphisms of drug-metabolizing enzymes CYP1A1, GSTT and GSTP contributed to the development of diffuse large B-cell lymphoma risk the Saudi Arabian population. Leuk Lymphoma. 2008, 49: 122-129. 10.1080/10428190701704605.

    Article  Google Scholar 

  55. Gemignani F, Landi S, Szeszenia-Dabrowska N, Zaridze D, Lissowska J, Rudnai P, Fabianova E, Mates D, Foretova L, Janout V, Bencko V, Gaborieau V, Gioia-Patricola L, Bellini1 I, Barale R, Canzian F, Hall J, Boffetta P, Hung RJ, Brennan P: Development of lung cancer before the age of 50: the role of xenobiotic metabolizing genes. Carcinogenesis. 2007, 28: 1287-1293. 10.1093/carcin/bgm021.

    Article  CAS  Google Scholar 

  56. Yang XR, Pfeiffer PM, Goldstein AM: Influence of glutathione-S-transferase (GSTM1, GSTP1, GSTT1) and cytochrome p450 (CYP1A, CYP2D6) polymorphisms on numbers of basal cell carcinomas (BCCs) in families with the naevoid basal cell carcinoma syndrome. J Med Genet. 2006, 43: e16-10.1136/jmg.2005.035006.

    Article  Google Scholar 

  57. De Roos AJ, Gold LS, Wang S, Hartge P, Cerhan JR, Cozen W, Yeager M, Chanock S, Rothman N, Severson RK: Metabolic gene variants and risk of non-Hodgkin's lymphoma. Cancer Epidemiol Biomarkers Prev. 2006, 15: 1647-1653. 10.1158/1055-9965.EPI-06-0193.

    Article  CAS  Google Scholar 

  58. Wahner AD, Glatt CE, Bronstein JM, Ritz B: Glutathione S-transferase mu, omega, pi, and theta class variants and smoking in Parkinson's disease. Neurosci Lett. 2007, 413: 274-278. 10.1016/j.neulet.2006.11.053.

    Article  CAS  Google Scholar 

  59. Melén E, Nyberg F, Lindgren CM, Berglind N, Zucchelli M, Nording E, Hallberg J, Svartengren M, Morgenstern R, Kere J, Bellander T, Wickman M, Pershagen G: Interactions between glutathione S-transferase P1, tumor necrosis factor, and traffic-related air pollution for development of childhood allergic disease. Environ Health Perspect. 2008, 116: 1077-1084.

    Article  Google Scholar 

  60. Mordukhovich I, Wilker E, Suh H, Wright R, Sparrow D, Vokonas PS, Schwartz : Black carbon exposure, oxidative stress genes, and blood pressure in a repeated measures study. Environ Health Perspect. 2009, 117: 1767-1772.

    CAS  Google Scholar 

  61. Ahn J, Albanes D, Berndt SI, Peters U, Chatterjee N, Freedman ND, Abnet CC, Huang W, Kibel AS, Crawford ED, Weinstein SJ, Chanock SJ, Schatzikin A, Hayes RB: Vitamin D-related genes, serum vitamin D concentrations and prostate cancer risk. Carcinogenesis. 2009, 30: 769-776. 10.1093/carcin/bgp055.

    Article  CAS  Google Scholar 

  62. Raimondi S, Johansson H, Maisonneuve P, Gandini S: Review and meta-analysis on vitamin D receptor polymorphisms and cancer risk. Carcinogenesis. 2009, 30: 1170-1180. 10.1093/carcin/bgp103.

    Article  CAS  Google Scholar 

  63. McCullough ML, Bostick RM, Mayo TL: Vitamin D gene pathway polymorphisms and risk of colorectal, breast, and prostate cancer. Annu Rev Nutr. 2009, 29: 111-132. 10.1146/annurev-nutr-080508-141248.

    Article  CAS  Google Scholar 

  64. Wang TJ, Pencina MJ, Booth SL, Jacques PF, Ingelsson E, Lanier K, Benjamin EJ, D'Agostino RB, Wolf M, Vasan RS: Vitamin D deficiency and risk of cardiovascular disease. Circulation. 2008, 117: 503-511. 10.1161/CIRCULATIONAHA.107.706127.

    Article  CAS  Google Scholar 

  65. Quan F, Korneluk RG, Tropak MB, Gravel RA: Isolation and characterization of the human catalase gene. Nucleic Acids Res. 1986, 14: 5321-5335. 10.1093/nar/14.13.5321.

    Article  CAS  Google Scholar 

  66. Mueller S, Riedel HD, Stremmel W: Direct evidence for catalase as the predominant H2O2 removing enzyme in erythrocytes. Blood. 1997, 90: 4973-4978.

    CAS  Google Scholar 

  67. Ahn J, Nowell S, McCann SE, Yu J, Carter L, Lang NP, Kadlubar FF, Ratnasinghe LD, Ambrosone CB: Associations between catalase phenotype and genetype: modification by epidemiologic factors. Cancer Epidemiol Biomarkers Prev. 2006, 15: 1217-1222. 10.1158/1055-9965.EPI-06-0104.

    Article  CAS  Google Scholar 

  68. Góth L, Rass P, Páy A: Catalase enzyme mutations and their association with diseases. Mol Diagn. 2004, 8: 141-149.

    Article  Google Scholar 

  69. Tamagawa E, Bai N, Morimoto K, Yatera E, Zhang X, Xing L, Li Y, Laher I, Sin DD, Man SFP, van Eeden SF: Particulate matter exposure induces persistent lung inflammation and endothelial dysfunction. Am J Physiol Cell Mol Physiol. 2008, 295: 79-85. 10.1152/ajplung.00048.2007.

    Article  Google Scholar 

  70. Sarnat JA, Brown KW, Schwartz J, Coull BA, Koutrakis P: Ambient gas concentrations and personal particulate matter exposures: implications for studying the health effects of particles. Epidemiology. 2005, 16: 385-395. 10.1097/01.ede.0000155505.04775.33.

    Article  Google Scholar 

Download references


This work was supported by the National Institute of Environmental Health Sciences grants ES014663, ES 15172, and ES-00002, by U.S. Environmental Protection Agency grant R832416 and USDA Contract 58-1950-7-707. The Normative Aging Study is supported by the Cooperative Studies Program/Epidemiology Research and Information Center of the U.S. Department of Veterans Affairs, and is a component of the Massachusetts Veterans Epidemiology Research and Information Center. It is partially supported by Harvard-NIOSH ERC Pilot (T42 OH008416).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Cizao Ren.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CR was responsible for study design, data analyses, result interpretation and manuscript writing. JS was responsible for study design, data collection and result interpretation. Other coauthors participated in the study design, data collection and result interpretation. All authors read and approved the final manuscripts.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Ren, C., Vokonas, P.S., Suh, H. et al. Effect modification of air pollution on Urinary 8-Hydroxy-2'-Deoxyguanosine by genotypes: an application of the multiple testing procedure to identify significant SNP interactions. Environ Health 9, 78 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Heart Rate Variability
  • Multiple Test Procedure
  • Oxidative Stress Pathway
  • Normative Aging Study
  • Significant Effect Modification