The COPD multi-dimensional phenotype: A new classification from the STORICO Italian observational study
Raffaele Antonelli Incalzi aff001; Giorgio Walter Canonica aff002; Nicola Scichilone aff003; Sara Rizzoli aff004; Lucia Simoni aff004; Francesco Blasi aff005;
Authors place of work:
University Biomedical Campus of Rome, Rome, Italy
aff001; Personalized Medicine Asthma & Allergy Clinic Humanitas University Humanitas research Hospital, Rozzano, Milan, Italy
aff002; DIBIMIS, University of Palermo, Piazza delle Cliniche, Palermo, Italy
aff003; Medineos Observational Research, Modena, ltaly
aff004; Department of Pathophysiology and Transplantation, University of Milan, Internal Medicine Department, Respiratory Unit and Cystic Fibrosis Adult Center Fondazione IRCCS Cà Granda Maggiore Hospital, Milan, Italy
Published in the journal:
PLoS ONE 14(9)
This paper is aimed to (i) develop an innovative classification of COPD, multi-dimensional phenotype, based on a multidimensional assessment; (ii) describe the identified multi-dimensional phenotypes.
An exploratory factor analysis to identify the main classificatory variables and, then, a cluster analysis based on these variables were run to classify the COPD-diagnosed 514 patients enrolled in the STORICO (trial registration number: NCT03105999) study into multi-dimensional phenotypes.
The circadian rhythm of symptoms and health-related quality of life, but neither comorbidity nor respiratory function, qualified as primary classificatory variables. Five multidimensional phenotypes were identified: the MILD COPD characterized by no night-time symptoms and the best health status in terms of quality of life, quality of sleep, level of depression and anxiety, the MILD EMPHYSEMATOUS with prevalent dyspnea in the early-morning and day-time, the SEVERE BRONCHITIC with nocturnal and diurnal cough and phlegm, the SEVERE EMPHYSEMATOUS with nocturnal and diurnal dyspnea and the SEVERE MIXED COPD distinguished by higher frequency of symptoms during 24h and worst quality of life, of sleep and highest levels of depression and anxiety.
Our results showed that properly collected respiratory symptoms play a primary classificatory role of COPD patients. The longitudinal observation will disclose the discriminative and prognostic potential of the proposed multidimensional phenotype.
Trial registration number: NCT03105999, date of registration: 10th April 2017.
Medicine and health sciences – Pulmonology – Chronic obstructive pulmonary disease – Dyspnea – Diagnostic medicine – Signs and symptoms – Pathology and laboratory medicine – Public and occupational health – Physical activity – Biology and life sciences – Physiology – Physiological processes – Coughing – Chronobiology – Circadian rhythms – Anatomy – Body fluids – Mucus – Research and analysis methods – Mathematical and statistical techniques – Statistical methods – Factor analysis – Simulation and modeling – Clustering algorithms – Physical sciences – Mathematics – Statistics – Applied mathematics – Algorithms
Chronic Obstructive Pulmonary Disease (COPD) is an umbrella definition encompassing a variety of clinical and pathophysiological conditions. Selected clinical phenotypes, namely the bronchitic, the emphysematous and the asthma-like, have been repeatedly recognized, but disease pattern not always univocally conforms to one of them . Indeed, the clinical presentation frequently resembles a mixture of the classical phenotypes. Furthermore, comorbidity and age itself largely shape the clinical and health status of COPD patients. Indeed, a senescence-associated secretory phenotype has been identified and might affect the clinical expression of COPD in the elderly . Unfortunately, large pharmacological randomized clinical trials (RCTs) could not shed light on phenotypic variability due to the stringent selecting criteria. Further complicating this issue is the lack of attention to the circadian rhythm of respiratory symptoms in RCTs. On the other hand, several observational studies point out at circadian symptom variability as a correlate of disease severity  or a distinctive clinical trait  or a correlate of derangement of selected health status dimensions . Furthermore, nocturnal symptoms are reported by the majority of COPD patients  and dramatically increased with age in a COPD population over 65 . However, circadian rhythm variability in COPD is only marginally and partially explored in selected RCT with the only intent of providing 24 hours pharmacological coverage and, then, symptoms control. Thus, how circadian rhythm of symptoms affects the clinical patterns and therapeutic needs of COPD has not been the object of interest. Paradoxically, the growing interest in comorbidity of COPD and its role as a determinant of health status seems to overcome the research on the respiratory dimension of COPD.
The STORICO (STudio Osservazionale sulla caratteRizzazione dei sIntomi delle 24 ore nei pazienti con broncopneumopatia cronica ostruttiva, Observational study on characterization of 24-h symptoms in patients with COPD) study assessed the association between circadian rhythm of symptoms and several measures of health status . Thus, it represents the ideal framework to verify whether an in depth characterization of symptoms may allow define clinical variants of COPD with distinctive classificatory and discriminatory properties. Thus, the present paper is aimed to:
develop an innovative classification of COPD based on a multidimensional assessment (named multi-dimensional phenotype) which takes into consideration primarily circadian rhythm of symptoms, but also demographic characteristics, health related quality of life, respiratory function and comorbidities.
describe the patients belonging to the so-identified multi-dimensional phenotypes with respect to body mass index, respiratory parameters, previous exacerbations, level of physical activity, quality of sleep, level of depression and anxiety and ongoing COPD pharmacological therapies.
The final objective of such an attempt is to investigate the discriminatory properties and potential clinical implications of this alternative classificatory method.
Materials and methods
The STORICO study (trial registration number: NCT03105999) enrolled subjects aged ≥50, current/ex-smokers, with a diagnosis of stable COPD according to the GOLD 2014. The study was approved by the ethical committee of the coordinating center (Fondazione Toscana G. Monasterio Pisa, Italy) and was conducted in accordance with the Declaration of Helsinki and the Good Clinical Practices guidelines for observational studies, complying with all requirements of local regulations. Patients provided written, informed consent before study participation. Subjects included, study design and methodology of the STORICO study are extensively described elsewhere .
Patients with available information about early-morning, day- and night-time COPD symptoms and clinical phenotype at enrollment were considered evaluable for analyses at baseline; among these, the ones with available variables for (factor and cluster) analyses (below detailed in Methods paragraph) were analyzed in this paper.
STORICO is an Italian observational cohort multicentre currently ongoing study conducted in 40 pneumology centers. The study lasted from February 2016 (first subject first visit) to June 2018 (last subject last visit); three visits were planned (baseline, 6- and 12-months follow up).
At enrollment patients completed the Night-time, Morning and Day-time Symptoms of COPD questionnaire  (hereafter named “symptoms questionnaire”) covering the frequency and severity of COPD symptoms (breathlessness, coughing, bringing up phlegm or mucus, chest tightness, chest congestion and wheezing) during each part of the day (night-time, early- and day-time). Linguistic validation of the questionnaire in Italian was performed by the authors to ensure accurate translation and a clear understanding of the questionnaire itself .
Health Related Quality of Life was evaluated by the St. George’s Respiratory Questionnaire (SGRQ), in its Symptoms, Activity component and Impact scores [9–10] ranging between 0 (no impairment) and 100 (highest impairment), lower scores corresponding to better health.
Anxiety and depression levels were evaluated through the Hospital Anxiety and Depression Scale (HADS) , a total score (ranging 0–42) and anxiety and depression subscales scores (ranging 0–21) were computed, with higher scores indicating more distress.
The impact on sleep due to respiratory disease was assessed with the total score of COPD and Asthma Sleep Impact Scale (CASIS) , ranging 0–100, with higher scores indicating greater sleep impairment.
Physical activity was assessed by means of the categorical score (low, medium, high physical activity) of the International Physical Activity Questionnaire (IPAQ)  on patients aged 15–69 years. Spirometry was performed according to the recommendations of the American Thoracic Society (ATS) and the European Respiratory Society (ERS) and lung function measurements were done with patients either standing or sitting with the nose clipped after at least 10 minutes rest .
Presence of relevant comorbidities according to clinical judgment, body mass index (BMI), spirometry functional assessment and occurrence of exacerbations in the 5 years before baseline were recorded too as far as ongoing pharmacological therapies for COPD (long-acting beta-agonists (LABA), long-acting muscarinic antagonist (LAMA), inhaled corticosteroid (ICS) and any combination, Others).
In order to identify the multi-dimensional phenotypes (m-phenotypes), a multi-step approach was followed.
STEP 1—Exploratory factor analysis
An exploratory factor analysis (EFA) was performed to find independent latent constructs (factors), not directly measurable and influencing responses on observed variables. EFA is a variable reduction technique which does not impose any preconceived structure on the outcome  and the observed variables included in the model are a linear combination of the underlying factors.
The following items of the symptoms questionnaire were included in the factor analysis: presence/absence of feeling short of breath or breathless (items 18a, 26a, 7a), cough (items 18b, 26b, 7b), bringing up phlegm or mucus (items 18c, 26c, 7c), each symptom evaluated during the 24-hours (i.e. in the early-morning, day-time and night-time). The presence of the symptom was coded as 1 and the absence as 0.
The value of the FEV1% of the predicted, of (symptoms, activity, impact) SGRQ scores, demographic variables (age as continuous variable and gender coded as 1 for males and 2 for females) and presence/absence of relevant comorbidities (cardiac ischemic disease, arterial hypertension, heart failure, atrial fibrillation, diabetes, osteoporosis, depression, kidney insufficiency) were also included in the model (presence was coded as 1 and absence as 0). Orthogonal VARIMAX rotation was applied; factors having an eigenvalue > 1.0 [16–17] and individual variables with higher-than-0.5 loadings on retained factors were retained .
STEP 2 –Definition of classificatory variables
STEP 1 brought to identification of n factors. Then, the variables included in the n-th factor were evaluated and, whenever possible, combined into a new classificatory variable (a classificatory variable for each factor was defined).
STEP 3 –Cluster analysis
A cluster analysis with classificatory variables mentioned at STEP 2 as input was performed. Average linkage was chosen as clustering method and average distance between clusters equal to 0.70 was taken as reference to cut the dendrogram. Quality of obtained clusters was evaluated by means of Semipartial R-squared (i.e. the loss of homogeneity due to combining two clusters to form a new cluster) and R-square (i.e. the proportion of variance accounted for by the clusters).
Once identified, the different m-phenotypes were described with respect to BMI, spirometric parameters, number of previous exacerbations, quality of sleep, level of physical activity, depression and anxiety and therapies for COPD at enrollment by means of median, 25th and 75th percentile for quantitative variables and absolute/relative frequency for categorical ones within each class of m-phenotype. Moreover, although SGRQ scores were included in EFA, SGRQ scores of m-phenotypes were compared to provide a comprehensive description of health status.
Analysis of variance and Kruskal–Wallis test by ranks (on means and medians respectively) and Fisher exact test (for categorical variables) were used to compare variables vs m-phenotype (variable in 5 classes). Upon statistical significance of these tests, Mann-Whitney test on medians and Chi-square or Fisher exact tests were then applied to compare variables among specific pairs of m-phenotypes. Alpha (with Bonferroni correction) was set to 0.0004 considering the total number of performed tests.
Statistical analysis was performed using SAS v9.4 and Enterprise Guide v7.1.
Among the 683 COPD patients enrolled in the STORICO study, 606 (88.7%) subjects (age 71.4±8.2 years, 75.1% males) were deemed evaluable for the analysis at baseline; 92 subjects were then excluded because they had missing information on variables analyzed in the factor and cluster analyses. Violations causing exclusion are shown in Fig 1.
So, 514 (mean±SD age: 71.4± 8.0 years) were deemed evaluable for the analyses here described.
Identification of m-phenotypes
STEP 1—Exploratory factor analysis
As results of the factors analysis, two factors accounting for 82% of the total variability were retained. In Table 1, the factor loadings of individual variables are shown. Factor 1 covers 58.6% of variability and its components are bringing up phlegm or mucus and cough at any time of the day. Factor 2, explaining 23.2% of variability, is composed by breathlessness at any time of the day and by SGRQ scores. A third factor explained a reduced amount of variability (9.5%) and no variables had a factor loading >0.5; so, it was not retained (output of factor analysis is reported in S1 Table and S1 Fig).
STEP 2 –definition of classificatory variables
The items of symptoms questionnaire contributing to each factor (shown in grey in Table 1), were combined into 2 new 3-levels classificatory variables, namely “cough and/or bringing up phlegm or mucus” and “breathlessness”. The levels of them had increasing severity and were defined based on the frequency of occurrence of the symptoms: (i) never during 24 hours, (ii) in the early-morning and/or day-time, but never in night and (iii) in night and in the early-morning and/or day-time. The distribution of patients according to classificatory variables is shown in Table 2.
As SGRQ scores were retained in the second factor, they were described too within the levels of “breathlessness” classificatory variable (Table 2). The (symptoms, activity, impacts) scores significantly increased with level of severity of “breathlessness” (p-value < .0001 for scores in all pairwise comparisons among “breathlessness” categories).
STEP 3 –cluster analysis
The cluster analysis run on previous mentioned classificatory variables brought to five well defined clusters (Semipartial R-squared = 0.0260, R-square = 0.828; dendrogram is shown in S2 Table and S2 Fig). Such clusters are thereafter called m-phenotypes and they are below described.
M-phenotypes and circadian rhythm of symptoms
As shown in Fig 2 the m-phenotypes were characterized with regard to the circadian rhythm of symptoms as follows:
No night-time symptoms (MILD COPD), n = 190. In this m-phenotype, 88 (46.3%) and 64 (33.7%) patients presented cough and/or bringing up phlegm or mucus in the early-morning and day-time, respectively. Breathlessness was reported by 18 (9.5%) and 46 (24.2%) patients in the early-morning and day-time, respectively. Patients belonging to this phenotype have no night symptoms.
Prevalent dyspnea in the early-morning and day-time, not by night (MILD EMPHYSEMATOUS), n = 94. Patients in this cluster suffered mainly from breathlessness in the early-morning (n = 89, 94.7%) and day-time (n = 93, 98.9%). Diurnal cough or bringing up phlegm or mucus was less frequent (63.8% in the early-morning, 59.6% in the day-time), and night symptoms were distinctly rare (only breathlessness in 5.3% patients).
Nocturnal and diurnal cough and phlegm, no nocturnal dyspnea (SEVERE BRONCHITIC), n = 132. These patients suffered from cough and/or bringing up phlegm or mucus, at night (n = 132, 100.0%), in the early-morning (n = 115, 87.1%) and day-time (n = 109, 82.6%). Early-morning dyspnea is present in 47 (35.6%) patients and day-time in 65 (49.2%) patients.
Nocturnal and diurnal dyspnea (SEVERE EMPHYSEMATOUS), n = 18. At variance with MILD EMPHYSEMATOUS, they presented breathlessness also during night (n = 18, 100%).
Higher frequency of symptoms, both cough and breathlessness, during 24h (SEVERE MIXED COPD), n = 80. These patients had the highest frequency of symptoms during 24 hours: 100% of them had cough and breathlessness during night, 77 (96.3%) and 66 (82.5%) had these symptoms in the early-morning and 74 (92.5%) and 73 (91.3%) in the day-time.
M-phenotypes vs BMI, respiratory parameters and number of exacerbations
In Table 3 BMI, spirometric parameters at enrollment and previous exacerbations are described. Patients with MILD COPD had a lower median (25th - 75th percentile) number of exacerbations (1.0 (0.0–3.0)) than patients with SEVERE MIXED COPD (3.0 (1.0–5.0); p-value < .0001). No other statistically significant differences emerged in analyzed parameters.
M-phenotypes vs quality of life, quality of sleep, level of physical activity, anxiety and depression
Quality of life, as expressed by SGRQ scores, was better in MILD COPD than in patients of other m-phenotypes (p-values tests MILD COPD vs each of the other m-phenotypes < .0001) (Table 4). The SEVERE MIXED COPD m-phenotype had higher SGRQ symptoms score than the other phenotypes (p-values < .0001) and higher SGRQ activity and impacts scores than in MILD EMPHYSEMATOUS or SEVERE BRONCHITIC patients (p-values < .0001).
The patients in MILD COPD m-phenotype had better quality of sleep than patients of other m-phenotypes (p-values < .0001) (Table 4). A negative grading of quality of sleep was evident from MILD COPD to MILD EMPHYSEMATOUS and, then, SEVERE BRONCHITIC, SEVERE EMPHYSEMATOUS and SEVERE MIXED COPD.
The patients with MILD COPD had lower HADS total, anxiety and depression scores than MILD EMPHYSEMATOUS or SEVERE MIXED COPD patients (p-values < .0001) and had lower total and depression scores than SEVERE BRONCHITIC patients (p-values < .0001) (Table 4). Differences between HADS scores of patients with SEVERE MIXED COPD vs MILD EMPHYSEMATOUS or SEVERE BRONCHITIC patients were also significant (p-values < .0001).
Low physical activity, as expressed by the IPAQ score, largely prevailed in the patients with SEVERE MIXED COPD and SEVERE EMPHYSEMATOUS m-phenotypes. However, differences between m-phenotypes did not reach significance (Table 4).
M-phenotypes vs pharmacological therapies for COPD
Pharmacological therapy did not differ significantly between m-phenotypes (Table 5).
We found that the circadian rhythm of the three main respiratory symptoms and an index of disease-related health status could classify the vast majority of patients, whereas comorbidity, FEV1, age and gender added little to the classificatory power of the model. Indeed, based on these components of the latent factors, it was possible to identify five multidimensional phenotypes with distinctive pattern of disease (Fig 2): the MILD COPD, the MILD EMPHYSEMATOUS, the SEVERE BRONCHITIC, the SEVERE EMPHYSEMATOUS and the SEVERE MIXED COPD.
Taken together, these findings point out that the “core business” of COPD, the clinical one, can distinguish patients with different disease profiles. The fact that two factors, extracted through a parsimonious model, could explain 82% of intrinsic variability testifies to the strength of the classificatory procedure and to the overall quality of the clustering based mainly on clinical ground.
Consistent with the classical distinction of bronchitic and emphysematous phenotypes, dyspnea and productive cough contributed to define two well distinguished dimensions of the disease. Accordingly, attempts at optimally caring the individual COPD patient should have individually tailored goals: relieving bronchial obstruction might not be the main objective of the care in patients burdened mainly or exclusively with productive cough, whereas relieving dyspnea could be the core therapeutic goal in another group of patients, and the timing of symptoms further adds to the decision making . This seemingly obvious finding has important practical implications. Indeed, FEV1 or the six minute walked distance should not be considered universal effect measures in randomized clinical trials. While these indexes have great merit as overall indexes of disease severity, they might not catch the core expression of COPD and, then, the main expected benefit of therapy in an important proportion of patients. The findings of the current study strongly advocate for the design of multi-outcome or patient-tailored RCTs in the perspective of precision medicine. Pursuing the patient centered objective of care is also expected to increase adherence to therapy. Indeed, the surprisingly high rate of non-adherence to inhaled medication, including inappropriate inhaler use by COPD patients  as well as the usually high rate of enrollees lost in RCT due to withdrawn consent [21–22] might reflect the perception of the care as missing the own needs by many patients.
The circadian rhythm of symptoms played a major classificatory role. Indeed, MILD COPD and MILD EMPHYSEMATOUS patients never/very unfrequently experience nighttime symptoms, which, instead, are always present in SEVERE EMPHYSEMATOUS, SEVERE BRONCHITIC and SEVERE MIXED COPD m-phenotypes. Several evidences support the role of circadian rhythm of symptoms in defining the clinical spectrum of COPD. Indeed, bronchial obstruction is strictly related to nocturnal symptoms  and wheezing has been reported to be the most common of them [23–24]. Interestingly, night time symptoms qualified as the main correlate of depression in the Assess study . However, Physicians significantly underestimate the impact of COPD on the patient’s ability to get up in the morning and on sleep . Unfortunately, only in the last few years attention has been paid to nocturnal symptoms, but mainly to empirically modulate the pharmacological therapy rather than to rigorously understand COPD heterogeneity.
Selected innovative models have recently been proposed to classify COPD patients. They variably rely upon comorbidities and complications of COPD [25–26] or on a comprehensive respiratory function assessment . However, they do not rate symptoms according to a standardized procedure. This likely explains why symptoms played a major classificatory role in ours and not in these models. Thus, two different phylosophies, respectively centered on the non respiratory and the respiratory dimension of COPD, found previous model and the present one. Only the longitudinal observation will allow compare discriminative and prognostic properties of these models. On the other hand, classificatory models including analytes such as cytokines or genetic markers are intended to define idyotypes and not phenotypes and, thus, are not comparable with our prettily clinical and, then, phenotypic model . Also the Spanish clinical classification of COPD patients differs from ours because it simply aims at assessing the perceived prevalence of predefined clinical phenotypes . Given patients were on regular topical therapy for COPD, the important prevalence of night time symptoms and severely impaired health status testifies to the poor quality of the overall pharmacological approach. Analogously, in a European survey from primary care centers two thirds of COPD patients complained of dyspnea despite current treatment ; an even higher percentage of subjects suffered from chronic cough, and this was similar across all severity COPD stages. Otherwise, the important clinical differences among clusters suggest that many patients are expected to gain the greatest benefit from non- pharmacological measures contrasting phlegm and mucus. Indeed, being recognized as a multidimensional and heterogeneous condition, COPD should be the object of personalized care.
The fact that the SGRQ contributed to categorize COPD patients further confirms that health status is something different from and, thus, may not be merely expressed by symptoms. Interestingly, dyspnea associates with SGRQ in latent factor 2, as if dyspnea were the main correlate of perceived health status. This partly reflects the structure of the SGRQ, which largely relates the perceived status to physical activities which are obviously limited by dyspnea. Furthermore, this finding stresses the need of including an index of health status among the instruments rating COPD severity and the effects of therapeutic interventions.
Limitations of the study deserve consideration. First, we chose the classificatory variables with the aim of limiting recall bias. This is the reason why we selected the Night-time, Morning and Day-time Symptoms of COPD questionnaire, with a recall period of a week. Indeed, we based the classification on variables reflecting the patient’s condition at the time of enrollment in order to prevent any classificatory bias linked to uncertain definition of variables requiring an accurate recall, such as frequency of exacerbations. Second, the definition of a clinical variant or subtype of COPD should be founded on the clinical pattern in the absence of any therapy. However, it is almost impossible to recruit patients with severe COPD free from pharmacological treatment, and those untreated likely represent a poorly representative sample. Third, the observed clinical variants of COPD should not be considered true phenotypes. Instead, phenotypes are intrinsically stable, and only the follow-up will assess cluster stability. Thus, the full spectrum of classificatory properties of our model will emerge only at the end of the follow-up phase of the study, which is currently ongoing. Finally, our patients had mean age of 71.4 years; thus, present results might not fully apply to the growing fraction of COPD patients classified as old (age over 75 years), which pose distinctive diagnostic and therapeutic problems .
Despite these limitations, our study emphasizes the clinical meaning of the circadian rhythm of COPD symptoms. It suggests that a classificatory approach based on the respiratory clinical dimension of COPD succeeds in identifying classes of patients with distinctive features. In an era characterized by a growing interest in extra pulmonary features of COPD, which obviously are worthy of consideration, our study redirects the attention of physicians to the elementary truth that COPD is a primarily respiratory disease and an accurate clinical rating of patient’s status is the first task required to any physician caring for a COPD patient.
S1 Table [doc]
Factor analysis output (eigenvalues of the reduced correlation matrix).
1. Lange P, Halpin DM, O'Donnell DE, MacNee W. Diagnosis, assessment, and phenotyping of COPD: beyond FEV1. Int J Chron Obstruct Pulmon Dis. 2016 Feb 19;11 Spec Iss:3–12.
2. Kumar M, Seeger W, Voswinckel R. Senescence-associated secretory phenotype and its possible role in chronic obstructive pulmonary disease. Am J Respir Cell Mol Biol. 2014 Sep;51(3):323–33. doi: 10.1165/rcmb.2013-0382PS 25171460
3. Price D Proceedings of the IV REVIEW: NIGHT-TIME SYMPTOMS IN COPD A. AGUSTI ET AL. 192 VOLUME 20 NUMBER 121 EUROPEAN RESPIRATORY REVIEW World Asthma and COPD Forum (Paris, France, April 30 to May 3, 2011)
4. Partridge MR, Karlsson N, Small IR. Patient insight into the impact of chronic obstructive pulmonary disease in the morning: an internet survey. Curr Med Res Opin. 2009 Aug;25(8):2043–8. doi: 10.1185/03007990903103006 19569976
5. Miravitlles M, Worth H, Soler Catalauña JJ, Price D, De Benedetto F, Roche N, et al. Observational study to characterize 24-hour COPD symptoms and their relationship with patient-reported outcomes; results fromthe ASSESS study. Respir Res. 2014 Oct;21; 15: 122. doi: 10.1186/s12931-014-0122-1 25331383
6. Price D, Small M, Milligan G, Higgins V, Garcia Gil EG, Estruch J. Impact of night-time symptoms in COPD: a real-world study in five European countries. Int J Chron Obstruct Pulmon Dis. 2013; 8:595–603. doi: 10.2147/COPD.S48570 24348032
7. Bellia V, Catalano F, Scichilone N, Incalzi RA, Spatafora M, Vergani C et al. Sleep disorders in the elderly with and without chronic airflow obstruction: the SARA study. Sleep. 2003 May 1;26(3):318–23. 12749552
8. Canonica GW, Blasi F, Scichilone N, Simoni L, Zullo A, Giovannetti C et al. on behalf of STORICO study group Characterization of circadian COPD symptoms by phenotype: Methodology of the STORICO observational study. Eur J Inter Med. 2017 Sep; 43: 62–68.
9. Jones PW, Quirk FH, Baveystock CM, Littlejohns P. A self-complete measure for chronic airflow limitation. The St George's respiratory questionnaire. Am Rev Respir Dis. 1992 June;145 (6):1321–7. 1595997
10. Meguro M, Barley EA, Spencer S, Jones PW. Development and validation of an improved COPD-specific version of the St George's respiratory questionnaire. Chest. 2007 Aug;132 (2):456–63. 17646240
11. Zigmond AS, Snaith RP. The hospital anxiety and depression scale. Acta Psychiatr Scand 1983;67:361–70. 6880820
12. Pokrzywinski RF, Meads DM, McKenna SP, Glendenning GA, Revicki DA. Development and psychometric assessment of the COPD and Asthma Sleep Impact Scale (CASIS). Health Qual Life Outcomes 2009;7:98. doi: 10.1186/1477-7525-7-98 19968881
13. Ainsworth BE, Bassett DR Jr, Strath SJ, Swartz AM, O'Brien WL, Thompson RW et al. Comparison of three methods for measuring the time spent in physical activity. Med Sci Sports Exerc 2000;32:S457–64. 10993415
14. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A et al. Standardisation of spirometry. Eur Respir J 2005; 26: 319–338 16055882
15. Child D. (1990). The essentials of factor analysis, second edition. London: Cassel Educational Limited.
16. Guttman L. (1953) Image Theory for the Structure of Quantitative Variables, Psychometrica, 18, 277–296.
17. Kaiser H.F., Rice J. Little Jiffy, Mark IV, Educational and Psychological Measurement. 1974; 34 (1), 111–117.
18. Truong Y., McColl R. Intrinsic motivations, self-esteem, and luxury goods consumption. Journal of Retailing and Consumer Services.2011; 18 (6): 555–561
19. Singh D, Miravitlles M, Vogelmeier C. Chronic Obstructive Pulmonary Disease Individualized Therapy: Tailored Approach to Symptom Management. Adv Ther. 2017 Feb;34(2):281–299. doi: 10.1007/s12325-016-0459-6 27981495
20. Van Boven JFM, Lavorini F, Dekhuijzen PNR, Blasi F, Price DB, Viegi G Urging Europe to put non-adherence to inhaled respiratory medication higher on the policy agenda: a report from the First European Congress on Adherence to Therapy. Eur Respir J. 2017 May 19;49(5).
21. Calverley PM, Anderson JA, Celli B, Ferguson GT, Jenkins C, Jones PW et al. Salmeterol and fluticasone propionate and survival in chronic obstructive pulmonary disease. N Engl J Med. 2007 Feb 22;356(8):775–89. 17314337
22. Tashkin DP, Celli B, Senn S, Burkhart D, Kesten S, Menjoge S et al. A 4-year trial of tiotropium in chronic obstructive pulmonary disease. N Engl J Med. 2008 Oct 9;359(15):1543–54. doi: 10.1056/NEJMoa0805800 18836213
23. Kessler R, Partridge MR, Miravitlles M, Cazzola M, Vogelmeier C, Leynaud D et al. Symptom variability in patients with severe COPD: a pan-European cross-sectional study. Eur Respir J. 2011 Feb;37 (2):264–72. doi: 10.1183/09031936.00051110 21115606
24. Partridge MR, Karlsson N, Small IR. Patient insight into the impact of chronic obstructive pulmonary disease in the morning: an internet survey. Curr Med Res Opin. 2009 Aug;25 (8):2043–8. doi: 10.1185/03007990903103006 19569976
25. Burgel PR, Paillasseur J-L, Janssens W, Piquet J, ter Riet G, Garcia-Aymerich J et al A simple algorithm for the identification of clinical COPD phenotypes Eur Respir J 2017; 50: 1701034 doi: 10.1183/13993003.01034-2017 29097431
26. Pikoula M, Quint JK, Nissen F, Hemingway H, Smeeth L, Denaxas S Identifying clinically important COPD sub-types using data-driven approaches in primary care population based electronic health records BMC Med Inform Decis Mak. 2019; 19: 86. doi: 10.1186/s12911-019-0805-0 30999919
27. Augustin IML, Spruit MA, Houben-Wilke S, Franssen FME, Vanfleteren LEGW, Gaffron S et al. The respiratory physiome: Clustering based on a comprehensive lung function assessment in patients with COPD. PLoS One. 2018 Sep 12;13(9):e0201593. doi: 10.1371/journal.pone.0201593 30208035
28. Agusti A, Bel E, Thomas M, Vogelmeier C, Brusselle G, Holgate S et al. Treatable traits: toward precision medicine of chronic airway diseases European Respiratory Journal 2016 47: 410–419; doi: 10.1183/13993003.01359-2015 26828055
29. Calle Rubio M, Casamor R, Miravitlles M. Identification and distribution of COPD phenotypes in clinical practice according to Spanish COPD Guidelines: the FENEPOC study. Int J Chron Obstruct Pulmon Dis. 2017 Aug 9;12:2373–2383. doi: 10.2147/COPD.S137872 28848338
30. Jones PW, Brusselle G, Dal Negro RW, Ferrer M, Kardos P, Levy ML et al. Health-related quality of life in patients by COPD severity within primary care in Europe. Respir Med. 2011 Jan;105(1):57–66. doi: 10.1016/j.rmed.2010.09.004 20932736
31. Corsonello A, Scarlata S, Pedone C, Bustacchini S, Fusco S, Zito A et al. Treating COPD in Older and Oldest Old Patients. Curr Pharm Des. 2015;21(13):1672–89. 25633118