Chronic Obstructive Pulmonary Disease Subtypes. Transitions over Time

Abstract

Background

Although subtypes of chronic obstructive pulmonary disease are recognized, it is unknown what happens to these subtypes over time. Our objectives were to assess the stability of cluster-based subtypes in patients with stable disease and explore changes in clusters over 1 year.

Methods

Multiple correspondence and cluster analysis were used to evaluate data collected from 543 stable patients included consecutively from 5 respiratory outpatient clinics.

Results

Four subtypes were identified. Three of them, A, B, and C, had marked respiratory profiles with a continuum in severity of several variables, while the fourth, subtype D, had a more systemic profile with intermediate respiratory disease severity. Subtype A was associated with less dyspnea, better health-related quality of life and lower Charlson comorbidity scores, and subtype C with the most severe dyspnea, and poorer pulmonary function and quality of life, while subtype B was between subtypes A and C. Subtype D had higher rates of hospitalization the previous year, and comorbidities. After 1 year, all clusters remained stable. Generally, patients continued in the same subtype but 28% migrated to another cluster. Together with movement across clusters, patients showed changes in certain characteristics (especially exercise capacity, some variables of pulmonary function and physical activity) and changes in outcomes (quality of life, hospitalization and mortality) depending on the new cluster they belonged to.

Conclusions

Chronic obstructive pulmonary disease clusters remained stable over 1 year. Most patients stayed in their initial subtype cluster, but some moved to another subtype and accordingly had different outcomes.

About the Authors:

Cristóbal Esteban

* E-mail: [email protected]

Affiliations Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain, Health Services Research on Chronic Patients Network (REDISSEC), Galdakao, Bizkaia, Spain

Inmaculada Arostegui

Affiliations Health Services Research on Chronic Patients Network (REDISSEC), Galdakao, Bizkaia, Spain, Department of Applied Mathematics, Statistics and Operative Research, University of the Basque Country (UPV/EHU), Leioa, Bizkaia, Spain, Basque Center for Applied Mathematics (BCAM), Leioa, Bizkaia, Spain

Myriam Aburto

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Javier Moraza

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

José M. Quintana

Affiliations Health Services Research on Chronic Patients Network (REDISSEC), Galdakao, Bizkaia, Spain, Research Unit, Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Amaia García-Loizaga

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Luis V. Basualdo

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Amaia Aramburu

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Susana Aizpiri

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Ane Uranga

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Alberto Capelastegui

Affiliation: Respiratory Department. Hospital Galdakao-Usansolo, Galdakao, Bizkaia, Spain

Introduction

Grouping patients by their characteristics has long been an aspiration in chronic obstructive pulmonary disease (COPD) [1]. In recent years, with stronger evidence of heterogeneity in COPD [2], the level of interest has increased, and this has seen the creation of groups [3], phenotypes [4], and severity scores [5] in an attempt to determine prognosis, improve the process of care and/or move towards precision medicine in the treatment of the disease. One of these approaches has involved identifying characteristics to classify these patients into clusters. The goal of clustering COPD patients into subtypes may be related to prediction of outcomes (prognosis) or to establishing more effective treatment or care strategies. The information available in this area is scarce but important, some studies focusing on stable COPD patients [6–9] and others on exacerbated patients [10,11].

Current information available about COPD phenotypes in stable COPD using clusters has identified several clusters, but there are three main phenotypes that consistently appear in all the different studies. A group with a high level of respiratory symptoms, severe obstruction, lower exercise capacity and lung density, and a low level of severe comorbidities (mainly diabetes and heart disease), this group being labeled the “severe respiratory phenotype” by Garcia-Aymerich et al. [7] The second phenotype, labeled “systemic COPD,” is associated with considerable comorbidities (cardiovascular diseases and diabetes) but moderate-to-severe airway obstruction. The third phenotype had a lower level of dyspnea, a lower rate of comorbidities, and mild-to-moderate airway obstruction. The aforementioned three clusters are in concordance with those described by Burgel et al. [8] and by Rennard et al. [9]. These clusters are related to important outcomes such as hospitalization [7, 9] and mortality [7–9].

To the best of our knowledge, this process has not been taken further, and to do so we need to know the evolution of these COPD clusters over time. Combining multivariate techniques to analyze clinical data obtained in a prospective cohort study could provide validation of the structure of stable COPD phenotypes and transitions between phenotypes over time. The objective of our study was to assess the stability of cluster-based COPD subtypes in stable patients and to identify any changes in the individuals´ subtypes, that is, transitions between COPD clusters, over 1 year of follow-up.

Methods

Participants and data collection

We recruited patients being treated for COPD between January 2003 and January 2004 at one of five outpatient respiratory clinics run by the Respiratory Service of Galdakao-Usansolo Hospital. Patients were consecutively included in the study if they had been diagnosed with COPD for at least 6 months and had been stable for 6 weeks before enrollment. Other inclusion criteria were forced expiratory volume in one second (FEV1) <80% of the predicted value and a FEV1/forced vital capacity ratio <70%.

Patients were not eligible for the study if they had been diagnosed with asthma, they had another important respiratory disease or cancer, or they were suffering from psychiatric or neurological problems that might prevent effective collaboration. The protocol was approved by the Ethics and Research Committees of Hospital Galdakao-Usansolo (03/09062005). Each candidate patient was given detailed information about the study and all those included provided written informed consent.

Study protocol

Sociodemographic variables and smoking habits were recorded. The level of dyspnea was established using modified Medical Research Council scale (mMRC) [12]. Comorbidities were determined by reviewing patients' medical records and summarized in the Charlson comorbidity index [13]. Health-related quality of life (HRQoL) was assessed using the Spanish validated version of the Saint George's Respiratory Questionnaire (SGRQ) [14,15]. The level of physical activity, defined as the time patients spent walking during their leisure time, was classified into four categories: low, moderate, high, and very high corresponding respectively to engaging in light physical activity such as walking for less than 2 hours/week, for 2 to 4 hours/week, more than 4 hours/week and more than 4 hours/week plus playing sports or doing some other type of intense physical activity [16,17]. Complete pulmonary function tests included forced spirometry, bronchodilator testing, and body plethysmography, as well as measurements of carbon monoxide diffusing capacity, and respiratory muscle strength. These tests were performed following the standards of the Spanish Society of Respiratory Medicine and Thoracic Surgery (SEPAR) [18]. For theoretical values, we considered those established by the European Community for Steel and Coal [19].

Two 6-minute walking tests were performed according to American Thoracic Society guidelines [20]. Peripheral muscle strength was evaluated in terms of handgrip strength [21] and extension knee force and shoulder abduction using a hand-held dynamometer [22].

Follow-up

One year after inclusion in the study, survivors were interviewed and repeated the aforementioned assessments. No interventions were performed related to this study, and the research team did not take part in patients' routine treatment that mainly include long acting bronchodilators and inhaled corticosteroids, or the treatment of any exacerbations.

During this 2-year follow-up period, patient medical records and the hospital database on hospitalizations were reviewed. Vital status was established by reviewing medical records, examining the hospital database and public death registries. Deaths were considered confirmed if the name, sex, and date of birth on the record matched those of the participant.

Statistical analysis

We combined multiple correspondence analysis (MCA) and cluster analysis to analyze clinical data on all participants included in the study [23,24]. Appendix in the Supporting Information shows a brief description on MCA and cluster analysis, including interpretation of the results (S1 File). MCA was performed twice: first, baseline characteristics were used to define the phenotype structure, and second, data for the same cohort 1 year later were added to the baseline data in order to identify any transitions between phenotypes. Variables included in the MCA, along with the corresponding categories, are listed in the Online Data Supplement, (S1 Table). Cluster analysis was performed using the relative position of the categories given by the MCA factors, the health components initially provided by the MCA being used to identify groups of individuals [25]. The number of groups was selected based on minimum loss of inertia [26]. Participants were classified using the health components at baseline, identifying cluster-based COPD phenotypes. Reclassification after 1 year was performed with both baseline and follow-up health components, showing the transition of participants across the COPD groups during the follow-up period. Phenotype structure stability was assessed by comparing the individual results from the two time points. Phenotypes were validated by assessing their relationships with key outcomes such as mortality, hospitalization, quality of life and COPD multidimensional severity scores. Association of cluster with outcomes was performed by chi-squared test for mortality, hospitalizations and GOLD classification; whereas analysis of variance was performed for HRQoL and severity scores (HADO-AH, BODE and ADO). Kaplan-Meier survival curves stratified by clusters were calculated. All statistical analyses were performed using R v2.13.0 [27].

Results

The cohort included 543 patients, of whom 96% were male. The cohort had a mean age of 68.3 (SD ± 8.3) years, body mass index of 28.3 (SD ± 4.4) kg/m², mMRC dyspnea score of 2.4 (SD ± 0.9), post-bronchodilator FEV1% of 55% (SD ± 13.3) of the predicted value, and 6-minute walking distance of 408.9 m (SD ± 92.4).

MCA analysis

Results from the MCA showed that 73.1% of the variance in the original variables was explained by four factors when analysis was performed at baseline and a similar figure of 73.3% was obtained for analysis of the variables at 1 year. Results were almost identical for the two periods of time, showing the high stability of the variables that contribute to the loading factors obtained from the MCA. Patient groupings can be shown graphically considering the two main factors (Fig 1). Categories that significantly contributed to the first factor were age < 65 years, FEV1% ≥ 50%, hand strength ≥ 40 kg, 6-minute walking distance ≥ 425 m, a very high level of physical activity and dyspnea score of 0–1 on the negative side; and FEV1% < 30%, hand strength < 30 kg, 6-minute walking distance < 350 m, a low-to-moderate level of physical activity and dyspnea score ≥ 2 on the positive side. Therefore, we have labeled this first factor “respiratory condition” (graphically represented by the horizontal axis, going from a mild impairment on the negative part of the axis to a severe impairment on the positive part).

[Figure omitted. See PDF.]

Fig 1. Map of clusters and distribution of patients.

Map created by the first and second components derived from the MCA is shown at the center. Four circles at the sides show how patients move between clusters after one year of follow-up. Relative positions of the subjects in the planes are represented by different colors, depending on the subtype provided by the cluster analysis. Definition of the axes is suggested based on information provided in appendix Table A1. The horizontal axis, first component, can be defined as an index of the respiratory conditions of the patient, milder (left side) vs. more severe (right side). The vertical axis, second component, can be defined as an index of the systemic status, worse (bottom) vs. better (top).

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

The second factor was based mainly on body mass index ≥ 30 kg/m², Charlson comorbidity score > 3 and presence of comorbidities on the negative side and by body mass index < 25 kg/m², Charlson comorbidity score 0–1 and no comorbidities on the positive side. Therefore, we have labeled the second factor “systemic condition”, basically related to number of comorbidities (graphically represented by the vertical axis, with presence of comorbidities on the negative part of the axis and absence of comorbidities on the positive part). Contributions from each variable to the first and second health components is shown in detail in the Supporting Information material (S1 Table). The third and fourth factors were a mixture of many of the original variables, and hence we decided not to give them a specific label.

Clustering

In the clustering process (Fig 2), patients were first classified based equally on the respiratory and systemic conditions (first and second factors). At that point, it was the comorbidity status that helped discriminate between two subtypes of patients (second factor). Subsequently, participants with fewer comorbidities were split again into two clusters based on their respiratory impairment (first factor). Third and fourth factors contribute at all three splitting points, although to a lesser extent. The four subtypes (A, B, C, and D) provide a typology of stable COPD patients.

[Figure omitted. See PDF.]

Fig 2. Partial dendrogram obtained from the cluster analysis.

The dendogram represents the results from the cluster analysis performed with the four components obtained from the multiple correspondence analysis. The graphical display includes an easy interpretation of the partition and a brief description of the resulting clusters.

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

Subtype A included younger participants, with less dyspnea, higher diffusing capacity of the lung for carbon monoxide (DLCO), muscle strength, exercise capacity, and level of physical activity, better HRQoL and lower Charlson comorbidity scores (Tables 1 and 2). Those with subtype C had the heaviest smoking history, highest dyspnea scores, poorest pulmonary function and HRQoL, and lowest muscle strength, exercise capacity, and level of physical activity. However, this subtype was associated with a similar Charlson comorbidity score to subtype A. In relation to the aforementioned variables, patients in subtype B are between subtypes A and C, and had very similar characteristics to subtype D regarding respiratory condition. However, apart from more overall comorbidities, patients in subtype D had a higher rate of hospitalization the year before, and lower exercise capacity and level of physical activity than those in subtype B (p<0.05). As regard as treatment, the only statistically significant difference between clusters was that cluster A shows a lower use of long-acting bronchodilator medication (LABA and/or LAMA) than the rest of the clusters. Our cluster classification matched the Global Strategy for Diagnosis, Management, and Prevention of COPD (GOLD) classification for subtypes A and C. In contrast, participants with subtype D were distributed in similar proportions across every GOLD category, and those with subtype B were also spread across the GOLD categories (though somewhat less evenly). With respect to the relationship of our clusters and multidimensional severity scales, the cluster classification showed a good association with all the multidimensional scales studied. However, the BODE index allowed to establish a different prognosis in every cluster, whereas neither the ADO nor the HADO-AH were able to differentiate between the cluster B and D prognosis (Table 3).

[Figure omitted. See PDF.]

Table 1. Distribution of the main clinical and functional variables related to clusters at baseline.

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

[Figure omitted. See PDF.]

Table 2. Distribution of the comorbidity variables related to clusters at baseline.

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

[Figure omitted. See PDF.]

Table 3. Association of cluster-based classification to GOLD classification and COPD severity score (HADO-AH, BODE and ADO).

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

Clusters and patients transition and outcomes

Performing cluster analysis for the two periods 1 year apart showed stability of the initial subtypes. When this analysis was performed first with baseline health components and later with both, baseline and follow-up assessments, we observed that movement of participants from one group to another was associated with changes in variables related to their respiratory condition (pulmonary function, exercise capacity, and physical activity) and also, to a lesser extent, the fact that baseline variables placed them on the border between two subtypes and hence small changes in some of variables analyzed moved them from one cluster to another. Transition of patients between clusters over time is shown graphically in Fig 1 and numerically in Table 4.

[Figure omitted. See PDF.]

Table 4. Mortality rate, hospitalizations and HRQoL for patients with COPD exacerbation based on the clusters at baseline (1 year period) and based on the clusters transition for a 1 year period.

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

More complete information about characteristics of patients in the clusters over time is provided in the Supporting Information material (S2–S5 Tables).

Regarding the relation of the four clusters with specific outcomes, significant differences were found for mortality rate, HRQoL scores and hospitalization due to COPD (Table 4). On inclusion in the study, HRQoL was poorest for subtype C, and best for subtype A, with subtypes B and D showing intermediate impairment with similar scores, but significantly different from the two other subtypes. Mortality and hospitalization rates during the first year of follow-up tended to grow from subtype A to subtype D. Regarding as transition of patients between clusters at 1 year, in general those who remained in the same cluster maintained their HRQoL score, whereas those who changed subtype had better or poorer HRQoL depending on their new subtype. Similar trend in results was observed regarding as mortality and hospitalization rates during the next year. The clearest example was patients that moved from subtype D to C with significantly higher mortality and hospitalization rates of 40% and 60% and significantly poorer HRQoL than those remaining at cluster D (Table 4). Moreover, significant differences between survival curves were observed between clusters A and C; A and D; and B and D (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Kaplan-Meier estimate of the survival function during the one year period of follow-up stratified by cluster.

Log-rank test for homogeneity (p < 0.001). Significant differences adjusted for multiple comparisons (Bonferroni) were observed between clusters A and C; A and D; and B and D.

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

Discussion

Knowing that in previous important studies [7–9] a longitudinal follow-up was carried out, to our knowledge, this is the first time that temporal evolution of COPD patient clusters has been analyzed, and the following were our main findings. 1) The general characteristics of the clusters remained stable over 1 year. 2) Most patients stayed in their initial cluster, but there was a migration of patients between clusters A, B and C. In general, this migration was towards a more severe cluster, although some patients moved to clusters indicating an improvement in their clinical condition. 3) Our cluster classification is in concordance with previously established COPD clusters. 4) Our clusters showed a clear relationship with different relevant COPD outcomes, and patient migration between clusters was associated with changes in outcomes.

Generally speaking, our four subgroups are based on characteristics that are in concordance with previous cluster analysis studies [7–9]. The identification of clusters with very similar characteristics from different cohorts [7,9] consolidates these subtypes for COPD patients, despite the fact that in one study a considerable number of patients (more than 50%) were not included in any specific cluster [9]. Of the four subtypes we have identified, three (A, B, and C) are suitable for describing a growing severity of the respiratory disease (going from mild to severe) and its expression in other spheres like age, muscle strength, physical activity performed, exercise capacity and HRQoL. Our fourth subtype (D) was very similar to the intermediate subtype of the respiratory group (B) from the point of view of severity of the respiratory disease but with a significantly higher rate of comorbidities. Specifically, the factors that differentiated this “more systemic” subtype from the “intermediate respiratory” subtype are a history of more hospitalizations (i.e., higher rate in the previous year) for COPD exacerbation (higher in the systemic subtype), and the distance covered in the 6-minute walking test (lower in the systemic subtype), as well as, evidently, the rate of comorbidities.

The distribution of our subtypes in relation to the GOLD categories shows that the “mild respiratory” subtype (A) were placed mainly in GOLD group A, and the “severe respiratory” subtype (C) mainly in GOLD group D. On the other hand, the “intermediate respiratory” subtype (B) showed a slight tendency to be placed in GOLD group A, but like our “more systemic” subgroup (D) these patients were distributed fairly evenly across all GOLD groups. This could be interpreted as the GOLD combined assessment not being able to reliably categorize COPD patients with an intermediate level of respiratory disease. Other authors found patients with higher rates of comorbidities to be GOLD groups B or D [28], but this was not the case for our patients. In fact, as mentioned above, patients from subtype D were uniformly distributed across the four GOLD groups. Therefore, it seems that the GOLD categorization does not properly take into account the importance of comorbidities.

We did not find an association between the level of impairment of pulmonary function and the rate of comorbidities, unlike in a previous study [29], patients with the highest level of comorbidities in our cohort being mainly grouped in a cluster with an intermediate level of pulmonary deterioration. When we studied the relationship between our cluster classification and three multidimensional severity scales, we found it to be most strongly correlated with the BODE index, [5], in fact it was the only multidimensional severity scale that was able to capture significant differences between clusters B and D.

There is growing evidence of the importance of comorbidities in COPD and this is an issue that should be considered in the evaluation of patients with this disease. In fact, comorbidity-based clusters have been described [30], but we need to know the impact of these different clusters on outcomes. Moreover, comorbidities have been included in predictive mortality models and scores in COPD [31] and were able to improve the predictive capacity of the BODE index [32].

Our subtype cluster classification is supported by its relationship with powerful outcomes namely mortality and hospitalization for COPD exacerbation. The “more systemic” subtype and the “respiratory severe” subtype showed the highest mortality rates during the 1 year of follow-up after inclusion in the study, of 16.8% and 12.4% respectively.

Reflecting on the previously described patterns, what is particularly interesting in our study is that cluster assignment was not static over the 1-year follow-up. Rather, although most patients stayed in their initial cluster, notable movement was seen across the three respiratory clusters (subtypes A, B, and C). During the 1 year follow-up period, 134 patients (28%) moved to another cluster. Specifically, 81 (17%) migrated to a more severe cluster and 53 (11%) to a less severe cluster, that is, far from negligible movement. Those movements were associated with changes in variables of pulmonary function (FEV1 and diffusion lung capacity for carbon monoxide), exercise capacity, muscle strength and physical activity, and it is plausible that this reflects the concept of disease activity, which would also be changeable over time. Previously, Casanova et al. reported that 27% of patients of a cohort had statistically significant changes in FEV1% and/or BODE index over a year [33], though these changes only reflected deterioration in those variables. In turn, in our study, patients who migrated to other clusters changed their initial prognosis in terms of mortality and hospitalization to those of the cluster to which they had moved.

Twenty-six patients of the “more systemic” cluster (D) migrated to the “intermediate (B) and severe (C) respiratory” clusters maintaining a high level of comorbidity, and this could be interpreted as a weakness of the model. However, those who the model placed in the more “severe respiratory” cluster (C) experienced a very substantial impairment in pulmonary function and exercise capacity over the 1-year follow-up. These changes would explain the change in cluster, as well as their hospitalizations and mortality rates shooting up. In other words those initially in the “more systemic” cluster (D) who move to the “severe respiratory” cluster (C) show a marked worsening in their prognosis. On the other hand, those who migrate from the “more systemic” subtype (D) to the “intermediate respiratory” cluster (B) showed an improvement in exercise capacity and physical activity and, while Charlson comorbidity scores remained steady, the specific comorbidities were very different, those who stayed in the “more systemic” cluster (D) having a higher level of cardio-vascular comorbidity. This pattern could mean that not all comorbidities have the same importance and probably a ranking of impact on prognosis should be established for different comorbidities.

In general terms, some variables repeatedly appeared to be associated with migration between clusters, whatever the direction of the change (better or worse). This could be because changes in 6-minute walking distance over the year were associated with outcomes such as mortality independently of changes in FEV1 [34], while changes in physical activity have been related to hospitalizations and changes in HRQoL [35,36], and muscle strength predicted mortality in a cohort with moderate-to-severe COPD [37]. However, the design of our study does not enable us to establish whether changes in these variables are the cause or the consequence of the migration of patients between clusters.

This study has several limitations. The clinical condition of the three respiratory clusters (A, B, and C) must be understood as a continuum. Although this continuum could be initially attributable to age, in our opinion there are other patients’ conditions that weigh more than age. This assertion is based on the fact that several conditions, not only respiratory, but also muscle strength, physical activity performed or exercise capacity had higher contribution to the first factor than age. Moreover, this continuum means that some patients lie on the border between clusters. Accordingly, small changes in their clinical condition would partially explain changes in classification, although this is not the main reason for the migration between clusters. Further, given the relatively small total number of patients for whom we had complete data for the second clustering process, relatively small numbers of patients migrated to other clusters and hence some results showed interesting tendencies but were not statistically significant. Our cohort had a small proportion of women, because the smoking trend of women in our country, there after the generalizability of the result to other cohort could be controversial. Getting a CT scan of every patient would have enriched our clustering process. In its place we have used the CO transfer as a subrogated marker of emphysema. In fact the CO transfer was one of the variables most frequently involved in the clustering and in the changes of the patient over the follow-up.

On the other hand, the types of analysis performed are a strength of the study. Multivariate techniques combined with cluster analysis have been widely used to differentiate groups of individuals, including for COPD [6–9]. Cluster analysis has previously been used in the literature to classify individuals into groups and also to assess phenotypes in patients with airway diseases [38,39]. Some studies in patients with airway diseases have combined factor and cluster analysis [5]. A strong assumption of factor analysis is that variables included must be continuous and normally distributed. MCA provides an alternative to factor analysis for categorical variables, without any distributional assumptions, and hence was more suitable for our study, that includes patient characteristics that were dichotomous (yes/no) or had three or four discrete possibilities, such as FEV1%, physical activity, and mMRC dyspnea score. MCA has been used previously with similar objectives to ours, including in relation to COPD, when continuous and qualitative explanatory variables were both candidate variables for a model [11]. Thus, we used a methodology previously shown to be useful for the creation of subtypes of individuals and describe their typology for transforming the information contained in the original categorical variables into two factors, which were interpreted as disease components.

In summary, our data allow us to demonstrate the stability of COPD clusters over time and understand the course of COPD patients over a short period of time. During this period, most patients continued to show the same phenotype, but some of them changed their situation, improving or worsening, their condition and therefore their outcomes. These changes should make us consider that we have to improve the characterization of COPD patients and identify these “moving” patients to try to avoid or slow progression of the disease and improve patients’ condition.

Supporting Information

[Figure omitted. See PDF.]

S1 File. Multiple Correspondence Analysis (MCA) and Cluster Analysis (CA).

https://doi.org/10.1371/journal.pone.0161710.s001

(DOCX)

S1 Table. Absolute contribution of each variable by category to the first and second factors of the MCA.

¹ Indicates the contribution of each category/ level of the variable to that factor or health component, higher value indicates higher contribution. ² Sign indicates the side of the map where each category of each variable is located (referring to Fig 1). The most important contributions along with their signs are highlighted.

https://doi.org/10.1371/journal.pone.0161710.s002

(DOCX)

S2 Table. Distribution of the main variables related to patient’s COPD at baseline for patients in cluster A (n = 164) and evolution in a one year period including cluster transition.

*5 patients (3%) were lost during the follow-up. Confidence intervals for groups of less than 5 individuals were not calculated (showed as --). Mean (95% CI) for continuous variables and n (%) for categorical variables. Dyspnea (mMRC): modified Medical Research Council Dyspnea Scale. RV: residual volume. DLCO: diffusion lung capacity for carbon monoxide. VA: alveolar volume.

https://doi.org/10.1371/journal.pone.0161710.s003

(DOCX)

S3 Table. Distribution of the main variables related to patient’s COPD at baseline for patients in cluster B (n = 195) and evolution in a one year period including cluster transition.

*10 patients (5.1%) were lost during the follow-up. Confidence intervals for groups of less than 5 individuals were not calculated (showed as --). Mean (95% CI) for continuous variables and n (%) for categorical variables. Dyspnea (mMRC): modified Medical Research Council Dyspnea Scale. RV: residual volume. DLCO: diffusion lung capacity for carbon monoxide. VA: alveolar volume.

https://doi.org/10.1371/journal.pone.0161710.s004

(DOCX)

S4 Table. Distribution of the main variables related to patient’s COPD at baseline for patients in cluster C (n = 89) and evolution in a one year period including cluster transition.

*7 patients (7.9%) were lost during the follow-up. Confidence intervals for groups of less than 5 individuals were not calculated (showed as --). Mean (95% CI) for continuous variables and n (%) for categorical variables. Dyspnea (mMRC): modified Medical Research Council Dyspnea Scale. RV: residual volume. DLCO: diffusion lung capacity for carbon monoxide. VA: alveolar volume.

https://doi.org/10.1371/journal.pone.0161710.s005

(DOCX)

S5 Table. Distribution of the main variables related to patient’s COPD at baseline for patients in cluster D (n = 95) and evolution in a one year period including cluster transition.

*4 patients (4.2%) were lost during the follow-up. Confidence intervals for groups of less than 5 individuals were not calculated (showed as --). Mean (95% CI) for continuous variables and n (%) for categorical variables. Dyspnea (mMRC): modified Medical Research Council Dyspnea Scale. RV: residual volume. DLCO: diffusion lung capacity for carbon monoxide. VA: alveolar volume.

https://doi.org/10.1371/journal.pone.0161710.s006

(DOCX)

Acknowledgments

We thank Julia G Fenn for editing the manuscript and Ministerio de Economía y Competitividad del Gobierno de España (MTM2013-40941-P), Departamento de Educación, Política Lingüística y Cultura del Gobierno Vasco (IT620-13)

Author Contributions

1. Conceptualization: CE IA MA JM JMQ.

2. Formal analysis: CE IA JMQ.

3. Funding acquisition: CE IA AC.

4. Investigation: MA JM AG-L AA LVB SA AU.

5. Methodology: CE IA JMQ.

6. Project administration: CE.

7. Supervision: CE.

8. Visualization: CE IA JMQ.

9. Writing – original draft: CE IA JMQ.

10. Writing – review & editing: CE IA MA JM JMQ AG-L AA LVB SA AU AC.

Citation: Esteban C, Arostegui I, Aburto M, Moraza J, Quintana JM, García-Loizaga A, et al. (2016) Chronic Obstructive Pulmonary Disease Subtypes. Transitions over Time. PLoS ONE 11(9): e0161710. https://doi.org/10.1371/journal.pone.0161710

References

1. Dornhorst AC. Respiratory insufficiency (Frederick Price Memorial Lecture). Lancet 1955; 1: 1185–1187.

2. Agusti A, Calverley PM, Celli B, Coxson HO, Edwards LD, Lomas DA, et al; Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) investigators. Characterisation of COPD heterogeneity in the ECLIPSE cohort. Respir Res 2010;11:122. pmid:20831787

3. Vestbo J, Hurd SS, Agustí AG, Jones C, Vogelmeier A, Anzueto P, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med 2013; 187: 347–365. pmid:22878278

4. Miravitlles M, Soler-Cataluña JJ, Calle M, Soriano JB. Treatment of COPD by clinical phenotypes: putting old evidence into clinical practice. Eur Respir J 2013; 41:1252–1256. pmid:23060631

5. Celli BC, Cote CG, Marín JM, Casanova C, Montes de Oca M, Mendez RA, et al. The body-mass index, airflow obstruction, dyspnea, and exercise capacity index in chronic obstructive pulmonary disease. N Engl J Med 2004; 350: 1005–1012. pmid:14999112

6. Burgel PR, Paillasseur JL, Caillaud D, Tillie-Leblond I, Chanez P, Escamilla R, et al. Clinical COPD phenotypes: a novel approach using principal component and cluster analyses. Eur Respir J 2010; 36: 531–539. pmid:20075045

7. Garcia-Aymerich J, Gómez FP, Benet M, Farrero E, Basagaña X, Gayete A, et al. Identification and prospective validation of clinically relevant chronic obstructive pulmonary disease (COPD) subtypes. Thorax 2011; 66: 430–437. pmid:21177668

8. Burgel PR, Paillasseur JL, Peene B, Dusser D, Roche N, Coolen J, et al. Two Distinct Chronic Obstructive Pulmonary Disease (COPD) Phenotypes Are Associated with High Risk of Mortality. PLOS One 2012; 7:e51048. pmid:23236428

9. Rennard SI, Locantore N, Delafont B, Tal-Singer R, Silverman EK, Vestbo J, et al; Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints. Identification of Five Chronic Obstructive Pulmonary Disease Subgroups with Different Prognoses in the ECLIPSE Cohort Using Cluster Analysis. Ann Am Thorac Soc 2015; 12: 303–312. pmid:25642832

10. Bafadhel M, McKenna S, Terry S, Mistry V, Reid C, Haldar P, et al. Acute exacerbations of chronic obstructive pulmonary disease: identification of biologic clusters and their biomarkers. Am J Respir Crit Care Med 2011; 184: 662–671. pmid:21680942

11. Arostegui I, Esteban C, García-Gutierrez S, Bare M, Fernández-de-Larrea N, Briones E, et al; IRYSS-COPD Group. Subtypes of Patients Experiencing Exacerbations of COPD and Associations with Outcomes. PLoS ONE 2014; 9(6): e98580. pmid:24892936

12. Mahler DA, Ward J, Waterman LA, McCusker C, Zuwallack R, Baird JC. Patient-reported dyspnea in COPD reliability and association with stage of disease. Chest 2009; 136: 1473–1479. pmid:19696126

13. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis 1987; 40: 373–383. pmid:3558716

14. Jones PW, Quirk FH, Baveystock CM. The St George's Respiratory Questionnaire. Respir Med 1991: 85(Suppl B): 25–31. pmid:1759018

15. Ferrer M, Alonso J, Prieto L, Plaza V, Monsó E, Marrades R, et al. Validity and reliability of the St George's Respiratory Questionnaire after adaptation to a different language and culture: the Spanish example. Eur Respir J 1996; 9: 1160–1166. pmid:8804932

16. Garcia-Aymerich J, Lange P, Benet M, Schnohr P, Antó JM. Regular physical activity reduces hospital admission and mortality in chronic obstructive pulmonary disease: a population based cohort study. Thorax 2006;61: 772–778. pmid:16738033

17. Schnohr P, Scharling H, Jensen JS. Changes in leisure time physical activity and risk of death: An observational study of 7000 men and women. Am J Epidemiol 2003;158: 639–644. pmid:14507599

18. Manual SEPAR de Procedimientos. Módulo 4. Procedimientos de evaluación de la función pulmonar-II. Barcelona, Spain: Publicaciones Permanyer; 2004.

19. Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, Peslin R, Yernault JC. Lung volumes and forced ventilatory flows. Report working party standardization of lung function test, European Community for Steel and Coal. Official statement of the European Respiratory Society. Eur Respir J 1993;16:5–40.

20. ATS statement: guidelines for the six-minute walk test. Am J Respir Crit Care Med 2002; 166:111–117. pmid:12091180

21. Mathiowetz V, Kashman N, Volland G, Weber K, Dowe M, Rogers S. Grip and pinch strength: normative data for adults. Arch Phys Med Rehabil 1985; 66:69–74. pmid:3970660

22. Bohannon RW. Responsiveness of measurements of knee extension force obtained by hand-held dynamometry. Isokinet Exerc Sci 2009;17:169–172.

23. Benzécri JP. Statistical analysis as a tool to make patterns emerge from data. In: Watanabe S. Methodologies of Pattern Recognition. New York: Academic Press. 1969; 35–74.

24. Greenacre M. Correspondence analysis in medical research. Stat Methods Med Res 1992; 1: 97–117. pmid:1341654

25. Jambu M. Classification Automatique pour l’Analyse des Donne´s. Tome1: Méthodes et Algoritmes. Paris: Dunod; 1978.

26. Ward JH. Hierarchical grouping to optimize an objective function. J Am Stat Assoc 1963; 58: 236–244.

27. R development Core Team. R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2010; Available: http://cran.r-project.org/.

28. Lange P, Marott JL, Vestbo J, Olsen KR, Ingebrigtsen TS, Dahl M, et al. Prediction of the clinical course of chronic obstructive pulmonary disease, using the new GOLD classification: a study of the general population. Am J Respir Crit Care Med 2012;186: 975–981. pmid:22997207

29. Mannino DM, Thorn D, Swensen A, Holguin F. Prevalence and outcomes of diabetes, hypertension and cardiovascular disease in COPD. Eur Respir J 2008; 32: 962–969. pmid:18579551

30. Vanfleteren LE, Spruit MA, Groenen M, Gaffron S, van Empel VP, Bruijnzeel PL, et al. Clusters of comorbidities based on validated objective measurements and systemic inflammation in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med 2013;187: 728–735 pmid:23392440

31. Almagro P, Cabrera FJ, Diez J, Boixeda R, Alonso-Ortiz MB, Barreiro B, et al. Comorbidities and short-term prognosis in patients hospitalized for acute exacerbation of COPD: the EPOC en Servicios de medicina interna (ESMI) study. Chest 2012;142: 1126–1133. pmid:23303399

32. Divo M, Cote C, de Torres JP, Casanova C, Marin JM, Pinto-Plata V, et al. BODE Collaborative Group. Comorbidities and risk of mortality in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med 2012;186: 155–161. pmid:22561964

33. Casanova C, de Torres JP, Aguirre-Jaíme A, Pinto-Plata V, Marin JM, Cordoba E, et al. The progression of chronic obstructive pulmonary disease is heterogeneous: the experience of the BODE cohort. Am J Respir Crit Care Med. 2011;184:1015–1021. pmid:21836135

34. Pinto-Plata VM, Cote C, Cabral H, Taylor J, Celli BR. The 6-min walk distance: change over time and value as a predictor of survival in severe COPD. Eur Respir J 2004; 23: 28–33. pmid:14738227

35. Esteban C, Quintana JM, Aburto M, Moraza J, Egurrola M, Pérez-Izquierdo J, et al. Impact of changes in physical activity on health-related quality of life among patients with COPD. Eur Respir J 2010; 36: 292–300. pmid:20075059

36. Esteban C, Arostegui I, Aburto M, Moraza J, Quintana JM, Aizpiri S, et al. Influence of changes in physical activity on frequency of hospitalization in chronic obstructive pulmonary disease. Respirology 2014;19: 330–338. pmid:24483954

37. Swallow EB, Reyes D, Hopkinson NS, Man WD, Porcher R, Cetti EJ, et al. Quadriceps strength predicts mortality in patients with moderate to severe chronic obstructive pulmonary disease. Thorax 2007;62:115–120. pmid:17090575

38. 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. pmid:18480428

39. Weatherall M, Travers J, Shirtcliffe PM, Marsh SE, Williams MV, Nowitz MR, et al. Distinct clinical phenotypes of airways disease defined by cluster analysis. Eur Respir J 2009; 34: 812–818. pmid:19357143