Skip to content


You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Pulmonary Medicine

Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Persistent asthma phenotype related with late-onset, high atopy, and low socioeconomic status in school-aged Korean children

  • Eun Lee1,
  • Si Hyeon Lee2,
  • Ji-Won Kwon3,
  • Young-Ho Kim4,
  • Jisun Yoon5,
  • Hyun-Ju Cho5,
  • Song-I Yang6,
  • Young-Ho Jung7,
  • Hyung Young Kim8,
  • Ju-Hee Seo9,
  • Hyo Bin Kim10,
  • So Yeon Lee5,
  • Ho-Jang Kwon11 and
  • Soo-Jong Hong5Email author
BMC Pulmonary MedicineBMC series – open, inclusive and trusted201717:45

Received: 19 October 2016

Accepted: 16 February 2017

Published: 23 February 2017



Treatment guidelines for asthma have been established based on asthma severity; there are limitations in the identification of underlying pathophysiology and prediction of prognosis in heterogeneous phenotypes of asthma. Although the complex interactions between environmental and genetic factors affect the development and progression of asthma, studies on asthma phenotypes considering environmental factors are limited. This study aimed to identify asthma phenotypes using latent class analysis including environmental factors in school-age children.


We included 235 children (6–8 years) with parent-reported, physician-diagnosed asthma from the Children’s HEalth and Environmental Research (CHEER) study, which is a 4-year prospective follow-up study with 2-year intervals. At every survey, pulmonary function tests, methacholine challenge tests and blood tests with questionnaire were conducted.


Four asthma phenotypes were identified. Cluster 1 (22% of children) was characterized by high prevalence of atopy and mild symptoms; subjects in cluster 2 (17%) consisted of less atopy and normal lung function, but intermittent troublesome; cluster 3 (29%) experienced late-onset atopic troublesome asthma with decreased lung function in combination with low socioeconomic status; and cluster 4 was associated with early-onset and less-atopic infrequent asthma.


Late-onset, high atopy, and low socioeconomic status are associated with troublesome persistent asthma phenotype in school-age children. Environmental factors might be implicated in the clinical heterogeneity of asthma. Asthma phenotypes considering diverse factors might be more helpful in the identification of asthma pathogenesis and its prevention.


AsthmaAtopyChildrenPhenotypeSocioeconomic status


Asthma is characterized by airway hyperresponsiveness and reversible airway obstruction, which results from complex interactions between genetic and environmental factors [1]. Wheezing in early life is usually associated with viral respiratory tract infection; however, frequent recurrence in these groups causes a considerable healthcare burden [2]. Owing to the clinical and pathophysiological heterogeneity of asthma, there have been diverse attempts to phenotype asthma [37]. However, the actual clinical application of phenotyping asthma was based on symptom severity [8]. Classification of asthma is important to predict the prognosis and response to treatment and understand its underlying mechanisms; however, studies on asthma phenotypes using comprehensive factors are limited.

Epidemiologic studies have considered instability of asthma phenotypes over time [3, 4, 7] and identified associated risk factors for each asthma phenotype [9]. However, such studies considering environmental factors for the classification of asthma phenotypes are limited, although environmental factors play an important role as causative and provocative factors for asthma [10]. Additionally, few studies artificially classified asthma phenotypes according to onset age and/or persistence of wheezing, despite diverse factors affecting the development and progression of asthma [9, 11]. Environmental factors affect the development of airway during antenatal and/or postnatal period through epigenetic or non-epigenetic mechanisms [12, 13]. Although persistent troublesome wheeze is known to be associated with atopy, bronchial hyperresonsiveness (BHR), and decreased pulmonary function compared with controlled asthma [6, 14], these studies described phenotypes of preschool wheezing based on temporal patterns of symptoms combined with symptom control. Another study showed that sex is a risk factor for different asthma phenotypes [15]. However, these studies have limitations in identification of the underlying pathology according to phenotypes. No study has shown how environmental factors combined with genetic factors reflected in parental history of allergic diseases and demographic factors are associated with the heterogeneity of asthma, which would be helpful to investigate the origin of asthma.

The aim of this study was to identify asthma phenotypes based on environmental factors and clinical characteristics of asthma in school-age children using latent class analysis (LCA), which is an unsupervised statistical method for classification [16]. Additionally, we investigated the prognosis of each asthma phenotype and changes in biomarkers during the prospective follow-ups.


Study population

This study was a part of Children’s HEalth and Environmental Research (CHEER) study, which aimed to evaluate the effect of environmental factors on children’s health [1618]. The study population was enrolled from January 2005 to December 2006, and they were followed-up for 4 years with 2-year intervals until December 2010 in randomly selected nationwide elementary schools in Korea [16, 17]. This cohort study consisted of parent-completed questionnaire, skin prick tests (SPTs), blood tests for total serum IgE levels and eosinophil percentages, pulmonary function tests, and methacholine challenge tests at every survey. The Korean version of the International Study of Asthma and Allergies in Childhood (ISAAC) questionnaire was used to assess the presence of allergic diseases [1921]. Out of the total 2491 children at the time of enrollment in CHEER study, those (mean age ± standard deviation, 7.5 ± 0.9 years) with parent-reported, physician-diagnosed asthma in lifetime were included in the present study. The study protocols were reviewed and approved by the institutional review board (IRB) of Ulsan University College of Medicine (IRB no. 2006–0081). Written informed consent was obtained from the parents or guardians of the study participants Parents or guardians completed the questionnaire.

Skin prick tests (SPTs)

The SPTs were performed for 13 common inhalant allergens and 4 food allergens [17]. Atopy was defined as at least one positive SPT response.

Measurement of total serum IgE levels and blood eosinophil percentages

Total serum IgE levels were measured by using fluorescent enzyme immunoassay (AutoCAP system, Phadia AB, Uppsala, Sweden). The lower limit of total serum IgE detection was 2 kUA/L. Blood eosinophil percentages were evaluated using NE-8000 system (Sysmex, Kobe, Japan) from peripheral blood samples.

Lung function test and BHR

Spirometry was performed according to the American Thoracic Society (ATS)/European Respiratory Society (ERS) guidelines by using a portable microspirometer (Microspiro HI-298, Chest Corporation, Tokyo, Japan) [17, 22, 23]. Details of BHR measurement are provided in elsewhere [24]. Subjects were considered to have BHR to methacholine when their methacholine PC20 was <8 mg/mL. The new development of BHR was defined as conversion from a negative BHR at enrollment to a positive BHR at the 4-year follow-up survey.

Cluster analysis

Asthma phenotypes were identified by using LCA. LCA has been employed previously to identify phenotypes of asthma, atopic dermatitis (AD), and allergic rhinitis (AR) [16, 25, 26] because this statistical method can classify distinctly based on the observed heterogeneity in total population without bias. Several variables were used to define asthma phenotypes: demographic characteristics, namely, sex, body mass index (BMI), and socioeconomic status (SES) using monthly house hold income ( and maternal education level; and environmental factors, namely, prior exposure to environmental tobacco smoke (ETS), parental history of allergic diseases (AD, AR, and/or asthma), onset age of the first wheezing, number of absence days from school and nocturnal awakening days owing to asthma attacks, number of asthma attacks in the previous 12 months, number of days with treatment for asthma in the previous 12 months at the time of enrollment, and SPT, blood tests, and lung function tests at the time of enrollment. The FEV1 % predicted and FEV1/FVC were dichotomized as ≥80% and <80% and for FEF25–75% % predicted as ≥65% and <65%. Based on the Bayesian information criteria, a model with four latent classes was selected as the best model for CHEER data [27].

Statistical analysis

To evaluate the differences in characteristics and parameters among the four asthma phenotypes, we used the t test for independent comparisons between two groups, Fisher exact test for categorical variables, or Kruskal–Wallis test for continuous variables. Multiple testing was conducted using the Bonferroni correction when needed. P < 0.05 was considered significant. All statistical analyses were performed using SAS for Windows (version 9.2).


Characteristics of study population

The characteristics of the study population are shown in Table 1. Mean age of the enrolled children was 7.5 years, and 63.4% of the study population were boys. The prevalence of parental history of allergic diseases, atopy, and physician-diagnosed AR and AD in lifetime was significantly higher in the included subjects compared to those excluded due to the absence of physician-diagnosed asthma history.
Table 1

Characteristics of study population


Participants in 4 clusters

Excluded in LCA

P value


235/2491 (9.4%)

2256/2491 (90.6%)



149/235 (63.4%)

1148/2256 (50.9%)


Age (years), mean ± SD

7.5 ± 0.9

7.7 ± 1.1


BMI (kg/m2), mean ± SD

16.95 ± 2.4

16.87 ± 2.6


Parental allergic diseases, yes

97/227 (42.7%)

582/2114 (27.5%)


Maternal educational levels, ≤ High school

133/222 (59.9%)

1261/1967 (64.1%)


Parental monthly income, < 3000 USD

152/230 (66.1%)

1501/2132 (70.4%)


ETS exposure, yes

106/227 (46.7%)

910/2004 (45.4%)


Atopy on skin prick tests

104/234 (44.4%)

583/2254 (25.9%)


Total serum IgE (IU/ml), mean [range]

454 [3.80–4630]

250.80 [1.40–6934]


Blood eosinophils (%), mean [range]

4.67 [0.50–39.10]

3.57 [0–22.60]


AR diagnosis, ever

107/234 (45.7%)

462/2081 (22.2%)


AD diagnosis, ever

99/233 (42.5%)

642/2075 (30.9%)


FEV1 % predicted

102.4 ± 13.77

106.1 ± 12.96



88.86 ± 6.06

90.86 ± 5.54


FEF25–75% % predicted

91.92 ± 21.88

99.91 ± 22.64


BHR (PC20 < 8 mg/ml)

80/221 (36.2%)

408/2204 (18.5%)


Data are presented as number (%) or mean ± standard deviation

Definition of abbreviations: AD atopic dermatitis, AR allergic rhinitis, BHR methacholine PC20 < 8mg/ml, BMI, body mass index, ETS environmental tobacco smoke, FEF 25–75% Forced expiratory flow at 25–75% of forced vital capacity, FEV 1 Forced expiratory volume in 1 s, FVC Forced vital capacity, LCA latent class analsis, SD standard deviation, USD US dollars

Cluster 1 (atopic mild asthma)

Fifty-one children (21.7%) were grouped into cluster 1, which was characterized by atopic mild asthma (Table 2). Atopy was observed in 46 children (90.2% of the children in this cluster). Onset age of the first wheezing in this cluster ranged widely from ≤ 1 year to > 6 years of age. None of the children in this cluster complained of absence days from school and nocturnal awakening days due to asthma attacks in the previous 12 months at the time of enrollment. Forty-two children (82.4%) had no asthma attacks in the previous 12 months at the time of enrollment. However, BHR was observed for 44.9% of the children (n = 22/49) and > 100 IU/mL of total serum IgE levels for 98.0% (n = 50/51).
Table 2

Results of latent class analysis


Cluster 1

Cluster 2

Cluster 3

Cluster 4

Number, n (%)

51 (21.7%)

40 (17.0%)

29 (12.3%)

115 (49.1%)

Sex, Male

28/51 (54.9%)

30/40 (75.0%)

14/29 (48.3%)

77/115 (67.0%)

Body mass index


43/51 (84.3%)

31/39 (79.5%)

29/29 (100.0%)

88/111 (79.3%)


8/51 (15.7%)

1/39 (2.6%)

0/29 (0.0%)

16/111 (14.4%)


0/51 (0.0%)

7/39 (17.9%)

0/29 (0.0%)

7/111 (6.3%)

Family history of allergic diseases, yes

28/50 (56.0%)

14/38 (36.8%)

16/26 (61.5%)

39/113 (34.5%)

Monthly income, < 3000 USD

32/49 (65.3%)

21/40 (52.5%)

26/28 (92.9%)

73/113 (64.6%)

Maternal educational level

 Less than high school

30/46 (65.2%)

15/36 (41.7%)

21/29 (72.4%)

67/111 (60.4%)

Exposure to ETS, yes

12/50 (24.0%)

20/38 (52.6%)

17/28 (60.7%)

57/111 (51.4%)

Atopy, yes

46/51 (90.2%)

18/40 (45.0%)

27/29 (93.1%)

13/114 (11.4%)

Onset age of first wheezing

  < 1 year

10/44 (22.7%)

5/35 (14.3%)

0/28 (0.0%)

18/103 (17.5%)

 1–3 years

9/44 (20.5%)

9/35 (25.7%)

6/28 (21.4%)

46/103 (44.7%)

 3–6 years

11/44 (25.0%)

14/35 (40.0%)

17/28 (60.7%)

29/103 (28.2%)

  ≥ 6 years

14/44 (31.8%)

7/35 (20.0%)

5/28 (17.9%)

10/103 (9.7%)

Number of absence days from school due to wheezinga


51/51 (100.0%)

24/40 (60.0%)

18/29 (62.1%)

115/115 (100.0%)

 1–3 days

0/51 (0.0%)

15/40 (37.5%)

5/29 (17.2%)

0/115 (0.0%)

 4–6 days

0/51 (0.0%)

0/40 (0.0%)

4/29 (13.8%)

0/115 (0.0%)

  ≥ 7 days

0/51 (0.0%)

1/40 (2.5%)

2/29 (6.9%)

0/115 (0.0%)

Number of nocturnal awakening days due to wheezinga


51/51 (100.0%)

13/40 (32.5%)

9/29 (31.0%)

115/115 (100.0%)

  < 1 /week

0/51 (0.0%)

24/40 (60.0%)

15/29 (51.7%)

0/115 (0.0%)

  ≥ 1 /week

0/51 (0.0%)

3/40 (7.5%)

5/29 (17.2%)

0/115 (0.0%)

Number of asthma attacka


42/51 (82.4%)

0/40 (0.0%)

0/29 (0.0%)

108/115 (93.9%)


9/51 (17.6%)

26/40 (65.0%)

11/29 (37.9%)

7/115 (6.1%)


0/51 (0.0%)

8/40 (20.0%)

7/29 (24.1%)

0/115 (0.0%)

  ≥ 13

0/51 (0.0%)

6/40 (15.0%)

11/29 (37.9%)

0/115 (0.0%)

Asthma treatmenta, yes

19/51 (37.3%)

32/40 (80.0%)

27/29 (93.1%)

21/115 (18.3%)

Total serum IgE ≥100 IU/mL

50/51 (98.0%)

24/40 (60.0%)

28/28 (100.0%)

30/111 (27.0%)

Blood eosinophil ≥4%

28/51 (54.9%)

24/40 (60.0%)

23/29 (79.3%)

29/112 (25.9%)

BHR <8 mg/ml

22/49 (44.9%)

12/38 (31.6%)

21/26 (80.8%)

25/108 (23.1%)

FEV1 < 80% predicted

7/47 (14.9%)

0/39 (0.0%)

4/29 (13.8%)

1/113 (0.9%)

FEV1/FVC <80% predicted

1/49 (2.0%)

0/40 (0.0%)

11/29 (37.9%)

8/114 (7.0%)

FEF25–75% <65% predicted

5/29 (17.2%)

0/30 (0.0%)

8/19 (42.1%)

4/90 (4.4%)

ain the previous 12 months

Definition of abbreviations: BHR methacholine PC20 < 8mg/ml, ETS environmental tobacco smoke, FEF 25–75% Forced expiratory flow at 25–75% of forced vital capacity, FEV 1 Forced expiratory volume in 1 s, FVC Forced vital capacity, USD US dollars

Cluster 2 (less atopic, troublesome asthma)

Forty children (17.0% of the enrolled children) were assigned to cluster 2, characterized by normal lung function with moderately troublesome asthma with high SES, reflected by monthly household income and maternal educational levels. Obesity was observed in 17.9% of the children (n = 7/39) in this cluster. This cluster had the lowest prevalence of low monthly income (n = 21/40, 52.5%) and lower levels of maternal educational levels (n = 15/36, 41.7%). All the children in this cluster experienced wheezing more than once in the previous 12 months at the time of enrollment. Sixteen children (40%) of this cluster experienced absence days from school due to asthma attacks in the previous 12 months, and 27 children (67.5%) experienced nocturnal awakening days due to asthma attacks in the previous 12 months at the time of enrollment. However, all the children had normal lung function. Out of 38, 12 children (31.6%) in this cluster showed positive BHR.

Cluster 3 (atopic, persistently troublesome asthma)

Twenty–nine children (12.3%) were grouped in this cluster. Children in this cluster experienced their first wheezing after 1 year of age. All children had normal BMI. Further, the highest prevalence of low SES was observed. The quality of life due to asthma attacks was poorest in this cluster; all the children experienced wheezing in the previous 12 months, and 27 children (93.1%) were treated for asthma attacks in the previous 12 months at the time of enrollment; 11 children (37.9%) were reported to have absence days from school due to asthma attacks and 20 children (69.0%) experienced nocturnal awakening owing to asthma attacks in the previous 12 months at the time of enrollment. The total serum IgE levels were >100 IU/mL for all the children, and blood eosinophil percentage was >4% in 23 children (79.3%). Out of 26 children, 21 (80.8%) showed positive BHR, and four of the total children showed <80% of FEV1 % predicted.

Cluster 4 (less-atopic mild asthma)

One hundred fifteen children (48.9%) were included in this cluster, which was the largest. The prevalence of atopy and parental history of allergic diseases including AD, AR, and asthma was the lowest (n = 39/113, 34.5%; n = 13/114, 11.4%, respectively). Seven children (6.1%) experienced asthma attacks, although 21 children (18.3%) were treated due to asthma exacerbation during the previous 12 months at the time of enrollment. However, none of the children in this cluster were reported to have any absence days from school or nocturnal awakening days due to asthma attacks in the previous 12 months at the time of enrollment. The prevalence of ≥100 IU/ml of total serum IgE levels (n = 30/111; 27.0%) and ≥4% of blood eosinophil percentages (n = 29/112; 25.9%) was the lowest in this cluster. Additionally, the prevalence of BHR was the lowest (n = 25/108; 23.1%), and most of the children in this cluster showed normal lung function.


The prevalence of parent-reported, physician-diagnosed AR in lifetime was the highest in cluster 4 (60.0%), followed by clusters 2 (57.5%) and 1 (56.0%) (Table 3). The prevalence of parent-reported, physician-diagnosed AD in lifetime was the highest in cluster 3 (55.2%), followed by clusters 2 (50.0%) and 4 (39.5%). The prevalence of parent-reported, physician-diagnosed food allergy in lifetime was the highest in cluster 2 (15.4%), followed by cluster 4 (12.4%), although no significant differences were observed among the clusters.
Table 3

Comorbidities of allergic diseases at the time of enrollment according to asthma clusters


Cluster 1

Cluster 2

Cluster 3

Cluster 4

P value

Allergic rhinitis

28/50 (56.0)

23/40 (57.5)

00/29 (34.5)

46/115 (60.0)


Atopic dermatitis

18/50 (36.0)

20/40 (50.0)

16/29 (55.2)

45/114 (39.5)


Food allergy

2/49 (4.1)

6/39 (15.4)

2/28 (7.1)

14/113 (12.4)



Mean values of the total serum IgE and blood eosinophil percentages showed a decreasing pattern as the children grew in each cluster (Table 4). The total serum IgE levels and blood eosinophil percentages were the highest in cluster 3, followed by clusters 2, 4, and 1. This decreasing pattern was observed even at the final survey.
Table 4

Total serum IgE levels and blood eosinophil percentages at the time of enrollment and follow-up


At the time of enrollment


Cluster 1

Cluster 2

Cluster 3

Cluster 4

P value

Cluster 1

Cluster 2

Cluster 3

Cluster 4

P value

Total serum IgE levels (IU/mL)

mean (N)

746.8 (51)

306.1 (40)

1238.0 (28)

175.0 (111)


378.9 (35)

148.7 (22)

447.7 (17)

77.2 (68)


95% CI











Blood eosinophil, %

mean (N)

5.7 (51)

5.3 (40)

6.2 (29)

3.6 (112)


4.4 (35)

3.8 (22)

5.4 (17)

3.1 (69)


95% CI











Definition of abbreviations: CI confidence intervals, N number

The prevalence of sensitization to indoor allergens was the highest in cluster 3 (n = 25/29; 86.2%), followed by clusters 1 (n = 42/51; 82.4%), 2 (n = 17/40; 42.5%), and 4 (n = 11/114; 9.7%) at the time of enrollment (Fig. 1). The statistical differences in the prevalence of sensitization were observed between clusters 2 and 3 and between clusters 1 and 4. The same pattern was observed at the final survey with no statistical differences between clusters 2 and 3. Cluster 3 showed the highest prevalence of sensitization to outdoor allergens (n = 8/29; 27.6%) at the time of enrollment without statistical differences.
Fig. 1

Sensitization according to asthma phenotypes. a Sensitizations to indoor allergens both at the time of enrollment and follow-up. b Sensitizations to outdoor allergens both at the time of enrollment and follow-up. $ P value < 0.005, compared with cluster 1. # P value < 0.005, compared with cluster 2. & P value < 0.005, compared with cluster 3

Lung function

At the time of enrollment, mean values of FEV1 % predicted in cluster 3 (mean, 94.4; 95% confidence intervals [CI], 87.9–101.0) were significantly lower than that in cluster 1 (mean, 98.3; 95% CI, 94.0–102.5), compared with clusters 2 (mean, 107.0; 95% CI, 103.0–111.1) and 4 (mean, 104.6; 95% CI, 102.4–106.8) (Fig. 2). This pattern in FEV1 % predicted was also observed at the final survey. Mean levels of FEV1/FVC % predicted in cluster 3 (mean, 83.3; 95% CI, 80.7–85.9) were significantly lower compared with those in cluster 2 (mean, 89.9; 95% CI, 88.4–91.5) at the time of enrollment. However, no significant differences in FEV1/FVC % predicted between clusters 2 (mean, 88.9; 95% CI, 85.9–92.0) and 3 (mean, 84.4; 95% CI, 77.4–91.3) were observed, although FEV1/FVC % predicted was the lowest in cluster 3 at the final survey.
Fig. 2

Changes in pulmonary function according to asthma clusters. a FEV1 % predicted at the time of enrollment. b FEV1 % predicted at the time of follow-up. c FEV1/FVC % predicted at the time of enrollment. d FEV1/FVC % predicted at the time of follow-ups. e FEF25–75% % predicted at the time of enrollment. f FEF25–75% % predicted at the time of follow-ups. * P value < 0.05

The levels of FEF25–75% % predicted were lowest in cluster 3. The levels of FEF25–75% % predicted were significantly lower in cluster 3 compared with those in cluster 2 both at the time of enrollment and at the final survey. However, no significant differences in the levels of FEF25–75% % predicted were observed between clusters 1 and 4 both at the time of enrollment and at the final survey.


The prevalence of new-onset BHR during the prospective follow-ups was the highest in cluster 3 (n = 3/5; 60.0%), followed by clusters 1 (n = 8/21; 38.1%) and 4 (n = 5/64; 7.8%) (Table 5). The persistence rate of BHR was the highest in cluster 1 (n = 12/20; 60.0%), followed by cluster 3 (n = 7/15; 46.7%). The remission rate of BHR was the highest in cluster 4 (n = 20/22; 90.9%), followed by cluster 2 (n = 6/7; 85.7%). The prevalence of asthma attacks during the previous 12 months at follow-ups was the highest in cluster 3 (n = 13/24; 54.2%), followed by cluster 2 (n = 10/26; 38.5%).
Table 5

Prognosis of asthma phenotypes according to asthma clusters


Cluster 1

Cluster 2

Cluster 3

Cluster 4

P value

New-onset BHR, n (%)

8/21 (38.1)

1/17 (5.9)

3/5 (60.0)

5/64 (7.8)*


Persistence of BHR, n (%)

12/20 (60.0)

1/7 (14.3)

7/15 (46.7)

2/22 (9.1)*


Remission of BHR, n (%)

8/20 (40.0)

6/7 (85.7)

8/15 (53.3)

20/22 (90.9)*


Wheezing in the previous 12 months at follow-up, n (%)

4/23 (17.4)

10/26 (38.5)

13/24 (54.2)

3/25 (12.0)†


New-onset AR, n (%)

9/18 (50.0)

6/9 (66.7)

7/15 (46.7)

15/60 (25.0)


New sensitization on SPT

13/42 (31.0)

4/26 (15.4)

6/23 (26.1)

27/91 (29.7)


Definition of abbreviations: AR allergic rhinitis, BHR methacholine PC20 < 8 mg/ml, SPT skin prick test

* P value < 0.005, compared with cluster 1

P value < 0.005, compared with cluster 3


In the present prospective general-population based study, we identified four distinct asthma clusters in school-age children (Additional file 1: Figure S1). More frequent asthma attacks in the previous 12 months at the time of enrollment were observed in two asthma phenotypes (clusters 2 and 3) among the four clusters. However, the two clusters differed in SES; cluster 3 included children with lower SES, whereas prevalence of high SES was higher in cluster 2. Additionally, cluster 3 had a higher prevalence of atopy combined with an increased prevalence of BHR and decreased lung function. The other two clusters (clusters 1 and 4) included less frequent asthma attacks; however, they differed in atopy prevalence. The higher prevalence of atopy and ≥100 IU/mL of total serum IgE was observed in cluster 1. However, 17.6% of the children in cluster 1 had experienced asthma attacks in the previous 12 months at the time of enrollment; positive BHR was observed in about 45%. These findings suggest that asthma phenotypes in children can be determined based on atopy and SES combined with symptom severity. The comprehensive consideration of diverse factors might contribute to the identification of the origin of asthma. To our knowledge, this study is the first to identify asthma phenotypes using severity, atopy burden, and SES in school-age children.

Guidelines for asthma treatment were suggested based on asthma severity [8]. However, this classification lacks in the reflection of heterogeneity in asthma phenotypes, especially in children with unpredictable prognosis as they grew. To identify the phenotypic heterogeneity over time in children with asthma, strenuous efforts have been made via the epidemiologic studies [6, 9, 23, 28, 29]. The previous studies classified wheezing phenotype according to onset age of wheezing and/or persistence of asthma symptoms, especially in preschool children [29, 30]. Another study has identified wheezing phenotypes based on symptom patterns such as episodic viral or multiple-trigger wheeze in preschool children [31]. Although these phenotypes have contributed to the prediction of prognosis in preschool children with wheezing, the previous studies could not identify the underlying complex pathologies in each phenotype; because associated factors were partially considered for the heterogeneity of asthma.

We included BMI, parental history of allergic diseases, lung function, atopy, asthma symptoms, and SES for LCA to better understand the complex interactions between environmental and genetic factors according to asthma phenotypes and to identify the underlying mechanism. SES, one of the differentiating points in two troublesome asthma phenotypes in the present study, has been identified to have an inverse association with asthma regardless of age [32]. In another study, low-income levels were associated with the development and persistency of asthma [3335]. Although the mechanisms underlying the association between SES and asthma have not been clearly identified, this factor might affect the pathology of asthma via complex interactions with other factors in some phenotypes of asthma, as identified in this study. This might be partially explained by inconsistent associations determined between SES and asthma [36]. The association between SES and asthma can be partially explained by exposure degree to allergens, gene by environmental interactions, nutrition, and psychosocial factors [3739]. These factors are considered to affect the immune development through balance between Th1 and Th2 immune response [40, 41], which might be linked to alteration in microbiota in humans [4244]. Alterations of colonization in gut or airway might affect the development and regulation of immune system with shaping of airway structures and gene expression through gene-environment interactions; therefore, they might contribute to the development and severity of asthma [4548]. Also, low SES might affect the control of asthma in the aspect of adherence or compliance to asthma treatment. However, higher levels of total serum IgE as well as eosinophil (%) both at the time of enrollment and follow-up might imply the different pathophysiologies underlying more severe asthma, as shown in cluster 3.

The previous studies have identified that atopy, late-onset wheeze, and/or BHR are associated with persistent or severe asthma in children [6, 11, 30]. However, in another study, multiple trigger wheezing is associated with decreased pulmonary function regardless atopy or current wheezing [5]. These contrasting findings suggest that atopy alone might not explain the severity and heterogeneity in asthma; however, it can differently affect the prognosis of asthma depending on phenotypes. In the present study, cluster 1 was characterized by high prevalence of atopy similar to cluster 3, in which atopy was highly prevalent with troublesome and persistent asthma; however, children in cluster 1 had mild symptoms of asthma. Based on the previous and present studies [5, 14], certain phenotypes of asthma are influenced by complex combinations that affect interactions between innate and adaptive immunity [49, 50]; therefore, this might contribute to the understanding of the mechanisms and the development of individualized treatment according to asthma phenotypes.

The total serum IgE levels and blood eosinophil percentages had similar patterns with degrees of atopy and BHR, rather than asthma severity, similar to the previous study [51]. The assessments of these biomarkers help to identify asthma phenotypes with partial explanation for the underlying pathology in these phenotypes. However, high atopy is not equal to persistence or severity of asthma, as reflected in cluster 1 in the present and previous studies [4, 5]; it would be more helpful to consider the more complex factors in the prediction of prognosis in each phenotype of asthma.

In contrast to the previous studies on wheezing phenotypes in preschool children [3, 6, 11, 14], a wide overlap in onset-age of the first wheezing was observed across four asthma phenotypes in the present study. Overlapping features are commonly encountered in the results of LCA and this might practically reflect the real phenomenon in the classification of asthma. In the previous studies, late-onset asthma was classified as an independent phenotype regardless of persistency of asthma [6, 14]. This might suggest that onset-age of the first wheezing does not contribute to the differences in underlying mechanisms according to asthma phenotypes.

Discrepancy in frequency of asthma symptoms and treatment in the previous 12 months at the time of enrollment across asthma phenotypes might be attributable to the early application of asthma treatment in the presence of asthma like symptoms such as cough. The prevalence of this discrepancy was highest in cluster 4, which was characterized by less atopic mild asthma phenotype. This classification might help to predict the treatment behavior of asthma in each phenotype, and thereby, facilitate the targeted therapy and patient education based on asthma phenotypes. This might be useful to reduce the socioeconomic burden owing to patient morbidity and mortality in asthma [52].

There are some limitations in this study. Retrospective recall of childhood events has limitations in accurate assessment; however, this study included a prospective follow-up with 2-year intervals for 4 years in school-age children selected from a nationwide general population. Physician-diagnosed asthma was also based on parent reports. However, the prevalence of asthma diagnosis in the present study was 9.4% (n = 235/2491), which is compatible with the worldwide prevalence of asthma (6–9%) [53]. Follow-ups were relatively short. The assessment on treatment pattern and responses to treatment was lacking. However, disease burden was assessed using questions on the frequency of nocturnal awakening days and absence days from school owing to asthma attacks; this information can indirectly reflect the diseases course. Moreover, we enrolled nationwide school-age children with asthma; therefore, the results of this study can be generalized. Although this study included several environmental factors for cluster analysis, further asthma phenotypes including index representatives of air pollution are needed.


We identified four asthma phenotypes in school-age children based on severity of asthma, SES, and atopy. Persistently troublesome asthma was divided into two phenotypes based on SES with atopy. Mild, intermittent asthma was divided into two phenotypes according to the degree of atopy. These findings suggest the importance of factors that affect SES in the classification of asthma phenotypes to identify the mechanisms underlying the heterogeneous asthma phenotypes. Consideration of diverse factors including environmental factors might be more helpful to identify the underlying pathogenesis according to asthma phenotype, prevent progression of asthma, and elucidate targeted therapy.



Atopic dermatitis


Allergic rhinitis


Bronchial hyperresponsiveness


Children’s HEalth and Environmental Research


Environmental tobacco smoke


Food allergy


Forced expiratory flow at 25–75% of forced vital capacity


Forced expiratory volume in 1 second


Forced vital capacity


Cumulative methacholine dose that caused a 20% reduction in FEV1


Socioeconomic status


Skin prick test.



This study was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), which is funded by the Ministry of Science, ICT and Future Planning (NRF-2014R1A2A1A10050687), and partially supported by the Korean Ministry of the Environment, Republic of Korea.

Availability of data and materials

The data supporting the findings of the manuscript are available on request to the first author (Eun Lee).

Authors’ contributions

EL, SHL, JWK and SJH contributed to the analysis and interpretation of data and drafting of the manuscript. EL, JWK, YHK, JY, HJC, SIY, YHJ, JHS, HBK, SYL, HJK, and SJH contributed to the collection of data and the thorough review of the manuscript. SJH contributed to the conception of this study, its design, and interpretation of the data as co-corresponding authors. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable (The manuscript does not contain any individual persons’ data).

Ethics approval and consent to participate

The study protocol was approved by the institutional review board of the University of Ulsan College of Medicine (IRB no. 2006–0081). Written informed consent was obtained from the parents or guardians of all participants.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Department of Pediatrics, Chonnam National University Hospital
Asan Institute for Life Sciences, University of Ulsan College of Medicine
Department of Pediatrics, Seoul National University Bundang Hospital
Department of Pediatrics, Gyeongsang National University Changwon Hospital
Department of Pediatrics, Childhood Asthma and Atopy Center, Environmental Health Center, Asan Medical Center, University of Ulsan College of Medicine
Department of Pediatrics, Hallym University Sacred Heart Hospital
Department of Pediatrics, CHA University School of Medicine
Department of Pediatrics, Pusan National University Yangsan Hospital
Department of Pediatrics, Dankook University College of Medicine
Department of Pediatrics, Sanggye Paik Hospital, Inje University College of Medicine
Department of Preventive Medicine, Dankook University College of Medicine


  1. Brusasco V, Crimi E, Pellegrino R. Airway hyperresponsiveness in asthma: not just a matter of airway inflammation. Thorax. 1998;53(11):992–8.View ArticlePubMedPubMed CentralGoogle Scholar
  2. Stevens CA, Turner D, Kuehni CE, Couriel JM, Silverman M. The economic impact of preschool asthma and wheeze. Eur Respir J. 2003;21(6):1000–6.View ArticlePubMedGoogle Scholar
  3. Fitzpatrick AM, Teague WG, Meyers DA, Peters SP, Li X, Li H, Wenzel SE, Aujla S, Castro M, Bacharier LB, et al. Heterogeneity of severe asthma in childhood: confirmation by cluster analysis of children in the National Institutes of Health/National Heart, Lung, and Blood Institute Severe Asthma Research Program. J Allergy Clin Immunol. 2011;127(2):382–9. e381-313.View ArticlePubMedGoogle Scholar
  4. Haldar P, Pavord ID, Shaw DE, Berry MA, Thomas M, Brightling CE, Wardlaw AJ, Green RH. Cluster analysis and clinical asthma phenotypes. Am J Respir Crit Care Med. 2008;178(3):218–24.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Sonnappa S, Bastardo CM, Wade A, Saglani S, McKenzie SA, Bush A, Aurora P. Symptom-pattern phenotype and pulmonary function in preschool wheezers. J Allergy Clin Immunol. 2010;126(3):519–26. e511-517.View ArticlePubMedGoogle Scholar
  6. Savenije OE, Granell R, Caudri D, Koppelman GH, Smit HA, Wijga A, de Jongste JC, Brunekreef B, Sterne JA, Postma DS, et al. Comparison of childhood wheezing phenotypes in 2 birth cohorts: ALSPAC and PIAMA. J Allergy Clin Immunol. 2011;127(6):1505–12. e1514.View ArticlePubMedGoogle Scholar
  7. Pite H, Gaspar A, Morais-Almeida M. Preschool-age wheezing phenotypes and asthma persistence in adolescents. Allergy Asthma Proc. 2016;37(3):231–41.View ArticlePubMedGoogle Scholar
  8. Boulet LP, FitzGerald JM, Reddel HK. The revised 2014 GINA strategy report: opportunities for change. Curr Opin Pulm Med. 2015;21(1):1–7.View ArticlePubMedGoogle Scholar
  9. Rusconi F, Galassi C, Corbo GM, Forastiere F, Biggeri A, Ciccone G, Renzoni E. Risk factors for early, persistent, and late-onset wheezing in young children. SIDRIA Collaborative Group. Am J Respir Crit Care Med. 1999;160(5 Pt1):1617–22.View ArticlePubMedGoogle Scholar
  10. Duffy DL, Mitchell CA, Martin NG. Genetic and environmental risk factors for asthma: a cotwin-control study. Am J Respir Crit Care Med. 1998;157(3 Pt1):840–5.View ArticlePubMedGoogle Scholar
  11. Sherriff A, Peters TJ, Henderson J, Strachan D. Risk factor associations with wheezing patterns in children followed longitudinally from birth to 3(1/2) years. Int J Epidemiol. 2001;30(6):1473–84.View ArticlePubMedGoogle Scholar
  12. Li YF, Langholz B, Salam MT, Gilliland FD. Maternal and grandmaternal smoking patterns are associated with early childhood asthma. Chest. 2005;127(4):1232–41.PubMedGoogle Scholar
  13. Breton CV, Byun HM, Wenten M, Pan F, Yang A, Gilliland FD. Prenatal tobacco smoke exposure affects global and gene-specific DNA methylation. Am J Respir Crit Care Med. 2009;180(5):462–7.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Belgrave DC, Simpson A, Semic-Jusufagic A, Murray CS, Buchan I, Pickles A, Custovic A. Joint modeling of parentally reported and physician-confirmed wheeze identifies children with persistent troublesome wheezing. J Allergy Clin Immunol. 2013;132(3):575–83. e512.View ArticlePubMedGoogle Scholar
  15. Tse SM, Rifas-Shiman SL, Coull BA, Litonjua AA, Oken E, Gold DR. Sex-specific risk factors for childhood wheeze and longitudinal phenotypes of wheeze. J Allergy Clin Immunol. 2016;138(6):1561-1568:e6.Google Scholar
  16. Lee E, Lee SH, Kwon JW, Kim YH, Cho HJ, Yang SI, Jung YH, Kim HY, Seo JH, Kim BJ, et al. Atopic dermatitis phenotype with early onset and high serum IL-13 is linked to the new development of bronchial hyperresponsiveness in school children. Allergy. 2016;71(5):692–700.View ArticlePubMedGoogle Scholar
  17. Kim BJ, Seo JH, Jung YH, Kim HY, Kwon JW, Kim HB, Lee SY, Park KS, Yu J, Kim HC, et al. Air pollution interacts with past episodes of bronchiolitis in the development of asthma. Allergy. 2013;68(4):517–23.View ArticlePubMedGoogle Scholar
  18. Kim BJ, Lee SY, Kwon JW, Jung YH, Lee E, Yang SI, Kim HY, Seo JH, Kim HB, Kim HC, et al. Traffic-related air pollution is associated with airway hyperresponsiveness. J Allergy Clin Immunol. 2014;133(6):1763–5. e1762.View ArticlePubMedGoogle Scholar
  19. Fukutomi Y, Taniguchi M. Sensitization to fungal allergens: Resolved and unresolved issues. Allergol Int. 2015;64(4):321–31.View ArticlePubMedGoogle Scholar
  20. Eder W, Ege MJ, von Mutius E. The asthma epidemic. N Engl J Med. 2006;355(21):2226–35.View ArticlePubMedGoogle Scholar
  21. Asher MI, Keil U, Anderson HR, Beasley R, Crane J, Martinez F, Mitchell EA, Pearce N, Sibbald B, Stewart AW, et al. International Study of Asthma and Allergies in Childhood (ISAAC): rationale and methods. Eur Respir J. 1995;8(3):483–91.View ArticlePubMedGoogle Scholar
  22. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, van der Grinten CP, Gustafsson P, et al. Standardisation of spirometry. Eur Respir J. 2005;26(2):319–38.View ArticlePubMedGoogle Scholar
  23. Kim BJ, Lee SY, Kim HB, Lee E, Hong SJ. Environmental changes, microbiota, and allergic diseases. Allergy, Asthma Immunol Res. 2014;6(5):389–400.View ArticleGoogle Scholar
  24. Kim BJ, Kwon JW, Seo JH, Kim HB, Lee SY, Park KS, Yu J, Kim HC, Leem JH, Sakong J, et al. Association of ozone exposure with asthma, allergic rhinitis, and allergic sensitization. Ann Allergy Asthma Immunol. 2011;107(3):214–9. e211.View ArticlePubMedGoogle Scholar
  25. Siroux V, Basagana X, Boudier A, Pin I, Garcia-Aymerich J, Vesin A, Slama R, Jarvis D, Anto JM, Kauffmann F, et al. Identifying adult asthma phenotypes using a clustering approach. Eur Respir J. 2011;38(2):310–7.View ArticlePubMedGoogle Scholar
  26. Lee E, Lee SH, Kwon JW, Kim Y, Cho HJ, Yang SI, Jung YH, Kim HY, Seo JH, Kim BJ, et al. A rhinitis phenotype associated with increased development of bronchial hyperresponsiveness and asthma in children. Ann Allergy Asthma Immunol. 2016;117(1):21–8. e21.View ArticlePubMedGoogle Scholar
  27. Hagenaars JA, McCutcheon AL. Applied latent class analysis. Cambridge: Cambridge University Press; 2002.View ArticleGoogle Scholar
  28. Cowan K, Guilbert TW. Pediatric asthma phenotypes. Curr Opin Pediatr. 2012;24(3):344–51.View ArticlePubMedGoogle Scholar
  29. Henderson J, Granell R, Heron J, Sherriff A, Simpson A, Woodcock A, Strachan DP, Shaheen SO, Sterne JA. Associations of wheezing phenotypes in the first 6 years of life with atopy, lung function and airway responsiveness in mid-childhood. Thorax. 2008;63(11):974–80.View ArticlePubMedPubMed CentralGoogle Scholar
  30. Martinez FD, Wright AL, Taussig LM, Holberg CJ, Halonen M, Morgan WJ. Asthma and wheezing in the first six years of life. N Engl J Med. 1995;332(3):133–8.View ArticlePubMedGoogle Scholar
  31. Brand PL, Baraldi E, Bisgaard H, Boner AL, Castro-Rodriguez JA, Custovic A, de Blic J, de Jongste JC, Eber E, Everard ML, et al. Definition, assessment and treatment of wheezing disorders in preschool children: an evidence-based approach. Eur Respir J. 2008;32(4):1096–110.View ArticlePubMedGoogle Scholar
  32. Gong T, Lundholm C, Rejno G, Mood C, Langstrom N, Almqvist C. Parental socioeconomic status, childhood asthma and medication use--a population-based study. PLoS One. 2014;9(9):e106579.View ArticlePubMedPubMed CentralGoogle Scholar
  33. Kozyrskyj AL, Kendall GE, Jacoby P, Sly PD, Zubrick SR. Association between socioeconomic status and the development of asthma: analyses of income trajectories. Am J Public Health. 2010;100(3):540–6.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Lindbaek M, Wefring KW, Grangard EH, Ovsthus K. Socioeconomical conditions as risk factors for bronchial asthma in children aged 4–5 years. Eur Respir J. 2003;21(1):105–8.View ArticlePubMedGoogle Scholar
  35. Mielck A, Reitmeir P, Wjst M. Severity of childhood asthma by socioeconomic status. Int J Epidemiol. 1996;25(2):388–93.View ArticlePubMedGoogle Scholar
  36. Rona RJ. Asthma and poverty. Thorax. 2000;55(3):239–44.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Lanphear BP, Kahn RS, Berger O, Auinger P, Bortnick SM, Nahhas RW. Contribution of residential exposures to asthma in us children and adolescents. Pediatrics. 2001;107(6):E98.View ArticlePubMedGoogle Scholar
  38. Chen E, Hanson MD, Paterson LQ, Griffin MJ, Walker HA, Miller GE. Socioeconomic status and inflammatory processes in childhood asthma: the role of psychological stress. J Allergy Clin Immunol. 2006;117(5):1014–20.View ArticlePubMedGoogle Scholar
  39. Trompette A, Gollwitzer ES, Yadava K, Sichelstiel AK, Sprenger N, Ngom-Bru C, Blanchard C, Junt T, Nicod LP, Harris NL, et al. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat Med. 2014;20(2):159–66.View ArticlePubMedGoogle Scholar
  40. Gern JE, Visness CM, Gergen PJ, Wood RA, Bloomberg GR, O’Connor GT, Kattan M, Sampson HA, Witter FR, Sandel MT, et al. The Urban Environment and Childhood Asthma (URECA) birth cohort study: design, methods, and study population. BMC Pulm Med. 2009;9:17.View ArticlePubMedPubMed CentralGoogle Scholar
  41. Celedon JC, Milton DK, Ramsey CD, Litonjua AA, Ryan L, Platts-Mills TA, Gold DR. Exposure to dust mite allergen and endotoxin in early life and asthma and atopy in childhood. J Allergy Clin Immunol. 2007;120(1):144–9.View ArticlePubMedPubMed CentralGoogle Scholar
  42. Rook GA, Raison CL, Lowry CA. Microbial ‘old friends’, immunoregulation and socioeconomic status. Clin Exp Immunol. 2014;177(1):1–12.View ArticlePubMedPubMed CentralGoogle Scholar
  43. Hanski I, von Hertzen L, Fyhrquist N, Koskinen K, Torppa K, Laatikainen T, Karisola P, Auvinen P, Paulin L, Makela MJ, et al. Environmental biodiversity, human microbiota, and allergy are interrelated. Proc Natl Acad Sci U S A. 2012;109(21):8334–9.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Fujimura KE, Demoor T, Rauch M, Faruqi AA, Jang S, Johnson CC, Boushey HA, Zoratti E, Ownby D, Lukacs NW, et al. House dust exposure mediates gut microbiome Lactobacillus enrichment and airway immune defense against allergens and virus infection. Proc Natl Acad Sci U S A. 2014;111(2):805–10.View ArticlePubMedGoogle Scholar
  45. Huang YJ, Nelson CE, Brodie EL, Desantis TZ, Baek MS, Liu J, Woyke T, Allgaier M, Bristow J, Wiener-Kronish JP, et al. Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma. J Allergy Clin Immunol. 2011;127(2):372–81. e371-373.View ArticlePubMedGoogle Scholar
  46. Sekhon HS, Keller JA, Proskocil BJ, Martin EL, Spindel ER. Maternal nicotine exposure upregulates collagen gene expression in fetal monkey lung. Association with alpha7 nicotinic acetylcholine receptors. Am J Respir Cell Mol Biol. 2002;26(1):31–41.View ArticlePubMedGoogle Scholar
  47. Huang YJ, Nariya S, Harris JM, Lynch SV, Choy DF, Arron JR, Boushey H. The airway microbiome in patients with severe asthma: associations with disease features and severity. J Allergy Clin Immunol. 2015;136(4):874–84.View ArticlePubMedPubMed CentralGoogle Scholar
  48. Holt PG. The mechanism or mechanisms driving atopic asthma initiation: The infant respiratory microbiome moves to center stage. J Allergy Clin Immunol. 2015;136(1):15–22.View ArticlePubMedGoogle Scholar
  49. Holt PG, Strickland DH. Interactions between innate and adaptive immunity in asthma pathogenesis: new perspectives from studies on acute exacerbations. J Allergy Clin Immunol. 2010;125(5):963–72. quiz 973–964.View ArticlePubMedGoogle Scholar
  50. Holt PG, Sly PD. Interaction between adaptive and innate immune pathways in the pathogenesis of atopic asthma: operation of a lung/bone marrow axis. Chest. 2011;139(5):1165–71.View ArticlePubMedGoogle Scholar
  51. Kim BS, Jin HS, Kim HB, Lee SY, Kim JH, Kwon JW, Kim BJ, Yu J, Yoo S, Hong SJ. Airway hyperresponsiveness is associated with total serum immunoglobulin E and sensitization to aeroallergens in Korean adolescents. Pediatr Pulmonol. 2010;45(12):1220–7.View ArticlePubMedGoogle Scholar
  52. Long AA. The burden of asthma and improving patient outcomes. Am J Manag Care. 2011;17 Suppl 3:S75–81.PubMedGoogle Scholar
  53. Asher MI, Montefort S, Bjorksten B, Lai CK, Strachan DP, Weiland SK, Williams H. Worldwide time trends in the prevalence of symptoms of asthma, allergic rhinoconjunctivitis, and eczema in childhood: ISAAC phases one and three repeat multicountry cross-sectional surveys. Lancet. 2006;368(9537):733–43.View ArticlePubMedGoogle Scholar


© The Author(s). 2017