Skip to main content
  • Research article
  • Open access
  • Published:

Coronary microvascular function in male physicians with burnout and job stress: an observational study

A Correction to this article was published on 16 February 2024

This article has been updated

Abstract

Background

As a professional group, physicians are at increased risk of burnout and job stress, both of which are associated with an increased risk of coronary heart disease that is at least as high as that of other professionals. This study aimed to examine the association of burnout and job stress with coronary microvascular function, a predictor of major adverse cardiovascular events.

Methods

Thirty male physicians with clinical burnout and 30 controls without burnout were included. Burnout was assessed with the Maslach Burnout Inventory and job stress with the effort-reward imbalance and overcommitment questionnaire. All participants underwent myocardial perfusion positron emission tomography to quantify endothelium-dependent (cold pressor test) and endothelium-independent (adenosine challenge) coronary microvascular function. Burnout and job stress were regressed on coronary flow reserve (primary outcome) and two additional measures of coronary microvascular function in the same model while adjusting for age and body mass index.

Results

Burnout and job stress were significantly and independently associated with endothelium-dependent microvascular function. Burnout was positively associated with coronary flow reserve, myocardial blood flow response, and hyperemic myocardial blood flow (r partial = 0.28 to 0.35; p-value = 0.008 to 0.035). Effort-reward ratio (r partial =  − 0.32 to − 0.38; p-value = 0.004 to 0.015) and overcommitment (r partial =  − 0.30 to − 0.37; p-value = 0.005 to 0.022) showed inverse associations with these measures.

Conclusions

In male physicians, burnout and high job stress showed opposite associations with coronary microvascular endothelial function. Longitudinal studies are needed to show potential clinical implications and temporal relationships between work-related variables and coronary microvascular function. Future studies should include burnout and job stress for a more nuanced understanding of their potential role in cardiovascular health.

Peer Review reports

Background

Occupational burnout can be described as a condition characterized by a sense of exhaustion, detachment from one’s work, and decreased productivity resulting from prolonged exposure to unmanaged job stress [1]. This conceptual framing of burnout is bolstered, for instance, by a recent systematic review of longitudinal studies, which demonstrates that high job demands and negative job attitudes are predictors of occupational burnout with effect sizes ranging from small to medium [2]. With a prevalence of about 50%, burnout is high in physicians [3] who also report higher job stress compared with the general population [4, 5]. Physician burnout is associated with poor mental health, including depression and suicidal ideation, reduced productivity, early job retirement, and lower patient satisfaction [6]. Similar negative effects on physician well-being and patient care have been found to result from increased job stress [7, 8]. Accumulating evidence also shows associations of both burnout and job stress with an increased risk of somatic diseases [9, 10], particularly coronary heart disease (CHD) [11, 12]. Specifically, among 8838 seemingly healthy employees in Israel, elevated initial levels of burnout were indicative of an elevated CHD risk after a follow-up of 3.4 years [11]. This association remained significant even after adjusting for demographics, cardiovascular risk factors, job stress, and depression [11]. Despite their medical background, physicians have at least as high a risk of cardiovascular disease (CVD) as other professionals [13, 14].

Explanations for the observed associations of burnout and job stress with an increased CHD may include unhealthy lifestyle behaviors and psychobiological changes like autonomic and neuroendocrine dysfunction and low-grade inflammation [15, 16]. However, the underlying pathophysiology of CHD risk in physicians with burnout and job stress is unknown. Moreover, examining the independent associations of these related constructs may lead to a more nuanced understanding of their impact on cardiovascular health. For instance, one meta-analysis found that the cortisol awakening response was negatively associated with burnout, but positively so with job stress [17].

With respect to cardiovascular health, the coronary microcirculation plays a pivotal role in modulating coronary blood flow to meet the oxygen demands of the heart [18]. Coronary microvascular dysfunction is expressed as reduced coronary flow reserve (CFR) which can be non-invasively assessed through positron emission tomography (PET) myocardial perfusion imaging [19]. Reduced CFR serves as an indicator of coronary artery disease severity and has been prospectively associated with an increased risk of major adverse cardiovascular events and all-cause mortality [20]. To the best of our knowledge, there have been no studies that have explored the correlation between both burnout and job stress with CFR. Therefore, the aim of this study was to examine the independent association of burnout and job stress with CFR in male physicians.

Methods

Participants

The local ethics committee Zurich (BASEC-Nr. 2018–01974) approved the study, and written informed consent was obtained from all participants. Between 09/2019 and 12/2021, we recruited male physicians in Switzerland to investigate the relation of burnout to cardiovascular health. Table 1 provides names and explains in straightforward language key parameters that were measured in this study. Various channels, including hospitals, clinics, physician associations, and direct email communication, were used to reach potential participants. Interested physicians were provided with study details and objectives through text/flyer and had the option to contact the study management. The study aimed to enroll 60 participants who were assigned to either the burnout group or the healthy control group, with 30 participants in each group. The sample size was determined based on a previous study from our institution, showing that n = 23 patients with obstructive sleep apnea in each of two groups were required to detect a group difference (p < 0.05) in adenosine-induced hyperemic MBF of clinical significance (0.75 ± 0.85 mL/min/g) [21]. To our knowledge, there is no data available in the literature to perform a formal power analysis regarding CFR, which was the primary outcome variable of this study.

Table 1 Key parameter explanations

Participants were screened and assigned to the burnout or healthy control group using the MBI-Human Services Survey (MBI-HSS) [22] and the Patient Health Questionnaire (PHQ)-9 [23] via telephone interview, with additional inclusion and exclusion criteria reviewed. To determine the cutoff points for group assignment, we referred to a previous systematic review on physician burnout [2]. For the clinical burnout group, we used emotional exhaustion (EE) ≥ 27 and/or depersonalization (DP) ≥ 10 (with a minimum EE score ≥ 20) as cutoff points. For the control group, the cutoff points were EE < 16 and DP < 7. The personal accomplishment (PA) subscale was not used for group assignment as it develops largely independently of EE and DP [24]. In addition, for the burnout group, a PHQ-9 score of ≤ 14 was required, indicating at most moderate depressive symptoms, experienced stress at work, and undue exhaustion for at least 6 months prior to enrollment in the study [25]. For the control group, a PHQ-9 score of ≤ 10 was required, indicating at most mild depressive symptoms [23]. We sought to match controls with participants with burnout in terms of age in the range of 5 years, body mass index (BMI) in the range of 5 kg/m2, and family history of early CVD in first-degree relatives (< 55 years in men, < 65 years in women).

Both groups had additional inclusion criteria, which required non-smoking for at least 5 years and age between 28 and 65 years, as in Switzerland most physicians do not practice medicine for at least 6 months before age 28, and 65 is the official retirement age. Exclusion criteria for both groups included a previous episode of clinical depression or burnout; known heart disease; familial hypercholesterolemia; type I or type II diabetes; stage II hypertension; renal insufficiency; any known active serious disease; intake of lipid-lowering, antihypertensive, or antidiabetic drugs; BMI ≥ 35 kg/m2; chronic risky alcohol consumption (≥ 4 standard drinks per day); allergy to iodine-containing contrast media; contraindications for adenosine, beta-blockers, or nitrates; a medication that interferes with blood biomarker levels (e.g., corticosteroids, anticoagulants); and the decision to waive information about clinically relevant cardiac imaging findings.

Positron emission tomography for myocardial blood flow assessment

Participants were imaged with the latest generation PET-CT scanner (Discovery MI, GE Healthcare, Waukesha, WI, USA) using 13N-ammonia as a flow tracer at the Department of Nuclear Medicine, University Hospital of Zurich. We used adenosine infusion and the cold pressor test (CPT) as pharmacological and sympathetic stressors. MBF was assessed at rest, during adenosine-induced hyperemia (reflecting primarily endothelium-independent vasodilation), and in response to standardized CPT (reflecting primarily endothelium-dependent vasodilation).

Beginning with the intravenous bolus administration of 13N-ammonia, serial dynamic and static PET images were acquired for 20 minutes (min); 10 min later, adenosine (at a continuous rate of 140 μg/kg/min) was administered intravenously for 6 min. Three minutes into the adenosine infusion, a second dose of 13N-ammonia was injected and images were recorded in the same acquisition sequence. Ten minutes later, a CPT was performed by immersing the participant’s right foot and lower half of the calf in a container filled with ice water (4° C) for 2 min [26]. A third dose of 13N-ammonia was administered after 60 s into the CPT, and PET image acquisition was performed. Heart rate, blood pressure, and a 12-lead electrocardiogram were recorded at baseline and throughout the adenosine infusion and CPT. MBF at rest (corrected for rate-pressure product) and during adenosine-stress and CPT was calculated using commercially available software (QPET 2017.7, Cedars-Sinai Medical Center, Los Angeles, CA, USA). Each participant’s global CFR (primary outcome) was calculated as the ratio of stress to rest absolute MBF for the whole left ventricle. As secondary outcomes, we also examined MBF response (i.e., the difference between stress and rest MBF (ΔMBF)) that is not influenced by MBF at rest, and hyperemic MBF (i.e., the maximal flow achieved during hyperemia) [27].

Psychometric assessment

For the analyses, we used burnout and depressive symptom scores that were obtained on the same day as MBF measurements, because these may be more closely associated with MBF than scores obtained through the previous telephone interview. One participant lacked this information, so we used the symptom scores collected during the telephone screening.

Burnout

We used a validated German version of the 22-item MBI-HSS, which includes three subscales assessing EE (nine items), DP (five items), and PA (eight items) [22]. Participants rated each item on a scale ranging from 1 (“never”) to 7 (“daily”). While EE examines the sensation of exhaustion and energy depletion that arise from work, DP assesses a cynical and detached attitude towards work/patients, and PA measures feelings of competence and successful job performance as a physician. In our sample, Cronbach’s α was 0.95, 0.87, and 0.77 for the EE, DP, and PA subscales, respectively.

Job stress

To assess job stress, we used the validated German short form of the Effort-Reward Imbalance (ERI) and overcommitment (OC) questionnaire, consisting of three items regarding work effort, seven items concerning work-related rewards, and six items pertaining to OC to work [28]. The ERI and OC model posits that job stress with consequences for CVD [6, 7] stems from two factors: a lack of reciprocity in terms of high cost and low gain, reflected by the ER ratio (extrinsic component), and OC to high achievement (intrinsic component) [29]. Moreover, employees who report both a high level of extrinsic and intrinsic job stress are thought to be at the highest risk of CVD [30]. Each questionnaire item is rated on a 4-point Likert scale ranging from 1 (“strongly disagree”) to 4 (“strongly agree”). Higher scores on the effort scale (range 3–12) indicate a greater degree of stress caused by high work effort, while lower scores on the reward scale (range 7–28) reflect more stressful experiences due to inadequate rewards. Higher scores on the OC scale (range 6–24) indicate more stress caused by OC. The ER ratio is calculated using a correction factor that accounts for the unequal number of items of the effort and reward scores. A higher ER ratio indicates greater job stress. In our sample, Cronbach’s α was 0.76, 0.77, and 0.81 for the effort, reward, and OC scales, respectively.

Depressive symptoms

We used the German version of the PHQ-9 questionnaire to measure depressive symptoms in the previous 2 weeks [31]. The scale comprises nine items that are rated by participants on a 4-point Likert scale between 0 (“not at all”) and 3 (“nearly every day”). Higher total scores (range 0–27) indicate greater severity of depressive symptoms. Cronbach’s α for the PHQ-9 total scale was 0.79 in our sample.

Cardiovascular risk and health behavior assessment

We employed the Systematic Coronary Risk Evaluation (SCORE)2 risk prediction algorithm to estimate the 10-year risk of fatal and non-fatal CVD in individuals younger than 70 years without prior CVD or diabetes [32]. Body mass index (BMI) was calculated by dividing measured weight (in kilograms) by measured height (in meters) squared (kg/m2). Physical activity was assessed by asking participants how many times (range 0–7) in an average week they engaged in sports activities that caused them to sweat. Intake of alcohol was determined by the number of standard drinks consumed per week.

Data analysis

Data were analyzed using IBM SPSS Statistics for Windows, Version 29.0 (Armonk, NY: IBM Corp) with a two-tailed significance level of p < 0.05. The expectation–maximization algorithm was used to substitute the few missing items (3.4%) in the ERI and OC questionnaire. Because of the non-normal distribution, all measurements of microvascular function were transformed to statistical normality using a two-step procedure described elsewhere [33], preserving the original series mean and standard deviation. Independent samples t-test was used to compare health characteristics between the burnout group and the control group. Pearson correlation analysis was calculated to express zero-order correlations between two variables.

Linear regression analyses were modeled with CRF, ΔMBF, and hyperemic MBF as dependent variables and burnout (yes/no) and job stress as independent variables, adjusting for age and BMI. All independent variables were entered in one block. We included the “constant” (also known as the “intercept”) in all regression models, which represents the expected value of the dependent variable when all independent variables are set to zero. The inclusion of the constant ensures that the model is “unbiased” as the mean of the residuals (i.e., the differences between the observed and predicted values) is exactly zero, thereby contributing to the statistical rigor of the analysis. Adjustment for multiple comparisons was carried out for the multivariable associations of burnout and both job stress measures (ER ratio and OC) with CFR, which was the predefined primary outcome variable of our study. Since we assessed CFR using two different methods, significance was considered if the independent associations of burnout, ER ratio, or OC with endothelium-dependent CFR or endothelium-independent CFR was p < 0.025. Given that ER ratio and OC tap into distinct aspects of job stress (i.e., extrinsic vs. intrinsic component), we did not apply further p-value adjustment for these two measures. All other associations examined were exploratory and aimed to confirm and identify additional relationships between work-related variables and measures of coronary microvascular function, including ΔMBF and hyperemic MBF, supporting findings for CFR, and generating new hypotheses. For job stress, two models were calculated, one for the ER ratio and one for the OC score. Incorporating these models facilitates an assessment of result robustness across two distinct job stress measures, enhancing the statistical rigor of the analysis. When burnout (yes/no) or the ER ratio showed a significant independent association with a particular measure of microvascular function, additional analyses were performed to examine whether EE, DP, PA, effort, or reward would show separate associations with that same measure.

Two complementary analyses were conducted with interaction terms between burnout and job stress and between the ER ratio and OC, respectively, additionally included in the models. The first analysis examined if the relationship between job stress (either ER ratio or OC) and microvascular function would be moderated by the presence or absence of clinical burnout. The second analysis tested the theoretical assumption that physicians who have both high ER ratio (extrinsic work stress) and OC scores (intrinsic work stress) have a greater reduction in microvascular function.

Due to high multicollinearity between the PHQ-9 score and burnout (by study design), and between age and the SCORE2, depressive symptoms and the SCORE2 were not considered as covariates to prevent inaccurate coefficient estimates. Specifically, the Pearson correlation coefficient yielded a value of 0.711 when examining the relationship between burnout and the PHQ-9 score. Additionally, when included in the multivariable model, the PHQ-9 score exhibited a variance inflation factor (VIF) of 2.761. We adjusted for age instead of SCORE2 because, despite our attempts to match groups on age, the burnout group was significantly younger, whereas both groups had a similar SCORE2 that already accounts for age. All VIFs were below 2.5 in multivariable models, indicating no concern for multicollinearity. For instance, in their respective multivariable models, VIFs were 1.642 for burnout and 1.760 for the ER ratio, and 1.602 for burnout and 1.620 for OC. One case was omitted for all models with the ER ratio because Mahalanobis distance indicated an influential outlier. Effect sizes of associations are expressed as correlation coefficients, where r = 0.1, 0.3, and 0.5 reflect small, medium, and large effects, respectively.

Results

Participant characteristics

Table 2 shows that physicians with burnout were significantly younger than their counterparts without burnout (p = 0.012), whereas the other cardiovascular risk factors studied did not show a significant group difference (p-values > 0.05). In contrast, psychometric scores differed significantly between the two groups. Physicians with burnout reported higher levels of job stress in terms of both higher ER ratio (p = 2 × 10−6) and overcommitment (p = 2 × 10−6). Due to the study design, physicians with burnout also reported higher levels of burnout symptoms (p-values ≤ 7 × 10−5) and depressive symptoms (p = 5 × 10−10) than those without burnout.

Table 2 Characteristics of the 60 study participants

Myocardial blood flow at rest

Table 3 shows zero-order correlations with MBF at rest. Both higher CVD risk (r = 0.270, p = 0.037) and higher reward (r = 0.314, p = 0.014) were significantly associated with higher MBF at rest, whereas higher job stress in terms of a higher ER ratio was significantly associated with lower MBF at rest (r =  − 0.335, p = 0.009). In regression analysis with age, BMI, burnout, and ER ratio in the model, ER ratio showed an independent association with MBF at rest (unstandardized coefficient B =  − 0.108, 95% confidence interval (CI) =  − 0.215, − 0.001; p = 0.048); this was not observed for burnout (p = 0.18), age (p = 0.27), and BMI (p = 0.27). Neither burnout (p = 0.74) nor OC (p = 0.89) was associated with MBF at rest in regression analysis with age (p = 0.084), BMI (p = 0.15), burnout, and OC as independent variables.

Table 3 Zero-order correlations between variables of interest (n = 60)

Endothelium-dependent coronary microvascular function

There were several significant associations of burnout and job stress with CFR, ΔMBF, and hyperemic MBF in the CPT.

Zero-order correlations

As can be seen in Table 3, there were significant inverse correlations between job stress and measures of microvascular function of mostly medium effect size (r between 0.27 and 0.41). A higher ER-ratio was significantly correlated with lower CRF (r =  − 0.270, p = 0.037), lower ΔMBF (r =  − 0.322, p = 0.012), and lower hyperemic MBF (r =  − 0.408, p = 0.001). Likewise, higher OC was significantly correlated with lower CRF (r =  − 0.261, p = 0.044), lower ΔMBF (r =  − 0.270, p = 0.037), and lower hyperemic MBF (r =  − 0.298, p = 0.021). Of the two ERI components, higher efforts showed a significant association with lower microvascular function for all three measures (p-values < 0.023), but reward did not (p-values > 0.083). Burnout (yes/no) showed no significant correlation with any measure of microvascular function (p-values > 0.58).

Independent associations

As shown in Table 4 and Fig. 1, burnout (yes/no) was significantly and independently associated with higher CRF (r partial = 0.325, p = 0.014), higher ΔMBF (r partial = 0.351, p = 0.008), and higher hyperemic MBF (r partial = 0.334, p = 0.012), controlling for age, BMI, and the ER ratio. Likewise, burnout was significantly and independently associated with higher CRF (r partial = 0.345, p = 0.009), higher ΔMBF (r partial = 0.332, p = 0.012), and higher hyperemic MBF (r partial = 0.281, p = 0.035), controlling for age, BMI, and OC.

Table 4 Associations of burnout and job stress with endothelium-dependent microvascular function (cold pressor test)
Fig. 1
figure 1

The individual data points and mean values (95% confidence interval) of the primary outcome coronary flow reserve (CFR) and the secondary outcomes myocardial blood flow (MBF) response and hyperemic MBF in response to the cold pressure test between burnout and control physicians. All analyses were adjusted for age and body mass index and either the effort-reward ratio (A, C, and E) or overcommitment scores (B, D, and F). The figure was generated using R version 4.2.0

In contrast, as shown in Table 4 and Fig. 2, higher job stress, measured by the ER ratio, was significantly and independently associated with lower CRF (r partial =  − 0.322, p = 0.015), lower ΔMBF (r partial =  − 0.368, p = 0.005), and lower hyperemic MBF (r partial =  − 0.381, p = 0.004). Job stress measured by the OC score was also significantly and independently associated with lower CRF (r partial =  − 0.366, p = 0.005), lower ΔMBF (r partial =  − 0.346, p = 0.008), and lower hyperemic MBF (r partial =  − 0.303, p = 0.022). The majority of adjusted correlation coefficients exceeded 0.3, with these independent associations indicating a medium effect size.

Fig. 2
figure 2

Partial regression plots with fit line (95% confidence interval) for the associations between the effort-reward ratio (A, C, and E) and overcommitment scores (B, D, and F) with the primary outcome coronary flow reserve (CFR) and the secondary outcomes myocardial blood flow (MBF) response and hyperemic MBF in response to the cold pressure test among all participants. All analyses were adjusted for age, body mass index, and burnout (yes/no)

Follow-up analyses with individual burnout dimension revealed that higher EE was significantly and independently associated with greater microvascular function regarding all three measures (r partial between 0.280 and 0.329, p-values between 0.013 and 0.037). These associations held true after adjusting for age, BMI, and either the ER ratio (Table 5, Model 1) or OC score (Table 5, Model 2). Similar significant and independent associations with greater microvascular function were observed for higher DP (r partial between 0.261 and 0.317, p-values between 0.016 and 0.049), except for the association between DP and CFR (p = 0.051) (Table 5, Model 1); PA showed no significant association with microvascular function measures (p-values > 0.09). Of the ERI components, higher efforts showed significant and independent associations with lower microvascular function with all three measures (Table 5, Models 3 and 4: r partial between − 0.414 and − 0.372, p-values between 0.001 and 0.005), but reward did not (p-values > 0.32).

Table 5 Associations of burnout and job stress dimensions with endothelium-dependent microvascular function (cold pressor test)

Endothelium-independent microvascular function

With the adenosine challenge, one participant in the burnout group, but none in the control group, had CFR < 2. Zero-order correlations (Table 3) showed that higher job stress, measured by the ER ratio, was significantly correlated with lower ΔMBF (r =  − 0.323, p = 0.012) and lower hyperemic MBF (r =  − 0.406, p = 0.001). Similarly, a higher OC score was also correlated with lower ΔMBF (r =  − 0.270, p = 0.037) and lower hyperemic MBF (r =  − 0.299, p = 0.020). Of the ERI components, the effort, but not the reward scale, showed a significant correlation with ΔMBF (r =  − 0.351, p = 0.006) and hyperemic MBF (r =  − 0.411, p = 0.001). These correlations were mostly of medium effect size. However, as can be seen in Table 6, none of the significant zero-order correlations maintained significance in multivariable regression models (p-values > 0.05). Burnout (yes/no) showed no significant independent association with any measure of endothelium-independent coronary microvascular function (p-values > 0.69).

Table 6 Associations of burnout and job stress with endothelium-independent microvascular function (adenosine challenge)

Complementary analyses

The interaction between the ER ratio and OC showed no independent association with any measure of microvascular function (p-values > 0.16). Likewise, the interaction between burnout (yes/no) and the ER ratio and the interaction between burnout (yes/no) and OC were not significant for any measure of microvascular function (p-values > 0.20).

Discussion

We found that both burnout and job stress showed significant and independent, albeit opposing, associations with endothelium-dependent coronary microvascular function in male physicians. Specifically, burnout exhibited positive associations with endothelium-dependent CFR, the primary outcome variable in our study, MBF response, and hyperemic MBF. In turn, job stress, defined by the ER ratio and OC, showed negative associations with endothelium-dependent CFR, MBF response, and hyperemic MBF. The magnitude of the observed effects suggests that these associations may be of clinical relevance. When examining the individual dimensions of burnout and job stress, EE, DP, and efforts were significantly and independently associated with coronary microvascular endothelial function. To put it simply, coronary microvascular endothelial function was better in male physicians versus those without burnout, but reduced in those with higher job stress.

To our knowledge, endothelial function has not previously been investigated in burnout, whereas, consistent with our findings, paramedics with high job strain were shown to have endothelial dysfunction measured by peripheral artery tonometry [34]. Our finding of and association of endothelial-dependent coronary microvascular dysfunction with higher ER ratio or OC aligns with meta-analyses showing a link between various aspects of job stress, including ERI, and an increased CHD risk [10, 14, 35,36,37]. While the ER ratio has been extensively studied in relation to CHD risk, OC has received comparatively less attention [12, 38]. However, OC has been associated with subclinical atherosclerosis in Chinese workers [39] and an increased CHD risk in patients referred for coronary angiography [40] and restenosis in men following coronary angioplasty [41]. We therefore interpret that job stress resulting both extrinsically from the work environment and intrinsically from an exhaustive pattern of coping with demanding work-related situations could be associated with poor cardiovascular health in male physicians. High effort was significantly associated with microvascular endothelial function whereas low reward was not. In contrast, we found no evidence for the theoretical concept [30] that the association between higher ER ratio and lower microvascular endothelial function would be strongest in male physicians scoring high on OC. This concurs with the majority of empirical investigations in which OC was not found to be a moderator of the association between ERI and various health outcomes [42].

Regarding clinical implications, effective psychosocial interventions to manage job stress in physicians [43] could be tested for a potential benefit on microvascular endothelial dysfunction. Such research may be informed by a previous trial showing that 12 weeks of meditation practice improved PET-assessed CFR in CHD patients [44]. Another study in patients with stable CHD showed that the combination of usual care and weekly 1.5-h stress management sessions for 16 weeks resulted in greater improvement in endothelial function (i.e., flow-mediated dilation) than usual care alone [45].

As burnout has been associated with an increased risk of CHD [11], the observation that burnout was related to increased endothelium-dependent coronary microvascular function seems counterintuitive. Although speculative, we offer some biobehavioral explanations for interpreting the finding of apparently improved endothelium-dependent coronary microvascular function in the burnout group. Based on a previous study in which male patients with burnout exhibited higher systolic blood pressure during psychosocial stress than healthy male controls [46], there could be an increased sympathetic reaction to the CPT in our male physicians with burnout. This hypothesis gains support from the observation that resting MBF was not different between physicians with burnout and those without, yet significantly lower in physicians with a higher ER ratio. At least during the acute phase, the burnout syndrome could possibly have a protective effect on cardiovascular health against high levels of stress via inhibiting even higher levels of stress and effort. Specifically, sympathetic activity has been associated with increased low-grade inflammation [47] and there is a sequence leading first to endothelial proinflammatory activation, which precedes endothelial dysfunction in atherogenesis [48]. Therefore, systemic inflammation in burnout [15] due to sympathetic activity [46] might be associated with endothelial activation rather than dysfunction in our participants. Longitudinal studies could investigate whether the positive relationship between burnout and coronary microvascular endothelial function changes to a negative one with a longer duration of burnout symptomatology. This relationship could be further explored, particularly with a focus on investigating the role of the sympathetic nervous system and low-grade inflammation. This may particularly be the case for individuals with high EE and high DP, both of which were associated with microvascular function in our study and an increased CVD risk in a large population from Finland [49].

In statistical terms, the observation that the effect of burnout on coronary microvascular endothelial function only emerged in the fully adjusted model suggests a suppressor effect. The inclusion of job stress likely improved the model’s predictive power by accounting for the confounding relationship between burnout and job stress, thereby revealing a “true” association between burnout and microvascular function. We found no evidence that burnout modified the association between job stress and microvascular endothelial function. By all means, examining the independent associations of burnout and job stress in the same study could provide a more comprehensive understanding of the complex relationships between work-related stressors and workers’ cardiovascular health. It should be noted, however, that CFR and hyperemic MBF and MBF response without adjustment for covariates did not show significant group differences. This could be due to the small sample size or to the fact that burnout does not really affect coronary microvascular function.

We found no significant independent association of burnout and job stress with endothelial-independent coronary microvascular function as assessed through an adenosine challenge. The selective impairment of endothelium-dependent regulation of CFR, MBF response, and hyperemic MBF indicates a specific dysfunction in the endothelium’s ability to respond to sympathetic stimuli (e.g., the CPT used in our study), and release vasodilators, like nitric oxide. The endothelium plays a crucial role in the mental stress-induced coronary microvascular response [50], and impaired endothelium-dependent coronary vasodilation is linked to exercise-induced myocardial ischemia [51]. Therefore, in individuals with high job stress, there could be increased vulnerability to acute coronary events when there is a temporary rise in sympathetic activity due to, for instance, acute emotional stress or intense physical exertion [52,53,54].

While we made efforts to match physicians with burnout to controls by age, it is worth noting that the burnout group turned out to be younger, consistent with findings from large population-based studies showing an inverse association between burnout and age [55, 56]. This age difference may be partly explained by the “healthy worker effect” [55, 56], a phenomenon observed when employees with burnout leave their jobs or professions at a younger age than those who remain in the workforce without health issues. It is possible that older physicians have developed better strategies for coping with job-related stress, or there may be differences in training among different cohorts of physicians.

The use of PET to quantify microvascular function with established techniques and validated questionnaires to limit measurement bias, and the inclusion of a well-defined group of male physicians with and without clinical burnout were strengths of our study. However, this study also has several noteworthy limitations. The calculation of the target sample size was not based on power analysis for the primary study outcome. The criteria for matching between the burnout and control groups, along with the exclusion criteria related to numerous health conditions, potentially affecting both cardiovascular health and workplace-related variables, aimed to minimize residual confounding. However, the sample size allowed only a limited number of covariates to be included in the multivariable analyses, so residual confounding cannot be excluded, resulting in some level of uncertainty in our findings. For instance, the inclusion of inflammatory biomarkers could have provided additional valuable information. Therefore, the study was likely underpowered and the results should be considered hypothesis-generating. Moreover, while we made adjustment for multiple comparisons for the primary outcome CFR, the additional examination of multiple exposures and multiple outcomes in our study was exploratory and raises the possibility that some p-values became significant by chance alone, so these results should be interpreted with caution. The study employed a cross-sectional design, aimed to provide a snapshot of current relationships between work-related variables and microvascular function, facilitating the identification of patterns of association. However, the deliberate choice of this design limits the possibility of inferring causality or mechanisms and ruling out reverse causality. Existing literature and conceptual considerations substantiate the idea that job stress and burnout are risk factors for CVD development. Notably, our study participants were free of CVD and reported experiencing undue exhaustion and job stress for a minimum of 6 months before the study. While these inclusion criteria may alleviate concerns about reverse causality, it does not preclude the possibility that subclinical CVD contributed to exhaustion, requiring increased efforts to maintain job performance. To enhance the robustness of our findings, we emphasize the need for future longitudinal studies to examine the temporal relationships between burnout, job stress, and coronary microvascular function. The specificity of the study population increased the study’s internal validity, but the numerous exclusion criteria limited the sample size affecting the generalizability of the results. We only included male physicians without a history of CVD, so it is unclear whether the findings generalize to other occupational groups, patients with established CVD, and particularly female physicians. This is important because, in response to acute stress, pre-menopausal women show both smaller autonomic responses [57] and greater endothelial dysfunction [58] than men of the same age.

Conclusions

The present study provides some evidence that burnout and job stress show opposite associations with coronary microvascular endothelial function in male physicians. Longitudinal studies with larger sample sizes are preferably needed to determine the transferability of our findings to other populations, including female physicians, and their clinical relevance for targeted interventions to reduce CHD risk in the working population.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Change history

Abbreviations

BMI:

Body mass index

CFR:

Coronary flow reserve

CHD:

Coronary heart disease

CPT:

Cold pressor test

CVD:

Cardiovascular disease

DP:

Depersonalization

EE:

Emotional exhaustion

ERI:

Effort-reward imbalance

MBF:

Myocardial blood flow

MBI-HSS:

Maslach Burnout Inventory-Human Services Survey

OC:

Overcommitment

PA:

Personal accomplishment

PET-CT:

Positron emission tomography computed tomography

PHQ:

Patient Health Questionnaire

SCORE:

Systematic Coronary Risk Evaluation

VIF:

Variance inflation factor

References

  1. WHO. ICD-11: International Classification of Diseases (11th Revision). https://icd.who.int/browse11/l-m/en#/http://id.who.int/icd/entity/129180281. Page last Accessed 18 Nov 2023.

  2. Shoman Y, El May E, Marca SC, Wild P, Bianchi R, Bugge MD, et al. Predictors of occupational burnout: a systematic review. Int J Environ Res Public Health. 2021;18(17):9188.

    PubMed  PubMed Central  Google Scholar 

  3. Rotenstein LS, Torre M, Ramos MA, Rosales RC, Guille C, Sen S, et al. Prevalence of burnout among physicians: a systematic review. JAMA. 2018;320(11):1131–50.

    PubMed  PubMed Central  Google Scholar 

  4. Fuss I, Nübling M, Hasselhorn HM, Schwappach D, Rieger MA. Working conditions and work-family conflict in German hospital physicians: psychosocial and organisational predictors and consequences. BMC Public Health. 2008;8:353.

    PubMed  PubMed Central  Google Scholar 

  5. Le Huu P, Bellagamba G, Bouhadfane M, Villa A, Lehucher MP. Meta-analysis of effort-reward imbalance prevalence among physicians. Int Arch Occup Environ Health. 2022;95(3):559–71.

    PubMed  Google Scholar 

  6. West CP, Dyrbye LN, Shanafelt TD. Physician burnout: contributors, consequences and solutions. J Intern Med. 2018;283(6):516–29.

    CAS  PubMed  Google Scholar 

  7. Tsutsumi A, Kawanami S, Horie S. Effort-reward imbalance and depression among private practice physicians. Int Arch Occup Environ Health. 2012;85(2):153–61.

    PubMed  Google Scholar 

  8. Loerbroks A, Weigl M, Li J, Angerer P. Effort-reward imbalance and perceived quality of patient care: a cross-sectional study among physicians in Germany. BMC Public Health. 2016;16:342.

    PubMed  PubMed Central  Google Scholar 

  9. Salvagioni DAJ, Melanda FN, Mesas AE, González AD, Gabani FL, Andrade SM. Physical, psychological and occupational consequences of job burnout: a systematic review of prospective studies. PLoS ONE. 2017;12(10): e0185781.

    PubMed  PubMed Central  Google Scholar 

  10. Ganster DC, Rosen CC. Work stress and employee health: a multidisciplinary review. J Manage. 2013;39(5):1085–122.

    Google Scholar 

  11. Toker S, Melamed S, Berliner S, Zeltser D, Shapira I. Burnout and risk of coronary heart disease: a prospective study of 8838 employees. Psychosom Med. 2012;74(8):840–7.

    PubMed  Google Scholar 

  12. Dragano N, Siegrist J, Nyberg ST, Lunau T, Fransson EI, Alfredsson L, et al. Effort-reward imbalance at work and incident coronary heart disease: a multicohort study of 90,164 individuals. Epidemiology. 2017;28(4):619–26.

    PubMed  PubMed Central  Google Scholar 

  13. Hegde SK, Vijayakrishnan G, Sasankh AK, Venkateswaran S, Parasuraman G. Lifestyle-associated risk for cardiovascular diseases among doctors and nurses working in a medical college hospital in Tamil Nadu, India. J Fam Med Prim Care. 2016;5(2):281–5.

    Google Scholar 

  14. Jardim TV, Sousa AL, Povoa TR, Barroso WS, Chinem B, Jardim PC. Comparison of cardiovascular risk factors in different areas of health care over a 20-year period. Arq Bras Cardiol. 2014;103(6):493–501.

    PubMed  PubMed Central  Google Scholar 

  15. Bayes A, Tavella G, Parker G. The biology of burnout: causes and consequences. World J Biol Psychiatry. 2021;22(9):686–98.

    PubMed  Google Scholar 

  16. Kivimäki M, Steptoe A. Effects of stress on the development and progression of cardiovascular disease. Nat Rev Cardiol. 2018;15(4):215–29.

    PubMed  Google Scholar 

  17. Chida Y, Steptoe A. Cortisol awakening response and psychosocial factors: a systematic review and meta-analysis. Biol Psychol. 2009;80(3):265–78.

    PubMed  Google Scholar 

  18. Vancheri F, Longo G, Vancheri S, Henein M. Coronary Microvascular Dysfunction. J Clin Med. 2020;9(9):2880.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Ziadi MC. Myocardial flow reserve (MFR) with positron emission tomography (PET)/computed tomography (CT): clinical impact in diagnosis and prognosis. Cardiovasc Diagn Ther. 2017;7(2):206–18.

    PubMed  PubMed Central  Google Scholar 

  20. Kelshiker MA, Seligman H, Howard JP, Rahman H, Foley M, Nowbar AN, et al. Coronary flow reserve and cardiovascular outcomes: a systematic review and meta-analysis. Eur Heart J. 2022;43(16):1582–93.

    PubMed  Google Scholar 

  21. Schwarz EI, Schlatzer C, Stehli J, Kaufmann PA, Bloch KE, Stradling JR, et al. Effect of CPAP Withdrawal on myocardial perfusion in OSA: a randomized controlled trial. Respirology. 2016;21(6):1126–33.

    PubMed  Google Scholar 

  22. Büssing A, Perrar KM. Die Messung von Burnout. Untersuchung einer deutschen Fassung des Maslach Burnout Inventory (MBI-D). [Measuring burnout: a study of a German version of the Maslach Burnout Inventory (MBI-D)]. Diagnostica. 1992;38(4):328–53.

    Google Scholar 

  23. Kroenke K, Spitzer RL, Williams JB. The PHQ-9: validity of a brief depression severity measure. J Gen Intern Med. 2001;16(9):606–13.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Lee RT, Ashforth BE. A meta-analytic examination of the correlates of the three dimensions of job burnout. J Appl Psychol. 1996;81(2):123–33.

    CAS  PubMed  Google Scholar 

  25. Jönsson P, Österberg K, Wallergård M, Hansen ÅM, Garde AH, Johansson G, et al. Exhaustion-related changes in cardiovascular and cortisol reactivity to acute psychosocial stress. Physiol Behav. 2015;151:327–37.

    PubMed  Google Scholar 

  26. Siegrist PT, Gaemperli O, Koepfli P, Schepis T, Namdar M, Valenta I, et al. Repeatability of cold pressor test-induced flow increase assessed with H(2)(15)O and PET. J Nucl Med. 2006;47(9):1420–6.

    PubMed  Google Scholar 

  27. Prior JO, Schindler TH, Facta AD, Hernandez-Pampaloni M, Campisi R, Dahlbom M, et al. Determinants of myocardial blood flow response to cold pressor testing and pharmacologic vasodilation in healthy humans. Eur J Nucl Med Mol Imaging. 2007;34(1):20–7.

    PubMed  Google Scholar 

  28. Siegrist J, Wege N, Pühlhofer F, Wahrendorf M. A short generic measure of work stress in the era of globalization: effort-reward imbalance. Int Arch Occup Environ Health. 2009;82(8):1005–13.

    PubMed  Google Scholar 

  29. Siegrist J. Effort-reward imbalance at work and cardiovascular diseases. Int J Occup Med Environ Health. 2010;23(3):279–85.

    MathSciNet  PubMed  Google Scholar 

  30. van Vegchel N, de Jonge J, Bosma H, Schaufeli W. Reviewing the effort-reward imbalance model: drawing up the balance of 45 empirical studies. Soc Sci Med. 2005;60(5):1117–31.

    PubMed  Google Scholar 

  31. Gräfe K, Zipfel S, Herzog W, Löwe B. Screening psychischer Störungen mit dem “Gesundheitsfragebogen für Patienten (PHQ-D)“. Diagnostica. 2004;50(4):171–81.

    Google Scholar 

  32. SCORE2 working group and ESC Cardiovascular risk collaboration. SCORE2 risk prediction algorithms: new models to estimate 10-year risk of cardiovascular disease in Europe. Eur Heart J. 2021;42(25):2439–54.

    Google Scholar 

  33. Templeton GF. A two-step approach for transforming continuous variables to normal: implications and recommendations for IS research. Commun Assoc Inf Syst. 2011;28(1):4.

    Google Scholar 

  34. Wong IS, Ostry AS, Demers PA, Davies HW. Job strain and shift work influences on biomarkers and subclinical heart disease indicators: a pilot study. J Occup Environ Hyg. 2012;9(8):467–77.

    PubMed  Google Scholar 

  35. Kivimäki M, Nyberg ST, Batty GD, Fransson EI, Heikkilä K, Alfredsson L, et al. Job strain as a risk factor for coronary heart disease: a collaborative meta-analysis of individual participant data. Lancet. 2012;380(9852):1491–7.

    PubMed  PubMed Central  Google Scholar 

  36. Virtanen M, Nyberg ST, Batty GD, Jokela M, Heikkilä K, Fransson EI, et al. Perceived job insecurity as a risk factor for incident coronary heart disease: systematic review and meta-analysis. BMJ. 2013;347: f4746.

    PubMed  PubMed Central  Google Scholar 

  37. Taouk Y, Spittal MJ, LaMontagne AD, Milner AJ. Psychosocial work stressors and risk of all-cause and coronary heart disease mortality: a systematic review and meta-analysis. Scand J Work Environ Health. 2020;46(1):19–31.

    CAS  PubMed  Google Scholar 

  38. Eddy P, Wertheim EH, Kingsley M, Wright BJ. Associations between the effort-reward imbalance model of workplace stress and indices of cardiovascular health: a systematic review and meta-analysis. Neurosci Biobehav Rev. 2017;83:252–66.

    PubMed  Google Scholar 

  39. Xu W, Zhao Y, Guo L, Guo Y, Gao W. Job stress and coronary heart disease: a case-control study using a Chinese population. J Occup Health. 2009;51(2):107–13.

    PubMed  Google Scholar 

  40. Xu W, Hang J, Cao T, Shi R, Zeng W, Deng Y, et al. Job stress and carotid intima-media thickness in Chinese workers. J Occup Health. 2010;52(5):257–62.

    PubMed  Google Scholar 

  41. Joksimovic L, Siegrist J, Meyer-Hammer M, Peter R, Franke B, Klimek WJ, et al. Overcommitment predicts restenosis after coronary angioplasty in cardiac patients. Int J Behav Med. 1999;6(4):356–69.

    CAS  PubMed  Google Scholar 

  42. Siegrist J, Li J. Associations of extrinsic and intrinsic components of work stress with health: a systematic review of evidence on the effort-reward imbalance model. Int J Environ Res Public Health. 2016;13(4):432.

    PubMed  PubMed Central  Google Scholar 

  43. Clough BA, March S, Chan RJ, Casey LM, Phillips R, Ireland MJ. Psychosocial interventions for managing occupational stress and burnout among medical doctors: a systematic review. Syst Rev. 2017;6(1):144.

    PubMed  PubMed Central  Google Scholar 

  44. Bokhari S, Schneider RH, Salerno JW, Rainforth MV, Gaylord-King C, Nidich SI. Effects of cardiac rehabilitation with and without meditation on myocardial blood flow using quantitative positron emission tomography: a pilot study. J Nucl Cardiol. 2021;28(4):1596–607.

    PubMed  Google Scholar 

  45. Blumenthal JA, Sherwood A, Babyak MA, Watkins LL, Waugh R, Georgiades A, et al. Effects of exercise and stress management training on markers of cardiovascular risk in patients with ischemic heart disease: a randomized controlled trial. JAMA. 2005;293(13):1626–34.

    CAS  PubMed  Google Scholar 

  46. de Vente W, van Amsterdam JG, Olff M, Kamphuis JH, Emmelkamp PM. Burnout is associated with reduced parasympathetic activity and reduced HPA axis responsiveness, predominantly in males. Biomed Res Int. 2015;2015: 431725.

    PubMed  PubMed Central  Google Scholar 

  47. Tawakol A, Ishai A, Takx RA, Figueroa AL, Ali A, Kaiser Y, et al. Relation between resting amygdalar activity and cardiovascular events: a longitudinal and cohort study. Lancet. 2017;389(10071):834–45.

    PubMed  PubMed Central  Google Scholar 

  48. Zhang J. Biomarkers of endothelial activation and dysfunction in cardiovascular diseases. Rev Cardiovasc Med. 2022;23(2):73.

    MathSciNet  PubMed  Google Scholar 

  49. Kajantie E, Phillips DI. The effects of sex and hormonal status on the physiological response to acute psychosocial stress. Psychoneuroendocrinology. 2006;31(2):151–78.

    CAS  PubMed  Google Scholar 

  50. Honkonen T, Ahola K, Pertovaara M, Isometsä E, Kalimo R, Nykyri E, et al. The association between burnout and physical illness in the general population–results from the Finnish Health 2000 Study. J Psychosom Res. 2006;61(1):59–66.

    PubMed  Google Scholar 

  51. Hammadah M, Kim JH, Al Mheid I, Samman Tahhan A, Wilmot K, Ramadan R, et al. Coronary and peripheral vasomotor responses to mental stress. J Am Heart Assoc. 2018;7(10): e008532.

    PubMed  PubMed Central  Google Scholar 

  52. Zeiher AM, Krause T, Schächinger V, Minners J, Moser E. Impaired endothelium-dependent vasodilation of coronary resistance vessels is associated with exercise-induced myocardial ischemia. Circulation. 1995;91(9):2345–52.

    CAS  PubMed  Google Scholar 

  53. Mostofsky E, Penner EA, Mittleman MA. Outbursts of anger as a trigger of acute cardiovascular events: a systematic review and meta-analysis. Eur Heart J. 2014;35(21):1404–10.

    PubMed  PubMed Central  Google Scholar 

  54. Kang DO, Eo JS, Park EJ, Nam HS, Song JW, Park YH, et al. Stress-associated neurobiological activity is linked with acute plaque instability via enhanced macrophage activity: a prospective serial 18F-FDG-PET/CT imaging assessment. Eur Heart J. 2021;42(19):1883–95.

    CAS  PubMed  Google Scholar 

  55. Ahola K, Honkonen T, Virtanen M, Aromaa A, Lönnqvist J. Burnout in relation to age in the adult working population. J Occup Health. 2008;50(4):362–5.

    PubMed  Google Scholar 

  56. Cunradi CB, Chen MJ, Lipton R. Association of occupational and substance use factors with burnout among urban transit operators. J Urban Health. 2009;86(4):562–70.

    PubMed  PubMed Central  Google Scholar 

  57. Willich SN, Lewis M, Löwel H, Arntz HR, Schubert F, Schröder R. Physical exertion as a trigger of acute myocardial infarction. Triggers and Mechanisms of Myocardial Infarction Study Group. N Engl J Med. 1993;329(23):1684–90.

    CAS  PubMed  Google Scholar 

  58. Martin EA, Tan SL, MacBride LR, Lavi S, Lerman LO, Lerman A. Sex differences in vascular and endothelial responses to acute mental stress. Clin Auton Res. 2008;18(6):339–45.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors express their gratitude to the nuclear medicine technicians, Josephine Trinckauf, Corina Weyermann, and their entire team, for the excellent technical support.

Funding

The study was financially supported by an institutional grant from the University of Zurich, Switzerland, to RvK.

Author information

Authors and Affiliations

Authors

Contributions

RvK, APP, MP, PAK, and RBB conceptualized, designed, and supervised the study. MP, SAH, DCB, AAG, RRB, CZH, and APP collected the data. RvK, AR, and APP performed the data analysis. RvK, SAH, and APP wrote the manuscript. MP, SAH, CG, DCB, AAG, RRB, CZH, and APP were responsible for recruiting patients and recording clinical data. All authors revised the manuscript for important intellectual contents and read and approved the final manuscript for submission.

Corresponding author

Correspondence to Roland von Känel.

Ethics declarations

Ethics approval and consent to participate

The local ethics committee Zurich (BASEC-Nr. 2018–01974) approved the study, and written informed consent was obtained from all participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

von Känel, R., Princip, M., Holzgang, S.A. et al. Coronary microvascular function in male physicians with burnout and job stress: an observational study. BMC Med 21, 477 (2023). https://doi.org/10.1186/s12916-023-03192-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-023-03192-z

Keywords