Skip to main content

Polymorphisms in the glutathione pathway modulate cystic fibrosis severity: a cross-sectional study



Cystic fibrosis (CF) clinically manifests with various levels of severity, which are thought to be modulated by mutations in the cystic fibrosis transmembrane conductance regulator gene (CFTR), modifier genes, and the environment. This study verified whether polymorphisms in modifier genes associated with glutathione (GSH) metabolism influence CF severity.


A cross-sectional study of 180 CF patients was carried out from 2011 to 2012. We analyzed CFTR mutations, polymorphisms (GSTM1 and GSTT1 deletions, GSTP1 + 313A > G, GCLC-129C > T, and GCLC-3506A > G) in modifier genes and CF clinical severity as assessed by 28 clinical and laboratory variables.


Significant associations were found between modifier gene polymorphisms and particular phenotypes or genotype changes. These included GCLC-129C > T with a higher frequency of the Pseudomonas aeruginosa mucoid to CC genotype (p = 0.044), and GCLC-3506A > G with a higher frequency of the no-mucoid P. aeruginosa (NMPA) to AA genotype (p = 0.012). The GSTT1 deletion was associated with a higher frequency of the NMPA to homozygous deletion (p = 0.008), GSTP1 + 313A > G with a minor risk of osteoporosis (p = 0.036), and patient age ≤ 154 months (p = 0.044) with the AA genotype. The Bhalla score was associated with GCLC-3506A > G (p = 0.044) and GSTM1/GSTT1 deletion polymorphisms (p = 0.02), while transcutaneous hemoglobin oxygen saturation levels were associated with GSTT1 deletions (p = 0.048).


CF severity is associated with polymorphisms in GSH pathways and CFTR mutations.

Peer Review reports


Cystic fibrosis (CF) presents with broad phenotypic variability, even in patients with identical mutations in the causative gene, cystic fibrosis transmembrane conductance regulator (CFTR) [1]. Explanations for this include environmental factors [2], medical management [3], nutritional status [4], emotional maladjustments [5], socioeconomic status [3], CFTR mutations [1], and modifier genes [1, 3, 6]. In this context, CF modifier genes have been studied with the aim of increasing chlorine transport and/or controlling pulmonary inflammation and infection [69].

Our group studied CF severity in association with several modifier genes including polymorphisms in the genes: MBL-2, TGF-β1, CD14[10], ACE[11], ADRB2[12], TCF7L2[13], ADRA2A[14], COX-2[15] and IFRD1[16]. These polymorphisms were associated with clinical variables including lung and digestive disease.

Glutathione (GSH) is a tripeptide composed of L-cysteine, L-glutamic acid, and glycine. It is a crucial part of the intracellular defense system, which protects the epithelium against the injuries and inflammation [17] common to CF that are caused by oxidation [18]. As polymorphisms can alter the GSH metabolic pathway, genetic variations of this pathway have previously been studied in association with CF [1921].

The glutathione S-transferase (GST) family of enzymes comprises proteins with distinct genetic origins that form a detoxification system, which protects the human body against electrophilic compounds and oxidative stress [22]. The GST protein is responsible for combining compounds that cause oxidative stress with GSH. It is therefore possible that GST polymorphisms are involved in CF severity [18, 22], especially with regard to pulmonary disease.

Genetic variants of the GST genes include glutathione S-transferase mu 1 (GSTM1) located on chromosome 1p13.3, and glutathione S-transferase theta 1 (GSTT1) on chromosome 22q11.23 [23], which both exhibit polymorphic deletions [22, 24]. The null GST allele does not encode a GST protein, so homozygous genotypes are associated with increased CF clinical severity [25, 26]. The glutathione S-transferase pi gene (GSTP1) on chromosome 11q13 [23] is associated with xenobiotic metabolism and susceptibility to cancer and other diseases [22]. Its most commonly studied polymorphism is an A → G base exchange at the +313 position (substituting isoleucine by valine at codon 105) [27].

The glutamate-cysteine ligase, catalytic subunit gene (GCLC) on chromosome 6p12 [23] encodes the catalytic subunit of glutamate-cysteine ligase (GCL), which is the first limiting enzyme in GSH synthesis [28]. The GCL holoenzyme is a heterodimer of approximately 104 kDa composed of catalytic-GCLC and regulatory-GCLR subunits [18]. The -129C > T and -3506A > G polymorphisms of GCLC are located in the promoter region and are responsible for reduced production of GSH [18, 28].

Of these genes, GSTP1 is associated with hepatic disease [19] and infection [20], GSTM1 with greater CF clinical severity [21], GSTT1 with no CF clinical variables, while GCLC has not been previously studied in relation to CF. However, as the action of the GSH protein is closely related to that of CFTR [29], it is conceivable that GCLC and GST polymorphisms influence CF severity [1921, 26, 30]. This study therefore aimed to determine whether genetic polymorphisms in the GSH metabolic pathway are associated with CF severity under different phenotypes of the disease.


This cross-sectional study was conducted in a university center for CF care between 2011 and 2012. Two hundred and fifteen patients were selected for the study, of which 35 were excluded for not signing the consent form or because of a lack of clinical data for statistical analysis. CF diagnosis was confirmed if levels of chloride in the sweat exceeded 60 mEq/L and by CFTR mutation screening when possible. CF patients, with no identified CFTR mutation or with one CFTR mutation screened, were classified as CF disease, considering: (i) all patients had levels of chloride in the sweat exceeded 60 mEq/L; (ii) CF clinical symptoms were diagnosed in all patients as: chronic obstructive pulmonary disease, bacteria in sputum, spirometry with obstruction values for forced expiratory volume in the first second (FEV1%), associated comorbidities (i.e. osteoporosis, nasal polyps, diabetes mellitus and pancreatic insufficiency); (iii) the dosage of active CFTR in epithelium via rectal biopsy was performed - all patients included had abnormal values for biopsy – absence of active CFTR was found; (iv) nasal potential was realized in some patients - all values were changed – but the comparison was not performed, taking into account a control standard curve, being an inconclusive data. By this method was possible to exclude Cystic Fibrosis Related Diseases.

No patients were diagnosed by a neonatal screening test. Patient DNA was obtained by phenol-chloroform extraction and 50 ng/mL was used for analysis as evaluated by a GE NanoVue™ Spectrophotometer (GE Healthcare Biosciences, Pittsburgh, PA, USA).

Clinical variables

Several clinical variables were employed, including Shwachman-Kulczycki, Kanga and Bhalla clinical scores [31]; body mass index (BMI) [for patients older than 19 years, the BMI = weight/(height)2 formula was used, while remaining patients used the WHO ANTHRO program (children 0–5 years of age) or the WHO ANTHRO PLUS program (children 5–19 years of age)]; patient’s age (≤154 and >154 months); time to diagnosis (≤24 and >24 months); time of first clinical symptoms (digestive: ≤3 and >3 months; pulmonary: ≤6 and >6 months); time to first colonization by Pseudomonas aeruginosa (≤31 and >31 months); bacteria in the respiratory airways: mucoid P. aeruginosa and no mucoid P. aeruginosa, Achromobacter xylosoxidans, Burkholderia cepacia and Staphylococcus aureus - the positive status was evaluated considering chronic infection (patients in whom more than 50% of the preceding 12 months was culture positive) + intermittent infection (patients with less than 50% of cultures positive). A patient was negative considering as free of bacterium (when no bacterium was grown from samples in the previous 12 months, despite a history of prior colonization) + never infected (patients in whom the bacterium) has never been cultured, i.e. this consensus was formulated for P. aeruginosa, but in our data was used for all bacteria [32]; transcutaneous hemoglobin oxygen saturation (SpO2) and spirometry variables.

Spirometry was performed in patients older than seven years of age with the CPFS/D spirometer (MedGraphics, Saint Paul, MN, USA) and data were recorded using the PF BREEZE software version 3.8B for Windows 95/98/NT [33]. The following variables were included: forced vital capacity [FVC(%)]; forced expiratory volume in the first second [FEV1(%)], the ratio between FEV1 and FVC(%) [FEV1/FVC(%)]; and forced expiratory flow between 25 and 75% of the FVC [FEF25-75%]. The data was analyzed considering international curves values for spirometry tests [34, 35].

The comorbidities analyzed were nasal polyps, osteoporosis, meconium ileus, diabetes mellitus, and pancreatic insufficiency. This study was approved by the Institutional Ethics Committee from the Faculty of Medical Sciences, University of Campinas (#528/2008), and all included patients or their parents signed a consent form before beginning the study.

CFTRmutation identification

CFTR mutation identification was performed by polymerase chain reaction (PCR) for F508del and the fragment-length polymorphism method for G542X, R1162X, R553X, G551D, and N1303K mutations. Some CF mutations were identified by sequencing or Multiplex Ligation-dependent Probe Amplification (MLPA) analysis: S4X, 2183A > G, 1717-G > A, and I618T. A MegaBace1000® sequencer (GE Healthcare Biosciences) was used for sequencing and MLPA.

The CFTR genotype was used as a correction factor for statistical analysis. All class I, II or III mutations, but not class IV mutations (P205S and R334W), identified were included in statistical analysis.

Identification of polymorphisms associated with GSH metabolic pathway genes

Polymorphism identification was carried out using PCR analysis. For GSTM1 and GSTT1 genes, a multiplex PCR reaction was performed using the CYP1A1 gene as an internal amplification control [36]. GCLC-129C > T, -3506A > G [18, 28] and GSTP1 + 313A > G [27] polymorphisms were identified by PCR followed by enzymatic digestion.

Statistical analysis

Statistical analysis was performed using Statistical Package for Social Sciences (SPSS) software version 21.0 (SPSS Inc., Chicago, IL, USA), Epi Info version 6.0 [37] and R version 2.12 (Comprehensive R Archive Network, 2011). GPower software [38] was used to calculate the statistical power, which was required to be above 80% for analysis.

Statistical tests included the analysis of variance (ANOVA) and the chi-square (χ2) test (Odds Ratio -OR) for GSTM1 and GSTT1, and the t-test and Fisher’s exact test for GCLC-129C > T, GCLC-3506A > G and GSTP1 + 313A > G polymorphisms. To avoid spurious data caused by the performance of multiple tests [39], the significance level (α) was adjusted by the Bonferroni correction (α corrected = 0.05/number of tests → 0.05/4 = 0.0125). The value of α was corrected considering clinical marker analysis of the same group of patients, taking into account, the CFTR mutation genotype.

Data distribution showing a high standard deviation was analyzed in groups distributed according to median value. Variables that were adjusted by median to short (more severe) and longtime were patient’s age, time to diagnosis, onset of pulmonary and digestive symptoms, and time to the first isolation of P. aeruginosa.

Analyses were performed of four cohorts: (i) all patients with CF (n = 180); (ii) patients with no identified CFTR mutation (n = 44); (iii) patients with an identified mutant CFTR allele (Class I, II and/or III) (n = 51); and (iv) patients with two identified CFTR mutations (Class I, II and/or III) (n = 85). For (ii) and (iii) groups, a second analysis was performed. In this case, CF patients with pancreatic sufficiency (PS) were excluded. Patients with mutations Class I, II and III for CFTR gene have severe disease, strongly associated with pancreatic insufficiency (PI). Excluding PI patients was a method to associated different CFTR mutation groups with no atypical CF – associated less severe mutation (Class IV, V and VI). After exclusion, we have in (ii) and (iii) groups, respectively, 35 and 43 CF patients.

Results and discussion

One of the most intriguing aspects of CF is that patients with the same CFTR genotype can present with phenotypic differences [40]. At our CF center, all patients receive free medication provided by the state, have a similar socioeconomic status, share similar Class I, II and/or III mutations, receive support from the Cystic Fibrosis Association (, and there are no severe cases of malnutrition. This therefore makes our sample more phenotypically homogeneous for studies involving gene modulation characteristics.

Variations in CF severity can be associated with a modifier gene, such as those associated with oxidative stress [1921, 26, 30] that are related to chronic obstructive pulmonary disease (COPD) [18]. The COPD pathophysiology is similar, in some aspects, to CF in that it involves cellular responses, inflammatory mediators, and oxidative stress [41]. However, there is no mention in the scientific literature of GCLC polymorphisms as clinical modulators of CF severity, and is necessary new studies to illuminate about GST genes and CF severity.

One of the main functions of GSH is to detoxify xenobiotics and their metabolites, and this function is dependent on GST proteins. The GST gene family has been linked with several diseases [22], as GSTM1, GSTT1 and GSTP1 polymorphisms were found to be associated with cancer, drugs, chemotherapy resistance [42], and respiratory diseases such as asthma [30]. For example, expression of the variant form of GSTP1 (where isoleucine is substituted for valine at codon 105) results in lower enzymatic activity, which is a risk factor for the development of cancer and pulmonary diseases such as CF [43].

The effects of GSTM1, GSTT1 and GSTP1 polymorphisms on spirometry were previously investigated in 1,940 children (aged 8–11 years) [44]. The null GSTM1 genotype was associated with a decrease in annual FVC(%) and FEV1(%) gain; likewise, homozygosity for the GSTP1 allele was linked with slower spirometric gain for the same markers. The GSTM1 and GSTP1 genotypes therefore appear to be associated with spirometric evolution, and could increase the severity of diseases of pulmonary obstruction, depending on the genotype and gene combination.

Table 1 shows the GCLC, GSTM1, GSTT1 and GSTP1 polymorphism distribution according to genotype in the present study. The -129C > T polymorphism in the promoter region of GCLC stimulates different responses to oxidative stress by decreasing GSH production and reducing cellular antioxidant capacity [45]. In the present study, it was associated with a higher frequency of the mucoid P. aeruginosa to CC genotype for GCLC-129C > T polymorphism in patients with one CFTR mutation identified (Table 2; p = 0.044). This association may be related to the lower GCLC protein expression in CC genotypes, which reduces circulating GSH levels. The T allele is also associated with increased GSH expression, as described in protein expression studies on cardiovascular disease [45]. In the cited literature, we only found one study that associated the GCLC polymorphism with CF severity. In this previous study, the GAG micro-satellite GCLC polymorphism was analyzed in 440 CF patients, and CFTR mutations of lower gravity and highest number of GAG repeats in the GCLC gene were associated with higher values of FEV1(%) [20].

Table 1 Distribution of GCLC , GSTM1, GSTT1 and GSTP1 polymorphisms
Table 2 Polymorphisms in modifier genes associated with categorical variables of cystic fibrosis severity

The GCLC-3506A > G polymorphism is not in Hardy-Weinberg equilibrium as shown in Table 3, which also shows the complete genotypic characteristics of GCLC, GSTM1, GSTT1, and GSTP1 polymorphisms and CFTR mutations in CF patients with regard to chromosomal position, polymorphism location within the gene, and minor allele frequency. GCLC-3506A > G was associated with a higher frequency of the no mucoid P. aeruginosa to AA genotype, and with a lower frequency of the no mucoid P. aeruginosa to AG + GG genotype group in patients with one CFTR mutation identified (Table 2; p = 0.012) and higher Bhalla score values (without taking CFTR mutation into account; p = 0.044).

Table 3 Genotyping of GCLC , GSTM1 , GSTT1 , and GSTP1 polymorphisms and CFTR mutations

The Bhalla score is associated with an impairment of the pulmonary parenchyma structure and higher values characterize major changes in thoracic tomography. Unexpectedly, we also found that the greatest expression of the A allele in the GCLC-3506A > G polymorphism did not protect against no mucoid P. aeruginosa colonization. However, protection against lung deterioration was evident when we considered the Bhalla score. This score was also associated with GSTM1/GSTT1 deletions (p = 0.02), with a lower frequency of heterozygous compared with homozygous deletions. Moreover, the GSTT1 deletion was found to be associated with SpO2 values (p = 0.048; Table 4).

Table 4 Polymorphisms in modifier genes associated with numerical variables of cystic fibrosis severity

GCLC haplotype analysis for GCLC-129C > T and GCLC-3506A > G showed association for A. xylosoxidans and CC + AA genotypes (OR = 17.9; CI95% = 2.781-411.6; Table 5).

Table 5 GCLC-129C > T and GCLC-3506A > G haplotype polymorphisms in modifier genes associated with categorical variables of cystic fibrosis severity

The present study found that the AA genotype of the GSTP1 + 313A > G polymorphism was associated with a low risk of osteoporosis (p = 0.036; with two CFTR mutations identified) as a protective factor and with young age ≤ 154 months (p = 0.044; without taking the CFTR gene into account) as a risk factor. The G allele, however, is responsible for increased GSTP1 expression. The presence of osteoporosis is influenced by several different factors, including mutations in the CFTR gene, the environment, modifier genes, and increased life expectancy [46]. In this context, in our data, the A allele is protective against osteoporosis, and is increased among young patients with unresolved CFTR mutation genotype. The osteoporosis frequency is shown in Table 6.

Table 6 Age’s distribution and osteoporosis among CFTR mutation groups

The role of the GSTP1 polymorphism in CF hepatic disease has previously been analyzed [19]. The authors noted that CFTR protein expression was limited in liver epithelium; however, recent discoveries indicate that CFTR modulates the transport of GSH, creating a dysfunction in the antioxidant defense [47]. Of the liver detoxifying enzymes, GST plays a major role in protection against oxidative stress. The impact of GSTM1 and GSTP1 was also previously assessed in 106 CF patients where it was verified that the frequency of the GG genotype for the GSTP1 + 313A > G polymorphism was significantly higher in CF patients with hepatic disease. This genotype was associated with an eight-fold increase in hepatic disease risk in patients younger than six years of age. These findings suggest that the identification of this polymorphism may have prognostic and awareness values for the treatment of CF patients with hepatic disease.

Considering the importance of the glutathione transport versus CFTR protein-mediated, patients with residual CFTR protein expression would have better performance in the extracellular oxidative stress response being favorable for the passage of GSH to the outside by residual CFTR activity. However, CF patients with two mutations screened in CFTR gene have principally alternate routes for the passing of GSH. Even taking into account that the most of GSH is transferred to the external environment via CFTR, in cases of residual CFTR (mutations Class IV, V and VI) would be modified slightly in relation to the GSH activity, since it is known that under 5% of CFTR expression occurs for minor severe CFTR mutations Classes, and approximately 65% of GSH passage occurs via CFTR, we had a percentage response to the GSH presence in external environment of at most 3.25% in cases of residual CFTR. Considering this factor, the analysis excluding the presence of PI patients enables better grouping of patients and optimizes the response of the associations found in our study.

Most studies analyzing the GSTP1 gene related it to cancer and other diseases [22, 30, 44, 48]. For example, the AA genotype of the GSTP1 + 313A > G polymorphism was shown to offer protection against asthmatic symptoms [22]. Indeed, the GSTP1 polymorphism was not previously found to affect pulmonary function in CF patients [30]. In an analysis of different genes involved in GST, there were no differences in GST activity and antioxidant levels observed between CF patients and controls. However, GST activity was lower in P. aeruginosa-infected CF children with severe clinical symptoms, as was the frequency of the GSTP1 + 313A > G polymorphism AA genotype in uninfected (75%) compared with infected (33%) children [21]. It is possible that GST activity and GSTP1 genotype play an important role in P. aeruginosa infection in CF patients. In support of this, the G allele of the GSTP1 gene appears to be associated with an increased risk of severe pulmonary disease [21]. However, in a previous investigation into GSTM1 and GSTP1 polymorphisms in patients with CF and COPD, no significant associations were found between GSTM1 activity and pulmonary disease severity. An analysis of genotypic combinations for GSTM1 and GSTP1 polymorphic loci showed that changes in GSTP1 activities produced adverse effects in patients with COPD. Although GSTM1 gene deletions may not themselves be implicated in pathogenesis, they may aggravate the disease in combination with GSTP1 polymorphisms. Perhaps the strongest performance for the GSTP1 gene in CF may result from the primary expression of this GST in the airways [48].

The present study showed that the homozygous deletion in GSTT1 was a no mucoid P. aeruginosa risk factor in the no CFTR mutation group (p = 0.008; Table 2) and a protective factor for low values of SpO2. GSTT1 expression is likely to act in the inflammatory response of the pulmonary parenchyma. As chronic airway infection by no mucoid P. aeruginosa is associated with greater clinical severity [49], the GSTT1 polymorphism may be associated with the presence of P. aeruginosa through different mechanisms, including a low antioxidant response leading to further pulmonary degradation and the formation of a favorable environment for no mucoid P. aeruginosa colonization or infection.

The mechanism of gene action that determines which bacteria can colonize the lungs of CF patients is not fully understood. Similarly, it is also unclear which microorganisms are risk factors for the disease. Therefore, confirmation of a gene acting as modulator of an important metabolic pathway, such as GSH, may open up novel ways to identify the genetic factors that determine the severity of pulmonary disease. Future pharmacogenetic studies could then use this knowledge to provide new CF therapies.

Many previous studies have revealed that polymorphisms of GSTM1 and GSTT1 are associated with cancer [22, 24, 50, 51], but few have been conducted in CF. Fifty-three children with CF were studied by Hull and Thomson [26], of which 26 with the GSTM1 null allele had a significantly lower Shwachman-Kulczycki score. This supports the hypothesis that inflammation in CF contributes to tissue injury. Indeed, GSTM1 null alleles can be a risk factor for pulmonary diseases in individuals with a reduced ability to deal with oxidants. There is also evidence that a high level of oxidative stress in the lungs of CF patients is caused by the release of reactive oxygen species by neutrophils [26]. In the present study, we found that expression of only one allele of GSTM1 and GSTT1 polymorphisms was associated with a low Bhalla score in patients with no CFTR mutation identified.

An interesting aspect was the high frequency of PS patients. The presence of PS occurred at exactly 20% of the sample. However, there was no difference distribution between the groups of patients with CF taking into account CFTR mutations groups (p = 0.621). Patients with two mutations identified in CFTR gene had 22.36% (19/85) of PS, values close to the other groups of patients [one identified mutation and no mutation identified with, respectively, 15.7% (8/51) and 20.5% (9/44)].

The PI is an important clinical marker of CF and is considered associated with the severity of disease and severe CFTR mutations (Class I, II and/or III). Studies considering populations of patients with CF, as performed by the Cystic Fibrosis Foundation give the prevalence of PI ranging from 5-10%. In our study, the high prevalence of PI may be associated with: (i) presence of higher frequency of mutations Class IV, V and/or VI, (ii) presence of modifier genes acting on the symptom of the disease, (iii) high miscegenation could be a protective factor for PI, (iv) environmental factor as an unknown protector.

The PI was used in statistical analysis as factor correction for no determination of CFTR mutation in CF groups with no or one CFTR mutation screened. After the patient exclusion to statistical analysis, all the previous positive associations were negative, except for GSTM1 null allele. The null allele was associated as protector factor for onset of digestive symptoms (OR = 0.134; CI = 0.023-0.606; Table 2).

One important aspect considered was the age. Before the statistical analysis, the age was considered between the CFTR mutations groups (p = 0.854). The same occurred for CFTR mutations groups + insufficiency pancreatic (p = 0.940) (Table 6). No positive association was find considering age.

The divergent immune response is associated with multiple factors that denote the CF complexity such as the multigenic response, environmental influences, and interaction between airway microorganisms [49, 52]. Clinically severe patients may have high initial inflammatory response, characterizing CF as a disease where inflammation occurs prior to infection [53]. Polymorphisms in genes that are involved in inflammation may be a risk factor for early severity of the disease [1], and patients with airways colonized by bacteria suffer early clinical deterioration and high levels of airway inflammation [54].

For the same population, a first study taking into account the same polymorphisms and clinical variables was performed. The previous data analyzed the genetic interaction among GST and GCLC polymorphisms, CFTR mutations and clinical markers. The data showed an interaction of GSTM1 and GSTT1 genes deletion, GSTP1* + 313A > G, and CFTR mutations (p = 0.008) and Bhalla clinical score by multifactor dimensionality reduction test. The Bhalla score is a computed tomography, which measures pulmonary involvement, therapeutic effects and selection of patients for transplantation, which detects anatomical changes of the lung parenchyma. The data published showed a first step to understand the complex mechanisms associated with the CF severity and modifier genes [55].

In the present study, we studied a CF population with complex clinical characteristics. By considering the different possible groupings of polymorphisms and clinical variables (Table 7) in relation to the CFTR gene, we performed various association studies. Supplementary data for GCLC-129C > T, GCLC-3506A > G, GSTM1 gene deletion, GSTT1 gene deletion, GSTM1/GSTT1 gene deletions and GSTP1* + 313A > G are shown in Tables 8, 9, 10, 11, 12, 13 and 14. Further multicenter studies should be conducted to verify the influence of modifier genes in different CFTR genotypes.

Table 7 Patient characteristics (n = 180)
Table 8 GCLC-129C > T polymorphism associated with CF clinical variables as distributed by CFTR mutation
Table 9 GCLC-3506A > G polymorphism in association with CF clinical variables as distributed by CFTR mutation
Table 10 GCLC-129C > T and GCLC-3506A > G polymorphisms by haplotype in association with CF clinical variables as distributed by CFTR mutation
Table 11 GSTM1 deletion polymorphism in association with CF clinical variables as distributed by CFTR mutation
Table 12 GSTT1 deletion polymorphism in association with CF clinical variables as distributed by CFTR mutation
Table 13 GSTM1/GSTT1 deletion polymorphism in association with CF clinical variables as distributed by CFTR mutation
Table 14 GSTP1 + 313A > G polymorphism in association with CF clinical variables as distributed by CFTR mutation

Study limitations: (i) CFTR mutation with no complete screening; (ii) short population of CF patients; (iii) spirometry test performed by transversal method and did no performed longitudinally; (iv) no measure of GSH activity or GST and GCLC proteins, taking into account the sample collection limitation in our center and time to process all data. Study highlights the data by: (i) one CF center collection – considering an admixed population, the CF patients from one center minimizes miscegenation factors. Another fact, is the similar environmental and the same access to treatment; (ii) high number of clinical markers evaluated provides better association and characterization of modifier genes action; (iii) complete CF diagnosis performed by different methods.


Our results show that, although a monogenic disease, CF is heavily influenced in its clinical characteristics, evolution and severity by polymorphisms in modifier genes. Nevertheless, there is still a long way before the dynamics of polymorphisms in genes active in the GSH metabolic pathway and involved in detoxification in CF are fully understood.

Another fact is the prevalence of PS and PI that should be considered in all studies in the future, being associated with different phenotype and genotype.



Cystic fibrosis


Cystic fibrosis transmembrane regulator


Glutamate-cysteine ligase, catalytic subunit


Glutathione S-transferase


Glutathione S-transferase mu 1


Glutathione S-transferase tetha 1


Glutathione S-transferase pi 1


No mucoid Pseudomonas aeruginosa


Transcutaneous hemoglobin oxygen saturation


Forced expiratory volume in 1 second


Forced vital capacity


Forced expiratory flow 25–75%


Body mass index


World health organization


Multiplex ligation-dependent probe amplification


Polymerase chain reaction


Statistical package for social science for windows


Pancreatic sufficiency


Pancreatic insufficiency


Chronic obstructive pulmonary disease.


  1. Cutting GR: Modifier genes in Mendelian disorders: the example of cystic fibrosis. Ann N Y Acad Sci. 2010, 1214: 57-69. 10.1111/j.1749-6632.2010.05879.x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Collaco JM, Blackman SM, McGready J, Naughton KM, Cutting GR: Quantification of the relative contribution of environmental and genetic factors to variation in cystic fibrosis lung function. J Pediatr. 2010, 157: 802-807. 10.1016/j.jpeds.2010.05.018.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Simmonds NJ, MacNeil SJ, Cullinan P, Hodson ME: Cystic fibrosis and survival to 40 years: a case–control study. Eur Respir J. 2010, 36: 1277-1283. 10.1183/09031936.00001710.

    Article  CAS  PubMed  Google Scholar 

  4. Matel JL, Milla CE: Nutrition in cystic fibrosis. Semin Respir Crit Care Med. 2009, 30: 579-586. 10.1055/s-0029-1238916.

    Article  PubMed  Google Scholar 

  5. Casier A, Goubert L, Theunis M, Huse D, Baets F, Matthys D, Crombez G: Acceptance and well-being in adolescents and young adults with cystic fibrosis: a prospective study. J Ped Psychol. 2011, 36: 476-487. 10.1093/jpepsy/jsq111.

    Article  Google Scholar 

  6. Stanke F, Becker T, Kumar V, Hedtfeld S, Becker C, Cuppens H, Tamm S, Yarden J, Laabs U, Siebert B, Fernandez L, Macek M, Radojkovic D, Ballmann M, Greipel J, Cassiman JJ, Wienker TF, Tümmler B: Genes that determine immunology and inflammation modify the basic defect of impaired ion conductance in cystic fibrosis epithelia. J Med Genet. 2011, 48: 24-31. 10.1136/jmg.2010.080937.

    Article  CAS  PubMed  Google Scholar 

  7. Merlo CA, Boyle MP: Modifier genes in cystic fibrosis lung disease. J Lab Clin Med. 2003, 141: 237-242. 10.1067/mlc.2003.29.

    Article  CAS  PubMed  Google Scholar 

  8. Slieker MG, Sanders EAM, Rijkers GT, Ruven HJT, Van Der Ent CK: Disease modifying genes in cystic fibrosis. J Cyst Fibros. 2005, 4: 7-13.

    Article  CAS  PubMed  Google Scholar 

  9. Davies JC, Griesenbach U, Alton E: Modifier genes in cystic fibrosis. Pediatr Pulmonol. 2005, 39: 383-391. 10.1002/ppul.20198.

    Article  CAS  PubMed  Google Scholar 

  10. Faria EJ, Faria IC, Ribeiro JD, Ribeiro AF, Hessel G, Bertuzzo CS: Association of MBL2, TGF-beta1 and CD14 gene polymorphisms with lung disease severity in cystic fibrosis. J Bras Pneumol. 2009, 35 (4): 334-342. 10.1590/S1806-37132009000400007.

    Article  PubMed  Google Scholar 

  11. Marson FAL, Bertuzzo CS, Hortencio TD, Ribeiro JD, Bonadia LC, Ribeiro AF: The ACE gene D/I polymorphism as a modulator of severity of cystic fibrosis. BMC Pulm Med. 2012, 12: 41-10.1186/1471-2466-12-41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Marson FAL, Bertuzzo CS, Ribeiro AF, Ribeiro JD: Polymorphisms in ADRB2 gene can modulate the response to bronchodilators and the severity of cystic fibrosis. BMC Pulm Med. 2012, 12 (1): 50-10.1186/1471-2466-12-50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Furgeri DT, Marson FAL, Ribeiro AF, Bertuzzo CS: Association between the IVS4G > T mutation in the TCF7L2 gene and susceptibility to diabetes in cystic fibrosis patients. BMC Res Notes. 2012, 5: 561-10.1186/1756-0500-5-561.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Marson FAL, Rezende LM, Furgeri DT, Ribeiro AF, Ribeiro JD, Bertuzzo CS: ADRA2A is a Cystic Fibrosis Modifier Gene. Int J Genet. 2013, 5 (1): 125-131.

    Article  CAS  Google Scholar 

  15. Marson FAL, Marcelino ARB, Ribeiro AF, Ribeiro JD, Bertuzzo CS: COX-2 Gene Polymorphisms: Genetic Determinants of Cystic Fibrosis Comorbidities. Int J Genet. 2013, 5 (1): 132-138.

    Article  CAS  Google Scholar 

  16. Marson FAL, Marcelino ARB, Rezende LM: The IFRD1 (57460C > T polymorphism) gene: a negative report in cystic fibrosis clinical severity. J Mol Genet Med. 2013, in press

    Google Scholar 

  17. Liu S, Li B, Zhou Y, Zhong N, Ran P: Genetic analysis of CC16, OGG1 and GCLC polymorphisms and susceptibility to COPD. Respirology. 2007, 12: 29-33. 10.1111/j.1440-1843.2006.00958.x.

    Article  CAS  PubMed  Google Scholar 

  18. Galli F, Battistoni A, Gambari R, Pompella A, Bragonzi A, Pilolli F, Iuliano L, Piroddi M, Dechecchi MC, Cabrini G, Working Group on Inflammation in Cystic Fibrosis: Oxidative stress and antioxidant therapy in cystic fibrosis. Biochim Biophys Acta. 2012, 1822 (5): 690-713. 10.1016/j.bbadis.2011.12.012.

    Article  CAS  PubMed  Google Scholar 

  19. Henrion-Caude A, Flamant C, Roussey M, Housset C, Flahault A, Fryer AA, Chadelat K, Strange RC, Clement A: Liver disease in pediatric patients with cystic fibrosis is associated with glutathione S-transferase P1 polymorphism. Hepatology. 2002, 36: 913-917.

    Article  CAS  PubMed  Google Scholar 

  20. McKone EF, Shao J, Frangolias DD, Keener CL, Shephard CA, Farin FM, Tonelli MR, Pare PD, Sandford AJ, Aitken ML, Kavanagh TJ: Variants in the glutamate-cysteine-ligase gene are associated with cystic fibrosis lung disease. Am J Resp Crit Care Med. 2006, 174: 415-429. 10.1164/rccm.200508-1281OC.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Feuillet-Fieux MN, Nguyen-Khoa T, Loriot MA, Kelly M, De Villartay P, Sermet I, Verrier P, Bonnefont JP, Beaune P, Lenoir G, Lacour B: Glutathione S-transferases related to P. aeruginosa lung infection in cystic fibrosis children: preliminary study. Clin Biochem. 2009, 42: 57-63. 10.1016/j.clinbiochem.2008.09.116.

    Article  CAS  PubMed  Google Scholar 

  22. Strange RC, Spiteri MA, Ramachandran S, Fryer AA: Glutathione-S-transferase family of enzymes. Mutat Res. 2001, 482: 21-26. 10.1016/S0027-5107(01)00206-8.

    Article  CAS  PubMed  Google Scholar 

  23. National Center for Biotechnology Information. []

  24. Amorim LMF, Rossini A, Mendonça G, Lotsch P, Simão TA, Gallo CM, Pinto L: CYP1A1, GSTM1, and GSTT1 polymorphisms and breast cancer risk in Brazilian women. Cancer Lett. 2002, 181: 179-186. 10.1016/S0304-3835(02)00058-7.

    Article  Google Scholar 

  25. Baranov VS, Ivaschenko T, Bakay B, Aseev M, Belotserkovskaya R, Baranova H, Malet P, Perriot J, Mouraire P, Baskakov VN, Savitskyi GA, Gorbushin S, Deyneka SI, Michnin E, Barchuck A, Vakharlovsky V, Pavlov G, Shilko VI, Guembitzkaya T, Kovaleva L: Proportion of the GSTM1 0/0 genotype in some Slavic populations and its correlation with cystic fibrosis and some multifactorial diseases. Hum Genet. 1996, 97: 516-520. 10.1007/BF02267078.

    Article  CAS  PubMed  Google Scholar 

  26. Hull J, Thomson AH: Contribution of genetic factors other than CFTR to disease severity in cystic fibrosis. Thorax. 1998, 53: 1018-1021. 10.1136/thx.53.12.1018.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Harries LW, Stubbins MJ, Forman D, Howard GCW, Wolf CR: Identification of genetic polymorphisms at the glutathione S-transferase Pi locus and association with susceptibility to bladder, testicular and prostate cancer. Carcinogenesis. 1997, 18: 641-644. 10.1093/carcin/18.4.641.

    Article  CAS  PubMed  Google Scholar 

  28. Koide S, Kugiyama K, Sugiyama S, Nakamura S, Fukushima H, Honda O, Yoshimura M, Ogawa H: Association of polymorphism in glutamate-cysteine ligase catalytic subunit gene with coronary vasomotor dysfunction and myocardial infarction. J Am Coll Cardiol. 2003, 41: 539-545.

    Article  CAS  PubMed  Google Scholar 

  29. Hudson VM: Rethinking cystic fibrosis pathology: the critical role of abnormal reduced glutathione (GSH) transport caused by CFTR mutation. Free Radic Biol Med. 2001, 1530: 1440-1461.

    Article  Google Scholar 

  30. Collaco JM, Cutting GR: Update on gene modifiers in cystic fibrosis. Curr Opin Pulm Med. 2008, 14: 559-566. 10.1097/MCP.0b013e3283121cdc.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Santos CIS, Ribeiro JD, Ribeiro AF, Hessel G: Critical analysis of scoring systems used in the assessment of cystic fibrosis severity: state of the art. J Bras Pneumol. 2004, 30: 286-298.

    Google Scholar 

  32. Döring G, Flume P, Heijerman H, Elborn JS, Consensus Study Group: Treatment of lung infection in patients with cystic fibrosis: current and future strategies. J Cyst Fibros. 2012, 11 (6): 461-479. 10.1016/j.jcf.2012.10.004.

    Article  PubMed  Google Scholar 

  33. American Thoracic Society (ATS).,

  34. Polgar G, Promadhat V: Pulmonary function testing in children: techniques and Standards. 1971, Philadelphia, PA: WB Saunders Company

    Google Scholar 

  35. Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, Culver BH, Enright PL, Hankinson JL, Ip MS, Zheng J, Stocks J: ERS Global Lung Function Initiative: multi-ethinic values for spirometry for the 3-95yr age range: the global lung function 2012 equations. Eur Respir J. 2012, 40: 1324-1343. 10.1183/09031936.00080312.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Abdel-Rahman SZ, El-Zein RA, Anwar WA, Au WW: A multiplex PCR procedure for polymorphic analysis of GSTM1 e GSTT1 genes in population studies. Cancer Lett. 1996, 107: 229-233. 10.1016/0304-3835(96)04832-X.

    Article  CAS  PubMed  Google Scholar 

  37. Dean AG, Dean JA, Coulombier D, Brendel KA, Smith DC, Burton AH, Dicker RC, Sullivan K, Fagan RF, Arner TG: EPI Info, Version 6.0: a word processing database and statistics program for epidemiology on microcomputers. 1994, Atlanta, Georgia: Center of Disease Control and Prevention

    Google Scholar 

  38. Faul F, Erdfelde E, Lang AG, Buchner A: GPower 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Meth. 2007, 39: 175-191. 10.3758/BF03193146.

    Article  Google Scholar 

  39. Drăghici S: Data analysis tools for DNA microarrays. 2003, New York: Chapman & Hall/CRC

    Book  Google Scholar 

  40. Drumm ML, Ziady AG, Davis PB: Genetic variation and clinical heterogeneity in cystic fibrosis. Annu Rev Pathol. 2012, 7: 267-282. 10.1146/annurev-pathol-011811-120900.

    Article  CAS  PubMed  Google Scholar 

  41. Ghezzi P: Role of glutathione in immunity and inflammation in the lung. Int J Gen Med. 2011, 4: 105-113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Hayes JD, Pulford DJ: The glutathione S-transferase supergene family: regulation of GST and the contribution of the isoenzymes to cancer chemoprotection and drug resistance. Crit Rev Biochem Mol Biol. 1995, 30: 445-600. 10.3109/10409239509083491.

    Article  CAS  PubMed  Google Scholar 

  43. Gresner P, Gromadzinska J, Wasowicz W: Polymorphism of selected enzymes involved in detoxification and biotransformation in relation to lung cancer. Lung Cancer. 2007, 57: 1-25. 10.1016/j.lungcan.2007.02.002.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  45. Campolo J, Penco S, Bianchi E, Colombo L, Parolini M, Caruso R, Sedda V, Patrosso MC, Cighetti G, Marocchi A, Parodi O: Glutamate-cysteine ligase polymorphism, hypertension, and male sex are associated with cardiovascular events. Biochemical and genetic characterization of Italian subpopulation. Am Heart J. 2007, 154: 1123-1130. 10.1016/j.ahj.2007.07.029.

    Article  CAS  PubMed  Google Scholar 

  46. Döring G, Conway SP: Osteoporosis in cystic fibrosis. J Pediatr. 2008, 84: 1-3. 10.1590/S0021-75572008000100001.

    Article  Google Scholar 

  47. Gould NS, Min E, Martin RJ, Day BJ: CFTR is the primary known apical glutathione transporter involved in cigarette smoke-induced adaptive responses in the lung. Free Radic Biol Med. 2012, 52 (7): 1201-1206. 10.1016/j.freeradbiomed.2012.01.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Korytina GF, Yanbaeva DG, Victorova TV: Polymorphism of glutathione S-transferase M1 and P1 in patients with cystic fibrosis and chronic respiratory diseases. Russ J Genet. 2004, 40: 314-320.

    Article  CAS  Google Scholar 

  49. Park JE, Yung R, Stefanowicz D, Shumansky K, Akhabir L, Durie PR, Corey M, Zielenski J, Dorfman R, Daley D, Sandford AJ: Cystic fibrosis modifier genes related to Pseudomonas aeruginosa infection. Genes Immun. 2011, 12: 370-377. 10.1038/gene.2011.5.

    Article  CAS  PubMed  Google Scholar 

  50. Mitrunen K, Jourenkova N, Kataja V, Eskelinen M, Kosma VM, Benhamou S, Vainio H, Uusitupa M, Hirvonen A: Glutathione S-transferase M1, M3, P1, and T1 genetic polymorphisms and susceptibility to breast cancer. Cancer Epidemiol Biomarkers Prev. 2001, 10: 229-236.

    CAS  PubMed  Google Scholar 

  51. Goloni-Bertollo EM, Biselli JM, Corrêa LCL, Maníglia JV, Rossit ARB, Ruiz MT, Pavarino-Bertelli EC: Avaliação da influência da nulidade dos genótipos GSTT1 e GSTM1 na carcinogênese em cabeça e pescoço. Rev Assoc Med Bras. 2006, 52: 365-368. 10.1590/S0104-42302006000500028.

    Article  PubMed  Google Scholar 

  52. Guss AM, Roeselers G, Newton ILG, Young CR, Klepac-Ceraj V, Lory S, Cavanaugh CM: Phylogenetic and metabolic diversity of bacteria associated with cystic fibrosis. ISME J. 2011, 5: 20-29. 10.1038/ismej.2010.88.

    Article  PubMed  Google Scholar 

  53. Balough K, McCubbin M, Weinberger M, Smits W, Ahrens R, Fick R: The relationship between infection and inflammation in the early stages of lung disease from cystic fibrosis. Pediatr Pulmonol. 1995, 20: 63-70. 10.1002/ppul.1950200203.

    Article  CAS  PubMed  Google Scholar 

  54. Sagel SD, Gibson RL, Emerson J, McNamara S, Burns JL, Wagener JS, Ramsey BW, Inhaled Tobramycin in Young Children Study Group; Cystic Fibrosis Foundation Therapeutics Development Network: Impact of Pseudomonas and Staphylococcus infection on inflammation and clinical status in young children with cystic fibrosis. J Pediatr. 2009, 154: 183-188. 10.1016/j.jpeds.2008.08.001.

    Article  PubMed  Google Scholar 

  55. Marson FAL, Bertuzzo CS, Secolin R, Ribeiro AF, Ribeiro JD: Genetic interaction of GSH metabolic pathway genes in cystic fibrosis. BMC Med Genet. 2013, 10 (14): 60-

    Article  Google Scholar 

Pre-publication history

Download references


We thank Luciana Cardoso Bonadia, Taís Daiene Russo Hortencio, Kátia Cristina Alberto Aguiar, Aline Gonçalves and Simoni Avansini for assistance in data collection and organization of ideas, Rodrigo Secolin for reviewing the use of English, and Maria Angela Ribeiro for spirometry analysis. Frauk Stanke, Andreas Hector and Harriet Corvol to manuscript correction. Margarida do Amaral and Karl Kunzelmann to provide CFTR-activity analysis.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Fernando Augusto de Fernando Augusto de Lima Marson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FALM contributed to the study conception and design, acquired, analyzed and interpreted the data, drafted the manuscript and revised it for intellectual content. CSB carried out the molecular genetic studies and drafted the manuscript. AFR drafted the manuscript and revised it for intellectual content. JDR approved the manuscript for publication. All authors read and approved the final manuscript.

Rights and permissions

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

Reprints and Permissions

About this article

Cite this article

Marson, F.A.d.F.A.d.L., Bertuzzo, C.S., Ribeiro, A.F. et al. Polymorphisms in the glutathione pathway modulate cystic fibrosis severity: a cross-sectional study. BMC Med Genet 15, 27 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: