Next Article in Journal
Selection Criteria for Determination of Optimal Reconstruction Method for Cu-64 Trastuzumab Dosimetry on Siemens Inveon PET Scanner
Previous Article in Journal
Intra-Vitreal Administration of Microvesicles Derived from Human Adipose-Derived Multipotent Stromal Cells Improves Retinal Functionality in Dogs with Retinal Degeneration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Disease-Specific Comorbidity Clusters in COPD and Accelerated Aging

1
CIRO+, Centre of Expertise for Chronic Organ Failure, 6085 NM Horn, The Netherlands
2
Department of Respiratory Medicine, MUMC+, Maastricht University Medical Centre, 6229 ER Maastricht, The Netherlands
3
Department of Respiratory Medicine, AZ Sint-Lucas, 9000 Ghent, Belgium
4
NUTRIM School of Nutrition and Translational Research in Metabolism, Maastricht University Medical Centre+, 6229 ER Maastricht, The Netherlands
5
Viscovery Software GmbH, 1130 Vienna, Austria
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2019, 8(4), 511; https://doi.org/10.3390/jcm8040511
Submission received: 17 March 2019 / Revised: 31 March 2019 / Accepted: 11 April 2019 / Published: 14 April 2019
(This article belongs to the Section Pulmonology)

Abstract

:
Background: Patients with chronic obstructive pulmonary disease (COPD) often suffer from multiple morbidities, which occur in clusters and are sometimes related to accelerated aging. This study aimed to assess the disease specificity of comorbidity clusters in COPD and their association with a biomarker of accelerated aging as a potential mechanistic factor. Methods: Body composition, metabolic, cardiovascular, musculoskeletal, and psychological morbidities were objectively evaluated in 208 COPD patients (age 62 ± 7 years, 58% males, FEV1 50 ± 16% predicted) and 200 non-COPD controls (age 61 ± 7 years, 45% males). Based on their presence and severity, the morbidities were clustered to generate distinct clusters in COPD and controls. Telomere length in circulating leukocytes was compared across the clusters. Results: (co)morbidities were more prevalent in COPD patients compared to controls (3.9 ± 1.7 vs. 2.4 ± 1.5, p < 0.05). A “Psychologic” and “Cachectic” cluster were only present in the COPD population. “Less (co)morbidity”, “Cardiovascular”, and “Metabolic” clusters were also observed in controls, although with less complexity. Telomere length was reduced in COPD patients, but did not differ between the (co)morbidity clusters in both populations. Conclusions: Two COPD-specific comorbidity clusters, a “Cachectic” and “Psychologic” cluster, were identified and warrant further studies regarding their development. Accelerated aging was present across various multimorbidity clusters in COPD.

1. Introduction

Chronic obstructive pulmonary disease (COPD) is considered a complex, heterogeneous, and multicomponent condition [1]. Comorbidities are prevalent in COPD, impact on morbidity and mortality [2], and appear in various patterns [3]. Using a systematic and objective assessment, Vanfleteren et al. [4] identified five comorbidity clusters (a “Less Comorbidity”, “Cardiovascular”, “Metabolic”, “Psychologic”, and “Cachectic” cluster). Their identification suggested that differential pathophysiological pathways underlie these clusters. At a given genetic background, chronic diseases may develop and progress (at different speeds in various combinations) in response to common risk factors, such as smoking, alcohol, diet, pollution, and physical inactivity [5]. These environmental factors and intrinsic mechanisms may lead to accelerated aging [6]. Reduced telomere length, a surrogate marker for accelerated aging, was reported in circulating leukocytes in COPD [7], myocardial infarction, stroke, and diabetes [8]. While, some resemblance in comorbidity associations between COPD patients and (unmatched) controls was shown [9], the co-existence of multiple comorbidities in COPD patients and aging controls has not been explored by a unifying cluster analysis, nor by objective measurements. We hypothesized that comorbidity clustering in patients with COPD might differ from that of a non-COPD elderly population, and might relate to accelerated aging. We aimed to compare and validate the clustering of (co)morbidities in a sample of COPD and non-COPD elderly control subjects, and secondarily, to compare telomere length between the different (co)morbidity clusters in both COPD and non-COPD subjects.

2. Experimental Section

2.1. Study Design and Population

The current cross-sectional analysis is part of the “Individualized COPD Evaluation in relation to Ageing” (ICE Age) study [10] conducted between December 2010 and August 2016 in a tertiary care pulmonary rehabilitation center in the southeastern region of the Netherlands. Inclusion and exclusion criteria have been published before (Appendix A) [10]. In summary, stable COPD patients with moderate to severe airflow limitation and elderly controls without COPD or a history of debilitating chronic disease were included. All participants provided written informed consent prior to study participation. The study was conducted in accordance with the Declaration of Helsinki and good clinical practice guidelines. It was approved by the local ethics review board and registered on www.isrctn.com.

2.2. Assessments

Details of patient characterization, (co)morbidity assessment, and references of the (co)morbidity definitions are provided in Appendix A. Body composition, metabolic, cardiovascular, musculoskeletal, and psychological morbidities were objectively evaluated by the following measurements [4]: body mass index (BMI); fat-free mass index (FFMI by dual X-ray absorptiometry (DEXA); bone mineral density (BMD by DEXA) at the hip, lumbar spine, and whole body; pulse wave velocity (PWV); carotid intima-media thickness (c-IMT); peripheral systolic and diastolic blood pressure; fasting plasma glucose and serum insulin concentration to calculate the homeostasis model assessment method (HOMA) index; plasma creatinine concentration to calculate the glomerular filtration rate (GFR, using the simplified Modification-of-Diet-in-Renal-Disease equation); plasma triglyceride and high-density lipoprotein (HDL) levels; and Hospital Anxiety and Depression Scale (HADS). Telomere length in circulating leukocytes was determined by a monochrome multiplex quantitative PCR-procedure.

2.3. Definitions

Thirteen (co)morbidities were identified based on predefined cut-offs, mostly corresponding to our previous cluster analysis [4]. In summary: obesity was defined as BMI ≥ 30 kg/m2; hyperglycemia as plasma glucose concentration ≥ 5.6 mmol/L; insulin resistance as HOMA index ≥ 4.29/4.43 for women/men; dyslipidemia as plasma triglyceride concentration > 1.7 mmol/L or HDL < 1.29/1.03 mmol/L for women/men; hypertension as systolic blood pressure ≥ 140 mmHg or diastolic blood pressure ≥ 90 mmHg; arterial stiffness as PWV > 10 m/s; atherosclerosis as c-IMT > 0.9 mm; underweight as BMI < 21 kg/m2; low muscle mass as FFMI < 14.62/17.05 kg/m2 for women/men; osteoporosis as a T score < −2.5/−2.8 in women/men; renal impairment as estimated GFR < 60 mL/min/1.73 m2); and anxiety or depression as a HADS score ≥ 10 points.

2.4. Statistical Analysis

Statistical analyses were performed using Viscovery SOMine 7.1 by Viscovery Software GmbH (www.viscovery.net; Vienna, Austria). Self-organizing maps (SOMs) were used to create an ordered representation for the presence of COPD and (co)morbidities. The SOM method can be viewed as a non-parametric regression technique that converts multidimensional data spaces into lower dimensional abstractions. It generates a non-linear representation of the data distribution and allows the user to identify homogenous data groups visually. Patients were primarily ordered by the presence/absence of COPD and secondarily by their overall similarity concerning their present (co)morbidities and also by the degree of its presence given by parameters from which the (co)morbidities were calculated. Based on the created SOM model, clusters were generated using the SOM–Ward Cluster algorithm that applies the classical hierarchical method of Ward on top of the SOM topology. Summary variables are presented as mean ± standard deviation for quantitative variables (except pack-years: median (interquartile range)), and percentage for discrete variables. Comorbidities, clinical characteristics, and telomere length were compared between: (1) The COPD and control group; (2) each (co)morbidity cluster and its corresponding group (COPD or control group); (3) each parallel (co)morbidity cluster in COPD and controls, using the integrated two-sided t-test with a confidence of 95%. Post-hoc in SPSS 19.0, additional binary logistic regression models were run within each group (unadjusted; adjusted for age; adjusted for age and sex), using the cluster as a dependent variable and telomere length as an independent variable.

3. Results

3.1. General Characteristics and (Co)Morbidities

Two-hundred-and-eight patients with COPD (62 ± 7 years) and 200 non-COPD controls (61 ± 7 years) were included (Figure 1). Patients with COPD were more often male and had a more extensive smoking history (Table 1). On average, both groups were slightly overweight, but controls had a lower cardiovascular risk. The COPD subjects had moderate to severe airflow limitation, moderately impaired diffusion capacity, increased static lung volumes, and impaired health status. Except for hypertension, atherosclerosis, and dyslipidemia, most (co)morbidities were more prevalent in patients with COPD (Figure 2a, baseline and missing data: Appendix A, Table A1). On average, patients with COPD had more (3.9 ± 1.7 versus 2.4 ± 1.5, p < 0.05, Figure 2b) and more overlapping (Figure 3) (co)morbidities than controls.

3.2. (Co)Morbidity Clusters

Five comorbidity clusters were identified in the COPD patients, and three (co)morbidity clusters in the non-COPD group (Figure 4).

3.2.1. Comorbidity Clusters in COPD

Cluster C1 (“Less Comorbidity” cluster in COPD; n = 46 (22%)) had significantly fewer morbidities compared to the whole COPD group (Table 2). Cluster C2 (“Cardiovascular”, n = 27 (13%)) had the highest prevalence of hypertension, arterial stiffness, and osteoporosis, but a lower prevalence of metabolic features. Cluster C3 (“Metabolic”, n = 56 (27%)) had a higher prevalence of obesity, hyperglycemia, insulin resistance, and dyslipidemia, but also cardiovascular features. Cluster C4 (“Psychologic”, n = 40 (19%)) had the highest number of subjects with anxiety and depression. Cluster C5 (“Cachectic”, n = 39 (19%)) had the highest portion of subjects with underweight and low muscle mass.
Functional exercise capacity, dyspnea scores, and frequency of exacerbations in the year prior to the study were not different among the comorbidity clusters (Table 3). However, patients in cluster C1 (“Less Comorbidity”) had less pack-years, less cardiovascular risk, and a better health-related quality of life. Cluster C2 (“Cardiovascular”) included older patients who were less frequently active smokers. Patients in cluster C3 (“Metabolic”) had a lower degree of airflow limitation, and static hyperinflation, and better diffusion capacity, while Framingham-defined cardiovascular risk was increased. Cluster 4 (“Psychologic”) had a higher proportion of women, a higher amount of pack-years, and a worse health status. Patients in cluster C5 (“Cachectic”) were more frequently active smokers, had worse pulmonary function with a higher degree of airflow limitation and static hyperinflation, but a lower Framingham-defined cardiovascular risk.

3.2.2. (Co)Morbidity Clusters in Controls

Cluster X1 (“Less (co)morbidity” cluster in controls) was a large cluster (n = 97 of 200 non-COPD controls (49%)) and was characterized by a clearly fewer number of (co)morbidities (Table 4). In cluster X2 (“Cardiovascular”, n = 34 (17%)), hypertension and arterial stiffness was present in all subjects. Cluster X3 (“Metabolic”, n = 69 (35%)) had the highest portion of subjects with obesity, hyperglycemia, insulin resistance, and dyslipidemia.
Cluster X1 (“Less (co)morbidity”) was a younger, predominantly female cluster with less pack-years, and a lower cardiovascular risk (Table 5). Cluster X2 (“Cardiovascular”) was older with a higher Framingham-defined cardiovascular risk. Cluster X3 (“Metabolic”) also had a high cardiovascular risk and accumulated more pack-years.

3.2.3. Comparing (Co)Morbidity Clusters in COPD and Controls

Details of the comparison of the parallel (co)morbidity clusters “Less (Co)morbidity”, “Cardiovascular”, and “Metabolic” are shown in Table 4 and Table 5 (footnote 2 of the table). In summary, COPD subjects in the “Less (Co)morbidity” cluster (C1) had more (co)morbidities than those in the parallel control cluster (X1), including hyperglycemia, dyslipidemia, obesity, arterial stiffness, and renal impairment. The presence of (co)morbidities in the “Cardiovascular” clusters (C2, X2) was mainly comparable, except for a higher frequency of osteoporosis and less dyslipidemia in the COPD cluster C2. Both cluster C3 and X3 had a high prevalence of metabolic diseases. The “Metabolic” COPD cluster (C3) was older, more overweight, and had more (co)morbidities than cluster X3, in particular insulin resistance, hypertension, and arterial stiffness.

3.3. Telomere Length

Telomere length was reduced in the patients with COPD compared to controls (Table 1). The mean telomere length of the (co)morbidity clusters did not differ significantly between the COPD clusters, nor between the non-COPD controls clusters (Figure 5 and Appendix A, Table A2), with similar results after adjusting for age and gender (Appendix A, Table A3 and Table A4). As for the parallel (co)morbidity clusters in COPD and controls, telomere length was comparable between both “Less (Co)morbidity” clusters C1 and X1, and significantly reduced in cluster C2 versus X2, while the difference failed to reach statistical significance for cluster C3 versus X3 (10.22 ± 1.90 versus 10.82 ± 1.45 kbp, p = 0.05) (Figure 5 and Appendix A, Table A2).

4. Discussion

This study showed that a “Psychologic” and “Cachectic” comorbidity cluster were specifically related to COPD, and that “Cardiovascular”, “Metabolic”, and “Less-(co)morbidity” clusters were also present in aging controls, while it validates the previously identified comorbidity clusters [4]. Accelerated aging, measured by reduced telomere length, was present across various multimorbidity clusters in COPD, but was not limited to COPD-specific clusters.
Using a similar statistical approach, but a slightly different assessment of comorbidities, the current study identified a predominantly similar comorbidity clustering, compared to that in the CIROCO cohort [4], except for a few intriguing differences. Osteoporosis was now most prevalent in the “Cardiovascular” instead of the “Cachectic” cluster. This might be driven by the inclusion of arterial stiffness in the clustering model, which in turn relates to osteoporosis [11]. With respect to the clinical characterization, we confirmed the worse quality of life in the “Psychologic” cluster, the higher amount of smokers in the “Cachectic”, and former smokers in an older “Cardiovascular” cluster. In both studies, females tended to be more represented in the “Psychologic” or “Cachectic” cluster.
A systematic review of studies that derived COPD phenotypes using unsupervised methods identified two recognizable phenotypes with poor longitudinal health outcomes across multiple studies [3]. A phenotype with moderate respiratory disease, obesity, cardiovascular, and metabolic comorbidities corresponds to the “Metabolic” cluster. Conversely, a phenotype with poor nutritional and health status, severe respiratory disease, and few cardiovascular comorbidities corresponds to the “Cachectic” cluster [3]. Indeed, underweight is prevalent in COPD, often co-occurs with low muscle mass [4], and is related to a more severe disease status [12]. Low body weight itself is also an independent risk factor for developing COPD [13].
The formation of emphysematous lesions might be promoted by cellular senescence which infers with tissue regeneration [14]. In turn, alveolar type II and endothelial cells from emphysematous lungs show shortened telomeres associated with increased cell senescence [15]. In a study in mice, it was reported that short telomeres lower the thresh-hold of cigarette smoke-induced damage, which implicated telomere length as a genetic susceptibility factor in emphysema [16]. On the other hand, emphysema severity is also associated with arterial stiffness [17] and shorter telomere length predict arterial stiffness [18]. On the other site of the phenotypic spectrum, telomere attrition has also been linked to type 2 diabetes [8]. These associations between telomere length and different comorbidities and emphysema, a widespread hallmark lesion in COPD [19], might explain the comparable presence of accelerated aging in the various comorbidity clusters in COPD. Also, in a subgroup of the BODE cohort, no differences in pulmonary function, nor BMI were found across telomere tertiles [20].
Anxiety and depression often co-occur [21,22] and are prevalent in COPD, especially in females [23]. Symptoms, poor health-related quality of life, behavioral factors, and social isolation in COPD contribute to the presence of mood disorders [24]. While the current study confirms the presence of an independent (predominantly female) “Psychologic” cluster [4], anxiety and depression are also associated to each other in non-COPD subjects [9]. The lower prevalence of mood disorders in the control population should be considered as a cause for the absence of a “Psychologic” cluster in the control population.
Uniquely, this cluster study included elderly controls, allowing for a meaningful comparison of their—rigorously and objectively assessed—comorbidity profiles with those of patients with COPD, although the control population was, on average, one year younger. Most of them suffered from several morbidities, which underlines the importance of including a control population when studying multimorbidity from the perspective of an index disease [25]. Whereas three corresponding (co)morbidity profiles were identified in controls, these clusters were more complex in COPD. First, a higher number of comorbidities was seen in the “Less Comorbidity” COPD cluster, emphasizing the clinical need to consider comorbidities in all patients with COPD. Second, osteoporosis was more prevalent in the “Cardiovascular” COPD cluster, which could be linked to the use of oral corticosteroids in patients with COPD. Third, subjects in the “Metabolic” COPD cluster were older, more obese, and had more insulin resistance and cardiovascular comorbidity, which warrants future pathophysiologic research on cardiometabolic disease in COPD patients.
Considerations: Most COPD patients were recruited in a tertiary care pulmonary rehabilitation setting, which limits the external validity of the current findings. Reflecting clinical reality, patients with COPD had a more substantial smoking history than their controls. While influenced by smoking, reduced telomere length in COPD compared to smoking and non-smoking controls was consistently reported [7,20,26,27]. Telomere attrition is only one hallmark in the complex concept of aging [28], and the study was not powered to exclude differences in telomere attrition between the clusters.

5. Conclusions

This study confirms previously identified clusters of objectively identified multiple comorbidities in a well-characterized cohort of patients with COPD and adds to this the comparison to a non-COPD elderly population, and the evaluation of telomere length between clusters, as a measure of biological aging. The current findings underline the complexity and heterogeneity of patients with COPD. Interestingly, exacerbation frequency, dyspnea, and functional exercise performance were similar amongst COPD comorbidity clusters, emphasizing again that detailed phenotyping is needed beyond the aforementioned clinical outcomes. In addition, a subgroup with low nutritional status and more severe lung function impairment, and a well-represented subgroup with psychological conditions are more specifically related to the presence of COPD and are important subgroups for targeted mechanistic and intervention studies. Reduced telomere length, is present across various comorbidity clusters in COPD, but did not seem a major determinant of (co)morbidity clustering itself.

Author Contributions

Conceptualization and methodology, all authors; formal analysis and investigation, F.J.J.T., N.R., S.G., and L.E.G.W.V.; resources and funding acquisition, E.P.A.R., E.F.M.W.; data curation F.J.J.T., S.G., and E.P.A.R.; writing—original draft preparation, all authors; writing—review and editing, F.J.J.T., F.M.E.F., M.A.S., and L.E.G.W.V.; visualization, F.J.J.T. and S.G.; supervision, F.M.E.F., E.F.M.W., and L.E.G.W.V.; project administration, F.J.J.T. and E.P.A.R.

Funding

This research was partly funded by the Lung Foundation Netherlands, previously Dutch Asthma Foundation, grant number 3.2.09.049 and an unrestricted grant from the Weijerhorst Foundation.

Acknowledgments

We thank the patients with COPD and the non-COPD controls who volunteered to participate in this study; Trineke Hofstra; Ans Suntjens, Marco Akkermans, Linda Op ‘t Veld, Koen Stakenborg, Jos Peeters, Martijn Cuijpers; Annie van de Kruijs, RN; Irma Timmermans, RN, Miriam Groenen, and Riny van Kessel, RN for planning and performing all the tests at CIRO+; Renske Krijgsman and Claudia Geijselaers for technical assistance with laboratory analyses at the Department of Respiratory Medicine; and J.M. Dekker for providing data of the HOMA index in the Hoorn Study.

Conflicts of Interest

F.J.J.T., F.M.E.F., N.R., S.G., E.P.A.R., and L.E.G.W.V. declare no conflict of interest. M.A.S. reports personal fees from Boehringer Ingelheim AstraZeneca and GSK outside the submitted work. D.J.A.J. reports personal fees from Almirall, Chiesi, GSK, Boehringer Ingelheim, and Mayne Pharma outside the submitted work. E.F.M.W. reports personal fees from Board membership (Nycomed, Boehringer), grants (AstraZeneca, GSK), personal fees from payment for lectures including service on speakers’ bureaus (AstraZeneca, GSK, Novartis, Chiesi) outside the submitted work. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results

Appendix A

Appendix A.1. Subject Enrollment

Inclusion and exclusion criteria have been published before [10]. COPD patients as well as non-COPD controls, 45–75 years of age were included. Following exclusion criteria were defined: any kind of carcinogenic pathology <5 years before study participation, other chronic lung diseases or previous lung surgery, alpha-1-antitrypsin deficiency, chronic use of oral corticosteroids >10 mg prednisolone/day, uncertainty about the willingness or ability of the subject to comply with the protocol requirements, and participation in a study involving investigational or marketed products concomitantly or <4 weeks prior to study entry. Inclusion criteria for patients with COPD were: moderate to severe airflow limitation (post-bronchodilator forced expiratory volume of 1s (FEV1) < 80% predicted, FEV1/forced vital capacity (FVC) < 70%) [2], and clinical stability, defined as absence of respiratory tract infection or exacerbation for <4 weeks before study entry. Control subjects should not have COPD or other diagnosis of debilitating chronic disease. Patients with COPD were recruited on admission to a tertiary care pulmonary rehabilitation center or from the same southeastern region of the Netherlands. Smoking and never-smoking (<1 pack-year) non-COPD controls were also recruited from this region.

Appendix A.2. Assessments

Appendix A.2.1. General Characteristics

At study entry, subjects’ demographics, smoking status, and Framingham 10-year risk for cardiovascular disease [29] were documented. In subjects with COPD, the number of exacerbations in the past year, dyspnea (Medical Research Council (MRC) grade) [30]), six-minute walk distance [31], and health related quality of life (St. George’s Respiratory Questionnaire (SGRQ) [32]) were also assessed. Baseline post-bronchodilator spirometry using a standardized spirometer, measurement of static lung volumes, and carbon monoxide transfer factor lung function (Masterlab®, Jaeger, Würzburg, Germany) was performed in all subjects.

Appendix A.2.2. Assessments of (Co)Morbidities

The following clinically relevant variables were assessed:
  • Body mass index (BMI = body weight in kg/(height in m)2. Height was measured to the nearest 0.1 cm. Body weight was assessed to the nearest 0.1 kg after emptying the bladder and with the subjects standing barefoot and wearing light indoor clothing.
  • Fat-free mass index (FFMI) = fat free mass in kg (=lean mass + bone mineral content by dual X-ray absorptiometry (DEXA) scan (Lunar Prodigy, GE Healthcare)) divided by (height in m)2 [33].
  • Bone mineral density (BMD by DEXA scan) at the hip, lumbar spine, and whole body [34].
  • Anxiety and depression were assessed by using the Hospital Anxiety and Depression Scale (HADS), a validated and reliable measurement instrument used widely in medically ill patients to screen for clinically relevant symptoms of anxiety and/or depression [35]. The HADS is a self-administered questionnaire, consisting of 14 questions. The HADS is divided into an anxiety subscale (HADS-A) and a depression subscale (HADS-D), both containing seven items. Total scores for each subscale can range from 0 (optimal) to 21 (worst) points.
  • Pulse wave velocity (PWV) assessment methods were similar in the CIROCO and ICE Age studies and have been published before [36]. Radial artery waveforms were recorded with a high-fidelity micromanometer (Millar Instruments, Houston, TX, USA). The APWV was measured by recording ECG-gated carotid and femoral artery waveforms. The Sphygmocor systems software helped to assure the quality of the pulse wave measurement. A detailed screen showed 10 s of recorded and analyzed waveforms which can be examined to assess overall consistency of the waveforms. In addition, a detailed report helped to interpret the consistency of the waveforms during the 10-s measurement. Only when the pulse height and diastolic variation was less than 5%, the average pulse height was more than 100 units, the augmentation index was less than 50%, T1 was between 80 and 150 ms, and the overall “quality index” above 80 points, data were retained. After marking the exact location measurements, it was repeated three times for securing reproducibility. A measurement was accepted when it was reproducible three times with minimal variation as judged by the biomedical technologist. The retained APWV was the mean of the three measurements. Shortest distances from manubrium to the marked location on the femoral artery (via the navel) were measured. Wave transit time was calculated by the system software, using the R-wave of the simultaneously recording electrocardiography as reference frame. The APWV was determined by dividing the distance between the two recording sites by the wave transit time.
  • Carotid intima-media thickness (c-IMT) was assessed using high-resolution B-mode ultrasound with a 10-MHz linear transducer (Art. LabEsoatePicus, Pie-medical Netherlands/Italy). The ultrasound device was connected to an acquisition modem with an automatic boundary detection system (Art.LabExoatePicus). Carotid IMT was thus quantified semi-automatically, reducing the interobserver variability [37]. With the patient in a supine position, measurements of the carotid IMT were performed throughout 10-mm segments across the bifurcation free of visual plaques. The probe was moved to obtain measurements of the left common carotid artery at 4 angles (180°, 150°, 120°, and 90°). For each measured segment, mean and maximum IMT values were acquired automatically throughout the 10-mm vessel length. The average of segmental maximum carotid values was determined as carotid IMT per patient.
  • Peripheral systolic and diastolic blood pressure measurements were performed early in the morning after 15 min in a resting supine position. A small head pillow was accepted. Peripheral blood pressure was measured three times with intervals of 5 min.
  • Laboratory analysis: A venous blood sample was collected in the fasted state. Serum, plasma, and peripheral blood mononuclear samples were stored at −80 °C prior to their analysis. Glucose, creatinine, high-density lipoprotein (HDL), and triglycerides were determined in all subjects. Fasting plasma glucose and serum insulin concentration were used to calculate the homeostasis model assessment method (HOMA index: fasting serum insulin (μU/mL) × fasting plasma glucose (mmol/L)/22.5), plasma creatinine concentration to calculate the estimated glomerular filtration rate (eGFR) using the simplified Modification of Diet in Renal Disease (MDRD) equation [38].

Appendix A.2.3. Telomere Length

Telomere length of blood leucocytes was determined by a monochrome multiplex quantitative PCR procedure as described by Cawthon [39] with DNA of two reference cell lines with known telomere length [40] included in each run. Primers used for the telomeres were: Telg (ACACTAAGGTTTGGGTTTGGGTTTGGGTTTGGGTTAGTGT), and Telc (TGTTAGGTATCCCTATCCCTATCCCTATCCCTATCCCTAACA). For the single copy reference gene beta-globin the primers used were: hbgu (CGGCGGCGGGCGGCGCGGGCTGGGCGGcttcatccacgttcaccttg), and hbgd (GCCCGGCCCGCCGCGCCCGTCCCGCCGgaggagaagtctgccgtt).

Appendix A.3. Definitions

  • COPD and grades of severity were classified according to the GOLD document [41].
  • OBESITY was defined as a BMI equal to or above 30 kg/m2 according to the WHO international classification of BMI [42].
  • HYPERGLYCEMIA was defined as fasting glucose level equal to or above 5.6 mmol/L according to the American Diabetes Association [43].
  • INSULIN RESISTANCE. The homeostasis model assessment method (HOMA 1 IR index) is a simple and validated method to assess insulin sensitivity [44]. It is also an independent predictor of cardiovascular disease in patients with type 2 diabetes mellitus [45]. It was calculated as follows: HOMA index = fasting serum insulin (μU/mL) × fasting plasma glucose (mmol/L)/22.5. Cut-off values depended on geographic variation and gender, and defining a local population, percentile-based cut-off is recommended for the HOMA index and/or hyperinsulinemia [46,47]. Therefore, we used the 75th percentile of the HOMA 1 IR of the Hoorn study, a large Dutch cohort study in men and women aged 50 to 75 years was used to define insulin resistance [48]. We based the cut-off on the entire population of the Hoorn study without exclusion of patients with diabetes mellitus and/or prevalent cardiovascular disease and used separate values for males and females. These data were provided by personal communication with JM Dekker. Insulin resistance was defined as a HOMA 1 IR equal to or above 4.43 for males and 4.29 for females.
  • DYSLIPIDEMIA was defined as a triglyceride level above 1.7 mmol/L or a HDL cholesterol level below 1.03 mmol/L in males or below 1.29 mmol/L in females [49].
  • HYPERTENSION was defined as a systolic blood pressure equal to or above 140 mmHg or diastolic pressure equal to or above 90 mmHg, in line with the European Society of Hypertension (ESH)/European Society of Cardiology (ESC) guidelines [50].
  • ARTERIAL STIFFNESS was defined as pulse wave velocity (PWV) > 10 m/s, as proposed by the expert consensus document of the Artery Society, the European Society of Hypertension Working Group on Vascular Structure and Function, and the European Network for Noninvasive Investigation of Large Arteries [51]. This value is at the upper part of the second quartile in the Framingham Heart Study and represents in this general population with a mean age of 63 years about 4% risk for a first major cardiovascular event within the next 8 years [52].
  • ATHEROSCLEROSIS was defined as a maximal carotid intima media thickness (IMT) > 0.9 mm, which is according to the ESH/ESC guidelines compatible with subclinical organ damage [50].
  • UNDERWEIGHT was defined as a BMI lower than 21 kg/m2. The purpose of a BMI cut-off point is to identify within a population, the proportion of people with a high risk of an undesirable health state. Hence, a BMI lower than 21 kg/m2, which has shown to be independently associated with one-year mortality in COPD, is an appropriate cut-off [42,53].
  • LOW MUSCLE MASS was defined as a fat free mass index (FFMI) less than 17.05 kg/m2 for men or less than 14.62 kg/m2 for women, according to the lowest 10th percentile of the general population as reported in the Copenhagen City Heart Study [54], which is also supported by the Nutritional assessment and therapy in COPD statement of the European Respiratory Society statement [33].
  • OSTEOPOROSIS was assessed by DEXA scanning and established if the lowest T score of the lumbar spine or hip was lower than −2.5 for females and lower than −2.8 for males, according to Dutch guidelines [55].
  • RENAL IMPAIRMENT was defined as an eGFR less than 60 mL/min/1.73 m2, which was estimated, using the simplified Modification of Diet in Renal Disease (MDRD) equation [38], corresponding with stage 3 chronic kidney disease according to the National Kidney Foundation Kidney Disease Outcome Quality Initiative (NKF KDOQI) guidelines [56].
  • ANXIETY was assessed, using the Hospital Anxiety and Depression Scale (HADS). A score equal to or greater than 10 points on the items related to anxiety, was considered to indicate symptoms of anxiety [35].
  • DEPRESSION was also assessed, using the HADS. A score equal to or greater than 10 points on the items related to depression, was considered to indicate symptoms of depression [35].

Appendix A.4. Baseline Data

Table A1. Comparison of baseline data between COPD and controls.
Table A1. Comparison of baseline data between COPD and controls.
COPDAll Controls
Plasma glucose level, mmol/L6.2 ± 1.65.9 ± 0.9 1
Systolic blood pressure, mmHg145 ± 22143 ± 20
Diastolic blood pressure, mmHg84 ± 984 ± 10
Pulse wave velocity, m/s10.6 ± 2.78.8 ± 2.2 1
Plasma triglyceride level, mmol/L1.5 ± 0.81.4 ± 1.0
Plasma HDL level, mmol/L1.6 ± 0.51.6 ± 0.5
cIMT, µm836 ± 181800 ± 156 1
HOMA index3.4 ± 6.02.0 ± 1.4 1
eGFR, mL/min/1.73 m273.6 ± 15.877.6 ± 13.6 1
BMI, kg/m226.9 ± 5.426.9 ± 3.4
FFMI, kg/m217.6 ± 2.618.1 ± 3.1
Lowest T score−1.7 ± 1.0−0.9 ± 1.1 1
HADS anxiety, score6.6 ± 4.23.5 ± 2.7 1
HADS depression, score6.0 ± 3.92.0 ± 2.2 1
Abbreviations: HDL = high density lipoprotein; cIMT= carotid intima media thickness; HOMA = homeostasis model assessment method; GFR = glomerular filtration rate; BMI = body mass index; FFMI = fat-free mass index; HADS = Hospital Anxiety and Depression Scale. Summary variables are presented as mean ± standard deviation for quantitative variables, and percentage for discrete variables. 1 p < 0.05 versus COPD group. Missing data for subjects with COPD: 26 subjects with COPD lacked an acceptable quality pulse wave velocity measurement; 15 subjects lacked an acceptable quality carotid intima-media thickness measurement, 11 did not fill out the entire hospital anxiety and depression scale, 1 subject did not give a blood sample, for 6 additional patients insulin levels could not be measured, and 4 subjects did not have a dual-energy X-ray absorptiometry. Missing data for controls: 9 control subjects lacked an acceptable quality pulse wave velocity measurement; 8 subjects lacked an acceptable quality carotid intima-media thickness measurement; insulin level could not be measured for one subject; 1 did not fill out the entire hospital anxiety and depression scale.

Appendix A.5. Telomere Length

Table A2. Telomere length across the (co)morbidity clusters in COPD and controls.
Table A2. Telomere length across the (co)morbidity clusters in COPD and controls.
C1C2C3C4C5X1X2X3
Subjects, n (%)46
(22%)
27
(13%)
56
(27%)
40
(19%)
39
(19%)
97
(49%)
34
(17%)
69
(35%)
Telomere length, kbp10.71 ± 1.9910.27 ± 1.9010.22 ± 1.9010.20 ± 1.3110.44 ± 1.7310.92 ± 1.5411.23 ± 1.32 110.82 ± 1.45 2
Abbreviations: C = COPD cluster. X = control cluster. kbp = kilobase pairs. Summary variables are presented as mean ± standard deviation for quantitative variables and percentage for discrete variables. The mean telomere length of the (co)morbidity clusters did not differ significantly between the COPD clusters, nor between the non-COPD controls clusters. 1 p < 0.05 versus the parallel COPD cluster. 2 p = 0.05 versus the parallel COPD cluster.
Table A3. Telomere length and morbidity clusters in COPD: binary logistic regression.
Table A3. Telomere length and morbidity clusters in COPD: binary logistic regression.
Crude ORAge Adj. ORAge–Sex Adj. OR
C11.15 (0.95–1.40)1.15 (0.94–1.39)1.15 (0.95–1.40)
C20.96 (0..77–1.21)0.99 (0.78–1.25)0.99 (0.78–1.26)
C30.94 (0.79–1.12)0.95 (0.80–1.13)0.95 (0.80–1.14)
C40.94 (0.77–1.14)0.93 (0.76–1.12)0.91 (0.75–1.11)
C51.03 (0.84–1.25)1.01 (0.83–1.23)1.01 (0.82–1.23)
Odds ratios (ORs) from binary logistic regression models run within the COPD group, using the cluster as a dependent variable and telomere length as an independent variable: 1. unadjusted; 2. adjusted for age; 3. adjusted for age and sex. Odds ratios (95% confidence interval) are reported.
Table A4. Telomere length and morbidity clusters in controls: binary logistic regression.
Table A4. Telomere length and morbidity clusters in controls: binary logistic regression.
Crude ORAge Adj. ORAge–Sex Adj. OR
X10.98 (0.81–1.19)0.99 (0.82–1.20)1.01 (0.83–1.23)
X21.19 (0.91–1.54)1.22 (0.91–1.63)1.21 (0.91–1.61)
X30.92 (0.75–1.12)0.92 (0.76–1.13)0.91 (0.74–1.11)
Odds ratios (ORs) from binary logistic regression models performed within the control group, using the cluster as a dependent variable and telomere length as an independent variable: 1. unadjusted; 2. adjusted for age; 3. adjusted for age and sex. Odds ratios (95% confidence interval) are reported.

References

  1. Agusti, A.G. COPD, a multicomponent disease: Implications for management. Respir. Med. 2005, 99, 670–682. [Google Scholar] [CrossRef]
  2. Vogelmeier, C.F.; Criner, G.J.; Martinez, F.J.; Anzueto, A.; Barnes, P.J.; Bourbeau, J.; Celli, B.R.; Chen, R.; Decramer, M.; Fabbri, L.M.; et al. Global Strategy for the Diagnosis, Management, and Prevention of Chronic Obstructive Lung Disease 2017 Report: GOLD Executive Summary. Eur. Respir. J. 2017. [Google Scholar] [CrossRef]
  3. Pinto, L.M.; Alghamdi, M.; Benedetti, A.; Zaihra, T.; Landry, T.; Bourbeau, J. Derivation and validation of clinical phenotypes for COPD: A systematic review. Respir. Res. 2015, 16, 50. [Google Scholar] [CrossRef] [PubMed]
  4. Vanfleteren, L.E.; Spruit, M.A.; Groenen, M.; Gaffron, S.; van Empel, V.P.; Bruijnzeel, P.L.; Rutten, E.P.; Op ‘t Roodt, J.; Wouters, E.F.; Franssen, F.M. 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. [Google Scholar] [CrossRef] [PubMed]
  5. Clini, E.M.; Beghé, B.; Fabbri, L.M. Chronic Obstructive Pulmonary Disease Is Just One Component of the Complex Multimorbidities in Patients with COPD. Am. J. Respir. Crit. Care Med. 2013, 187, 668–671. [Google Scholar] [CrossRef] [PubMed]
  6. Ito, K.; Barnes, P.J. COPD as a disease of accelerated lung aging. Chest 2009, 135, 173–180. [Google Scholar] [CrossRef] [PubMed]
  7. Savale, L.; Chaouat, A.; Bastuji-Garin, S.; Marcos, E.; Boyer, L.; Maitre, B.; Sarni, M.; Housset, B.; Weitzenblum, E.; Matrat, M.; et al. Shortened telomeres in circulating leukocytes of patients with chronic obstructive pulmonary disease. Am. J. Respir. Crit. Care Med. 2009, 179, 566–571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. D’Mello, M.J.; Ross, S.A.; Briel, M.; Anand, S.S.; Gerstein, H.; Pare, G. Association between shortened leukocyte telomere length and cardiometabolic outcomes: Systematic review and meta-analysis. Circ. Cardiovasc. Genet. 2015, 8, 82–90. [Google Scholar] [CrossRef] [PubMed]
  9. Divo, M.J.; Casanova, C.; Marin, J.M.; Pinto-Plata, V.M.; de-Torres, J.P.; Zulueta, J.J.; Cabrera, C.; Zagaceta, J.; Sanchez-Salcedo, P.; Berto, J.; et al. COPD comorbidities network. Eur. Respir. J. 2015, 46, 640–650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Rutten, E.P.A.; Gopal, P.; Wouters, E.F.M.; Franssen, F.M.E.; Hageman, G.J.; Vanfleteren, L.E.; Spruit, M.A.; Reynaert, N.L. VArious mechanistic pathways representing the aging process are altered in copd. Chest 2016, 149, 53–61. [Google Scholar] [CrossRef] [PubMed]
  11. Sabit, R.; Bolton, C.E.; Edwards, P.H.; Pettit, R.J.; Evans, W.D.; McEniery, C.M.; Wilkinson, I.B.; Cockcroft, J.R.; Shale, D.J. Arterial stiffness and osteoporosis in chronic obstructive pulmonary disease. Am. J. Respir. Crit. Care Med. 2007, 175, 1259–1265. [Google Scholar] [CrossRef] [PubMed]
  12. Divo, M.J.; Cabrera, C.; Casanova, C.; Marin, J.M.; Pinto-Plata, V.; de-Torres, J.P.; Zulueta, J.; Zagaceta, J.; Sanchez-Salcedo, P.; Berto, J. Comorbidity Distribution, Clinical Expression and Survival in COPD Patients with Different Body Mass Index. Chronic Obstr. Pulm. Dis. 2014, 1, 229–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Harik-Khan, R.I.; Fleg, J.L.; Wise, R.A. Body mass index and the risk of COPD. Chest 2002, 121, 370–376. [Google Scholar] [CrossRef]
  14. Aoshiba, K.; Nagai, A. Senescence Hypothesis for the Pathogenetic Mechanism of Chronic Obstructive Pulmonary Disease. Proc. Am. Thorac. Soc. 2009, 6, 596–601. [Google Scholar] [CrossRef]
  15. Tsuji, T.; Aoshiba, K.; Nagai, A. Alveolar cell senescence in patients with pulmonary emphysema. Am. J. Respir. Crit. Care Med. 2006, 174, 886–893. [Google Scholar] [CrossRef] [PubMed]
  16. Alder, J.K.; Guo, N.; Kembou, F.; Parry, E.M.; Anderson, C.J.; Gorgy, A.I.; Walsh, M.F.; Sussan, T.; Biswal, S.; Mitzner, W.; et al. Telomere length is a determinant of emphysema susceptibility. Am. J. Respir. Crit. Care Med. 2011, 184, 904–912. [Google Scholar] [CrossRef] [PubMed]
  17. McAllister, D.A.; Maclay, J.D.; Mills, N.L.; Mair, G.; Miller, J.; Anderson, D.; Newby, D.E.; Murchison, J.T.; Macnee, W. Arterial Stiffness is Independently Associated with Emphysema Severity in Patients with COPD. Am. J. Respir. Crit. Care Med. 2007, 176, 1208–1214. [Google Scholar] [CrossRef]
  18. Benetos, A.; Okuda, K.; Lajemi, M.; Kimura, M.; Thomas, F.; Skurnick, J.; Labat, C.; Bean, K.; Aviv, A. Telomere length as an indicator of biological aging: The gender effect and relation with pulse pressure and pulse wave velocity. Hypertension 2001, 37, 381–385. [Google Scholar] [CrossRef]
  19. Coxson, H.O.; Dirksen, A.; Edwards, L.D.; Yates, J.C.; Agusti, A.; Bakke, P.; Calverley, P.M.; Celli, B.; Crim, C.; Duvoix, A.; et al. The presence and progression of emphysema in COPD as determined by CT scanning and biomarker expression: A prospective analysis from the ECLIPSE study. Lancet Respir. Med. 2013, 1, 129–136. [Google Scholar] [CrossRef]
  20. Cordoba-Lanus, E.; Cazorla-Rivero, S.; Espinoza-Jimenez, A.; de-Torres, J.P.; Pajares, M.J.; Aguirre-Jaime, A.; Celli, B.; Casanova, C. Telomere shortening and accelerated aging in COPD: Findings from the BODE cohort. Respir. Res. 2017, 18, 59. [Google Scholar] [CrossRef]
  21. Baty, F.; Putora, P.M.; Isenring, B.; Blum, T.; Brutsche, M. Comorbidities and burden of COPD: A population based case-control study. PLoS ONE 2013, 8, e63285. [Google Scholar] [CrossRef] [PubMed]
  22. Chubachi, S.; Sato, M.; Kameyama, N.; Tsutsumi, A.; Sasaki, M.; Tateno, H.; Nakamura, H.; Asano, K.; Betsuyaku, T. Identification of five clusters of comorbidities in a longitudinal Japanese chronic obstructive pulmonary disease cohort. Respir. Med. 2016, 117, 272–279. [Google Scholar] [CrossRef] [PubMed]
  23. Maurer, J.; Rebbapragada, V.; Borson, S.; Goldstein, R.; Kunik, M.E.; Yohannes, A.M.; Hanania, N.A.; for the ACCP Workshop Panel on Anxiety and Depression in COPD. Anxiety and depression in COPD: Current understanding, unanswered questions, and research needs. Chest 2008, 134, 43S–56S. [Google Scholar] [CrossRef]
  24. Yohannes, A.M.; Alexopoulos, G.S. Depression and anxiety in patients with COPD. Eur. Respir. Rev. 2014, 23, 345–349. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Barnett, K.; Mercer, S.W.; Norbury, M.; Watt, G.; Wyke, S.; Guthrie, B. Epidemiology of multimorbidity and implications for health care, research, and medical education: A cross-sectional study. Lancet 2012, 380, 37–43. [Google Scholar] [CrossRef]
  26. Albrecht, E.; Sillanpaa, E.; Karrasch, S.; Alves, A.C.; Codd, V.; Hovatta, I.; Buxton, J.L.; Nelson, C.P.; Broer, L.; Hagg, S.; et al. Telomere length in circulating leukocytes is associated with lung function and disease. Eur. Respir. J. 2014, 43, 983–992. [Google Scholar] [CrossRef] [PubMed]
  27. Boyer, L.; Chouaïd, C.; Bastuji-Garin, S.; Marcos, E.; Margarit, L.; Le Corvoisier, P.; Vervoitte, L.; Hamidou, L.; Frih, L.; Audureau, E.; et al. Aging-Related Systemic Manifestations in COPD Patients and Cigarette Smokers. PLoS ONE 2015, 10, e0121539. [Google Scholar] [CrossRef] [PubMed]
  28. Lopez-Otin, C.; Blasco, M.A.; Partridge, L.; Serrano, M.; Kroemer, G. The hallmarks of aging. Cell 2013, 153, 1194–1217. [Google Scholar] [CrossRef]
  29. D’Agostino, R.B., Sr.; Vasan, R.S.; Pencina, M.J.; Wolf, P.A.; Cobain, M.; Massaro, J.M.; Kannel, W.B. General cardiovascular risk profile for use in primary care: The Framingham Heart Study. Circulation 2008, 117, 743–753. [Google Scholar] [CrossRef]
  30. Fletcher, C.M.; Elmes, P.C.; Fairbairn, A.S.; Wood, C.H. The significance of respiratory symptoms and the diagnosis of chronic bronchitis in a working population. Br. Med. J. 1959, 2, 257–266. [Google Scholar] [CrossRef]
  31. Hernandes, N.A.; Wouters, E.F.; Meijer, K.; Annegarn, J.; Pitta, F.; Spruit, M.A. Reproducibility of 6-minute walking test in patients with COPD. Eur. Respir. J. 2011, 38, 261–267. [Google Scholar] [CrossRef] [PubMed]
  32. Jones, P.W.; Quirk, F.H.; Baveystock, C.M.; Littlejohns, P. A self-complete measure of health status for chronic airflow limitation. The St. George’s Respiratory Questionnaire. Am. Rev. Respir. Dis. 1992, 145, 1321–1327. [Google Scholar] [CrossRef] [PubMed]
  33. Schols, A.M.; Ferreira, I.M.; Franssen, F.M.; Gosker, H.R.; Janssens, W.; Muscaritoli, M.; Pison, C.; Rutten-van Molken, M.; Slinde, F.; Steiner, M.C.; et al. Nutritional assessment and therapy in COPD: A European Respiratory Society statement. Eur. Respir. J. 2014, 44, 1504–1520. [Google Scholar] [CrossRef] [PubMed]
  34. Graat-Verboom, L.; Spruit, M.A.; van den Borne, B.E.; Smeenk, F.W.; Wouters, E.F. Whole-Body versus Local DXA-Scan for the Diagnosis of Osteoporosis in COPD Patients. J. Osteoporos. 2010, 2010, 640878. [Google Scholar] [CrossRef] [PubMed]
  35. Zigmond, A.S.; Snaith, R.P. The hospital anxiety and depression scale. Acta Psychiatr. Scand. 1983, 67, 361–370. [Google Scholar] [CrossRef] [PubMed]
  36. Vanfleteren, L.E.; Spruit, M.A.; Groenen, M.T.; Bruijnzeel, P.L.; Taib, Z.; Rutten, E.P.; Op ‘t Roodt, J.; Akkermans, M.A.; Wouters, E.F.; Franssen, F.M. Arterial stiffness in patients with COPD: The role of systemic inflammation and the effects of pulmonary rehabilitation. Eur. Respir. J. 2014, 43, 1306–1315. [Google Scholar] [CrossRef] [PubMed]
  37. Rossi, A.C.; Brands, P.J.; Hoeks, A.P. Automatic recognition of the common carotid artery in longitudinal ultrasound B-mode scans. Med. Image Anal. 2008, 12, 653–665. [Google Scholar] [CrossRef]
  38. Levey, A.S.; Greene, T.; Kusek, J.W.; Beck, G.J. A simplified equation to predict glomerular filtration rate from serum creatinine. J. Am. Soc. Nephrol. 2000, 11, 155A. [Google Scholar]
  39. Cawthon, R.M. Telomere length measurement by a novel monochrome multiplex quantitative PCR method. Nucleic Acids Res. 2009, 37, e21. [Google Scholar] [CrossRef]
  40. de Vos-Houben, J.M.; Ottenheim, N.R.; Kafatos, A.; Buijsse, B.; Hageman, G.J.; Kromhout, D.; Giltay, E.J. Telomere length, oxidative stress, and antioxidant status in elderly men in Zutphen and Crete. Mech. Ageing Dev. 2012, 133, 373–377. [Google Scholar] [CrossRef]
  41. Vestbo, J.; Hurd, S.S.; Agusti, A.G.; Jones, P.W.; Vogelmeier, C.; Anzueto, A.; Barnes, P.J.; Fabbri, L.M.; Martinez, F.J.; Nishimura, M.; 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. [Google Scholar] [CrossRef] [PubMed]
  42. Consultation, W.E. Appropriate body-mass index for Asian populations and its implications for policy and intervention strategies. Lancet 2004, 363, 157–163. [Google Scholar] [CrossRef]
  43. American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care 2010, 33, S62–S69. [Google Scholar] [CrossRef] [PubMed]
  44. Bonora, E.; Targher, G.; Alberiche, M.; Bonadonna, R.C.; Saggiani, F.; Zenere, M.B.; Monauni, T.; Muggeo, M. Homeostasis model assessment closely mirrors the glucose clamp technique in the assessment of insulin sensitivity: Studies in subjects with various degrees of glucose tolerance and insulin sensitivity. Diabetes Care 2000, 23, 57–63. [Google Scholar] [CrossRef]
  45. Bonora, E.; Formentini, G.; Calcaterra, F.; Lombardi, S.; Marini, F.; Zenari, L.; Saggiani, F.; Poli, M.; Perbellini, S.; Raffaelli, A.; et al. HOMA-Estimated Insulin Resistance Is an Independent Predictor of Cardiovascular Disease in Type 2 Diabetic Subjects: Prospective data from the Verona Diabetes Complications Study. Diabetes Care 2002, 25, 1135–1141. [Google Scholar] [CrossRef] [Green Version]
  46. Gayoso-Diz, P.; Otero-Gonzalez, A.; Rodriguez-Alvarez, M.X.; Gude, F.; Garcia, F.; De Francisco, A.; Quintela, A.G. Insulin resistance (HOMA-IR) cut-off values and the metabolic syndrome in a general adult population: Effect of gender and age: EPIRCE cross-sectional study. BMC Endocr. Disord. 2013, 13, 47. [Google Scholar] [CrossRef]
  47. Balkau, B.; Charles, M.A.; Drivsholm, T.; Borch-Johnsen, K.; Wareham, N.; Yudkin, J.S.; Morris, R.; Zavaroni, I.; van Dam, R.; Feskins, E.; et al. Frequency of the WHO metabolic syndrome in European cohorts, and an alternative definition of an insulin resistance syndrome. Diabetes Metab. 2002, 28, 364–376. [Google Scholar]
  48. Dekker, J.M.; Girman, C.; Rhodes, T.; Nijpels, G.; Stehouwer, C.D.; Bouter, L.M.; Heine, R.J. Metabolic syndrome and 10-year cardiovascular disease risk in the Hoorn Study. Circulation 2005, 112, 666–673. [Google Scholar] [CrossRef]
  49. Alberti, K.G.; Zimmet, P.; Shaw, J. The metabolic syndrome—A new worldwide definition. Lancet 2005, 366, 1059–1062. [Google Scholar] [CrossRef]
  50. Mancia, G.; De Backer, G.; Dominiczak, A.; Cifkova, R.; Fagard, R.; Germano, G.; Grassi, G.; Heagerty, A.M.; Kjeldsen, S.E.; Laurent, S.; et al. 2007 Guidelines for the Management of Arterial Hypertension: The Task Force for the Management of Arterial Hypertension of the European Society of Hypertension (ESH) and of the European Society of Cardiology (ESC). J. Hypertens. 2007, 25, 1105–1187. [Google Scholar] [CrossRef] [PubMed]
  51. Van Bortel, L.M.; Laurent, S.; Boutouyrie, P.; Chowienczyk, P.; Cruickshank, J.K.; De Backer, T.; Filipovsky, J.; Huybrechts, S.; Mattace-Raso, F.U.; Protogerou, A.D.; et al. Expert consensus document on the measurement of aortic stiffness in daily practice using carotid-femoral pulse wave velocity. J. Hypertens. 2012, 30, 445–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Mitchell, G.F.; Hwang, S.J.; Vasan, R.S.; Larson, M.G.; Pencina, M.J.; Hamburg, N.M.; Vita, J.A.; Levy, D.; Benjamin, E.J. Arterial stiffness and cardiovascular events: The Framingham Heart Study. Circulation 2010, 121, 505–511. [Google Scholar] [CrossRef] [Green Version]
  53. Celli, B.R.; Cote, C.G.; Marin, J.M.; Casanova, C.; Montes de Oca, M.; Mendez, R.A.; Pinto Plata, V.; Cabral, H.J. The Body-Mass Index, Airflow Obstruction, Dyspnea, and Exercise Capacity Index in Chronic Obstructive Pulmonary Disease. N. Engl. J. Med. 2004, 350, 1005–1012. [Google Scholar] [CrossRef] [PubMed]
  54. Vestbo, J.; Prescott, E.; Almdal, T.; Dahl, M.; Nordestgaard, B.G.; Andersen, T.; Sorensen, T.I.; Lange, P. Body mass, fat-free body mass, and prognosis in patients with chronic obstructive pulmonary disease from a random population sample: Findings from the Copenhagen City Heart Study. Am. J. Respir. Crit. Care Med. 2006, 173, 79–83. [Google Scholar] [CrossRef] [PubMed]
  55. Reumatologie, N.V.v. Richtlijn osteoporose en fractuurpreventie, 3e herziening. In Hoofdstuk 8: Diagnostiek en Behandeling van Osteoporose bij Mannen; NVR: Utrecht, The Netherlands, 2011. [Google Scholar]
  56. Kdoqi. KDOQI Clinical Practice Guidelines and Clinical Practice Recommendations for Diabetes and Chronic Kidney Disease. Am. J. Kidney Dis. 2007, 49, S12–S154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Subject enrollment. COPD = chronic obstructive pulmonary disease.
Figure 1. Subject enrollment. COPD = chronic obstructive pulmonary disease.
Jcm 08 00511 g001
Figure 2. Frequencies of (co)morbidities. (a) Frequencies of (co)morbidities in COPD and controls: * More prevalent in the COPD group (p < 0.05); (b) Number of (co)morbidities per subject: mean number of (co)morbidities: 3.9 ± 1.7 in COPD vs. 2.4 ± 1.5 in controls (p < 0.05).
Figure 2. Frequencies of (co)morbidities. (a) Frequencies of (co)morbidities in COPD and controls: * More prevalent in the COPD group (p < 0.05); (b) Number of (co)morbidities per subject: mean number of (co)morbidities: 3.9 ± 1.7 in COPD vs. 2.4 ± 1.5 in controls (p < 0.05).
Jcm 08 00511 g002
Figure 3. Co-existence of (co)morbidities in COPD and controls. The frequencies of co-existing (co)morbidities (%, ranging from 0 to 100%, horizontal row) in the presence of every index morbidity (vertical column) in COPD (left) and controls (right) are shown. These are descriptive data and do not offer a statistical comparison. For interpretation, the cells are color filled: red indicating the highest values, blue the lowest, and yellow for those in between.
Figure 3. Co-existence of (co)morbidities in COPD and controls. The frequencies of co-existing (co)morbidities (%, ranging from 0 to 100%, horizontal row) in the presence of every index morbidity (vertical column) in COPD (left) and controls (right) are shown. These are descriptive data and do not offer a statistical comparison. For interpretation, the cells are color filled: red indicating the highest values, blue the lowest, and yellow for those in between.
Jcm 08 00511 g003
Figure 4. (Co)morbidity clusters in COPD and controls. Panels generated using Viscovery (Viscovery Software GmbH, Vienna, Austria). The Viscovery program primarily ordered all subjects on the map based on the presence or absence of COPD (second panel from the left in upper row) and secondarily by their overall similarity concerning the prevalent (co)morbidities and the degree (severity) of these (co)morbidities. The more subjects are comparable in terms of their (co)morbidity profile, the closer they are on the map. Contrarily, the more they differ, the further they are away from each other. Within the COPD patients, 5 clusters were identified (C1–C5), while there were only 3 in the controls (X1–X3). Within the panels showing the specific comorbidities, a red color indicates the presence of a specific comorbidity, while a blue color represents the absence of this condition. All panels represent the same map and a subject is always on the same place on the map. If the same node on the map colors red for “Low Muscle Mass” and “Underweight”, but blue for “Insulin Resistance”, a subject in that node, has a low muscle mass and is underweight, but does not have insulin resistance.
Figure 4. (Co)morbidity clusters in COPD and controls. Panels generated using Viscovery (Viscovery Software GmbH, Vienna, Austria). The Viscovery program primarily ordered all subjects on the map based on the presence or absence of COPD (second panel from the left in upper row) and secondarily by their overall similarity concerning the prevalent (co)morbidities and the degree (severity) of these (co)morbidities. The more subjects are comparable in terms of their (co)morbidity profile, the closer they are on the map. Contrarily, the more they differ, the further they are away from each other. Within the COPD patients, 5 clusters were identified (C1–C5), while there were only 3 in the controls (X1–X3). Within the panels showing the specific comorbidities, a red color indicates the presence of a specific comorbidity, while a blue color represents the absence of this condition. All panels represent the same map and a subject is always on the same place on the map. If the same node on the map colors red for “Low Muscle Mass” and “Underweight”, but blue for “Insulin Resistance”, a subject in that node, has a low muscle mass and is underweight, but does not have insulin resistance.
Jcm 08 00511 g004
Figure 5. Telomere length across (co)morbidity clusters in COPD and controls. Scatter plots of telomere length across (co)morbidity clusters in COPD and controls, representing mean ± standard deviation. C1–5: COPD. X1–3: controls. n.s. = not significant. * Indicates significantly (p < 0.05) reduced telomere length in the COPD cluster compared to its parallel control cluster. For cluster C3 vs. X3: p = 0.05.
Figure 5. Telomere length across (co)morbidity clusters in COPD and controls. Scatter plots of telomere length across (co)morbidity clusters in COPD and controls, representing mean ± standard deviation. C1–5: COPD. X1–3: controls. n.s. = not significant. * Indicates significantly (p < 0.05) reduced telomere length in the COPD cluster compared to its parallel control cluster. For cluster C3 vs. X3: p = 0.05.
Jcm 08 00511 g005
Table 1. Subject characteristics.
Table 1. Subject characteristics.
COPDControls
Number of subjects, n208200
Male, n (%)5845 1
Age, years62 ± 761 ± 7 1
BMI, kg/m226.9 ± 5.426.9 ± 3.4
Former smoker, %6555 1
Current smoker, %3411 1
Pack-years, n43 (29–58)6 (0–20) 1
Fram.10-yr. risk ≥ 30%, %3219 1
FEV1, % predicted50 ± 16119 ± 15 1
FEV1/FVC41 ± 1279 ± 5 1
ITGV, % predicted144 ± 33100 ± 18 1
KCO, % predicted69 ± 22101 ± 16 1
MRC dyspnea grade3.0 ± 1.0NA
SGRQ, total score55.0 ± 16.5NA
6MWD, m458 ± 131NA
≥2 ex. previous year, n (%)45NA
Telomere length, kbp10.37 ± 1.7810.94 ± 1.47 1
Abbreviations: BMI = body mass index; Fram.10-yr. risk = Framingham 10-year risk of overall cardiovascular disease; FEV1 = Forced expiratory volume in 1 s; FVC = Forced vital capacity. ITGV = Intrathoracic gas volume. KCO = Transfer factor for carbon monoxide per alveolar volume. MRC = Medical Research Council. SGRQ = St. George’s Respiratory Questionnaire. 6MWD = 6-min walking distance. ex. = Exacerbations. kbp = Kilobase pairs. NA = Not applicable. Summary variables are presented as mean ± standard deviation for quantitative variables, with exception of pack-years (median and interquartile range), and percentage for discrete variables. 1 p < 0.05 versus COPD group.
Table 2. Prevalence of comorbidities across the comorbidity clusters in COPD.
Table 2. Prevalence of comorbidities across the comorbidity clusters in COPD.
C1C2C3C4C5
Cluster descriptionLess ComorbidityCardiovascularMetabolicPsychologicCachectic
Subjects, n (%)46 (22%)27 (13%)56 (27%)40 (19%)39 (19%)
Comorbidities, n2.5 ± 1.2 13.4 ± 0.95.0 ± 1.3 14.8 ± 1.9 13.4 ± 1.4 1
Hypertension, %20 196 186 15528 1
Arterial Stiffness, %27 192 178 15329 1
Atherosclerosis, %331647 1426 1
Hyperglycemia, %8048 182 16754 1
Insulin resistance, %5 10 152 1110 1
Obesity, %282254 1250 1
Dyslipidemia, %260 179 14415 1
Renal Impairment, %22221428 13 1
Osteoporosis, %733 171524
Low muscle mass, %0 1116 13 180 1
Underweight, %0 10 10 15 180 1
Anxiety, %2 18 19 195 114
Depression, %90 19 159 114
Summary variables are presented as mean ± standard deviation for quantitative variables and percentage for discrete variables. 1 p < 0.05 versus the COPD group.
Table 3. Clinical characteristics across the (co)morbidity clusters in COPD.
Table 3. Clinical characteristics across the (co)morbidity clusters in COPD.
C1C2C3C4C5
Cluster descriptionLess ComorbidityCardiovascularMetabolicPsychologicCachectic
Subjects, n (%)46 (22%)27 (13%)56 (27%)40 (19%)39 (19%)
Male, %63676443 151
Age, years61 ± 765 ± 7 164 ± 661 ± 860 ± 6
BMI, kg/m227.2 ± 4.026.9 ± 3.531.1 ± 4.9 127.5 ± 3.819.9 ± 2.0 1
Waist circumference, cm102 ± 12102 ± 12 1112 ± 14102 ± 1582 ± 8 1
Former smoker, %5989 1755849 1
Current smoker, %3911 1264049 1
Pack-years, n39 (27–45) 140 (27–45)45 (32–57)58 (31–84) 143 (28–59)
Fram.10yr. risk ≥ 30%, %18 13760 12610 1
FEV1, % predicted51 ± 1649 ± 1656 ± 13 148 ± 1540 ± 15 1
FEV1/FVC41 ± 1240 ± 1046 ± 11 140 ± 1136 ± 12 1
ITGV, % predicted139 ± 28149 ± 30128 ± 30 1144 ± 30171 ± 32 1
KCO, % predicted67 ± 2268 ± 1780 ± 24 166 ± 1858 ± 18 1
MRC dyspnea grade3.0 ± 0.82.9 ± 1.12.9 ± 1.23.2 ± 0.93.0 ± 1.1
SGRQ, total score48.3 ± 15.0 156.1 ± 18.753.6 ± 15.166.2 ± 13.7 152.2 ± 16.2
6MWD, m445 ± 168450 ± 159462 ± 121466 ± 100465 ± 116
≥2 ex. previous year, %4633444853
Abbreviations: BMI = body mass index; Fram.10-yr. risk = Framingham 10-year risk of overall cardiovascular disease; FEV1 = Forced expiratory volume in 1 s; FVC = Forced vital capacity. ITGV = Intrathoracic gas volume. KCO = Transfer factor for carbon monoxide per alveolar volume. MRC = Medical Research Council. SGRQ = St. George’s Respiratory Questionnaire. 6MWD = 6-min walking distance. ex. = Exacerbations. kbp = Kilobase pairs. NA = Not applicable. Summary variables are presented as mean ± standard deviation for quantitative variables, with exception of pack-years (median and interquartile range), and percentage for discrete variables. 1 p < 0.05 versus COPD group.
Table 4. Prevalence of (co)morbidities across the (co)morbidity clusters in controls.
Table 4. Prevalence of (co)morbidities across the (co)morbidity clusters in controls.
X1X2X3
Cluster descriptionLess (Co)morbidityCardiovascularMetabolic
Subjects, n (%)97 (49%)34 (17%)69 (35%)
Comorbidities, n1.4 ± 1.0 1,23.8 ± 0.8 13.1 ± 1.3 1,2
Hypertension, %44 1,2100 1462
Arterial Stiffness, %1 1,2100 15 1,2
Atherosclerosis, %17 1,22832
Hyperglycemia, %46 1,26874 1
Insulin resistance %1 1019 1,2
Obesity, %0 1,22136 1
Dyslipidemia, %0 1,235 278 1
Renal Impairment, %4 21210
Osteoporosis, %96 23
Low muscle mass %763
Underweight, %530
Anxiety, %134
Depression, %1 203
Summary variables are presented as mean ± standard deviation for quantitative variables and percentage for discrete variables. 1 p < 0.05 versus the control group. 2 p < 0.05 versus the parallel COPD cluster.
Table 5. Clinical characteristics across the (co)morbidity clusters in controls.
Table 5. Clinical characteristics across the (co)morbidity clusters in controls.
X1X2X3
Cluster descriptionLess (co)morbidityCardiovascularMetabolic
Subjects, n (%)97 (49%)34 (17%)69 (35%)
Comorbidities, n1.4 ± 1.0 1,23.8 ± 0.8 13.1 ± 1.3 1,2
Male, %34 1,25954
Age, years60 ± 6 166 ± 6 160 ± 6 2
BMI, kg/m225.4 ± 2.4 1,227.4 ± 3.828.7 ± 3.5 1,2
Waist circumference, cm89 ± 10 1,295 ± 10 298 ± 11 1,2
Former smoker, %5362 254 2
Current smoker, %8 2915
Pack-years, n5 (0–16) 1,25 (0–22) 210 (0–25) 1,2
Fram.10yr. risk ≥30%, %9 133 129 1,2
FEV1, % predicted119 ± 15 2122 ± 16 2118 ± 15 2
FEV1/FVC78 ± 5 278 ± 5 279 ± 5 2
ITGV, % predicted106 ± 20 1,298 ± 17 294 ± 14 1,2
KCO, % predicted99 ± 15 2100 ± 14 2105 ± 18 1,2
Abbreviations: BMI = body mass index; Fram.10-yr. risk = Framingham 10-year risk of overall cardiovascular disease; FEV1 = Forced expiratory volume in 1 s; FVC = Forced vital capacity. ITGV = Intrathoracic gas volume. KCO = Transfer factor for carbon monoxide per alveolar volume. Summary variables are presented as mean ± standard deviation for quantitative variables, and percentage for discrete variables. 1 p < 0.05 versus the control group. 2 p < 0.05 versus the parallel COPD cluster.

Share and Cite

MDPI and ACS Style

Triest, F.J.J.; Franssen, F.M.E.; Reynaert, N.; Gaffron, S.; Spruit, M.A.; Janssen, D.J.A.; Rutten, E.P.A.; Wouters, E.F.M.; Vanfleteren, L.E.G.W. Disease-Specific Comorbidity Clusters in COPD and Accelerated Aging. J. Clin. Med. 2019, 8, 511. https://doi.org/10.3390/jcm8040511

AMA Style

Triest FJJ, Franssen FME, Reynaert N, Gaffron S, Spruit MA, Janssen DJA, Rutten EPA, Wouters EFM, Vanfleteren LEGW. Disease-Specific Comorbidity Clusters in COPD and Accelerated Aging. Journal of Clinical Medicine. 2019; 8(4):511. https://doi.org/10.3390/jcm8040511

Chicago/Turabian Style

Triest, Filip J. J., Frits M. E. Franssen, Niki Reynaert, Swetlana Gaffron, Martijn A. Spruit, Daisy J. A. Janssen, Erica P. A. Rutten, Emiel F. M. Wouters, and Lowie E. G. W. Vanfleteren. 2019. "Disease-Specific Comorbidity Clusters in COPD and Accelerated Aging" Journal of Clinical Medicine 8, no. 4: 511. https://doi.org/10.3390/jcm8040511

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop