Relationship between job burnout and somatic diseases: a network analysis

Job burnout, characterized by feelings of exhaustion, cynicism and reduced personal efficacy, has been defined as a risk state for the development of diseases, but its association with somatic diseases is underexplored. Study participants were 5671 respondents (mean age 44.1 years, range 18–70; 38.6% men) to an online survey of job burnout, using a mobile health web application. Respondents provided data on sociodemographic characteristics, symptoms of burnout, measured with the Maslach Burnout Inventory-General Survey, depressive symptoms, measured with the Profile of Mood States, and 11 categories of somatic diseases. Adjusting for age, sex, educational level, depressive symptoms, and all disease categories included, network analysis showed a significant association of high exhaustion with “high blood pressure” and a category of “other chronic somatic diseases”. Low personal efficacy showed a significant association with “chronic lung diseases”. In network analysis, clinically relevant depressive symptoms were also significantly associated with “high blood pressure”, “other chronic somatic diseases” and “skin diseases”. Logistic regression analysis confirmed these associations. Burnout dimensions were significantly associated with an increased risk for somatic diseases, independently of sociodemographic factors and clinically relevant depressive symptoms. This relationship seems particularly evident for the exhaustion dimension of job burnout.


Scientific Reports
| (2020) 10 14 . Data were collected with Burnout Protector, a previously described mobile health web application 13 . Age below 18 was the only exclusion criterion. We did not consider data of respondents above 70 years, the age of the latest official retirement in Switzerland. We identified participants which answered the survey twice or more by their assigned user identification. We included only their initial answer in the analysis. This study was carried out in accordance with the Swiss Human Research Ordinance (i.e., under strict confidentiality and privacy, with coding of health-related personal data). Ethical authorization was not required, as this study does not fall within the scope of the Swiss Human Research Act. Nevertheless, all data was collected in accordance with scientific standards specified by Switzerland's Federal Act on Research Involving Human Beings 15 and in compliance with the national ethical guidelines on the handling of personal data in the field of medicine issued by the Federal Data Protection Commissioner of Switzerland 16 . All respondents provided informed consent to their anonymized data being used for research 13 . Measures. Respondents indicated their age, sex, and highest level of education. Somatic diseases were inquired with the closed-ended question "Are you suffering from any of the following somatic diseases?" referring to 11 categories. These were high blood pressure, high cholesterol, diabetes, stroke, myocardial infarction, chronic heart failure, chronic lung disease, chronic liver disease, chronic renal disease, chronic skin disease, and a residual category of "other chronic somatic disease". Participants could select the response option "other chronic somatic disease" if they suffered from a chronic disease which was not among the ten specified categories on the shortlist.
The 16-item MBI-General Survey (MBI-GS) was applied to assess EX (5 items), CYN (5 items) and PE (6 items) 17,18 . Each item is rated on a 7-point Likert scale from 0 ("never") to 6 ("daily"). In accordance with a previous population-based study from Finland, showing associations of both the total score of the MBI-GS and each of its dimensions with physical health outcomes, we used a mean score of 3.50 or higher across all items of a subscale to define high EX, high CYN and high lack of PE 8 . This categorization means that the three burnout dimensions are in the severe range (3.50-6) when symptoms of either EX, CYN or lack of PE are experienced about once a week or daily 8 . In a previous study on outpatients with work-related psychological problems, the MBI-GS showed the highest sum of sensitivity (78%) and specificity (48%) with a cutoff point of 3.50 on the EX subscale, which was recommended to minimize false negative burnout cases 19 . Moreover, individuals working under highly demanding conditions who scored above 3.50 points on the EX subscale of the MBI-GS showed greater physiological dysregulation, evidenced by lower heart rate variability, higher cortisol levels and higher blood pressure, than their counterparts with lower EX scores 10 . A 14-item subscale of the Profile of Mood States (POMS) was used to assess the severity of depressive symptoms 20,21 . Each item is rated on a 5-point Likert scale from 0 ("not at all") to 4 ("extremely"). Clinically relevant depressive symptoms were defined by a score of 14 or higher 22 . Statistical analysis. For network analysis and following logistic regression, categorical variables for high vs. low scores of each burnout dimension, clinically relevant vs. non-relevant depressive symptoms and high versus low educational level were used 13 . Age was dichotomized at a threshold of 45 years, as EX was shown to be predictive of increased all-cause mortality in employees under 45 years of age 23 . In all analyses, the included variables were dichotomous. The network analysis was performed using the eLASSO procedure 24 . Conceptually, for each node in the network, a regularized logistic regression is performed, with the given node as dependent variable (outcome) and all other included variables as independent variables (predictors) 24 . The chosen eLASSO procedure is a L1-regularization regression procedure based on LASSO (least absolute shrinkage and selection operator) 25 and adapted to analyze networks of dichotomous variables. In such a procedure, the penalization causes small estimates to shrink to zero. Hence, adjustment of individual p values is not necessary (for more details see 24 ). In a resulting network, nodes and edges represent the variables and the associations between variables, respectively, whereas the associations correspond to parameter estimates of a logistic regression model, whose magnitude can be transformed to odds ratio. In a second step, we investigated the associations of interest revealed by network analysis using multivariable, non-regularized logistic regression. In these models, the dependent variable (outcome) was a somatic disease of interest and all other variables also included in the network analysis were again used as independent variables. We tested no interaction terms. The variance inflation factor for all included variables in all tested models was below 3, thus indicating no multicollinearity. For the logistic regressions, all variables were entered collectively. The event rate of the included variables, their coding, the detailed results of the analyses (incl. goodness-of-fit indices and results of the significance testing of the individual parameters) are presented in the supplementary materials (Supplementary Tables S1-S4). All analyses were carried out in R (version 3.6.1) 26 . Estimation and visualization of the network was carried out using R-Package Bootnet (version 1.2.4) 27 . The analytic file ("Syntax") is provided as supplementary information. Table 1 shows the characteristics of the 5671 study participants. Figure 1 shows the results of the network analysis. The estimated network was sparse. Most edges were between somatic diseases, and between psychological variables, whereas only few edges were found between somatic diseases and psychological variables. Among them, high EX was associated with high BP and other chronic diseases, and low PE was associated with chronic lung diseases. Figure 1 additionally shows that clinical depressive symptoms were associated with high BP, skin diseases  Table S4).

Discussion
Performing network analysis, we found significant associations of high EX with high BP and a residual category of other chronic somatic diseases, and of low PE with chronic lung diseases, but no associations with high CYN, in respondents to an online survey of job burnout. These associations were adjusted for sociodemographic factors, clinically relevant depressive symptoms, and other somatic diseases and confirmed by logistic regression analysis. More precisely, using a previously proposed cut-off of 3.50 for each burnout dimension 8 , we found high EX to be associated with a 40% and 33% increased relative risk of high BP and other chronic somatic diseases, respectively. Low PE was significantly associated with a 47% increased relative risk of chronic lung diseases. Among other associations, we also observed a significant association of clinical depressive symptoms with high EX on the one hand and with high BP, skin diseases and other chronic somatic diseases on the other. These observations confirm an overlap between depression and burnout 11,12 , but also suggest partially independent associations of depression and burnout with somatic diseases. The overlap between high EX and clinical depressive symptoms on the one hand and their somatic outcomes on the other is perhaps not so surprising, since feelings of exhaustion are part of depressive symptomatology.
Our findings concur with results from some previous studies, but with partly different definitions of burnout as a single score or cut-off scores for high EX, high CYN and low PE than those we used in our study. The advantages and disadvantages of different burnout definitions and measurements have been acknowledged elsewhere 3,28 , and were not part of our study, but must be kept in mind when comparing results between different studies. For instance, high EX, defined by the same cut-off as in our study, was associated with elevated BP in individuals working under highly demanding conditions 10 . Moreover, the total score of the School Burnout Inventory (EX at school, CYN toward the meaning of school, sense of inadequacy at school) was associated with elevated 24-h ambulatory BP in female students 29 . A burnout total score consisting of EX and CYN symptoms was associated with self-reported treatment for hypertension in South African employees 30 . In 5000 civil servants from China, a higher burnout total score (sum of emotional exhaustion, depersonalization, and reduced accomplishment of a 15-item MBI-Human-Services Survey) was associated with a higher risk of chronic diseases, including hypertension, coronary heart disease, stroke, and diabetes 31 .
We did not find an association of burnout dimensions with CVD or cardiovascular risk factors, except high BP. Limited statistical power could explain this, as the prevalence of CVD was low in our sample. The diseases categorized under the residual category "other chronic somatic diseases" are subject to speculation, as this was a www.nature.com/scientificreports/ residual category which could be selected if participants had a disease form a category which was not specified on the shortlist. For instance, in Turkish physicians, peptic ulcers have previously been associated with high emotional exhaustion, defined by a cut-off of 2.0 on the emotional exhaustion subscale of the 22-item MBI-Human Services Survey 32 . Moreover, as the category of other chronic somatic disease showed a strong relationship with liver disease, and to a lesser extent with pulmonary and skin diseases, one option could be inflammatory bowel disease or systemic diseases with organ manifestations. Unlike a previous study applying a clinical health examination, applying the same cut-off as used in our study 8 , we found the association of low PE with self-reported chronic lung diseases to be significant. Behavioral and biological mechanisms might be involved in the association between burnout and somatic diseases. For instance, high EX, defined by the same cut-off as in our study, has been associated with low physical activity and heavy alcohol drinking 33 . Also, one explanation for the association between low PE and chronic lung disease could be smoking, which has been shown to be more frequent in employees in the health care sector with low PE than in those with high PE 34 . Greater total frequency and total severity of burnout symptoms, measured with a 23-item early version of the MBI-Human Services Survey, were associated with increased allostatic load, denoting chronic stress-induced biological burden leading to diseases via neuroendocrine, immune and autonomic pathways 35 . However, of the three burnout dimensions, only the emotional exhaustion symptom frequency score was significantly associated with allostatic load, whereas depersonalization and lack of personal accomplishment were not 35 . In contrast, we found no significant association between high CYN and any somatic disease category. The absence of this relationship remains to some extent speculative. For instance, a  www.nature.com/scientificreports/ previous study found an association between depersonalization and alcohol abuse in health care workers 34 . However, since our sample was only 44 years old on average, excessive alcohol consumption may not yet have led to manifest somatic disease in a significant proportion of those with high CYN. Another explanation could be limited statistical power, as only 10.4% of our study participants reported high CYN. The use of network analysis and the large sample size are strengths of our study, which has also its limitations. An investigation of users of a web application, including medical doctors, limits the generalizability of our findings. In addition, we did not validate our results in a second, independent dataset, so the external validity of our results is currently unknown. There was insufficient occupational information regarding job content and environment of study participants; however, it has been argued that the dimensions of the MBI-GS, used in our study, can be applied across all work domains 28 . Self-reported diseases are less accurate than clinical diagnosis. The network model does not allow for unidirect associations. For instance, whereas it is possible that EX has led to high BP through increased sympathetic activity, inflammation, accompanying high BP, could also have triggered sickness behavior including feelings of exhaustion and depressive mood. There could be residual confounding through unmeasured health behaviors. There is no definite consensus on cut-offs that may best define clinically relevant symptoms of burnout and depression with the MBI and POMS, respectively, and results might vary across different psychometric instruments. Taken together, our network analysis provides evidence that burnout dimensions are significantly associated with an increased risk of somatic diseases. The observed associations were independent of sociodemographic factors and clinically relevant depressive symptoms and more apparent for high EX compared with the other dimensions of the burnout experience.

Data availability
The datasets analysed during the current study are available from the corresponding author on reasonable request. www.nature.com/scientificreports/