Association between Glucose-6-Phosphate Dehydrogenase Deficiency and Asthma.

BACKGROUND: Among the determinants contributing to the pathogenesis of asthma, antioxidant genetic factors play a leading role. Glucose-6-phosphate dehydrogenase (G6PD) is an enzyme that is competent to detoxify free radicals. Although a relationship between G6PD deficiency and asthma has been previously reported, the literature is still scanty. In this study, we test this hypothesis in a large cohort of patients from Sardinia, Italy.
METHODS: A retrospective case-control study was performed using data from 11,829 clinical records of outpatients referred to a teaching hospital for a medical visit. In total, 455 cases (asthma-positive) and 11,374 controls (asthma-negative) were compared for G6PD status using multivariable analysis, adjusting for all covariates.
RESULTS: Overall, G6PD deficiency was detected in 11.2% of study participants and was associated with an increased risk of asthma (odds ratio (OR) 1.63; 95% confidence interval (CI) 1.27-2.10). Additional variables significantly associated with asthma were female sex (OR 1.66; 95% CI 1.34-2.06), overweight/obesity (OR 1.56; 95% CI 1.27-1.92), smoking (OR 1.44; 95% CI 1.449-3.963), and high socioeconomic status (OR 1.40; 95% CI 1.16-1.70), whereas age was inversely related with asthma (OR 0.49; 95% CI 0.39-0.61).
CONCLUSIONS: Our study shows that G6PD deficiency is an independent risk for asthma. These findings suggest that G6PD should be assessed in asthmatic patients for better risk stratification.

1. Introduction

Asthma is an inflammatory non-communicable disease of the small airways that affects more than 330 million people worldwide [1]. The prevalence rate shows differences between countries and studies [2], ranging from 0.2% in China to 21.0% in Australia [3]. Although asthma is considered the most common chronic disease among children, it can affect adults as well. As with other inflammatory diseases, asthma is more prevalent in adult females compared to males, although the gender disparity is reversed in childhood, where asthma is more frequent among boys than among girls [4]. The physiopathology of asthma is characterized by the persistent inflammation of small airways with eosinophilic infiltration in the mucosal lamina propria, leading to airway hyperresponsiveness to a wide variety of exogenous and endogenous stimuli [5]. Consequently, patients suffer from wheezing and cough, sometimes undergoing a progressive lifetime disability.

Apart from mucosal eosinophils, additional immune cells participate in long-lasting disease as well. Exposure to allergens triggers the production of specific IgE antibodies, leading to the overexpression of Th2 type T-cell response. The IgE bound to high-affinity receptors on dendritic cells facilitate allergen internalization; once inside the dendritic cell, the processing of allergens by cathepsin S and the subsequent selection of peptides loaded onto and presented by HLA molecules (MHC class II) are fundamental steps for these cells to act as antigen-presenting cells to T lymphocytes [6]. Presentation of a selected antigen peptide to the T-cell receptor induces sensitization and a subsequent immune response to the specific allergen [7].

Among the risk factors involved in the development of asthma, dietary, environmental, and genetic factors that are capable of reducing cell antioxidant capacity by increasing tissue vulnerability to oxidative stress raise the disease risk. For instance, selenium deficiency lowers red cell glutathione peroxidase activity and is associated with increased asthma risk [8]. Similarly, low dietary intakes of vitamins C and E appear to increase the risk as well [9,10]. Additional predisposing factors are family history, perinatal factors such as maternal smoking [11], age, diet and vitamin D deficiency [12], presence of other atopic diseases, maternal and childhood exposure to medications such as acetaminophen [13], proton pump inhibitors [14], certain antibiotics, and exposure to bacteria and bacterial products early in life. Moreover, abdominal obesity has been suggested to raise the risk of developing asthma [15]. The major environmental risk factor is a smoking habit, which seems to alter mucosal proteins [16]. Although the pathogenesis of asthma remains elusive, several genetic factors can interact with environmental risk factors to enhance the inflammation and tissue damage caused by oxidative stress, promoting the disorder.

Interestingly, several years ago, a few studies reported an increased risk for developing asthma in subjects carrying a glucose-6-phosphate dehydrogenase (G6PD) defect [17,18]. However, for several decades, the topic was shelved, and, until now, the literature on this association is scarce.

G6PD is known as the rate-limiting enzyme of the pentose phosphate pathway [19]. The enzyme produces nicotinamide adenine dinucleotide phosphate (NADPH), the essential cofactor in oxidoreductive metabolism, maintaining a high ratio of reduced/oxidized glutathione (GSH) and acting as a substrate for NADPH oxidase (NOX) and nitric oxide synthase (NOS). G6PD deficiency is the most common enzyme defect worldwide, affecting 400–500 million people [20]. Individuals harboring enzyme deficiency are generally asymptomatic but, in certain circumstances, may show episodic anemia. Upon exposure to fava beans, infections, or medicines with a high redox potential, including non-steroidal anti-inflammatory drugs, older red cells, where the G6PD deficiency gradient is higher, undergo hemolysis. Moreover, in newborns with G6PD deficiency, decreased bilirubin excretion may result in jaundice. More recently, it has been reported that G6PD deficiency may exert a proinflammatory effect depending on the animal or human model and the tissue involved [21].

On the Sardinian island in Italy, G6PD deficiency is very common, e.g., around 10–12% of the general population [22]. On the other hand, data according to the Gender Environment Interactions in Respiratory Diseases (GEIRD) indicate a relative risk of asthma in Italy of 6.6 (95% CI 6.1–7.1) and, more specifically, for Sassari (northern Sardinia), a relative risk of 7.3 (95% CI 6.0–8.8) [23]. An additional epidemiological study in southern Italy, including Sardinia and Sicily, reported a prevalence of asthma ranging between 2.6% and 3.2%, making the island the ideal model to test the G6PD/asthma association [24].

Based on these premises, we attempted to evaluate the association between asthma and G6PD deficiency in a defined population of northern Sardinia.

2. Materials and Methods

2.1. Study Design

This was a retrospective case–control, single-center study that took advantage of the availability of the clinical records of adult outpatients referred to a teaching hospital of northern Sardinia (Department of Internal Medicine, University of Sassari), Italy, from January 2002 to December 2019. Patients were referred by their family physicians and/or specialists for any reason. Patients with asthma were considered, and the controls were patients without asthma.

2.2. Eligibility Criteria

Records of personal information such as sex, age, smoking habits, socioeconomic status, anthropometric parameters (body height and weight), all signs and symptoms, treatments, and an accurate medical history were considered eligible for the analysis. More specifically, for the purpose of the study, the presence of asthma was retrieved from the computerized database among previously diagnosed conditions. Moreover, the diagnosis of asthma retrieved from patient’s chart was double-checked by matching the treatment (for example, inhaled corticosteroid or mixed inhaled corticosteroid, long-acting β2-agonist on maintenance or on demand, or short-acting β2-agonist upon need), according to the Global Initiative for Asthma (GINA) consensus [25]. Availability of G6PD status was considered the major inclusion criteria. Data from each patient were collected using a standard form for the entire study period, and each visit was supervised by the same attending physician.

2.3. Exclusion Criteria

Incomplete records or records belonging to patients younger than 18 years were not included in the analysis. In the case of multiple visits for the same patient within the given time period, only the most recent was considered for the analysis.

2.4. Diagnostic Criteria

G6PD status. Since Sardinia is a high-prevalence region for G6PD deficiency, the subjects are routinely assessed for G6PD status, especially before exposure to certain medications. The G6PD activity of whole blood was measured using a standard method, as previously described [26,27]. Less than 10% residual activity was used to define the presence of G6PD total deficiency, while residual activity ranging from 10% to 50% was used to define partial deficiency. Molecular analysis was not available for G6PD-deficient patients.

Asthma diagnosis. Asthma is a chronic obstructive condition characterized by both the inflammatory component and airway obstruction, which may be spontaneous or induced by different triggers, different from chronic obstructive pulmonary disease (COPD), which typically shows a fixed bronchial obstruction. The measure of pulmonary function assessed by spirometry is considered the gold standard to diagnose asthma [28]. According to the European Respiratory Society and American Thoracic Society (ERS/ATS) guidelines, a ratio between forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) less than 0.7 is the cut-off used to diagnose an obstructive respiratory pattern, while the GINA consensus suggests a cut-off of 0.8 [25]. To establish if the bronchial obstruction is permanent or temporary, pharmacological reversibility testing is used. The test is performed by repeating spirometry after inhalation of a rapid onset β2-agonist such as salbutamol [25,28]. The test is considered positive if there are FEV1 increases of 200 mL or FEV1 or FVC increases of 12% from the basal value [29]. A history of wheezing, cough, and chest tightness could be suggestive of the disorder, especially when there is a positive family history of asthma.

2.5. Ethical Considerations

The protocol was approved by the local Ethics Committee (Comitato di Bioetica, Azienda Ospedaliero-Universitaria di Sassari, Italy) (Protocol no. 3004/CE, 2016).

2.6. Statistical Analysis

Distributions of age, sex, socioeconomic status (SES), body mass index (BMI), and smoking habit were descriptively compared between cases and controls. More specifically, age was recoded into a binary variable by splitting at the age of 60 years. Socioeconomic status was estimated using current or past occupation and divided into four categories, from class I (the highest) to class IV (the lowest), as previously reported [30], and expressed as a binary variable: (i) high SES (classes I and II) and (ii) low SES (classes III and IV). BMI was calculated using the formula of weight (kg)/height (m2), and overweight/obesity was defined as BMI > 25 kg/m2. In relation to smoking habits, patients were stratified as never or current/former smokers. Patients with total and partial G6PD deficiency were pooled together.

Univariate and multivariable logistic regression was used to examine the association of asthma with G6PD status by calculating odds ratios (ORs) and their 95% confidence intervals (CI) using the Wald formula: 95% CI = OR1±(β/SE). The main effect of variables was assessed first by entering the covariates one by one and then simultaneously (the adjusted model). In addition, patients with asthma were stratified according to the GINA classification, based on asthma symptoms and therapy, in order to ascertain the frequency of G6PD deficiency in each subgroup (intermittent, mild persistent, moderate persistent, and severe persistent asthma).

All statistical analyses were carried out using SPSS statistical software (version 22.0, Chicago, IL, USA), and two-sided p values lower than 0.05 were considered statistically significant.

3. Results

A total of 11,829 clinical records of patients (7323 females, 61.9%) who underwent a medical visit were available for the analysis. Descriptive statistics of the studied population are reported in Table 1.

Table 1

Characteristics of 11,829 study participants with and without asthma.

Covariates	Asthma (Cases)	No Asthma (Controls)	p-Value
Age, n (%)
<60	322 (4.8)	6352 (95.2)
≥60	133 (2.6)	5022 (97.4)	<0.0001
Sex, n (%)
Female	315 (4.3)	7008 (95.7)
Male	140 (3.1)	4366 (96.9)	0.001
SES 1, n (%)
Low	230 (3.3)	6807 (96.7)
High	225 (4.7)	4567 (95.3)	<0.0001
BMI 2, kg/m2
<25	211 (3.6)	5652 (96.4)
≥25	244 (4.1)	5722 (95.9)	0.166
Smoke
No	316 (3.5)	8602 (96.5)
Yes	139 (4.8)	2772 (95.2)	0.003
G6PD 3 status
Normal	378 (3.6)	10,132 (96.4)
Deficiency	77 (5.8)	1242 (94.2)	<0.0001

1 Socioeconomic status; 2 body mass index; 3 glucose-6-phosphate dehydrogenase.

The prevalence of asthma among study participants was 3.8%, according to a previous report [24]. Asthmatic patients were significantly younger than non-asthmatics (asthma: 49.1 ± 17.1 vs. non-asthmatics: 55.2 ± 18.0 years; p < 0.0001). As expected, the prevalence of asthma was greater among females compared with males (4.3 vs. 3.1, p = 0.001).

Individuals with a low SES were fewer in the asthmatic group (3.3% vs. 4.7%, p < 0.0001). Overweight/obesity (≥25 kg/m2) was more common among asthmatics, although this difference did not reach statistical significance (Table 1). The proportion of current or former smokers was significantly greater among asthmatics (4.8% vs. 3.5%, p = 0.003). Based on the laboratory data retrieved from the participants’ clinical records, 1319 patients (11.2%) were partially or totally deficient in G6PD, with greater frequency among females (13.4%) than among males (7.7%), in line with the sex-linked trait. Interestingly, the proportion of asthma sufferers was higher among G6PD-deficient patients than among patients with normal enzyme activity (5.8% vs. 3.6%, p < 0.0001) in both males (6.4%) and females (4.3%).

The proportion of asthma decreased in patients in relation to age, as shown in Figure 1.

Figure 1

Asthma prevalence among the study participants, stratified by age groups.

In Table 2, the risk of asthma occurrence is listed according to sex and age in G6PD-normal and -deficient patients.

Table 2

Risk of asthma according to sex and age in G6PD-normal and -deficient individuals.

	Cases (Asthma)		Controls (No Asthma)		OR § (95% CI †)
	G6PD #Normal	G6PDDeficient	G6PDNormal	G6PDDeficient
Sex, n (%)
Female	253	62 (19.6)	6097	911 (12.9)	1.64 (1.23–2.18) **
Male	124	16 (11.4)	4035	331 (7.6)	1.57 (0.92–2.68)
Age, n (%)
<60	271	51 (15.8)	5629	723 (11.4)	1.47 (1.08–2.00) *
≥60	106	27 (20.3)	4503	519 (10.3)	2.21 (1.43–3.40) **
Total patients	377	78 (17.1)	10,132	1242 (10.9)	1.69 (1.31–2.17) **

§ OR = odds ratio; † CI = confidence interval; # G6PD = glucose-6-phosphate dehydrogenase; * p < 0.05, ** p < 0.01.

The unadjusted risk of asthma was significantly higher than the unity among G6PD-deficient patients (1.69, 95% CI 1.31–2.17).

The increased risk of asthma among G6PD-deficient subjects was observed in both sexes, although statistical significance was reached in only the female subgroup. More specifically, an increasing trend was observed across the various age groups; the trend reached significance in the 60–79 years group (ANOVA, p < 0.0001) (Figure 2).

Figure 2

Asthma risk in G6PD-deficient patients compared to non-deficient patients according to age decades.

Remarkably, although the absolute risk of asthma was inversely related to age, in subjects with G6PD deficiency compared to normal subjects, the risk of asthma increased in an age-dependent manner (Figure 2) and reached a maximum after age 60 (OR 2.21, 95% CI 1.43–3.40) (Table 2).

Table 3 reports the univariate and multivariate logistic regression analysis for the association of G6PD status with asthma, adjusting for covariates. G6PD deficiency significantly increased the risk in the unadjusted model (1.69, 95% CI 1.31–2.17), remaining significant after adjusting for sex, age, SES, BMI, and smoking (1.63, 95% CI 1.27–2.10).

Table 3

Logistic regression analysis for G6PD status and other variables potentially associated with the risk of asthma.

Covariates	UnadjustedOR ‡s and 95% CI	AdjustedORs and 95% CI
G6PD # status
Normal	Ref	Ref
Deficiency	1.69 (1.31–2.17) **	1.63 (1.27–2.10) **
Age, yrs
<60	Ref	Ref
≥60	0.52 (0.42–0.64) **	0.49 (0.39–0.61) **
Sex
Male	Ref	Ref
Female	1.40 (1.14–1.72) **	1.66 (1.34–2.06) **
SES §
Low	Ref.	Ref.
High	1.46 (1.21–1.76) **	1.40 (1.16–1.70) **
BMI, kg/m2
<30	Ref.	Ref.
≥30	1.14 (0.94–1.38)	1.56 (1.27–1.92) **
Smoking
No	Ref.	Ref.
Yes	1.36 (1.11–1.67) **	1.44 (1.17–1.77) **

# G6PD = glucose-6-phosphate dehydrogenase, § SES = socioeconomic status, ‡ OR = odds ratio; ** p < 0.01.

Table 4 reports the G6PD status according to the severity of the asthma phenotype. Compared with study participants without asthma, there was a trend to higher frequencies of G6PD deficiency among participants with the highest severity asthma levels, which was statistically significant for moderate persistent asthma.

Table 4

The severity of asthma based on the Global Initiative for Asthma (GINA) classification according to glucose-6-phosphate dehydrogenase (G6PD) status.

Severity of Asthma According to GINA # Guidelines	Drugs Used to Treat Asthma	G6PD-NormalNo. (%)	G6PD-DeficientNo. (%)
Intermittent asthma	Low dose ICS §-formoterol as needed, with rapid onset LABA ¶ as needed, or low-dose ICS whenever SABA $ used	206 (85.5)	35 (14.5)
Mild persistent asthma	Daily low dose ICS with SABA or low dose ICS-formoterol as needed, or low-dose ICS plus SABA $ concomitantly as needed, or LTRA‡ daily and SABA $ as needed	74 (84.1)	14 (15.9)
Moderate persistent asthma	Low dose ICS-LABA as maintenance and reliever therapy, or low-dose ICS plus LTRA daily, +/− SABA as needed	84 (76.4)	26 (23.6) **
Severe persistent asthma	Medium or high-dose ICS-LABA daily and SABA as needed, or high dose of ICS plus tiotropium, or LTRA +/− short course of oral glucocorticoids +/− add-on therapy (e.g., tiotropium, zileuton, anti-IgE, anti-IL-5, anti-IL-5R, anti-IL-4R), or oral glucocorticoids, or addiction of biologics	14 (87.5)	2 (12.5)

# GINA, Global Initiative for Asthma; § ICS, inhaled corticosteroid; ¶ LABA, long-acting beta agonists; $ SABA, short-acting beta agonists ‡ LTRA, leukotriene receptor antagonists; ** p < 0.01.

Since steroids can inhibit G6PD activity, at least in animal models [31], there is a theoretical possibility that asthmatic patients may have an enzyme deficiency due to therapy. However, as can be seen in Table 4, the frequency of G6PD deficiency in patients with the most severe form of asthma, exposed to systemic steroids, was similar to that of the normal population.

4. Discussion

In the present study, conducted in a cohort of subjects referred to the Department of Medicine of Sassari, Italy, a significant association was found between inherited G6PD deficiency and the risk of asthma. This association, first reported over three decades ago [17,18], for many years, did not attract enough interest from investigators for two probable main reasons. First, the populations with the highest prevalence of G6PD deficiency live mostly in low-resource countries and, therefore, are less exposed to stressful environmental factors, such as pollution [32]. A second reason may be that until very recently, the enzyme was considered crucial only for red blood cells, which depend entirely on PPP for the maintenance of reducing capacity, while other cell types possess alternative sources of NADPH [33]. In reality, the expression of the enzyme is high in leukocytes and other immune cells [34], and, recently, G6PD deficiency has been found to be involved in several disorders beyond blood diseases, including sepsis [35], neonatal hyperbilirubinemia [36], and cardiovascular disease [37], especially in the elderly [38].

Among Sardinians, there is a prevalence of about 10–12% [22], essentially determined by a founder mutation (G6PD Med, S188F) in more than 95% of deficient subjects. The residual activity of G6PD Med is less than 10% in hemizygous males and homozygous females and between 10% and 50% in heterozygous females. Therefore, this population is an ideal model to test the hypothesis of the G6PD/asthma association with an adequate degree of statistical power.

The results obtained in this study provide compelling evidence that subjects with inherited G6PD deficiency are prone to develop asthma more frequently than subjects with normal enzyme activity, regardless of sex, SES, smoking habit, and excess weight. The increased risk was observed especially in older subjects, reaching statistical significance in the 60–79 years age group. The association probably persists in later ages as well, but the small number of older asthmatics in the examined sample cohort did not allow this hypothesis to be tested with sufficient statistical power. This finding suggests that G6PD-deficient individuals may experience progressively impaired antioxidant defense, putting them at increased risk to develop asthma in the presence of predisposing triggers, including aging. The same association may exist in children, although, in this study, it was not evaluated.

The main pathogenetic mechanism of asthma is due to chronic inflammation, largely driven by the innate immune response. Reactive oxygen species (ROS) are powerful mediators of the inflammatory response, and the increased production or decreased neutralization level of ROS plays an important role in the pathogenesis of asthma.

In a study aimed at identifying differentially expressed proteins in the serum of children with and without asthma, G6PD was found to be significantly downregulated in asthmatic children [39], suggesting that enzyme deficiency may increase ROS, triggering kinase pathways that are able to facilitate viral replication and, in turn, aggravating airways injury. An additional study performed on a pediatric cohort using the microarray profiling technique identified the G6PD gene among the first three differentially expressed genes out of the 274 genes associated with asthma [40]. G6PD deficiency could increase the risk of infections, including those of the respiratory tract [41], that are implicated in the etiology of asthma. It is known that in vitro cells with G6PD knockdown are more susceptible to viral infections [42], and during the recent COVID-19 pandemic, G6PD deficiency was seen as an aggravating factor in the clinical picture and prognosis [37]. It can, therefore, be hypothesized that children with G6PD deficiency experience viral and/or bacterial infections more frequently, paving the road for subsequent hyperreactivity.

In addition to infections, a key role in the pathogenesis of asthma is played by oxidative stress; this is supported by experimental and clinical studies, for example, those focused on exposure to strong oxidants such as ozone [43,44]. G6PD activity increases following intravenous administration of ozone in humans [45,46] as well as in alveolar macrophages in mouse models exposed to ozone [47]. Interestingly, Varghese et al. demonstrated that enzyme deficiency influences metabolic fluxes and pulmonary hypertension in a mouse line [48].

A major function of G6PD is to regenerate the reduced thiol GSH consumed during the neutralization of ROS [49]. In experimental asthma models, the GSH is significantly decreased compared to controls, supporting the pivotal role of intracellular antioxidative mechanisms in this disorder [50]. Additionally, in human asthma, the level of reduced GSH is significantly increased, indicating a repairing function [51]. Thus, it is reasonable to speculate that G6PD-deficient subjects may have insufficient production of NADPH to maintain adequate GSH stores in asthma patients (Figure 3).

Figure 3

The putative role of G6PD deficiency in the pathogenesis of asthma. Nitric oxide (NO) production relaxes the smooth muscle of airways and vessels, inducing bronchodilation and vasodilation. G6PD deficiency induces bronchoconstriction by lowering NO generation by nitric oxide synthase (NOS). In addition, G6PD deficiency causes reduced glutathione (GSH) depletion, thus decreasing antioxidant defense and increasing the generation of free radicals.

Furthermore, G6PD deficiency causes the depletion of nitric oxide (NO) [52]. This compound, being both a major signaling molecule and a free radical, has a dual effect: at low levels, it has mainly an airway-muscle-relaxing role, and, at high levels, it activates proinflammatory mechanisms; therefore, in asthma, the effect of NO depletion is complex, and, in severe G6PD, a deficiency may be beneficial as well as adverse. It can be conjectured that the chronic depletion of NO, resulting from G6PD deficiency, may affect the basal bronchodilator tone sustained by this molecule [53]. More importantly, NO depletion manifests its effect, especially in the elderly, where there is a progressive impairment of antioxidant mechanisms [54]; perhaps this may partly explain why in our study, the magnitude of the association of G6PD deficiency with asthma was higher in patients after age 60 (Figure 3).

Finally, G6PD deficiency is associated with a proinflammatory state or the exacerbation of a proinflammatory state [21]. Recently, an increasing number of studies have reported that G6PD deficiency alters specific cytokine pathways involved in several disorders. For example, an in vitro model based on human hepatocarcinoma HepG2 cells, in which G6PD was inactivated by siRNA, revealed an increased production of neutrophil chemoattractant interleukin-8 (IL-8), which was mirrored by an enhanced production of reactive oxygen species [55]. At the molecular level, the mechanism has been investigated in detail and involves the activation of the NF-κB signaling pathway. Increased serum levels of IL-8, as well IL-4, IL-5, and IL-9 have been reported in asthma hyperreactivity following Th2 lymphocyte release [56]. Therefore, G6PD-deficient subjects may experience over-expression of IL-8 and increased chemotaxis of eosinophils and other inflammatory cells, which are found in abundance in sputum or bronchoalveolar lavage fluid. Moreover, since TGF-β secretion was found to be increased in G6PD-deficient macrophages [57], this can be an aggravating factor of the inflammatory state (Figure 3).

The association of asthma with the G6PD gene, which maps on the X chromosome (OMIM 305900), could theoretically be due to a linkage with a predisposing gene located near the G6PD locus. The most interesting candidate is probably the cysteinyl leukotriene receptor 1 (CYSLTR1) gene, located on chromosome Xq21.1 [58], which has been implicated in several allergic disorders, including asthma [59]. However, the relative distance between the two loci does not totally exclude the possibility of recombination events, making this hypothesis weak and implicitly strengthening a direct pathogenetic role of the G6PD gene itself.

Some limitations of this study need to be mentioned. First, being a retrospective study, we cannot exclude that potential confounders may be differently distributed between cases and controls. However, given the large number of participants, we are confident that the extent of such a bias would be minimized. An additional limitation of our study may be the lack of molecular genotyping for patients with biochemically identified G6PD deficiency, precluding the possibility, in the case of females, of distinguishing homozygotes and heterozygotes and, therefore, testing the association of asthma with total or partial G6PD deficiency separately. However, from a clinical standpoint, the lack of a more detailed analysis would not have changed the overall findings.

5. Conclusions

In conclusion, in our study, we found that G6PD deficiency was associated with an increased risk of asthma, and this risk remained significant after adjusting for well-known traditional variables such as female gender, overweight/obesity, and high SES. These findings suggest that G6PD should be assessed in patients with asthma for better risk stratification.