Beta2 Adrenergic Receptor (ADRβ2) Haplotype Pair (2/4) Is Associated with Severe Asthma

Abstract

Background

β2 adrenergic receptor (ADRβ2) polymorphisms including ADRβ2+46G>A have been reported to cause adverse outcomes in mild asthmatics. The extent to which ADRβ2 polymorphisms and in particular their haplotypes contribute to severe asthma is unknown.

Objective

To determine the association of ADRβ2 polymorphisms and haplotypes with asthma severity.

Methods

Caucasians (n = 2979) were genotyped for 11 ADRβ2 polymorphisms. The cohort (mean age 39.6, 60% female) included 2296 non-asthmatics, 386 mild asthmatics, 172 moderate asthmatics and 125 severe asthmatics. Haplotype frequency and haplotype pair for each subject was determined using the PHASE algorithm.

Results

The three asthmatic cohorts were comparable in age and gender but were distinguishable from each other in terms of symptoms, spirometry, medication use and health care utilisation (p <0.001). None of the polymorphisms showed a genotypic or allelic association with asthma diagnosis or severity. Nine haplotypes were identified and no association was found with asthma diagnosis or severity per se. Haplotype pair 2/4 was associated with asthma severity (Trend Test, OR 1.42, p = 0.0008) but not with asthma per se. Prevalence of haplotype pair 2/2 appeared to decrease with asthma severity (Trend Test, OR 0.78, p = 0.067). Two new haplotypes were identified, occurring exclusively in asthmatics at a frequency of ≥ 1%. In addition, a positive association between carriage of ADRβ2 +523*C and increased risk of atopy was discovered.

Conclusions

ADRβ2 haplotype pair 2/4 is associated with severe asthma and is consistent with findings of poor bronchodilator response in mild asthmatics who are also haplotype 2/4.

About the Authors:

Li Ping Chung

* E-mail: [email protected]

Affiliations Molecular Genetics and Inflammation Unit, Lung Institute of Western Australia and Centre for Asthma, Allergy and Respiratory Research, University of Western Australia, Crawley, Nedlands, Western Australia, Australia, School of Medicine and Pharmacology, University of Western Australia, Crawley, Western Australia, Australia, Department of Respiratory Medicine, Royal Perth Hospital, Perth, Western Australia, Australia

Svetlana Baltic

Affiliations Molecular Genetics and Inflammation Unit, Lung Institute of Western Australia and Centre for Asthma, Allergy and Respiratory Research, University of Western Australia, Crawley, Nedlands, Western Australia, Australia, School of Medicine and Pharmacology, University of Western Australia, Crawley, Western Australia, Australia

Manuel Ferreira

Affiliation: The Queensland Institute of Medical Research, Brisbane, Queensland, Australia

Suzanna Temple

Affiliations Molecular Genetics and Inflammation Unit, Lung Institute of Western Australia and Centre for Asthma, Allergy and Respiratory Research, University of Western Australia, Crawley, Nedlands, Western Australia, Australia, School of Medicine and Pharmacology, University of Western Australia, Crawley, Western Australia, Australia

Grant Waterer

Affiliations School of Medicine and Pharmacology, University of Western Australia, Crawley, Western Australia, Australia, Department of Respiratory Medicine, Royal Perth Hospital, Perth, Western Australia, Australia

Philip J. Thompson

Affiliations Molecular Genetics and Inflammation Unit, Lung Institute of Western Australia and Centre for Asthma, Allergy and Respiratory Research, University of Western Australia, Crawley, Nedlands, Western Australia, Australia, School of Medicine and Pharmacology, University of Western Australia, Crawley, Western Australia, Australia, Department of Respiratory Medicine, Sir Charles Gairdner Hospital, Nedlands, Western Australia, Australia

Abstract

Background

β2 adrenergic receptor (ADRβ2) polymorphisms including ADRβ2+46G>A have been reported to cause adverse outcomes in mild asthmatics. The extent to which ADRβ2 polymorphisms and in particular their haplotypes contribute to severe asthma is unknown.

Objective

To determine the association of ADRβ2 polymorphisms and haplotypes with asthma severity.

Methods

Caucasians (n = 2979) were genotyped for 11 ADRβ2 polymorphisms. The cohort (mean age 39.6, 60% female) included 2296 non-asthmatics, 386 mild asthmatics, 172 moderate asthmatics and 125 severe asthmatics. Haplotype frequency and haplotype pair for each subject was determined using the PHASE algorithm.

Results

The three asthmatic cohorts were comparable in age and gender but were distinguishable from each other in terms of symptoms, spirometry, medication use and health care utilisation (p <0.001). None of the polymorphisms showed a genotypic or allelic association with asthma diagnosis or severity. Nine haplotypes were identified and no association was found with asthma diagnosis or severity per se. Haplotype pair 2/4 was associated with asthma severity (Trend Test, OR 1.42, p = 0.0008) but not with asthma per se. Prevalence of haplotype pair 2/2 appeared to decrease with asthma severity (Trend Test, OR 0.78, p = 0.067). Two new haplotypes were identified, occurring exclusively in asthmatics at a frequency of ≥ 1%. In addition, a positive association between carriage of ADRβ2 +523*C and increased risk of atopy was discovered.

Conclusions

ADRβ2 haplotype pair 2/4 is associated with severe asthma and is consistent with findings of poor bronchodilator response in mild asthmatics who are also haplotype 2/4.

Figures

Figure 3

Figure 1

Table 1

Table 2

Table 3

Table 4

Table 5

Table 6

Table 7

Figure 2

Figure 3

Figure 1

Table 1

Table 2

Citation: Chung LP, Baltic S, Ferreira M, Temple S, Waterer G, Thompson PJ (2014) Beta2 Adrenergic Receptor (ADRβ2) Haplotype Pair (2/4) Is Associated with Severe Asthma. PLoS ONE 9(4): e93695. https://doi.org/10.1371/journal.pone.0093695

Editor: Tesfaye Mersha Baye, Cincinnati Children's Hospital Medical center, United States of America

Received: November 20, 2013; Accepted: March 5, 2014; Published: April 1, 2014

Copyright: © 2014 Chung et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: The funding for this report was derived from the Lung Institute of Western Australia as part of a PhD project. There were no external sources of funding that would otherwise affect the study design, data collection and analysis, decision to publish or presentation of the manuscript.

Competing interests: All authors have declared that no competing interests exist.

Background

Asthma is a common inflammatory airways disease affecting 12% of adults in the Australian population. Approximately 10% of these have severe persistent disease and remains poorly controlled despite maximal therapy. Response to medications varies significantly between individuals. Although medication adherence is a significant issue, it is estimated that up to 70% of variability in therapeutic responses to pharmacotherapy is genetically based.[1]–[3]

β2 adrenergic receptor (ADRβ2) agonists are the most commonly prescribed and used inhaled asthma therapy. Their efficacy is partially dependent upon their molecular conformation and properties and the density, structure and conformation of the ADRβ2 on the cell surface. The expression/regulation of ADRβ2 is altered by several single nucleotide polymorphisms (SNP) within the promoter, coding and the 3′UTR domains. [4]–[6] ADRβ2+46*G (Gly16) results in increased agonist-stimulated receptor down-regulation while ADRβ2+79*G (Glu27) is highly resistant.[7] ADRβ2+491*T (Ile164) alters receptor binding affinity with reduced activation in vitro [8] and blunted in vivo responses to terbutaline in healthy and cardiac failure subjects. [9], [10]

Clinical studies in mild asthmatics have shown differential clinical response attributable to genetic variation. ADRβ2+46*A homozygotes are more likely to bronchodilate with salbutamol compared to individuals homozygous for the G allele.[11], [12] In contrast, prolonged or frequent use of short or long acting ADRβ2 agonist (SABA, LABA respectively) in mild asthmatics homozygous for ADRβ2+46*A is associated with adverse outcomes.[13]–[17]

Larger pharmacogenetic studies on ADRβ2 polymorphisms and asthma have almost exclusively involved those with mild or moderate disease while the significance of genetic variation in severe asthma remains unclear.[18]–[21] Three studies that examined the relationship between ADRβ2 polymorphisms and severe asthma produced contradictory results.[22]–[24] These include an association between severe asthma and ADRβ2+79*G (odds ratio 1.91) [22]; a fall in FEV1 and tachyphylaxis with formoterol in ADRβ2+46*G homozygotes [23]; and no association between ADRβ2+46 or +79 genotype in fatal or near fatal asthma.[24] However, severe asthma is a complex disease and it is unlikely explained by single genotypes (ADRβ2+46 and +79) per se but rather requires analysis of rarer variants and haplotypes.

Linkage disequilibrium of the ADRβ2 gene dictates that alleles at different locations of the ADRβ2 gene with opposing in-vitro effects are commonly co-inherited as a haplotype and the outcome of their interaction is hard to predict. Receptors with ADRβ2+79*G polymorphism are highly resistant to agonist-induced downregulation in vitro yet co-inheritance of ADRβ2+46*G favours ADRβ2+46*G induced down-regulation.[25] It is therefore critical that ADRβ2 haplotypes are assessed in addition to individual genotypes.[26] Each individual inherits two copies of the ADRβ2 gene derived from the maternal and paternal copy of chromosome 5. This in effect means that ADRβ2 gene is inherited, not simply as genotypes or single haplotype, but as a haplotype pair (Figure 1) which makes it even more difficult to ascertain the functional and clinical significance of ADRβ2 polymorphisms.

[Figure omitted. See PDF.]

Figure 1. Schematic representation of ADRβ2 genotype, haplotype and haplotype pair.

https://doi.org/10.1371/journal.pone.0093695.g001

We hypothesized that ADRβ2 polymorphisms or haplotypes associated with altered receptor expression/regulation are more prevalent in severe asthma. As such our aim was to determine the frequency of selected SNPs and identify and compare the common haplotypes and haplotype pairs in mild, moderate and severe asthma and with healthy individuals.

Methods

Overview of Study Design

Asthmatics and healthy controls were recruited from the Lung Institute of Western Australia and the Queensland Institute of Medical Research. All participants were extensively genotyped for the ADRβ2 gene. Statistical analyses were performed to determine the relationship between ADRβ2 gene and asthma severity and presence of atopy by looking at the frequency of various genotypes, haplotypes as well as haplotype pairs in the relevant patients groups.

Subject Recruitment

Lung Institute of Western Australia (LIWA) Asthma Genetic Cohort. Non-asthmatic controls (n = 200) were recruited through contacting individuals randomly selected from the metropolitan telephone directory. Mild, moderate and severe asthmatic cohorts were similarly recruited but enriched by recruiting from specialty asthma clinics and hospital databases.

Asthma (n = 683) was defined as having a respiratory physician-diagnosis of asthma based upon clinical history, spirometry, salbutamol reversibility and/or positive methacholine challenge tests. Mild and moderate asthma were defined as published previously using a number of criteria including lung function (FEV1 % predicted, use of oral steroids in the last 12 months, daily inhaled corticosteroid dose, weekly frequency of use of rescue medication, weekly frequency of daytime symptoms, weekly frequency of nocturnal awakening due to asthma, unplanned visits to doctor in the last 12 months and hospital admissions for asthma in the last 12 months.[27], [28] Patients were classified as having mild asthma if they have normal lung function (FEV1 greater than 80% predicted), frequency of daytime or night-time symptoms and use of rescue medication use less than twice per week with no unplanned medical visits or hospitalisations in the previous 12 months while treated with low dose inhaled corticosteroid therapy. Severe refractory asthma patients had to meet the American Thoracic Society (ATS) Criteria [29] and to have experienced at least 2 asthma-related hospital admissions within any 12 month period between January 2003 and February 2010. Patients with symptoms and treatment intermediate to mild and severe classification were labelled as moderate asthmatics.

All participants completed a questionnaire which included an in depth description of their symptoms and medication use. They were reviewed by a respiratory physician who confirmed their diagnosis and severity. All participants (including healthy controls) underwent spirometry in accordance with the ATS guidelines as well as allergen skin prick testing to 5 common aeroallergens: cat, dog, house dust mite, mould mix and grass pollen mix. Atopy was defined as positive skin reaction (wheal diameter ≥ 3 mm) to at least one allergen. The study was approved by the Human Research Ethics Committees at Sir Charles Gairdner and Royal Perth Hospitals. All subjects provided written consent to all study procedures including allergy and genetic testing.

Queensland Institute of Medical Research (QIMR) Non-Asthmatics from Genome Wide Association Study (GWAS).

To improve the power to detect a true genetic association, we expanded the non-asthmatic LIWA cohort to include 2096 unrelated Caucasians ascertained by QIMR and included in a recent GWAS of asthma.[30] Of these, 1248 were considered asthma-free based on a negative response to the questions “How often have you had asthma?”, “Have you ever had asthma?”, “Has a doctor ever diagnosed you as suffering from asthma?”. There was no information available on the asthma status of the remaining 848 individuals, but these were nonetheless included as unselected controls to improve power. There were no significant differences in SNP allele frequencies between the QIMR and LIWA control groups, thus supporting our decision to merge both datasets.

Genotyping and Determination of ADRβ2 Haplotype and Haplotype Pairs

Genomic DNA was extracted from buffy coat by QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). A total of 11 SNPs of the ADRβ2 gene (Table 1) were genotyped using Taqman Genotying Assay (Applied Biosystems, Foster City, CA, USA). Unequivocal results from Taqman Assays were confirmed by PCR and direct sequencing. The QIMR_610 k cohort was genotyped using the Illumina 610K Assay (Illumina, San Diego, CA, USA) as described previously.[30] A subgroup of the LIWA cohort were also included in the GWAS (n = 645), thereby allowing quality control assessment of ADRβ2 genotyping by both genotyping techniques. The QIMR_610k cohort has been genotyped for numerous other genes related to asthma and atopy and these results are not yet available for public access but may be shared upon request.

[Figure omitted. See PDF.]

Table 1. ADRβ2 Polymorphisms Genotyped.§

https://doi.org/10.1371/journal.pone.0093695.t001

ADRβ2 haplotype and haplotype pairs for each individual patient were inferred from the genotype data using the PHASE program. (v2.0). The PHASE prediction algorithm was run several times to ensure that the result was reproducible. The number and characteristics of haplotypes with frequency > 1% remained consistent. For each individual, only the most likely HP (predicted with over 80% accuracy) was considered in the final analysis.

Statistical Analysis

Patient characteristics are expressed as mean with standard deviation (SD), unless otherwise stated. Demographic data and clinical parameters were compared between all groups using ANOVA and Tukey's test or where pooled data was used, Student's t-test for continuous data or with Chi square or Fishers' exact test (where appropriate) for categorical data.

Genotypic and allelic frequencies were calculated for each patient group stratified by asthma severity or atopy. Univariate analysis was then performed using chi-square tests and logistic regression analysis using SPSS v18.0 (IBM, Armonk, CA, USA). Power analysis showed that a difference of 10% in allelic and genotypic frequency could be detected, with 80% power and Bonferroni-corrected p-value of 0.05.

Frequency of the 5 most common HPs were calculated for mild, moderate and severe asthma group. The differences were tested using the odds ratio method where only a more stringent p value of 0.01 was considered statistically significant (Bonferroni-corrected). The prevalence of the HPs in the three asthma cohorts were also assessed using a logistic regression trend test performed on SAS v9.2 (SAS Institute Inc, Cary, CA, USA) to further support the association with asthma severity.

Results

Subjects

The overall demographic and clinical characteristics of the study cohort are summarised in Table 2. The LIWA cohort included a total of 883 unrelated Caucasians (125 severe, 172 moderate, 386 mild and 200 non-asthmatics) aged between 18 and 85. These four study cohorts were similar in age but there was a greater proportion of females (p <0.001) and a higher prevalence of atopy (p < 0.001) in the pooled asthmatic cohort compared with the non-asthmatics.

[Figure omitted. See PDF.]

Table 2. Characteristics of Study Participants.

https://doi.org/10.1371/journal.pone.0093695.t002

The severe, moderate and mild asthmatics were distinguishable in all clinical parameters including frequency of symptoms, lung function, medication use and health care utilisation (p<0.0001 for ANOVA and all subsequent pairwise analysis between severity groups, Table 2).

The QIMR_610K dataset (n = 2096) included 1248 non-asthmatics with mean age 31 (62% female) and 848 individuals (55.9% female, age unknown) with no available asthma information that were considered as unselected controls to improve power.[27]

Genotype and Allelic Analysis

All participants were successfully genotyped for all SNPs (see Table 3). There was complete concordance between the Illumina and Taqman Assay results (n = 269). The allelic prevalence of ADRβ2+46*A and ADRβ2+79*C in the entire study cohort was 0.36 and 0.56, respectively. ADRβ2+491*T was only seen in the heterozygous form in 13 asthmatics (8 mild, 2 moderate, 4 severe) and 37 non-asthmatics. ADRβ2 SNPs frequencies were comparable for all study groups and were within Hardy-Weinberg equilibrium (see Table 4). ADRβ2 polymorphisms did not correlate with lung function for the whole group or when stratified by asthma severity (see Table 5). Despite lack of association between ADRβ2 genotypes and asthma, there was a strong correlation between carriage of the C allele of SNP +523 and atopy (OR 2.05, CI 1.33-4.48 for +523CC or AC genotypes compared with +523AA, p = 0.03; Table 6).

[Figure omitted. See PDF.]

Table 3. Genotypic and allelic frequency of ADRβ2 polymorphisms stratified by asthma diagnosis and severity.

https://doi.org/10.1371/journal.pone.0093695.t003

[Figure omitted. See PDF.]

Table 4. Hardy Weinberg equilibrium value (χ2) for each ADRβ2 Polymorphisms.

https://doi.org/10.1371/journal.pone.0093695.t004

[Figure omitted. See PDF.]

Table 5. Lung function stratified by ADRβ2 genotypes.

https://doi.org/10.1371/journal.pone.0093695.t005

[Figure omitted. See PDF.]

Table 6. ADRβ2+523 Genotypes are Associated with Atopy (n = 2979).

https://doi.org/10.1371/journal.pone.0093695.t006

Haplotype Analysis

The haplotype frequencies are summarised in Table 7. Nine of the ten Drysdale haplotypes [25] were identified with haplotype 2, 4 and 6 being the most common with frequencies of 43.9%, 34.5% and 15.0%, respectively. Two additional haplotypes were identified occurring at a frequency >1%, denoted as 2′ and 4′ as significant homology with Drysdale haplotype 2 and 4 existed. The differences consisted of three 5′UTR SNPs in linkage disequilibrium (D′≥0.97, r2 = 1 ADRβ2-367, −47, 20) for 2′ and one 5′UTR SNP (ADRβ2-1023) for 4′. These two new haplotypes were found exclusively in asthmatics (see Table 4). While numbers were too small for further statistical analysis there was an overall trend for a higher frequency of pooled novel haplotypes in severe asthmatics compared with moderate, mild or non-asthmatics (p = 0.07).

[Figure omitted. See PDF.]

Table 7. Haplotype frequencies stratified by asthma diagnosis and severity.

https://doi.org/10.1371/journal.pone.0093695.t007

Overall, the most common HP was 2/4 (32.0%), followed by 2/2 (19.3%), 2/6 (14.8%), 4/4 (12.2%), 4/6 (10.9%) and 6/6 (3.0%). There was no association between any HPs and asthma diagnosis per se (Figure 2a). However, the frequency of HP2/4 increased with asthma severity (Figure 2b). The prevalence of HP2/4 was significantly higher in the severe asthma group (42.4%) compared with mild (27.7%, p = 0.0008 with OR 2.40, 95%CI 1.34, 3.11). The association was strengthened by a positive Trend test for mild, moderate and severe asthma (Mild OR 1.0; Moderate OR 1.42. 95%CI 1.16, 1.74; Severe OR 2.01, 95%CI 1.64, 2.46; p =  0.0008). Consequently, there was an absolute reduction in frequency of HP2/2 and 2/6, namely 6.8% and 9.7% in severe compared with mild asthmatics respectively. Although pairwise comparisons of the frequency of HP2/2 between mild, moderate and severe asthmatics were not significant, there was a trend for reduced prevalence with increasing asthma severity (Trend Test, OR 0.76, 95%CI 0.59, 1.02, p = 0.06, Figure 1b). In comparison, although the prevalence of HP2/6 in severe asthmatics was significantly different on pairwise comparison against mild or moderate asthmatics (OR 2.82, 95%CI 1.31, 6.06 p = 0.008; OR 2.98, 95%CI 1.32, 6.80, p = 0.008), there was no observable trend with asthma severity. The HP4/4 frequency was similar across the 3 cohort groups (14.5%, 13.2% and 11.6% for mild, moderate and severe asthmatics respectively, See Figure 1). There was no statistical difference in any clinical parameters when HP2/4 was compared with other HP for the asthmatic cohort as a whole or within each severity group. Within the severe asthma group, HP2/4 was not associated with steroid resistance.

[Figure omitted. See PDF.]

Figure 2. Comparison of prevalence of ADRβ2 haplotype pairs with (A) asthma diagnosis (n = 2979) and (B) asthma severity (n = 683).

The most common haplotype pair 2/4 was associated with severe asthma (Trend test, OR 1.42, p = 0.0008).

https://doi.org/10.1371/journal.pone.0093695.g002

The association with atopy was also present but not significantly strengthened in the relevant haplotype analysis inclusive of the +523CC genotype. Prevalence of atopy was also similar amongst the 6 most common haplotype pairs. (Figure 3).

[Figure omitted. See PDF.]

Figure 3. Prevalence of Atopy Similar Amongst Common ADRβ2 Haplotype Pairs (n = 2979).

https://doi.org/10.1371/journal.pone.0093695.g003

Discussion

This is the largest cohort of asthmatics and healthy controls to be studied in detail with respect to their ADRβ2 haplotypic features. Given the heterogeneity in asthma pathogenesis and clinical manifestations, not unexpectedly, ADRβ2 polymorphisms were not associated with asthma diagnosis or severity when individual SNPs or haplotypes were examined in isolation. However the frequency of the most common haplotype pair 2/4 was significantly higher in severe asthmatics with a raised odds ratio compared with mild asthmatics. The potential role for this ADRβ2 haplotype pair in the pathogenesis of severe asthma was further supported by a significant, positive trend in the prevalence of haplotype pair 2/4 with increasing asthma severity. In addition, carriage of ADRβ2+523*C was associated with increased risk of atopy.

ADRβ2 Haplotype Pair Is Associated With Asthma Severity

The haplotypes in this study were statistically inferred via the PHASE algorithm which is a well-recognised, time-efficient method for haplotype reconstruction using unphased genotype data derived from population samples such as ours. The accuracy and predictive power of PHASE algorithm was enhanced given our large sample size, the tight linkage disequilibrium between ADRβ2 SNPs, successful genotyping for the entire cohort at each SNP locus and by performing the algorithm multiple times. In the vast majority of our cohort (99%), the final haplotype pair for each individual was predicted with over 98% specificity. Each of Drysdale's haplotypes are distinct and specifically defined by the subset of SNPs genotyped in this study.

Using a subset of 11 ADRβ2 polymorphisms similar to those used by Drysdale et al [25], we identified nine Caucasian haplotypes whose frequency was similar to published data.[4], [25] Importantly two new haplotypes (2′ and 4′) were identified exclusively in asthmatics and possibly represent a rarer subset of Drysdale's haplotypes 2 and 4 contributing to asthma severity.

While previous studies have investigated the relationship between ADRβ2 haplotypes and mild asthma [4], [25], [31], this is the first study to examine haplotypes and severe asthma. In agreement with previous studies, there was no association between single ADRβ2 haplotypes and mild asthma. Five HPs (2/4, 2/2, 4/4, 2/6, 4/6) accounted for 80% of our subjects which is similar to that reported in mild or moderate asthma. [31] More importantly, HP2/4 was associated with severe asthma in the current study with an associated decreased frequency of HP2/2 and 2/6. Whilst it may be attractive to undertake a meta-analysis of our severe asthmatic subgroup and other studies on severe asthma [22]–[24], this is unfortunately unrealistic as these studies used varying definitions of asthma severity and genotyping was limited to only three ADRβ2 SNPs.

The association between HP2/4 with asthma severity was unexpected but significantly substantiated by the findings from formal trend tests across the 3 asthma severity groups. It is worth noting that haplotype 2 and haplotype 4 are heterogeneous at the main functional locus as well as many others and it is therefore rather difficult to know how or which functional effects would be specifically driven when inherited as a haplotype pair. ADRβ2 preferentially exists as dimerised pair in vivo. Whether the coinheritance of haplotype 2 and 4 affects receptor dimerisation and therefore function of the combined receptor contributing to asthma pathogenesis needs to be considered and further explored.

There was no correlation between ADRβ2 haplotypes and any individual clinical parameter used to define asthma severity. This is not surprising given that asthma is heterogeneous and classification of disease severity can be difficult when the various parameters used in defining severity are inter-related in many ways. In this study, there were no association between the inheritance of ADRβ2 haplotypes and diagnosis of mild or moderate asthma. Inflammation plays a critical role in the pathogenesis of asthma but there is increasing evidence for considerable variability in the pattern of inflammation between individuals which most likely contributes to phenotypic differences in disease severity and treatment responses. Detailed examination of the role of ADRβ2 haplotypes in the pathogenesis of asthma is beyond the scope of this study. However, it is unlikely any one gene can fully account for the pathogenesis of severe asthma but ADRβ2 haplotypes may explain some of the therapeutic resistance seen in severe asthma. Mild asthmatics homozygous for ADRβ2+46*A bronchodilate to SABA[11] but have reduced PEFR and FEV1 and increased risk of exacerbation if SABA or LABA are used regularly.[13], [14], [16], [17] Contradictory findings on therapeutic responses in larger prospective studies [19]–[21] may be due to a lack of haplotypic analysis. The few studies that attempted to address the relationship between haplotypes and various clinical characteristics of asthma including pharmacological response in asthmatics are inconclusive. [32]–[37] Furthermore, the haplotypes in these analyses were simply defined by 2–4 SNPs which is often inadequate for more detailed distinction between ADRβ2 haplotypes defined by Drysdale and colleagues.

Data on the effect of HPs on acute bronchodilator response (BDR) is limited and controversial. Taylor et al [31] showed no association between haplotype and BDR to inhaled salbutamol in mild to moderate asthmatics weaned off LABA, SABA and ICS. In contrast, Drysdale et al [25] reported a significant association between Drysdale haplotype two and a positive response to ADRβ2 agonist. The response for HP2/4 was intermediate while HP4/4 showed a much weaker response. [25] Bronchodilator response was similarly high in HP2/6 and statistically no different to HP2/2. Whether Drysdale's apriori selection of patients with a positive BDR has led to selection bias is unclear. However, the implication that BDR can alter depending on prior treatment and HPs takes on much greater importance in severe asthma particularly with the current finding of a higher frequency of the less favourable HP2/4 and lower frequency of the more responsive HP2/2 and 2/6. [25] All of our asthmatics were selected based upon chronic disease status and not acute response to ADRβ2 agonist and therefore similar to population studies by Drysdale el al. It is possible that the loss of protective, more responsive haplotypes 2 and 6 in severe asthmatics leads to a dynamic increase in the frequency of the most common haplotype pair in any given population, thus explaining the apparent association between HP2/4 and severe asthma.

To further enhance the current study data the functional significance of ADRβ2 haplotypes with respect to in vitro receptor characteristics and response to different ADRβ2 agonists needs addressing. Panebra et al reported variable baseline ADRβ2 expression and agonist-induced downregulation between ADRβ2 haplotypes. [38] Haplotype 2 demonstrated higher basal receptor expression while haplotype 4 showed greater downregulation following exposure to the ADRβ2 agonists. Whether such observation translates to clinically important differences in acute bronchodilator response and tachyphylaxis respectively requires further investigation. Furthermore, the relative effect of each haplotype when coinherited as haplotype pair 2/4 is unknown. This would be best performed in a prospective pharmacogenetic study complementing in vitro studies to characterise the effect of haplotype pairs on receptor expression, function and regulation as well as pharmacological response. Such studies would provide important insights into whether individuals with severe asthma might benefit from alternative treatment strategies that minimise excessive use ADRβ2 agonists.

Association Between ADRβ2 +523CC Genotype and Atopy

The strong correlation between ADRβ2+523CC genotype and atopy is a novel finding. The precise mechanism for this is unclear and is beyond the scope of the current study. It is most likely that ADRβ2 +523C acts as a tag or proxy for the actual genetic variation associated with atopy. This would not be surprising given that the ADRβ2 gene is located on chromosome 5q, a region where multiple candidate gene polymorphisms has been associated with atopy and/or IgE production. These genes include interleukin 4 (IL4) [39], [40], interleukin 13 (IL13) [41], monocyte differentiation antigen CD14 [42]–[45] and serine protease inhibitor kazal 5 (SPINK5) [46], [47]. Unfortunately, there is no published data on whether polymorphisms of these candidate genes are in linkage disequilibrium with any of the ADRβ2 polymorphisms. Woszczek et al demonstrated significantly higher total IgE levels and prevalence of atopy in individuals homozygous with the ADRβ2 haplotype -47T/-20T/+46A/+79C/+252G. [48] Although their genotyping was not inclusive of SNP+523, results from the current study together with those reported in the literature [4], [25] suggest that the +523*C allele is likely to be in tight linkage equilibrium with the above haplotype.

ADRβ2+523C>A is a synonymous polymorphism in which both alleles produce the same protein without any changes in amino acid sequence. It was long assumed that synonymous SNPs are “silent” and inconsequential, as the primary sequence of the protein is retained. However, it has been shown in other disease states that synonymous polymorphisms can predict disease susceptibility if it is in linkage with another polymorphism of the same or different gene that contributes to disease pathogenesis. One such example is the association between the immune-related GTPase family M

(IRGM) gene and Crohn's inflammatory bowel disease. [49]–[52]

Conclusions

Despite widespread use of ADRβ2 agonists in the treatment of asthma, there is persistent background concern about their safety and efficacy. Functional polymorphisms of ADRβ2 potentially generate differential treatment response and suboptimal outcomes for asthmatics treated with regular SABA or LABA. Most studies have investigated only a limited number of ADRβ2 polymorphisms in mild or moderate asthmatics. However, it is the overall haplotype that best predicts receptor characteristics and treatment response. This study identified 13 new ADRβ2 haplotypes and demonstrated, for the first time, that the ADRβ2 haplotype pair 2/4 is associated with asthma severity. Further functional and clinical studies are needed to determine the role of ADRβ2 polymorphisms/haplotypes both in the development of severe asthma and in influencing therapeutic responses to ADRβ2 agonists. The positive association between ADRβ2+523CC genotype and atopy is novel. This most likely reflects its linkage with another candidate gene or causal polymorphism.

Acknowledgments

We thank Scott Gordon, Dixie Statham, Ann Eldridge, Marlene Grace, Lisa Bowdler, Steven Crooks, David Smyth, Harry Beeby, Anjali K. Henders, Patrick Danoy, Dale R. Nyholt, Peter M. Visscher, Pamela A. Madden, Andrew C. Heath, David L. Duffy, Matthew A. Brown, Grant W. Montgomery and Nicholas G. Martin, who were part of the original team that ascertained, genotyped and processed SNP data for the twins and families who participated in the various QIMR cohorts.

Author Contributions

Conceived and designed the experiments: LPC GW PT. Performed the experiments: LPC. Analyzed the data: LPC MF. Contributed reagents/materials/analysis tools: LPC MF SB ST. Wrote the paper: LPC MF GW PT.

Citation: Chung LP, Baltic S, Ferreira M, Temple S, Waterer G, Thompson PJ (2014) Beta2 Adrenergic Receptor (ADRβ2) Haplotype Pair (2/4) Is Associated with Severe Asthma. PLoS ONE 9(4): e93695. https://doi.org/10.1371/journal.pone.0093695

References

1. Drazen JM, Silverman EK, Lee TH (2000) Heterogeneity of therapeutic responses in asthma. Br Med Bull 56: 1054–1070.

2. Evans WE, McLeod HL (2003) Pharmacogenomics - drug deposition, drug targets, and side effects. N Engl J Med 348: 538–549.

3. Weiss ST, Litonjua AA, Lange C, Lazarus R, Liggett SB, et al. (2006) Overview of the pharmacogenetics of asthma treatment. The Pharmacogenomics Journal 6: 311–326.

4. Hawkins GA, Tantisira K, Meyers DA, Ampleford EJ, Moore WC, et al. (2006) Sequence, haplotype, and association analysis of ADRbeta2 in a multiethnic asthma case-control study. Am J Respir Crit Care Med 174: 1101–1109.

5. Kobilka BK, Frielle T, Dohlman HG, Bolanowski MA, Dixon RA, et al. (1987) Delineation of the intronless nature of the genes for the human and hamster beta 2-adrenergic receptor and their putative promoter regions. J Biol Chem 262: 7321–7327.

6. Reihsaus E, Innis M, MacIntyre N, Liggett SB (1993) Mutations in the gene encoding for the beta 2-adrenergic receptor in normal and asthmatic subjects. Am J Respir Cell Mol Biol 8: 334–339.

7. Green SA, Turki J, Innis M, Liggett SB (1994) Amino-terminal polymorphisms of the human beta 2-adrenergic receptor impart distinct agonist-promoted regulatory properties. Biochemistry 33: 9414–9419.

8. Green SA, Cole G, Jacinto M, Innis M, Liggett SB (1993) A polymorphism of the human beta 2-adrenergic receptor within the fourth transmembrane domain alters ligand binding and functional properties of the receptor. J Biol Chem 268: 23116–23121.

9. Barbato E, Penicka M, Delrue L, Van Durme F, De Bruyne B, et al. (2007) Thr164Ile polymorphism of the beta2-adrenergic receptor negatively modulates cardiac contractility: implications for prognosis in patients with idiopathic dilated cardiomyopathy. Heart 93: 856–861.

10. Brodde OE, Buscher R, Tellkamp R, Radke J, Dhein S, et al. (2001) Blunted cardiac responses to receptor activation in subjects with Thr164Ile beta(2)-adrenoceptors. Circulation 103: 1048–1050.

11. Martinez FD, Graves PE, Baldini M, Solomon S, Erickson R (1997) Association between genetic polymorphisms of the beta2-adrenoceptor and response to albuterol in children with and without a history of wheezing. J Clin Invest 100: 3184–3188.

12. Silverman EK, Kwiatkowski DJ, Sylvia JS, Lazarus R, Drazen JM, et al. (2003) Family-based association analysis of beta2-adrenergic receptor polymorphisms in the childhood asthma management program. J Allergy Clin Immunol 112: 870–876.

13. Israel E, Chinchilli VM, Ford JG, Boushey HA, Cherniack R, et al. (2004) Use of regularly scheduled albuterol treatment in asthma: genotype-stratified, randomised, placebo-controlled cross-over trial. Lancet 364: 1505–1512.

14. Israel E, Drazen JM, Liggett SB, Boushey HA, Cherniack RM, et al. (2001) Effect of polymorphism of the beta(2)-adrenergic receptor on response to regular use of albuterol in asthma. Int Arch Allergy Immunol 124: 183–186.

15. Palmer CN, Lipworth BJ, Lee S, Ismail T, Macgregor DF, et al. (2006) Arginine-16 beta2 adrenoceptor genotype predisposes to exacerbations in young asthmatics taking regular salmeterol. Thorax 61: 940–944.

16. Taylor DR, Drazen JM, Herbison GP, Yandava CN, Hancox RJ, et al. (2000) Asthma exacerbations during long term beta agonist use: influence of beta(2) adrenoceptor polymorphism. Thorax 55: 762–767.

17. Wechsler ME, Lehman E, Lazarus SC, Lemanske RF Jr, Boushey HA, et al. (2006) beta-Adrenergic Receptor Polymorphisms and Response to Salmeterol. Am J Respir Crit Care Med 173: 519–526.

18. Anderson HR, Ayres JG, Sturdy PM, Bland JM, Butland BK, et al.. (2005) Bronchodilator treatment and deaths from asthma: case-control study. BMJ 330: 117–.

19. Bleecker ER, Postma DS, Lawrance RM, Meyers DA, Ambrose HJ, et al. (2007) Effect of ADRB2 polymorphisms on response to longacting beta2-agonist therapy: a pharmacogenetic analysis of two randomised studies. Lancet 370: 2118–2125.

20. Bleecker ER, Yancey SW, Baitinger LA, Edwards LD, Klotsman M, et al. (2006) Salmeterol response is not affected by beta2-adrenergic receptor genotype in subjects with persistent asthma. J Allergy Clin Immunol 118: 809–816.

21. Wechsler ME, Kunselman SJ, Chinchilli VM, Bleecker E, Boushey HA, et al. (2009) Effect of [beta]2-adrenergic receptor polymorphism on response to longacting [beta]2 agonist in asthma (LARGE trial): a genotype-stratified, randomised, placebo-controlled, crossover trial. The Lancet 374: 1754–1764.

22. Holloway JW, Dunbar PR, Riley GA, Sawyer GM, Fitzharris PF, et al. (2000) Association of beta2-adrenergic receptor polymorphisms with severe asthma. Clin Exp Allergy 30: 1097–1103.

23. Tan S, Hall IP, Dewar J, Dow E, Lipworth B (1997) Association between beta 2 adrenoceptor polymorphisms and susceptibility to bronchodilator desensitisation in moderately severe stable asthmatics. Lancet 350: 995–999.

24. Weir TD, Mallek NOA, Sandford AJ, Bai TR, Awadh N, et al. (1998) beta 2-Adrenergic Receptor Haplotypes in Mild, Moderate and Fatal/Near Fatal Asthma. Am J Respir Crit Care Med 158: 787–791.

25. Drysdale CM, McGraw DW, Stack CB, Stephens JC, Judson RS, et al. (2000) Complex promoter and coding region beta 2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness. Proc Natl Acad Sci U S A 97: 10483–10488.

26. Chung L, Waterer G, Thompson P (2010) Pharmacogenetics of beta2 adrenergic receptor gene polymorphisms, long acting beta agonists and asthma. Clin Exp Allergy in press.

27. Ferreira MA, Matheson MC, Duffy D, Marks GM, Hui J, et al. (2011) Identification of IL6R and chromosome 11p13.5 as risk loci for asthma. Lancet 378: 1006–1014.

28. Kedda M-A, Shi J, Duffy D, Phelps S, Yang I, et al. (2004) Characterization of two polymorphisms in the leukotriene C4 synthase gene in an Australian population of subjects with mild, moderate, and severe asthma. J Allergy Clin Immunol 113: 889–895.

29. Proceedings of the ATS Workshop on Refractory Asthma. Current Understanding, Recommendations, and Unanswered Questions. Am J Respir Crit Care Med 162: 2341–2351.

30. Wan YI, Shrine NRG, Soler Artigas M, Wain LV, Blakey JD, et al. (2012) Genome-wide association study to identify genetic determinants of severe asthma. Thorax 67: 762–768.

31. Taylor DR, Epton MJ, Kennedy MA, Smith AD, Iles S, et al. (2005) Bronchodilator Response in Relation to {beta}2-Adrenoceptor Haplotype in Patients with Asthma. Am J Respir Crit Care Med 172: 700–703.

32. Qui Y, Zhang X, Qin Y, Yin K, Zhang D (2010) Beta(2)-adrenergic receptor haplotype/polymorphisms and asthma susceptibility and clinical phenotype in a Chinese Han population. Allergy Asthma Proc 31: 91–97.

33. Manoharan A, Anderson W, Lipworth BJ (2013) Influence of beta(2)-adrenergic receptor polymorphism on methacholine hyperresponsiveness in asthmatics patients. Ann Allergy Asthma Immunol 110: 161–164.

34. Carroll C, Sala K, Zucker A, Schramm C (2013) Beta2-adrenergic receptor haplotype linked to intubation and mechanical ventilation in children with asthma. J Asthma 49: 563–568.

35. Moore P, Williams S, Gebretsadit T, Jiang L, Minton P, et al. (2008) beta(2) adrenergic receptor promoter haplotype influences spirometric response during an acute asthma exacerbation. Clin Transl Sci 1: 155–161.

36. Wilson A, Gray R, Hall I, Lipworth BJ (2006) The effect of beta2-adrenoceptor haplotypes on bronchial hyper-responsiveness in patients with asthma. Allergy 61: 254–259.

37. Torjussen T, Munthe-Kaas M, Mowinckel P, Carlsen K, Undlien D, et al. (2013) Childhood lung function and the association with beta2-adrenergic receptor haplotypes. Acta Paediatr 102: 727–731.

38. Panebra A, Wang WC, Malone MM, Pitter DRG, Weiss ST, et al. (2010) Common ADRβ2 Haplotypes Derived from 26 Polymorphic Sites Direct β2-Adrenergic Receptor Expression and Regulation Phenotypes. PLoS ONE 5: e11819.

39. Beghé B, Barton S, Rorke S, Peng Q, Sayers I, et al. (2003) Polymorphisms in the interleukin-4 and interleukin-4 receptor α chain genes confer susceptibility to asthma and atopy in a Caucasian population. Clin Exp Allergy 33: 1111–1117.

40. Kabesch M, Tzotcheva I, Carr D, Hofler C, Weiland SK, et al. (2003) A complete screening of the IL4 gene: Novel polymorphisms and their association with asthma and IgE in childhood. J Allergy Clin Immunol 112: 893–898.

41. Heinzmann A, Jerkic S-P, Ganter K, Kurz T, Blattmann S, et al. (2003) Association study of the IL13 variant Arg110Gln in atopic diseases and juvenile idiopathic arthritis. J Allergy Clin Immunol 112: 735–739.

42. Bučková D, Hollá LI, Schüller M, Znojil V, Vácha J (2003) Two CD14 promoter polymorphisms and atopic phenotypes in Czech patients with IgE-mediated allergy. Allergy 58: 1023–1026.

43. Leung TF, Tang NLS, Sung YM, Li AM, Wong GWK, et al. (2003) The C−159T polymorphism in the CD14 promoter is associated with serum total IgE concentration in atopic Chinese children. Pediatr Allergy Immunology 14: 255–260.

44. O'Donnell AR, Toelle BG, Marks GB, Hayden CM, Laing IA, et al. (2004) Age-specific Relationship between CD14 and Atopy in a Cohort Assessed from Age 8 to 25 Years. Am J Respir Crit Care Med 169: 615–622.

45. Woo JG, Assa'ad A, Heizer AB, Bernstein JA, Khurana Hershey GK (2003) The -159 C>T polymorphism of CD14 is associated with nonatopic asthma and food allergy. J Allergy Clin Immunol 112: 438–444.

46. Moffatt M (2004) SPINK5: A gene for atopic dermatitis and asthma. Clin Exp Allergy 34: 325–327.

47. Kabesch M, Carr D, Weiland SK, Von Mutius E (2004) Association between polymorphisms in serine protease inhibitor, kazal type 5 and asthma phenotypes in a large German population sample. Clin Exp Allergy 34: 340–345.

48. Woszczek G, Borowiec M, Ptasinska A, Kosinski S, Pawliczak R, et al. (2005) Beta2-ADR haplotypes/polymorphisms associate with bronchodilator response and total IgE in grass allergy. Allergy 60: 1412–1417.

49. Deretic V (2011) Autophage as an innate immune paradign: expanding the scope and repertoire of pattern recognition receptors. Curr Opin Immunol: Nov 14 [Epub ahead of print].

50. Brest P, Corcelle E, Cesaro A, Chargui A, Belaid A, et al. (2010) Autophagy and Crohn's disease: at the crossroads of infection, inflammation, immunity, and cancer. Curr Mol Med 10: 486–502.

51. McCarroll SA, Huett A, Kuballa P, Chilewski SD, Landry A, et al. (2008) Deletion polymorphism upstream of IRGM associated with altered IRGM expression and Crohn's disease. Nat Genet 40: 1107–1112.

52. Parkes M, Barrett JC, Prescott NJ, Tremelling M, Anderson CA, et al. (2007) Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility. Nat Genet 39: 830–832.

53. Westland R, van Veen A, Jansen HM, Jonkers RE, Wierenga EA (2004) Limited impact of multiple 5' single-nucleotide polymorphisms on the transcriptional control of the human beta 2-adrenoceptor gene. Immunogenetics 56: 625–630.

54. McGraw DW, Forbes SL, Kramer LA, Liggett SB (1998) Polymorphisms of the 5' leader cistron of the human beta2-adrenergic receptor regulate receptor expression. J Clin Invest 102: 1927–1932.