ORIGINAL RESEARCH article

Front. Public Health, 16 June 2025

Sec. Occupational Health and Safety

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1595550

Abnormal intrinsic functional hubs and connectivity in nurses with occupational burnout: a resting-state fMRI study

Jian-Ping Liu,&#x;Jian-Ping Liu1,2Si-Yu Gu&#x;Si-Yu Gu3Chun-Mei Song&#x;Chun-Mei Song4Hu-Cheng YangHu-Cheng Yang5Yang ShiYang Shi6Yu-Fang GuYu-Fang Gu1Shu-Fang Wang
Shu-Fang Wang2*Ying-Zhu Chen
Ying-Zhu Chen7*
  • 1School of Nursing and Public Health, Yangzhou University, Yangzhou, China
  • 2Department of Neurology, Yancheng Clinical Medical College of Yangzhou University, Yancheng Third People’s Hospital, Yancheng, China
  • 3Department of Radiology, Affiliated Hospital 6 of Nantong University, Yancheng Third People’s Hospital, Yancheng, China
  • 4Department of Disinfection Supply Center, Affiliated Hospital 6 of Nantong University, Yancheng Third People’s Hospital, Yancheng, China
  • 5Department of Radiology, Binhai Maternal and Child Health Hospital, Yancheng, China
  • 6Department of Neurology, Yangzhou Wutaishan Hospital of Jiangsu Province, Teaching Hospital of Yangzhou University, Yangzhou, China
  • 7Department of Geriatrics, Northern Jiangsu People’s Hospital Affiliated to Yangzhou University, Yangzhou, China

Background: Occupational burnout is a significant problem among nurses, linked to negative outcomes. Understanding its neurobiological basis is crucial, yet remains limited.

Methods: Resting-state functional magnetic resonance imaging (rs-fMRI) data were acquired from 40 female nurses with occupational burnout and 40 healthy controls. Degree centrality (DC) was calculated to identify functional hubs, and subsequent functional connectivity (FC) analysis was performed. Group differences in DC and FC were statistically compared. Their correlations with Maslach Burnout Inventory-Human Services Survey (MBI-HSS) scores were assessed, and a classification model was built using DC and FC features to distinguish between burnout and control groups.

Results: The burnout group showed significantly decreased DC in bilateral precuneus and reduced FC between left precuneus and right medial orbitofrontal cortex (mOFC) compared to the healthy control group. These neuroimaging markers correlated with clinical burnout dimensions: precuneus DC negatively associated with emotional exhaustion and depersonalization, while precuneus-mOFC connectivity positively correlated with personal accomplishment. A linear discriminant analysis model combining DC and FC measures achieved 85% classification accuracy (sensitivity 80%, specificity 90%) in distinguishing burnout from controls.

Conclusion: These findings identify the precuneus and its mOFC connectivity as key neural substrates of occupational burnout, suggesting disrupted integration between self-referential processing and reward/emotion regulation systems. Our results advance understanding of burnout’s neurobiological mechanisms and demonstrate the potential of neuroimaging markers for objective burnout assessment.

1 Introduction

Burnout is defined as a psychological state arising from chronic emotional or interpersonal stressors in the workplace (1). It manifests as an emotional and behavioral disorder caused by prolonged exposure to high occupational stress, particularly prevalent among nurses in the healthcare sector (2). Globally, approximately 30% of nurses experience occupational burnout (3), leading to increased turnover rates (4), economic losses (5), and compromised patient safety (6). The COVID-19 pandemic has exacerbated this crisis, intensifying psychological strain and accelerating burnout rates (7). Despite its severe societal and clinical implications, the neurobiological mechanisms underlying occupational burnout in nurse remain poorly understood, hindering early diagnosis and targeted interventions.

Over the past two decades, resting-state functional magnetic resonance imaging (rs-fMRI) has emerged as a powerful tool for investigating brain function alterations. By measuring blood oxygen level-dependent (BOLD) signals, rs-fMRI reveals spontaneous brain activity and the synchronization of neural networks in the absence of specific tasks (8, 9). Functional connectivity (FC), which assesses the temporal correlation of neural signals between brain regions, has been widely used to explore the intrinsic interactions within brain networks (10, 11). Previous studies have applied FC to investigate occupational burnout, identifying disrupted connectivity patterns (1215); however, these efforts often rely on seed-based regions selected from prior literature, limiting their ability to capture the full scope of burnout-related neurobiological changes. Degree centrality (DC), a graph-based metric, offers an unbiased approach by quantifying the importance of brain regions as functional hubs within the connectome, reflecting their connectivity strength across the whole brain (16, 17). Combining DC and FC analyses may overcome this limitation by first identifying hub regions and then mapping their connectivity patterns, a strategy proven effective in neuropsychiatric disorders (17, 18).

This study aimed to investigate brain functional alterations in nurses with burnout using rs-fMRI. First, we employed DC analysis to identify abnormal functional hubs associated with burnout. Subsequently, these hubs served as regions of interest (ROIs) for whole-brain FC analysis to elucidate connectivity changes specific to burnout. Additionally, correlation analyses explored the relationship between these brain alterations and burnout severity, while machine learning techniques, integrating DC and FC features, was used to develop a classification model for distinguishing burnout cases from controls.

2 Methods

2.1 Subjects

This study recruited female nurses from Yancheng Clinical Medical College of Yangzhou University as subjects, with data collection conducted from September 2024 to December 2024. Inclusion criteria: the occupational burnout group: (1) female, aged 20–40 years; (2) right-handed; (3) according to the burnout norms of Chinese nurses (19), the critical value of burnout was determined as all three-dimensional scores exceeding critical values (emotional exhaustion (EE) ≥ 27 points, depersonalization (DP) ≥ 8 points, and personal accomplishment (PA) ≤ 24 points) and the healthy control group: (1) female, aged 20–40 years; (2) right-handed; (3) all three-dimension scores below critical values (EE < 27 points, DP < 8 points, PA > 24 points). Exclusion criteria: all subjects with any of the following conditions were excluded: (1) endocrine, neurological, or psychiatric disorders or other primary diseases; (2) pregnant or lactating women; (3) history of drug dependence, smoking, or alcohol consumption; (4) adverse reactions during scanning leading to termination of the experiment or contraindications to MRI scanning; (5) data collection failure during scanning or unclear images; (6) MRI images showing organic brain lesions; (7) other serious physical illnesses. Based on the inclusion and exclusion criteria, 80 subjects were ultimately selected, with 40 in the occupational burnout group and 40 in the heathy control group. The two groups were matched in terms of age and years of education. This study strictly adhered to the ethical principles of the Declaration of Helsinki and has received approval from the Ethics Committee of the Yancheng Clinical Medical College of Yangzhou University (2024–82) and obtained informed consent from all subjects involved.

Prior to MRI scanning, general information and clinical data were collected, including age, years of education, body mass index (BMI), Beck Anxiety Inventory (BAI) (20), Beck Depression Inventory-II (BDI-II) (21), and Maslach Burnout Inventory-Human Services Survey (MBI-HSS) scale. MBI-HSS scale assesses burnout across three dimensions in service industries: EE (9 items): emotional exhaustion from work; DP (5 items): depersonalized responses to care recipients; PA (8 items, reverse-scored): feelings of competence and achievement. Items are rated on a 7-point Likert scale (0: never to 6: daily) (22). Within our sample, the Cronbach’s α coefficients for the dimensions of EE, DP, and PA in the MBI-HSS were 0.970, 0.965, and 0.960, respectively.

2.2 MRI data acquisition

Rs-fMRI and structural 3D-T1-weighted images were acquired using a 3.0 T MRI scanner with a 24-channel head coil (Discovery 750w, GE, United States) at Yancheng Clinical Medical College of Yangzhou University. Parameters included: rs-fMRI: repetition time (TR)/ echo time (TE) = 3,000/35 ms, 128 volumes, field of view (FOV) = 24 cm × 24 cm, Slice thickness = 5.0 mm, and voxel size = 3.75 × 3.75 × 4 mm; structural 3D-T1: TR = 750 ms, TE = 2.8 ms, FOV = 24 cm × 24 cm, Slice thickness = 1.0 mm, number of slices = 152, flip angle = 15°, and voxel size = 0.5 × 0.5 × 1 mm.

2.3 Rs-fMRI preprocessing

Preprocessing was performed using the DPABI 8.2 software, including: (1) removal of the first 10 time points; (2) slice timing correction; (3) realignment for head motion; (4) exclusion of participants with maximum displacement > 3 mm or rotation > 3°; (5) spatial normalization to the standard Montreal Neurological Institute (MNI) space achieved through the Dartel alignment method; (6) linear regression to reduce errors; (7) regression of nuisance covariates; and (8) band-pass filtering (0.01–0.1 Hz). All participants included in the final analysis met these head motion criteria.

2.4 Total intracranial volume (TIV) extraction

Brain structural 3D-T1-weighted images were preprocessed using SPM12 and CAT12, including bias field correction, skull stripping, alignment to MNI template, and segmentation into gray matter, white matter, and cerebrospinal fluid. TIV was extracted for all participants.

2.5 DC analysis

DC was calculated using DPABI 8.2, computing the DC value for each voxel in the brain. Pearson correlation coefficients were used to estimate functional connectivity between all pairs of gray matter voxels. A threshold of r > 0.25 was used to derive the adjacency matrix, followed by conversion of individual voxelwise DC values into a z-score map. Subsequently, the DC maps obtained were smoothed spatially using a 6-mm full width at half-maximum (FWHM) Gaussian kernel.

2.6 FC analysis

Using the AAL 90 template, brain regions with significant DC differences between the burnout and control groups were selected as ROIs for whole-brain FC analysis. The average time series for each ROI was calculated, and Fisher’s z-transformation was applied to obtain normally distributed z-score maps.

2.7 Statistical analysis

SPSS 27.0 software was used for statistical analysis. First, normality tests were conducted for age, years of education, BMI, TIV, BAI, BDI-II, and MBI-HSS. For normally distributed measurement data, independent sample t-tests were used, while non-parametric tests were applied for skewed distribution data. A p-value < 0.05 was considered to indicate significant between-group differences.

Two-sample t-tests were conducted to compare differences in FC and DC between the burnout group and the control group, with age, years of education, TIV, BAI, and BDI-II as covariates using SPM12 in MATLAB (R2020b). Multiple comparison corrections used cluster-level False Discovery Rate (FDR), with voxel-level p < 0.001 and cluster-level p < 0.05. DC and FC values from ROIs (showing significant group differences) were extracted for subsequent correlation and machine learning analyses.

2.8 Correlation analysis

SPSS 27.0 software was used to conduct Spearman correlation analysis between the extracted brain DC and FC signal values with significant between-group differences and burnout-related scales (EE, DP, PA). A p-value < 0.05 was considered to indicate significant correlation.

2.9 Classification model construction

This study used DC values and FC values from brain regions with significant differences obtained through two-sample t-tests and cluster-level FDR correction, as well as combined DC and FC values, as input features for Linear discriminant analysis (LDA) to construct classification models for burnout and healthy populations. The predictive performance of the models was evaluated using receiver operating characteristic (ROC) curves, calculating the area under curve (AUC), sensitivity, specificity, and accuracy, and comparing the performance of these three models in the validation set. Leave-one-out cross-validation (LOOCV) was used to validate the predictive ability and stability of the classification models. Subsequently, permutation tests were conducted to evaluate the significance of classification accuracy and AUC. During the permutation tests, group labels were randomly shuffled 5,000 times, and classification was performed on the newly generated datasets, calculating the classification accuracy and AUC each time. The p-value for classification accuracy and AUC was the number of times the classification accuracy or AUC in the 5,000 random cases exceeded the true value, divided by 5,000. Results were considered significant when the p-values for both accuracy and AUC were less than 0.05.

3 Results

3.1 Demographic and clinical characteristics

The study included 40 participants in occupational burnout group and healthy control group. Demographic and clinical characteristics of both groups are summarized in Table 1. No statistically significant differences were observed between the occupational burnout and healthy control groups in age (median [interquartile range, IQR]: 33 [27–37] vs. 33.5 [28–36] years; p = 0.973), years of education (median [IQR]: 16 [16–16] vs. 16 [16–16]; p = 0.724), or body mass index (BMI; median [IQR]: 21.89 [19.64–23.81] vs. 20.81 [19.55–23.63]; p = 0.544). The occupational group exhibited markedly higher anxiety levels on the BAI (median [IQR]: 27 [23.25–33] vs. 24 [22–26]; p = 0.009) and greater depressive symptoms on the BDI-II (median [IQR]: 9.5 [5.5–13.75] vs. 4 [0–9.5]; p = 0.002). Scores on MBI-HSS subscales further distinguished the groups. The occupational group reported significantly higher EE (median [IQR]: 36.0 [30.5–40.00] vs. 13.5 [9–17]; p < 0.001) and DP (median [IQR]: 13.0 [10.0–16.75] vs. 3 [0–5]; p < 0.001), alongside significantly lower PA (median [IQR]: 20.0 [17.25–22.0] vs. 38.5 [30–43.75]; p < 0.001) compared to the healthy control group.

Table 1
www.frontiersin.org

Table 1. Demographic information and clinical data.

3.2 Group differences in DC

As shown in Table 2 and Figure 1, significant differences in DC were observed between the occupational burnout group and the healthy control group. After applying cluster-level FDR correction (voxel threshold: p < 0.001 and cluster threshold: p < 0.05), the occupational burnout group exhibited significantly reduced DC values in the left precuneus (t = −4.37, p < 0.001) and the right precuneus (t = −4.55, p < 0.001) compared to the healthy control group.

Table 2
www.frontiersin.org

Table 2. Group differences in DC between nurses with occupational burnout and HCs.

Figure 1
www.frontiersin.org

Figure 1. Group differences in DC between female nurses with occupational burnout and HCs. Compared to HCs, nurses with burnout showed a significant decrease in DC in the bilateral PCUN. DC, degree centrality; HCs, healthy controls; PCUN, precuneus.

3.3 Group differences in FC

Seed-based whole-brain FC analysis was conducted using the DC-differentiated brain regions as seeds. As summarized in Table 3 and illustrated in Figure 2, the occupational burnout exhibited significantly reduced FC values between the left precuneus and the right medial orbitofrontal cortex (mOFC, t = −4.57, p < 0.001) compared to the healthy control group.

Table 3
www.frontiersin.org

Table 3. Group differences in FC between nurses with occupational burnout and HCs.

Figure 2
www.frontiersin.org

Figure 2. Group differences in FC between female nurses with occupational burnout and HCs. Compared to HCs, FC between the PCUN. L and mOFC. R is significantly reduced in nurses with burnout. FC, functional connectivity; HCs, healthy controls; PCUN, precuneus; mOFC, medial orbitofrontal cortex; L, left; R, right.

3.4 Correlation analysis

As illustrated in Figures 3A,B, EE was negatively correlated with DC values in the left precuneus (r = −0.46, p = 0.003) and right precuneus (r = −0.36, p = 0.023). Additionally, DP was negatively correlated with DC values in the left precuneus (Figure 3C) (r = −0.348, p = 0.028). As shown in Figure 3D, PA was positively correlated with FC values between the left precuneus and the right mOFC (r = 0.378, p = 0.016).

Figure 3
www.frontiersin.org

Figure 3. The relationship between DC values and FC values with clinical variables in the occupational burnout group. (A) EE was negatively correlated with DC values in the left precuneus. (B) EE was negatively correlated with DC values in the right precuneus. (C) DP was negatively correlated with DC values in the left precuneus. (D) PA was positively correlated with FC values between the left precuneus and the right orbital part of the superior frontal gyrus. DC, degree centrality; FC, functional connectivity; EE, emotional exhaustion; DP, depersonalization; PA, reduced personal accomplishment.

3.5 Classification model analysis

To evaluate the classification performance, we employed LOOCV LDA model to build models distinguishing occupational burnout individuals from healthy controls. Input features for these models were DC and FC values extracted from brain regions exhibiting significant group differences. As shown in Table 4 and Figure 4, the LDA model incorporating both DC and FC values achieved the highest performance (sensitivity = 0.800, specificity = 0.900, accuracy = 0.850, and AUC = 0.902), significantly outperforming models using only DC values (sensitivity = 0.800, specificity = 0.675, accuracy = 0.738, and AUC = 0.851) or only FC values (sensitivity = 0.575, specificity = 0.675, accuracy = 0.625, and AUC = 0.657).

Table 4
www.frontiersin.org

Table 4. The ability of an LDA model to differentiate between nurses with occupational burnout and HCs based on different input features.

Figure 4
www.frontiersin.org

Figure 4. The ability of an LDA model to differentiate between nurses with burnout and HCs based on different input features. (A) Use the DC values of the bilateral PCUN as input features. (B) Use the FC values between the PCUN. L and the mOFC. R as input features. (C) Use both the DC values of the bilateral precuneus and the FC values between the PCUN. L and the mOFC. R as input features. LDA, linear discriminant analysis; HCs, healthy controls; DC, degree centrality; FC, functional connectivity; PCUN, precuneus; mOFC, medial orbitofrontal cortex; L, left; R, right; ROC, receiver operating characteristic; TPR, true positive rate; FPR, false positive rate.

4 Discussion

This study investigated the neural correlates of occupational burnout in female nurses by examining alterations in DC and corresponding FC using rs-fMRI. Our findings revealed significant reductions of DC within the bilateral precuneus, and decreased FC between the left precuneus and the right mOFC in the occupational burnout group compared to healthy control group. Furthermore, we observed significant correlations between these neuroimaging markers and clinical burnout dimensions and demonstrated the potential of combining DC and FC measures to effectively classify individuals with occupational burnout.

The decreased DC in the precuneus in the occupational burnout group was a key finding observed in our study. The precuneus is a central hub within the default mode network (DMN), widely recognized for its involvement in a range of higher-order cognitive functions, including self-referential processing, introspection, episodic memory recall, and maintaining awareness (2325). Reduced DC in this region suggests a disruption in its functional integration and potentially altered processing of internal states in individuals experiencing burnout. This aligns with the core features of burnout, which include emotional exhaustion and a sense of detachment from work, both of which may be linked to altered self-awareness and introspective abilities (2628). Previous studies have also implicated the precuneus in stress-related conditions and emotional regulation (2931), further supporting our findings. The negative correlations between precuneus DC and both emotional exhaustion and depersonalization dimensions of the MBI-HSS further reinforced this interpretation, suggesting that lower precuneus centrality is associated with higher levels of burnout severity in these domains.

Beyond the hub properties of the precuneus, the observed reduction in FC between the left precuneus and the right mOFC offers deeper insights into the neurobiological pathways affected by burnout. The mOFC is firmly established as a pivotal area for processing reward, regulating emotions, and guiding decision-making processes (3235). The precuneus and mOFC are anatomically and functionally interconnected regions (23, 32), and their coordinated activity is likely essential for integrating self-referential information with emotional and regulatory processes. The decreased FC between these regions in burnout may indicate a disruption in this crucial communication pathway, potentially explaining why burnout patients often experience symptoms such as emotional exhaustion, reduced sense of work value, cognitive function decline, and decreased responsiveness to work-related rewards (2628, 36, 37). Interestingly, we found a positive relationship between this precuneus-mOFC connectivity and personal accomplishment, a key facet of burnout. This suggests that stronger functional integration within this neural circuit may be indicative of greater professional self-efficacy and enhanced resilience against burnout. It is conceivable that robust connectivity here facilitates a more effective integration of self-perception with the emotional and motivational systems pertinent to work, thereby providing a buffer against the development of burnout (26).

Moving beyond the neurobiological findings, the performance of our classification model points towards potential clinical utility. The superior performance of the LDA model when incorporating both DC and FC values underscores the complementary nature of these measures in capturing the neurobiological underpinnings of burnout. The high accuracy, sensitivity, and specificity achieved by the combined model, particularly with LOOCV validation, suggest that these rs-fMRI-derived metrics hold promise as objective biomarkers. Such biomarkers could be invaluable in aiding the identification and diagnosis of occupational burnout, especially considering the inherent subjectivity of current burnout assessment tools and the recognized need for more objective measures within both clinical and occupational health contexts.

While the present study offers valuable perspectives on the neural mechanisms associated with occupational burnout, it is important to acknowledge certain limitations. Firstly, this study focuses solely on female nurses aged 20–40, recruited from a single center, potentially limiting the generalizability of our findings. Future multicenter studies should incorporate more diverse samples to validate these results in a broader population and varied institutional settings. Secondly, the cross-sectional nature of our study design prevents us from establishing definitive causal relationships between the observed brain functional changes and the development of burnout. Longitudinal investigations are essential to clarify whether these neuroimaging alterations precede the onset of burnout or emerge as a consequence of prolonged occupational stress. Thirdly, while DC and FC provided a valuable initial window into neural function, future studies could benefit from exploring a wider array of neuroimaging measures. Examining cerebral blood flow, as well as structural and functional network organization, could offer a more comprehensive and nuanced understanding of the neural underpinnings of burnout. Fourthly, a significant limitation of our study is the lack of control over menstrual cycle phases or the use of hormonal contraceptives, as our research focused solely on female nurses. Given that hormonal fluctuations can impact emotional regulation and rs-fMRI signals, future studies should systematically consider these factors to gain a more comprehensive understanding of their potential effects on brain function in the context of burnout. Fifthly, while LOOCV was employed for model validation and is suitable for the current sample size, its performance on this dataset does not guarantee generalizability to entirely new, independent populations. Future research should rigorously test these neuroimaging-based classifiers on larger, external datasets, ideally using independent held-out test sets, to establish their robustness and potential for broader clinical application. Finally, while our sample size was adequate for detecting statistically significant group differences, replication of these findings in larger cohorts would further strengthen the robustness and generalizability of our conclusions.

5 Conclusion

In summary, this study identifies reduced DC in the precuneus and its decreased FC with the mOFC as key neural substrates of occupational burnout, indicating impaired integration between self-referential processing and reward/emotion regulation systems. Furthermore, the significant correlations observed with clinical burnout scales, coupled with the robust diagnostic accuracy of our integrated DC-FC model, underscore the considerable potential of neuroimaging biomarkers for the objective assessment of burnout. These findings contribute to a better understanding of the pathophysiological mechanisms underlying burnout and may pave the way for the development of targeted neuromodulation therapies.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Ethics Committee of the Yancheng Clinical Medical College of Yangzhou University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

J-PL: Data curation, Formal analysis, Visualization, Writing – original draft, Writing – review & editing. S-YG: Data curation, Investigation, Project administration, Writing – original draft. C-MS: Data curation, Investigation, Project administration, Writing – original draft. H-CY: Data curation, Formal analysis, Writing – review & editing. YS: Data curation, Formal analysis, Writing – review & editing. Y-FG: Data curation, Formal analysis, Writing – review & editing. S-FW: Conceptualization, Methodology, Supervision, Validation, Writing – review & editing. Y-ZC: Conceptualization, Methodology, Supervision, Validation, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by Yancheng Science and Technology Bureau (YCBK2024084), the College-local Collaborative Innovation Research Project of Jiangsu Medical College (20229150), and the Nantong University Clinical Medicine Special Project (2024JY022, 2024JY023, 2024JY024, and 2024LZ003).

Acknowledgments

The authors gratefully acknowledge the female nurses who participated in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Verhavert, Y, Deliens, T, Van Cauwenberg, J, Van Hoof, E, Matthys, C, de Vries, J, et al. Associations of lifestyle with burnout risk and recovery need in Flemish secondary schoolteachers: a cross-sectional study. Sci Rep. (2024) 14:3268. doi: 10.1038/s41598-024-53044-w

PubMed Abstract | Crossref Full Text | Google Scholar

2. Chen, Z, Leng, J, Pang, Y, He, Y, Heng, F, and Tang, L. Demographic, occupational, and societal features associated with burnout among medical oncology staff members: cross-sectional results of a Cancer Center in Beijing, China. Psychooncology. (2019) 28:2365–73. doi: 10.1002/pon.5230

PubMed Abstract | Crossref Full Text | Google Scholar

3. Ge, MW, Hu, FH, Jia, YJ, Tang, W, Zhang, WQ, and Chen, HL. Global prevalence of nursing burnout syndrome and temporal trends for the last 10 years: a meta-analysis of 94 studies covering over 30 countries. J Clin Nurs. (2023) 32:5836–54. doi: 10.1111/jocn.16708

PubMed Abstract | Crossref Full Text | Google Scholar

4. Wang, QQ, Lv, WJ, Qian, RL, and Zhang, YH. Job burnout and quality of working life among Chinese nurses: a cross-sectional study. J Nurs Manag. (2019) 27:1835–44. doi: 10.1111/jonm.12884

PubMed Abstract | Crossref Full Text | Google Scholar

5. Rudman, A, Arborelius, L, Dahlgren, A, Finnes, A, and Gustavsson, P. Consequences of early career nurse burnout: a prospective long-term follow-up on cognitive functions, depressive symptoms, and insomnia. EClinicalMedicine. (2020) 27:100565. doi: 10.1016/j.eclinm.2020.100565

PubMed Abstract | Crossref Full Text | Google Scholar

6. Lee, M, and Cha, C. Interventions to reduce burnout among clinical nurses: systematic review and meta-analysis. Sci Rep. (2023) 13:10971. doi: 10.1038/s41598-023-38169-8

PubMed Abstract | Crossref Full Text | Google Scholar

7. Ma, Y, Faraz, NA, Ahmed, F, Iqbal, MK, Saeed, U, Mughal, MF, et al. Curbing nurses' burnout during COVID-19: the roles of servant leadership and psychological safety. J Nurs Manag. (2021) 29:2383–91. doi: 10.1111/jonm.13414

PubMed Abstract | Crossref Full Text | Google Scholar

8. Lu, Q, Jia, Z, and Gu, H. Association between brain resting-state functional activities and migraine: a bidirectional mendelian randomization study. Sci Rep. (2024) 14:23901. doi: 10.1038/s41598-024-74745-2

PubMed Abstract | Crossref Full Text | Google Scholar

9. Tahmi, M, Kane, VA, Pavol, MA, and Naqvi, IA. Neuroimaging biomarkers of cognitive recovery after ischemic stroke. Front Neurol. (2022) 13:923942. doi: 10.3389/fneur.2022.923942

PubMed Abstract | Crossref Full Text | Google Scholar

10. Puvogel, S, Blanchard, K, Casas, BS, Miller, RL, Garrido-Jara, D, Arizabalos, S, et al. Altered resting-state functional connectivity in hiPSCs-derived neuronal networks from schizophrenia patients. Front Cell Dev Biol. (2022) 10:935360. doi: 10.3389/fcell.2022.935360

PubMed Abstract | Crossref Full Text | Google Scholar

11. Zhao, F, Chen, Z, Rekik, I, Lee, SW, and Shen, D. Diagnosis of autism Spectrum disorder using central-moment features from low- and high-order dynamic resting-state functional connectivity networks. Front Neurosci. (2020) 14:258. doi: 10.3389/fnins.2020.00258

PubMed Abstract | Crossref Full Text | Google Scholar

12. Jovanovic, H, Perski, A, Berglund, H, and Savic, I. Chronic stress is linked to 5-HT(1A) receptor changes and functional disintegration of the limbic networks. NeuroImage. (2011) 55:1178–88. doi: 10.1016/j.neuroimage.2010.12.060

PubMed Abstract | Crossref Full Text | Google Scholar

13. Chen, M, Su, Q, Zhao, Z, Li, T, Yao, Z, Zheng, W, et al. Rich Club reorganization in nurses before and after the onset of occupational burnout: a longitudinal MRI study. J Magn Reson Imaging. (2024) 60:1918–31. doi: 10.1002/jmri.29288

PubMed Abstract | Crossref Full Text | Google Scholar

14. Bai, L, Ji, GJ, Song, Y, Sun, J, Wei, J, Xue, F, et al. Dynamic brain connectome and high risk of mental problem in clinical nurses. Hum Brain Mapp. (2021) 42:5300–8. doi: 10.1002/hbm.25617

PubMed Abstract | Crossref Full Text | Google Scholar

15. Dong, Y, Wu, X, Zhou, Y, and Qiu, K. Differences in functional activity and connectivity in the right Frontoparietal network between nurses working Long-term shifts and fixed day shifts. J Integr Neurosci. (2024) 23:9. doi: 10.31083/j.jin2301009

PubMed Abstract | Crossref Full Text | Google Scholar

16. Yang, L, Zhang, L, Liu, Y, Liu, J, Li, K, and Cai, J. The different impacts of pain-related negative emotion and trait negative emotion on brain function in patients with inflammatory bowel disease. Sci Rep. (2024) 14:23897. doi: 10.1038/s41598-024-75237-z

PubMed Abstract | Crossref Full Text | Google Scholar

17. Cañete-Massé, C, Carbó-Carreté, M, Peró-Cebollero, M, Cui, SX, Yan, CG, and Guàrdia-Olmos, J. Abnormal degree centrality and functional connectivity in down syndrome: a resting-state fMRI study. Int J Clin Health Psychol. (2023) 23:100341. doi: 10.1016/j.ijchp.2022.100341

PubMed Abstract | Crossref Full Text | Google Scholar

18. Shan, A, Zhang, H, Gao, M, Wang, L, Cao, X, Gan, C, et al. Aberrant voxel-based degree centrality and functional connectivity in Parkinson's disease patients with fatigue. CNS Neurosci Ther. (2023) 29:2680–9. doi: 10.1111/cns.14212

PubMed Abstract | Crossref Full Text | Google Scholar

19. Ye, Z, Luo, H, and Jiang, A. Diagnostic standard and norms of Maslach Burnout Inventory for nurses in Hangzhou. Chin J Nurs. (2008) 43:207. doi: 10.3761/j.iSSN.0254-1769.2008.03.005

Crossref Full Text | Google Scholar

20. Beck, AT, Epstein, N, Brown, G, and Steer, RA. An inventory for measuring clinical anxiety: psychometric properties. J Consult Clin Psychol. (1988) 56:893–7. doi: 10.1037/0022-006X.56.6.893

PubMed Abstract | Crossref Full Text | Google Scholar

21. Beck, AT, Steer, RA, and Brown, OK. Beck depression inventory manual, 2nd Edn. San Antonio, TX: Psychological Corporation. (1996).

Google Scholar

22. Maslach, C, Jackson, SE, and Leiter, MP. Maslach burnout inventory manual. Edina, MN, USA: Consulting Psychologists (1996).

Google Scholar

23. Zhang, S, and Li, CS. Functional connectivity mapping of the human precuneus by resting state fMRI. NeuroImage. (2012) 59:3548–62. doi: 10.1016/j.neuroimage.2011.11.023

PubMed Abstract | Crossref Full Text | Google Scholar

24. Cavanna, AE, and Trimble, MR. The precuneus: a review of its functional anatomy and behavioural correlates. Brain. (2006) 129:564–83. doi: 10.1093/brain/awl004

PubMed Abstract | Crossref Full Text | Google Scholar

25. Dadario, NB, and Sughrue, ME. The functional role of the precuneus. Brain. (2023) 146:3598–607. doi: 10.1093/brain/awad181

PubMed Abstract | Crossref Full Text | Google Scholar

26. Bayes, A, Tavella, G, and Parker, G. The biology of burnout: causes and consequences. World J Biol Psychiatr. (2021) 22:686–98. doi: 10.1080/15622975.2021.1907713

PubMed Abstract | Crossref Full Text | Google Scholar

27. Dam, A. A clinical perspective on burnout: diagnosis, classification, and treatment of clinical burnout. Eur J Work Organ Psy. (2021) 30:732–41. doi: 10.1080/1359432x.2021.1948400

Crossref Full Text | Google Scholar

28. Tavella, G, Hadzi-Pavlovic, D, and Parker, G. Burnout: redefining its key symptoms. Psychiatry Res. (2021) 302:114023. doi: 10.1016/j.psychres.2021.114023

PubMed Abstract | Crossref Full Text | Google Scholar

29. Messina, I, Grecucci, A, and Viviani, R. Neurobiological models of emotion regulation: a meta-analysis of neuroimaging studies of acceptance as an emotion regulation strategy. Soc Cogn Affect Neurosci. (2021) 16:257–67. doi: 10.1093/scan/nsab007

PubMed Abstract | Crossref Full Text | Google Scholar

30. Lee, D, Kim, W, Lee, JE, Lee, J, Kim, YT, Lee, SK, et al. Changes in intrinsic functional brain connectivity related to occupational stress of firefighters. Psychiatry Res. (2022) 314:114688. doi: 10.1016/j.psychres.2022.114688

PubMed Abstract | Crossref Full Text | Google Scholar

31. Li, X, Qureshi, MNI, Laplante, D, Elgbeili, G, Jones, SL, Long, X, et al. Atypical brain structure and function in young adults exposed to disaster-related prenatal maternal stress: project ice storm. J Neurosci Res. (2023) 101:1849–63. doi: 10.1002/jnr.25246

PubMed Abstract | Crossref Full Text | Google Scholar

32. Zald, DH, McHugo, M, Ray, KL, Glahn, DC, Eickhoff, SB, and Laird, AR. Meta-analytic connectivity modeling reveals differential functional connectivity of the medial and lateral orbitofrontal cortex. Cereb Cortex. (2014) 24:232–48. doi: 10.1093/cercor/bhs308

PubMed Abstract | Crossref Full Text | Google Scholar

33. Rolls, ET. Emotion, motivation, decision-making, the orbitofrontal cortex, anterior cingulate cortex, and the amygdala. Brain Struct Funct. (2023) 228:1201–57. doi: 10.1007/s00429-023-02644-9

PubMed Abstract | Crossref Full Text | Google Scholar

34. Henssen, A, Zilles, K, Palomero-Gallagher, N, Schleicher, A, Mohlberg, H, Gerboga, F, et al. Cytoarchitecture and probability maps of the human medial orbitofrontal cortex. Cortex. (2016) 75:87–112. doi: 10.1016/j.cortex.2015.11.006

PubMed Abstract | Crossref Full Text | Google Scholar

35. Elliott, R, Dolan, RJ, and Frith, CD. Dissociable functions in the medial and lateral orbitofrontal cortex: evidence from human neuroimaging studies. Cereb Cortex. (2000) 10:308–17. doi: 10.1093/cercor/10.3.308

PubMed Abstract | Crossref Full Text | Google Scholar

36. Edú-Valsania, S, Laguía, A, and Moriano, J. Burnout: a review of theory and measurement. Int J Environ Res Public Health. (2022) 19:1780. doi: 10.3390/ijerph19031780

PubMed Abstract | Crossref Full Text | Google Scholar

37. Gavelin, HM, Domellöf, M, Åström, E, Nelson, A, Launder, NH, Stigsdotter-Neely, A, et al. Cognitive function in clinical burnout: a systematic review and meta-analysis. Work Stress. (2021) 36:86–104. doi: 10.1080/02678373.2021.2002972

Crossref Full Text | Google Scholar

Glossary

AAL - anatomical automatic labeling

AUC - area under curve

BAI - Beck Anxiety Inventory

BDI-II - Beck Depression Inventory-II

BMI - body mass index

BOLD - blood oxygen level-dependent

DC - degree centrality

DMN - default mode network

DP - depersonalization

EE - emotional exhaustion

FC - functional connectivity

FDR - False Discovery Rate

FOV - field of view

FPR - false positive rate

FWHM - full width at half-maximum

HCs - healthy controls

IQR - interquartile range

L - left

LDA - linear discriminant analysis

LOOCV - leave-one-out cross-validation

MBI-HSS - Maslach Burnout Inventory-Human Services Survey

MNI - Montreal Neurological Institute

mOFC - medial orbitofrontal cortex

PA - personal accomplishment

PCUN - precuneus

R - right

ROC - Receiver Operating Characteristic

ROIs - regions of interest

rs-fMRI - resting-state functional magnetic resonance imaging

TE - echo time

TIV - total intracranial volume

TPR - true positive rate

TR - repetition time

3D-T1 - structural 3D-T1-weighted

Keywords: burnout, resting-state functional MRI, degree centrality, functional connectivity, precuneus, medial orbitofrontal cortex

Citation: Liu J-P, Gu S-Y, Song C-M, Yang H-C, Shi Y, Gu Y-F, Wang S-F and Chen Y-Z (2025) Abnormal intrinsic functional hubs and connectivity in nurses with occupational burnout: a resting-state fMRI study. Front. Public Health. 13:1595550. doi: 10.3389/fpubh.2025.1595550

Received: 18 March 2025; Accepted: 29 May 2025;
Published: 16 June 2025.

Edited by:

Masafumi Yoshimura, Faculty of Rehabilitation Kansai Medical University, Japan

Reviewed by:

Yanzhen Zhang, University of California, Irvine, United States
Francisco Palencia-Sánchez, Pontifical Javeriana University, Colombia

Copyright © 2025 Liu, Gu, Song, Yang, Shi, Gu, Wang and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ying-Zhu Chen, eXpjaGVuZHJAMTYzLmNvbQ==; Shu-Fang Wang, MjcwMDEyODQzQHFxLmNvbQ==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.