key: cord-0737535-1t0nk08o
authors: Tsou, Meng-Ting
title: Association of 5-item Brief Symptom Rating Scale scores and health status ratings with burnout among healthcare workers
date: 2022-05-03
journal: Sci Rep
DOI: 10.1038/s41598-022-11326-1
sha: dfc10186ea2e4f87e03bb47b2f9789a23abf0596
doc_id: 737535
cord_uid: 1t0nk08o

This cross-sectional study evaluated and quantified the possible association of psychological symptoms and health status ratings on the burnout of healthcare workers from a tertiary medical center. Demographic data were obtained through a questionnaire survey. We evaluated their psychological symptoms using a 5-item Brief Symptom Rating Scale (BSRS-5) and burnout was measured using the Chinese version of the Maslach Burnout Inventory–Health Services Survey. The study began in Nov. 2018 and ended in Nov. 2020. A total of 2813 participants (men = 296, 10.5%; women = 2517, 89.5%) completed the questionnaires between December 2018 and March 2019. The response rate and complete rate were 71.7% and 93.2%. The multivariate analysis showed that, as the BSRS-5 level added, the odds ratio (OR) of burnout increased (BSRS-5 scores 6–9, 10–14, and 15–20; OR = 1.83, 3.23, and 9.15, separately; p value < 0.05). Overall, men (≥ 30 years of age) and women staffs with BSRS-5 scores ≥ 6; women with longer working hours (more than 46 h/week), men and women (≥ 30 years of age) working night shifts, and poor health status ratings were highly associated with burnout. The findings highlight the importance of screening for the BSRS-5 scores and health status ratings level for healthcare professionals at high risk of burnout, especially men ≥ 30 years of age and women with stressful working conditions.

The potential participants were 4457 persons. 1262 staffs refused to joint it. Response rate of survey was 65.3% (doctor), 88.2% (nurse), and 54.1% (non-doctor/nurse). A total of 187 doctors, 1860 nurses, and 970 non-doctor/ nurse participants were invited after performing the exclusion procedure (see Fig. 1 ), and 110, 1758, and 945 cases from the three groups completed the survey. The actual response rates were 58.8% (110/187) in the doctor group, 94.5% (1758/1860) in the nurse group, and 97.4% (945/970) in the non-doctor/nurse group. In total of 2813 participants were included in this study, implying a sufficient statistical power. Reliability of MBI-HSS and BSRS-5 questionnaire for all participants, stratified by sex and age was performed in Table 1 . The results showed the internal consistency coefficients (Cronbach's α) of the MBI-HSS, and BSRS-5 were 0.79 (acceptable) and 0.89 (good), respectively (22) (23) (24) for all participants. Stratifying by sex and age, the Cronbach's α of MBI-HSS was more than 0.8 (good) for 60-79 years among both men and women, and between 0.7 and 0.8 (acceptable) for less than 60 years. The Cronbach's α of BSRS-5 score was more than 0.9 (excellent) for men among different age groups, and between 0.8 and 0.9 (good) for women among different age groups.

This study included a total of 2582 non-burnout and 231 burnout participants. The prevalence of men in the burnout group was significantly higher than in the non-burnout group (15.22% vs. 10 .19%, p = 0.02; Table 2 ). The college education level and exercise habits of the non-burnout group were significantly higher than those of the burnout group. Table 2 shows no statistical difference in terms of department and work style between the two groups. The seniority of the burnout group was less than ten years. A total of 45.8% of the non-burnout group had more than ten years of experience. The working hours of the burnout group exceeded those of the non-burnout group-approximately 48.5% of the burnout group reported working for more than 46 h per week. The BSRS-5 score ≥ 6 subgroup proportion was higher in the burnout group than in the non-burnout group (76.62% vs. 32.18%, p < 0.001). Sleep duration was shorter for the burnout group. Poor self-assessed physical health status and inadequate self-assessed sleep time were found in the burnout group (p < 0.001).

The participants' data was divided into four categories of BSRS-5 scores, as shown in Table 3 ; younger age and women were found predominantly in the groups with higher BSRS-5 scores (p < 0.001 and p = 0.003). We also identified a significant trend of association between less exercise, drinking habit, and above college level education and higher BSRS-5 scores ( Table 3) .

The seniority of the higher BSRS-5 scores was primarily 4-10 years. We also found a significant decreasing trend of association between seniority of more than 10 years and higher BSRS-5 scores. Similarly, working more than 46 h per week and night shifts were related to higher BSRS-5 scores. There was a significant trend of association between shorter sleep duration, poor self-assessed physical health status, poor self-assessed sleep quality, and higher burnout prevalence and higher BSRS-5 scores (p < 0.001). Medium to high effect size test in ANOVA/chi square test was in self-assessment physical health status and burnout. Medium to high effect size test in Cochran-Armitage Trend test was exercise, drinking, education, and seniority; high effect size test in Cochran-Armitage Trend test was department, working hours/week, work style, s self-assessment physical health status, self-assessment of sleep time, and burnout. Figure 2 shows the relationship between the three domains of the MBI-HSS scores and the four categories of the BSRS-5 scores (p < 0.05). Table 4A shows that EE and DP have a good correlation with anxiety, depression, and total BSRS-5 scores (r = 0.5-0.75, p < 0.001); however, PA was irrelevant by Spearman's correlation. Table 4B showed that EE and DP had good correlation with all items except poor sleep quality among men less than 30 years, DP had good correlation with anxiety, inferiority, and total BSRS-5 score among men more than 30 years. EE and DP had a good correlation with anxiety, depression, hostility, and total BSRS-5 scores among women with different age group. Supplement Table 2 showed the similar results by Pearson correlation. The results showed a higher prevalence of burnout among 20-39 years old women with a BSRS-5 score > 6 or the health status ratings measure ≥ 1 ( www.nature.com/scientificreports/ The logistic regression analysis was performed to estimate the crude and multivariable-adjusted odds (OR) ratios and identify factors associated with burnout ( Table 6 ). The BSRS-5 scores, age, sex, and factors (working characters and health status ratings that were significantly associated with burnout in crude analysis were included in the multiple logistic regression model. The BSRS-5 scores (compared to 0-5, OR of 6-9, 10-14, and 15-20 = 1.83, 3.23, 9.15, respectively; p < 0.01) and being a man (OR = 1.88, 95% confidence interval [CI] = 1.15-3.09) were still associated with burnout, whereas the work characteristics that remained significant included a seniority of 2-4 years and 4-10 years (OR = 1.79, 1.69, p < 0.05), poor self-assessed health status (OR = 3.31, 95% CI 1.84-5.97), and poor self-assessed sleep quality (OR = 1.78, 95% CI 1.08-2.23), and were significantly associated with burnout. 

To the best of our knowledge, this is the first study to demonstrate that the severity of psychological symptoms, measured by the BSRS-5, can significantly associated with burnout among different sex and age, and that this effect persists after controlling for common potential confounders, including age, sex, smoking, alcohol consumption, exercise, and chronic diseases (HTN, DM, and hyperlipidemia). We also performed the stratified analysis to identify the association of high workload and poor health status ratings with burnout among middle-aged men and women.

The result showed that the prevalence of burnout among the entire sample population, men, and women was 8.2%, 11.8% and 7.8%, respectively (p = 0.04). Compared with other studies 1,2,4,5 , the possible reasons for the lower incidence of burnout found in this study included different departments including in this study, not only physicians and nurses. Second, the participants in this study were generally younger than those in other studies (median 38.4 years) [19] [20] [21] . Finally, 86.8% of the participants in this study possessed a higher level of education (higher than university) as compared to participants in other studies [19] [20] [21] . A previous study described a high risk of burnout mostly among women with low education levels 22 .

Burnout, SRH and psychologic symptoms are not consistent with age and gender in HCWs in different departments 12, 13, 23, 24 . Previous study had clarified that age and sex are a potential modulator of health risks caused by psychosocial work conditions 25 . Surveys of physicians from different specialties found that burnout was self-reported had found that women were more likely to exhibit burnout symptoms than were men [26] [27] [28] . Age was an independent variable associated with burnout, with younger physicians at a higher risk of burnout than older physicians 29 . In previous studies mentioned that it was also worth noting that age may influence burnout more than gender did 27 .

Most of previous studies are based on different indicators to discuss age or sex stratification 13, 21, 24, 29 . Few studies have combined psychological symptoms, health status ratings, age, and sex into groups to explore the correlation with burnout. This study shows that whether men and women are over 30 with a high BSRS-5 score or poor health status ratings is higher than that of people under 30 years old. This echoes with previous studies of individual indicators, these previous results show that poor SRH increases with age, and the level of burnout among employees aged 30-40 is higher than that of other age groups 11, 13 .

In terms of sex, previous studies proposed that burnout symptoms vary considerably according to the different life stages of working men and women 28, 30 . Their research found that two groups including aged 20-35 years and aged 55 years and above were particularly susceptible to burnout, and that different psychosocial conditions needed to be considered. Therefore, in our study, we included different psychosocial factors (such as BSRS-5 and health status ratings) that caused different levels of burnout among men and women of different age groups.

Studies have mentioned that physical and mental diseases affect psychological symptoms and burnout 15, 31 . We excluded participants with uncontrolled chronic diseases and mental and sleep disorders. However, we found that some under-treatment chronic diseases (HTN Recent research has stated that poor sleep quality is common, and associated with increased burnout, among clinicians delivering care to COVID-19 patients 32 . In this model, self-reported poor sleep quality was significantly associated with burnout (OR = 4.13, 95% CI 2.33-7.32, p < 0.05), whereas short sleep duration (< 6 h) was not (OR = 0.73, 95% CI 0.41-1.30, p = 0.28). Interestingly, we found that sleep quality was more strongly related 13 . Our study shows that middle-aged/elderly HCWs with longer working hours and night shifts exhibit greater physical and psychological reactions than younger staff, especially among women. Among the work characteristics, longer working hours (more than 46 h/week) for women and night shift jobs for both men and women aged 30 years or above were found to be the most common predictors of high burnout. These findings support the results of a previous study conducted in Taiwan and other countries 13,33-35 . Table 3 . Demographic and clinical characteristics of study participants across different level of BSRS-5. a p < 0.05 versus BSRS-5 score ≦ 5 group. b p < 0.05 versus BSRS-5 score 6-9 group. c p < 0.05 versus BSRS-5 score 10-14 group in the Bonferroni post hoc comparisons. BSRS-5: 5-item Brief Symptom Rating Scale, n = numbers, SD = Standard deviation, HTN = hypertension, DM = diabetes mellitus, BMI = body mass index, WC = waist circumference. 37, 38 . This may be due to the influence of the work characteristics.

Although burnout, anxiety, and depression have many overlapping issues, Fischer et al. focused on 715 critical care clinicians in Brazil and, using appropriate statistical methods, they found that burnout was statistically distinct from anxiety and depression 15 . These findings suggest that health professionals at high risk of stress need to be screened for both burnout and clinical psychological symptoms to provide timely and efficient treatment. Even the BSRS-5 is not a generic questionnaire and is not commonly used in burnout studies. The reasons why we chose it for the measurement tool first: it contains only five items and can be administered to subjects to complete independently and quickly. This questionnaire is already recommended to the general public to detect psychological symptoms and suicide prevention in Taiwan [39] [40] [41] . Second: good reliability of MBI-HSS and BSRS-5 was checked among all participants in our study (Table 1) . This study provides additional evidence that may be suitable for the evaluation of psychological symptoms in HCWs with burnout syndrome.

The strengths of our study are as follows: first, we use simple BSRS-5 score and health status ratings problems to reveal the association and evaluate the burnout problem among HCWs. We also performed age and sex stratification to facilitate burnout prevention in high-risk HCWs of a particular age and sex. Second, we excluded participants with untreated chronic diseases and mental and sleep disorders to prevent any interference with the results for burnout.

Alternatively, this study had some limitations. First, nearly 90% of the participants were women (men = 296, 10.5%; women = 2517, 89.5%). Although the chi-square test showed significant differences between the proportions of men and women for the burnout categories and levels of the BSRS-5 scores, there were some restrictions in the case of men when stratifying analysis by age and sex because of the small sample size (men: < 30 years of age, n = 69; ≥ 30 years of age, n = 227). Due to the absence of measurement of social desirability, in the future, data should be collected from multiple hospitals and regions, and studies should include more male participants. Second, this study was cross-sectional, and it was insufficient to confirm the causes and effects of psychological symptoms and burnout. Third, this study did not discuss family-related factors, including the pressure and time spending of taking care of children, organizing housework, etc., which have made the gender-related differences between work and non-work stressors. Fourth, we evaluated health status ratings by asking two questions including two items (physical health condition and sleep quality) because of the limitation related to the total number of items in the questionnaire; we did not use the more common "Perceived Wellness Survey, " which comprises 36 items 42 . Finally, this study was a voluntary response sample, and there was a problem of self-selection bias.

The BSRS-5 score is associated with the scores of at least two domains of the Chinese version of MBI-HSS for HCWs, whereas mild physical illness (such as HTN, DM, and hyperlipidemia) may not influence the results in the case of some participants. Researchers and health practitioners should pay attention to the influence of age Socio-demographic data and past medical history. Information was collected using structured questionnaires, designed according to the instructions of the Institute of Labor, Occupational Safety, Ministry of Labor (Survey of Perceptions of Safety and Health in the Work Environment in 2016 Taiwan; ILOSH105-M309) 45 . The following data were collected: socioeconomic status, work characteristics, job burnout syndrome status, and mental health status. Analysis of the internal consistency reliability of each part of the scale showed Cronbach's α to be between 0.82 and 0.93, indicating good reliability. The age group was divided to young adult (20-39 years old), middle-aged adult (40-59 years old), and old adult (60-79 years old) 46 . According to distribution, we divided into < 30 Y/O and ≧ 30 Y/O". The presence of hypertension (HTN), diabetes mellitus (DM), and hyperlipidemia was defined as a previous diagnosis of diseases or current medication use for these diseases.

including age, body height, body weight, body mass index, and waist circumference (WC), were recorded. Height was measured to the nearest 0.01 cm using a standard stadiometer. Weight was measured in light clothes to the nearest 0.01 kg using a set of standard calibrated electronic scales. The WC was measured using a constanttension tape at the midpoint between the lowest rib and the upper point of the iliac crest, and at the end of normal expiration. Anthropometric measurements such as height, weight, and WC were examined and recorded by trained nurses who were blinded to the patients' information. Table 5 . The association between level of BSRS-5 score and health status ratings with the prevalence of burnout by using chi square test stratified by age and gender. BSRS-5 = 5-item Brief Symptom Rating Scale. Definition: Health status ratings included 1. self-assessment physical health status (good, common, and poor); and 2. self-assessment sleep quality (good and poor). 0 means the answers included self-assessment physical health status (good or moderate) and self-assessment sleep quality (good). 1 means the answers included selfassessment physical health status (good or moderate) and self-assessment sleep quality (poor); self-assessment physical health status (poor) or self-assessment sleep quality (good). 2 means the answers included selfassessment physical health status (poor) and self-assessment sleep quality (poor). www.nature.com/scientificreports/ Work characteristics. Basic work-related issues included department and seniority. Information about working hours and shift jobs was derived from the work status provided by participants during the month before the survey. A night shift was described as that which encompassed midnight; work shifts that began in the morning, afternoon, or evening but did not extend to midnight were classified as day shifts. The participants were categorized into the following: day shift (day fixed/rotation) and night shift (night fixed/rotation) 47 . 48 . Three domains of the MBI-HSS were evaluated: emotional exhaustion (EE), depersonalization (DP), and personal accomplishment (PA) (Supplement Table 1 ). The MBI-HSS includes 22 questions in total, and the frequency of occurrence is scored using a 7-point Likert scale from 0 (never felt) to 6 (felt every day). According to the classification standard of previous study 48 , if EE score ≧ 27, DP score ≧ 13, and PA ≦ 31 was categorized as high job burnout; if EE score is between 17 and 26, DP score is between 7 and 12, and PA score is between 32 and 38, it is considered to be moderate job burnout; if EE score ≦ 16, DP score ≦ 6, and PA score ≧ 39 are considered low job burnout. The psychological work demands, EE and DP showed a direct relation to burnout, whereas PA showed a negative association 4 . The internal consistency coefficients (Cronbach's α) of the EE, DP, and PA dimensions were 0.77, 0.83, and 0.82, respectively 48 .

Assessment of psychological symptoms: the 5-item Brief Symptom Rating Scale. The BSRS-5 contains five items of psychological symptoms and is commonly used for screening psychological disorders. It is available in Taiwan with excellent validity and reliability 40,41 (Supplement Table 1) . It is a self-administered questionnaire that measures anxiety (feeling tense or high-strung), depression (feeling depressed or in a low mood), hostility (feeling easily annoyed or irritated), interpersonal sensitivity (feeling inferior to others), and additional symptoms (having trouble falling asleep) during the past week. Responses are rated on a scale of 0-4 (0 = not at all; 1 = a little bit; 2 = moderately; 3 = quite a bit; and 4 = extremely). The total score ranges from 0 to 20 39 . When a score of ≥ 6 was used as the cutoff point for psychiatric cases, the BSRS-5 was divided into four groups: "no symptoms" (0-5), "mild" (6-9), "moderate" (10) (11) (12) (13) (14) , and "severe" (> 15) 40 . The internal consistency was analyzed using Cronbach's α coefficient: the BSRS-5 coefficients ranged from 0.77 to 0.90 (20) , and the test-retest reliability coefficient was 0.82. Using a score of ≥ 6 as the cutoff point for psychiatric cases, the BSRS-5 accurate classification rate was 76.3% (sensitivity 78.9%; specificity 74.3%; positive predictive value 69.9%; and negative predictive value, 82.3%) 40 .

Measure of health status: health status ratings. Health status ratings was assessed using two-item questions for physical health and sleep quality issues. The physical issue question was: "In general, how is your health?" it had three possible answers ranging: "good, " "moderate, " or "poor. " In case of the sleep quality issue, which was mentioned in the literature 6, 13 : "How about your sleep quality?" which had two possible answers: "good" and "poor" 13 . Table 7 . The association between BSRS-5, working character, and perception of wellness and burnout by using multivariate logistic regression analysis stratified by age and sex. Adjusted factors: body mass index, waist circumference, personal habits, chronic diseases, and education level. www.nature.com/scientificreports/ The health status ratings levels were divided into answers of 0, 1, and 2; their definitions are as follows: 0 means that the answers included self-assessed physical health status as good or moderate and self-assessed sleep quality as good; 1 means the answers included self-assessed physical health status as common or moderate or self-assessed sleep quality as poor; self-assessed physical health status as poor or self-assessed sleep quality as good; and 2 means the answers included self-assessed health status as poor and self-assessed sleep quality as poor.

Statistical analysis. In this study, the findings of screening for the BSRS-5 scores and health status ratings level for healthcare professionals at high risk of burnout was primary outcome; and stratification analysis for age and sex was second outcome. Data are summarized as means and standard deviations for continuous variables or frequencies and percentages for categorical variables. The distribution of burnout and non-burnout syndrome according to basic characteristics was analyzed using chi-square tests for categorical variables and Student's t-tests for continuous variables. We conducted an analysis of variance and chi-square tests for evaluating differences among variables; the Cochran-Armitage test for trend and the Bonferroni post hoc comparisons were performed for the participants, who were divided into four groups according to BSRS-5 score intervals of 0-5, 6-9, 10-14, and 15-20. Effect size, including Cohen's W for chi-square tests and eta squared for ANOVA, was performed as a quantitative measure of the magnitude of the experimental effect 43 . The larger the effect size the stronger the relationship between two variables. Based on the rules of thumb 49 , we categorized the magnitudes of effect sizes into small, medium, and large.

Spearman's correlation and Pearson correlation were performed between the five items of the BSRS-5 and the three domains of burnout 39, 48, 50 . Pearson correlation measures the "linear" correlation. The three domains of MBI-HSS and BSRS-5 in this study are "non-linear" variables, which is the main reason for adopting Spearman's correlation in this study. Stratifying by sex and age was performed. The univariate logistic regression analyses were performed to investigate the possible associations between burnout and the BSRS-5, work-related factors, psycho-physical symptoms, and chronic disease status. Thereafter, the multivariate logistic regression analyses were performed in order to adjust for age, sex, and other variables that showed significant associations in the univariate analyses. Finally, we evaluated the association between the BSRS-5, work characteristics, perception of wellness, and burnout in the stratum of age and sex using the logistic regression analysis. All analyses were performed using SPSS 22 (IBM Corp., Armonk, NY, USA) for Windows. All tests were 2-tailed and p values of < 0.05 were considered statistically significant.

SABA Study Investigators and the Asian Critical Care Clinical Trials Group. Professional burnout among physicians and nurses in Asian intensive care units: A multinational survey

Burnout in anesthesiology, critical care, and perioperative care clinicians in Asia: A developing picture

Minor mental disorders in Taiwanese healthcare workers and the associations with psychosocial work conditions

Nurse burnout in Taiwan

Surveys of burnout among physicians in Taiwan

BSRS-5 (5-item Brief Symptom Rating Scale) scores affect every aspect of quality of life measured by WHOQOL-BREF in healthy workers

Is precarious employment damaging to self-rated health? Results of propensity score matching methods, using longitudinal data in South Korea

Occupational burnout and medically certified sickness absence: A population-based study of Finnish employees

Burnout levels and self-rated health prospectively predict future long-term sickness absence: A study among female health professionals

Burnout and risk of regional musculoskeletal pain-a prospective study of apparently healthy employed adults

Psychological symptoms among hospital nurses in Taiwan: A cross sectional study

Burnout and metabolic syndrome among different departments of medical center nurses in Taiwan-cross-sectional study and biomarker research

The influence of age on the distribution of self-rated health, burnout and their associations with psychosocial work conditions

Workplace relationships impact self-rated health: A survey of Swedish municipal health care employees

Association of burnout with depression and anxiety in critical care clinicians in Brazil

The relationship between Burnout, depression, and anxiety: A systematic review and meta-analysis

Women and gender in research on work and family stress

Occupational and nonwork stressors among female physicians in Taiwan: A single case study

Job stress and burnout in hospital employees: Comparisons of different medical professions in a regional hospital in Taiwan

Prevalence of burnout syndrome in emergency nurses: A meta-analysis

Burnout and metabolic syndrome in female nurses: An observational study

Burnout among women: Associations with demographic/socio-economic, work, lifestyle and health factors

Age and sex differences in burnout, career satisfaction, and well-being in US neurologists

effort-reward imbalance at work, and burnout: Findings from a developing Caribbean country

Age as factor in the relation between work and mental health: Results of the longitudinal TAS survey

Medscape national physician burnout and depression

Burnout and psychological distress among pediatric critical care physicians in the United States

Gender-based differences in burnout: Issues faced by women physicians. NAM Perspectives. Discussion Paper

Physician burnout: Contributors, consequences, and solutions

Do age and gender contribute to workers' burnout symptoms?

Chronic diseases, age and gender: Examining the contribution to burnout symptoms in a sample of 2075 Canadian workers

Associations of sleep quality and burnout in clinicians during the COVID-19 pandemic

Long working hours and burnout in health care workers: Non-linear dose-response relationship and the effect mediated by sleeping hours-a cross-sectional study

Work hours and difficulty in leaving work on time in relation to work-to-family conflict and burnout among female workers in Taiwan

Exploring the role of shift work in the self-reported health and wellbeing of long-term and assisted-living professional caregivers in Alberta

Self-rated health, work characteristics and health related behaviours among nurses in Greece: A cross sectional study

Social inequalities in self-rated health by age: Cross-sectional study of 22,457 middle-aged men and women

What determines Self-Rated Health (HEALTH STATUS RATINGS)? A cross-sectional study of SF-36 health domains in the EPIC-Norfolk cohort

Validity of the five-item Brief Symptom Rating Scale among subjects admitted for general health screening

Development and verification of validity and reliability of a short screening instrument to identify psychiatric morbidity

The five-item Brief-Symptom Rating Scale as a suicide ideation screening instrument for psychiatric inpatients and community residents

The conceptualization and measurement of perceived wellness: Integrating balance across and within dimensions

Population of the administrative area of the Republic of China

Statistical Power Analysis for the Behavioral Sciences 2nd edn

Survey of perception of safety and health in the work environment in 2016 Taiwan; ILOSH105-M309. Institute of Labor, Occupational Safety, Ministry of Labor

Classification of age groups based on facial features

Night shift and rotating shift in association with sleep problems, burnout and minor mental disorder in male and female employees

Occupational stress, health and occupational burnout among dentists: A study of clinical dentists in Kaohsiung

Available at: FAQ/effectSize-CBU statistics Wiki (cam.ac.uk) Accessed

Jmp For Basic Univariate And Multivariate Statistics: A Step-by-Step Guide

Thanks go to all the hospital staff who participated in the study and cooperated in completing the questionnaires.

M.T.T. wrote the main manuscript text and prepared the table and figures. All authors reviewed the manuscript.

The author declares no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-11326-1.Correspondence and requests for materials should be addressed to M.-T.T.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.