Abstract

Oxidative stress resulting from an increased amount of reactive oxygen species and an imbalance between oxidants and antioxidants plays an important role in the pathogenesis of asthma. The present study tested the hypothesis that genetic susceptibility to allergic and nonallergic variants of asthma is determined by complex interactions between genes encoding antioxidant defense enzymes (ADE). We carried out a comprehensive analysis of the associations between adult asthma and 46 single nucleotide polymorphisms of 34 ADE genes and 12 other candidate genes of asthma in Russian population using set association analysis and multifactor dimensionality reduction approaches. We found for the first time epistatic interactions between ADE genes underlying asthma susceptibility and the genetic heterogeneity between allergic and nonallergic variants of the disease. We identified GSR (glutathione reductase) and PON2 (paraoxonase 2) as novel candidate genes for asthma susceptibility. We observed gender-specific effects of ADE genes on the risk of asthma. The results of the study demonstrate complexity and diversity of interactions between genes involved in oxidative stress underlying susceptibility to allergic and nonallergic asthma.

1. Introduction

Bronchial asthma (BA) is a common chronic inflammatory disease of the airways characterized by variable and recurring symptoms, reversible airflow obstruction, and bronchospasm [1]. There is a considerable body of evidence demonstrating that asthma is a multifactorial disease which results from complex interactions between susceptibility genes of small-to-modest effects and equally important environmental factors [2, 3].

In the recent years, the relationships between common genetic variants and BA risk are being reported with rapidly increasing frequency. Large-scale genome-wide association studies (GWAS) have been recently done to look for asthma susceptibility genes in ethnically diverse populations of the world [4, 5]. However, the findings obtained by GWAS were limited by the strongest associations of a few number of genetic variants that achieved genome-wide significance level. In addition, the genome-wide associations are found to be quite difficult to interpret with respect to disease pathogenesis [47]. Meanwhile, hundreds to thousands of genetic markers associated with a disease risk are not interpreted because they have not achieved the genome-wide significance level, thus accounting for the “missing heritability” of complex diseases [8]. This means that GWAS approach is powerless in the detecting genes of small-to-modest effects which represent a polygenic background of multifactorial disease. Additionally, genetic diversity of human populations [9], heterogeneity of disease pathogenesis [10], and especially a complexity of gene-gene interactions [11, 12] complicate our opportunities in unraveling the molecular mechanisms underlying complex diseases including asthma.

It is widely agreed that the expression of a disease phenotype may not accurately be predicted from the knowledge of the effects of individual genes because of complex nonlinear interactions between genes, including epistatic and additive interactions [8, 13]. To address this issue, several powerful data-mining approaches have been developed to identify susceptibility genes involved in such complex interactions [1416]. One of them, multifactor dimensionality reduction (MDR) method, was developed to reduce the dimensionality of multilocus information to improve the ability to detect genetic combinations that confer disease risk in relatively small samples [16, 17]. With set association analysis (SAA), contributions from multiple SNPs are combined by forming a sum of single-marker statistics, which results in a single genome-wide test statistic with high power [14].

An important task for a genetic epidemiologist utilizing a candidate gene approach is the selection of appropriate genes and SNPs for testing the disease association. Compared with studying individual genes, the inferences derived from a hypothesis-driven candidate pathway study are enhanced by allowing global conclusions about the involvement of entire biochemical pathway to the pathogenesis of disease [18]. Following this approach in our previous study [19], we pointed out the potential relevance of toxicogenomic mechanisms of BA in the modern world and proposed that genes for xenobiotic-metabolizing enzymes would be the most appropriate candidate genes for asthma susceptibility whose effects on the disease risk can be associated with exposure to air pollution. Due to the fact that air pollutants are the sources of reactive oxygen species (ROS), genes involved in oxidative stress can potentiate harmful effects of xenobiotics on the respiratory system [20, 21].

It is well known that oxidative stress resulting from an increased amount of ROS and an imbalance between oxidants and antioxidants plays a role in the molecular mechanisms underlying BA [2224]. We have demonstrated that genes of antioxidant defense enzymes (ADE), such as glutamate cysteine ligase (GCLM) [25], glutathione peroxidase (GPX1) [26], catalase (CAT) [27], myeloperoxidase (MPO) [28], NADPH oxidase (CYBA, p22phox subunit) [29, 30], NAD(P)H: quinone oxidoreductase type 1 (NQO1) [31] and microsomal epoxide hydrolase (EPHX1) [19], are important determinants of genetic susceptibility to asthma in Russians. In the present study, we tested the hypothesis that genetic susceptibility to both allergic and nonallergic asthma is determined by complex interactions between genes involved in oxidative stress. We performed for the first time a comprehensive analysis of genomic interactions between 34 ADE genes and 12 other candidate genes in order to identify gene-gene interactions in redox homeostasis underlying polygenic mechanisms of BA.

2. Materials and Methods

2.1. Study Population

The study protocol was approved by the Ethical Review Committee of Kursk State Medical University, and written informed consent was obtained from each participant before the study. The participants comprised a total of 429 unrelated individuals (215 patients with asthma and 214 healthy controls); all are ethnically Russians from Central Russia (mainly from the Kursk region). All study subjects were recruited from the Division of Pulmonology at the Kursk Regional Clinical Hospital between 2003 and 2004. Asthma was diagnosed by qualified pulmonologists on the basis of the WHO criteria, as described previously [19, 32]. The mean age of the patients with asthma (94 men and 121 women) was 43.3 years (ranging from 16 to 67 years), and the mean age of the healthy subjects (105 men and 109 women) was 41.3 years (ranging from 17 to 84 years). Skin prick tests were conducted and total serum IgE levels were determined in all study subjects. Patients with positive skin prick tests and high level of total IgE were defined as patients with allergic asthma (64 men and 92 women). Asthmatics who showed either negative skin prick test results (wheal size: <5 mm) or a normal total IgE level (<0.35 IU) were considered to be patients with nonallergic asthma (29 men and 27 women). Data on allergic status were not available for three asthmatics. A strong positive family history of asthma was found in the case group (40.1%) in comparison with controls (6.7%).

2.2. Selection of Candidate Genes

The candidate genes for this study were selected according to the guidelines for genetic association studies proposed by Cooper and coauthors [33]. We used the following criteria to select ADE genes and their genetic polymorphisms satisfying our study’s purposes: (1) enzymes should represent key players involved in the regulation of redox processes; (2) enzymes should cover all biochemical pathways of redox homeostasis entirely, including enzymes possessing antioxidant activity (GPX1, SOD2, CAT, GSTM1, etc.) and those with prooxidant activity (i.e., ROS-generating enzymes such as CYBA, MPO, and CYP2E1); (3) enzymes should be expressed in the lung and/or airways (the expression patterns of the selected ADE genes in human tissues and organs are shown in Supplementary Material available online at http://dx.doi.org/10.1155/2014/708903); (4) SNPs should be functionally significant, whenever possible; and (5) minor allele frequency should be more than 5%. Following these criteria, 34 polymorphisms of 24 ADE genes have been selected from published literature and public databases.

2.3. DNA Extraction and Genotyping

Genomic DNA of all study participants was isolated from 5–10 mL of peripheral blood samples, collected in K3-EDTA tubes by venipuncture, and maintained at −20°C until processed. Twenty-five of the selected gene polymorphisms had been genotyped in our previous studies [19, 2531]. In the present study, another nine ADE gene polymorphisms such as GPX2 (rs17880492), GPX3 (rs2070593) GPX4 (rs713041), GSR (rs2551715), SOD2 (rs4880), SOD3 (rs2536512), PRDX1 (rs17522918), TXNRD1 (rs1128446), and FMO3 (rs2266782) have been genotyped. Additionally, we genotyped 12 polymorphisms of 9 candidate genes of asthma such as TNF (rs1800629), IL1B (rs16944), IL3 (rs40401 and rs31480), IL5 (rs2069812), CSF2RB (rs131840), IL9 (rs2069885), SCGB1A1 (rs11549442), and SERPINA1 (rs17580, rs143370956, and rs11568814). Majority of them have been reported to be associated with the risk of asthma and/or asthma-related phenotypes in Russians [3437]. A complete list of 46 studied SNPs is given in Table 1. Genotyping of the selected polymorphisms was done using restriction fragment length polymorphism assays according to the published protocols (genotyping protocols are available upon request). All of the genotyping was done blindly to the case-control status and the repeatability test was conducted for the 5% of total subjects, resulting in a 100% concordance rate.

2.4. Statistical Analysis

The concordance of genotypes prevalence in patients with asthma and healthy controls with values expected under Hardy-Weinberg equilibrium was assessed by Pearson’s chi-square test. The association between ADE gene polymorphisms and asthma was examined with binary logistic regression analysis with calculation of odds ratios (OR) and 95% confidence intervals (CI). The statistical calculations were done using Statistica for Windows (v8.0) software package (StatSoft; Tulsa, OK, USA). The statistical significance was established at the level. Bonferroni correction for values was applied in cases when multiple tests were performed.

Two bioinformatic approaches, SAA and MDR, were applied for the analysis of gene-gene interactions. The principle of SAA is described in detail elsewhere [14, 55] and implemented in a statistical program SUMSTAT (http://linkage.rockefeller.edu/ott/sumstat.html/). Briefly, the method combines the information derived from measurements of allelic/genotype association and departure from Hardy-Weinberg equilibrium into a single, genome-wide statistic. The markers with high Hardy-Weinberg disequilibrium (HWD) values in the control group are trimmed and are not considered for further analysis. For the remaining markers, effects of allelic/genotype association with disease and HWD values are then combined into a single Sum statistic [56]. values reported by the program were calculated by permutations. The number of permutation tests was set at 10000.

MDR is a flexible nonparametric and genetic model free method for analysis of high-order nonlinear or nonadditive gene-gene interactions [15]. The method has been proposed to overcome limitation of logistic regression which deals with many factors simultaneously and fails to characterize epistatic models in the absence of main effects, due to the hierarchical model-building process leading to an increase in type II errors and decreased power [17]. The MDR method uses a constructive induction algorithm that converts two or more variables such as SNPs into a single attribute. In particular, SNPs are pooled into high and low risk group, effectively reducing the multifactor prediction from dimension to one dimension. Best models for each locus combination are selected by repeating the analysis for up to 10 seeds after shuffling the order of individuals and applying 10-fold cross-validation each time. Average of cross-validation consistency (CVC) together with training and test accuracy is calculated for each locus combination. CVC is defined as the number of times a particular interaction model is selected across 10 cross-validation datasets. We performed statistical calculations using MDR software (http://www.multifactordimensionalityreduction.org/). Statistical significance of the best models selected for each SNP combination was determined using 1000-fold permutation testing. The significance of the final MDR model was determined empirically by 1000 permutations using the Monte-Carlo procedure implemented into the MDRpt software (http://sourceforge.net/projects/mdr/). values for CVC were considered statistically significant at ≤0.05 levels. To visualize and interpret the results obtained from MDR, we used interaction dendrograms.

Both the SAA and MDR are limited by the identification of a few number of high penetrance interacting genes, whereas a larger portion of genes of low-to-moderate effects remain out of the analysis. To address this issue, we performed post hoc comparisons of two-locus genotype combinations (only for those SNPs which were found in gene-gene interaction models obtained by SAA and/or MDR methods) between the case and control groups to look for the genotype combinations which determine the risk of asthma. The observed associations were adjusted for multiple tests using Bonferroni procedure.

3. Results

3.1. Allele and Genotype Frequencies in Asthmatics and Controls

Allele and genotype frequencies of the studied genes are shown in Tables 2 and 3, respectively. After adjusting for multiple tests, the only statistically significant association was found between the IL5 C-703T polymorphism and BA. The −703CC genotype was found to be associated with the risk of allergic asthma (OR = 0.44; 95% CI 0.29–0.67; ). In gender-specific analysis, this association was seen in both men (OR = 0.50; 95% CI 0.27–0.94; ) and women (OR = 0.40; 95% CI 0.22–0.70; ) but did not reach a statistical significance after Bonferroni correction for multiple tests (). No association of this genotype was found with nonallergic asthma in both sexes.

3.2. Gene-Gene Interactions in Asthma Revealed by Set Association Analysis

Taking into account the polygenic basis of asthma, it was an important task to investigate high-order gene-gene interactions using specialized bioinformatics approach called set association analysis which captures the simultaneous effects of multiple genes and achieves a global view of gene action and interaction [14, 55]. For trimming, we considered values of departures from HWE (HWD values) exceeding the 99th percentile of chi-square (, ) in the control group [14]. There were no HWD values larger than 6.6 in the control subjects, so the trimming procedure to our dataset was avoided. To calculate single-locus test statistics, we used the difference in the distribution of genotypes for the th SNP between cases and controls (chi-square test for 2 × 3 tables). We tested up to sums () in allergic and nonallergic asthma, separately in men and women. When we added more SNPs to the , values tended to increase; that is, adding additional markers introduces noise to the . In allergic asthma, the smallest values were obtained for a sum of 5 SNPs (the GPX1 P198L, CAT −21A>T, EPHX1 H139R, GCLM −588C>T, and IL5 C-703T) for men and a sum of 3 SNPs (the EPHX1 Y113H, NQO1 R139W, and IL5 C-703T) for women (Figure 1). After 1000 permutation tests, the global significance levels of 0.0042 for men and 0.0001 for women were obtained. In nonallergic asthma, the smallest significance levels appeared for 5 SNPs (the GCLM −588C>T, GSR T>C, CAT −21A>T, CYBA −930A>G, and EPHX1 Y113H) in men and for 2 SNPs (the EPHX1 Y113H and GPX2 G>A) in women (Figure 2). The global significance levels of 0.0003 for men and 0.0001 for women were obtained.

3.3. Modeling for Gene-Gene Interactions in Asthma Using MDR Method

The MDR method was used for a purpose of modeling gene-gene interactions underlying allergic and nonallergic asthma in men and women. Firstly, we used an exhaustive search algorithm to evaluate all interactions among all possible subsets of the polymorphisms. Table 4 shows the cross-validation consistency and the prediction error for gene-gene interactions (from two- to four-locus interactions) obtained from MDR analysis in both allergic and nonallergic asthma. The only statistically significant (empirical ) three-locus model involving interactions between EPHX1 Y113H, IL5 C-703T, and GPX1 P198L loci was discovered. The model had a minimum prediction error of 40.9 and a maximum cross-validation consistency of 50% in allergic asthma in women . None of the rest -locus models in both allergic and nonallergic asthma showed a statistical significance in the MDR analysis using an exhaustive search algorithm, thereby motivating us to apply a forced search algorithm for further MDR analyses in order to build the best -locus models in men with allergic asthma and in both sexes with nonallergic asthma. Following this approach, we obtained one statistically significant (empirical ) four-locus model comprising interactions between CAT −21A>T, GPX2 G>A, GSR T>C, and IL5 C-703T in men with allergic asthma. The model had a minimum prediction error of 26.8 and a maximum cross-validation consistency of 100%.

Figure 3 shows the dendrograms illustrating high-order gene-gene interactions between the ADE loci in the pathogenetic variants of asthma in men and women. According to the figure, there is a strong difference in the structure of gene-gene interactions between men and women. In particular, synergistic interaction effect was found between the GSR T>C and IL5 C-703T loci in men. Moreover, the CAT −21A>T and GPX2 G>A gene polymorphisms had a strong antagonistic effect on the risk of allergic asthma in men. On the contrary, the hierarchical cluster analysis of the MDR data in women showed that the CYBA 640A>G, GPX4 C718T, and PON2 S311C gene polymorphisms have a strong synergistic interaction effect on the risk of allergic asthma. The EPHX1 Y113H and IL5 C-703T SNPs had a moderate antagonistic effect on the allergic asthma risk in women. Also, a relatively independent effect of the GPX1 P198L gene polymorphism on the risk of allergic asthma was seen.

On the next step, a forced search algorithm was applied to analyze all possible n-locus interactions in nonallergic asthma. The best 4-locus model including GPX1 P198L, GPX3 G/A, CYBA −930A>G, and FMO3 E158K polymorphisms was found in men (empirical ). The model had a minimum prediction error of 26.1 and a cross-validation consistency of 100%. The forced MDR analysis performed in women revealed a model including GPX1 P198L, GPX2 G>A, EPHX1 Y113H, and IL5 C-703T polymorphisms with a minimum prediction error of 28.1 and a cross-validation consistency of 100% (empirical ).

3.4. Post Hoc Association Analysis of Two-Locus Genotype Combinations

Then, we performed a post hoc comparison of genotype frequencies between the case and control groups with a focus on those ADE genes which were present in gene-gene interaction models obtained using SAA and MDR methods. Ten and nine two-locus combinations were found to be associated with allergic asthma in men and women, respectively (Table 5). However, only one genotype combination GPX4 718TC × CYBA 640AG achieved statistically significant inverse association with the risk of allergic asthma in women after adjustment for multiple tests (OR = 0.37; 95% CI 0.20–0.71; ). Twelve and five two-locus genotype combinations were found to be associated with the risk of nonallergic asthma in men and women, respectively (Table 6). Four two-locus genotype combinations showed statistically significant associations with nonallergic asthma in men after Bonferroni correction for multiple comparisons: GPX1 198PL × CAT −21AA (OR = 11.45; 95% CI 2.49–52.66; ), GSR TT × GCLM −588CT (OR = 11.58; 95% CI 3.07–43.72; ), CAT −21AA × CYBA −930GG (OR = 15.64; 95% CI 2.44–100.3; ), and GCLM −588CT × CYBA −930GG (OR = 6.71; 95% CI 2.5–17.96; ). One genotype combination EPHX1 113HH × IL5 −703CC showed a significant association with increased risk of nonallergic asthma in women (OR = 8.58; 95% CI 2.43–30.26; ).

4. Discussion

4.1. A Summary of the Study Findings

The main purpose of our study was to investigate a comprehensive contribution of ADE genes to genetic susceptibility to allergic and nonallergic variants of BA. The single-locus analysis revealed that none of the ADE genes was associated with the risk of asthma. However, using two bioinformatics approaches, we found multilocus gene-gene interactions which are associated with the risk of allergic and nonallergic asthma in men and women in a gender-specific manner. Further, post hoc analysis allowed revealing two-locus combinations of genotypes which are significantly associated with allergic and nonallergic asthma in both sexes. A majority of the susceptibility genes identified in our study represented antioxidant defense enzymes. Moreover, interactions between ADE genes varied across the pathogenetic variants of asthma and were different in men and women suggesting both genetic heterogeneity and gender-specific genetic effects in the disease susceptibility.

4.2. Genetic Heterogeneity of Asthma and Complexity of Genomic Interactions Underlying the Disease

The observed differences in gene-gene interactions between allergic and nonallergic variants of asthma demonstrate a genetic heterogeneity of the disease, a situation in which the same or similar phenotype of a complex disorder is caused by different susceptibility genes [57]. It is well known that genetic heterogeneity is the general feature of many common diseases [58] and may be explained at least partially by genetic differences between human populations [9]. Bronchial asthma is a typical example for complex multifactorial disease being characterized by genetic heterogeneity [59]. In fact, the models of gene-gene interactions in the pathogenetic variants of asthma overlap only partially, thereby reflecting, on the one hand, possible differences in the molecular mechanisms of allergic and nonallergic asthma and, on the other hand, the existence of shared genes that determine common susceptibility to the disease. In particular, three ADE genes such as GSR, EPHX1, and GPX1 showed significant interaction in both variants of asthma in both men and women (except for the GPX1 gene in nonallergic asthma in women); therefore, they can be considered as common susceptibility genes to asthma. While the IL5 and PON2 genes showed an association only with allergic asthma, none of the studied genetic polymorphisms was found to be associated exclusively with the risk of nonallergic asthma.

The results of gene-gene interactions analysis are consistent with observations of other genetic studies which demonstrated an importance of ADE genes for asthma pathogenesis. In particular, we confirmed a potential role in the pathogenesis of asthma for CYBA and CAT genes that was associated with asthma in Czech [60] and Canadian [61] populations, respectively. In addition, the associations of asthma with the IL5 C-703T polymorphism in Russians from the city of Tomsk [35] and IL1B −511C>T in the Canadian Asthma Primary Prevention Study [62] have been successfully replicated in our study. Our study is consistent with the observation that glutathione-S-transferase genes M1, T1, and P1 alone and in combination with other ADE genes do not play a substantial role in the development of BA [63]. We also found for the first time genetic polymorphisms of the GSR and PON2 genes can be important determinants of susceptibility to asthma, but their associations need to be confirmed in independent populations. Further studies should also be focused on the analysis of gene-gene interactions to better understanding the role of ADE genes in asthma pathogenesis. In the study of Millstein et al. [16], interactions between the NQO1, MPO, and CAT genes have been identified in ethnically diverse cohorts of patients with childhood asthma, whereas marker-by-marker analysis did not reveal the associations of these genes with disease susceptibility. This means that marker-by-marker approach ignores the multigenic nature of BA and does not evaluate a complexity of interactions between susceptibility genes.

Comparing the results obtained by the three statistical approaches to the analysis of gene-gene interactions, we can say that, despite gender-specific effects of genotypes on the pathogenetic variants of the disease, each of the methods showed own uniqueness and efficacy in the detecting genes associated with asthma risk. In our point of view, the advantage of SAA method is its capacity in the identification of “gene dosage effects” of different sets of ADE genes on asthmatic phenotype. Meanwhile, MDR method, especially its cluster technique, was found to be powerful in the detecting high-order epistatic interactions between ADE genes and their synergic and antagonistic effects on the asthma risk. The variability in the structure of gene-gene interactions models across the pathogenetic variants of asthma can be partially explained by differences in bioinformatic approaches to the analysis of multiple genes. Apparently, a similarity in gene-gene interactions between the models obtained by the two different bioinformatic tools may be explained by strong effects of particular genes on the asthmatic phenotype. This means that strong phenotypic effects of the CAT, GPX1, GSR, GCLM, EPHX1, CYBA, and IL5 genes (they showed the similarity between the models) may be considered as major gene effects. And, finally, a post hoc comparative analysis of the frequencies of genotype combinations was useful in the detection of ADE genes with low or moderate effects on asthma as well as the unique combinations of ADE genotypes strongly associated with disease susceptibility. In particular, the PON2, GPX2, GPX3, GPX4, NQO1, FMO3, SOD2, and IL3 genes showed low or moderate effects on asthma risk and may represent a polygenic background of the disease susceptibility. Thus, the methods complemented each other and contributed to the understanding of the polygenic nature of asthma and complexity of gene-gene interactions underlying the asthmatic phenotypes. Due to the fact that there is no universal method for comprehensive analysis of gene-gene interactions in genomic epidemiology, it makes sense to use several methods, as it has been successfully applied in the present study.

The dendrograms obtained by MDR technique (Figure 3) clearly showed complex and hierarchic pattern of interactions between ADE genes constituting the polygenic basis of the pathogenetic variants of asthma. In particular, the GSR T>C and IL5 C-703T genes in men and the CYBA 640A>G, PON2 C311S, and GPX4 C718T loci in women had the highest degree of synergy in their interactions to determine the susceptibility to allergic asthma, whereas the highest degree of synergy in gene-gene interactions in nonallergic asthma in men was found for the CYBA −930A>G, FMO3 E158K, and GPX1 P198L loci. In contrast, a different degree of redundancy (antagonism) in gene-gene interactions was observed between the CAT −21A>T and GPX2 G>A loci in men and between the EPHX1 Y113H and IL5 loci C-703T in women with allergic asthma, as well as between the GPX1 P198L, GPX2 G>A, EPHX1 Y113H, and IL5 C-703T loci in women with nonallergic asthma. Interestingly, the EPHX1 Y113H and IL5 C-703T genes showed an antagonistic character of gene-gene interactions exclusively in women with both pathogenetic variants of asthma. Notably, the EPHX1 Y113H genotypes did show the association with asthma risk in single-locus analysis performed in our previous study ( ) [19]. A strong synergism or antagonism in the interaction between the ADE genes in determining different types of asthma may suggest that the gene-gene effect can be driven by their true interaction, rather than by the main effect from the distinct gene. These findings may indicate epistatic interactions of the ADE genes, a situation when the effect of one gene may not be disclosed if the effect of another gene is not considered [64].

A post hoc comparative analysis of the frequencies of genotype combinations in the study groups revealed two-locus combinations of the ADE genotypes which increase the risk of the development of asthma. We found relatively rare combinations of genotypes which gave the highest asthma risk estimates but were limited to small subgroups of subjects. In particular, frequencies of these genotype combinations varied from 1 to 9% among healthy controls and from 17 to 41% among patients with nonallergic asthma, whereas odds ratios for disease risk varied from 6.7 to 15.6. Moreover, there was an obvious excess of combinations of variant genotypes among asthma patients compared with healthy subjects, and these differences reached statistical significance after adjusting for multiple tests.

4.3. Genetic Variation in ADE Genes and Asthma Pathogenesis

Despite nonsignificant differences in the genotype distributions, the two-locus comparison of genotype frequencies between the study groups has shown that asthma patients more often than healthy subjects carry combinations of the genotypes which are known to determine a diminished activity of ADE towards ROS. This is supported by a number biochemical studies that observed massive generation of ROS and the insufficiency in antioxidant capacity in asthma [6568]. In this context, it is important to highlight that the studied ADE genes alone cannot account for the whole polygenic mechanisms underlying such biochemical abnormalities in asthma. The interactions between ADE genes that we have identified using different statistical methods make a mechanistic sense because these genes are collectively involved in the maintaining and regulation of redox homeostasis. Moreover, the integrated function of ADE genes in the lung and airways can promote a coordinated detoxification of xenobiotics-induced ROS, thus preventing oxidative stress which plays an important role in the pathogenesis of asthma [21, 24, 51, 69]. Importantly, ADE genes showed interactions with other asthma-related genes such as IL5 and IL1B which are responsible for the immunological mechanisms of asthma and allergy.

Based on the literature data demonstrating biochemical abnormalities in redox homeostasis in asthma and the results of our study, we assumed possible relationships between these abnormalities and ADE genes showed the associations with asthma in our study (the data are shown in Table 7). Changes in the activity of antioxidant defense enzymes such as glutathione peroxidases and catalase in whole blood, plasma, platelets, and bronchoalveolar lavage fluid have been reported by a number of biochemical studies, the findings which are in accordance with our results demonstrating the relationship between the genes for these enzymes and the risk of different pathogenetic variants of asthma. Briefly, an enhanced production of ROS by blood neutrophils, monocytes, and eosinophils found in asthma can be explained by the effects of functional polymorphisms in the gene encoding p22 phox subunit (CYBA) of NADPH oxidase. Genetic variation in the GSR and GCLM genes may be responsible for biochemical perturbations of glutathione metabolism such as an increased level of oxidized glutathione in asthma. Polymorphisms of the EPHX1 gene determine the increased activity of the enzyme, thus leading to the enhanced production of reactive semiquinones.

Although we did not perform biochemical investigations of antioxidant status, taking the observed association of asthma with the ADE genotypes and their functional significance into account, it is likely that an imbalance between oxidants and antioxidants detected in asthma can be directly related to genetically diminished capacity of ADE. Such an imbalance results in oxidative stress caused by an excessive production of ROS and/or by inadequate antioxidant defense leading to damage of airway epithelial cells and inflammation due to upregulation of redox-sensitive transcription factors and proinflammatory genes [22, 24]. We may also conclude that ADE genes seem to play a greater role in the development of nonallergic asthma than in allergic asthma.

4.4. Gender-Specific Effects of ADE Genes on Susceptibility to Bronchial Asthma

An important finding of our study was that polymorphisms of many ADE genes showed sex-specific associations with the development of asthma. For instance, the CAT −21A>T and GCLM −588C>T gene polymorphisms were associated with asthma susceptibility exclusively in men. In contrast, polymorphism 640A>G of the CYBA gene showed a relationship with asthma risk only in women. These findings demonstrate sexual dimorphism in genetic susceptibility to asthma, a phenomenon established for many complex human diseases [70, 71] including asthma [72]. It has been proposed that existing variation in regulatory elements of genes rather than differences in their structure in men and women may explain sex-specific genotype-phenotype interactions in complex traits [70, 73]. Sex-specific changes in age-related gene regulation can result in the difference in asthma susceptibility between the sexes [70]. We suggest that the mechanisms underlying gender-related specificity in the associations of ADE genes with BA found in our study are related to differential expression of redox-sensitive genes in men and women. Since estrogen was found to depress oxidative stress in mice [74], sex-steroid receptors might be an example of sex-specific trans-regulatory elements [75] for redox-sensitive genes which in turn may differentially respond to the inducers due to their functionally unequal polymorphic alleles. This means that ADE genes can function differently in men and women in some circumstances. This suggestion is supported by the finding of gender difference in both expression and activity of antioxidant enzymes demonstrated in animal studies [76, 77]. Therefore, it can be concluded that identification of gender-specific genetic variants of ADE, which contribute to the shift of redox homeostasis towards oxidative stress, will provide a better understanding of sex-specific regulation of ADE gene expression and differences in the molecular mechanisms of asthma in men and women.

4.5. Limitations of the Study

The study has limitations. Due to the relatively small sample sizes of the studied groups, the association analysis of two-locus genotype combinations was underpowered, especially after Bonferroni adjustment for multiple tests. Because of the limited sample size, we also cannot exclude the possibility that small effects of some ADE genes were not detected. Since BA is a multifactorial and genetically heterogeneous disease [78, 79], further studies with larger sample sizes with genotyping of more polymorphic variants of ADE genes are required for better understanding of the roles of these genes in asthma pathogenesis. Because we did not analyze expression profiles of the genes and biochemical parameters of redox homeostasis, both functional genomics and metabolomics studies are required to clarify the molecular mechanisms by which polymorphisms of ADE genes contribute to the development of BA. Since the risk of BA is determined by a complex interplay between genetic and environmental factors, further genetic studies should take into account environmental factors that may play a significant role in the etiology of the disease.

5. Conclusions

To the best of our knowledge, this is the first study investigating the associations between BA and 34 functionally significant polymorphic variants of ADE genes and 12 other candidate genes. So far, no genetic studies have reported a comprehensive evaluation of asthma susceptibility with a number of ADE genes at once. Methodological approaches used in this study were proved fruitful in uncovering the genetic architecture of complex interactions between genes involved in the regulation of redox homeostasis. This allowed finding for the first time that antioxidant defense enzymes genes are collectively involved in the molecular mechanisms of BA and can explain genetic heterogeneity between allergic and nonallergic variants of the disease. In particular, we found for the first time that the GSR and PON2 genes can be referred to as novel asthma susceptibility genes, but their associations need to be confirmed in independent populations. We also showed both complexity and diversity of gene-gene interactions in allergic and nonallergic asthma. Finally, we have discovered gender-specific effects of ADE genes for the risk of the pathogenetic variants of asthma. Altogether the study results provide strong evidence for the pathogenetic role of ADE genes in asthma. Our data on the relationship of the ADE genes and asthma are concordant with the results of a number of biochemical studies demonstrating the massive generation of ROS and the insufficiency in antioxidant capacity which have been implicated in pathogenesis of asthma.

Further studies focusing on the molecular mechanisms regulating redox homeostasis can provide more complete understanding of the role of the ADE genes in bronchial asthma and end up in the discovery of new drug targets for antioxidant treatment and prevention of the disease.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The study was supported in part by the Federal Targeted Program “Scientific and Scientific-Pedagogical Personnel of the Innovative Russia in 2009–2013.” The authors thank Drs. Mikhail Kozhuhov and Valery Panfilov from Kursk Regional Clinical Hospital for their invaluable help in recruiting and examining asthma patients for the study.

Supplementary Materials

On comparative expression profiles of 23 ADE genes in various cell types/tissues/organs in human

  1. Supplementary Materials