Temporomandibular dysfunction experience is associated with oral health-related quality of life: an Australian national study

Background There are very few studies of the association between temporomandibular dysfunction (TMD) and oral health-related quality of life (OHRQoL) in a representative sample from the Asia–Pacific region. Accordingly, we aimed to quantify the association of TMD with OHRQoL dimensions and overall measurement scores in a representative sample of Australian adults while accounting for a range of confounders, and statistically estimating whether TMD experience is meaningfully associated with OHRQoL. Method Australia’s National Survey of Adult Oral Health 2004–2006 data were used. The outcome variables were the Oral Health Impact Profile (OHIP-14) domains and overall scores while the main exposure was self-reported Diagnostic Criteria Question for TMD. The analysis accounted for confounders including oral health status obtained from the oral examination, demographics, socioeconomics, health behaviours and health including perceived stress subscales of the PSS-14. We conducted complex samples analysis while using Cohen’s f2 effect size to estimate whether the association is meaningful. Results TMD prevalence was 9.9% (95% CI: 8.4–11.6%) among 4133 Australian adults. TMD experience was associated with impairments to the seven OHIP-14 OHRQoL domains (P < .05) with higher impairments observed in physical pain (B = 0.82, 95% CI: .45–1.20, P < .001), psychological discomfort (B = 0.68, 95% CI: .29–1.06, P = .001) and psychological disability (B = 0.52, 95% CI: .20–.84, P = .001) in adjusted multivariate analyses. The difference in the mean OHIP-14 scores for those reporting TMD (Mean = 13.1, 95% CI: 12.0–14.0) and those who did not (Mean = 6.6, 95% CI: 6.0–6.8) was significant (t = 7.51, P < .001). In an adjusted multivariable model for OHIP-14 scores, TMD experience was associated with higher OHIP-14 scores (B = 3.34, 95% CI: 1.94–4.75, P < .001) where the Cohen’s f2 was .022. Further, perceived stress subscales: perceived distress and perceived control were associated with TMD experience and OHIP-14 scores (P < .05). Conclusion Lower OHRQoL was observed in Australian adults who reported TMD experience but with small clinical importance which might support considering TMD in regular dental care. The higher impairments observed in physical pain, psychological discomfort and psychological disability domains of OHRQL can help clinicians and researchers focus their attention on these domains. The confounding effect exhibited by the perceived stress subscale might support their inclusion in future TMD and OHRQoL research. Supplementary Information The online version contains supplementary material available at 10.1186/s12903-021-01773-z.


Introduction
Temporomandibular dysfunction (TMD) is a group of degenerative musculoskeletal disorders that affects the morphology and function of the masticatory system [1]. In Australia, TMD is experienced by almost 10% of Australian adults [1] where evidence suggests that TMD symptoms might be negatively associated with individual's perception of their physical, mental and social wellbeing about oral health [2][3][4] which is known as oral health-related quality of life (OHRQoL) [5]. While associations between components of oral health status such as caries experience and periodontal diseases, and OHRQoL have been investigated among the Australian adult population [6,7], there is a lack of assessment of the association between TMD experience and OHRQoL among the Australian adult population. Understanding the OHRQoL impairment associated with TMD experience might be important for oral health education and treatment-seeking considering a minority of TMD patients (3-7%) are seeking professional advice [8]. While OHRQoL is known to be a multidimensional concept [5], there are limitations for studies from Australia that assessed which OHRQoL dimensions might be impaired by TMD experience.
TMD experience might be associated with impairment to several dimensions of OHRQoL. For instance, orofacial pain is the most common symptom experienced by TMD patients which severely affects the individual's QoL [3,9]. Furthermore, it is observed in another oral health condition that severe pain might interact with other QoL aspects such as mood status and the ability to perform daily activities such as work or study [10]. Besides, TMD patients experience chewing difficulty which is known to be negatively associated with OHRQoL [11]. This chewing difficulty often results in a dietary modification to reduce symptoms. The dietary modifications also constitute part of the conservative management plans for TMD [12], which is negatively associated with food enjoyment-an important OHRQoL aspect. Moreover, TMD patients might experience an audible click or grating [12] which might affect the social life of TMD patients. Understanding TMD association with dimensions and the overall OHRQoL among a representative sample of the Australian adult population might be important for clinical practice and health services research considering that researchers who conducted systematic reviews of the association between TMD and OHRQoL did not include representative population-based samples [2,3]. Beyond investigating the association of TMD experience with OHRQoL, it is becoming necessary to capture whether TMD experience has a clinically relevant association with the OHRQoL measurement-a concept known as the Minimally Important Difference (MID).
There is a growing demand to estimate the Minimally Important Difference (MID) in measures of OHRQoL between different groups based on their disease status in the population [13]. Determining the clinical relevance of the association between TMD experience and OHRQoL is important for decision-making on whether a treatment is needed as well as to choose the intervention to be provided (conservative or invasive) based on evidence of whether the intervention improves OHRQoL [14,15]. Determination of MID is challenging and several methods are suggested [15]. While some researchers used the anchorage method to establish the MID points for specific intervention on OHRQoL measure [16], Tsakos et al. [17] suggested calculating the effect size (ES) of the mean difference between groups of disease/condition status in cross-sectional data to statistically estimate the MID of the associations with OHRQoL measures. The use of effect size to estimate MID is argued to be less biased and makes it easy to compare between studies with different outcome measures or measures of exposure [15]. Several measures of standardized effect size are available however, Cohen's f 2 is argued to be a useful measure of effect size estimation considering it could estimate the effect size from multivariable hierarchical models [18] which account for potential confounders. Among the important confounders for OHRQoL such as biomedical and socioeconomic factors that needed to be accounted for when examining the association between TMD and OHRQoL, there has been a demand to include psychological factors [2] in assessing such a relationship.
Perceived stress, which refers to the extent an individual feels stressed with life events, is an important psychological factor that affects some aspects of oral health status. For instance, perceived stress is associated with recurrent aphthous ulcer [19], mechanical wear of teeth due to bruxism [20], reporting of dry mouth [21], poor perceived oral health [22], and perceived (work and non-work related) stress [1,23]. While these variable associations are recognised in the literature, the earlier quantifications did not account for sampling biases adequately due to the sampling restrictions in those studies included in Dahlström and Carlsson [2] systematic review. National studies with representative samples can play an important role in providing for such quantification and help triangulate the association of TMD with OHRQoL. Thus, this study aimed to quantify the association between TMD experience and OHRQoL domains and overall measurement scores while accounting for various confounders as well as estimating the effects sizes for TMD experience association with OHRQoL measurement score for necessary across-study comparisons and statistical estimation of Minimally Important Difference (MID-S) for the association between TMD and OHIP-14.

Study design
We used the National Survey of Adult Oral Health (NSAOH) 2004-2006 wave [24] which is a cross-sectional three staged random stratified clustered representative sample of Australia's adult residents aged 15-years or over who resided in a household with access to a telephone line that was listed in Australia's "Electronic White Pages" [25] with their corresponding postcode. The first sampling stage selected the postcodes followed by the second sampling stage where households were selected from the selected postcodes. The third sampling stage selected the target person from the selected household. Stratification of the selected postcodes from all Australian states and territories was performed and included two strata: capital city and rest of the state, however, the Australian Capital Territory (ACT) was considered as a single stratum major city. Further details of NSAOH sampling are reported in Slade et al. [26]. The NSAOH data used in this study comprised three datasets: Computer Assisted Telephone Interview (CATI), Self-Complete Questionnaire and a Standardized Oral Epidemiological Examination.

Setting
Australia's residents in all states and territories aged 15-years or over who randomly invited to join the NSAOH. NSAOH data were collected between July 2004 and September 2006.

Participants
The randomly selected adult aged 15-years or older present in the household at the time of the CATI interview was invited to join the study to proceed with the interview and if declined, it was reported as a declined response. Participants who completed the CATI were asked if they were happy to receive a mailed questionnaire and participate in the oral epidemiological examination. Participants who agreed to complete other parts of the study were contacted later for the mailed questionnaire and the oral epidemiological examination. For the oral epidemiological examination, participants had a medical questionnaire to determine their eligibility for the epidemiological oral examination according to the adopted protocol [27] such as rheumatic fever, endocarditis, bleeding disorders, joint replacement within the past 3-months, etc.

Sample size, response rate and quality check
The target sample size of the NSAOH was calculated for key outcomes such as the Decayed, Missing and Filled teeth to detect, on a national level in the age-specific estimates, 10% reduction on the national level compared with 1988 survey at 80% power and was 7500 participants. Considering that TMD was not one of the outcomes considered for sample size calculation, we provided posthoc power analysis. The average response rate for the NSAOH 2004-06 by the selected postcode was 49.0% for the interview and 44% for the in-scope oral examination [26]. The non-response bias was investigated using the "population benchmark" and "small area socio-economic indicators" approaches and found to be small (less than 3%) for most oral health indicators [26]. Missing data were examined for key variables and was small (0-4.4%) which found to be at random and were not imputed. NSAOH data were weighted to account for sampling design and to match the age and gender distribution for the selected strata which were provided by the 2005 Australian Bureau of Statistics Estimated Residential Population Data [26]. Participants with complete records across the three datasets were included in this study.

Ethical considerations
Ethical approval for the NSAOH was obtained from the University of Adelaide Human Research Ethics Committee with approval number: H-001-2004. Informed consent was obtained from all NSAOH participants. For the CATI interview, a verbal consent was obtained from the participants before proceeding with the CATI interview due to the nature of the survey method, and was included in the ethics approval. Written consent was obtained from participant for the self-complete questionnaire and the standardized oral epidemiological examination. Participants found on their oral examination to have a concerning medical or dental condition were referred appropriately to health professionals.

Self-complete questionnaire
OHRQoL (outcome variable) measured by the Oral Health Impact Profile-Short form (OHIP- 14) The mailed self-completed questionnaire collected data on OHRQoL using the 14-item Oral Health Impact Profile-short-form (OHIP-14) [28], a shorter version of the original OHIP-49. The OHIP-14 consists of seven domains with two items that represent each domain. OHIP-14 domains include functional limitation, physical pain, psychological discomfort, physical disability, psychological disability, social disability and handicap. Items for the OHIP-14 are scored on a 5-point scale ranging from 0 for "Never" to 4 for "Very often". The reference period of the OHIP-14 is "over the past 12 months". The total OHIP-14 scores range from 0 to 56 with a higher score meaning a lower OHRQoL. OHIP-14 was the most commonly used generic measure for oral conditions associated with OHRQoL with evidence of reliability and validity across a wide range of socio-cultural contexts [29].

TMD experience (main exposure variable) measured by the Diagnostic Criteria Question for TMD
TMD experience was assessed using the Diagnostic Criteria Question for TMD, which was adapted from a Canadian study where the questionnaire is reported to have 73% sensitivity and 75% specificity in predicting clinical diagnosis for TMD [30]. The questionnaire ( Fig. 1) consists of seven items in two domains: pain (three questions) and symptomatic TMD (four questions) which represent functional disturbance. The criterion adopted for the presence of TMD was recording a positive response for one or more of the pain items and a positive response in one or more of the symptomatic TMD items which are in line with the clinical diagnosis of TMD used in previous research [31]. This research diagnostic criteria have been recommended by Sanders and Slade [32] and was used by researchers in a recent study [1].

Perceived stress measured by the 14-item Perceived Stress Scale (PSS-14) (confounder variable)
The Perceived Stress Scale (PSS-14) was developed by Cohen et al. [33] as a measure of subjective stress and consists of 14 items representing how unpredictable, uncontrollable, and overloaded participants feel in their life. The reference period of the PSS-14 is"last month". Items for the PSS-14 are scored on a 5-point Likert-like scales ranging from 0 for 'Never' to 4 for 'Very often' . The PSS-14 consists of two subscales: Perceived distress (negative subscale) and Perceived control (positive subscale). The perceived distress subscale is the sum score for items 1,2,3,8,11,12 and 14 where higher scores represent higher perceived distress. The perceived control subscale is the sum score of items 4,5,6,7,9,10 and 13 where higher score represent higher perceived control (coping).

Self-reported static malocclusion (confounder)
The NSAOH self-complete questionnaire asked the participants whether their 'Teeth feel do not fit properly together' whereas the response was yes/no. We have used this question as a proxy of static malocclusion.

Alcohol consumption (confounder)
The NSAOH participants were asked about the 'Days per week of alcohol drinking' and 'The number of standard drinks per day' .

The Computer Assisted Telephone Interview (CATI) data
The CATI collected demographic and socioeconomic characteristics (included as confounders or predictors in our analysis) such as age, gender, educational attainment, current employment status, annual household income, having private dental insurance, indigenous status and country of birth. Further, NSAOH collected information on chronic diseases experience such as diabetes and behaviours such as current smoking status.

Standardized oral epidemiological examination data
The standardized dental examinations were carried out by 30 dentists who received training to comply with the oral examination protocol [34]. The intra-examiner reliability was checked for five participants for each dentist with replicated examinations conducted by the principal survey examiner for teeth present, caries experience and periodontitis. The median inter-examiner reliability was 0.85. The oral examination was carried out using a selfilluminated dental mirror and a periodontal probe (with 2 mm marking). The oral examination included:

Number of missing teeth (confounder)
The number of missing permanent teeth for any reason and not replaced by a removable or fixed prosthesis was recorded starting from the upper right quadrant at the position of tooth 18 to 28 then the lower left quadrant at tooth position 38 to 48.

Number of teeth with untreated decay (confounder)
The NSAOH examiners visually inspected (using the self-illuminated dental mirror and periodontal probe gently used to check the surface texture) the five surfaces of the present permanent teeth (occlusal, mesial, distal, buccal and lingual) for decay in enamel, dentine or involved enamel and dentine and recorded as 'D' . Recording of teeth with untreated decay started from tooth 18 if present to 28 then 38 to 48. The number of teeth with untreated decay was computed by summing the number of present teeth with untreated decay lesion on one or more of the examined surfaces.

Number of filled teeth (confounder)
The oral examiners visually inspected tooth surfaces for fillings using the same procedure for recording decayed teeth. The number of filled teeth was computed by summing the number of present teeth with one filling or more.

Periodontal diseases (predictor)
Participants of the NSAOH were examined for recession and probing depth which were used to compute clinical attachment loss (CAL). Measurements were obtained from three points: mesiobuccal, mid-buccal and distobuccal for all teeth present except for third molars. Periodontitis status was assessed using the US Centre for Disease Control and Prevention case definition of periodontitis [35]. For this study, we created a binary variable (Yes/No) for the presence of 'moderate or severe periodontitis' where moderate periodontitis was defined as the presence of ≥ 2 sites of inter-proximal CAL of ≥ 4 mm, and severe periodontitis was defined as the presence of ≥ 2 sites of inter-proximal CAL of ≥ 6 mm.

Data analysis plan
Data were analysed using SPSS Statistics for Windows version 27 [36] and STATA 14 IC [37]. We conducted a Confirmatory Factor Analysis (CFA) to check the structural validity of the used DCQ for TMD using the generalized structural equation model (Additional file 1: Figure  S1 and Table S1 represented CFA of the two-dimensional DCQ for TMD model, and Figure S2 and Table S2 represented the CFA for the one-dimensional DCQ for TMD model). We found that the two-dimensional structure of the DCQ for TMD had slightly better goodness of fit statistics compared with the one-dimensional structure. To account for the sampling weights, strata and clusters of the NSAOH, a complex sample analysis module in SPSS were used in all conducted analyses [38]. The characteristics of the participants were obtained using complex sample descriptive statistics using percentages or mean of population estimates with a 95% confidence interval as well as a distribution analysis of OHIP-14 scores.
To explore factors associated with TMD experience as an exposure, we conducted bivariate associations with potential confounders using complex sample bivariate logistic regression models. Further, we conducted complex samples Independent Sample t-test for OHIP-14 domains and total scores and TMD experience as well as conducted a posthoc power analysis. Moreover, we computed Cronbach's alpha reliability for the used scales (OHIP-14 and PSS-14 subscales). We used adjusted multivariate complex samples linear regression models for OHIP-14 domains and total scores for the association with TMD experience. To examine the association between TMD experience and OHIP-14 scores while accounting for various confounders, multivariable blockwise complex sample linear regression models for OHIP-14 scores were conducted using consecutive blocks. Block 1: TMD experience, Block 2: Block 1 + components of oral health status (Moderate or severe periodontitis, number of decayed, missing and filled teeth, and whether teeth are not fit together properly), Block 3: Block 2 + demographics (sex, age group, birthplace and Indigenous status), Block 4: Block 3 + socioeconomics (educational attainment, current employment status, annual household income level and having private dental insurance) and Block 5: Block 4 + health and health behaviours (diabetes status, perceived stress subscales, current smoking status and alcohol consumption). To estimate the minimally important difference statistically (MID-S) of TMD association with OHIP-14 scores, we used Cohen's f 2 as an appropriate measure of effect size [18,39]. The Cohen's f 2 was calculated using the following equation: [18] where R 2 AB was the OHIP-14 variance explained in multivariable linear regression by TMD experience together with a set of confounders and R 2 A is the OHIP-14 variance explained by confounders in the multivariable linear regression (excluding the TMD experience). We used Cohen's f 2 index to determine the magnitude of the effect size where f 2 ≥ 0.02, f 2 ≥ 0.15, and f 2 ≥ 0.35 represent small, medium, and large effect sizes, respectively [18,39]. Supplementary analyses were presented in the Additional file 1: Supplementary analyses that included a confirmatory factor analysis for the DCQ for TMD, examination of regression residuals and an alternative complex samples linear regression models for log-transformed OHIP-14 overall and domains scores and included a confirmatory factor analysis for the DCQ for TMD, examination of regression residuals and an alternative complex samples linear regression models for log-transformed OHIP-14 overall and domains scores.

Discussion
In this study, we aimed to quantify, among a representative sample of Australian adults, the association between TMD and OHRQoL dimensions and overall measurement scores while accounting for the hierarchical structure of the data and the relevant confounding variables.
We found that TMD experience was associated with impairments to all seven OHRQoL domains measured by the OHIP-14 with higher impairments observed in the physical pain, psychological discomfort and psychological disability domains of the OHIP-14. Furthermore, on examining the association between TMD experience and the overall OHRQoL measured by OHIP-14 scores using multivariate regression models, TMD had a significant negative association with OHRQoL and this association remained significant while accounting for various confounders. However, when we estimated the minimally important difference of the association between TMD experience and the OHRQoL measure statistically using a standardized estimate of effect size, it was found to be small. Moreover, we found that higher levels of perceived distress and perceived control were associated with TMD experience and OHRQoL suggesting a confounding effect.
The principal finding of this study was that TMD experience was associated with impairments to all OHRQoL domains measured by the OHIP-14 suggesting the broad extent to which TMD impaired the participant's perception of their wellbeing. While it might not be surprising that TMD experience was associated with the physical pain domain [3] with a small clinical relevance observed, It was interesting to observe associations with psychological discomfort and psychological disability in relation to OHRQoL. This observation might explain why pharmacological treatment with non-steroidal anti-inflammatory drugs (NSAID) and selective serotonin reuptake inhibitors (SSRIs) showed effectiveness in improving OHRQoL in TMD patients [14]. Moreover, TMD experience was associated with poorer overall OHRQoL. This finding is consistent with smaller studies included in the Dahlström and Carlsson [2] and Oghli et al. [3] systematic reviews.
We attempted to estimate whether TMD experience was clinically relevant using a statistical approach to determine relevance for clinicians, health services providers and health funds to offer treatment and whether the intervention provided is effective [41]. While effect size of mean difference of OHIP-14 scores across groups of a disease or a condition in cross-sectional data is recommended by Tsakos et al. [17] as a statistical method to estimate MID, we argued that using Cohen's f 2 to estimate the effect size might provide a standardized effect size estimate of TMD experience association with OHIP-14 score while considering confounders for OHIP-14 scores. The standardized effect size estimate will enable researchers to compare our findings across studies and OHRQoL outcome measures for association with TMD experience as well as it is argued to be less biased [15]. When we assessed the effect size of TMD experience association with OHIP-14 scores, it appeared to be small suggesting a small MID-S. The observed small MID-S of TMD experience association with OHRQoL measure might be interpreted with caution considering the OHIP-14 is a generic OHRQoL measure and other conditionspecific OHRQoL measures such as the OHIP-TMD might show a different MID-S [42]. Accordingly, our findings might support screening for TMD when providing regular dental care similar to the current practice for caries and periodontitis experiences. Further, health services researchers might need to consider TMD as a component of oral health status when conducting OHRQoL research.
In our study, we included the psychological factor of perceived stress as it might be associated with oral health and oral health perception [19][20][21]. We found that increased perceived distress and perceived control subscales of the perceived stress scale were associated with TMD experience and OHRQoL suggesting a confounding effect. While it might be argued that OHIP-14 and PSS-14 are measuring psychological aspects, the measured psychological constructs by the two instruments are different which reduced any potential of bias. Accordingly, future research investigating TMD association with OHRQoL might need to consider psychological factors Table 2 OHIP-14 domains and total scores among the Australians and associations with TMD experience status 1 Analyses accounted for cluster and stratum used in NSAOH sampling strategy, as well as sampling weights to ensure representativeness of the estimates 2 Multivariate complex samples linear regression models for OHIP-14 domains and total score association with TMD experience status adjusted for oral health status (Moderate or severe periodontitis, number of decayed, missing and filled teeth, and whether teeth are not fit together properly), demographics (sex, age group, birthplace and Indigenous status), socioeconomics (educational attainment, current employment status, annual household income level and having private dental insurance) and health and health behaviours (diabetes status, perceived stress subscales, current smoking status. and alcohol consumption) 3 The provided f 2 is for TMD association with OHIP-14 domains and total score using the equation presented in the methods section     such as perceived stress in the analysis similar to the current practice of including biomedical and socioeconomic factors. One limitation of our study was the use of the diagnostic criteria question for TMD which varies from proposed TMD questions by the International Network of Orofacial Pain and Related Disorders Methodology [43]. Further, the NSAOH oral epidemiological examination protocol did not include clinical examination for TMD [34]. However, the used TMD questionnaire has been validated using clinical data in another study [30] and reported to have a high level of sensitivity and specificity with the clinical diagnosis of TMD [32] which encouraged its use in past studies [1,32]. It is also pertinent to note that the perceived impairment found in the measure of TMD and OHRQoL gives credence to the presence of a real impairment felt by the individuals with TMD. Also, we acknowledged the temporality of the observed associations between TMD and OHRQoL given the cross-sectional design of the study. While we argued that OHIP-14 scores have an acceptable normal distribution [40] for using the linear regression model, close findings were obtained via an alternative modelling approach using the complex samples linear regression for log-transformed OHIP-14 overall and domains scores (Additional file 1: Table 3 and Table 4). Our study carried limitations associated with secondary data analysis which limited our ability to use a current case definition of periodontitis however, this has a limited effect on our finding as periodontitis surveillance was not the focus of this study. Further, the chronic medical conditions included in our analysis was limited to diabetes since data on other chronic medical conditions were not collected as part of the NSAOH such as osteoarthritis. However, this unlikely to impact our findings considering that osteoarthritis is common in the older age group and might only impact a limited subgroup of TMD patients.
This study adds to the knowledge in the field by identifying the negative association between TMD and OHRQoL domains and overall measurement scores using a representative sample of the Australian adult population. Further, our analysis has adjusted for a range of biomedical, psychological and socioeconomic confounders that might not be included in a single study due to the limitation of such a quality dataset which adds to the reliability of our findings. Our findings contributed to clinical practice by highlighting TMD relevance, as a component of oral health status, to OHRQoL-an important outcome for healthcare services. In addition, our study was able to explain the reported effectiveness of pharmacological treatment of TMD using NSAIDs and SSRIs in improving OHRQoL in TMD patients considering our findings of higher observed associations of TMD experience with physical pain, psychological discomfort and psychological disability domains of OHRQoL. Besides, this study has revealed the confounding effect of psychological factors measured by perceived stress subscales in the relationship between TMD experience and OHRQoL which might support its inclusion in future TMD and OHRQoL research.

Conclusion
In a representative sample of the Australian adult population, TMD experience was associated with lower OHRQoL, and this association remained significant after accounting for biomedical, psychological, and socioeconomic confounders. This finding supports that clinicians and health services providers might need to consider TMD in the overall OHRQoL of patients receiving dental care since OHRQoL is an end outcome of the provided oral healthcare services. When we used the standardized effect size estimate to determine whether TMD experience is clinically relevant to our participants' OHRQoL in an adjusted analysis, we found that it was of a small clinical relevance. Considering that we used a generic OHRQoL measure, clinicians might need to consider individual variation in the perception of how TMD experience might be associated with the patient's OHRQoL during decision-making about the need to provide an intervention (conservative or invasive). Besides, we found that TMD experience was associated with impairments to all OHRQoL domains measured by the OHIP-14. Also, we found that higher impairments were observed in physical pain, psychological disability and psychological discomfort which might explain the reported effectiveness of NSAID and SSRIs in improving the OHRQoL in TMD patients. Further, we found that increased levels of perceived distress and perceived control subscales of the perceived stress scale were associated with TMD experience and OHRQoL suggesting a confounding effect that might support the inclusion of perceived stress in TMD and OHRQoL research similar to the practice with biomedical and socioeconomic confounders. and an alternative complex samples linear regression models for logtransformed OHIP-14 overall and domains scores.