Journal List > Allergy Asthma Immunol Res > v.9(6) > 1052693

TOOLS
Kim, Shin, Park, Uh, Chang, Bae, Cho, Park, Yoon, Choi, Kim, and Park: Clinical Characteristics of Exacerbation-Prone Adult Asthmatics Identified by Cluster Analysis
Original Article

Allergy, Asthma & Immunology Research 2017; 9(6): 483-490.

Published online: 8 August 2017

DOI: https://doi.org/10.4168/aair.2017.9.6.483

Clinical Characteristics of Exacerbation-Prone Adult Asthmatics Identified by Cluster Analysis

Mi-Ae Kim1, , Seung-Woo Shin2, , Jong-Sook Park2, 3, Soo-Taek Uh4, Hun Soo Chang5, Da-Jeong Bae5, You Sook Cho6, Hae-Sim Park7, Ho Joo Yoon8, Byoung Whui Choi9, Yong Hoon Kim4, Choon-Sik Park3

1Department of Pulmonology, Allergy and Critical Care Medicine, CHA Bundang Medical Center, CHA University, Seongnam, Korea.

2Genome Research Center for Allergy and Respiratory Diseases, Soonchunhyang University Bucheon Hospital, Bucheon, Korea.

3Division of Allergy and Respiratory Medicine, Department of Internal Medicine, Soonchunhyang University Bucheon Hospital, Bucheon, Korea.

4Division of Allergy and Respiratory Medicine, Department of Internal Medicine, Soonchunhyang University Seoul Hospital, Seoul, Korea.

5Department of Interdisciplinary Program in Biomedical Science Major, Soonchunhyang Graduate School, Bucheon, Korea.

6Department of Allergy and Clinical Immunology, University of Ulsan College of Medicine, Asan Medical Center, Seoul, Korea.

7Department of Allergy and Clinical Immunology, Ajou University School of Medicine, Suwon, Korea.

8Department of Internal Medicine, Hanyang University College of Medicine, Seoul, Korea.

9Department of Internal Medicine, Chung-Ang University College of Medicine, Seoul, Korea.

Correspondence to: Choon-Sik Park, MD, PhD, Division of Allergy and Respiratory Medicine, Department of Internal Medicine, Soonchunhyang University Bucheon Hospital, 170 Jomaru-ro, Wonmi-gu, Bucheon 14584, Korea. Tel: +82-32-621-5105; Fax: +82-32-621-5023; mdcspark@hanmail.net

Equal contributor:

Mi-Ae Kim and Seung-Woo Shin contributed equally to this paper.

Received 25 January 2017       Revised 18 April 2017       Accepted 14 May 2017

Copyright © 2017 The Korean Academy of Asthma, Allergy and Clinical Immunology • The Korean Academy of Pediatric Allergy and Respiratory Disease

(open-access):

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Purpose

Asthma is a heterogeneous disease characterized by various types of airway inflammation and obstruction. Therefore, it is classified into several subphenotypes, such as early-onset atopic, obese non-eosinophilic, benign, and eosinophilic asthma, using cluster analysis. A number of asthmatics frequently experience exacerbation over a long-term follow-up period, but the exacerbation-prone subphenotype has rarely been evaluated by cluster analysis. This prompted us to identify clusters reflecting asthma exacerbation.

Methods

A uniform cluster analysis method was applied to 259 adult asthmatics who were regularly followed-up for over 1 year using 12 variables, selected on the basis of their contribution to asthma phenotypes. After clustering, clinical profiles and exacerbation rates during follow-up were compared among the clusters.

Results

Four subphenotypes were identified: cluster 1 was comprised of patients with early-onset atopic asthma with preserved lung function, cluster 2 late-onset non-atopic asthma with impaired lung function, cluster 3 early-onset atopic asthma with severely impaired lung function, and cluster 4 late-onset non-atopic asthma with well-preserved lung function. The patients in clusters 2 and 3 were identified as exacerbation-prone asthmatics, showing a higher risk of asthma exacerbation.

Conclusions

Two different phenotypes of exacerbation-prone asthma were identified among Korean asthmatics using cluster analysis; both were characterized by impaired lung function, but the age at asthma onset and atopic status were different between the two.

Keywords: Asthma, cluster analysis, exacerbation

INTRODUCTION

Asthma is a heterogeneous disease with various respiratory symptoms, such as wheezing, coughing, and dyspnea, and is characterized by chronic airway inflammation and remodeling.1 Owing to its heterogeneity, asthma has been classified into several phenotypes, such as atopic or non-atopic, eosinophilic or neutrophilic, and severe or non-severe subtypes.23 These classifications have been applied in the implementation of appropriate environmental management, prediction of prognosis, and treatment strategies.345678 Recently, candidate molecular targeted therapeutics, such as lebrikizumab and mepolizumab, have been used in cases with specific phenotypes, such as elevated immunoglobulin E (IgE) levels910 and blood or sputum eosinophilia.111213 Further subgroup analyses of asthmatics have also been initiated using newly adopted clinical, laboratory, and molecular markers.2141516
Cluster analysis involves grouping a set of objects together that are more similar to each other than to objects in other groups.17 For the first time, Haldar et al.18 demonstrated differences in clinical responses to treatments among clusters. In addition, in the Severe Asthma Research Program cohort, 2 severe asthma phenotype clusters were identified: 1) long disease duration with childhood onset and atopic disease; and 2) female predominance with mainly later-onset disease and less atopy.19 Further cluster analyses using additional parameters such as inflammatory biomarkers and obesity in asthmatics have also been reported.2021
Previously, we found differences in temporal changes in forced expiratory volume in 1 second (FEV1) over a 12-month follow-up among clusters of asthmatics.22 An additional longitudinal analysis of the populations in this cluster study demonstrated consistent differences in asthma control test scores and systemic corticosteroid usage among the clusters.23 Even the patients with refractory asthma could be classified into 4 subtypes.24 However, few cluster analyses have evaluated the risk of exacerbation in terms of long-term prognosis. In the present study, we attempted to determine whether clusters reflect the prognosis of risk of asthma exacerbation over a follow-up period of at least 1 year.

MATERIALS AND METHODS

Selection of the study subjects

Clinical and demographic data from 1,843 asthmatics registered in an asthma cohort in Korea22 were analyzed retrospectively. These asthmatics were recruited from 3 general hospitals in different areas. Asthma was diagnosed by physicians based on the Global Initiative for Asthma guidelines.1 Among the 1,843 subjects, 628 who were regularly followed-up for over 1 year were included in the analysis, after excluding current smokers and ex-smokers of 10 pack-year or more. The protocol was approved by the ethics committee of Soonchunhyang University Bucheon Hospital (SCHBC_2014_07_028).
At the baseline visit, demographic information was collected, such as age at enrollment, sex, body mass index (BMI), age at asthma onset, asthma duration, and smoking amount. All patients underwent a standardized assessment, which included peripheral blood total and differential cell counts, serum total IgE levels, chest radiography, sputum cell analysis and skin prick tests with 24 common inhalant allergens (Bencard Co., Brentford, UK). Baseline FEV1, forced vital capacity (FVC), short acting bronchodilator-induced percent change in FEV1 (ΔFEV1), and provocation concentration of methacholine inducing a 20% decrease in FEV1 (PC20) methacholine levels (mg/mL) were measured. The total amount of inhaled steroids was calculated from electronic medical records and expressed as the equivalent dosage of fluticasone propionate per day. Exacerbation of asthma was defined as episodes of increased dyspnea, wheezing, or coughing with a decrease in the personal best FEV1 by more than 20%.25 Short-term treatment with systemic corticosteroids was administered when the subjects experienced moderate exacerbation; the amount administered is expressed as the equivalent dosage of prednisone per year. Frequent exacerbation was defined when more than 2 episodes occurred. Sputum was induced using isotonic saline containing a short-acting bronchodilator and was processed within 2 hours of collection as described previously.4 Slides were prepared by cytocentrifugation, and 500 cells were counted on each slide after staining with Diff-Quick (American Scientific Products, Chicago, IL, USA).

Cluster analysis

Uniform cluster analysis was applied to each population using a 2-tiered approach.26 In the first step, hierarchical cluster analysis was performed using Ward's method to generate a dendrogram for estimation of the number of likely clusters within the population. Cluster boundaries were defined by large differences between successive fusion levels.27 K-means cluster analysis was then performed to select the most representative model. The K-means algorithm was repeated 10 times in each dataset at random starting points to ensure repeatability and stability within each model using SPSS version 18 (IBM Corp., Chicago, IL, USA). The variables used for cluster modeling were selected based on their contribution to characterization of the asthma phenotype. All measurements were standardized using z scores for continuous variables and 0 or 1 for categorical variables. Initially, we selected 12 variables (age at enrollment, sex, BMI, age at asthma onset, asthma duration, smoking amount, atopy, peripheral blood eosinophils (%), total IgE level, FEV1, ΔFEV1, and PC20) that have been identified as determinants of asthma.2829

Statistical analysis

SPSS version 18 (IBM Corp.) was used for the principal components analysis (PCA), clustering, and discriminant analyses. To compare clinical parameters among clusters, analysis of variance, Kruskal-Wallis test, and χ2 test were applied to parametric continuous, non-parametric continuous, and categorical variables, respectively. In addition, the times of exacerbation appearance were compared among the clusters using Kaplan-Meier analysis and the log-rank test.

RESULTS

Clinical characteristics of the study population

A total of 628 subjects were regularly followed-up for more than 1 year. Among them, 109 subjects with missing PC20 data and 260 subjects who presented with acute exacerbation at the initial visit were excluded; the remaining 259 subjects were analyzed for clustering. Clinical and demographic characteristics are summarized in Table 1.

Hierarchical cluster analysis

To estimate the number of likely clusters, hierarchical analysis using Ward's method was performed to generate dendrograms using 12 variables (Supplementary Figs. 1 and 2) and showed that the subjects could be classified into 4 clusters. The hierarchical analysis was repeated using the 12 variables, revealing that 4 clusters are likely optimal for discriminating the subjects; thus, subsequent K-means cluster analysis was performed using the 4 clusters. The heap map and hierarchical clustering of these 12 variables are presented in Fig. 1.

Characteristics of the 4 clusters

Age, age at asthma onset, BMI, atopic status, FEV1, ΔFEV1, and PC20 were markedly different among the 4 clusters (Table 1). Cluster 1 (C1) comprised the youngest patient group with early-onset asthma. Individuals in this cluster had well-preserved FEV1 with a high rate of positive skin prick tests. Thus, C1 was defined as early-onset atopic asthma with preserved lung function. Cluster 2 (C2) comprised the oldest patient group with an age at onset of 56 years. These individuals had the lowest rate of positive skin prick tests and decreased FEV1. Thus, C2 was defined as late-onset non-atopic asthma with impaired lung function. Cluster 3 (C3) showed the lowest FEV1 and PC20 with the highest rate of positive skin prick tests, while cluster 4 (C4) comprised older patients with the highest FEV1 and a lower rate of atopy. Thus, C3 was defined as early-onset atopic asthma with severely impaired lung function and C4 as late-onset non-atopic asthma with well-preserved lung function. The age at asthma onset was different among age groups; early-onset asthma was found in the younger groups (C1 and C3), while late-onset asthma was found in the older groups (C2 and C4). Atopic status was distinguishable between age groups, showing higher rates in the early-onset groups (C1 and C3) than in the late-onset groups (C2 and C4).
We also performed a post-hoc analysis (Table 1) of variables with P values <0.05. All pairs of clusters differed significantly in terms of age and FEV1. C1 represented the youngest group and C2 the oldest. C3 showed the worst FEV1 while C4 showed the best. In early-onset asthmatics (C1 and C3), the distinguishable parameters were FEV1 and PC20, similar variables to the late-onset groups (C2 and C4), FEV1, and ΔFEV1.

Use of corticosteroids and exacerbation frequency among the clusters

The dosages of inhaled steroids, number of subjects experiencing exacerbation, and the annual number of exacerbations showed inter-individual differences (Table 2). More than 80% of the subjects in each cluster used inhaled steroid. The dosage of inhaled steroid was lower in C1 than other clusters (Fig. 2A), while the proportion of subjects using high-dose inhaled steroid was not significantly different. There was no significant difference in the proportion or dosage of systemic steroid among clusters (Fig. 2B).
The exacerbation rate over 1 year was significantly higher in C2 and C3 than in C1 and C4 (Table 2). The difference was still apparent in a post-hoc analysis. Furthermore, the annual number of exacerbations was higher in C2 and C3 than in C1 and C4 (Fig. 2C). The proportion of subjects with frequent exacerbations (more than 2 times per year) was higher in C2 and C3 than in the other clusters (Table 2). Thus, we defined C2 and C3 as the exacerbation-prone clusters.
Additionally, correlations between the annual rate of exacerbation and the dosage of inhaled or systemic corticosteroids used were analyzed. The annual rate of exacerbation was significantly correlated with the dosage of inhaled corticosteroids in all study subjects (Supplementary Fig. 3).

Comparison of exacerbation risk

To compare the exacerbation risk among the 4 clusters, Kaplan-Meier analysis, the log-rank test, and Cox regression analyses were applied. The cumulative rate of the first and second attack was different among the 4 clusters in the Kaplan-Meier analysis (Fig. 3). C3 showed a significantly increased risk of the first exacerbation compared with C1 (odds ratio [OR], 17.569; P<0.001) and C4 (OR, 12.664; P<0.001). C3 also showed an increased risk of a second exacerbation compared with C1 and C4. No significant differences in the first and second exacerbation were observed between C2 and C3 (Supplementary Table 1).
Next, Cox regression analysis was performed to compare the risk of exacerbation among clusters after adjusting for other clinical parameters. To identify multicollinearity, correlation analysis was performed using 12 variables (Supplementary Table 2). The result of the Cox regression analysis using these 12 variables and cluster groups are presented in Supplementary Table 3. Asthma duration and FEV1 were significantly associated with exacerbation, while no significant associations were found in any of the clusters.
To reveal the differences in inflammatory patterns among the clusters, total and differential cell counts in peripheral blood (n=259) and sputum (n=86) were compared. There were no differences in the cellular fraction or total cell counts among the clusters, except for white blood cell and neutrophil counts in peripheral blood (Supplementary Table 4).

DISCUSSION

In the present study, we performed a longitudinal analysis of the clinical outcomes of asthma subphenotypes defined by cluster analysis. We included only patients who were followed-up for longer than 1 year to analyze the annual exacerbation rate and exacerbation risk within 1 year. Four clusters among patients with adult asthma were identified in this study. Among them, C2 and C3 comprised patients with lower lung function and higher exacerbation risk. We defined these clusters as exacerbation-prone asthma. The C2 and C3 clusters showed contrasting clinical characteristics; C2 showed late-onset (age 50s) and non-atopic asthma, while C3 showed early-onset (age 20s) and atopic asthma. To the best of our knowledge, this study is the first to analyze the exacerbation risk of asthma clusters using 1-year follow-up data.
The clinical characteristics of the exacerbation-prone clusters in our study differed somewhat from those in other studies. First, obesity and the inflammatory patterns in the airway may not have been factors that aggravated asthma in our clusters. A recent cohort study by Amelink et al.21 surveying 200 patients with adult-onset asthma showed the following exacerbation-prone clusters: 1) obese females with frequent symptoms, high healthcare utilization, and low sputum eosinophil counts; and 2) a severe eosinophilic inflammation-predominant group with persistent airflow limitation. Obese asthma patients showed a different expression profile of clinical and inflammatory biomarkers and reduced glucocorticoid receptor expression, which may mediate glucocorticoid insensitivity.20 However, the exacerbation-prone subjects in our study did not have a higher BMI or eosinophilic inflammation, which is in contrast to studies conducted in Caucasian populations.18192021 We decided that BMI was not a significant parameter to define asthma clusters, because Asian asthmatics usually have a normal BMI range compared with Western populations. Few Asian asthmatics have BMIs over 30 kg/m2, and all of our subjects were in the range of 21.9–26.2 kg/m2. Second, the asthma duration was not significantly different among clusters. Instead, there were 2 age groups among the 4 clusters: older subjects with an age at asthma onset of over 50 years in C2 and C4 and younger subjects with an age at asthma onset in their 20s in C1 and C3. In other studies, asthmatics were in their late 30s to late 40s and had a long asthma duration of approximately 30 years in Moore et al.'s study,19 or they were in their mid- to late 40s and had early-onset asthma with a duration of approximately 25 years in the study by Haldar et al.'s study.18 We believe these differences originated from the relatively short asthma duration (4.7–10.0 years) of our cohort.
When comparing the present results with those of a previous cluster analysis performed in South Korea,2224 there are apparent differences in the demographic parameters of study populations and variables used for clustering. In the present study, we excluded current smokers and ex-smokers of 10 pack-year or more and patients presented with acute exacerbation at the initial visit while the previous study did not. The main result of that study22 is the presence of 4 asthma subtypes: 1) smoking asthma; 2) severe obstructive asthma; 3) early-onset atopic asthma; and 4) late-onset mild asthma using 6 variables (FEV1, BMI, age at onset, atopic status, smoking history, and history of hospital use due to exacerbation) in 724 asthmatics of the Cohort for Reality and Evolution of Adult Asthma in Korea (COREA) and 1,843 asthmatics of Soonchunhyang University Asthma Genome Research Centre (SCH) cohorts. However, the present study analyzed a complete set of 12 variables and exacerbation frequency with asthma control medications' usage for 1 year follow-up period. As the result, only 259 patients are included in the present analysis among the 1,843 subjects of SCH cohort. In another cluster analysis,24 86 refractory asthma patients were selected from SCH cohort on the basis of the American Thoracic Society's criteria.30 This study used 5 variables (age at the onset of asthma, BMI, FEV1/FVC ratio, amount of cigarette smoking, and PC20) in the cluster analysis and mean smoking amount was 8.4 pack-year, while the present study excluded current smokers and ex-smokers of more than 10 pack-year. In addition, only 39 subjects were included in both 2 studies. Accordingly, the present study seems to present new cluster, C2, which the previous studies of ours could not, on the basis of exacerbation frequency and asthma medication usages. C2 comprised the oldest patient group with an age at onset of 56 years. These individuals had the lowest rate of positive skin prick tests and decreased FEV1. Thus, C2 was defined as late-onset non-atopic asthma with impaired lung function.
Our study has several limitations caused by the difficulty in collecting clear follow-up data in practice. First, we did not record adherence at clinic visits, making it difficult to assess whether the study subjects visited a hospital for every occurrence of exacerbation. Second, we did not check the remaining dosage of inhaled steroids after use. Because some patients did not use steroids even when they had been prescribed them, the dosages of inhaled steroids calculated in our study might be overestimates. Third, the subjects were recruited from only 3 hospitals (Seoul, Bucheon, and Cheonan hospital), which may have contributed to a selection bias. Fourth, there is a chance that mild asthma was excluded, because only patients followed-up over 1 year were included in the analysis. Fifth, the proportion of female patients was higher than that in other cohorts, so we could not analyze the effects of sex on asthma exacerbation. However, we believe that our study was valuable in terms of the racial differences and the data completeness including the exact exacerbation-free duration in all subjects. We showed different clinical characteristics among exacerbation-prone asthmatics in an Asian population and analyzed complete exacerbation data from all subjects included in the cluster analysis. Further studies should be performed to define the parameters or pathogeneses that contributed to the asthma exacerbation observed in patients in C2 and C3, in addition to FEV1.
In conclusion, 2 different phenotypes of exacerbation-prone asthma were identified in Korean asthmatics using cluster analysis, revealing somewhat different clinical characteristics from those of asthmatics in Western countries. Both of the phenotypes showed impaired lung function with different ages at asthma onset and atopic statuses, while no significant differences in BMI or eosinophilic inflammation were identified.

Figures and Tables

Fig. 1

A heap map. Hierarchical cluster analysis of 259 subjects using 12 parameters: age at enrollment, age at asthma onset, BMI, atopy, log total IgE, duration of asthma, ΔFEV1, FEV1 (% predicted), log bood eosinophils (%), sex, smoking history, and PC20. Red color indicates 100% and green color indicates 0%. The gradations of the 2 colors indicate the range from 0% to 100%. BMI, body mass index; IgE, immunoglobulin E; FEV1, forced expiratory volume in 1 second; ΔFEV1, short acting bronchodilator-induced percentage increase in FEV1; PC20, provocation concentration of methacholine inducing a 20% decrease in FEV1.

aair-9-483-g001
Fig. 2

Comparison of the corticosteroid dose and exacerbation frequency during follow-up. The dosages of inhaled corticosteroids (A) and systemic corticosteroids (B) are presented as fluticasone equivalents (µg/day) and prednisone equivalents (mg/year), respectively. The annual rate of exacerbation (C) was calculated as the total number of exacerbations divided by the follow-up duration (years). The subject numbers were 72 (C1), 48 (C2), 52 (C3), and 87 (C4). Box plots show the median values with interquartile ranges. P values were obtained using the Mann-Whitney U test. C1, cluster 1; C2, cluster 2; C3, cluster 3; C4, cluster 4.

aair-9-483-g002
Fig. 3

Kaplan-Meier plots of the accumulation risk of the first (A) and second (B) exacerbations. Y-axis, accumulation rate of exacerbation; X-axis, time to exacerbation. P values were obtained using the χ2 test.

aair-9-483-g003
Table 1

Clinical and demographic characteristics of the study subjects

aair-9-483-i001
Variable Total C1 C2 C3 C4 P value
Number of subjects 259 72 48 52 87
Age (year) 56 (18-88) 40 (18–71)*,†,‡ 68 (52–88)§,II 47 (22–79) 61 (46–84) 1.57E-42
Age at asthma onset (year) 41 (1–77) 27 (1–43)*,‡ 56 (37–74)§,II 28 (5–46) 51 (32–77) 3.28E-39
Asthma duration (year) 7.3 (1.1–25.2) 7.1 (1.1–16.7) 6.7 (1.1–19.8) 7.8 (1.2–17.1) 7.4 (1.2–25.2) 0.992
Female (%) 211 (81.5%) 52 (72.2%) 40 (83.3%) 42 (80.8%) 77 (88.5%) 0.07
BMI (kg/m2) 24.1 (15.9–35) 22.5 (16–33.8)*,‡ 24.2 (18.3–35) 23.5 (15.9–33.1) 25 (19.9–33.7) 1.95E-04
Never-smoker/ex-smoker 238/21 (8.1%) 63/9 (12.5%) 43/5 (10.4%) 49/3 (5.8%) 83/4 (4.6%) 0.258
Smoking amount (pack-year) 0.0 (0–9) 0 (0–9) 0 (0–6) 0 (0–3) 0 (0–5) 0.213
Skin prick test (%) 131 (50.8%) 44 (62.0%)*,‡ 18 (37.5%)§ 34 (65.4%) 35 (39.8%) 0.002
Total IgE (IU/L) 151 (2–3564) 162 (4–3564) 124.5 (4–3182) 164 (10–2374) 136 (2–1734) 0.855
FEV1 (% predicted) 77 (18–140) 83 (65–140)*,†,‡ 57.5 (18–71)§,II 52.5 (24–68) 87 (73–116) 1.04E-39
ΔFEV1 (%) 7 (0–55) 6.5 (0–39)* 12 (0–55)II 10.5 (0–45) 6 (0–34) 0.015
PC20 (mg/mL) 1.1 (0–25) 0.9 (0–25)† 1.4 (0–25)§ 0.3 (0–25) 1.5 (0–25) 0.005
Blood eosinophils (%) 4.8 (0.0–38.7) 5.2 (0.1–22.9) 4.9 (0–36) 3.7 (0.1–27.7) 5 (0.1–38.7) 0.136

BMI: body mass index, FEV1: forced expiratory volume in 1 second, ΔFEV1: short acting bronchodilator-induced percentage increase in FEV1, PC20: provocation concentration of methacholine inducing a 20% decrease in FEV1. All values are presented as medians with ranges. The normality of the distribution was evaluated using the Shapiro-Wilk test. The Kruskal-Wallis test was used to compare variables among the groups, and P values are presented in the right hand column. Post hoc analyses using Mann-Whitney U tests were performed for comparisons of two groups. The chi-squared test was used to compare qualitative variables. *C1 vs C2 (P<0.05); †C1 vs C3 (P<0.05); ‡C1 vs C4 (p<0.05); §C2 vs C3 (p<0.05); ∥C2 vs C4 (P<0.05); C3 vs C4 (P<0.05).

Table 2

Dosages of the corticosteroids used and exacerbation rates in the study subjects

aair-9-483-i002
Parameter Total C1 C2 C3 C4 P value
Subject numbers 259 72 48 52 87
Subjects using inhaled steroids (%) 214 (82.6) 61 (84.7) 42 (87.5) 39 (75.0) 72 (82.8) 0.373
Dosage of inhaled steroids (µg/day)* 441.7 (20.5–2630.1) 287.6 (41–863)§,II 493.1 (20.5–1397.2)†† 494.1 (82.1–1109.5) 431.5 (41–2630.1) 0.002
 0–249 55 (21.2) 25 (34.7)§,II,¶ 3 (6.3)**,†† 10 (19.2) 17 (19.5) 0.002
 250–499 81 (31.2) 21 (29.2) 19 (39.6) 12 (23.1) 29 (33.3) 0.322
 ≥500 78 (30.1) 15 (20.8) 20 (41.7) 17 (32.7) 26 (29.9) 0.104
0.009
Number of subjects using systemic steroids (%) 67 (25.8) 256 (30-4934) 17 (35.4) 17 (32.7) 20 (23.0) 0.102
Dosage of systemic steroids (mg/year)† 441 (8–4934) 5 (45.4) 280 (90–3885) 530 (8–4340) 497 (70–4645) 0.736
 90–360 26 (10.3) 2 (18.1) 9 (52.9) 4 (25) 8 (42.1) 0.355
 361–999 22 (8.4) 4 (36.3) 5 (29.4) 8 (50) 7 (36.8) 0.110
≥1,000 15 (5.7) 3 (17.6) 4 (25) 4 (21) 0.120
0.570
Exacerbation rate over 1 year (%) 143 (55.2) 34 (47.2)§,II 37 (77.1)†† 39 (75.0)‡‡ 33 (37.9) 1.00E-06
Annual number of exacerbations‡ 1 (1-11) 1 (1–8)§,II 2 (1–10)†† 2 (1–11)‡‡ 1 (1–5) 5.30E-05
Subjects having exacerbations
 0 116 (44.7) 38 (52.7)§,II 11 (22.9)†† 13 (25)‡‡ 54 (62) 1.00E-6
 1 81 (31.2) 29 (40.27) 15 (31.3) 14 (26.9) 23 (26.4) 0.249
 2 41 (15.8) 3 (4.16)§,II 12 (25.0)†† 19 (36.5)‡‡ 7 (8.0) 8.21E-7
 ≥3 21 (8.1) 2 (2.77)§,II 10 (20.8)†† 6 (11.5) 3 (3.4) 8.81E-4
9.91E-10

*Dosages of inhaled corticosteroid are presented as fluticasone equivalents (µg/day), calculated as the total amount of inhaled corticosteroids during the first year after the initial visit. †Dosages of systemic corticosteroid are presented as prednisone equivalents (mg/year), calculated as the total amount of prednisone during the first year after the initial visit, *, †, ‡ All values are presented as medians with ranges. Other data are presented as the numbers of subjects (% of the total subjects or subtotal in each cluster). The Kruskal-Wallis test was used to compare the variables among groups, and P values are presented in the right hand columns. Post hoc analysis using the Mann-Whitney U test was performed to compare two groups. §C1 vs C2 (P<0.05); ∥C1 vs C3 (P<0.05); C1 vs C4 (P<0.05); **C2 vs C3 (P<0.05); ††C2 vs C4 (P<0.05); ‡‡C3 vs C4 (P<0.05).

ACKNOWLEDGMENTS

The study was supported by a research grant of Soonchunhyang University to Choon-Sik Park and by a grant of the Korean Health Technology Research and Development Project, Ministry of Health and Welfare, Republic of Korea (2014-ER7401-00). Clinical and laboratory data are provided from a biobank of Soonchunhyang University Bucheon Hospital.

Notes

There are no financial or other issues that might lead to conflict of interest.

References

1. Global Initiative for Asthma. Global strategy for asthma management and prevention [Internet]. [place unknown]: Global Initiative for Asthma;2016. 2017 January 25. Available from: http://www.ginasthma.org/.
2. Kim MA, Shin YS, Pham le D, Park HS. Adult asthma biomarkers. Curr Opin Allergy Clin Immunol. 2014; 14:49–54.
3. Lemière C, Ernst P, Olivenstein R, Yamauchi Y, Govindaraju K, Ludwig MS, et al. Airway inflammation assessed by invasive and noninvasive means in severe asthma: eosinophilic and noneosinophilic phenotypes. J Allergy Clin Immunol. 2006; 118:1033–1039.
4. Choi JS, Jang AS, Park JS, Park SW, Paik SH, Park JS, et al. Role of neutrophils in persistent airway obstruction due to refractory asthma. Respirology. 2012; 17:322–329.
5. Jang AS, Lee JH, Park SW, Park JS, Kim DJ, Park CS. Risk factors related to fixed airway obstruction in patients with asthma after anti-asthma treatment. Ann Allergy Asthma Immunol. 2007; 99:408–412.
6. Park CS. Eosinophilic bronchitis, eosinophilia associated genetic variants, and notch signaling in asthma. Allergy Asthma Immunol Res. 2010; 2:188–194.
7. Haldar P, Pavord ID. Noneosinophilic asthma: a distinct clinical and pathologic phenotype. J Allergy Clin Immunol. 2007; 119:1043–1052.
8. Fajt ML, Wenzel SE. Development of new therapies for severe asthma. Allergy Asthma Immunol Res. 2017; 9:3–14.
9. Bousquet J, Cabrera P, Berkman N, Buhl R, Holgate S, Wenzel S, et al. The effect of treatment with omalizumab, an anti-IgE antibody, on asthma exacerbations and emergency medical visits in patients with severe persistent asthma. Allergy. 2005; 60:302–308.
10. Bousquet J, Siergiejko Z, Swiebocka E, Humbert M, Rabe KF, Smith N, et al. Persistency of response to omalizumab therapy in severe allergic (IgE-mediated) asthma. Allergy. 2011; 66:671–678.
11. Ortega HG, Liu MC, Pavord ID, Brusselle GG, FitzGerald JM, Chetta A, et al. Mepolizumab treatment in patients with severe eosinophilic asthma. N Engl J Med. 2014; 371:1198–1207.
12. Corren J, Lemanske RF Jr, Hanania NA, Korenblat PE, Parsey MV, Arron JR, et al. Lebrikizumab treatment in adults with asthma. N Engl J Med. 2011; 365:1088–1098.
13. Nair P, O’Byrne PM. Measuring eosinophils to make treatment decisions in asthma. Chest. 2016; 150:485–487.
14. Hanania NA, Wenzel S, Rosén K, Hsieh HJ, Mosesova S, Choy DF, et al. Exploring the effects of omalizumab in allergic asthma: an analysis of biomarkers in the EXTRA study. Am J Respir Crit Care Med. 2013; 187:804–811.
15. Izuhara K, Ohta S, Ono J. Using Periostin as a biomarker in the treatment of asthma. Allergy Asthma Immunol Res. 2016; 8:491–498.
16. Wardzyńska A, Makowska JS, Pawełczyk M, Piechota-Polańczyk A, Kurowski M, Kowalski ML. Periostin in exhaled breath condensate and in serum of asthmatic patients: relationship to upper and lower airway disease. Allergy Asthma Immunol Res. 2017; 9:126–132.
17. Wikipedia. Cluster analysis [Internet]. Los Angeles (CA): Wikipedia;2017. cited year month day. Available from: http://en.wikipedia.org/wiki/Cluster_analysis.
18. Haldar P, Pavord ID, Shaw DE, Berry MA, Thomas M, Brightling CE, et al. Cluster analysis and clinical asthma phenotypes. Am J Respir Crit Care Med. 2008; 178:218–224.
19. 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–323.
20. Sutherland ER, Goleva E, King TS, Lehman E, Stevens AD, Jackson LP, et al. Cluster analysis of obesity and asthma phenotypes. PLoS One. 2012; 7:e36631.
21. Amelink M, de Nijs SB, de Groot JC, van Tilburg PM, van Spiegel PI, Krouwels FH, et al. Three phenotypes of adult-onset asthma. Allergy. 2013; 68:674–680.
22. Kim TB, Jang AS, Kwon HS, Park JS, Chang YS, Cho SH, et al. Identification of asthma clusters in two independent Korean adult asthma cohorts. Eur Respir J. 2013; 41:1308–1314.
23. Park SY, Baek S, Kim S, Yoon SY, Kwon HS, Chang YS, et al. Clinical significance of asthma clusters by longitudinal analysis in Korean asthma cohort. PLoS One. 2013; 8:e83540.
24. Jang AS, Kwon HS, Cho YS, Bae YJ, Kim TB, Park JS, et al. Identification of subtypes of refractory asthma in Korean patients by cluster analysis. Lung. 2013; 191:87–93.
25. Kim S, Kim Y, Lee MR, Kim J, Jung A, Park JS, et al. Winter season temperature drops and sulfur dioxide levels affect on exacerbation of refractory asthma in South Korea: a time-trend controlled case-crossover study using Soonchunhyang asthma cohort data. J Asthma. 2012; 49:679–687.
26. Ball GH, Hall DJ. A clustering technique for summarizing multivariate data. Behav Sci. 1967; 12:153–155.
27. Everitt BS. Cluster analysis. 3rd ed. New York (NY): John Wiley;1993.
28. Rosi E, Ronchi MC, Grazzini M, Duranti R, Scano G. Sputum analysis, bronchial hyperresponsiveness, and airway function in asthma: results of a factor analysis. J Allergy Clin Immunol. 1999; 103:232–237.
29. Miranda C, Busacker A, Balzar S, Trudeau J, Wenzel SE. Distinguishing severe asthma phenotypes: role of age at onset and eosinophilic inflammation. J Allergy Clin Immunol. 2004; 113:101–108.
30. Proceedings of the ATS workshop on refractory asthma: current understanding, recommendations, and unanswered questions. American Thoracic Society. Am J Respir Crit Care Med. 2000; 162:2341–2351.

SUPPLEMENTARY MATERIALS

Supplementary Table S1

Pair-wise comparisons of the first (A) and the second (B) exacerbations among the four clusters

Supplementary Table S2

Correlation coefficients between the variables.

Supplementary Table S3

Cox's regression analysis of variables from the first attack (A) to the second attack (B)

Supplementary Table S4

Cellular profiles of the peripheral blood and sputum in the four clusters

Supplementary Fig. 1

Dendrogram obtained using the 12 parameters.

Supplementary Fig. 2

Hierarchical clustering using the 12 parameters.

Supplementary Fig. 3

Correlations between the total number of exacerbations and corticosteroid dosage. (A) Inhaled corticosteroids (Spearman's rho=0.197, P=0.019) and (B) systemic corticosteroids (Spearman's rho=0.034, P=0.684).
TOOLS
ORCID iDs

You Sook Cho
https://orcid.org/0000-0001-8767-2667

Hae-Sim Park
https://orcid.org/0000-0003-2614-0303

Choon-Sik Park
https://orcid.org/0000-0003-2977-0255

Similar articles