Sleep-disordered breathing (SDB) has been associated with right-sided heart dysfunction and adverse cardiovascular outcomes. Longitudinal data are sparse in terms of understanding the prognostic implications of right ventricular remodeling in SDB on cardiovascular risk. We therefore investigated the predictive value of right-sided cardiac functional alterations on incident heart failure (HF) or death in SDB. Patients with SDB who underwent echocardiography within 1 month of index polysomnogram from January 2002 to July 2011 with normal left ventricular ejection fraction were included. Cox proportional prognostic hazard models predicting HF or death were used. Of a potential 375 subjects, 202 fulfilled the inclusion criteria (58 ± 14 years; 50% men). Subjects were followed for 3.1 ± 2.4 years with a total of 34 (16.8%) developing HF or death. Right ventricular end-systolic area (hazard ratio [HR] 1.3, 95% CI 1.01 to 1.6, p = 0.038), pulmonary vascular resistance (PVR; HR 1.4, 95% CI 1.1 to 1.7, p = 0.005) and also left atrial volume index (HR 1.7, 95%, CI 1.3 to 2.3, p <0.001) and E/A ratio (HR 1.4, 95% CI 1.1 to 1.7, p <0.001), were predictive of HF or death. Patients with increased PVR had significantly shorter event-free survival than without increased PVR (p = 0.04). In sequential Cox models, a model based on clinical data and left ventricular ejection fraction (χ 2 , 5.4) was improved by left atrial volume index (χ 2 , 12.7; p = 0.011) and further increased by PVR (χ 2 , 19.7; p = 0.015). In conclusion, right-sided heart dysfunction provides important prognostic information in SDB and may aid in identifying those at highest risk to target for closer follow-up.
It is recognized that in those with severe heart failure (HF) with reduced ejection fraction (EF), right-sided heart dysfunction contributes to increased mortality even independent of left ventricular (LV) function. Although right-sided heart abnormalities in sleep-disordered breathing (SDB) have been described in physiologic experiments and some epidemiologic work, existing data are limited by small sample sizes and an inability to characterize the longitudinal relations of right-sided heart dysfunction as it relates to development of morbidity and mortality in SDB. We therefore chose to investigate the relation of right-sided heart function and the development of adverse cardiovascular outcomes in SDB in those with baseline preserved EF in a clinic-based cohort. We postulated that in those with SDB, indexes of right-sided heart dysfunction might predict incident HF or death.
Methods
This is a longitudinal observational study of adult patients residing in Ohio and diagnosed with SDB by a sleep physician from April 2002 to July 2011 as a result of a minimum 1 overnight polysomnogram performed within the Cleveland Clinic Foundation. Inclusion criteria were as follows: (1) uninterrupted care within the Cleveland Clinic Foundation at least 6 months before and during the follow-up period after the polysomnogram and (2) completion of an echocardiographic examination within 1 month of the baseline polysomnogram recorded in the Cleveland Clinic prospective echocardiogram registry. Exclusion criteria included EF of <50%, patients with a history of HF, patients with moderate or more valvular regurgitation, those residing outside Ohio, and patients without an echocardiographic examination within 1 month of the polysomnogram evaluation. The entry date in the study was set at the date of polysomnogram, and all patients were followed up until incident HF or death occurred or until the last date of a physical examination performed by a physician within the Cleveland Clinic Foundation system up to September 2013. Institutional review board approval was obtained, and data were collected using multiple sources. Sleep evaluation included a sleep history and a routine overnight polysomnogram attended by a sleep technologist, reviewed by a board certified sleep physician. Apnea/hypopnea was defined by a minimum 50% decrease from the baseline airflow pressure signal excursions, lasting at least 10 seconds and with at least 3% desaturations from preevent baseline or/and association with an arousal (minimum 3-second awaking). Apnea/Hypopnea Index (AHI) was defined as the total number of respiratory events per hour of sleep, and SDB was defined based on an AHI ≥5 events/hour of sleep.
Comprehensive baseline transthoracic echocardiograms were performed with commercially available systems. Measurements performed for clinical purposes and available in the registry included left ventricular EF (LVEF), LV, and left atrial (LA) volumes (biplane Simpson method). Standard echocardiographic measurements of the right-sided heart function were made in accordance with current guidelines. Right ventricular (RV) fractional area change was defined using the formula ([end-diastolic area − end-systolic area]/end-diastolic area) × 100. Tricuspid annular plane systolic excursion was measured as the distance of systolic movement of the junction between the tricuspid valve and the RV free wall. Systolic pulmonary artery pressure was estimated from the maximal continuous-wave Doppler velocity of the tricuspid regurgitation jet using systolic transtricuspid pressure gradient calculated by the modified Bernoulli equation and the addition of estimated right atrial pressure. An index of pulmonary vascular resistance (PVR) was derived by dividing the maximal velocity of the tricuspid regurgitation jet by the RV outflow tract velocity − time integral. All measurements were performed and averaged over 3 cardiac cycles.
Baseline and follow-up clinical data were collected using multiple sources, and in conformity with definitions and recommendations published by the American College of Cardiology/American Heart Association in 2013: data elements and definitions and/or with 2010 American College of Cardiology/American Heart Association HF Performance Measurement Set. Clinical data available in the prospective echocardiographic registry were merged with data obtained during the baseline polysomnogram evaluation, by manual abstraction from comprehensive physician sleep study reports, and with socioeconomic and clinical characteristics recorded in Epic through intrainstitutional automated electronic data fields capturing all the aspects of care at any time before and during the follow-up period of the enrolled patients. Comprehensive assessments of co-morbidities at baseline were used as independent covariates for calculating the Charlson Comorbidity Index.
A composite outcome was calculated using incident HF or death. HF was collected from follow-up clinical data, and mortality was determined using telephone visits and the Social Security Death Index obtained from the Ohio Department of Health. HF was defined as meeting stage C or D criteria of classification of HF.
Standard statistical software packages (SAS version 9.3; SAS Institute, Inc., Cary, North Carolina and SPSS version 20.0; SPSS Inc., Chicago) were used for data analysis. “Time zero” was date of polysomnogram. The Kolmogorov–Smirnov test was performed to test the normality of continuous baseline cohort characteristics further summarized as mean ± SD or medians (interquartile range). Frequencies and percentages were calculated for the categorical data. Missing values (<15% and only covariables) were imputed using fivefold multiple imputations with the Markov chain Monte Carlo technique. Bootstrap bagging for variable selection and regression coefficients and their variance-covariance matrix for the final model were subsequently estimated for each imputed data set and further combined using the method of Rubin et al to produce the final estimates. Crude cumulative incidence function estimates of composite HF or death risk were plotted as a step function using the method of Fine and Gray. Cox proportional univariate and multivariate prognostic hazard models predicting HF and death were built. Sequential Cox models were performed to determine the incremental prognostic benefit of LA volume index (LAVi) and PVR over clinical data and LVEF, with incremental prognostic value being defined by a significant increase in global chi-square. Median values of LAVi and PVR were used to divide patients into 2 equal groups for Kaplan–Meier analysis, with survival compared using a 2-sided log-rank test.
Results
Of 375 consecutive patients, identified from billing records, with an echocardiographic examination done within 1 month of an overnight in-lab sleep study from April 2002 to July 2011, 285 patients were diagnosed by a sleep physician with SDB. Of those, 70 patients were excluded due to an HF co-morbidity, and an additional 13 patients were excluded because of reduced LVEF, more than mild valvular disease, or an insufficient follow-up period (<6 months). Baseline characteristics of the final cohort (n = 202) are included in Table 1 . They show an expected prevalence of risk factors and preserved LV function. Eighty-one patients (40%) had severe SDB (AHI ≥30 events/hour). Table 2 demonstrates the correlation between SDB severity (assessed by AHI) and echocardiographic characteristics of interest. Of note, there was no significant correlation between AHI and echocardiographic characteristics except for LV mass index which had a weak correlation with AHI ( r = 0.19, p = 0.007).
Variables | All |
---|---|
Age (years) | 57.7±13.8 |
Women | 101(50.0%) |
Married | 128 (63.4%) |
White | 148 (73.3%) |
Black | 423 (20.8%) |
Body surface area (m 2 ) | 2.3±0.4 |
Body mass index (kg/m 2 ) | 37.1±11.9 |
Systolic blood pressure (mmHg) | 129±17 |
Diastolic blood pressure (mmHg) | 77 ±45 |
Coronary artery disease risk factors | |
Hypertension | 128 (63.4%) |
Dyslipidemia | 119 (58.9%) |
Diabetes mellitus | 82(40.6%) |
Ever smoked | 100 (49.5%) |
Coronary artery disease | 36 (17.8%) |
Cerebrovascular disease | 23 (11.4%) |
Cardiomyopathy | 12 (5.9%) |
Chronic obstructive pulmonary disease | 51(25.3%) |
Obesity | 138 (68.3%) |
Overweight/Obesity | 182 (90.1%) |
Atrial Fibrillation and/or Flutter | 27 (13.4%) |
Depression | 54 (26.7%) |
Cancer ( No metastasis) | 45 (22.3%) |
Charlson co-morbidity index | 2.0±1.7 |
Medication | |
Angiotensin converting enzyme inhibitor | 56 (27.7%) |
Angiotensin II receptor blocker | 35 (17.3%) |
ß-blockers | 88 (42.9%) |
Calcium channel blocker | 61 (30.2%) |
Aspirin | 72 (35.6%) |
Insulin | 20 (9.9%) |
Apnea hypopnea index | 34.6±28.3 |
Arousal index (arousals/hour of sleep) | 34.5±20.8 |
Mean overnight oxygen saturation | 92±7 |
Minimum overnight oxygen saturation | 80±8 |
Periodic leg movement index | 3.1±8.5 |
Total pts Severe Sleep Disordered Breathing | 89 (44.1%) |
Total pts with Central Apnea | 2 (1.0%) |
Use of continuous positive airway pressure | 52 (25.7%) |
Left ventricular end-diastolic volume index (ml/m 2 ) | 43±13 |
Left ventricular end-systolic volume index (ml/m 2 ) | 17±7 |
Left ventricular ejection fraction (%) | 64±13 |
Left atrial volume index (ml/m 2 ) | 25±10 |
Left ventricular mass index (g/m 2 ) | 91±28 |
Trans-mitral early diastolic wave (cm/s) | 77±22 |
Trans-mitral early diastolic deceleration time (msec) | 220±60 |
Trans-mitral atrial filling wave (cm/s) | 76±23 |
Early diastolic / atrial filling ratio | 1.1±0.6 |
Right ventricular end diastolic area (cm 2 ) | 20.5±5.2 |
Right ventricular end systolic area (cm 2 ) | 11.6±3.6 |
Right ventricular fractional area change (%) | 44±7 |
Right atrial cavity area (cm 2 ) | 16.3±4.3 |
Tricuspid annular plane systolic excursion (cm ) | 2.1±0.4 |
Systolic pulmonary artery pressure (mmHg ) | 34.2±12.2 |
Pulmonary vascular resistance (Wood Units) | 1.6±0.5 |
Variables | Correlation | p value |
---|---|---|
Left ventricular ejection fraction (%) | -0.06227 | 0.3787 |
Left atrial volume index (ml/m 2 ) | -0.07753 | 0.2728 |
Left ventricular mass index (g/m 2 ) | 0.18785 | 0.0074 |
Left ventricular end-systolic volume index (ml/m 2 ) | 0.02603 | 0.7131 |
Left ventricular end-diastolic volume index (ml/m 2 ) | 0.03048 | 0.6667 |
Trans-mitral early diastolic wave (cm/s) | -0.08205 | 0.2457 |
Trans-mitral atrial filling wave (cm/s) | 0.00765 | 0.914 |
Early diastolic/atrial filling ratio | -0.0468 | 0.5084 |
Trans-mitral early diastolic deceleration time (msec) | -0.04273 | 0.546 |
Right ventricular end diastolic area (cm 2 ) | 0.11398 | 0.1063 |
Right ventricular end systolic area (cm 2 ) | 0.09695 | 0.1699 |
Right ventricular fractional area change (%) | 0.01714 | 0.8087 |
Tricuspid annular plane systolic excursion (cm) | -0.02357 | 0.7391 |
Pulmonary vascular resistance (Wood Units) | 0.00959 | 0.8923 |
All 202 patients were followed up over a period of 3.1 ± 2.4 years (range, 0.5 to 8.4 years). Thirty-four patients (16.8%) reached the composite end point defined as HF or death. There were no documented noncardiac deaths during follow-up. Incident HF occurred in 27 patients (13.4%), and 7 patients died (3.5%) during the study period. In our cohort, event-free survival rate was 98% (4 events total) and 97.5% (5 events in 2 years) at 1 and 2 years, respectively. Hazard ratios (HRs) of the relevant parameters before and after adjustment for age are listed in Table 3 . Interestingly, there was no association between AHI/use of continuous positive airway pressure and clinical outcome shown on Cox proportional hazards models. In multivariable Cox proportional hazard models (after adjustment for age and AHI, Table 4 ), RV end-systolic area (HR 1.3, p = 0.038) and PVR (HR 1.4, p = 0.005) were associated with the composite end point, independent of age and AHI.
Variables | Univariate Hazard ratio (95% CI) | p | Age ∗ adjusted Hazard Ratio (95% CI) | p |
---|---|---|---|---|
Apnea hypo index | 1.0 (0.99-1.01) | 0.787 | 1.0 (0.99-1.01) | 0.713 |
Use of continuous positive airway pressure | 1.9 (0.9-3.7) | 0.075 | 1.8 (0.9-3.6) | 0.102 |
Left ventricular ejection fraction | 0.8 (0.6-1.1) | 0.110 | 0.6 † (0.5-0.9) | 0.010 |
Left atrial volume index | 1.8 (1.4-2.4) | <0.001 | 1.7(1.3- 2.3) | <0.001 |
Left ventricular mass index | 1.5 (1.1-2.1) | 0.007 | 1.5(1.1- 2.0) | 0.016 |
Left ventricular end-systolic volume index, | 1.0 (0.7-1.4) | 0.912 | 1.1 (0.8- 1.5) | 0.625 |
Left ventricular end-diastolic volume index, | 1.1 (0.8-1.5) | 0.467 | 1.2 (0.9-1.7) | 0.219 |
Trans-mitral early diastolic wave, | 1.5(1.2-2.0) | 0.002 | 1.5(1.1-2.0) | 0.005 |
Trans-mitral atrial filling wave, | 0.7(0.5-1.1) | 0.128 | 0.7(0.5- 0.99) | 0.041 |
Early diastolic / atrial filling ratio | 1.5 (1.2- 1.8) | <0.001 | 1.4 (1.2- 1.7) | <0.001 |
Trans-mitral early diastolic deceleration time | 1.1( 0.8- 1.5) | 0.685 | 1.0 ( 0.7- 1.5) | 0.875 |
Right ventricular end systolic area | 1.2 (0.99 -1.5) | 0.059 | 1.3 (1.02- 1.6) | 0.035 |
Right ventricular end diastolic area | 1.2 ( 0.9-1.6) | 0.142 | 1.3 (0.97- 1.7) | 0.084 |
Right ventricular fractional area change | 0.8 (0.6- 1.04) | 0.095 | 0.8 (0.6- 1.1) | 0.118 |
Tricuspid annular plane systolic excursion | 0.8 ( 0.5-1.1) | 0.098 | 0.7 ( 0.5-1.1) | 0.092 |
Pulmonary vascular resistance, | 1.3(1.1- 1.6) | 0.006 | 1.4 (1.1-1.7) | 0.005 |