The CRTh2 polymorphism rs533116 G > A associates with asthma severity in older females

Background: CRTh2 is G protein coupled receptor for prostaglandin D2 (PGD)₂ expressed by immune cells that drive type 2 inflammation such as CD4⁺ T cells (Th2), eosinophils and group 2 innate lymphoid cells (ILC2) as well as structural cells including smooth muscle and epithelium. CRTh2-expressing cells are increased in the blood and airways of asthmatics and severe asthma is characterized by increased activity of the PGD₂-CRTh2 pathway. The CRTh2 single nucleotide polymorphism (SNP) rs533116 G > A is associated with development of asthma and increased Th2 cell differentiation.

Objective: To examine whether CRTh2 rs533116G > A associates with asthma severity. Since severe asthma is more common in females than males, we performed a sex-stratified analysis.

Methods: Clinical data from asthmatics (n = 170) were obtained from clinic visits and chart review. Asthma severity was assessed according to ERS/ATS guidelines. Peripheral blood cells were characterized by flow cytometry and qRT-PCR. Genotyping was performed by TaqMan assay.

Results: Older females (≥45 years) homozygous for minor A allele of rs533116 were more likely to have severe asthma, lower FEV₁, a higher prescribed dose of inhaled corticosteroid and more type 2 inflammation than females carrying GA or GG genotypes. Comparing females and males with the AA genotype also revealed that women had more type 2 inflammation.

Conclusions and significance: The polymorphism CRTh2 rs533116 G > A associates with severe asthma and type 2 inflammation in older females. This study reveals a gene-sex-aging interaction influencing the effect of CRTh2 on asthma severity.

Introduction

Asthma is a heterogeneous disease with many endotypes. Type 2 high asthma is the most common form and is mediated by the cytokines IL-4, IL-5, and IL-13 (1). These cytokines are produced primarily by Th2 (T-helper 2) cells and ILC2s (group 2 innate lymphoid cells) and together promote type 2 inflammation and infiltration of other important inflammatory cells such as eosinophils and basophils (2–4). Th2 cells (5–7) and ILC2 (4) are defined by the expression of CRTh2 (chemoattractant homologous receptor expressed by Th2 cells). Other inflammatory cells including eosinophils and basophils and structural cells such as smooth muscle and epithelium also express CRTh2 (8–10).

While both Th2 cells and ILC2 produce type 2 cytokines, these cells are considered to play inherently different roles in immune responses. Th2 cells are differentiated and activated by allergen-primed dendritic cells (11) and are memory T cells that circulate between the lymph nodes and periphery providing surveillance and allergen specific responses (12, 13). ILC2 can be found in the blood, but at 1,000-fold lower numbers than Th2 cells (14); this is considered to be due to their role as innate cells residing primarily at mucosal sites. ILC2 are activated by IL-25 and IL-33, cytokines released from the airway epithelium following exposure to allergens (15) or other environmental factors such as cigarette smoke or viral and/or bacterial infections [reviewed in (16)]. Both Th2 cells and ILC2 have been shown to be increased in the blood and airways of severe compared to mild/moderate asthmatics (14, 17, 18). CRTh2 is a G protein coupled receptor for prostaglandin D₂ (PGD)₂ (19, 20) and activation of CRTh2 mediates chemotaxis (19) as well as type 2 cytokine expression (21). Since PGD₂ is produced by both allergen-activated mast cells (22) and microbial-activated macrophages (23), the PGD₂-CRTh2 pathway is considered to play a role in contributing to both day to day asthma symptoms as well as exacerbations [reviewed in (24, 25)].

Genetic variations in CRTh2 have been associated with development of asthma and other allergic phenotypes (26, 27). We reported that the minor allele of rs533116 G > A, within a CRTh2 enhancer region, was associated with higher levels of CRTh2 expression on circulating CD4⁺ T cells and eosinophils and in vitro differentiated Th2 cells (28). We recently observed that the proportion of CD4⁺CRTh2⁺ T cells (Th2 cells) are increased in severe asthma (17). Others have shown more CRTh2-expressing cells and higher PGD₂ levels in bronchial alveolar lavage and epithelial brushings from severe asthmatics (18). To date, no studies have investigated whether this polymorphism is associated with clinical characteristics of asthma including severity of the disease.

In this study, we investigated whether the genotype status of rs533116 is associated with severe asthma and/or indices of asthma severity. In light of our emerging understanding of the importance of sex differences in health and disease, coupled with reports that females are more likely to suffer from severe asthma and to have more severe symptoms (29), we additionally performed a sex-stratified analysis. We found that the homogeneous minor allele of rs533116 (AA genotype) associated with asthma severity and degree of type 2 inflammation exclusively in females. Collectively, our results suggest that genetic variation in CRTh2 may play a sex-specific role in asthma severity.

Materials and methods

Subjects

The institutional ethics review boards of the University of Alberta and Western University approved this study. Patients were recruited and consented from the tertiary care Asthma Clinics at the University of Alberta, Edmonton Alberta and London Health Sciences Centre/St. Joseph’s Health Care, London, Ontario. Inclusion criteria were age >18 years and a physician diagnosis of asthma. Severe asthma was defined as asthma that requires treatment with high dose inhaled corticosteroids plus a second controller (and/or systemic corticosteroids) to prevent it from becoming “uncontrolled” or that remains “uncontrolled” despite this therapy. Inadequate symptom control was defined by any one of the following: (i) airflow limitation (FEV₁ < 80% predicted after bronchodilator medication withheld); (ii) Asthma Control Questionnaire (ACQ) > 1.5; (iii) frequent severe exacerbations (two or more systemic corticosteroids bursts (>3 days each); (iv) serious exacerbation (at least one hospitalization or ICU stay in the previous year) (30).

Genotyping

Using the Wizard^® Genomic DNA Purification Kit (Promega, Madison, WI, USA), DNA was isolated from peripheral blood mononuclear cells or whole blood following manufacturer instructions. CRTh2 rs533116 [g. −6,391 bp (G > A), upstream of the translation start site] genotyping was performed using TaqMan^® allelic discrimination assay (Applied Biosystems, Foster City, CA, USA).

Quantitative real time polymerase chain reaction

Whole blood (2 ml) was collected in PAX gene tubes and total RNA was isolated using the PAXgene Blood RNA Kit (PreAnalytiX, Qiagen, BD, Mississauga, ON, Canada) according to the manufacturer’s instruction. Reverse transcription (RT) reactions were performed using 1 μg of RNA, oligo-dT primers and Superscript II Reverse Transcriptase (Invitrogen, Burlington, ON, Canada). TaqMan gene expression assays for CRTh2 (Hs00173717_ml) and GATA3 (Hs00231122_ml) were quantified using pre-designed Taqman assays (Applied Biosystems, Carlsbad, CA, USA), 19 μl of TaqMan gene expression master mix and 1 μl of cDNA. GAPDH mRNA was used as an internal control (housekeeping control) and was quantified by a custom ⁶FAM-labeled TAMRA probe (5′-AAA TCC CAT CAC CAT CTT CCA GGA GCG A-3′) (Applied Biosystems) and the following primers GAPDH-F (5′-CAA GGCT GAG AAC GGG AAG-3′) and GAPDH-R (5′-GCA AAT GAG CCC CAG CCT T-3′). The PCR protocol consisted of 10 min at 95°C followed by 40 cycles of 30 s at 95°C and 60 s at 60°C. All samples were run in triplicate. Differences in Ct values of the gene of interest and the house keeping gene GAPDH were used to calculate delta Ct (ΔCt). Relative fold changes (RFC) were then calculated using the 2^–ΔΔCt algorithm. All ΔCt were subtracted from one subjects’ ΔCt (lowest Ct value) to calculate ΔΔCt and fold increase then calculated using the ΔΔCt as a negative exponent to the base of 2 (2^–ΔΔCt).

Flow cytometry

Determining the proportion of circulating cells expressing CRTh2 was performed as we described previously (17, 31). Briefly, whole peripheral blood was stained using antibodies against CD4, CCR3, CRTh2, or isotype control and fixed with paraformaldehyde (2%). Flow cytometry analysis was performed using a BD LSR II Flow Cytometer, with gates set in accordance with the profiles of the isotype control and/or negative control beads. The proportion of Th2 cells was determine as the proportion of CD4⁺CRTh2⁺ cells/peripheral white blood cells (pWBC). CRTh2⁺eosinophils were identified as the proportion of high side scatter (SSc^high) cells expressing CCR3 and CRTh2/WBC. Results were analyzed using FlowJo^® (TreeStar, Ashland OR, USA).

Statistical analyses

Difference in the mean value of phenotypic clinical characteristics was compared using the Mann-Whitney U-test or independent t-test (based on normality test) for continuous variables and Fisher exact t-test for categorical variables. Significant departures of genotype frequency from Hardy-Weinberg equilibrium were determined by Fisher exact t-test. For associations of the CRTh2 rs533116 polymorphism (dominant and recessive model) with severe asthma, odds ratios were calculated by binary logistic regression adjusted for sex, age and/or BMI according to analysis and p-values determined by Chi-Square. Statistical analyses were performed using SPSS (version 28, Chicago, IL, USA). Statistical significance was set at P < 0.05.

Results

Asthma population

Table 1 represents the clinical and demographic data of the population with physician diagnosed asthmatics (n = 170) recruited from specialty asthma clinics, separated into patients with severe or mild/moderate asthma determined by ERS/ATS guidelines (30). There were no differences in age, BMI or sex distribution between the two groups. Severe asthmatics had lower forced expiratory volume in one second (FEV_1;% predicted), ratio of FEV₁/forced vital capacity (FVC) and were taking higher doses of inhaled corticosteroids (ICS). Those in the severe asthma group were more likely to have FEV₁ < 80% predicted and to have needed oral corticosteroids, indicating poorer symptom control (Table 1). A sex-stratified analysis showed that males had lower FEV₁ and FEV₁/FVC than females and in older subjects (>45 years) that males were more likely to require oral corticosteroids (39%) than females (17%, p = 0.014; Supplementary Table 1).

TABLE 1

Table 1. Clinical characteristics of study subjects.

CRTh2 rs533116 G > A is associated with severe asthma in older females

We previously reported CRTh2 rs533116 G > A to be associated with susceptibility to develop asthma (28, 32). In this study, after confirming that genotype frequencies within this asthma cohort exhibited Hardy Weinberg equilibrium (Supplementary Table 2), we assessed the relationship between this polymorphism and diagnosis of severe asthma. Examining the entire population, we found no significant associations with severe asthma, although those homozygous for minor A allele (AA genotype) showed a trend for significance in females (OR = 1.48 [95% CI = 0.86–2.52]), but not males (OR = 0.87 [95% CI = 0.46–1.63]; Supplementary Table 3). Recently, the probability of having severe asthma was shown to increase with each year of life until age 45 (33). We therefore restricted our further analyses to older asthmatics. In females 45 years or older, the AA genotype was significantly associated with having severe asthma (OR = 2.50 [95% CI = 1.26–4.98], p = 0.009), while no association was observed in males (OR = 0.84 [95% CI = 0.38–1.86]; p = 0.660; Supplementary Table 4 and Figure 1).

FIGURE 1

Figure 1. The CRTh2 rs533116 G > A polymorphism increases risk of severe asthma in older females. Asthmatics with mild/moderate (n = 109) and severe asthma (n = 61) were genotyped and odds ratio (OR) and 95% CI were determined. Females 45 years and older carrying the AA genotype were at increased risk of having severe asthma compared to those of GG/GA genotype. Analysis using binary logistic regression were adjusted for age, sex and BMI.

What clinical parameters of severity are influenced by CRTh2 rs533116 G > A?

Since this polymorphism was associated with severe asthma, we next assessed its relationship to clinical outcomes of asthma severity. Examining ICS usage, we found similarly to severe asthma, that females 45 years and older homozygous for the A allele (AA genotype) were significantly more likely to be prescribed more than 1,600 μg/day of ICS (budesonide equivalent) than those of GG/GA genotype (Figure 2A). In contrast, this association was not observed in males (Figure 2B).

FIGURE 2

Figure 2. The CRTh2 rs533116 G > A polymorphism is associated with ICS usage. Asthmatics were genotyped and assessed for the likelihood of being prescribed high dose ICS (>1600 μg/day, budesonide equivalent), an indicator of asthma severity. (A) Older females (AA = 11, GG/GA = 68) and (B) older males (AA = 8, GG/GA = 31) were assessed for the influence of the A allele on this outcome. Statistical differences determined Mann-Whitney U-test.

Another critical determinant of asthma severity is lung function. Unlike analysis for severe asthma and ICS usage, we found no evidence that homozygosity for A allele (AA genotype) was associated with having lower FEV₁ compared to those with GG/GA. However, examining those with GG and AA genotypes showed significantly lower FEV₁ in the combined population 45 years and older (p = 0.005) as well as older females (p = 0.023), though not males. As such, we examined this polymorphism under the dominant model and found that older asthmatics with GA or AA genotype had FEV₁ levels lower than GG in the combined group (p = 0.001), females (p = 0.004) and males (p = 0.036; Supplementary Table 5). We also examined association of this polymorphism with symptom control, assessing its association with having FEV₁ below the normal range. This analysis revealed, with the dominant and genotypic models, that older females (Figure 3A) and older males (Figure 3B) carrying the A allele were more likely to have FEV₁ < 80% than those with the GG. As such, these results suggest both GA and AA genotypes influence FEV₁. Other measures of symptom control were only available in subjects recruited from the University of Alberta (n = 59, Supplementary Table 6) and showed no significant differences across genotypes. Although there was a promising trend for CRTh2 rs533116 AA to coincide with increased likelihood of exacerbation in females (Supplementary Figure 1), suggesting this should be studied further. Only a few subjects (n = 5) needed > 3 cycles of oral corticosteroid/year. Collectively, the nature by which CRTh2 rs533116 G > A associates with asthma severity appears to depend on biological sex and age.

FIGURE 3

Figure 3. The CRTh2 rs533116 G > A polymorphism is associated with lung function. Asthmatics were genotyped and assessed for the likelihood of having FEV₁ below the normal range of 80%, an indicator of asthma severity. (A) Older females (AA = 11, GA = 41, GG = 27) and (B) older males (AA = 8, GA = 14, GG = 17) were assessed for the influence of A allele on this outcome. Statistical differences determined Mann-Whitney U-test.

Does CRTh2 rs533116G > A influence type 2 inflammation?

To assess whether CRTh2 rs533116G > A influences the degree of type 2 inflammation, the subset recruited from the University of Alberta were extensively characterized (n = 59; Supplementary Table 6). This analysis showed that females of AA genotype had more CD4⁺ T cells, eosinophils, GATA3 and CRTh2 mRNA than those carrying GA/GG genotype, but only a trend for elevated levels of circulating Th2 cells. No significant differences in type 2 inflammation were observed across genotypes in males (Table 2). These findings were similar even when those on oral corticosteroids were removed from the analysis (females, n = 1; males, n = 4; data not shown). Since a gradual increase in type 2 inflammation with age has been reported (34), we specifically examined older asthmatics. This analysis revealed that in older females carrying the AA genotype coincided with having significantly more Th2 cells (Figure 4A) compared to those females with GG or GA genotypes. Moreover, in older asthmatics the level of GATA3 mRNA (Figure 4B) and proportion of Th2 cells (Figure 4C) were higher in females than males of AA genotype. These results suggest this polymorphism enhances type 2 inflammation in females, but not males, and that the effect is more evident with age.

TABLE 2

Table 2. Association of CRTh2 rs533116 G > A with type 2 inflammation in all asthmatics.

FIGURE 4

Figure 4. The CRTh2 rs533116 G > A polymorphism is associated with increased type 2 inflammation. A subset of asthmatics were characterized for type 2 inflammation. (A) Influence of the AA genotype on the proportion of Th2 cells (CD4⁺CRTh2⁺/peripheral white blood cells) determined by flow cytometry in older females (45 years and older; AA = 3, GG/GA = 13). (B) Level of GATA3 mRNA from whole blood and (C) proportion of Th2 cells in older females (n = 3) and males (n = 3) carrying the AA genotype. Statistical differences determined Mann-Whitney U-test.

Discussion

Development of chronic inflammatory diseases involve gene-by-environment interactions (35). Although many studies have focused on external environmental exposures, such as air pollution (36), the internal environment generated by biological sex is perhaps the most fundamental and complex of interacting factors. Males and females have different gene complements that shape the internal environment and this sex-specific genetic architecture contributes to human disease (37). Age is another factor potentially influencing the impact of genetic variation, with physiologic changes occurring throughout the lifespan (38). The prevalence of asthma is known to differ with age and sex. In childhood, higher prevalence is observed in males than in females, while asthma risk increases in girls after puberty (39). This shift is attributed to internal changes in hormonal milieu (40). By adulthood, females are twice as likely as males to have a severe asthma diagnosis and severe asthma symptoms (29, 41).

Despite the plethora of studies examining genetics of asthma and the progress toward personalized medicine (42), sex-specific analyses are often not considered. Mersha et al. showed the importance of stratifying by sex as 55% of the genetic variants identified in sex-specific analyses were not found in the combined analysis and effect sizes were often larger (43). Though some sex-specific risk alleles for asthma (44), asthma-related quantitative traits IgE and FEV₁ (45) and regulatory pathways have now been identified (46, 47), a deeper understanding is still needed. Further, age-related changes influencing asthma risk and asthma severity throughout life may also exhibit sex-specific effects. Lung function, for instance, peaks ∼20 years (48) and then declines with age (49). Type 2 inflammation was recently shown to increase with age, particularly in the fourth decade (34).

Our study revealed a gene-sex-aging relationship between CRTh2 533116 G > A and various parameters related to asthma severity. Indeed, we found that in older females the AA genotype was associated with having a diagnosis of severe asthma, increased likelihood of taking high dose ICS and having FEV₁ (% predicted) below the normal range (<80%). Many aspects of type 2 inflammation were also increased in females, but not males with the AA genotype, though the proportion of Th2 cells was only significantly higher in those older than 45 years. The mechanism underlying why the AA genotype associates with increased type 2 inflammation only in females remains unclear but we hypothesize that this may be related to their stronger propensity to develop type 2 responses (50–52). Estrogen has been shown to induce IL-4 and IL-13 expression in murine models of asthma (53, 54), while serum estradiol during the luteal phase was positively correlated with airway levels of IL-5 (55). In males, the AA genotype did not influence risk of having severe asthma, coincide with higher ICS usage or type 2 inflammation, which in turn may be due to androgen levels; androgens have been shown to suppress type 2 inflammation (50–52). In our study males were more likely to be taking systemic corticosteroids, which could blunt detection of type 2 responses. Removing these subjects from the analysis, however, gave similar results and so it seems unlikely the lack of association of AA with increased type 2 inflammation in males is due to corticosteroid usage. Although peripheral detection of type 2 inflammation (blood eosinophil counts) has become the standard for determining eligibility for anti-type 2 medications (56), airway assessment is considered more accurate (57). As such, additional studies examining association of CRTh2 rs533116 AA and airway levels of CRTh2⁺ cells in males and females are still needed to fully understand its influence on type 2 inflammation. Intriguingly, there may be sex differences in the role of peripheral vs. lung cell contribution to type 2 inflammation, since androgens directly reduced ILC2-mediated (50, 51), but indirectly reduced Th2 cell-mediated type 2 inflammation (52).

In terms of the influence of the CRTh2 rs533116 G > A polymorphism on FEV₁, the combined analysis of males and females showed a trend for significance for the AA genotype to influence risk for lower FEV₁ (recessive model, p = 0.055) but, in contrast to diagnosis of severe asthma, no effect was observed in sex-stratified analysis. Instead, we observed a graded gene dose effect of 1 and 2 copies of the A allele (GA and AA genotypes) on FEV₁. A number of studies now report preference for assessing genetic association using a genotypic model as they can be more sensitive than the additive model, even when there is an intermediate heterozygote effect (58, 59). Association between the FTO (Fat mass and obesity-associated) genetic variant rs1421085 and severe obesity, for example, was higher using a comparison of homozygous genotypes than the additive model (60). Using this approach, we found the AA genotype was associated with increased likelihood of having FEV₁ below the normal range (<80%) than those of GG genotype, in both females and males. Taken together, our results suggest the influence of CRTh2 rs533116 G > A on asthma severity emerges in females as they age, when declining lung function (49) interacts with genotype- and/or age-mediated increases in type 2 inflammation (34). The fact that in males this polymorphism associates with lower lung function, but not type 2 inflammation and asthma severity, suggests the potential for sex differences in asthma etiology. Moreover, it highlights the importance of future work assessing whether CRTh2 rs533116 G > A directly influences lung function per se and its impact on other lung conditions.

PGD₂ release from mast cells serves as a chemotactic factor drawing CRTh2-expressing inflammatory cells to allergen-exposed tissues (8, 22, 61, 62). On Th2 cells, PGD₂ activation of CRTh2 mediates production of IL-4, IL-5 and IL-13, amplifying the type 2 response (21). Our previous study showed association of CRTh2 rs533116 AA with higher eosinophil and T cell expression of CRTh2 in young adults with self-reported asthma (28). Here we report this association persists in physician diagnosed asthma. Whole genome discovery of regulatory regions identified H3K4me1 binding to the CRTh2 rs533116 locus, indicating this polymorphism resides within an enhancer region [HaploReg¹; (63)]. Transcription factor binding site analysis also revealed that the CRTh2 rs533116 G allele (but not A allele) contains a putative NFAT site. Mechanistically, this could result in increased CRTh2 expression since we previously showed this transcription factor reduces CRTh2 transcription (64). If so, CRTh2 rs533116AA could drive unfettered transcription resulting in high level CRTh2 mRNA and more CRTh2-expressing cells.

Differences in the influence of this polymorphism on type 2 inflammation and FEV₁ could be the result of genetic effects on different cell types. CRTh2 is expressed on airway smooth muscle and the peripheral nervous system (65). A direct role for CRTh2 in regulating asthma-related changes in lung function would be in line with a report of association between other CRTh2 polymorphisms and lower PC₂₀ (26). In silico analysis of the CRTh2 rs533116 G > A locus suggests the A allele reduces binding of MEF2 [HaploReg, see text footnote 1; (63)], a transcription factor expressed by smooth muscle cells that regulates contraction (66). PGD₂ activation of CRTh2 induces myocyte migration (9) and fibrosis (67) and so an enhanced CRTh2 pathway may increase airway smooth muscle area and tissue remodeling associated with asthma-related lower lung function (68). Observed differences in genotypic effects on type 2 inflammation (recessive) and FEV₁ (dominant), may be due to chromatin accessibility at the CRTh2 locus differing within female and male cells and/or different cell types. Ultimately these nuclear environments may regulate the impact of CRTh2 rs533116 G > A.

Despite this novel finding of biological sex and age influencing the contribution of CRTh2 rs533116 G > A to asthma severity and type 2 inflammation, this study had some limitations. While data were acquired from two separate centers, it had to be pooled to achieve sufficient power to detect genotypic differences. We also lacked power to detect differences in some indices of symptom control, as number of exacerbations and cycles of oral corticosteroid/year were acquired only in a subset of subjects. Another important question is whether CRTh2 rs533116 G > A associates with a particular subtype of Type 2 high asthma. Previous cluster analysis identified two subtypes of female-dominant asthma characterized by moderate airflow obstruction (>68% predicted FEV₁) and type 2 inflammation (either IgE or airway eosinophils) but differing in age-of-onset (< or >40 years of age) (69). Unfortunately, we could not assess the influence of CRTh2 rs533116 AA on age-of-onset or severity of asthma within each subtype due to low sample size (U of A, n = 59, AA < 40 = 4; AA > 40 = 2). Since the PGD₂-CRTh2 pathway plays a role in both IgE-mediated and eosinophilic type 2 inflammation, CRTh2 rs533116AA may serve as a biomarker of severity for both early and late onset asthma. In light of these remaining questions, a follow up validation study including these outcomes and features of asthma is warranted and would extend our understanding of the current findings.

Severe asthma comprises about 5–10% of asthmatics, though this group accounts for more than 50% of the asthma related total health costs due to hospital admissions, use of emergency service and unscheduled physician visits (70–72). Our study reveals that CRTh2 rs533116 AA associates with having more type 2 inflammation, need for ICS and severe asthma diagnosis exclusively in females. As such, this variant may serve as a sex-specific biomarker for type 2 high severe asthma and help guide precision delivery of anti-Th2 therapy.

Data availability statement

The raw data supporting the results and conclusions of this article will be made available by the corresponding author on reasonable request.

Ethics statement

The studies involving human participants were reviewed and approved by University of Alberta Human Ethics Committee and Western University Human Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author contributions

NS acquired all flow cytometry and qRT-PCR data and performed genotyping of the U of A dataset and performed all statistical analyses. CM, CL, and RK oversaw recruitment of patients and acquisition of clinical data for the WU dataset. RK also oversaw genotyping of the WU dataset. HV oversaw recruitment and genotyping of all U of A patients and provided invaluable discussion for data interpretation and writing of the manuscript. LC conceived the project, coordinated the team, oversaw acquisition and analysis of the flow cytometry and qRT-PCR data, and wrote the final manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was made possible with grant funding to LC from CIHR (154849), CAAIF and Western University; RK from CIHR (154849) and the Wolfe Medical Research Chair in Pharmacogenomics; HV from CIHR (133475) and the GlaxoSmithKline-CIHR Research Chair in Airway Inflammation (143901). The CIHR/GSK Chair was established with support from GSK and CIHR given to the University of Alberta. The funding organizations are not involved in the selection of the Chair holder or any of the Chair activities.

Acknowledgments

We thank Angela Hillaby and Miranda Bowen for coordinating visits to the University of Alberta (UA) asthma clinic and their help with data collection. Flow cytometry was performed at the UA Facility, which receives financial support from the Faculty of Medicine & Dentistry and Canada Foundation for Innovation (CFI) awards. We would like to acknowledge Dr. Michaela Schedel, University Medicine Essen-Ruhrlandklinik for invaluable discussions and review of the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2022.970495/full#supplementary-material

Footnotes

1. ^ https://pubs.broadinstitute.org/mammals/haploreg

References

1. Woodruff PG, Modrek B, Choy DF, Jia G, Abbas AR, Ellwanger A, et al. T-helper type 2-driven inflammation defines major subphenotypes of asthma. Am J Respir Crit Care Med. (2009) 180:388–95. doi: 10.1164/rccm.200903-0392OC

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Wenzel SE. Asthma phenotypes: the evolution from clinical to molecular approaches. Nat Med. (2012) 18:716–25. doi: 10.1038/nm.2678

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Mosmann TR, Cherwinski H, Bond MW, Giedlin MA, Coffman RL. Two types of murine helper T cell clone. I. Definition according to profiles of lymphokine activities and secreted proteins. J Immunol. (1986) 136:2348–57.

Google Scholar

4. Mjösberg J, Spits H. Type 2 innate lymphoid cells-new members of the “type 2 franchise” that mediate allergic airway inflammation. Eur J Immunol. (2012) 42:1093–6. doi: 10.1002/eji.201242549

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Cosmi L, Annunziato F, Galli MIG, Maggi RME, Nagata K, Romagnani S, et al. CRTH2 is the most reliable marker for the detection of circulating human type 2 Th and type 2 T cytotoxic cells in health and disease. Eur J Immunol. (2000) 30:2972–9. doi: 10.1002/1521-4141(200010)30:10<2972::AID-IMMU2972>3.0.CO;2-#

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Messi M, Giacchetto I, Nagata K, Lanzavecchia A, Natoli G, Sallusto F, et al. Memory and flexibility of cytokine gene expression as separable properties of human T(H)1 and T(H)2 lymphocytes. Nat Immunol. (2003) 4:78–86. doi: 10.1038/ni872

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Nagata K, Tanaka K, Ogawa K, Kemmotsu K, Imai T, Yoshie O, et al. Selective expression of a novel surface molecule by human Th2 cells in vivo. J Immunol. (1999) 162:1278–86.

Google Scholar

8. Nagata K, Hirai H, Tanaka K, Ogawa K, Aso T, Sugamura K, et al. CRTH2, an orphan receptor of T-helper-2-cells, is expressed on basophils and eosinophils and responds to mast cell-derived factor(s). FEBS Lett. (1999) 459:195–9. doi: 10.1016/s0014-5793(99)01251-x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Parameswaran K, Radford K, Fanat A, Stephen J, Bonnans C, Levy BD, et al. Modulation of human airway smooth muscle migration by lipid mediators and Th-2 cytokines. Am J Respir Cell Mol Biol. (2007) 37:240–7. doi: 10.1165/rcmb.2006-0172OC

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Stinson SE, Amrani Y, Brightling CE. D prostanoid receptor 2 (chemoattractant receptor-homologous molecule expressed on TH2 cells) protein expression in asthmatic patients and its effects on bronchial epithelial cells. J Allergy Clin Immunol. (2015) 135:395–406. doi: 10.1016/j.jaci.2014.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Woodland DL, Kohlmeier JE. Migration, maintenance and recall of memory T cells in peripheral tissues. Nat Rev Immunol. (2009) 9:153–61. doi: 10.1038/nri2496

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Mackay CR, Marston WL, Dudler L. Naive and memory T cells show distinct pathways of lymphocyte recirculation. J Exp Med. (1990) 171:801–17. doi: 10.1084/jem.171.3.801

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mojtabavi N, Dekan G, Stingl G, Epstein MM. Long-lived Th2 memory in experimental allergic asthma. J Immunol. (2002) 169:4788–96. doi: 10.4049/jimmunol.169.9.4788

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Smith SG, Chen R, Kjarsgaard M, Huang C, Oliveria JP, O’Byrne PM, et al. Increased numbers of activated group 2 innate lymphoid cells in the airways of patients with severe asthma and persistent airway eosinophilia. J Allergy Clin Immunol. (2016) 137:75–86. doi: 10.1016/j.jaci.2015.05.037

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Hammad H, Lambrecht BN. Barrier epithelial cells and the control of type 2 immunity. Immunity. (2015) 43:29–40. doi: 10.1016/j.immuni.2015.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lambrecht BN, Hammad H. Allergens and the airway epithelium response: gateway to allergic sensitization. J Allergy Clin Immunol. (2014) 134:499–507. doi: 10.1016/j.jaci.2014.06.036

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Palikhe NS, Laratta C, Nahirney D, Vethanayagam D, Bhutani M, Vliagoftis H, et al. Elevated levels of circulating CD4(+) CRTh2(+) T cells characterize severe asthma. Clin Exp Allergy. (2016) 46:825–36. doi: 10.1111/cea.12741

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Fajt ML, Gelhaus SL, Freeman B, Uvalle CE, Trudeau JB, Holguin F, et al. Prostaglandin D2 pathway upregulation: relation to asthma severity, control, and TH2 inflammation. J Allergy Clin Immunol. (2013) 131:1504–12. doi: 10.1016/j.jaci.2013.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hirai H, Tanaka K, Yoshie O, Ogawa K, Kenmotsu K, Takamori Y, et al. Prostaglandin D2 selectively induces chemotaxis in T helper type 2 cells, eosinophils, and basophils via seven-transmembrane receptor CRTH2. J Exp Med. (2001) 193:255–61. doi: 10.1084/jem.193.2.255

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Monneret G, Gravel S, Diamond M, Rokach J, Powell WS. Prostaglandin D2 is a potent chemoattractant for human eosinophils that acts via a novel DP receptor. Blood. (2001) 98:1942–8. doi: 10.1182/blood.v98.6.1942

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Xue L, Gyles SL, Wettey FR, Gazi L, Townsend E, Hunter MG, et al. Prostaglandin D2 causes preferential induction of proinflammatory Th2 cytokine production through an action on chemoattractant receptor-like molecule expressed on Th2 cells. J Immunol. (2005) 175:6531–6. doi: 10.4049/jimmunol.175.10.6531

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lewis RA, Soter NA, Diamond PT, Austen KF, Oates JA, Roberts LJ, et al. Prostaglandin D2 generation after activation of rat and human mast cells with anti-IgE. J Immunol. (1982) 129:1627–31.

Google Scholar

23. Eguchi M, Kariya S, Okano M, Higaki T, Makihara S, Fujiwara T, et al. Lipopolysaccharide induces proinflammatory cytokines and chemokines in experimental otitis media through the prostaglandin D2 receptor (DP)-dependent pathway. Clin Exp Immunol. (2011) 163:260–9. doi: 10.1111/j.1365-2249.2010.04292.x

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Singh AM, Busse WW. Asthma exacerbations 2: aetiology. Thorax. (2006) 61:809–16. doi: 10.1136/thx.2005.045179

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Earl CS, An S, Ryan RP. The changing face of asthma and its relation with microbes. Trends Microbiol. (2015) 23:408–18. doi: 10.1016/j.tim.2015.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Huang JL, Gao PS, Mathias RA, Yao TC, Chen LC, Kuo ML, et al. Sequence variants of the gene encoding chemoattractant receptor expressed on Th2 cells (CRTH2) are associated with asthma and differentially influence mRNA stability. Hum Mol Genet. (2004) 13:2691–7. doi: 10.1093/hmg/ddh279

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Cameron L, Webster RB, Strempel JM, Kiesler P, Kabesch M, Ramachandran H, et al. Th2 cell-selective enhancement of human IL13 transcription by IL13-1112C>T, a polymorphism associated with allergic inflammation. J Immunol. (2006) 177:8633–42. doi: 10.4049/jimmunol.177.12.8633

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Campos Alberto E, Maclean E, Davidson C, Palikhe NS, Storie J, Tse C, et al. The single nucleotide polymorphism CRTh2 rs533116 is associated with allergic asthma and increased expression of CRTh2. Allergy. (2012) 67:1357–64. doi: 10.1111/all.12003

PubMed Abstract | CrossRef Full Text | Google Scholar

29. The ENFUMOSA Study Group. The ENFUMOSA cross-sectional European multicentre study of the clinical phenotype of chronic severe asthma. European network for understanding mechanisms of severe asthma. Eur Respir J. (2003) 22:470–7. doi: 10.1183/09031936.03.00261903

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Chung KF, Wenzel SE, Brozek JL, Bush A, Castro M, Sterk PJ, et al. International ERS/ATS guidelines on definition, evaluation and treatment of severe asthma. Eur Respir J. (2014) 43:343–73. doi: 10.1183/09031936.00202013

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Shrestha Palikhe N, Wu Y, Konrad E, Gandhi VD, Rowe BH, Vliagoftis H, et al. Th2 cell markers in peripheral blood increase during an acute asthma exacerbation. Allergy. (2021) 76:281–90. doi: 10.1111/all.14543

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Cameron L, Depner M, Kormann M, Klopp N, Illig T, Von Mutius E, et al. Genetic variation in CRTh2 influences development of allergic phenotypes. Allergy. (2009) 64:1478–85. doi: 10.1111/j.1398-9995.2009.02053.x

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Zein JG, Dweik RA, Comhair SA, Bleecker ER, Moore WC, Peters SP, et al. Asthma is more severe in older adults. PLoS One (2015) 10:e0133490. doi: 10.1371/journal.pone.0133490

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Peters MC, Kerr S, Dunican EM, Woodruff PG, Fajt ML, Levy BD, et al. Refractory airway type 2 inflammation in a large subgroup of asthmatic patients treated with inhaled corticosteroids. J Allergy Clin Immunol. (2019) 143:104–13.e14. doi: 10.1016/j.jaci.2017.12.1009

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Renz H, von Mutius E, Brandtzaeg P, Cookson WO, Autenrieth IB, Haller D. Gene-environment interactions in chronic inflammatory disease. Nat Immunol. (2011) 12:273–7. doi: 10.1038/ni0411-273

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Castro-Giner F, Kauffmann F, de Cid R, Kogevinas M. Gene-environment interactions in asthma. Occup Environ Med. (2006) 63:776–86. doi: 10.1136/oem.2004.019216

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ober C, Loisel DA, Gilad Y. Sex-specific genetic architecture of human disease. Nat Rev Genet. (2008) 9:911–22. doi: 10.1038/nrg2415

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wang K, Basu M, Malin J, Hannenhalli S. A transcription-centric model of SNP-age interaction. PLoS Genet. (2021) 17:e1009427. doi: 10.1371/journal.pgen.1009427

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Zein JG, Erzurum SC. Asthma is different in women. Curr Allergy Asthma Rep. (2015) 15:28. doi: 10.1007/s11882-015-0528-y

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Almqvist C, Worm M, Leynaert B. Impact of gender on asthma in childhood and adolescence: a GA2LEN review. Allergy. (2008) 63:47–57. doi: 10.1111/j.1398-9995.2007.01524.x

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Schatz M, Meckley LM, Kim M, Stockwell BT, Castro M. Asthma exacerbation rates in adults are unchanged over a 5-year period despite high-intensity therapy. J Allergy Clin Immunol Pract. (2014) 2:570–4.e1. doi: 10.1016/j.jaip.2014.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

42. El-Husseini ZW, Gosens R, Dekker F, Koppelman GH. The genetics of asthma and the promise of genomics-guided drug target discovery. Lancet Respir Med. (2020) 8:1045–56. doi: 10.1016/S2213-2600(20)30363-5

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Mersha TB, Martin LJ, Biagini Myers JM, Kovacic MB, He H, Lindsey M, et al. Genomic architecture of asthma differs by sex. Genomics. (2015) 106:15–22.

Google Scholar

44. Myers RA, Scott NM, Gauderman WJ, Qiu W, Mathias RA, Romieu I, et al. Genome-wide interaction studies reveal sex-specific asthma risk alleles. Hum Mol Genet. (2014) 23:5251–9. doi: 10.1093/hmg/ddu222

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Ober C, Pan L, Phillips N, Parry R, Kurina LM. Sex-specific genetic architecture of asthma-associated quantitative trait loci in a founder population. Curr Allergy Asthma Rep. (2006) 6:241–6. doi: 10.1007/s11882-006-0041-4

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Gautam Y, Afanador Y, Abebe T, López JE, Mersha TB. Genome-wide analysis revealed sex-specific gene expression in asthmatics. Hum Mol Genet. (2019) 28:2600–14.

Google Scholar

47. Schieck M, Schouten JP, Michel S, Suttner K, Toncheva AA, Gaertner VD, et al. Doublesex and mab-3 related transcription factor 1 (DMRT1) is a sex-specific genetic determinant of childhood-onset asthma and is expressed in testis and macrophages. J Allergy Clin Immunol. (2016) 138:421–31. doi: 10.1016/j.jaci.2015.12.1305

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Sharma G, Goodwin J. Effect of aging on respiratory system physiology and immunology. Clin Interv Aging. (2006) 1:253–60.

Google Scholar

49. Thomas ET, Guppy M, Straus SE, Bell KJL, Glasziou P. Rate of normal lung function decline in ageing adults: a systematic review of prospective cohort studies. BMJ Open (2019) 9:e028150. doi: 10.1136/bmjopen-2018-028150

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Cephus JY, Stier MT, Fuseini H, Yung JA, Toki S, Bloodworth MH, et al. Testosterone attenuates group 2 innate lymphoid cell-mediated airway inflammation. Cell Rep. (2017) 21:2487–99. doi: 10.1016/j.celrep.2017.10.110

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Laffont S, Blanquart E, Savignac M, Cénac C, Laverny G, Metzger D, et al. Androgen signaling negatively controls group 2 innate lymphoid cells. J Exp Med. (2017) 214:1581–92. doi: 10.1084/jem.20161807

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Fuseini H, Yung JA, Cephus JY, Zhang J, Goleniewska K, Polosukhin VV, et al. Testosterone decreases house dust mite-induced type 2 and IL-17A-mediated airway inflammation. J Immunol. (2018) 201:1843–54. doi: 10.4049/jimmunol.1800293

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Cai Y, Zhou J, Webb DC. Estrogen stimulates Th2 cytokine production and regulates the compartmentalisation of eosinophils during allergen challenge in a mouse model of asthma. Int Arch Allergy Immunol. (2012) 158:252–60. doi: 10.1159/000331437

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Lambert KC, Curran EM, Judy BM, Milligan GN, Lubahn DB, Estes DM, et al. Estrogen receptor alpha (ERalpha) deficiency in macrophages results in increased stimulation of CD4+ T cells while 17beta-estradiol acts through ERalpha to increase IL-4 and GATA-3 expression in CD4+ T cells independent of antigen presentation. J Immunol. (2005) 175:5716–23. doi: 10.4049/jimmunol.175.9.5716

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Semik-Orzech A, Skoczyński S, Pierzchała W. Serum estradiol concentration, estradiol-to-progesterone ratio and sputum IL-5 and IL-8 concentrations are increased in luteal phase of the menstrual cycle in perimenstrual asthma patients. Eur Ann Allergy Clin Immunol. (2017) 49:161–70. doi: 10.23822/eurannaci.1764-1489.09

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Rabe KF, FitzGerald JM, Bateman ED, Castro M, Pavord ID, Maspero JF, et al. Dupilumab is effective in patients with moderate-to-severe uncontrolled GINA-Defined type 2 asthma irrespective of an allergic asthma phenotype. J Allergy Clin Immunol Pract. (2022) (in press). doi: 10.1016/j.jaip.2022.06.036

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Koenderman L, Hassani M, Mukherjee M, Nair P. Monitoring eosinophils to guide therapy with biologics in asthma: does the compartment matter? Allergy. (2021) 76:1294–7. doi: 10.1111/all.14700

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Guindo-Martinez M, Amela R, Bonàs-Guarch S, Puiggròs M, Salvoro C, Miguel-Escalada I, et al. The impact of non-additive genetic associations on age-related complex diseases. Nat Commun. (2021) 12:2436. doi: 10.1038/s41467-021-21952-4

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Liu HM, Amela R, Bonàs-Guarch S, Puiggròs M, Salvoro C, Miguel-Escalada I, et al. Recessive/dominant model: alternative choice in case-control-based genome-wide association studies. PLoS One (2021) 16:e0254947. doi: 10.1371/journal.pone.0254947

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Wood AR, Tyrrell J, Beaumont R, Jones SE, Tuke MA, Ruth KS, et al. Variants in the FTO and CDKAL1 loci have recessive effects on risk of obesity and type 2 diabetes, respectively. Diabetologia. (2016) 59:1214–21. doi: 10.1007/s00125-016-3908-5

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Obata T, Nagakura T, Kanbe M, Masaki T, Maekawa K, Yamashita K. IgE-anti-IgE-induced prostaglandin D2 release from cultured human mast cells. Biochem Biophys Res Commun. (1996) 225:1015–20. doi: 10.1006/bbrc.1996.1287

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Wenzel SE, Westcott JY, Smith HR, Larsen GL. Spectrum of prostanoid release after bronchoalveolar allergen challenge in atopic asthmatics and in control groups. An alteration in the ratio of bronchoconstrictive to bronchoprotective mediators. Am Rev Respir Dis. (1989) 139:450–7. doi: 10.1164/ajrccm/139.2.450

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. (2012) 40:D930–4. doi: 10.1093/nar/gkr917

PubMed Abstract | CrossRef Full Text | Google Scholar

64. MacLean Scott E, Solomon LA, Davidson C, Storie J, Palikhe NS, Cameron L. Activation of Th2 cells downregulates CRTh2 through an NFAT1 mediated mechanism. PLoS One (2018) 13:e0199156. doi: 10.1371/journal.pone.0199156

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Jandl K, Heinemann A. The therapeutic potential of CRTH2/DP2 beyond allergy and asthma. Prostaglandins Other Lipid Mediat. (2017) 133:42–8. doi: 10.1016/j.prostaglandins.2017.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Firulli AB, Miano JM, Bi W, Johnson AD, Casscells W, Olson EN, et al. Myocyte enhancer binding factor-2 expression and activity in vascular smooth muscle cells. Association with the activated phenotype. Circ Res. (1996) 78:196–204.

Google Scholar

67. Ito H, Yan X, Nagata N, Aritake K, Katsumata Y, Matsuhashi T, et al. PGD2-CRTH2 pathway promotes tubulointerstitial fibrosis. J Am Soc Nephrol. (2012) 23:1797–809. doi: 10.1681/ASN.2012020126

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Pepe C, Foley S, Shannon J, Lemiere C, Olivenstein R, Ernst P, et al. Differences in airway remodeling between subjects with severe and moderate asthma. J Allergy Clin Immunol. (2005) 116:544–9.

Google Scholar

69. Moore WC, Meyers DA, Wenzel SE, Teague WG, Li H, Li X, et al. Identification of asthma phenotypes using cluster analysis in the Severe Asthma Research Program. Am J Respir Crit Care Med. (2010) 181:315–23. doi: 10.1164/rccm.200906-0896OC

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Cisternas MG, Blanc PD, Yen IH, Katz PP, Earnest G, Eisner MD, et al. A comprehensive study of the direct and indirect costs of adult asthma. J Allergy Clin Immunol. (2003) 111:1212–8.

Google Scholar

71. Serra-Batlles J, Plaza V, Morejón E, Comella A, Brugués J. Costs of asthma according to the degree of severity. Eur Respir J. (1998) 12:1322–6. doi: 10.1183/09031936.98.12061322

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Godard P, Chanez P, Siraudin L, Nicoloyannis N, Duru G. Costs of asthma are correlated with severity: a 1-yr prospective study. Eur Respir J. (2002) 19:61–7. doi: 10.1183/09031936.02.00232001

PubMed Abstract | CrossRef Full Text | Google Scholar