Article Text

Download PDFPDF

Original article
Polysomnographic phenotypes and their cardiovascular implications in obstructive sleep apnoea
Free
  1. Andrey V Zinchuk1,
  2. Sangchoon Jeon2,
  3. Brian B Koo3,
  4. Xiting Yan1,
  5. Dawn M Bravata4,
  6. Li Qin5,
  7. Bernardo J Selim6,
  8. Kingman P Strohl7,
  9. Nancy S Redeker2,
  10. John Concato1,8,
  11. Henry K Yaggi1
  1. 1 Department of Medicine, Yale University School of Medicine, New Haven, Connecticut, USA
  2. 2 Division of Acute Care/Health Systems, Yale School of Nursing, Yale University, New Haven, Connecticut, USA
  3. 3 Department of Neurology, Yale University, New Haven, Connecticut, USA
  4. 4 Departments of Neurology and Internal Medicine, Richard L. Roudenbush VA Medical Center, Indianapolis, Indiana, USA
  5. 5 Center for Outcomes Research and Evaluation, Yale University School of Medicine, New Haven, Connecticut, USA
  6. 6 Section of Pulmonary, Critical Care, and Sleep Medicine, Mayo Clinic, Rochester, Minnesota, USA
  7. 7 Section of Pulmonary, Critical Care, and Sleep Medicine, Case Western Reserve University, Cleveland, Ohio, USA
  8. 8 Clinical Epidemiology Research Center, Veterans Affairs Connecticut Healthcare System, West Haven, Connecticut, USA
  1. Correspondence to Dr Henry K Yaggi, Department of Medicine, Yale University School of Medicine, 300 Cedar Street, New Haven, CT 06443, USA; henry.yaggi{at}yale.edu

Abstract

Background Obstructive sleep apnoea (OSA) is a heterogeneous disorder, and improved understanding of physiologic phenotypes and their clinical implications is needed. We aimed to determine whether routine polysomnographic data can be used to identify OSA phenotypes (clusters) and to assess the associations between the phenotypes and cardiovascular outcomes.

Methods Cross-sectional and longitudinal analyses of a multisite, observational US Veteran (n=1247) cohort were performed. Principal components-based clustering was used to identify polysomnographic features in OSA’s four pathophysiological domains (sleep architecture disturbance, autonomic dysregulation, breathing disturbance and hypoxia). Using these features, OSA phenotypes were identified by cluster analysis (K-means). Cox survival analysis was used to evaluate longitudinal relationships between clusters and the combined outcome of incident transient ischaemic attack, stroke, acute coronary syndrome or death.

Results Seven patient clusters were identified based on distinguishing polysomnographic features: ‘mild’, ‘periodic limb movements of sleep (PLMS)’, ‘NREM and arousal’, ‘REM and hypoxia’, ‘hypopnoea and hypoxia’, ‘arousal and poor sleep’ and ‘combined severe’. In adjusted analyses, the risk (compared with ‘mild’) of the combined outcome (HR (95% CI)) was significantly increased for ‘PLMS’, (2.02 (1.32 to 3.08)), ‘hypopnoea and hypoxia’ (1.74 (1.02 to 2.99)) and ‘combined severe’ (1.69 (1.09 to 2.62)). Conventional apnoea–hypopnoea index (AHI) severity categories of moderate (15≤AHI<30) and severe (AHI ≥30), compared with mild/none category (AHI <15), were not associated with increased risk.

Conclusions Among patients referred for OSA evaluation, routine polysomnographic data can identify physiological phenotypes that capture risk of adverse cardiovascular outcomes otherwise missed by conventional OSA severity classification.

  • obstructive sleep apnea (OSA)
  • cluster analysis
  • phenotype
  • heterogeneity
  • cardiovascular diseases
  • mortality

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Key messages

What is the key question?

  • Can routine polysomnographic data and unsupervised learning methods be used to identify subgroups (‘phenotypes’) of obstructive sleep apnoea (OSA) patients, and are those subgroups associated with increased risk of adverse cardiovascular outcomes or death?

What is the bottom line?

  • In patients referred for OSA evaluation, we identified physiological phenotypes that captured risk of adverse cardiovascular outcomes otherwise missed by conventional OSA severity classification using the apnoea–hypopnoea index.

Why read on?

  • We show that physiologically distinct groups of patients (eg, with higher prevalence of periodic limb movements of sleep, arousals, hypoxia and so on) exist within each of the traditional OSA severity categories and that these groups may have implications for treatment and risk of cardiovascular events or death.

Introduction

Obstructive sleep apnoea (OSA) is increasingly recognised as a complex and heterogeneous syndrome.1 Categorising such conditions into smaller, more homogeneous categories (‘phenotypes’) can help advance understanding of pathophysiology, develop targeted treatments and improve both prognostication and risk stratification for clinical trials. Polysomnography enables collection of multiple types of physiological data, which results in numerous metrics by which patients with OSA can be characterised. Previously suggested OSA phenotypes were based primarily on occurrence of respiratory events during specific sleep states,2 3 body position4 or hypoxia.5 These prior studies used approaches that identified the phenotypes a priori, based on a single or few characteristics, and therefore could have missed more complex and (up to this point) unidentified OSA subgroups.6

The clinical relevance, including prognostic implications, of such polysomnogaphic phenotypes remains largely uncertain. OSA is a well-known independent risk factor for cardiovascular disease (CVD) and death.7–9 Much of the epidemiological work linking OSA and CVD outcomes, however, is focused primarily on the apnoea–hypopnoea index (AHI),10 which represents a small fraction of the physiological data generated by polysomnography. Recent studies suggest that the AHI alone may not fully capture the physiological processes (eg, hypoxaemic burden,11 sleep fragmentation12 or periodic limb movements of sleep (PLMS))13 implicated in development of adverse cardiovascular outcomes.

Cluster analysis, a type of unsupervised learning methodology, can integrate multiple characteristics (such as polysomnographic metrics) without a priori groupings to identify unique patient categories.14 In other heterogeneous disorders such as asthma, COPD and heart failure, phenotypes with different prognoses and treatment responsiveness were identified using cluster analysis.15–17 In OSA, research using unsupervised learning methods to identify patient subgroups based on polysomnographic metrics is scant.18 19 Associations of such subgroups with cardiovascular outcomes remain largely unexplored.

Our objective was to use routine polysomnographic data to capture OSA’s physiological heterogeneity in relation to clinically relevant cardiovascular outcomes. Specifically, we hypothesised that unique clusters (phenotypes) of patients could be identified by applying unsupervised learning methods to polysomnographic data and that the clusters would be differentially associated with risk of adverse cardiovascular outcomes (acute coronary syndrome (ACS), transient ischaemic attack (TIA), stroke or death.

Methods

A comprehensive description of methodological details is provided in the online supplementary file.

Supplementary file 1

Participants and measures

The Determining Risk of Vascular Events by Apnea Monitoring (DREAM) study consists of 2041 veterans from three Veterans Affairs Centers (West Haven, Connecticut; Cleveland, Ohio; and Indianapolis, Indiana) that underwent OSA evaluation between 2000 and 2004, with follow-up through 201220. At enrolment (ie, time of polysomnography), data were collected on demographic and anthropometric characteristics, cardiovascular risk factors, medical comorbidities and other measures, as detailed previously.20 Cardiovascular risk was ascertained using the Framingham risk score (FRS; including age, sex, smoking status, presence of diabetes, antihypertensive medication, systolic blood pressure, total cholesterol and high-density lipoprotein levels)21 and comorbidity burden by the Charlson comorbidity index (CCI).22 Polysomnographic scoring, variable definitions, classification and continuous positive airway pressure (CPAP) treatment adherence (regular vs not regular) are described in the online supplementary file. All polysomnographic scoring was performed at a single centre (West Haven, Connecticut). The primary outcome was the composite of incident ACS, TIA, stroke or death from any cause. Rationale for composite outcome selection and its adjudication are described previously20 and in the online supplementary file.

Polysomnographic variable selection and analytic sample

Sixty-five polysomnographic variables were sorted into four domains, representing the pathophysiological processes previously reported to be in the causal pathway between OSA and CVD: breathing disturbance, autonomic dysregulation, hypoxaemia and sleep disturbance.10 23 Redundant variables (eg, time in bed) were excluded if other included variables reflected similar physiological parameters (eg, total sleep time and efficiency). Categorical variables (eg, snoring), and variables with >5% missing data, were excluded due to cluster analysis requirements (see online supplementary file for details). Twenty-nine polysomnographic variables (see online supplementary table E1) were ultimately selected for variable reduction analysis.

From the source population (n=2041), the final analytic sample (n=1247) was obtained. Patients with missing data on any polysomnographic variables (n=482) were excluded due to cluster analysis requirements. Patients with inconsistently abstracted polysomnograms or scoring discrepancies (n=308), or those with <15 min of sleep (n=4) were also excluded.

Statistical analyses

A two-step variable reduction analysis was used to identify unique polysomnographic features within each domain. First, principal components-based clustering was used to identify groups of variables highly correlated with their own cluster and uncorrelated with others.24 Second, polysomnographic features were selected to retain >75% of total variance within each domain, using the least number of features that were judged to be clinically interpretable (see online supplementary figure E1 and table E2).

We employed K-means analysis to generate clusters, based on polysomnographic feature scores. Silhouette widths (a cluster quality measure)25 and clinical interpretability were used to select the number of clusters (detailed description in online supplementary file). Jaccard coefficient of bootstrap cluster samples was used to evaluate cluster stability.26

We assessed cross-sectional contrasts between clusters, using analysis of variance or X2 tests. The relationship between OSA clusters and primary outcome was evaluated with unadjusted Cox proportional hazards models and with models adjusting for FRS and CPAP use. The proportional hazards assumption for the model comparing each cluster to the ‘mild’ group was assessed using the cumulative hazard function: H(t)=−logS(t). Sensitivity analyses incorporated the clusters and CPAP use with (1) CCI and (2) home oxygen use at baseline (all sensitivity analyses are reported in the online supplementary file). Analyses were also performed with CVD-specific mortality as part of the composite outcome. We repeated the analysis using AHI severity categories (mild/none: AHI <15; moderate: 15≤ AHI<30; and severe: AHI ≥30), in place of cluster membership, to assess whether risk was captured by conventional OSA severity categories. To evaluate for differences between the analytic sample and excluded patients, we compared the baseline risk (FRS and CCI) and primary outcomes between these groups (online supplementary file). SAS V.9.3 was used for all analyses. Institutional review board at each centre approved the study (WHVA HIC #0001).

Results

The study population was predominantly male, with an average age of 58.3±11.7 years, mean AHI of 25.0±29.7 events/hour and 29% prevalence of CVD (table 1). The 29 polysomnographic variables were reduced to 17 features that preserved 83% of total variance (see online supplementary table E3) and were judged to be clinically interpretable. Patient-level cluster analysis resulted in seven OSA subgroups (clusters A - G) supplementary clusters A–G), labeled based on the most distinguishing features (table 2). Feature distributions are shown in online supplementary figure E1 and table E6. Given that OSA severity is conventionally expressed according to AHI-based categories (none/mild, moderate and severe) and to highlight heterogeneity within each category, cluster descriptions below are grouped by conventional OSA severity.

Table 1

Summary of patient characteristics in the analytic sample (n=1247) at baseline

Table 2

Description of and labels for the polysomnographic clusters based on distinguishing features

None/mild OSA

The largest (n=533) cluster was labelled ‘mild’, based on few disordered breathing events (median AHI (IQR) of 4 (9)), high sleep efficiency, normal REM sleep and near normal nocturnal oxygen saturation (table 3). The ‘PLMS’ cluster (n=119), with median AHI of 10 (19), was distinguished from others primarily by a markedly higher PLMS index (PLMI). Compared with the ‘mild’ cluster, those in ‘PLMS’ spent more sleep time with oxygen saturation of <90% (T<90%).

Table 3

Polysomnographic characteristics of OSA clusters (ranked by AHI, A–G) grouped according to domain (breathing disturbance, hypoxaemia, sleep architecture disturbance and autonomic dysregulation)*

Patients in the ‘mild’ cluster were among the youngest, with lowest comorbidity (CCI) and CPAP use (table 4). ‘PLMS’ cluster members were on average 8 years older than the cases in the ‘mild’ group and exhibited higher comorbidity, including hypertension and diabetes (table 4, online supplementary table E1).

Table 4

Baseline demographic, anthropometric, sleep and comorbidity characteristics for clusters grouped by conventional OSA severity designations (based on median cluster AHI: none/mild (Clusters A and B), moderate (Clusters C and D) and severe (Clusters E, F and G))

Moderate OSA

The median AHIs for ‘NREM and poor sleep’ (n=186) and ‘REM and hypoxia’ (n=168) clusters were equivalent, but the polysomnographic measures differed markedly. NREM-related apnoeas, hypopnoea-related and spontaneous arousals characterised ‘NREM and poor sleep’ cluster (table 3, online supplementary table E6). These patients exhibited poor sleep (55% efficiency; 41% of stage 1 with low REM sleep). Patients in the ‘REM and hypoxia’ cluster exhibited the highest frequency of REM-related apnoeas, lowest oxygen saturation nadir (72.8%) and an increased hypoxic burden (median (IQR) T<90% of 9% (23%), table 3).

Patients in ‘NREM and poor sleep’ and ‘REM and hypoxia’ clusters were more obese, exhibited higher CPAP or home oxygen use and had higher prevalence of hypertension and diabetes, in comparison with the ‘mild’ cluster (table 4, online supplementary table E7).

Severe OSA

Cluster labelled ‘hypopnea and hypoxia’ (n=75, median AHI of 44 (39)) showed predominance of hypopnoeas with desaturations and hypoxic burden (T<90% of 14% (36%)), and relatively preserved sleep architecture (table 3, online supplementary table E6). Dominance of apnoeas with arousals, but without desaturations, and markedly disturbed sleep architecture distinguished the ‘arousal and poor sleep’ cluster (n=42, median AHI of 68 (30)). Its hypoxic burden (T<90% of 0% (4%); oxygen saturation nadir of 87.1%) was similar to clusters with markedly lower AHIs. The ‘combined severe’ cluster (n=126, median AHI of 84 (28)) exhibited combined apnoeas (with desaturations and arousals), severe hypoxic burden, high arousal indices and high percentage of stage 1 sleep.

Patients in the ‘arousal and poor sleep’ cluster tended to be older, less sleepy and have lower CPAP use, compared with those in other clusters in the severe OSA category (table 4). This cluster also had the highest percentage of Black and Hispanic patients. Although baseline comorbidity burden was similar for clusters in the ‘severe OSA’ category, patients in ‘arousal and poor sleep’ cluster had a trend towards higher prevalence of myocardial infarction, atrial fibrillation and chronic lung disease in comparison to the ‘mild’ cluster (online supplementary table E7).

Clinical outcomes

After an average follow-up of 4.9±1.9 years, 271 patients experienced ACS, TIA, stroke or death from any cause (online supplementary table E8). Statistically significant differences in Kaplan-Meier curves for the primary outcome were evident according to patient cluster (log rank, p value <0.001, figure 1). In unadjusted analysis (table 5), cluster membership was associated with adverse cardiovascular events or death (p value 0.002). Compared with the ‘mild’ cluster, the ‘PLMS’ cluster exhibited a greater than twofold higher risk for the primary outcome (hazard ratio (HR) 2.36; 95% CI 1.61 to 3.46). Among the remaining clusters, only ‘REM and hypoxia’ did not exhibit a statistically higher risk of primary outcome (HR 1.34, 95% CI 0.91 to 1.99) than ‘mild’ cluster. Overall, regular CPAP use was associated with a 36% decreased risk of primary outcome (p value 0.0015, table 5) and this attenuation varied by cluster (figure 2). After adjustment with FRS for age, sex and baseline cardiovascular risk, the association with the primary outcome remained significant for ‘PLMS,’ ‘hypopnea and hypoxia’ and ‘combined severe’ clusters, with the HR of 1.79, 95% CI 0.97 to 3.29 for the ‘arousal and poor sleep’ group (table 6). Because majority of combined outcome events were deaths, sensitivity analyses using CCI in place of FRS, as well as the addition of home oxygen use were performed with similar results (online supplementary tables E9 and E10). Finally, a sensitivity analysis using CVD deaths versus all-cause mortality as part of the composite outcome resulted in equivalent risk associations for the clusters (data not shown).

Figure 1

Kaplan-Meier survival probability curves for risk of acute coronary syndrome, transient ischaemic attack, stroke or death for seven polysomnographic clusters. NREM, non-rapid eye movement; PLMS, periodic limb movements of sleep; REM, rapid eye movement.

Figure 2

Risk (odds ratio, OR) of ACS, TIA, stroke or death for those with regular versus not regular CPAP use for each cluster. ACS, acute coronary syndrome; CPAP, continuous positive airway pressure; NREM, non-rapid eye movement; PLMS, periodic limb movements of sleep; REM, rapid eye movement; TIA, transient ischaemic attack. p value is for Breslow-Day homogeneity of ORs test. OR and 95% CI reported for those with regular CPAP use versus not regular CPAP use.

Table 5

Cox hazard ratios (HR) for incident acute coronary syndrome, transient ischaemic attack, stroke or death for OSA clusters (n=1036*)

Table 6

Comparison of adjusted Cox proportional hazards models for primary outcome using cluster membership versus OSA severity as measured by AHI (n=1036*, see legend for analyses by none vs mild, moderate and severe groups)

Conventional OSA severity categories, based on AHI, were not associated with the primary outcome (p value of 0.49) in an adjusted model, regardless of the reference group specification (none (AHI <5) or none/mild (AHI <15), table 6).

Discussion

Using readily available polysomnographic data and unsupervised analytic methods (eg, cluster analysis), we identified multiple patient phenotypes with mean AHIs within each of the conventional, AHI-based OSA severity categories. These clusters captured different phenomena, including co-occurrence with PLMS, association of respiratory events with arousals, hypoxic load and specific sleep states. To highlight the clinical relevance of physiological heterogeneity portrayed by the clusters (but not to create a prognostic tool), we evaluated their associations with incident cardiovascular events or death. We found that cluster membership independently predicted adverse CVD or death, whereas conventional AHI-based severity categories did not.

PLMS were the main distinguishing characteristic between ‘mild’ and ‘PLMS’, the two clusters with median AHI <15. The high prevalence of PLMS in patients with OSA (28%–48%), especially those with AHI <20,27 supports our finding of a ‘PLMS’ cluster as a group with a distinct physiological disturbance in OSA. Although PLMS are common, little is known about factors associated with their origin in patients with OSA, aside from older age and increased comorbidity burden,27 also found in the ‘PLMS’ cluster. Evidence of the adverse impact of PLMS on the cardiovascular system is emerging. Basal autonomic regulation is shifted towards sympathetic predominance in OSA patients with PLMS, in comparison with those with OSA alone.28 PLMS are associated with daytime blood pressure elevations29 and progression to persistent atrial fibrillation.30 After controlling for baseline cardiovascular risk, we found that the ‘PLMS’ cluster was associated with a twofold increase of the combined outcome compared with the ‘mild’ cluster and that CPAP use may attenuate this risk (figure 2). In another OSA cohort, PLMI independently predicted incidence of CVD or death from any cause.13 Taken together with our findings, evidence suggests that PLMS are important for identification of patients at increased risk of CVD or death.

In the moderate OSA severity range (15≤ AHI<30), the polysomnographic characteristics of the ‘NREM and poor sleep’ cluster suggest that a lower arousal threshold (in addition to higher ventilatory loop gain) may result in fragmented sleep with less severe oxygen desaturation.31 Although we did not measure arousal threshold directly, this cluster exhibited a nadir oxygen saturation of >82.5%, a previously reported predictor of a low arousal threshold.32 In contrast, the ‘REM and hypoxia’ cluster, with predominance of REM-related respiratory events, severe hypoxia, but relatively non-interrupted sleep, was similar to REM-related OSA, as a previously recognised phenotype.2 33 In our study, the ‘NREM and poor sleep’ cluster was not associated with increased composite outcome, and to our knowledge, incidence of CVD or death in patients with polysomnographic characteristics similar to ‘NREM and poor sleep’ cluster has not been described. The ‘REM and hypoxia’ cluster was associated with higher (unadjusted) prevalence of hypertension (online supplementary table E7), but not with increased incidence of ACS, TIA, stroke or death. One prior study of REM-related OSA found a lower prevalence of hypertension in comparison with non-stage-specific OSA,2 whereas another found increased prevalent and incident hypertension in OSA patients with REM AHI >15.33 The differences in study design, outcomes and cohort characteristics (eg, predominantly male) make direct comparisons difficult.

Our findings suggest there is considerable pathophysiological variability among patients with AHI ≥30. Some patients may develop hypoxia at lower intensity of airway obstruction (‘hypopnea and hypoxia’), compared with others who may require frank apnoeas for desaturation (‘combined severe’). Patients with equivalent AHIs and conventional arousal indices may exhibit markedly different degrees of hypoxia (‘arousal and poor sleep’ vs ‘combined severe’). Some arousals precede airway opening, acting as a ‘safety valve’ to prevent worsening hypoxia/hypercapnia, whereas others occur after the airway opening.31 In patients who open their airways prior to or without an arousal, spontaneous arousals may potentiate ventilatory instability, interfere with progression to deeper sleep and increase cycling frequency (more frequent events with less hypoxia).31 Measuring factors such as arousal timing, apnoea and hypopnoea length, pre-event oxygen saturation, ventilatory loop gain and pharyngeal muscle responsiveness34 may help refine the differences between clusters and provide insight into these categories of risk.

The ‘hypopnea and hypoxia’ and ‘combined severe’ clusters in our study were each associated with increased risk for the composite outcome, in contrast to the AHI-based severity categories. Although many studies have reported an association between AHI and cardiovascular events or death,7–9 a recent systematic review shows that this association has not been found consistently.10 and, such as T<90% or a number of awakenings, have been identified as independent predictors.7 11 13 It is plausible that different physiological characteristics may play an important role in some patients, but not in others. For example, for patients in the ‘hypopnea and hypoxia’ and ‘combined severe’ clusters, the duration and depth of hypoxia (not simply desaturation frequency) may pose increased risk. In contrast, for other patients (‘arousal and poor sleep’), with more non-white and older individuals, the rapid cycling apnoeas associated with arousals and fragmented sleep, but without marked hypoxia, may be important. It is possible that when categorising patients using AHI alone, we may group patients with different physiological characteristics and risk profiles (online supplementary figure E3), potentially missing an opportunity for risk mitigation.

Our findings may have important implications for selecting patients for treatment and clinical trials in OSA. In current practice guidelines, for example, the recommendation for treatment of mild OSA (AHI <15) without symptoms is optional, and PLMS are not included in risk stratification.35 Focusing on AHI alone may miss a subgroup of patients at greater risk of CVD or death. CPAP use attenuated risk for the ‘PLMS’ cluster in our study. Limited evidence regarding CPAP treatment for unselected patients with AHI <15 suggests no impact on fatal or non-fatal cardiovascular events.36 Our findings may also serve as a guide for recruitment into clinical trials of OSA treatment, given that patients at greatest risk of adverse outcomes may derive the greatest benefit.

Among patients traditionally classified as severe OSA, our results suggest different physiological mechanisms of sleep-disordered breathing. As such, these patients may be better served with targeted treatment based on the predominant causative physiological factors (eg, reducing arousal propensity with sedatives and/or ventilatory drive with oxygen)34 rather than a singular approach of PAP therapy. In addition, not all patients with severe OSA may stand to benefit from treatment (figure 2). The impact of using such approaches on cardiovascular outcomes is not known and is an important area of future investigation.

Although cluster analysis has been used in OSA,18 19 37–40 only two studies18 19 included polysomnographic measures other than AHI. Joosten et al 18 identified six clusters in 1064 patients with mild or moderate OSA, similar to a priori defined polysomnographic OSA groups (REM-predominant, supine-isolated, supine-predominant and REM/supine overlap). Most metrics in the domains of sleep architecture, arousals and hypoxia were excluded. Analysis by Lacedonia et al, incorporating AHI, T<90% and oxygen desaturation index, identified three clusters in 198 patients with OSA.19 Two clusters with AHI >55 differed markedly in T<90%, akin to the ‘combined severe’ and ‘arousal and poor sleep’ groups in our analysis. None of the clusters differed in prevalence of CVD or diabetes. Neither study focused on or modelled clinical outcomes by cluster, limiting the prognostic implications of the findings.

Several limitations of our study warrant consideration. Data from a single prior study were analysed (although from multiple VA medical centres), and as with most cluster analyses, our results are best used to generate hypotheses. Validation of our findings in other populations to identify distinguishing polysomnographic features is needed. Although our overall analytic sample was relatively large, because individual groups in cluster analyses may be small, we elected (a priori) not to perform split-cohort validation to avoid ‘missing’ potentially clinically relevant but small patient categories (eg, ‘hypopnea and hypoxia’, n=75). Similarly, better understanding of the cluster’s implications is needed. For example, whether the ‘PLMS’ cluster represents restless leg syndrome with comorbid OSA, or a unique entity of OSA with PLMS, as well as its significance to cardiovascular outcomes and death, needs to be explored. Given an extensive number of polysomnographic variables included in our study (n=29), and requirements of no missing data for cluster analysis, we anticipated a considerable proportion of participants from the source populations would be excluded from our analytic sample. A comparison of demographic factors, baseline cardiovascular risk comorbidities and primary outcomes between these groups revealed no differences (online supplementary table E11). In addition, our data did not include measures of symptoms, intermediate OSA phenotypes (eg, craniofacial structure and physiological traits) or cardiovascular risk biomarker data. The current study was also not designed to address the objective adherence to or effectiveness of CPAP on cardiovascular outcomes among the clusters. Nonetheless, we demonstrated increased risk of adverse events despite regular use of CPAP, suggesting that for some groups, risk was present even with treatment. Finally, the lack of gender and ethnic diversity in our cohort may limit generalisability of our findings to women or minorities.

Strengths of our study include use of a structured, data-driven approach to ascertain clinically interpretable polysomnographic features within each pathophysiological domain of OSA, and rigourous quantitative methodology (ie, cluster analysis), integrating multiple polysomnographic features to identify unique patient subgroups. Our approach is easily translatable across centres, as it is based on standard polysomnography and scoring. Assessing the clusters’ incident cardiovascular outcomes or death, while adjusting for baseline cardiovascular risk using validated measures, is another strength. We also included patients with a full range of sleep apnoea severity, given that some with AHI <5 have been previously shown to have increased risk of cardiovascular events.33 Our results serve to emphasise the need for multidimensional OSA characterisation to improve prognostication and trial design, especially as large, multidomain datasets from multiple study populations become available.

In summary, our study adds new evidence that polysomnographic heterogeneity exists within traditional OSA severity categories and that readily available data with unsupervised learning methods can be used to classify patients into subgroups reflecting different pathophysiological processes. Some of these subgroups capture the risk of adverse cardiovascular events or death otherwise missed by conventional AHI categorisation, suggesting that additional measures—such as PLMS or respiratory events with hypoxia or arousals—may be important for risk stratification and treatment selection. Future studies are needed to evaluate whether unsupervised learning methods can be applied in other populations and multidimensional datasets to better understand pathophysiology, improve prognostication and identify targeted treatments for patients with OSA.

Acknowledgments

The authors are grateful to US Veterans, staff at the Clinical Epidemiology Research Center at Connecticut’s Veterans Affairs Medical Center and faculty and fellows of Yale’s Section of Pulmonary, Critical Care and Sleep Medicine who provided valuable feedback in formulating this work.

References

Footnotes

  • Contributors All authors contributed substantially to the conception, design, analysis or interpretation of the data in this study. AVZ, SJ and HKY had full access to all the data in the study and take responsibility for the integrity and the accuracy of the data analyses. AVZ, SJ, BBK, NSR, JC and HKY were involved in the interpretation of data. AVZ, BBK and HKY drafted the manuscript and all authors revised it critically for important intellectual content. All authors gave final approval of this version to be submitted.

  • Funding This work was supported by the VA Clinical Science Research and Development (CSR&D) Merit Review Program (CSRDS07), VA Cooperative Studies Program, NRSA Institutional training Grant from the NIH (5T32HL07778), Yale Center for Investigating Sleep Disturbance in Acute and Chronic Conditions (P20NRO14126) and Robert E. Leet and Clara Guthrie Patterson Trust Fellowship Program in Clinical Research, Bank of America, N.A., Trustee. JC is supported by the VA Cooperative Studies Program.

  • Competing interests None declared.

  • Ethics approval West Haven Veterans Affairs Institutional Review Board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Presented at Part of this work was presented at as a poster discussion presentation: Zinchuk AV, Jeon S, Koo BB, et al. Clinically Relevant Phenotypes of Obstructive Sleep Apnea. Am J Respir Crit Care Med 2016;193:A6380.

Linked Articles