Abstract
Chronic exposure to ozone (O3) can cause changes in lung function that may reflect remodelling of small airways. It is likely that antioxidant enzyme function affects susceptibility to O3. The aim of the present study was to determine whether polymorphisms in antioxidant enzyme (GSTM1, GSTP1 and NQO1) genes affect the risk of lung function changes related to chronic exposure to O3.
In total, 210 young adults who participated in a previous study, which showed a relationship between lifetime exposure to O3 and decreased lung function, were genotyped. Multivariable linear regression was used to model sex-specific associations between genotypes and O3-related lung function changes, adjusting for height, weight, lifetime exposure to nitrogen dioxide and particles with a 50% cut-off aerodynamic diameter of 10 μm, and self-identified race/ethnicity.
The GSTM1-null/NQO1 Pro187Pro-combination genotype was significantly associated with increased risk of an O3-related decrease in mean forced expiratory flow between 25–75% of forced vital capacity in females (parameter estimate±se -75±35 mL·s−1), while the GSTP1 Val105 variant genotypes were significantly associated with greater risk of an O3-related decrease in mean forced expiratory flow at 75% of forced vital capacity in males (-81±31 mL·s−1). GSTM1-null status was not significantly associated with any O3-related changes in lung function in either sex.
The current authors conclude that the effects of antioxidant enzyme gene polymorphisms on the risk of decreased lung function related to chronic exposure to ozone may be modified by sex-specific factors.
Ozone (O3), a major component of air pollution, is a potent oxidant gas that causes airway injury in human lungs 1. A large proportion of inhaled O3 (up to 90%) is absorbed in the respiratory tract along the entire tracheobronchial tree 2, with the greatest dose being delivered to the peripheral airways at the junction between the conducting and respiratory airways 3. O3 reacts with respiratory tract lining fluid constituents and cellular membrane components to generate lipid ozonation products (LOP) and reactive oxygen species (ROS), which in turn can cause oxidative damage to other biological molecules 4–6. Acute exposure to ambient levels of O3 can induce short-term lung function abnormalities and airway inflammation, while chronic exposure may lead to remodelling of the small airways, where deposition is the greatest 7–12.
To minimise the potential for oxidative injury, the human lung has an integrated system of antioxidant enzymes and expendable soluble molecules. This system includes several mechanisms by which ROS are converted to products that are further detoxified by other enzymes. If the oxidant burden is sufficiently great, ROS may overwhelm the antioxidant system leading to a state of “oxidative stress”, which is thought to contribute to the pathogenesis of a number of respiratory diseases 13–15. Although antioxidant defences are available to decrease oxidative stress in the airways, individuals differ in their ability to deal with an oxidant burden; such differences are, in part, determined genetically 16. This genetic variability may account for the considerable between-subject variability seen in both the lung function and airway inflammatory responses to O3 17, 18.
Glutathione S-transferase (GST) enzymes, a superfamily of dimeric phase-II metabolic enzymes, play an important role in the antioxidant defence system. GST enzymes catalyse the conjugation of toxic electrophilic molecules with glutathione and thereby protect cellular macromolecules from damage due to LOP and ROS. The specific GST enzymes that have been proposed as important in antioxidant defence are those of the mu (GSTM), theta (GSTT) and pi (GSTP) classes, each with functional polymorphisms that affect protein expression or function 19.
A common polymorphism in the GSTM1 gene locus, which exists in 30–50% of the general population 19, 20, involves a null allele and results in a complete lack of GSTM enzymatic function. Therefore, the GSTM1 null genotype would be expected to affect the individual's response to O3 exposure, possibly causing increased susceptibility to oxidative injury. Since 2001, the results of several field studies showed that the GSTM1 null genotype is associated with greater acute lung function response and/or respiratory symptom response to O3-induced oxidative stress 21–23. Other studies have suggested that the GSTM1 null genotype may play a significant role in the development of asthma in response to oxidative stress 24, 25. In two of the aforementioned studies, a polymorphism (Ser187) of a second antioxidant enzyme, NQO1 (reduced nicotinamide adenine dinucleotide (phosphate):quinone oxidoreductase 1), provided a protective effect among GSTM1 null subjects 21, 25.
Another gene of interest with regard to responses to oxidant pollutants is GSTP1, which is the most abundant GST in lung tissue and has a common A105G polymorphism that results in an Ile105Val amino acid substitution. In the Children's Health Study (CHS), children who were homozygotes for GSTP1 Val105 variant allele had a lower rate of respiratory infections than those with the GSTP1 Ile/Ile105 wild type but, somewhat surprisingly, they also had a slower rate of lung function growth 26, 27. A recent study showed that the GSTP1 Val/Val105 genotype was associated with increased O3-related respiratory symptoms 23. In contrast, in a small (n = 19) controlled exposure study of sensitised allergic rhinitic adults, in which nasal instillation of diesel exhaust particles (known to cause oxidative stress) enhanced specific allergic responses to ragweed, Gilliland et al. 28 showed that the GSTM1 null genotype increased susceptibility and the GSTP1 Val105 variant had a protective effect.
Considered together, these human studies provide suggestive evidence that polymorphisms of phase-II enzymes contribute to susceptibility to inhaled oxidant-induced toxicity. In a recent epidemiological study 11, the current authors’ group demonstrated an association between lifetime exposure to ambient O3 and decreased lung function parameters consistent with small airway remodelling. To determine whether the GSTM1 null, GSTM1 null/NQO1 homozygous Pro187 combination or GSTP1 Val105 variant genotypes had an effect on the observed relationship between lifetime exposure to O3 and decreased lung function, the subjects who participated in the previous study 11 were genotyped, and it was assessed whether these genotypes affected the risk of O3-induced lung function changes. The current authors selected these three genotypes for study on the basis of the previous literature reviewed above.
METHODS
The protocol for the study was approved by the Committee for the Protection of Human Subjects, University of California, Berkeley (UCB; CA, USA), and the Committee on Human Research, University of California, San Francisco, CA, USA. Written informed consent was obtained from all study participants once eligibility was established.
Study design
The overall design of the study has been previously presented in detail 11. Briefly, a convenience sample of 255 freshman undergraduates at UCB was recruited in three waves that began on April 10, 2000, February 12, 2001 and February 6, 2002. All waves ended in the first week of June. Subjects were studied between February–May when students from Los Angeles (LA; CA, USA) would not have been exposed to high summertime O3 concentrations.
Students were eligible based on the following criteria: 1) lifelong resident of the greater LA or San Francisco Bay (SF) area prior to enrolment at UCB; 2) lifetime never-smoker, 3) no history of chronic respiratory disease (history of asthma before age 12 yrs was permitted, provided that student had no symptoms and had not taken any medication at any time after age 12 yrs (n = 6)); and 4) no physical impairment that would hinder performance of spirometry. Location of all residences within the geographical boundaries for the study was confirmed by study personnel.
Ozone exposure assessment
A detailed description of the creation of lifetime cumulative O3 exposure for each subject has been previously reported 10. Briefly, lifetime residential history was reconstructed with a standardised questionnaire and air pollutant (O3, NO2 and particles with a 50% cut-off aerodynamic diameter of 10 μm (PM10)) concentrations were assigned for each month of life to each residential location. Air quality data were acquired from the California Air Resources Board (ARB; CD No. PTSD-02-017-CD), the Aerometric Information Retrieval System and by special requests to the ARB. Monthly mean measures of O3 were interpolated spatially from air quality monitoring stations to the residence locations with inverse distance weighting and a maximum of three monitoring stations for each interpolation (maximum interpolation radius of 50 km). The details and reliability of the exposure assignment method have been previously published 10, 29, 30. Briefly, two basic models are fitted in order to estimate lifetime pollutant (O3, NO2 and PM10) exposure. There was no significant difference in the association between lifetime O3 exposure and lung function between the two models. In the present study, the so-called “ecological” model was used, which omitted estimates of time spent outdoors and used only the residence-specific monthly average interpolated pollutant concentrations.
Subject characteristics
Of the 255 enrolled subjects, 226 had sufficient DNA available for genotyping. The 29 subjects who were not genotyped were mostly female. In total, 16 subjects did not self-identify with one of three main racial/ethnic groups (Asians/Pacific Islanders, Caucasians and Hispanics) and were excluded from further analysis. Thus, 210 subjects were used in the final analysis (table 1⇓). Of these, ∼43% were male and 60% were lifelong residents of the LA area. Most of the participants were Asians (54% of males and 60% of females) or Caucasians (39% of males and 28% of females). There were no significant sex differences in lifetime exposure estimates of pollutants (O3, PM10 and NO2). Although subjects who grew up in LA had higher median estimated lifetime exposures than those from the SF area, distributions between the two regions overlapped and represented a continuum of individual exposure.
Antioxidant enzyme genotyping
DNA was isolated from clot with a Qiamp Blood DNA Maxi kit (Qiagen Inc., Santa Clarita, CA, USA) in accordance with the manufacturer's instructions and stored at -80°C until use. Genotyping for the GSTM1 polymorphism was carried out following a previously reported protocol 31. The TaqMan real time PCR method was used to detect polymorphisms of GSTP1 (A105G) and NQO1 (C187T). Primers and probes for the single-nucleotide polymorphisms were custom-designed by Applied Biosystems Inc. (Foster City, CA, USA; see online data supplement for primer sequences). The reaction was carried out in TaqMan Universal Master Mix with a 7900 Real-Time PCR machine (Applied Biosystems). Quality assurance procedures included: assessment of randomly distributed blank samples; duplicates of randomly selected samples; manual calls assisting automated calling for Taqman analysis; and repeated additional analysis from independently isolated DNA samples from the same subjects. Assays were repeated for all low-confidence samples until a reliable call was obtained. The genotype frequencies for GSTM1, GSTP1, and GSTM1 null/NQO1 did not deviate from Hardy–Weinberg equilibrium.
Spirometry
Forced expiratory volumes were obtained in the sitting position with nose clip with a Collins Survey rolling seal spirometer (Warren E. Collins Co., Braintree, MA, USA) with two modifications to the standard criteria of the American Thoracic Society 32, details of which have been previously reported 33. Forced vital capacity (FVC), forced expiratory volume in one second (FEV1), mean forced expiratory flow between 25–75% of FVC (FEF25–75%) and mean forced expiratory flow at 75% of FVC (FEF75%) were recorded. The FEF25–75%/FVC ratio, an estimate of the reciprocal of the time constant of the lung 34 and a reflection of intrinsic airway size 33, was also calculated. This measure was used, in part, to control for racial/ethnic differences in airway size 11. Sex-specific models for each measure of lung function were fit based on height, weight and age as described previously 11, 31. There was no association between any measure of lung function and history of asthma before age 12 (n = 6) or history of second-hand tobacco smoke exposure (n = 34) 11.
Statistical analysis
Genotypes for antioxidant enzyme polymorphisms were coded as follows: GSTM1 wild type (positive) = 0 or null = 1 and GSTP1 Ile105 wild type (homozygous AA) = 0 or Val105 variant (heterozygous AG or homozygous GG) = 1. The GSTM1/NQO1 combination genotype included only GSTM1 null subjects and used the following coding: NQO1 wild type Pro187 (homozygous CC) = 0 or Ser187 variant (heterozygous CT or homozygous TT) = 1. Because it is well established that lung function differs between males and females 35, sex-specific multivariable linear regression was used to model lung function variables. Except for FEV1, natural logarithmic transformations of lung function variables were used. The initial model for each lung function variable included the subject's height, weight and race/ethnicity, and each genotype and was based on the optimal model among several tested, as previously described 11, 31. Final models included height, weight, race/ethnicity, genotype and lifetime exposures to O3, PM10 and NO2, as previously described 11, 31. The combined genotype, GSTM1·NQO1, was treated as single term, based on reported interactions between the null allele of GSTM1 and the serine polymorphism for NQO1. The FEF25–75%/FVC ratio was treated as an interaction term for reasons previously discussed 11. Finally, measurement error correction procedures were not used for the O3 effects, since the current authors have shown previously 11 that the coefficients for O3 and the O3–(FEF25–75%/FVC) interaction were not affected by such corrections.
RESULTS
Allele and genotype frequencies differed significantly among ethnicities (table 2⇓). For GSTM1, frequency of GSTM1 null status was significantly different for both Hispanics (60%) and Caucasians (57%) compared with Asians (44%). For NQO1, the NQO1 T allele (Ser187) was significantly lower in Caucasians (18%) than Asians (40%) and Hispanics (35%). For GSTP1, the G allele (Val105) was significantly more common in Hispanics (48%) than Caucasians (20%) and Asians (29%).
Effect of genotype on lung function
The effects of the three polymorphisms on lung function can be seen in table 3⇓ (see online data supplement for complete model parameter estimates). In males, GSTM1 null homozygosity was associated significantly with decreased FEF25–75% (-98 mL·s−1; 95% confidence interval (CI): -15– -181 mL·s−1). However, when the GSTM1 null genotype occurred in combination with NQO1 wild type (homozygous Pro187), this association was no longer significant. Similarly, GSTM1 null was associated significantly with decreased FEF75% (-133 mL·s−1; 95% CI: -27– -240 mL·s−1) while GSTM1 null/NQO1 wild-type combination genotype was not. For FEV1, no association between either the GSTM1 null or GSTM1 null/NQO1 wild-type combination genotype was found in males (data not shown).
In contrast to males, there were no significant changes in FEF25–75% and FEF75% for the GSTM1 null variant alone in females. However, a decrease in FEF25–75% (-136 mL·s−1; 95% CI: -29– -243 mL·s−1) and FEF75% (-125 mL·s−1; 95% CI: -2– -253 mL·s−1) was associated significantly with the GSTM1 null/NQO1 wild-type combination genotype. For FEV1, there was no association with either GSTM1 null or GSTM1 null/NQO1 wild type-combination genotype in females (data not shown).
Finally, effects of GSTP1 Val105 in male and females subjects also differed. The variant allele was marginally associated with changes in both FEF25–75% and FEF75% in both sexes. However, for males the variant allele was nonsignificantly associated with decreases in these flow measures, while for females the trend is for an increase. No associations were found between GSTP1 val105 and FEV1 in either males or females (data not shown).
Effect of genotype and lifetime O3 exposure on lung function
To explore the effect of genotype on risk of lung function changes due to chronic exposure to O3, the final models of the present study included lifetime exposure to O3, PM10 and NO2, and an interaction term for O3·(FEF25–75%/FVC), in addition to the variables included in the initial genotype-only models. When this approach was used in sex-specific models for GSTM1 null and GSTM1 null/NQO1 wild-type combination genotypes, the only significant association with risk of O3-related decreased lung function was observed in females for the combination genotype and FEF25–75% (table 4⇓).
The GSTP1 variant allele, however, was associated with greater risk of O3-related decreases in FEF25–75% (p<0.11) and FEF75% in males (p<0.04; table 4⇑) after adjustment for lifetime O3 exposure and its interaction with airway size. The magnitude of the effect can be estimated from the final model based on the male-specific 25th percentile FEF25–75%/FVC ratio and mean lifetime O3 exposure difference between subjects from LA and SF (17 ppb). For males who are homozygous GSTP1 Iso105 (wild type), the 17 ppb lifetime O3 exposure difference results in a 20 mL·s−1 (95% CI: -18– -22 mL·s−1) decrease in FEF75%. For males carrying the Val105 variant GSTP1 allele, the 17 ppb lifetime O3 exposure difference results in a 28 mL·s−1 (95% CI: -26– -30 mL·s−1) decrease in FEF75%. The magnitude of the combined effect of the GSTP1 variant allele and lifetime O3 exposure is almost 50% less for males with the median FEF25–75%/FVC ratio (i.e. larger airway size).
For females, the GSTP1 Val105 variant allele did not have a statistically significant effect on lifetime O3-related decreases in lung function.
DISCUSSION
In a previous study, the current authors showed that estimated lifetime exposure to ambient O3 in a cohort of adolescents was associated with reduced levels of lung function measures that reflect the function of small airways 11. It was found that, without consideration of the effect of O3, the male subjects of the cohort with the GSTM1 null genotype had lower lung function measures that reflect small airways function, compared with those without this genotype. The current authors did not find this same gene effect for female subjects of the cohort. However, when lifetime exposure to O3 was included in the models, no deleterious role for GSTM1 null was found on lung function in either sex, although the GSTM1 null/NQO1 wild-type combination genotype was associated with increased risk of O3-related decreases in FEF25–75% in females.
A novel finding of the present study is that the GSTP1 Val105 variant genotype was a risk factor for decreased lung function in association with lifetime exposure to O3 in males. Conversely, the data suggest that this genotype may have a protective effect in females. This sex difference in the effect of the GSTP1 Val105 variant genotype may help explain the finding of greater male sensitivity to O3-induced lung function changes, which was previously reported by the current authors 11. Although the possible mechanism for such sex-specific modification can only be speculated upon, GSTP1 is known to have sex-specific patterns of expression 36, 37.
Previously, the current authors have shown that the deleterious effect of O3 on lung function was dependent on intrinsic airway size (measured by the FEF25–75%/FVC ratio) 33, 34, with a more deleterious effect of O3 on lung function occurring in subjects with smaller airway size 11. In the present study, in a model that includes antioxidant enzyme genotypes as well as the FEF25–75%/FVC ratio, the results show that the deleterious effect of O3 on lung function remains dependent on airway size.
Gilliland et al. 27 have studied a large group of subjects from the CHS in Southern California and found that non-Hispanic white children with GSTM1 null genotype had a lower rate of lung function growth. The results from male subjects of the present cohort are consistent with those of the CHS study. However, the analysis presented by the CHS investigators was not stratified by sex. Therefore, the role of sex in modifying the effects of enzyme genotypes on growth of lung function warrants further investigation.
Several studies have suggested an association between O3-induced airway oxidative injury and certain antioxidant enzyme genetic polymorphisms in nonasthmatic subjects, specifically with the GSTM1 null alone and the GSTM1 null/NQO1 wild-type combination genotypes. In a small field study, Bergamaschi et al. 21 showed an association between the O3 level in ambient air and decrements in lung function and changes in plasma CC16 in individuals with the GSTM1 null/NQO1 wild type-combination genotype only. Later, in a controlled exposure study, the same group of investigators showed a differential change in some biomarkers of oxidative stress after O3 exposure between subjects with the GSTM1 null/NQO1 wild-type combination genotype and those with other genotypes 38. The results of the current authors’ chronic exposure study also suggest that females with the combined GSTM1 null/NQO1 wild-type combination genotype have increased susceptibility to O3-related remodelling of the small airways.
The lack of concordant findings with regard to male subjects may be due to multiple differences between the Italian studies 21, 38 and the present study. Bergamaschi et al. 21 and Corradi et al. 38 studied the effects of acute O3 exposure while the current authors studied the effect of chronic lifetime exposure. The sample size of the present study is also much larger, potentially allowing sex-specific differences to be uncovered in the gene–environment interaction. In addition, the racial/ethnic composition of the study populations is most likely quite different, since the present population included Asian and Hispanic subjects, who were probably not represented in the Italian studies.
The genetic background of subjects from different self-identified racial/ethnic groups, which includes their genotypes for other antioxidant enzymes, probably plays an important role in determining their responses to O3 exposure. Population stratification, which can cause spurious associations in candidate gene-association studies, exists when the total population has been formed by admixture of two or more ancestral populations and when admixture proportions vary among individuals. If the risk of the outcome varies with ancestry proportions, then admixture will confound associations of the outcome with genotypes at any locus where allele frequencies vary between ancestral populations. Because genotype frequencies for the three candidate genes varied across racial/ethnic groups (table 2⇑), race/ethnicity was adjusted for in the regression models. However, inclusion of race/ethnicity in the models had little impact on the results (data not shown). Additionally, exclusion of the 16 subjects that did not self-identify as Asian, Caucasian or Hispanic did not significantly change the results of the analyses.
The current authors acknowledge several limitations of the study. First, while larger than many of the other studies that have assessed the effects of the GSTM1 null, GSTM1 null/NQO1 wild-type combination and GSTP1 Val105 genotypes on response to oxidant pollutants, the present study population is too small to definitely assess gene–environment interactions, especially if sex-specific modification is present. Secondly, other genes that were not studied are likely to play a role in determining susceptibility to chronic exposure to O3. Finally, although the current authors attempted to control for population stratification by including self-identified race/ethnicity in regression models, the use of genetic markers might have improved the ability to do so.
In conclusion, it was found that the GSTP1 Val105 variant genotype increases the risk of deleterious effects of chronic exposure to ozone on measures of lung function that reflect small airway remodelling in a group of healthy adolescent males. However, it was also found that this genotype may have a protective effect in their female counterparts. Unlike previous reports from smaller studies of acute exposures, the current authors did not find the GSTM1 null genotype in either sex or the GSTM1 null/NQO1 wild-type combination genotype in males to be associated with decreased lung function due to chronic exposure to ozone. However, it was found that the GSTM1 null/NQO1 wild-type combination genotype increases the risk of ozone-related loss of mean forced expiratory flow between 25–75% of forced vital capacity in females. The results of the present study suggest that the effects of antioxidant enzyme gene polymorphisms on the risk of decreased lung function related to chronic exposure to ozone may be modified by sex-specific factors.
Acknowledgments
The authors greatly appreciate the expert assistance of K. Beckman (Children’s Hospital of Oakland, Oakland, CA, USA) and M. Bastaki (University of California, Berkeley, CA, USA) in the genotyping, the dedication of research assistants L. Carlton, J. Murphy and S. Deamer (all University of California), and the participation of the University of California students who made the study possible.
Footnotes
-
This article has supplementary material accessible from www.erj.ersjournals.com
- Received December 8, 2006.
- Accepted June 29, 2007.
- © ERS Journals Ltd