Published ahead of print on May 11, 2006, doi:10.1164/rccm.200508-1281OC
© 2006 American Thoracic Society doi: 10.1164/rccm.200508-1281OC
Variants in the Glutamate-Cysteine-Ligase Gene Are Associated with Cystic Fibrosis Lung DiseaseDivision of Pulmonary and Critical Care Medicine, and Department of Environmental and Occupational Health Sciences and Center for Ecogenetics and Environmental Health, University of Washington, Seattle, Washington; and UBC McDonald Research Laboratories/iCAPTURE Center, Vancouver, British Columbia, Canada Correspondence and requests for reprints should be addressed to Terrance J. Kavanagh, Ph.D., Department of Environmental and Occupational Health Sciences, University of Washington Mail Box 354695, Seattle, WA 98195. E-mail: tjkav{at}u.washington.edu
Background: Chronic progressive lung disease is the most serious complication of cystic fibrosis (CF). Glutathione plays an important role in the protection of the CF lung against oxidant-induced lung injury. Objectives: We hypothesized that a polymorphism in a novel candidate gene that regulates glutathione synthesis might influence CF lung disease. Methods: In a cross-sectional study, subjects were recruited from CF clinics in Seattle and multiple centers in Canada. We tested for an association between CF lung disease and a functional polymorphism in the glutamate-cysteine ligase catalytic subunit (GCLC) gene. Multiple linear regression was used to test for association between polymorphisms of GCLC and severity of CF lung disease while adjusting for age, Pseudomonas aeruginosa infection, and cystic fibrosis transmembrane conductance regulator (CFTR) genotype. Analysis was repeated for patients with CF stratified by CFTR genotype. Measurements and Main Results: A total of 440 subjects with CF participated in the study (51% male; mean [± SD] age, 26 ± 11 yr; mean FEV1, 62 ± 28% predicted). In the total population, there was a trend toward an association between GCLC genotypes and CF lung disease (linear regression coefficient [SEM], 1.68 [1.0]; p = 0.097). In the stratified analysis, there was a highly significant association between GCLC genotype and CF lung function in subjects with a milder CFTR genotype (linear regression coefficient [SEM], 5.5 (1.7); p = 0.001). Conclusions: In patients with CF with a milder CFTR genotype, there is a strong association between functional polymorphisms of the GCLC gene and CF lung disease severity.
Key Words: CFTR genotype glutathione modifier genes Cystic fibrosis (CF) is an autosomal recessive genetic condition that is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene (1). CF is characterized by the clinical phenotype of recurrent respiratory infection, pancreatic insufficiency, and elevated sweat chloride. Median survival is into the mid-thirties, with over 90% of patients dying of progressive respiratory failure (2). The clinical course in CF is highly variable and is not fully explained by mutations in the CFTR gene (35). Family-based studies have suggested that host genetic factors that code outside of CFTR may modify the CF phenotype (6, 7). Identifying modifier genes that influence the development of lung disease is of primary importance in understanding the pathophysiology of CF lung disease and in identifying new therapeutic strategies to slow the progression of disease (8, 9). Glutathione (GSH) is a ubiquitous tripeptide that plays a crucial role in the protection of the lung from oxidant-induced lung damage (10). Oxidant-induced lung injury is an important mechanism of lung damage in CF (1114). Patients with CF have low levels of GSH in the epithelial lining fluid (15), suggesting that deficiency in GSH may influence the development of CF lung disease. We hypothesized that functional polymorphisms in a novel candidate gene that regulates GSH synthesis may modify the severity of lung disease in patients with CF.
The rate-limiting step in the production of GSH is the formation of We tested whether GCLC genotypes are associated with more severe CF lung disease. In addition, because certain CFTR genotypes are associated with a milder clinical phenotype (3, 5, 2125) and GSH secretion into the epithelial lining fluid may be influenced by CFTR function (2629), we hypothesized that this effect would be influenced by CFTR genotype. Some of the results presented in the manuscript have been presented in abstract form (30).
In a pilot study performed between 1999 and 2001, patients were recruited from the CF Clinic at the University of Washington, Seattle, WA. Based on interesting early results (30), the study population was subsequently expanded to include patients with CF recruited from Canadian CF Clinics in Vancouver, Victoria, Montreal, and Hamilton. All patients were diagnosed with CF based on established criteria (31). The study design was a cross-sectional study testing for an association between GCLC GAG TNR genotypes and severe CF lung disease. At time of entry into the study, the following clinical variables were collected on each subject using chart review: FEV1(% predicted), age, sex, presence or absence of Pseudomonas aeruginosa infection, pancreatic insufficiency status, and CFTR genotype. The best FEV1(% predicted) in the year before recruitment was used to define the pulmonary phenotype. P. aeruginosa infection was defined as present if the patient had at least one positive sputum sample in the year before recruitment. Pancreatic insufficiency was defined as present based on the need for pancreatic enzyme supplementation. CFTR genotypes were classified as "severe" and "mild" based on the effects of their CFTR mutations on CFTR production and clinical phenotype (5). Using the functional class system proposed by Tsui (32) and by Welsh and Smith (33), patients with both mutant alleles in Class I, Class II, or Class III were grouped together as a severe CFTR genotype, and patients with at least one mutant allele in Class IV, V, and Unclassified were considered a mild CFTR genotype. This classification system has been shown to be associated with significant differences in lung disease severity, with CFTR genotypes from Class IIII being associated with worse lung function than those with a Class IV, V, or Unclassified genotype (5). Microsatellite analysis of the GCLC TNR polymorphism is presented in the online supplement.
Statistical Analysis
Clinical characteristics of the total patient population and the Seattle and Canadian cohorts are shown in Table 1. Clinical and genetic data were available on a total of 440 patients with CF. The mean age of patients was 26.2 yr (range, 562 yr). Mean FEV1(% predicted) was 62% (range, 11128%). GCLC GAG TNR genotype frequencies are shown in Table 2 and were similar to that observed in the white population (20). GCLC GAG TNR genotype frequencies were in Hardy-Weinberg equilibrium.
Age, age squared (age2), the presence of P. aeruginosa infection, and CFTR genotype (severe versus mild) were associated with FEV1(% predicted) on univariate analysis (p < 0.2). On multivariate analysis, there was a trend toward an association between GCLC TNR polymorphism and FEV1(% predicted) after adjusting for age, age2, P. aeruginosa infection, and CFTR genotype (linear regression coefficient [SEM], 1.68 [1.0]; p = 0.097). In the subgroup with a mild CFTR genotype, there was a highly statistically significant association between GCLC TNR polymorphisms and FEV1(% predicted) (linear regression coefficient [SEM], 5.5 [1.7]; p = 0.001), explaining 4.3% of the variation in FEV1(% predicted). There was no association between GCLC TNR polymorphisms and CF lung disease in patients with a severe CFTR genotype (Table 3). Graphical interpretation of the multiple linear regression models demonstrating the effect of GCLC genotype on FEV1(% predicted) is shown in Figure 1.
Internal Validation of Results The decision to expand recruitment was based on pilot data found in the Seattle cohort. (30) This pilot study found a significant association between GCLC polymorphisms and CF lung disease in patients with a milder CFTR mutation (linear regression coefficient [SEM], 13.4 [2.8]; p < 0.001). To determine the replication validity of our pooled findings, we performed post-hoc analysis of the Canadian patients with CF with a mild CFTR genotype separately. The association was not as strong as that observed in the Seattle cohort but was in the same direction and approached statistical significance (linear regression coefficient [SEM], 3.93 [2.14]; p = 0.069).
The main findings of this study are that functional polymorphisms of the GCLC gene are associated with CF lung disease and that this effect is influenced by CFTR genotype. These observations are consistent with the known effect of GCLC polymorphisms on gene product and the biological role of GSH and CFTR in CF lung disease. In addition, the strength of association is large and highly statistically significant. We believe that this supplies compelling evidence that the GCLC gene is a genetic modifier of CF lung disease.
Patients with CF demonstrate substantial phenotypic variability (4, 5). Since the discovery of the CFTR gene that causes CF, over 1,000 mutations have been identified (34), yet only a small number of CFTR genotypes are associated with a milder clinical phenotype (3, 5, 2125). The effect of CFTR genotype on phenotype seems to relate to the functional effects of CFTR genotypes on CFTR production with very low levels of CFTR associated with a severe phenotype and intermediate levels of CFTR associated with milder CF (25, 35, 36). Despite this, there is evidence that the clinical course in patients with the same CFTR genotype is highly variable, and family-based studies have suggested that genetic factors that code outside of CFTR may modify the CF phenotype (6, 7). A number of possible CF modifier genes have been identified (3744), although a recent large, multicentered study found that, with the exception of variants in the TGF- Genes that influence the GSH pathway have been examined as potential modifiers of the CF phenotype, with polymorphisms in the glutathione-S-transferases M1, M3, and P1 being associated with a more severe CF phenotype (37, 43, 46). This is the first study to describe an association between a novel candidate gene that regulates GSH production and CF lung disease severity. These findings are supported by considerable evidence that GSH plays an important role in the pathophysiology of CF lung disease. Lung inflammation in patients with CF is characterized by an increase in oxidant production from activated neutrophils and macrophages on the respiratory epithelial surface (4749). This increase in oxidant production is thought to be one of the principal mechanisms through which the lung epithelium becomes damaged in CF (1114). In the healthy lung, one of the first lines of defense against oxidant stress is the extracellular antioxidant GSH, which is abundantly present in the epithelial lining fluid (ELF). In the CF lung, ELF levels of reduced and oxidized GSH are low, suggesting that in patients with CF, a relative deficiency of GSH is present (15). The reason for this seems to be related to the direct or indirect effects of CFTR on GSH transport out of epithelial cells. CFTR may be an important regulator of GSH release into the ELF (29), with deficiency in CFTR associated with lower levels of GSH (26, 27) and restoration of chloride channel function leading to increased GSH secretion (28). Our observation that GCLC TNR GAG polymorphisms have a significant influence on CF lung function only in patients with a mild CFTR genotype is consistent with the molecular biology of how CFTR influences the CF phenotype and GSH transport into the ELF. Patients with a severe CFTR genotype have very low levels (< 3%) of functioning CFTR protein, which leads to a more severe CF phenotype. Patients with a mild CFTR genotype have a milder clinical phenotype that is most likely due to the expression of increased amounts (513%) of functioning CFTR protein (25, 35, 36). Because CFTR influences the transport of GSH into the ELF, patients with CF with GCLC variants associated with greater GSH production, in the setting of more functioning CFTR, should have higher levels of GSH in the ELF and enhanced protection against oxidant-induced lung injury. Conversely, the severe reduction in functional CFTR observed in patients with CF with a severe CFTR genotype would result in low levels of GSH in the ELF irrespective of the effect of GCLC GAG TNR polymorphisms on GSH synthesis, which is consistent with our findings. Although the results of this study require confirmation at a molecular level, these findings may have important therapeutic implications. Our findings show that lung disease severity in patients with CF is influenced by factors that regulate lung GSH levels. It is therefore plausible that therapeutic interventions that increase lung GSH may reduce the decline in lung function seen in CF. Clinical trials that examine the role of supplemental GSH (or GSH precursors) in CF should consider stratification of enrolled patients by CFTR and GCLC genotype because the therapeutic effects may vary across different genetic subgroups. Our data also suggest that systemic therapies that increase epithelial cell GSH production could be of benefit to patients with CF with a mild CFTR genotype. For patients with CF with a severe CFTR genotype, our results suggest that systemic therapies that increase epithelial cell GSH levels may not be as beneficial, although inhaled therapy may be a more appropriate method of increasing GSH levels in the ELF of these patients (50). There are a number of limitations to our study. First, false-positive genetic association studies may occur as a result of chance, population stratification, or the presence of linkage disequilibrium between the allele of interest and the true disease-modifying allele. It is unlikely that chance or population stratification would fully explain our findings, given the strength of the observed association and the replication of findings in different centers (51). We are currently examining for linkage disequilibrium between GCLC TNR GAG repeats and other variants in the region of the GCLC gene, although, because the nearest gene (NM_021814) is > 140 kb downstream of the GCLC gene (52), it is likely that the GCLC gene is the true modifier gene. Second, we assumed a linear doseresponse relationship between GCLC GAG TNR genotype and CF lung function. This is based on earlier data that suggest that the GCLC 7 GAG TNR polymorphisms result in less GSH compared with the 8 GAG TNR and 9 GAG TNR (20). This has yet to be confirmed in patients with CF. Ongoing studies in our laboratory are examining GCLC mRNA and protein levels to fully understand the differences in transcriptional efficiency, message stability, or translational efficiency associated with GCLC GAG TNR alleles. Another concern is our method of classifying CFTR genotype. This was done a priori, with CFTR mutations associated with severely reduced CFTR production (members of the functional Class I, II, and III) grouped as "severe" and those that are likely to have more CFTR function (members of functional Class IV, V, and Unclassified) grouped as "mild." We believe this classification method is justified based on the knowledge that CFTR influences GSH transport and the substantial evidence that the CFTR genotypes we classified as mild are associated with a milder pulmonary phenotype (3, 5, 2125). Finally, our analysis was based on the hypothesis that low levels of CFTR leads to reduced GSH in the ELF, but there are alternative explanations. Recent studies suggest that CFTR expression may be reduced in the setting of increased oxidant stress (53, 54). Because a moderate amount of CFTR expression is associated with a milder CF phenotype (36), genetic variation in GCLC that leads to improved intracellular oxidant defense may lead to a milder CF phenotype as a result of increased CFTR expression. Alternatively, severe CFTR genotypes leading to worse clinical disease may result in lower GSH as a result of increased lung inflammation and impaired nutrition. Further studies are necessary to understand the interaction between GCLC, GSH, and CFTR as modifiers of the CF phenotype. In conclusion, this study found an association between variants in GCLC, a novel candidate gene that influences GSH production, and CF lung function. This effect was observed only in patients with a mild CFTR genotype. This is of clinical interest in terms of predicting which patients with CF will develop more severe lung disease and the potential role the GCLC genotype may play in the pharmacogenetics of GSH replacement therapy.
The authors thank the patients, their families, and the clinic staff who assisted in the completion of this study.
* These investigators contributed equally to this article. Supported by National Institute of Health K23 HL/70849-01 and the United States' Cystic Fibrosis Foundation (EFM), National Institutes of Health grants 1P21ES04696 and 1P30ES07033 (TJK), and the Canadian Cystic Fibrosis Foundation. This article has an online supplement, which is accessible from this issue's table of contents at www.atsjournals.org Originally Published in Press as DOI: 10.1164/rccm.200508-1281OC on May 11, 2006 Conflict of Interest Statement: None of the authors has a financial relationship with a commercial entity that has an interest in the subject of this manuscript. Received in original form August 18, 2005; accepted in final form May 11, 2006
This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||