Next Article in Journal
Implementing a Diabetes Education Program to Reduce Health Disparities in South Texas: Application of the RE-AIM Framework for Planning and Evaluation
Next Article in Special Issue
Endothelin-1 and LOX-1 as Markers of Endothelial Dysfunction in Obstructive Sleep Apnea Patients
Previous Article in Journal
Poor Oral Hygiene and High Levels of Inflammatory Cytokines in Saliva Predict the Risk of Overweight and Obesity
Previous Article in Special Issue
Comparison of Findings between Clinical Examinations and Drug-Induced Sleep Endoscopy in Patients with Obstructive Sleep Apnea Syndrome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Association between Temporomandibular Disorder and Sleep Apnea—A Nationwide Population-Based Cohort Study

1
Department of Dentistry, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
2
Department of Oral Hygiene, College of Dental Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
3
Division of Pulmonary and Critical Care Medicine, Department of Internal Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
4
Sleep Disorders Center, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
5
School of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
6
Department of Respiratory Care, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
7
Graduate Institute of Clinical Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
8
Department of Neurology, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(17), 6311; https://doi.org/10.3390/ijerph17176311
Submission received: 2 August 2020 / Revised: 24 August 2020 / Accepted: 28 August 2020 / Published: 30 August 2020
(This article belongs to the Special Issue Obstructive Sleep Apnea Syndrome: From Symptoms to Treatment)

Abstract

:
An increased incidence of temporomandibular disorders (TMD) among patients with sleep apnea (SA) has been reported. However, the association between TMD and SA has not been demonstrated in a large-scale study. This population-based cohort study with the Taiwan National Health Insurance (NHI) Research Database aimed to understand the association between SA and TMD. We identified adult patients with suspected SA (identified with diagnostic codes) and excluded those diagnosed with TMD prior to SA. Patients with SA diagnosis after polysomnography were also identified as probable SA patients. The index dates were the dates of their initial SA diagnosis. Ten control subjects were matched, by age and sex, to each SA patient, and were assigned the same index dates as the SA patients. In total, 10,408 suspected SA patients (including 4105 probable SA patients) matched to 104,080 control subjects (including 41,050 subjects matched to the probable SA patients) in this study. The TMD incidence rate was significantly higher in the SA patients than in the control subjects (2.8 vs. 1.0 per thousand-patient-year in probable SA patients vs. the corresponding control subjects, with an adjusted incidence rate ratio [95% confidence interval] = 2.5 [2.3–2.7], p < 0.0001). SA patients significantly showed a higher cumulative incidence of TMD than the corresponding control subjects (p < 0.0001). Multivariable Cox regression analysis revealed SA as an independent risk factor for the development of TMD (adjusted hazard ratio = 2.5 [1.7–3.7], p < 0.0001). In summary, this study confirmed an increased TMD incidence in the SA patients. While treating TMD patients, dentists should pay careful attention to the potential underlying SA.

1. Introduction

Temporomandibular disorders (TMDs) are musculoskeletal disorders characterized by persistent pain in the temporomandibular joint, masticatory muscles, and periauricular region [1]. Based on clinical evaluation, the prevalence of TMD was 6–93% in the general population [2]. It has been considered a complex and multifaceted disease process [3]. Biomechanical, neuromuscular, neurobiological, and biopsychosocial factors may contribute to the disorder [4]. Risk factors contributing to TMD includes age, genetic factors, sex, stress, anxiety, occlusion, poor posture, rheumatoid arthritis, and dysfunctional breathing [2,5,6]. TMD develops at a markedly higher rate in individuals with relatively poorer health status, whether in the form of other pain conditions, comorbid diseases, poor sleep quality, or cigarette smoking [7].
Being the most common type of sleep-disordered breathing, sleep apnea (SA) is usually characterized by repetitive shallowing and stopping of breaths during sleep, which is usually accompanied with intermittent hypoxia, arousal, and sleep fragmentation [8]. Polysomnography (PSG) is usually performed to establish the diagnosis of SA [9]. Most (>90%) cases are obstructive sleep apnea (OSA), while a small proportion of cases are central sleep apnea (CSA) [10,11]. SA has been associated with a variety of diseases, such as arrhythmia, hypertension, ischemic heart disease, stroke, metabolic syndrome, diabetes mellitus, chronic kidney disease, erectile dysfunction, hormonal dysfunction, and dementia [8,12,13,14].
Several previous studies have shown the association between sleep disorders and TMD [6,7,15], while only few of them discussed the association between SA and TMD. In the OPPERA (Orofacial Pain: Prospective Evaluation and Risk Assessment) cohort, which recruited subjects from 2006 to 2008 and followed them through 2012, signs/symptoms of OSA were associated with higher risk of incident TMD [6]. Because both SA and TMD are associated with several comorbidities, the association between SA and TMD might be very complex. A large-scale cohort study that takes comorbidities into the analyses is therefore needed. This large population-based cohort study aimed to confirm the association between SA and TMD by demonstrating a higher TMD incidence in SA patients.

2. Materials and Methods

2.1. Study Design

This study is a retrospective cohort study using a large nationwide population-based database.

2.2. Data Sources

Since March 1995, the Taiwan National Health Insurance (NHI) has covered outpatient care, dental services, hospital inpatient care, and drug prescriptions. In the population of about 23 million people, the coverage rate reached 99% in 2004 [8,16]. From the NHI Research Database (NHIRD) managed by the National Health Research Institutes in Taiwan, we obtained the Longitudinal Health Insurance Database 2005 (LHID2005), a cohort of one million randomly sampled subjects in the NHI system in 2005, for this study. For confidentiality, patient identification was encrypted in the database.

2.3. Ethical Consideration

The Institutional Review Board of Kaohsiung Medical University Hospital approved this study (KMUH-IRB-EXEMPT-20130034).

2.4. Sleep Apnea Cohorts

Patients with “suspected SA” were identified using an algorithm (Figure 1). Patients that received a diagnostic code of SA between 1995/03/01 and 2013/12/31 were identified initially. The International Classification of Diseases, 9th Revision, Clinical Modification (ICD-9-CM) codes of 780.51, 780.53, and 780.57 were used for the diagnosis of SA [8,17]. The ICD-9-CM codes of 524.6 was used to identify the diagnosis of TMD [18]. Each patient was followed from the date of first SA diagnosis (the index date). The exclusion criteria included: patients with insufficient (<1 year) washout periods (between NHI enrollment and the index date); patients with insufficient (<1 year) follow-up period; patients that received a TMD diagnosis before the index date; aged less than 18 or more than 90 years. Therefore, the “suspected SA cohort” consisted of SA patients identified with diagnostic codes from the database.
In the suspected SA patients, those who had retained their SA diagnosis after a PSG examination were further extracted and were defined as the “probable SA cohort.” Theoretically, these “probable SA patients” might have PSG-confirmed SA.

2.5. Control Cohorts

We randomly selected 10 age- and sex-matched control subjects for each suspected SA patient. The index dates of control subjects were assigned based on those of their corresponding SA patients. We applied the same exclusion criteria while selecting the control subjects in the matching process.

2.6. Definitions of Variables

The development of TMD was taken as the endpoint in the current study. Only those whose TMD diagnosis appeared at least three times in outpatient claims or at least one time in inpatient claims were considered TMD patients, to increase the accuracy of the diagnosis. A comorbidity was identified by the first appearance of its corresponding diagnostic codes (before the index date), and was confirmed by the appearance in at least three outpatient claims or an inpatient claim. We calculated the Charlson Comorbidity Index (CCI) score of each subject [19].
Each subject was followed from the index date to either the development of TMD, termination of the record due to death or withdrawal from NHI, or the end of 2013 (the end of study period), whichever came first.

2.7. Statistical Analysis

The current study comprised two study arms with the same statistical analyses: study arm A comparing suspected SA cohort and control A cohort; study arm B comparing probable SA cohort and control B cohort. Firstly, the demographic information was compared between SA cohorts and the corresponding control cohorts with Student’s t-test and Pearson χ2 test for continuous and categorical variables, respectively. The incidence rate of TMD was the number of TMDs occurred in the follow-up period divided by the total person-year. By calculating the incidence rate ratio (IRR), the incidence rates of TMD in the SA cohorts and the corresponding control cohorts were compared. With the assumption that the occurrence of TMD followed a Poisson probability distribution, the 95% confidence interval (CIs) of the IRRs were estimated. We also performed stratified analyses by classifying the subjects based on sex, age, residency, income level, or presence of a comorbidity. The cumulative incidence of TMD was estimated with the Kaplan–Meier method and compared with log-rank test. In order to further assess the effect of SA, multivariable Cox proportional hazards regression analyses, adjusting for sex, age, residency, income level, and comorbidities, were performed. Multivariable Cox regression analyses, built by a backward variable selection method, eliminating variables with a p-value of less than 0.05, were used to identify the risk factors for developing TMD in SA patients. The adjusted IRRs and hazard ratios (HRs) were adjusted for sex, age group, residency, income level, and the presence of comorbidities (except for the stratifying variable).
We used the Statistical Analysis System (SAS) system (version 9.4 for Windows, SAS Institute Inc., Cary, NC, USA) for extraction and computation of data, data linkage, processing and sampling, and statistical analysis. We set the statistical significance level at a two-sided p value of less than 0.05.

3. Results

A total of 10,408 suspected SA patients, including 4105 probable SA patients, were identified and were matched to control subjects (104,080 and 41,050 subjects in the control A and control B cohorts, respectively) by age and sex using the algorithm (Figure 1). Table 1 presents the baseline characteristics of the cohorts. Patients in the SA cohorts had more comorbidities than the subjects in the corresponding control cohorts (Table 1).
The suspected SA patients had a significantly higher TMD incidence rate than the subjects in the control A cohort (2.9 vs. 1.1 per thousand-patient-year, adjusted IRR [95% CI] = 2.0 [1.9–2.1], p < 0.0001). A significantly higher TMD incidence rate was also observed in patients with probable SA than in the subjects in the control B cohort (2.8 vs. 1.0 per thousand-patient-year, adjusted IRR [95% CI] = 2.5 [2.3–2.7], p < 0.0001) (Table 2). Stratified analyses showed that SA patients had significantly higher TMD incidence rates than the corresponding control subjects in all strata (Table 2).
The suspected SA cohort had a significantly higher cumulative TMD incidence than the control A cohort (p < 0.0001); the probable SA cohort also had a significantly higher cumulative TMD incidence than the control B cohort (p < 0.0001) (Figure 2). Consistent results were noted in stratified analyses except for female subjects in study arm B (Figure A1 in Appendix A).
Multivariable Cox proportional hazards regression analyses found that SA was an independent risk factor for the development of TMD (adjusted HR [95% CI] = 2.0 [1.6–2.6] and 2.5 [1.7–3.7] in study arm A and B, respectively; both p < 0.0001) (Figure 3 and Table A1 in Appendix A). Stratified analyses showed SA as an independent risk factor for incident TMD in nearly all strata (Figure 3).
Using multivariable Cox regression models, we tried to find the risk factors for developing TMD in SA patients. In suspected SA patients, the risk factors for incident TMD included peptic ulcer disease (adjusted HR [95% CI] = 2.0 [1.3–3.2], p = 0.0017) and cancer (adjusted HR [95% CI] = 2.4 [1.3–4.4], p = 0.0065) (Table 3). In probable SA patients, peptic ulcer disease was found as the only risk factor for incident TMD (adjusted HR [95% CI] = 3.1 [1.6–6.1], p = 0.0010).

4. Discussion

SA patients had a significantly higher TMD incidence than the control subjects in this large nationwide, population-based cohort study. The analyses in the study arm A, which included SA patients identified by only using the ICD-9-CM codes (i.e., suspected SA cohort), showed similar results as the analyses in the study arm B, which included SA patients diagnosed with SA after PSG (i.e., probable SA cohort). After being adjusted for sex, age, and comorbidities, SA remained an independent risk factor for the development of TMD.
This study is the first and largest long-term, population-based cohort study to survey TMD incidence in SA patients. In line with our finding, several previous studies have shown an association between SA and TMD [20]. Indeed, disturbed sleep may interfere with patients’ daily function, contributing to an increase in their sensitivity to pain [21]. Sleep disturbance may increase pain sensitivity in patients with chronic pain and may create a self-perpetuating cycle of sleep disruption, more pain, and depression. A bidirectional association between OSA and TMD has been suggested, which is evident in their prevalence rates: an increased prevalence of TMD in OSA patients [22] and an increased prevalence of OSA in TMD patients [15]. OSA has been reported as a common comorbidity in general chronic pain populations, with a pooled prevalence of about 37%, as well as in TMD patients, with a frequency of 28.6% [15]. Women with TMD present more respiratory effort-related arousals than healthy controls [23], predisposing them to the development of OSA.
As this study showed that SA patients had a higher risk of subsequent TMD, it is suggested that physicians treating SA patients should pay careful attention to the possible associated chronic pain, which affects the quality of sleep. Moreover, as the current study revealed that SA was an independent risk factor for the development of TMD, dentists must pay careful attention to the possible underlying SA while treating TMD patients so that appropriate referrals may be made to improve treatment outcomes.
Several mechanisms might be underlying the close relationship between OSA and TMD [24]. Firstly, inadequate and/or disrupted sleep in OSA patients might enhance pain sensitivity, contributing to hyperalgesia, an important feature found in many TMD patients [25]. Secondly, OSA patients usually display chronic intermittent hypoxemia, which increases the levels of inflammatory cytokines, contributing to the pathogenesis of multiple comorbidities [26]. As higher plasma levels of inflammatory cytokines were reported both in OSA and TMD [27,28], we postulated that OSA might contribute to the pathogenesis of TMD through enhancing systemic inflammation. It is therefore not surprising that SA patients have higher TMD incidence. Thirdly, malocclusion (misalignment of the mandible to the cranium) may prevent the airway from staying open during sleep, thus, positional therapy and intraoral protrusion devices are sometimes suggested to treat OSA [29]. However, using mandibular-advancement oral appliances as a treatment modality for OSA may also cause TMD [29,30,31,32,33]. Fourthly, TMD and sleep bruxism may concomitantly present in OSA patients. The bruxism episode index (BEI) positively correlated with AHI in patients with mild-to moderate OSA, whereas patients with severe OSA had lower BEI than those with mild-to-moderate OSA [34]. Fifthly, SA patients, mostly OSA, may have craniofacial configurations and/or muscle dysfunction that predispose them to the development of TMD.
A variety of risk factors for TMD have also been reported in the literature. Traditionally, female sex has been suggested as a risk factor for the development of TMD. For example, Contreras et al. found a higher TMD prevalence in women [35]. Sexual hormones, especially estrogen, might play an important role in the sensitivity of pain, even in the muscles of mastication, and therefore, might contribute to the pathogenesis of TMD [36]. A few studies also showed that the tolerance to pain varied in different phases during the menstrual cycle [36]. In the present study, the TMD prevalence was higher in women than in men in both study arm A and B. The multivariable Cox regression analyses models (Table A1 in the Appendix A) also suggested male sex as an independent protective factor (HR [95% CI] = 0.6 [0.5–0.7] and 0.5 [0.4–0.7] in study arm A and B, respectively). In addition, reduced 25-hydroxyvitamin D level is associated with decreased functional capability of musculoskeletal system [37]. Vitamin D deficiency might contribute to the development of TMD through increasing the parathyroid hormone [38]. Since vitamin D deficiency is more prevalent in women, and is usually associated with poorer musculoskeletal health in community-living people in Taiwan [39], it is not surprising that female sex might be a risk factor for TMD. In a previous study by Contreras [35], peptic ulcer disease has been associated with the presence of painful TMD in the univariate analysis, whereas this association was not shown in the multivariable model. In our current study, we found peptic ulcer disease to be an independent risk factor for incident TMD in SA patients.
Our study had several strengths. Firstly, our study is the largest research discussing the association between SA and TMD thus far. Secondly, a population-based study might minimize the risk of selection bias. Thirdly, this is a cohort study with a very long follow-up time. Fourthly, in contrast to the previous studies with cross-sectional designs, this study demonstrated a temporal trajectory and might therefore provide better evidence regarding the association between SA and TMD.
This study, however, still had a few limitations. Firstly, the diagnoses based on ICD-9-CM codes might be less reliable than those established with standardized diagnostic criteria in clinical trials. However, besides using SA patients identified with only the codes (i.e., suspected SA cohort) for analyses, we also used SA patients diagnosed with SA after PSG (i.e., probable SA cohort) for another set of diagnosis, which found similar results. Secondly, it is impossible to distinguish OSA and CSA with the ICD-9-CM codes. While OSA patients comprise more than 90% of SA patients, we believe that our findings represent the phenomenon in OSA patients. Thirdly, the information of some potential confounders, such as craniofacial structure, masticatory habits, tooth loss, and body mass index, was not available in the database. Matching the case and controls by age and sex might not eliminate the potential confounders. As SA patients had significantly higher prevalence of comorbidities than the control subjects, resulting from multiple factors, multivariable analyses adjusting for the comorbidities might mitigate the bias introduced by the confounders.

5. Conclusions

This large population-based cohort study showed a significantly higher TMD incidence in SA patients.

Author Contributions

Conceptualization, J.-H.W., K.-T.L., C.-Y.K., C.-H.C., J.-Y.C., J.-Y.H., and C.-Y.H.; Data curation, M.-J.T.; Formal analysis, J.-H.W. and M.-J.T.; Funding acquisition, M.-J.T.; Investigation, J.-H.W., K.-T.L., C.-Y.K., C.-H.C., J.-Y.C., and M.-J.T.; Methodology, J.-H.W., K.-T.L., C.-Y.K., C.-H.C., J.-Y.H., C.-Y.H., and M.-J.T.; Supervision, J.-Y.H., C.-Y.H., and M.-J.T.; Writing—original draft, J.-H.W. and M.-J.T.; Writing—review and editing, J.-H.W., K.-T.L., C.-Y.K., C.-H.C., J.-Y.C., J.-Y.H., C.-Y.H., and M.-J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from Kaohsiung Medical University Hospital [KMUH108–8R13; KMUH108-8M63].

Conflicts of Interest

All authors declare no conflict of interest.

Abbreviations

CCICharlson Comorbidity Index
CIconfidence interval
CSAcentral sleep apnea
HRhazard ratio
IRincident rate
IRRincidence rate ratio
NHINational Health Insurance
OSAobstructive sleep apnea
PSGpolysomnography
SAsleep apnea
TMDtemporomandibular disorders

Appendix A

Figure A1. The cumulative incidences of temporomandibular disorder in female subjects (a,b), male subjects (c,d), subjects ≤ 50 years of age (e,f), and subjects > 50 years of age (g,h). The analyses were taken in study arm A (suspected SA vs. control A) (a,c,e,g) and study arm B (probable SA vs. control B) (b,d,f,h). The red continuous lines and blue dashed lines show the cumulative incidence of temporomandibular disorder for the sleep apnea (SA) patients and the control subjects, respectively.
Figure A1. The cumulative incidences of temporomandibular disorder in female subjects (a,b), male subjects (c,d), subjects ≤ 50 years of age (e,f), and subjects > 50 years of age (g,h). The analyses were taken in study arm A (suspected SA vs. control A) (a,c,e,g) and study arm B (probable SA vs. control B) (b,d,f,h). The red continuous lines and blue dashed lines show the cumulative incidence of temporomandibular disorder for the sleep apnea (SA) patients and the control subjects, respectively.
Ijerph 17 06311 g0a1
Table A1. Multivariable Cox analysis of the contributing factors for the development of temporomandibular disorder.
Table A1. Multivariable Cox analysis of the contributing factors for the development of temporomandibular disorder.
VariablesStudy Arm A
(Suspected SA vs. Control A)
Study Arm B
(Probable SA vs. Control B)
HR [95% CI]p-ValueHR [95% CI]p-Value
SA patients vs. Control subjects2.0 [1.6–2.6]<0.00012.5 [1.7–3.7]<0.0001
Male vs. female0.6 [0.5–0.7]<0.00010.5 [0.4–0.7]0.0003
Age > 50 vs. Age ≤ 500.9 [0.8–1.2]0.60551.3 [0.9–1.8]0.1564
Residency
 (Northern Taiwan vs. Other areas)
1.0 [0.9–1.3]0.67751.2 [0.9–1.7]0.2368
Higher income (>NT$24,000)
 vs. lower income (≤NT$24,000)
1.2 [1.0–1.4]0.12761.4 [1.0–1.9]0.0611
Presence of underlying diseases:
 Heart disease0.6 [0.3–1.1]0.08400.8 [0.3–2.1]0.6785
 Peripheral vascular disease1.3 [0.6–2.6]0.52972.4 [1.0–5.5]0.0448
 Major neurological disorder1.1 [0.8–1.5]0.64641.4 [0.9–2.4]0.1701
 Chronic pulmonary disease1.2 [0.9–1.5]0.13041.1 [0.7–1.6]0.7958
 Connective tissue disease1.6 [1.0–2.6]0.06091.8 [0.8–3.8]0.1538
 Peptic ulcer disease1.7 [1.4–2.1]<0.00011.7 [1.2–2.5]0.0032
 Liver disease1.4 [1.1–1.7]0.00841.2 [0.8–1.7]0.4586
 Diabetes mellitus1.1 [0.8–1.5]0.53601.2 [0.7–1.9]0.4809
 Renal disease0.3 [0.1–0.7]0.00600.4 [0.1–1.2]0.1035
 Cancer2.0 [1.4–2.8]0.00012.3 [1.3–3.9]0.0031
Abbreviations: SA = sleep apnea; HR = hazard ratio; CI = confidence interval.

References

  1. Schiffman, E.; Ohrbach, R.; Truelove, E.; Look, J.; Anderson, G.; Goulet, J.P.; List, T.; Svensson, P.; Gonzalez, Y.; Lobbezoo, F.; et al. Diagnostic criteria for temporomandibular disorders (dc/tmd) for clinical and research applications: Recommendations of the international rdc/tmd consortium network* and orofacial pain special interest groupdagger. J. Oral Facial Pain Headache 2014, 28, 6–27. [Google Scholar] [CrossRef] [PubMed]
  2. Poveda Roda, R.; Bagan, J.V.; Diaz Fernandez, J.M.; Hernandez Bazan, S.; Jimenez Soriano, Y. Review of temporomandibular joint pathology. Part i: Classification, epidemiology and risk factors. Med. Oral Patol. Oral Cir. Bucal. 2007, 12, 292–298. [Google Scholar]
  3. Liao, C.H.; Chang, C.S.; Chang, S.N.; Lane, H.Y.; Lyu, S.Y.; Morisky, D.E.; Sung, F.C. The risk of temporomandibular disorder in patients with depression: A population-based cohort study. Community Dent. Oral Epidemiol. 2011, 39, 525–531. [Google Scholar] [CrossRef]
  4. Oral, K.; Bal Kucuk, B.; Ebeoglu, B.; Dincer, S. Etiology of temporomandibular disorder pain. Agri 2009, 21, 89–94. [Google Scholar]
  5. Lin, C.Y.; Chung, C.H.; Chu, H.Y.; Chen, L.C.; Tu, K.H.; Tsao, C.H.; Wu, Y.T.; Chien, W.C. Prevalence of temporomandibular disorders in rheumatoid arthritis and associated risk factors: A nationwide study in taiwan. J. Oral Facial Pain Headache 2017, 31, e29–e36. [Google Scholar] [CrossRef] [PubMed]
  6. Sanders, A.E.; Essick, G.K.; Fillingim, R.; Knott, C.; Ohrbach, R.; Greenspan, J.D.; Diatchenko, L.; Maixner, W.; Dubner, R.; Bair, E.; et al. Sleep apnea symptoms and risk of temporomandibular disorder: Oppera cohort. J. Dent. Res. 2013, 92, 70S–77S. [Google Scholar] [CrossRef]
  7. Slade, G.D.; Ohrbach, R.; Greenspan, J.D.; Fillingim, R.B.; Bair, E.; Sanders, A.E.; Dubner, R.; Diatchenko, L.; Meloto, C.B.; Smith, S.; et al. Painful temporomandibular disorder: Decade of discovery from oppera studies. J. Dent. Res. 2016, 95, 1084–1092. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, C.M.; Tsai, M.J.; Wei, P.J.; Su, Y.C.; Yang, C.J.; Wu, M.N.; Hsu, C.Y.; Hwang, S.J.; Chong, I.W.; Huang, M.S. Erectile dysfunction in patients with sleep apnea—A nationwide population-based study. PLoS ONE 2015, 10, e0132510. [Google Scholar] [CrossRef] [Green Version]
  9. Myers, K.A.; Mrkobrada, M.; Simel, D.L. Does this patient have obstructive sleep apnea? The rational clinical examination systematic review. JAMA 2013, 310, 731–741. [Google Scholar] [CrossRef]
  10. Chou, K.T.; Huang, C.C.; Chen, Y.M.; Su, K.C.; Shiao, G.M.; Lee, Y.C.; Chan, W.L.; Leu, H.B. Sleep apnea and risk of deep vein thrombosis: A non-randomized, pair-matched cohort study. Am. J. Med. 2012, 125, 374–380. [Google Scholar] [CrossRef]
  11. Chou, P.S.; Chang, W.C.; Chou, W.P.; Liu, M.E.; Lai, C.L.; Liu, C.K.; Ku, Y.C.; Tsai, S.J.; Chou, Y.H.; Chang, W.P. Increased risk of benign prostate hyperplasia in sleep apnea patients: A nationwide population-based study. PLoS ONE 2014, 9, e93081. [Google Scholar] [CrossRef] [Green Version]
  12. Gronewold, J.; Haensel, R.; Kleinschnitz, C.; Frohnhofen, H.; Hermann, D.M. Sleep-disordered breathing in hospitalized geriatric patients with mild dementia and its association with cognition, emotion and mobility. Int. J. Environ. Res. Public Health 2019, 16, 863. [Google Scholar] [CrossRef] [Green Version]
  13. Arzt, M.; Young, T.; Finn, L.; Skatrud, J.B.; Bradley, T.D. Association of sleep-disordered breathing and the occurrence of stroke. Am. J. Respir. Crit. Care Med. 2005, 172, 1447–1451. [Google Scholar] [CrossRef] [Green Version]
  14. Reichmuth, K.J.; Austin, D.; Skatrud, J.B.; Young, T. Association of sleep apnea and type ii diabetes: A population-based study. Am. J. Respir. Crit. Care Med. 2005, 172, 1590–1595. [Google Scholar] [CrossRef] [Green Version]
  15. Smith, M.T.; Wickwire, E.M.; Grace, E.G.; Edwards, R.R.; Buenaver, L.F.; Peterson, S.; Klick, B.; Haythornthwaite, J.A. Sleep disorders and their association with laboratory pain sensitivity in temporomandibular joint disorder. Sleep 2009, 32, 779–790. [Google Scholar] [CrossRef] [Green Version]
  16. Wu, C.Y.; Chen, Y.J.; Ho, H.J.; Hsu, Y.C.; Kuo, K.N.; Wu, M.S.; Lin, J.T. Association between nucleoside analogues and risk of hepatitis b virus-related hepatocellular carcinoma recurrence following liver resection. JAMA 2012, 308, 1906–1914. [Google Scholar] [CrossRef] [Green Version]
  17. Wu, K.L.; Kuo, C.Y.; Tsai, Y.C.; Hung, J.Y.; Sheu, C.C.; Yang, C.J.; Hsu, C.Y.; Wu, M.N.; Tsai, M.J. CHADS2, CHA2DS2ASc, and new ABCD scores predict the risk of peripheral arterial disease in patients with sleep apnea. J. Clin. Med. 2019, 8, 188. [Google Scholar] [CrossRef] [Green Version]
  18. Lin, S.L.; Wu, S.L.; Ko, S.Y.; Yen, C.Y.; Chiang, W.F.; Yang, J.W. Temporal relationship between dysthymia and temporomandibular disorder: A population-based matched case-control study in taiwan. BMC Oral Health 2017, 17, 50. [Google Scholar] [CrossRef] [Green Version]
  19. Deyo, R.A.; Cherkin, D.C.; Ciol, M.A. Adapting a clinical comorbidity index for use with icd-9-cm administrative databases. J. Clin. Epidemiol. 1992, 45, 613–619. [Google Scholar] [CrossRef]
  20. Collesano, V.; Segu, M.; Masseroli, C.; Manni, R. Temporomandibular disorders and sleep disorders: Which relationship? Minerva Stomatol. 2004, 53, 661–668. [Google Scholar]
  21. Moldofsky, H. Sleep and pain. Sleep Med. Rev. 2001, 5, 385–396. [Google Scholar] [CrossRef] [Green Version]
  22. Cunali, P.A.; Almeida, F.R.; Santos, C.D.; Valdrighi, N.Y.; Nascimento, L.S.; Dal’Fabbro, C.; Tufik, S.; Bittencourt, L.R. Prevalence of temporomandibular disorders in obstructive sleep apnea patients referred for oral appliance therapy. J. Orofac. Pain 2009, 23, 339–344. [Google Scholar]
  23. Dubrovsky, B.; Raphael, K.G.; Lavigne, G.J.; Janal, M.N.; Sirois, D.A.; Wigren, P.E.; Nemelivsky, L.V.; Klausner, J.J.; Krieger, A.C. Polysomnographic investigation of sleep and respiratory parameters in women with temporomandibular pain disorders. J. Clin. Sleep Med. 2014, 10, 195–201. [Google Scholar] [CrossRef] [Green Version]
  24. Talley, R.L. Tmj and osa are sisters. Cranio 2019, 37, 273–274. [Google Scholar] [CrossRef] [Green Version]
  25. Bair, E.; Gaynor, S.; Slade, G.D.; Ohrbach, R.; Fillingim, R.B.; Greenspan, J.D.; Dubner, R.; Smith, S.B.; Diatchenko, L.; Maixner, W. Identification of clusters of individuals relevant to temporomandibular disorders and other chronic pain conditions: The OPPERA study. Pain 2016, 157, 1266–1278. [Google Scholar] [CrossRef]
  26. Dewan, N.A.; Nieto, F.J.; Somers, V.K. Intermittent hypoxemia and osa: Implications for comorbidities. Chest 2015, 147, 266–274. [Google Scholar] [CrossRef] [Green Version]
  27. Park, J.W.; Chung, J.W. Inflammatory cytokines and sleep disturbance in patients with temporomandibular disorders. J. Oral Facial Pain Headache 2016, 30, 27–33. [Google Scholar] [CrossRef]
  28. Doufas, A.G.; Tian, L.; Padrez, K.A.; Suwanprathes, P.; Cardell, J.A.; Maecker, H.T.; Panousis, P. Experimental pain and opioid analgesia in volunteers at high risk for obstructive sleep apnea. PLoS ONE 2013, 8, e54807. [Google Scholar] [CrossRef]
  29. Randerath, W.J.; Verbraecken, J.; Andreas, S.; Bettega, G.; Boudewyns, A.; Hamans, E.; Jalbert, F.; Paoli, J.R.; Sanner, B.; Smith, I.; et al. Non-cpap therapies in obstructive sleep apnoea. Eur. Respir. J. 2011, 37, 1000–1028. [Google Scholar] [CrossRef] [Green Version]
  30. Marklund, M. Subjective versus objective dental side effects from oral sleep apnea appliances. Sleep Breath. 2020, 24, 111–117. [Google Scholar] [CrossRef] [Green Version]
  31. Merrill, R.L. Temporomandibular disorder pain and dental treatment of obstructive sleep apnea. Dent. Clin. North. Am. 2012, 56, 415–431. [Google Scholar] [CrossRef]
  32. Perez, C.V.; de Leeuw, R.; Okeson, J.P.; Carlson, C.R.; Li, H.F.; Bush, H.M.; Falace, D.A. The incidence and prevalence of temporomandibular disorders and posterior open bite in patients receiving mandibular advancement device therapy for obstructive sleep apnea. Sleep Breath. 2013, 17, 323–332. [Google Scholar] [CrossRef]
  33. Doff, M.H.; Veldhuis, S.K.; Hoekema, A.; Slater, J.J.; Wijkstra, P.J.; de Bont, L.G.; Stegenga, B. Long-term oral appliance therapy in obstructive sleep apnea syndrome: A controlled study on temporomandibular side effects. Clin. Oral Investig. 2012, 16, 689–697. [Google Scholar] [CrossRef]
  34. Martynowicz, H.; Gac, P.; Brzecka, A.; Poreba, R.; Wojakowska, A.; Mazur, G.; Smardz, J.; Wieckiewicz, M. The relationship between sleep bruxism and obstructive sleep apnea based on polysomnographic findings. J. Clin. Med. 2019, 8, 1653. [Google Scholar] [CrossRef] [Green Version]
  35. Contreras, E.F.R.; Fernandes, G.; Ongaro, P.C.J.; Campi, L.B.; Goncalves, D.A.G. Systemic diseases and other painful conditions in patients with temporomandibular disorders and migraine. Braz. Oral Res. 2018, 32, e77. [Google Scholar] [CrossRef] [Green Version]
  36. Ferreira, C.L.; Silva, M.A.; Felicio, C.M. Signs and symptoms of temporomandibular disorders in women and men. Codas 2016, 28, 17–21. [Google Scholar] [CrossRef]
  37. Bischoff-Ferrari, H.A.; Dietrich, T.; Orav, E.J.; Hu, F.B.; Zhang, Y.; Karlson, E.W.; Dawson-Hughes, B. Higher 25-hydroxyvitamin d concentrations are associated with better lower-extremity function in both active and inactive persons aged > or = 60 y. Am. J. Clin. Nutr. 2004, 80, 752–758. [Google Scholar] [CrossRef]
  38. Demir, C.Y.; Ersoz, M.E. Biochemical changes associated with temporomandibular disorders. J. Int. Med. Res. 2019, 47, 765–771. [Google Scholar] [CrossRef] [Green Version]
  39. Chen, C.H.; Liu, L.K.; Chen, M.J.; Lee, W.J.; Lin, M.H.; Peng, L.N.; Chen, L.K. Associations between vitamin d deficiency, musculoskeletal health, and cardiometabolic risk among community-living people in taiwan: Age and sex-specific relationship. Medicine 2018, 97, e13886. [Google Scholar] [CrossRef]
Figure 1. Algorithm to identify the study population. Abbreviations: LHID2005 = Longitudinal Health Insurance Database 2005; SA = sleep apnea; TMD = temporomandibular disorder; PSG = polysomnography; y/o = years old.
Figure 1. Algorithm to identify the study population. Abbreviations: LHID2005 = Longitudinal Health Insurance Database 2005; SA = sleep apnea; TMD = temporomandibular disorder; PSG = polysomnography; y/o = years old.
Ijerph 17 06311 g001
Figure 2. The cumulative incidences of temporomandibular disorder. The red continuous lines show the cumulative incidence of temporomandibular disorder for the sleep apnea patients; the blue dashed lines show the cumulative incidence of temporomandibular disorder for the control subjects. (a) Study arm A (the suspected sleep apnea cohort vs. control A cohort); (b) Study arm B (the probable sleep apnea cohort vs. control B cohort).
Figure 2. The cumulative incidences of temporomandibular disorder. The red continuous lines show the cumulative incidence of temporomandibular disorder for the sleep apnea patients; the blue dashed lines show the cumulative incidence of temporomandibular disorder for the control subjects. (a) Study arm A (the suspected sleep apnea cohort vs. control A cohort); (b) Study arm B (the probable sleep apnea cohort vs. control B cohort).
Ijerph 17 06311 g002
Figure 3. Stratified analyses of multivariable Cox analyses. (a) Study arm A (suspected SA vs. control A); (b) Study arm B (probable SA vs. control B). Adjusted hazard ratios (95% confidence intervals) of sleep apnea, adjusted for sex, age, residency, income level, and the presence of various comorbidities (except for the stratifying variable) are presented. Hazard ratio cannot be estimated because of small sample size. Abbreviations: SA = sleep apnea; CCI = Charlson Comorbidity Index.
Figure 3. Stratified analyses of multivariable Cox analyses. (a) Study arm A (suspected SA vs. control A); (b) Study arm B (probable SA vs. control B). Adjusted hazard ratios (95% confidence intervals) of sleep apnea, adjusted for sex, age, residency, income level, and the presence of various comorbidities (except for the stratifying variable) are presented. Hazard ratio cannot be estimated because of small sample size. Abbreviations: SA = sleep apnea; CCI = Charlson Comorbidity Index.
Ijerph 17 06311 g003
Table 1. Baseline characteristics of the study population.
Table 1. Baseline characteristics of the study population.
VariablesStudy Arm AStudy Arm B
Suspected SA Control Ap-ValueProbable SAControl B p-Value
N10,408104,080 410541,050
Sex, n (%)
 Female3823 (37%)38,230 (37%) 894 (22%)8940 (22%)
 Male6585 (63%)65,850 (63%) 3211 (78%)32,110 (78%)
Age (year), mean ± SD47.8 ± 14.947.8 ± 14.9 47.6 ± 13.347.6 ± 13.3
Age (year), n (%)
 ≤403531 (34%)35,310 (34%) 1298 (32%)12,980 (32%)
 40 < age ≤ 502574 (25%)25,740 (25%) 1139 (28%)11,390 (28%)
 >504303 (41%)43,030 (41%) 1668 (41%)16,680 (41%)
Residency <0.0001 <0.0001
 Northern Taiwan5605 (54%)52,331 (50%) 2330 (57%)20,771 (51%)
 Other areas4803 (46%)51,749 (50%) 1775 (43%)20,279 (49%)
Monthly income (NT$),
  median (IQR)
21,900
(1249–42,000)
21,900
(1249–36,300)
<0.000126,400
(1249–43,900)
21,900
(1249–42,000)
<0.0001
Monthly income (NT$), n (%) <0.0001 <0.0001
 ≤240006031 (58%)65,840 (63%) 2022 (49%)24,501 (60%)
 >240004377 (42%)38,240 (37%) 2083 (51%)16,549 (40%)
CCI score, mean ± SD1.7 ± 2.11.1 ± 1.7<0.00011.8 ± 21 ± 1.6<0.0001
CCI score, n (%) <0.0001 <0.0001
 =03419 (33%)54,800 (53%) 1287 (31%)21,585 (53%)
 =12642 (25%)23,116 (22%) 1039 (25%)9407 (23%)
 ≥24347 (42%)26,164 (25%) 1779 (43%)10,058 (25%)
Underlying diseases, n (%)
 Heart disease599 (6%)3030 (3%)<0.0001238 (6%)992 (2%)<0.0001
  Myocardial infarction157 (2%)819 (1%)<0.000169 (2%)314 (1%)<0.0001
  Congestive heart failure486 (5%)2467 (2%)<0.0001190 (5%)768 (2%)<0.0001
 Peripheral vascular disease204 (2%)1386 (1%)<0.000171 (2%)523 (1%)0.0146
 Major neurological disorder1306 (13%)7425 (7%)<0.0001547 (13%)2530 (6%)<0.0001
  Cerebral Vascular disease1230 (12%)6954 (7%)<0.0001522 (13%)2395 (6%)<0.0001
  Dementia136 (1%)870 (1%)<0.000149 (1%)217 (1%)<0.0001
  Hemiplegia111 (1%)891 (1%)0.02840 (1%)312 (1%)0.1365
 Chronic pulmonary disease3427 (33%)20,451 (20%)<0.00011407 (34%)7764 (19%)<0.0001
 Connective tissue disease360 (3%)1954 (2%)<0.0001131 (3%)665 (2%)<0.0001
 Peptic ulcer disease3633 (35%)23,157 (22%)<0.00011431 (35%)8976 (22%)<0.0001
 Liver disease2715 (26%)16,206 (16%)<0.00011188 (29%)6786 (17%)<0.0001
 Diabetes mellitus1544 (15%)10,849 (10%)<0.0001649 (16%)4212 (10%)<0.0001
 Renal disease593 (6%)3339 (3%)<0.0001238 (6%)1266 (3%)<0.0001
 Cancer668 (6%)4119 (4%)<0.0001236 (6%)1534 (4%)<0.0001
Abbreviations: SA = sleep apnea; NT$ = New Taiwan Dollar; CCI = Charlson Comorbidity Index; IQR = interquartile range.
Table 2. Incidence rate of temporomandibular disorder (TMD) after the index date in each group.
Table 2. Incidence rate of temporomandibular disorder (TMD) after the index date in each group.
VariablesStudy Arm A Study Arm B
Suspected SA Control ACrude
IRR [95% CI]
Adjusted
IRR [95% CI]
Probable SAControl BCrude
IRR [95% CI]
Adjusted
IRR [95% CI]
NTMDPYIRNTMDPYIRNTMDPYIRNTMDPYIR
Whole study population10,40818462,825.72.9104,080715634,918.61.12.3 [2.2–2.5] ***2.0 [1.9–2.1] ***41056623,204.82.841,050239232,345.41.03.1 [2.9–3.4] ***2.5 [2.3–2.7] ***
Stratified analyses
Sex
 Female38238823,241.73.838,230363236,432.31.51.7 [1.5–1.9] ***1.4 [1.3–1.6] ***894164716.43.489407647,813.41.61.5 [1.3–1.9] ***1.2 [1.0–1.5] *
 Male65859639,584.02.465,850352398,486.30.93.1 [2.9–3.3] ***2.7 [2.6–2.9] ***32115018,488.52.732,110163184,532.00.94.1 [3.8–4.5] ***3.4 [3.1–3.7] ***
Age
 ≤50610510038,558.22.661,050415389,170.41.12.3 [2.1–2.4] ***1.9 [1.8–2.1] ***24373714,314.62.624,370130143,405.10.92.5 [2.3–2.8] ***2.5 [2.2–2.8] ***
 >5043038424,267.63.543,030300245,748.21.22.4 [2.3–2.7] ***2.1 [2.0–2.3] ***1668298890.33.316,68010988,940.41.23.7 [3.3–4.1] ***2.6 [2.3–3.0] ***
Residents in
 Northern Taiwan560510333,394.53.152,331377317,999.41.22.3 [2.1–2.5] ***2.0 [1.8–2.1] ***23304013,652.22.920,771135116,551.11.23.1 [2.8–3.5] ***2.7 [2.4–3.0] ***
 Other areas48038129,431.22.851,749338316,919.21.12.4 [2.2–2.6] ***2.1 [1.9–2.2] ***1775269552.72.720,279104115,794.30.93.0 [2.6–3.4] ***2.3 [2.0–2.6] ***
Monthly income
 ≤NT$24,000603110136,318.02.865,840424398,259.61.12.5 [2.3–2.6] ***2.1 [2.0–2.3] ***202232111,75.72.924,501127135,874.10.93.8 [3.4–4.2] ***2.7 [2.4–3.0] ***
 NT$24,00043778326,507.83.138,240291236,659.01.22.2 [2.0–2.4] ***1.9 [1.8–2.1] ***20833412,029.22.816,54911296,471.31.22.5 [2.2–2.8] ***2.3 [2.1–2.6] ***
Comorbidity
 No (CCI score = 0)34195323,465.92.354,800307364,963.60.82.5 [2.3–2.7] ***2.4 [2.2–2.7] ***1287208271.12.421,585103132,653.40.83.1 [2.7–3.5] ***2.7 [2.4–3.1] ***
 Yes (CCI score ≥1)698913139,359.83.349,280408269,955.01.52.0 [1.8–2.1] ***1.9 [1.7–2.0] ***28184614,933.73.119,46513699,692.11.42.6 [2.4–2.9] ***2.4 [2.2–2.7] ***
The adjusted IRRs were calculated by multivariable analyses adjusting for sex, age, residency, income level, and the presence of various comorbidities (except for the stratifying variable). * p < 0.05; *** p < 0.0001. Abbreviations: SA = sleep apnea; N = number of patients; TMD = number of patients with temporomandibular disorder; PY = total patient-years; IR = incident rate (incidence of temporomandibular disorder per thousand-patient-years); IRR = incidence rate ratio; CI = confidence interval; CCI = Charlson Comorbidity Index.
Table 3. Multivariable Cox regression analyses of the factors contributing to temporomandibular disorder (TMD) in SA patients.
Table 3. Multivariable Cox regression analyses of the factors contributing to temporomandibular disorder (TMD) in SA patients.
VariablesSuspected SA PatientsProbable SA Patients
Maximal ModelReduced Model *Maximal ModelReduced Model *
HR [95% CI]p-ValueHR [95% CI]p-ValueHR [95% CI]p-ValueHR [95% CI]p-Value
Sex (male vs. female)1.0 [0.6–1.5]0.9027 1.4 [0.6–3.3]0.4233
Age > 50 (vs. Age ≤ 50)1.0 [0.6–1.7]0.8622 1.8 [0.8–3.8]0.1289
Residency in northern Taiwan
 (vs. in other areas)
1.0 [0.7–1.6]0.9741 1.3 [0.6–2.6]0.4782
Higher income (>NT$24,000)
 (vs. lower income)
1.0 [0.6–1.5]0.9131 0.9 [0.5–1.8]0.8036
Presence of underlying diseases
 (vs. absence of the diseases):
 Heart disease
 Peripheral vascular disease2.1 [0.7–6.9]0.2112 2.6 [0.6–11.5]0.1924
 Major neurological disorder0.6 [0.3–1.4]0.2387 0.8 [0.3–2.2]0.7341
 Chronic pulmonary disease1.0 [0.6–1.7]0.8444 0.7 [0.3–1.5]0.3687
 Connective tissue disease2.1 [0.9–4.9]0.0948 2.3 [0.7–8.0]0.1742
 Peptic ulcer disease2.0 [1.2–3.1]0.00562.0 [1.3–3.2]0.00172.7 [1.3–5.6]0.00813.1 [1.6–6.1]0.0010
 Liver disease1.4 [0.8–2.2]0.2191 1.2 [0.6–2.4]0.6526
 Diabetes mellitus1.4 [0.8–2.6]0.2388 1.7 [0.7–3.8]0.2247
 Renal disease0.3 [0.1–1.4]0.1267
 Cancer2.5 [1.3–4.7]0.00562.4 [1.3–4.4]0.00652.0 [0.7–5.5]0.1790
* Multivariable Cox regression model was built by a backward variable selection method, eliminating variables with the p value > 0.05. † Hazard ratio cannot be estimated because of small sample size. Abbreviations: SA = sleep apnea; HR = hazard ratio; CI = confidence interval.

Share and Cite

MDPI and ACS Style

Wu, J.-H.; Lee, K.-T.; Kuo, C.-Y.; Cheng, C.-H.; Chiu, J.-Y.; Hung, J.-Y.; Hsu, C.-Y.; Tsai, M.-J. The Association between Temporomandibular Disorder and Sleep Apnea—A Nationwide Population-Based Cohort Study. Int. J. Environ. Res. Public Health 2020, 17, 6311. https://doi.org/10.3390/ijerph17176311

AMA Style

Wu J-H, Lee K-T, Kuo C-Y, Cheng C-H, Chiu J-Y, Hung J-Y, Hsu C-Y, Tsai M-J. The Association between Temporomandibular Disorder and Sleep Apnea—A Nationwide Population-Based Cohort Study. International Journal of Environmental Research and Public Health. 2020; 17(17):6311. https://doi.org/10.3390/ijerph17176311

Chicago/Turabian Style

Wu, Ju-Hui, Kun-Tsung Lee, Chia-Yu Kuo, Chih-Hung Cheng, Jih-Yu Chiu, Jen-Yu Hung, Chung-Yao Hsu, and Ming-Ju Tsai. 2020. "The Association between Temporomandibular Disorder and Sleep Apnea—A Nationwide Population-Based Cohort Study" International Journal of Environmental Research and Public Health 17, no. 17: 6311. https://doi.org/10.3390/ijerph17176311

APA Style

Wu, J. -H., Lee, K. -T., Kuo, C. -Y., Cheng, C. -H., Chiu, J. -Y., Hung, J. -Y., Hsu, C. -Y., & Tsai, M. -J. (2020). The Association between Temporomandibular Disorder and Sleep Apnea—A Nationwide Population-Based Cohort Study. International Journal of Environmental Research and Public Health, 17(17), 6311. https://doi.org/10.3390/ijerph17176311

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