Sleep Apnea, Obstructive - Social Class - Disorders of Excessive Somnolence - Polysomnography
INTRODUCTION
Obstructive sleep apnea (OSA) is a disorder characterized by recurrent episodes of
partial or complete obstruction of the upper airway during sleep. Its increasing prevalence
is consistent with the frequency of established risk factors, such as older age and
obesity[1]. Functional outcomes, such as excessive daytime sleepiness (EDS) and poor concentration,
largely reversible with successful treatment, may affect occupational performance[2] and increase the risk of road traffic[3] and occupational accidents[4]. On the other hand, occupation and socio-economic status (SES) have been described
in recent research as independent prognostic factors for OSA presence[5], severity[6], cardiovascular consequences[7] and adherence to treatment[8], although controlling for potential confounders remains a matter for consideration.
The objectives of the current study were to investigate the relationship between OSA
severity and SES, measured by occupational activity, in a sample referred for sleep
study in a public hospital sleep laboratory and to corroborate whether possible differences
are explained by other known risk factors, such as age, sex, obesity, smoking and
alcohol consumption. We have hypothesized that lower class workers would have greater
OSA severity, mainly as a result of their greater exposures to unhealthy habits and
behaviors[9].
METHODS
A retrospective cross-sectional study was conducted in the setting of a public-hospital-based
sleep laboratory in Athens. All referred patients were evaluated by a sleep specialist
physician in the outpatient department, where they were asked to complete a questionnaire
on demographics, social history, sleep habits, nighttime and daytime symptoms and
other health problems, as well as the Greek version of the Epworth sleepiness scale
(ESS)[10].
Those judged as having high probability of OSA underwent a diagnostic polysomnography
(PSG) with a digital recording system (Alice 3, 4 and 5 Diagnostic Sleep System, Philips
Respironics) that monitors the following variables: electroencephalogram (two paracentral,
two frontal and two occipital leads), right and left electrooculograms, submental
and anterior tibial electromyograms, electrocardiogram, nasal (nasal pressure transducer)
and oronasal airflow (thermistor), respiratory effort (thoracoabdominal respiratory
inductance plethysmography belts) and pulse oximetry.
Sleep and respiratory event scoring was done manually in 30-second epochs according
to recommendations by American Academy of Sleep Medicine[11]. An apnea was scored when a reduction in the oronasal flow signal by ≥90% from baseline
was recorded for at least 10s. Apneas were classified as obstructive or mixed in the
presence of inspiratory effort in the entire event or a part of it respectively. Hypopnea
was scored when a reduction in the nasal pressure signal by ≥30% from baseline was
recorded for at least 10s, accompanied with a ≥3% oxygen desaturation from baseline
or an arousal. The study protocol was approved by the Ethics Committee of the National
School of Public Health.
The records of all patients who underwent diagnostic PSG between January 1 and December
31, 2015, were searched and analyzed. Titration, follow-up and studies with total
sleep time less than 2 hours were excluded. The diagnosis of OSA was made according
to the diagnostic criteria in the third edition of the International Classification
of Sleep Disorders[12], that require the presence of 15 or more predominantly obstructive respiratory events
per hour of sleep or the combination of 5 or more obstructive events per hour with
appropriate symptoms and comorbid conditions. Symptomatology and comorbidities were
assessed in each patient by their answers in the questionnaire. OSA severity was evaluated
on the basis of the following parameters: frequency of respiratory events with calculation
of the apnea-hypopnea index (AHI) and degree of subjective daytime sleepiness with
ESS score.
The European Socio-economic Classification (ESeC) was used as an indicator of SES.
It is a categorical social stratification scheme based on the idea that in market
economies, the market position and in particular the occupational division of labor
is fundamental to the production of social inequalities. Classification is performed
according to occupation, employment status (employer, manager, supervisor, employee,
self-employed) and the size of the organization[13].
Subjects were first given a 3-digit code, categorizing them in minor occupational
groups of the International Standard Classification of Occupations (ISCO-88)[14] based on their response about their current or last professional activity in the
questionnaire. This code was used to derive the ESeC ranking under the simplified
method, since information about employment status and size of organization was unavailable,
which has approximately 80% agreement with the full model. The classes were then collapsed
in the 3-class version, “salariat”, “intermediate” and “working class”, for analysis
purposes. Housewives and students were excluded from statistical analysis, since they
are not actively seeking work and their SES is influenced by the occupation and income
of other family members.
Variables that were considered potential confounders included sex, age, anthropometric
characteristics, such as body mass index (BMI), neck circumference and waist to hip
ratio (WHR), marital status, smoking, alcohol consumption and years of education.
The anthropometric data were measured by the unit staff before PSG, while the rest
were self-reported.
In descriptive statistical analysis, continuous variables were expressed as mean ±
standard deviation and categorical in the form of frequencies. Where continuous variables
were not normally distributed the median and interquartile range (IQR) were also given.
Frequency differences between categorical variables were analyzed by the chi-square
test, while the differences in means between continuous variables with analysis of
variance (ANOVA) or Kruskal-Wallis test, depending on the normality of data. For post
hoc analysis we applied the Bonferroni adjustment for pairwise comparisons to account
for type-1 errors.
To examine associations between several variables we used multiple linear regression
models with the backward stepwise method. Non-normally distributed dependent variables
were transformed to their square roots before performing linear regression, while
categorical independent variables were entered in the regression model in the form
of dummy variables. We also implemented ordinal regression to examine the probabilities
of independent variables to predict a higher severity category of AHI, applying the
usual cut-offs of 15 and 30 events/h to distinguish between mild, moderate and severe
OSA[12].
Missing data were handled with multiple imputation, assuming they were missing at
random, using the fully conditional specification method. Suitable transformations
were used to deal with non-normality and then complete data were transformed back
to their original scales before analysis. Five sets of data were created and the results
of linear or ordinal regression in each set were pooled to report an average estimate.
All significance tests were two-sided, with p-value < 0.05 being considered statistically significant. Analysis was conducted with
the statistical software package IBM(®) SPSS(®) Statistics Version 20.
RESULTS
Among 399 patients who underwent PSG during the study period, 53 were excluded for
not fulfilling the inclusion criteria and 36 for not having an OSA diagnosis. From
the rest, we could not codify 28 patients in neither of the 3 socio-economic classes
(housewives=20, students=2, long term unemployed=1, insufficient information=5), limiting
the final sample to 282 patients. 241 were males (85.5%) and 41 females (14.5%), while
the mean age of the sample was 54.61±12.17. The distribution of subjects in socio-economic
classes and the demographic, anthropometric and social characteristics of the sample
per class are presented in [Table 1]. Subjects of the intermediate class had significantly larger neck circumference
than those of the upper class (p=0.022). As expected, salariat and working class had more and less years of education
respectively than intermediate class (p<0.001). There were no significant differences between the 3 classes in terms of symptoms
and associated diseases, except from the fact that individuals of the upper class
had significantly more complaints of non-restorative sleep than the other groups (chi-square
10.561,p=0.005).
Table 1
Descriptive statistics of study population per socio-economic class with missing data
frequencies.
Variable
|
EseC 1 (Ν=99)
|
EseC 2 (Ν=70)
|
EseC 3 (Ν=113)
|
p
|
Sex
|
|
|
|
|
Male
|
84 (84.8)
|
63 (90)
|
94 (83.2)
|
0.436
|
Female
|
15 (15.2)
|
7 (10)
|
19 (16.8)
|
|
Age (years)
|
54.44±12.33
|
55.43±12.05
|
54.25±12.18
|
0.806
|
BMI (kg/m2)
|
33.43±5.93
|
35.41±6.88
|
34.68±6.62
|
0.127
|
Missing data
|
0 (0)
|
1 (1.4)
|
0 (0)
|
0.219
|
Neck circumference (cm)
|
42.26±3.48
|
43.61±3.59
|
42.95±3.67
|
0.025*
|
Median (IQR)
|
42 (40-44)
|
43 (41-47)
|
43 (41-46)
|
|
Missing data
|
0 (0)
|
1 (1.4)
|
2 (1.8)
|
0.430
|
WHR
|
1±0.1
|
1.01±0.05
|
1.01±0.08
|
0.376
|
Missing data
|
0 (0)
|
1 (1.4)
|
2 (1.8)
|
0.430
|
Marital status
|
|
|
|
|
Married
|
75 (75.8)
|
54 (77.1)
|
97 (85.8)
|
0.336
|
Single
|
21 (21.2)
|
13 (18.6)
|
16 (14.2)
|
|
Missing data
|
3 (3)
|
3 (4.3)
|
0 (0)
|
0.110
|
Smoking
|
|
|
|
|
Smoker
|
34 (34.3)
|
28 (40)
|
49 (43.4)
|
0.314
|
Former smoker
|
32 (32.3)
|
28 (40)
|
35 (31)
|
|
Non-smoker
|
32 (32.3)
|
14 (20)
|
28 (24.8)
|
|
Missing data
|
1 (1)
|
0 (0)
|
1 (0.9)
|
0.713
|
Alcohol consumption
|
|
|
|
|
Daily
|
4 (4)
|
11 (15.7)
|
13 (11.5)
|
0.053
|
Occasionally
|
45 (45.5)
|
20 (28.6)
|
43 (38.1)
|
|
Almost never
|
47 (47.5)
|
36 (51.4)
|
55 (48.7)
|
|
Missing data
|
3 (3)
|
3 (4.3)
|
2 (1.8)
|
0.603
|
Years of education
|
|
|
|
|
<7
|
0 (0)
|
3 (4.3)
|
18 (15.9)
|
<0.001
|
7-12
|
5 (5.1)
|
20 (28.6)
|
54 (47.8)
|
|
>12
|
88 (88.9)
|
39 (55.7)
|
23 (20.4)
|
|
Missing data
|
6 (6.1)
|
8 (11.4)
|
18 (15.9)
|
0.078
|
BMI=body mass index; WHR=waist to hip ratio; IQR=interquartile range
*statistical significant pairwise comparison only between EsEC 1 and EsEC 2 (p=0.022)
Categorical and missing data are presented as frequency count (%) and analyzed with
chi-square test
Continuous data are presented as mean ± standard deviation and analyzed with one-way
ANOVA, except from neck circumference presented also as median (IQR) and analyzed
with Kruskal-Wallis test
The mean AHI of the sample was 43.32±27.1 and 59.9% of the subjects had severe OSA,
based on an AHI ≥30. The mean ESS score was 9.16±5.16. There was a statistically significant
positive, but weak, correlation between AHI and ESS score (Spearman’s rho 0.236, p<0.001). The mean and median AHI and ESS scores and the frequencies of the severity
categories across socio-economic classes are shown in [Table 2]. Intermediate class was found to have the highest AHI, although the overall difference
between classes was marginally non-significant (p=0.059). In post-hoc analysis the above trend was only observed between higher and
intermediate class (p=0.075).
Table 2
Mean AHI and ESS scores and prevalence of AHI severity categories of study population
per socio-economic class with missing data frequencies.
Outcome variable
|
ESeC 1
|
ESeC 2
|
ESeC 3
|
p
|
ESS score
|
9.24±5.27
|
9.41±4.9
|
8.92±5.2
|
0.645
|
Median (IQR)
|
8.5 (5-12)
|
9 (6-13)
|
8 (5-12)
|
|
Missing data
|
1 (1)
|
2 (2.9)
|
4 (3.5)
|
0.485
|
AHI
|
39.81±27.01
|
50.51±30.43
|
41.94±24.26
|
0.059
|
Median (IQR)
|
31.7 (16.7-63.5)
|
45 (24-78.6)
|
39.3 (22-64.4)
|
|
AHI groups
|
|
|
|
|
<15
|
23 (23.2)
|
6 (8.6)
|
18 (15.9)
|
0.133
|
15-29.9
|
23 (23.2)
|
19 (27.1)
|
24 (21.2)
|
|
≥30
|
53 (53.5)
|
45 (64.3)
|
71 (62.8)
|
|
ESS=Epworth sleepiness scale; AHI=apnea-hypopnea index; IQR=interquartile range
Categorical and missing data are presented as frequency count (%) and analyzed with
chi-square test
Continuous data are presented as mean ± standard deviation, median (IQR) and analyzed
with Kruskal-Wallis test
A multiple linear regression analysis was performed to reveal statistically significant
predictors of the variance in ESS score and AHI ([Table 3]). To achieve normal distribution, both variables were transformed to their square
roots before analysis. Because of the linear association between socio-economic status
and years of education, the latter variable was excluded from analysis to avoid multicollinearity.
Table 3
Pooled estimated results of multiple linear regression analysis with the backward
stepwise method for the dependent variables √ESS and √AHI after multiple imputation
of missing data.
Independent variable
|
B Coefficient (SE)
|
p
|
|
Dependent variable: √ESS
|
|
Age (per year increase)
|
-0.01 (0.004)
|
0.020
|
Neck circumference (per cm increase)
|
0.06 (0.016)
|
<0.001
|
WHR (per unit increase)
|
-1.33 (0.718)
|
0.065
|
Alcohol consumption (vs almost never)
|
|
|
Daily
|
-0.34 (0.174)
|
0.051
|
Occasionally
|
-0.26 (0.111)
|
0.022
|
|
Dependent variable: √AHI
|
|
BMI (per unit increase)
|
0.06 (0.021)
|
0.004
|
Neck circumference (per cm increase)
|
0.2 (0.039)
|
<0.001
|
Marital status (married vs single)
|
0.74 (0.299)
|
0.014
|
Smoking (smoker vs non-smoker)
|
0.51 (0.228)
|
0.025
|
Socio-economic class (intermediate vs salariat)
|
0.45 (0.257)
|
0.082
|
ESS=Epworth sleepiness scale; AHI=apnea-hypopnea index; BMI=body mass index; WHR=waist
to hip ratio
Variables entered: sex, age, BMI, neck circumference, WHR, marital status, smoking
(two dummy variables), alcohol consumption (two dummy variables), socio-economic class
(two dummy variables)
Criterion for stepwise variable removal: p>0.1
For √ESS, age and occasional alcohol consumption were significant negative correlates
and neck circumference was significant positive correlate, explaining the variance
in ESS score by 10%. Respectively, BMI, neck circumference, being married and current
smoker were independent positive correlates for √AHI, the model accounting for 28%
of the variance. SES was not an independent predictor in both models (p>0.05), using salariat as the reference category. In the √AHI model, however, the
significance of the difference between higher and intermediate class remained over
the 90% confidence level (p=0.082).
We further performed an ordinal logistic regression analysis in order to reveal significant
predictors of OSA severity according to AHI category ([Table 4]). Since the frequency counts constantly increase from the lowest to the highest
severity category in our sample, we used the complementary log-log link function to
transform the cumulative probabilities. In this model, only neck circumference was
statistically significant (p=0.006), predicting a 14% increase in the odds of being in a higher severity category
for every 1 cm increase.
Table 4
Pooled estimated results of ordinal logistic regression analysis with the complementary
log-log link function for the AHI dependent variable ordered by severity category
after multiple imputation of missing data.
Independent variable
|
OR (CI 95%)
|
p
|
Sex
|
|
|
Male
|
1.54 (0.75-3.20)
|
0.243
|
Female
|
1*
|
|
Age (per year increase)
|
1.01 (0.99-1.03)
|
0.256
|
BMI (per unit increase)
|
1.04 (0.99-1.09)
|
0.089
|
Neck circumference (per cm increase)
|
1.14 (1.04-1.24)
|
0.006
|
WHR (per unit increase)
|
0.36 (0.03-4.63)
|
0.434
|
Marital status
|
|
|
Married
|
1.41 (0.88-2.27)
|
0.150
|
Single
|
1*
|
|
Smoking
|
|
|
Smoker
|
1.34 (0.81-2.21)
|
0.252
|
Former smoker
|
0.86 (0.53-1.40)
|
0.547
|
Non-smoker
|
1*
|
|
Alcohol consumption
|
|
|
Daily
|
0.79 (0.40-1.53)
|
0.483
|
Occasionally
|
0.92 (0.58-1.46)
|
0.714
|
Almost never
|
1*
|
|
Socio-economic class
|
|
|
Salariat
|
1*
|
|
Intermediate
|
1.15 (0.69-1.93)
|
0.588
|
Working class
|
1.16 (0.75-1.79)
|
0.500
|
BMI=body mass index; WHR=waist to hip ratio
*Reference category.
DISCUSSION
Our initial hypothesis that OSA patients from low socio-economic classes would present
with greater severity of respiratory events during sleep and subjective daytime sleepiness
was not confirmed. There seems to be a trend for higher AHI in intermediate class,
although it did not reach statistical significance neither in univariate nor in multivariate
analysis. Moreover, intermediate class in our sample had significant higher neck circumference
than upper class and that difference could partially be responsible for the observed
trend. The social factors that were found to independently predict the variance in
OSA severity were marital status and smoking for AHI variance and alcohol consumption
for ESS score variance.
The relationship between low SES and obesity[15], increased consumption of tobacco[16] and alcohol[17], consistently found in literature, was not observed in our sample. In fact, our
data show that intermediate class have more similarities in obesity scores and patterns
of social habits with working class than salariat, implying that they could also share
similar health risks. The last decade’s economic crisis in Greece has mostly affected
the urban middle class, shrinking its income and widening the gap between the wealthier
upper classes and the lower ones, resulting in lower self-rated health[18] and rising unmet needs for health care[19].
We did not find a statistically significant effect of SES on OSA severity, assessed
by ESS score and AHI, after applying multiple regression models to control for potential
confounders, such as age, sex, body habitus measurements and social factors. Ramsey
et al. studied 4042 OSA patients and comparing income categories in terms of AHI found
also no significant differences after adjustment for BMI in both sexes[6].
However, the trend for higher AHI in subjects of the intermediate class compared with
upper class, which was included in our final multiple regression model for AHI variance,
even with lower level of statistical significance (p<0.1), cannot be attributed solely to anthropometric differences. Possible explanations
are differences in referral patterns between classes, since upper class often has
better access to health care, or it could reflect differences between distinct occupations
in each class. The majority of subjects of the intermediate class in our sample were
office clerks (52.9%), being at most a sedentary occupation. In recent literature,
light activity or sedentary occupations have been associated with increased risk for
moderate to severe OSA[20].
Despite non-significant differences in AHI and ESS score between classes, patients
of the upper class complained significantly more for not obtaining restorative sleep
most of the nights than the other classes. It is possible that, since they are more
educated and thereby more cultured, they would recognize easier their symptoms and
their day-to-day variability.Results of a large US cross-sectional epidemiologic survey
also showed that individuals with the lowest educational attainment, particularly
immigrants, reported fewer sleep symptoms than the more educated groups or the native
born[21]. In a similar Brazilian survey, subjects with higher family income were more likely
to report the presence of any sleep complaint, as well as insufficient sleep, snoring
and bruxism, while those with lower income complained more about insomnia and superficial
sleep[22].
Older age and social alcohol drinking were protective factors for subjective daytime
sleepiness in our multiple regression model. Previous studies have reproduced the
same findings, using both subjective and objective measurements of EDS. Bixler et
al.[23] examined a large Pennsylvanian cohort from the general population and observed that
increasing age was associated with less subjective EDS, suggesting the presence of
unsatisfied sleep needs and depression in the young. Budhiraja et al.[24] recently showed that ESS score decreased and mean sleep latency in maintenance of
wakefulness test increased with advancing age in a multicenter OSA cohort, giving
the explanation of disrupted sleep homeostatic mechanisms with ageing. However, since
elderly individuals often consider their sleepiness normal and EDS was found to have
no impact on quality of life of elderly OSA patients[25], it is also possible that they seek less frequently medical assistance than younger
sleepy OSA patients.
Regarding alcohol consumption, Pack et al.[26] found in a sample of older adults that alcohol use reduced the risk for subjective
EDS, hypothesizing that awareness of the negative effect of alcohol on sleep gradually
leads to a decrease in its consumption. A similar result was obtained from a population
survey in US[27]; however, the authors discovered that the interaction between heavy alcohol drinking
and decreased sleep duration predicted increased EDS and considered sleep duration
to be a confounding factor. Despite the objectively evaluated detrimental effects
of alcohol consumption on sleep and daytime alertness in multiple studies, alcohol
users may still perceive its impact as beneficial and rate it accordingly, perhaps
due to differential expectations[28]. Further research using both subjective and objective measurements of EDS is required
to test this assumption in OSA patients. The fact that these risk factors account
for only a small percentage of ESS score variance in our sample highlights the multifactorial
nature of EDS.
Obesity, large neck circumference and smoking were independent risk factors for higher
AHI in our study, results consistent with previous research. Peppard et al.[29] showed that a 10% weight gain in 4 years predicted 32% increase in AHI. Neck circumference
has been recognized as better predictor of OSA severity than visceral obesity, especially
in non-obese patients[30], and smoking has been associated with upper airway inflammation and narrowing, worsening
OSA[31]. The finding that married patients had significantly higher AHI than singles was
also observed in another clinical-based study but not in community samples[32]. Since referral patterns between married and unmarried patients can substantially
differ, depending on the presence of a bed partner who witnesses the relevant symptoms
and behaviors, this relationship can be subjected to selection bias rather than represent
a true association. Moreover, the finding from our ordinal regression model that neck
circumference was the only significant correlate of the probability of being in a
higher severity category in terms of AHI implies that, unlike the social factors examined,
it could serve as a useful predictor in clinical practice, being able to identify
the most severe OSA cases.
Our study has several limitations. Because of the retrospective nature of our data,
classification in socio-economic classes was based in a single open-type question
about subjects’ most recent occupational activity. As a result, the entirety of description
varied greatly between individuals, allowing us to codify some of them in hierarchically
less detailed occupational group (major or sub-major) in ISCO-88. It is, however,
possible that this simplification could in some cases overestimate or underestimate
the positioning in socio-economic class in relation to subjects’ actual occupation.
In the same manner, there were no information about length of employment and former
occupational activities.
Since subjects in this study had not emerged from sampling of the general population,
where patients with less severe OSA are more likely to exist, caution must be taken
when attempting to generalise our results to the whole referent population. Furthermore,
the cross-sectional design lacks definite power in finding causative associations
between outcome and exposure, because they were assessed at the same time.
In conclusion, we have shown that SES has a minor effect on OSA severity. Intermediate
class patients tend to have worse OSA than upper class, although differences in certain
obesity indices were also noted. Further research with prospective studies is required
to test the effect of SES on OSA presence and severity. Already known risk factors,
such as obesity, large neck circumference and smoking, were found independent predictors
of severity of respiratory events at sleep, while the role of alcohol consumption
and marital status on OSA severity needs further clarification in future research.