Skip to main content

Main menu

  • Home
  • Current issue
  • ERJ Early View
  • Past issues
  • ERS Guidelines
  • Authors/reviewers
    • Instructions for authors
    • Submit a manuscript
    • Open access
    • Peer reviewer login
  • Alerts
  • Subscriptions
  • ERS Publications
    • European Respiratory Journal
    • ERJ Open Research
    • European Respiratory Review
    • Breathe
    • ERS Books
    • ERS publications home

User menu

  • Log in
  • Subscribe
  • Contact Us
  • My Cart
  • Log out

Search

  • Advanced search
  • ERS Publications
    • European Respiratory Journal
    • ERJ Open Research
    • European Respiratory Review
    • Breathe
    • ERS Books
    • ERS publications home

Login

European Respiratory Society

Advanced Search

  • Home
  • Current issue
  • ERJ Early View
  • Past issues
  • ERS Guidelines
  • Authors/reviewers
    • Instructions for authors
    • Submit a manuscript
    • Open access
    • Peer reviewer login
  • Alerts
  • Subscriptions

A simple algorithm for the identification of clinical COPD phenotypes

Pierre-Régis Burgel, Jean-Louis Paillasseur, Wim Janssens, Jacques Piquet, Gerben ter Riet, Judith Garcia-Aymerich, Borja Cosio, Per Bakke, Milo A. Puhan, Arnulf Langhammer, Inmaculada Alfageme, Pere Almagro, Julio Ancochea, Bartolome R. Celli, Ciro Casanova, Juan P. de-Torres, Marc Decramer, Andrés Echazarreta, Cristobal Esteban, Rosa Mar Gomez Punter, MeiLan K. Han, Ane Johannessen, Bernhard Kaiser, Bernd Lamprecht, Peter Lange, Linda Leivseth, Jose M. Marin, Francis Martin, Pablo Martinez-Camblor, Marc Miravitlles, Toru Oga, Ana Sofia Ramírez, Don D. Sin, Patricia Sobradillo, Juan J. Soler-Cataluña, Alice M. Turner, Francisco Javier Verdu Rivera, Joan B. Soriano, Nicolas Roche on behalf of Initiatives BPCO, EABPCO, Leuven and 3CIA study groups
European Respiratory Journal 2017 50: 1701034; DOI: 10.1183/13993003.01034-2017
Pierre-Régis Burgel
1University Paris Descartes (EA2511), Sorbonne Paris Cité, Paris, France
2Dept of Respiratory Medicine, Cochin Hospital, AP-HP, Paris, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: pierre-regis.burgel@cch.aphp.fr
Jean-Louis Paillasseur
3Effi-Stat, Paris, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wim Janssens
4Respiratory Division, University Hospital Gasthuisberg, K.U. Leuven, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jacques Piquet
5Dept of Respiratory Medicine, Le Raincy-Montfermeil Hospital, Montfermeil, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gerben ter Riet
6Dept General Practice – Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Judith Garcia-Aymerich
7ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Universitat Pompeu Fabra (UPF), CIBER Epidemiología y Salud Pública (CIBERESP), Barcelona, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Borja Cosio
8Unidad de Investigación, Servicio de Neumología, Hospital Universitario Son Espases, Palma de Mallorca, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Borja Cosio
Per Bakke
9Dept of Clinical Science, Faculty of Medicine and Dentistry, University of Bergen, Bergen, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Milo A. Puhan
10Epidemiology, Biostatistics und Prevention Institute (EBPI), University of Zurich, Zurich, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Arnulf Langhammer
11Dept of Public Health and General Practice, HUNT Research Centre, Norwegian University of Science and Technology, Levanger, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Inmaculada Alfageme
12Universidad de Sevilla, Seville, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pere Almagro
13Internal Medicine, Hospital Universitari Mutua de Terrassa, Universitat de Barcelona, Barcelona, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julio Ancochea
14Pneumology Service, La Princesa Institute for Health Research (IP), Hospital Universitario de la Princesa, Madrid, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bartolome R. Celli
15Brigham and Women's Hospital, Boston, MA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ciro Casanova
16Hospital Nuestra Señora de la Candelaria, Tenerife, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Juan P. de-Torres
17Clınica Universidad de Navarra, Pamplona, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marc Decramer
4Respiratory Division, University Hospital Gasthuisberg, K.U. Leuven, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrés Echazarreta
18Servicio de Neumonología Hospital San Juan de Dios de La Plata, Buenos Aires, Argentina
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cristobal Esteban
19Hospital Galdakao-Usansolo, Galdakao, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rosa Mar Gomez Punter
20Servicio de Neumología, Hospital Universitario La Princesa, Madrid, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
MeiLan K. Han
21University of Michigan, Ann Arbor, MI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ane Johannessen
22Centre for Clinical Research, Haukeland University Hospital, Bergen, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bernhard Kaiser
23Dept of Pulmonary Medicine, Paracelsus Medical University Hospital, Salzburg, Austria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bernd Lamprecht
24Dept of Pulmonary Medicine, General Hospital Linz (AKH), Linz, Austria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter Lange
25Section of Social Medicine, Dept of Public Health, Copenhagen University, Copenhagen, Denmark
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Linda Leivseth
26Centre for Clinical Documentation and Evaluation, Northern Norway Regional Health Authority, Tromso, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jose M. Marin
27Hospital Universitario Miguel Servet, Zaragoza, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jose M. Marin
Francis Martin
28Pneumologie, Centre Hospitalier de Compiègne, Compiègne, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pablo Martinez-Camblor
29Hospital Universitario Central de Asturias (HUCA), Oviedo, Spain
30Universidad Autónoma de Chile, San Miguel, Chile
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marc Miravitlles
31Pneumology Dept, Hospital Universitary Vall d'Hebron. CIBER de Enfermedades Respiratorias (CIBERES), Barcelona, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marc Miravitlles
Toru Oga
32Dept of Respiratory Care and Sleep Control Medicine, Graduate School of Medicine, Kyoto University, Kyoto, Japan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ana Sofia Ramírez
33Facultad de Medicina UASLP, San Luis Potosí, México
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ana Sofia Ramírez
Don D. Sin
34James Hogg Research Centre, University of British Columbia; Division of Respiratory Medicine, Dept of Medicine, St Paul's Hospital, Vancouver, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patricia Sobradillo
35Hospital Universitario Araba, Sede Txagorritxu, Vitoria, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Juan J. Soler-Cataluña
36Servicio de Neumología, Hospital Arnau de Vilanova, Valencia, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alice M. Turner
37Queen Elizabeth Hospital Research Laboratories, Birmingham, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Francisco Javier Verdu Rivera
38H.U. Son Espases, Palma de Mallorca, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joan B. Soriano
39Instituto de Investigación Hospital Universitario de la Princesa (IISP), Universidad Autónoma de Madrid, Madrid, Spain
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Joan B. Soriano
Nicolas Roche
1University Paris Descartes (EA2511), Sorbonne Paris Cité, Paris, France
2Dept of Respiratory Medicine, Cochin Hospital, AP-HP, Paris, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

This study aimed to identify simple rules for allocating chronic obstructive pulmonary disease (COPD) patients to clinical phenotypes identified by cluster analyses.

Data from 2409 COPD patients of French/Belgian COPD cohorts were analysed using cluster analysis resulting in the identification of subgroups, for which clinical relevance was determined by comparing 3-year all-cause mortality. Classification and regression trees (CARTs) were used to develop an algorithm for allocating patients to these subgroups. This algorithm was tested in 3651 patients from the COPD Cohorts Collaborative International Assessment (3CIA) initiative.

Cluster analysis identified five subgroups of COPD patients with different clinical characteristics (especially regarding severity of respiratory disease and the presence of cardiovascular comorbidities and diabetes). The CART-based algorithm indicated that the variables relevant for patient grouping differed markedly between patients with isolated respiratory disease (FEV1, dyspnoea grade) and those with multi-morbidity (dyspnoea grade, age, FEV1 and body mass index). Application of this algorithm to the 3CIA cohorts confirmed that it identified subgroups of patients with different clinical characteristics, mortality rates (median, from 4% to 27%) and age at death (median, from 68 to 76 years).

A simple algorithm, integrating respiratory characteristics and comorbidities, allowed the identification of clinically relevant COPD phenotypes.

Abstract

An algorithm integrating respiratory characteristics and comorbidities identifies clinical COPD phenotypes http://ow.ly/eSRp30fJPG5

Introduction

Airflow limitation is the hallmark of chronic obstructive pulmonary disease (COPD), and forced expiratory volume in 1 s (FEV1) has long been used as the main criterion for the characterisation of disease severity [1, 2]. Analyses of observational cohorts (e.g. the ECLIPSE cohort) have revealed that COPD patients with similar levels of FEV1 experience different degrees of disease burden, reflected by dyspnoea levels, exacerbation rates, health-related quality of life (HRQoL) impairment and exercise limitation [3]. Accordingly, the current classification of COPD proposed by the Global initiative for chronic Obstructive Lung Disease (GOLD) incorporates not only the FEV1, but also dyspnoea or HRQoL, and previous occurrence of COPD exacerbations and/or hospitalisation [1]. Although this classification is not fully evidence-based, it has the advantage of taking into account some of the clinical heterogeneity of COPD with the aim of predicting future risk and proposing corresponding treatment choices. A limitation of this classification is that it does not account for age, an important determinant of prognosis in patients with COPD [4]. Furthermore, the GOLD classification does not account for comorbidities, which can both be frequent and contribute to the prognosis [5–7].

Several groups have used cluster analyses to explore clinical heterogeneity in cohorts of patients with COPD [8–10]. These studies have identified consistent clinical COPD phenotypes at high risk of mortality, including 1) younger patients with severe respiratory disease, few cardiovascular comorbidities, and poor nutritional status; and 2) older patients with moderate respiratory disease, metabolic and cardiovascular comorbidities, and obesity [11]. They have also identified patients with mild disease and a good prognosis [12, 13]. However, all published studies had limitations related to a relatively small sample size and lack of further validation in independent samples [11, 13]. Furthermore, the results of cluster analyses are difficult to translate for use in daily practice, as they provide no tool for individual patient allocation in the identified phenotypes.

In the present study, our aim was to develop and validate an algorithm, based on easily available clinical data, to assign patients with COPD to clinically relevant phenotypes.

Methods

Overall design

Data from three French/Belgian COPD cohorts were used to identify clinical COPD phenotypes using cluster analysis. Classification and Regression Tree (CART) [14] analysis was then used to develop an algorithm to allocate individual COPD patients recruited in these French/Belgian cohorts to specific subgroups. This algorithm was further tested in an independent sample of patients with COPD, using data from the COPD Cohorts Collaborative International Assessment (3CIA) initiative [15].

COPD patient cohorts

The French/Belgian COPD cohorts are composed of three cohorts: the Initiatives BPCO cohort [8], the French College of General Hospital Respiratory Physicians (CPHG) cohort [16] and the Leuven cohort [12]. Patients within these cohorts had a diagnosis of COPD, based on post-bronchodilator FEV1/FVC<0.70, and were recruited in a stable state in university hospitals (Initiatives BPCO and Leuven cohorts) [8, 12], or at the time of hospitalisation for COPD exacerbations (CPHG cohort) [16], as previously described. The 3CIA initiative contains pooled individualised data from 22 cohorts of patients with COPD, who were recruited in publicly funded hospitals or in population-based studies [15]. All cohorts were approved by a local Ethics Committee and all subjects provided informed written consent.

Statistical analysis plan

First, COPD patients recruited in the French/Belgian cohorts were classified into subgroups, based on the results of cluster analysis of data obtained at inclusion in the cohorts. The clinical relevance of the identified subgroups was established by examining their association with 3-year all-cause mortality. Next, CARTs were used for the development of an algorithm, assigning COPD patients to the subgroups identified by cluster analysis. The clinical value of this algorithm was examined using 3-year all-cause mortality in the French/Belgian cohorts. Finally, the algorithm was tested for external validation using data from the 3CIA initiative database [15]. Mortality risks among subgroups were analysed using Kaplan–Meier curves and Cox models. The concordance probability estimate was used to evaluate the discriminatory power of classifications for mortality prediction. Data are presented as median (interquartile range, IQR) or n (%). Analyses were performed using SAS 9.2 (SAS Institute Inc., Cary, NC, USA) and Tanagra 1.4 (Lyon, France) software. Additional information on the methods used can be found in the online supplementary material.

Cluster analysis of the French/Belgian COPD cohorts

Variables were selected for inclusion in the cluster analysis, based on their previous association with future risk and prognosis in COPD patients [1, 6], and included age, body mass index (BMI), FEV1 (% predicted), modified Medical Research Council (mMRC) dyspnoea scale, number of exacerbations in the previous 12 months, and presence/absence of cardiovascular comorbidities (hypertension, coronary artery disease and/or left heart failure) and/or diabetes. Identification of subgroups of patients with COPD associated with survival was achieved using factor analysis for mixed data (FAMD) [17, 18], followed by classification of patients using Ward's agglomerative hierarchical cluster analysis [8, 12]. The clinical relevance of the identified subgroups was examined by comparing their all-cause mortality at 3 years, as previously described [8, 12]. These subgroups (phenotypes) were labelled using Roman numbers.

Development of an algorithm for assigning COPD patients to specific subgroups in the French/Belgian cohorts

The development of an algorithm to assign COPD patients to the subgroups identified by cluster analysis was achieved using CART analysis [14, 19], a non-parametric decision tree learning technique [19]. Variables included in this analysis were those selected for the cluster analysis (see above). Threshold values for these variables were based on those obtained by CART analysis and were slightly modified for improved practicality (see online supplement for a detailed explanation).

External validation of the algorithm

The algorithm established in the French/Belgian cohorts was then tested in an independent group of patients with COPD from the 3CIA database. Patients in this database (n=16 332) were considered eligible for the study if data necessary to apply the algorithm (age, BMI, FEV1% predicted, mMRC scale, presence/absence of cardiovascular comorbidities and diabetes) and information on vital status at 3 years were available. Patients with appropriate data (n=3651) were classified by the algorithm into the five classes described above (labelled using Arabic digits), and these classes were compared according to their clinical characteristics, all-cause mortality at 3 years and age at death.

Results

Patients and overall study design

The study design is presented in figure 1 and the characteristics of the patients with COPD at inclusion in the French/Belgian cohorts (n=2409 patients) and in the 3CIA database (n=3651 patients) are presented in supplementary table S1. Their 3-year all-cause mortality rates were 30.8% and 11.6%, respectively. Patients included in the French/Belgian cohorts were characterised by older age, more severe airflow limitation and higher rates of cardiovascular comorbidities and/or diabetes. Furthermore, 57% of patients in the French/Belgian cohorts were recruited at the time of hospitalisation for COPD exacerbations (as part of the CPHG cohort) [16].

FIGURE 1
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 1

Study design. Patients with chronic obstructive pulmonary disease (COPD) recruited in the French/Belgian cohorts were classified into subgroups (phenotypes), based on the results of a cluster analysis of clinical data obtained at inclusion in the cohorts. Next, classification and regression trees (CARTs) were used on the same data to determine the best variables and thresholds necessary for the development of an algorithm for assigning COPD patients to the subgroups identified by cluster analysis in the French/Belgian cohorts. This analysis led to the development of a simple algorithm for allocating patients with COPD into five classes. This algorithm was then tested for external validation using data from the 3CIA initiative database (n=16 332). This latter analysis was only possible in patients with available data (n=3651), i.e. with all the variables contained in the algorithm. In each analysis, the clinical relevance of the identified subgroups/classes was established by examining their association with 3-year all-cause mortality. BMI: body mass index; FEV1: forced expiratory volume in 1 s; mMRC: modified Medical Research Council.

Cluster analysis of the French/Belgian cohorts

Table 1 shows the five subgroups (labelled I to V) identified in the French/Belgian COPD cohorts using cluster analysis (see online supplementary tables S2–S6 and figure S1). Table 2 summarises the main descriptors of these subgroups, according to increasing rates of 3-year all-cause mortality. Subgroup V (mortality rate 2.5%) was characterised by mild respiratory disease and low rates of comorbidities. Subgroup II (mortality rate 21.8%) was characterised by moderate to severe respiratory disease and low rates of comorbidities. Subgroup III (mortality rate 30.0%) was generally characterised by an older age than that of subgroup II, with a high prevalence of comorbidities and obesity. Subgroup IV (mortality rate 47.0%) was characterised by very severe respiratory disease with low rates of cardiovascular comorbidities and diabetes. Subgroup I (mortality rate 50.9%) had less severe respiratory disease than subgroup IV, but was characterised by older age and very high rates of cardiovascular comorbidities and diabetes.

View this table:
  • View inline
  • View popup
TABLE 1

Characteristics and 3-year mortality in chronic obstructive pulmonary disease (COPD) patients (n=2409) recruited in the French/Belgian COPD cohort, according to the five subgroups identified using cluster analysis

View this table:
  • View inline
  • View popup
TABLE 2

Main descriptors of the five chronic obstructive pulmonary disease (COPD) phenotypes identified by cluster analysis in the French/Belgian COPD cohort#

Use of CART for the development of an algorithm to assign COPD patients to subgroups of patients, identified by cluster analysis in the French/Belgian cohorts

The CART analysis provided an algorithm that facilitated the assignment of up to 80% of the patients to the subgroups identified by cluster analysis (see online supplementary tables S7 and S8). This algorithm is presented in figure 2 and the clinical characteristics of patients, according to the five classes obtained by applying this algorithm, are presented in table 3. Kaplan–Meier survival curves by cluster analysis-defined subgroups (figure 3a) and CART-defined classes (figure 3b) showed comparable results. Concordance probability estimates were 0.61 (95% CI 0.59–0.63) for cluster analysis-defined subgroups and 0.60 (95% CI 0.58–0.62) for CART-defined classes, confirming that both methods had comparable discriminatory power for the identification of subgroups with different prognoses.

FIGURE 2
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 2

Algorithm developed by classification and regression tree (CART) analysis for the classification of chronic obstructive pulmonary disease (COPD) patients. Application to the French/Belgian and 3CIA cohorts. BMI: body mass index; FEV1: forced expiratory volume in 1 s; mMRC: modified Medical Research Council.

View this table:
  • View inline
  • View popup
TABLE 3

Characteristics and 3-year mortality rates in chronic obstructive pulmonary disease (COPD) patients recruited in the French/Belgian COPD cohorts, or in the 3CIA initiative database according to the five classes identified using the classification and regression tree (CART)-based algorithm

FIGURE 3
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 3

Kaplan–Meier analyses for assessing all-cause mortality at 3 years. a) French/Belgian chronic obstructive pulmonary disease (COPD) cohorts according to the five subgroups (phenotypes, Ph) identified by cluster analysis. b) French/Belgian COPD cohorts according to the five classes identified by classification and regression tree (CART) analysis. c) The 3CIA COPD cohort according to the five classes identified by the algorithm developed in the French/Belgian cohorts. All analyses, p<0.0001 (Log-rank test).

Evaluation of the algorithm using data from the 3CIA initiative database

The algorithm developed in the French/Belgian cohorts was then tested, using data obtained in COPD patients from the 3CIA database. Characteristics of the 3651 patients distributed into classes, according to this algorithm, are presented in table 3. Kaplan–Meier survival curves by classes are presented in figure 3c. The concordance probability estimate was 0.62 (95% CI 0.59–0.64).

Comparison of mortality rates among classes in the French/Belgian COPD cohorts versus the 3CIA database

Because 3-year mortality rates varied widely between French/Belgian COPD cohorts and the 3CIA database, we used Cox analysis to examine hazard ratios of mortality among patients in the five classes defined by our algorithm in both cohorts, respectively. Forest plots corresponding to these analyses are presented in figure 4. Although absolute rates of death were markedly higher in the French/Belgian cohorts, hazard ratios of mortality among the five classes were rather comparable in the French/Belgian cohorts and in the 3CIA initiative.

FIGURE 4
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 4

Relative mortality risks at 3 years among chronic obstructive pulmonary disease (COPD) patients in a) the French/Belgian COPD cohorts and b) the 3CIA initiative. COPD patients were classified into five classes according to the algorithm. Horizontal bars show hazard ratios and 95% confidence intervals of mortality risks between classes. For example, in the French/Belgian COPD cohorts, subjects in class 4 have a 23.2-fold (95% CI 10.2–52.7) increased risk of mortality when compared with subjects in class 5.

Distribution by GOLD grades of severity of airflow limitation [1] in patients who died during follow-up is presented in figure 5. When comparing classes with high rates of all-cause mortality, patients without cardiovascular comorbidities/diabetes (class 4) who died were predominantly in GOLD 4; whereas patients with cardiovascular comorbidities/diabetes (class 1) who died had less severe airflow limitation (predominantly GOLD 2 and 3). Comparable observations were made when comparing patients in class 2 versus class 3 (intermediate mortality rates).

FIGURE 5
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 5

Distribution of airflow limitation severity by Global Initiative for Chronic Obstructive Lung Disease (GOLD) grade at inclusion in the cohorts, in patients who died during follow-up. a) French/Belgian cohorts. b) 3CIA initiative. Data are presented as a percentage of the total number of deaths in each class. Absolute numbers of deaths (n) in each class are also presented.

Discussion

In the present study, we first performed cluster analysis in a pool of French/Belgian COPD cohorts, which identified five subgroups (phenotypes) of patients with different rates of all-cause mortality at 3 years and different ages at death. We then used CART analysis in this pool of French/Belgian cohorts to develop an algorithm that allowed allocation of patients into five classes that corresponded to the subgroups identified by cluster analysis. This simple algorithm was based on clinical variables (including cardiovascular comorbidities and/or diabetes and respiratory characteristics) that are routinely available in daily practice. Classification of COPD patients using this algorithm allowed the identification of subgroups of patients, which differed on 3-year all-cause mortality and age at death in the pool of French/Belgian cohorts, thereby providing internal validation of the approach. This method provided comparable results in patients included in the 3CIA initiative database, which contained an independent group of patients with COPD recruited in multinational cohorts, thereby providing external validation. The algorithm identifies clinical phenotypes that are relevant to the prognosis of patients with COPD, which could aid in the exploration of underlying pathophysiological mechanisms and development of novel strategies of care.

The algorithm described in the present study is the first to integrate comorbidities (cardiovascular diseases, diabetes and obesity) and age to more classical respiratory variables (FEV1 and dyspnoea) to improve the characterisation of patients with COPD. An important yield of this algorithm is to identify patients who belong to two subgroups with a poorer prognosis, i.e. classes 1 and 4; and to highlight the corresponding determinants, i.e. the severity of the respiratory component (as assessed by the degrees of lung function impairment and dyspnoea) and the presence of major cardiovascular comorbidities or metabolic risk factors (diabetes). These data confirm previous studies, which show that 1) cardiovascular and metabolic comorbidities contribute to worsening outcomes (e.g. mortality, hospitalisation and exacerbation) in patients with COPD [6, 20]; and 2) two very different phenotypes of COPD patients with a poor prognosis exist (those with severe respiratory disease, often occurring at a younger age; and those with multi-morbidities including cardiovascular and metabolic diseases, often characterised by an older age) [9, 12]. Importantly, this study extends previous data by studying larger numbers of patients (including larger numbers of women) recruited in multiple countries, and provides a simple algorithm that can be used in the clinic to classify patients. One notable characteristic of the algorithm is that it highlights the variables on which clinicians and researchers should focus during follow-up and treatment adaptation. Whether specific strategies need to be developed for all or some of the identified phenotypes now needs to be tested prospectively. Similarly, future studies should aim to determine whether these phenotypes are associated with specific biomarkers that reflect underlying pathophysiological mechanisms.

The main strengths of the present study were the application of exploratory statistical analyses complemented by clinical knowledge in large cohorts of patients, the validation of findings in an external pool of cohorts and the use of a robust variable (mortality) for validation. We also recognise that the present study has limitations. Our assessment of comorbidities was based on physician diagnoses that did not consider occult conditions, which reportedly occur in COPD patients [21]. To limit such underestimation of the impact of undiagnosed cardiovascular diseases, the definition of cardiovascular comorbidities was relatively loose and included hypertension (a risk factor for cardiovascular disease rather than a disease itself). This definition also corresponds to what occurs in real-life daily practice, where many patients do not benefit from systematic screening for cardiovascular comorbidities. Although COPD patients are at high risk for lung cancer, which is associated with a poor prognosis, patients with active lung cancer were generally excluded from the present cohorts, thus limiting our findings to COPD patients without active lung cancer. Specific causes of mortality were not available in the cohorts used in the present analyses, and the prognostic value of the phenotypes was confirmed using all-cause mortality. Previous studies have shown that causes of mortality in COPD populations differ between patients with mild versus severe airflow obstruction, with a higher relative weight of cancer and cardiovascular causes in patients with less severe airflow impairment, and more respiratory causes in those with more severe airflow impairment [22]. Among the patients who died, differences in the GOLD grades of airflow obstruction (see figure 4) between phenotypes with comparable survival rates (e.g. class 1 versus class 4 and class 2 versus class 3) suggest that patients with relatively high rates of cardiovascular comorbidities and/or diabetes (e.g. class 1 and 3) were more likely to die from extrapulmonary causes. Importantly, even if one of its purposes is to identify populations with different mortality rates, the algorithm is not intended to represent a prognostic index, as the determinants of a given prognosis might differ markedly between patients of a given group. The large difference in mortality rates between the two groups of cohorts largely relates to the fact that 57% of patients in the French/Belgian cohorts were recruited at the time of hospitalisation for a COPD exacerbation (CPHG cohort) [16], reflecting the prognostic impact of hospitalisations. Although hospitalisation appears to be an important prognostic factor, it should be considered a marker of disease severity rather than a phenotype per se. This was the basis for not including previous hospitalisation as a variable in the cluster analysis. However, COPD exacerbations (which are important in the characterisation of patients with COPD) [23] were included in the cluster analysis and the CART analysis. The finding that exacerbations were not retained in our final algorithm should not be misinterpreted, as exacerbations remain important events in the life of patients with COPD [24]; it merely reflects that non-hospitalised exacerbations were not significantly related to prognosis. The performance of classification trees could also be improved by the integration of biomarkers that reflect inflammatory (fibrinogen, white blood cell count, C-reactive protein, eosinophils, etc.) [25–27] and cardiovascular (brain natriuretic peptide, copeptin, pro-adrenomedullin etc.) [28] biological phenomena.

The field of COPD phenotypes was once considered ‘the future of COPD’ [29], but moving from exploratory research studies to the clinic has proven to be difficult. The algorithm described in the present study offers a new way of combining and hierarchising well-known prognostic criteria (including comorbidities, age and symptoms) to identify COPD phenotypes in the clinic. This approach could serve as a basis to develop phenotype-specific therapeutic strategies, by recruiting appropriate at-risk target populations in clinical trials. We speculate that our algorithm might also help in unravelling specific biological pathways that were previously missed, owing to the mixing of various phenotypes in the current classifications of COPD.

Supplementary material

Supplementary Material

Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.

Online supplement ERJ-01034-2017_supplement

Disclosures

Supplementary Material

I. Alfageme ERJ-01034-2017_Alfageme

P. Bakke ERJ-01034-2017_Bakke

P.-R. Burgel ERJ-01034-2017_Burgel

C. Casanova ERJ-01034-2017_Casanova

B. Celli ERJ-01034-2017_Celli

BG. Cosio ERJ-01034-2017_Cosio

M. Decramer ERJ-01034-2017_Decramer

A.L. Echazarreta ERJ-01034-2017_Echazarreta

M. Han ERJ-01034-2017_Han

W. Janssens ERJ-01034-2017_Janssens

P. Lange ERJ-01034-2017_Lange

J.M. Marin ERJ-01034-2017_Marin

M. Miravitlles ERJ-01034-2017_Mirevitlles

T. Oga ERJ-01034-2017_Oga

J.-L. Paillasseur ERJ-01034-2017_Paillasseur

A.S. Ramírez García Luna ERJ-01034-2017_Ramirez

N. Roche ERJ-01034-2017_Roche

D. Sin ERJ-01034-2017_Sin

J.J. Soler-Cataluña ERJ-01034-2017_SolerCataluña

J.B. Soriano ERJ-01034-2017_Soriano

A. Turner ERJ-01034-2017_Turner

Footnotes

  • This article has supplementary material available from erj.ersjournals.com

  • Support statement: The analyses reported here were supported by an unrestricted grant from Boehringer Ingelheim France, which played no role in study design, data collection, analysis and interpretation of data, writing of the manuscript nor decision to submit it for publication. The Initiatives BPCO study was supported by an unrestricted grant from Boehringer Ingelheim France and (until 2015) Pfizer. None of the funding sources of the individual trials were involved in any aspect of the 3CIA initiative, including the design, data collection and analysis, decision to publish, or preparation of the manuscript. P. Martinez-Camblor was supported by research grant MTM2011-23204 from the Spanish Ministerio de Ciencia e Innovación (FEDER support included). J. Garcia-Aymerich has a researcher contract from the Instituto de Salud Carlos III (CP05/00118), Ministry of Health, Spain. Funding information for this article has been deposited with the Crossref Funder Registry.

  • Conflict of interest: Disclosures can be found alongside this article at erj.ersjournals.com

  • Received May 19, 2017.
  • Accepted July 28, 2017.
  • Copyright ©ERS 2017

References

  1. ↵
    1. Vogelmeier CF,
    2. Criner GJ,
    3. Martinez FJ
    , et al. Global Strategy for the Diagnosis, Management, and Prevention of Chronic Obstructive Lung Disease 2017 Report: GOLD Executive Summary. Eur Respir J 2017; 49: 1700214.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Celli BR,
    2. Decramer M,
    3. Wedzicha JA
    , et al. An official American Thoracic Society/European Respiratory Society statement: research questions in COPD. Eur Respir J 2015; 45: 879–905.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Agusti A,
    2. Calverley PM,
    3. Celli B
    , et al. Characterisation of COPD heterogeneity in the ECLIPSE cohort. Respir Res 2010; 11: 122.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Puhan MA,
    2. Garcia-Aymerich J,
    3. Frey M
    , et al. Expansion of the prognostic assessment of patients with chronic obstructive pulmonary disease: the updated BODE index and the ADO index. Lancet 2009; 374: 704–711.
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. Divo M,
    2. Cote C,
    3. de Torres JP
    , et al. Comorbidities and risk of mortality in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med 2012; 186: 155–161.
    OpenUrlCrossRefPubMedWeb of Science
  6. ↵
    1. Mannino DM,
    2. Thorn D,
    3. Swensen A
    , et al. Prevalence and outcomes of diabetes, hypertension and cardiovascular disease in COPD. Eur Respir J 2008; 32: 962–969.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Fabbri LM,
    2. Boyd C,
    3. Boschetto P
    , et al. How to integrate multiple comorbidities in guideline development: article 10 in integrating and coordinating efforts in COPD guideline development. An official ATS/ERS workshop report. Proc Am Thorac Soc 2012; 9: 274–281.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Burgel PR,
    2. Paillasseur JL,
    3. Caillaud D
    , et al. Clinical COPD phenotypes: a novel approach using principal component and cluster analyses. Eur Respir J 2010; 36: 531–539.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Garcia-Aymerich J,
    2. Gomez FP,
    3. Benet M
    , et al. Identification and prospective validation of clinically relevant chronic obstructive pulmonary disease (COPD) subtypes. Thorax 2011; 66: 430–437.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Rennard S,
    2. Locantore N,
    3. Delafont B
    , et al. 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.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Pinto LM,
    2. Alghamdi M,
    3. Benedetti A
    , et al. Derivation and validation of clinical phenotypes for COPD: a systematic review. Respir Res 2015; 16: 50.
    OpenUrl
  12. ↵
    1. Burgel PR,
    2. Paillasseur JL,
    3. Peene B
    , et al. Two distinct chronic obstructive pulmonary disease (COPD) phenotypes are associated with high risk of mortality. Plos ONE 2012; 7: e51048.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Burgel PR,
    2. Paillasseur JL,
    3. Roche N
    . Identification of clinical phenotypes using cluster analyses in COPD patients with multiple comorbidities. Biomed Res Int 2014; 2014: 420134.
    OpenUrl
  14. ↵
    1. Breiman L,
    2. Friedman J,
    3. Olshen R
    , et al. Classification and regression trees. Monterey, CA, Wadsworth & Brooks, 1984.
  15. ↵
    1. Soriano JB,
    2. Lamprecht B,
    3. Ramirez AS
    , et al. Mortality prediction in chronic obstructive pulmonary disease comparing the GOLD 2007 and 2011 staging systems: a pooled analysis of individual patient data. Lancet Respir Med 2015; 3: 443–450.
    OpenUrl
  16. ↵
    1. Piquet J,
    2. Chavaillon JM,
    3. David P
    , et al. High-risk patients following hospitalisation for an acute exacerbation of COPD. Eur Respir J 2013; 42: 946–955.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Wikipedia
    . Factor analysis of mixed data. https://en.wikipedia.org/wiki/Factor_analysis_of_mixed_data Date last accessed: January 25, 2016. Date last updated: May 17, 2016.
  18. ↵
    1. Pagès J
    . Analyse factorielle de données mixtes. [Multiple Factor Analysis for Mixed Data]. Rev Statistique Appliquée 2004; 52: 93–111.
    OpenUrl
  19. ↵
    1. Wikipedia
    . Predictive analytics. https://en.wikipedia.org/wiki/Predictive_analytics#Classification_and_regression_trees Date last accessed: January 25, 2016. Date last updated: October 03, 2017.
  20. ↵
    1. Briggs A,
    2. Spencer M,
    3. Wang H
    , et al. Development and validation of a prognostic index for health outcomes in chronic obstructive pulmonary disease. Arch Intern Med 2008; 168: 71–79.
    OpenUrlCrossRefPubMedWeb of Science
  21. ↵
    1. Rutten FH,
    2. Cramer MJM,
    3. Grobbee DE
    , et al. Unrecognized heart failure in elderly patients with stable chronic obstructive pulmonary disease. Eur Heart J 2005; 26: 1887–1894.
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    1. Sin DD,
    2. Anthonisen NR,
    3. Soriano JB
    , et al. Mortality in COPD: role of comorbidities. Eur Respir J 2006; 28: 1245–1257.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Hurst JR,
    2. Vestbo J,
    3. Anzueto A
    , et al. Susceptibility to exacerbation in chronic obstructive pulmonary disease. N Engl J Med 2010; 363: 1128–1138.
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    1. Wedzicha JA,
    2. Seemungal TA
    . COPD exacerbations: defining their cause and prevention. Lancet 2007; 370: 786–796.
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    1. Duvoix A,
    2. Dickens J,
    3. Haq I
    , et al. Blood fibrinogen as a biomarker of chronic obstructive pulmonary disease. Thorax 2013; 68: 670–676.
    OpenUrlAbstract/FREE Full Text
    1. Thomsen M,
    2. Dahl M,
    3. Lange P
    , et al. Inflammatory biomarkers and comorbidities in chronic obstructive pulmonary disease. Am J Respir Crit Care Med 2012; 186: 982–988.
    OpenUrlCrossRefPubMedWeb of Science
  26. ↵
    1. George L,
    2. Brightling CE
    . Eosinophilic airway inflammation: role in asthma and chronic obstructive pulmonary disease. Ther Adv Chronic Dis 2016; 7: 34–51.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Stolz D,
    2. Meyer A,
    3. Rakic J
    , et al. Mortality risk prediction in COPD by a prognostic marker panel. Eur Respir J 2014; 44: 1557–1570.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Han MK,
    2. Agusti A,
    3. Calverley PM
    , et al. Chronic obstructive pulmonary disease phenotypes: the future of COPD. Am J Respir Crit Care Med 2010; 182: 598–604.
    OpenUrlCrossRefPubMedWeb of Science
View Abstract
PreviousNext
Back to top
View this article with LENS
Vol 50 Issue 5 Table of Contents
European Respiratory Journal: 50 (5)
  • Table of Contents
  • Index by author
Email

Thank you for your interest in spreading the word on European Respiratory Society .

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
A simple algorithm for the identification of clinical COPD phenotypes
(Your Name) has sent you a message from European Respiratory Society
(Your Name) thought you would like to see the European Respiratory Society web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
Citation Tools
A simple algorithm for the identification of clinical COPD phenotypes
Pierre-Régis Burgel, Jean-Louis Paillasseur, Wim Janssens, Jacques Piquet, Gerben ter Riet, Judith Garcia-Aymerich, Borja Cosio, Per Bakke, Milo A. Puhan, Arnulf Langhammer, Inmaculada Alfageme, Pere Almagro, Julio Ancochea, Bartolome R. Celli, Ciro Casanova, Juan P. de-Torres, Marc Decramer, Andrés Echazarreta, Cristobal Esteban, Rosa Mar Gomez Punter, MeiLan K. Han, Ane Johannessen, Bernhard Kaiser, Bernd Lamprecht, Peter Lange, Linda Leivseth, Jose M. Marin, Francis Martin, Pablo Martinez-Camblor, Marc Miravitlles, Toru Oga, Ana Sofia Ramírez, Don D. Sin, Patricia Sobradillo, Juan J. Soler-Cataluña, Alice M. Turner, Francisco Javier Verdu Rivera, Joan B. Soriano, Nicolas Roche
European Respiratory Journal Nov 2017, 50 (5) 1701034; DOI: 10.1183/13993003.01034-2017

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Share
A simple algorithm for the identification of clinical COPD phenotypes
Pierre-Régis Burgel, Jean-Louis Paillasseur, Wim Janssens, Jacques Piquet, Gerben ter Riet, Judith Garcia-Aymerich, Borja Cosio, Per Bakke, Milo A. Puhan, Arnulf Langhammer, Inmaculada Alfageme, Pere Almagro, Julio Ancochea, Bartolome R. Celli, Ciro Casanova, Juan P. de-Torres, Marc Decramer, Andrés Echazarreta, Cristobal Esteban, Rosa Mar Gomez Punter, MeiLan K. Han, Ane Johannessen, Bernhard Kaiser, Bernd Lamprecht, Peter Lange, Linda Leivseth, Jose M. Marin, Francis Martin, Pablo Martinez-Camblor, Marc Miravitlles, Toru Oga, Ana Sofia Ramírez, Don D. Sin, Patricia Sobradillo, Juan J. Soler-Cataluña, Alice M. Turner, Francisco Javier Verdu Rivera, Joan B. Soriano, Nicolas Roche
European Respiratory Journal Nov 2017, 50 (5) 1701034; DOI: 10.1183/13993003.01034-2017
del.icio.us logo Digg logo Reddit logo Technorati logo Twitter logo CiteULike logo Connotea logo Facebook logo Google logo Mendeley logo
Full Text (PDF)

Jump To

  • Article
    • Abstract
    • Abstract
    • Introduction
    • Methods
    • Results
    • Discussion
    • Supplementary material
    • Disclosures
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

Subjects

  • COPD and smoking
  • Tweet Widget
  • Facebook Like
  • Google Plus One

More in this TOC Section

Original articles

  • Identifying early PAH biomarkers in systemic sclerosis
  • Viable virus aerosol propagation by PAP circuit leak
  • Ambulatory management of secondary spontaneous pneumothorax
Show more Original articles

COPD

  • Airway smooth muscle area to predict steroid responsive in COPD
  • Lung volume reduction surgery versus endobronchial valves
  • Long-acting bronchodilator combination therapy and cardiovascular events in COPD
Show more COPD

Related Articles

Navigate

  • Home
  • Current issue
  • Archive

About the ERJ

  • Journal information
  • Editorial board
  • Press
  • Permissions and reprints
  • Advertising

The European Respiratory Society

  • Society home
  • myERS
  • Privacy policy
  • Accessibility

ERS publications

  • European Respiratory Journal
  • ERJ Open Research
  • European Respiratory Review
  • Breathe
  • ERS books online
  • ERS Bookshop

Help

  • Feedback

For authors

  • Instructions for authors
  • Publication ethics and malpractice
  • Submit a manuscript

For readers

  • Alerts
  • Subjects
  • Podcasts
  • RSS

Subscriptions

  • Accessing the ERS publications

Contact us

European Respiratory Society
442 Glossop Road
Sheffield S10 2PX
United Kingdom
Tel: +44 114 2672860
Email: journals@ersnet.org

ISSN

Print ISSN:  0903-1936
Online ISSN: 1399-3003

Copyright © 2023 by the European Respiratory Society