{'fraction_non_alphanumeric': 0.03447076310047121, 'fraction_numerical': 0.01976534581449037, 'mean_word_length': 5.103454931480409, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Predictors of Shift Work Sleep Disorder Among Nurses During the COVID-19 Pandemic: A Multicenter Cross-Sectional Study

Background: Nurses have a high incidence of shift work sleep disorder, which places their health and patient safety in danger. Thus, exploring the factors associated with shift work sleep disorder in nurses is of great significance in improving their sleep health, nursing personnel staffing, and scheduling during the COVID-19 pandemic. Objectives: The purpose of this study was to investigate the incidence of shift work sleep disorder during the COVID-19 pandemic and explore the factors associated with shift work sleep disorder in Chinese nurses. Methods: This was a multicenter cross-sectional study using an online survey. Stratified cluster sampling was used to include 4,275 nurses from 14 hospitals in Shandong, China from December 2020 to June 2021. Stepwise multivariate logistic regression analysis and random forest were used to identify the factors associated with shift work sleep disorder. Results: The prevalence of shift work sleep disorder in the sampled shift nurses was 48.5% during the COVID-19 pandemic. Physical fatigue, psychological stress, shift work more than 6 months per year, busyness during night shift, working more than 40 h per week, working more than four night shifts per month, sleeping more than 8 h before night shift, using sleep medication, irregular meals, and high-intensity physical activity were associated with increased odds of shift work sleep disorder. Good social support, good work-family balance, napping two or three times per week, resting more than one day after shifts, intervals of 8 days or more between shifts, and taking turns to rest during the night shift were associated with decreased odds of shift work sleep disorder. Conclusions: Shift work sleep disorder may be associated with scheduling strategies and personal behavior during the COVID-19 pandemic. To reduce the incidence of shift work sleep disorders in nurses, nursing managers should increase night shift staffing, extend rest days after shift, increase night shift spacing, and reduce overtime, and nurses need to seek more family and social support and control their sleep schedules and diet.


INTRODUCTION
The National Sleep Foundation defines shift work as a work that takes place outside the traditional daily schedule of 09:00-17:00 (1). Shift work is very common in the healthcare sector, in which nurses are usually involved. A previous longitudinal study showed that shift work schedules had no effect on lifestyle (2). However, shift work interferes with biological rhythms, and can lead to a wide variety of health problems, with sleep problems being the most common health problem (3). Thus, shift work sleep disorder (SWSD) is considered a chronic condition that is directly related to work schedule (4). Previous studies have shown that nurses are taking more frequent shifts, with heavier workloads and a sharp increase in fatigue and mental stress during the COVID-19 pandemic (5). Large numbers of patients, long working hours and the definition of new roles for healthcare workers in the context of the COVID-19 pandemic have led to insufficient breaks and sleep, and reduced job satisfaction (6). However, previous studies have lacked reports on the prevalence of SWSD for nurses during the COVID-19 pandemic. Therefore, nursing managers and nurses themselves should improve their awareness of the incidence and harm of SWSD during the COVID-19 pandemic, and pay more attention to the possible factors associated with SWSD in the physical, psychological, environmental, occupational, social, and behavioral aspects, so as to improve professional identity and happiness and ensure the health and safety of nurses on duty.
In 2014, the third edition of ICSD (ICSD-3) defined SWSD as a circadian rhythm disorder and updated the ICSD-2 criteria for this condition (7,8). Major changes from ICSD-2 include linking insomnia/sleepiness to a reduction in the total sleep time associated with work schedules. In addition, the minimum duration of symptoms increased from 1 month to at least 3 months (9). The incidence of SWSD among nurses is reported to be between 24.4 and 37.6% using ICSD-2 criteria, accounting for a quarter of the total number of nurses (10,11). In a Norwegian study using ICSD-3 criteria, 45.6% of nurses reported previously experiencing SWSD during 2 years of follow-up (4). However, a recent study in Finland reported that using ICSD-3 criteria led to a lower incidence of SWSD diagnosis compared to using ICSD-2 criteria (12). Thus far, few studies have investigated the prevalence of SWSD according to the latest criteria, and its prevalence according to country and occupation remains to be determined (13).
SWSD is associated with serious adverse consequences for individuals, employers, and society, and can damage health (1). In addition, there are many health problems that occur with SWSD, such as anxiety or depression (14), headaches (15), and impaired cognitive function (16). SWSD is associated with lower alertness, reduced attention, and increased likelihood of sleepiness-related accidents, all of which can reduce work performance, patient safety, and quality of life (13,(17)(18)(19). In addition, SWSD leads to loss of productivity and increased economic costs (20). The dissonance between endogenous circadian rhythms and sleep and work schedules of shift workers may be the mechanism responsible for SWSD (17). In addition, short rest periods (<11 h) between two shifts have also been reported to be associated with the presence of SWSD (21). Previous reviews have also pointed to the importance of investigating the association of shift work schedules and lifestyle factors with SWSD and exploring preventive measures for SWSD (1,20). Jang proposed interventions to manage SWSD, including adjusting the shift work schedules and working hours, providing shift workers with adequate sleep health education, and providing the opportunities to take naps and breaks during night shifts (20). However, specifics of the factors and interventions, including periods between two shifts, staffing, frequency and duration of naps, sleep duration before and after night shifts, and dietary adjustments, remain unclear (13,17).
Therefore, in this study we first aimed to measure the incidence of SWSD in nurses based on the latest diagnostic criteria. Second, we aimed to explore the comprehensive factors of SWSD in shift nurses from physical, psychological, environmental, occupational, social, and behavioral perspectives based on the Integrated Framework for Population Health Risk Management (22) to fill the gaps in extant research. Furthermore, random forest was used to rank the importance of the factors. Currently, nursing management emphasizes on making full use of strengths and avoiding weaknesses, and our study intended to find approaches that will aid in scientifically and effectively organizing shift schedules for nurses to reduce the incidence of SWSD, promote health, and improve the quality of life of shift nurses.

Design, Setting, and Participants
In this multicenter cross-sectional study, we used data from the first survey of Chinese nurses conducted between December 2020 to June 2021 from the Nurses' Health Cohort Study of Shandong. The cohort study mentioned above is an ongoing prospective cohort study of Chinese nurses conducted in Shandong province, China (registration number: ChiCTR2100043202). The study adopted a multistage sampling method for a more extensive sample coverage. The steps are as follows: (a) Shandong was categorized into eastern, western, southern, northern and central regions; (b) hospitals were selected by convenience sampling (at least one secondary and one tertiary hospital were selected from each region); and (c) participants were selected through cluster sampling.
Registered nurses with nurse qualification certificates who worked in shifts, and agreed to participate in this study were included in the present study. Retired nurses, refresher nurses, and student nurses; nurses who had been working for a duration of <6 months; and nurses who were on leave during the investigation were excluded from the present study.

Measures
The questionnaire consisted of the following four parts: demographic and lifestyle variables; SWSD assessment; shift work characteristics; standardized questionnaires.

Demographic and Lifestyle Variables
The basic demographic variables included sex, age, marital status, children, educational background, department, hospital level, professional title, and income. The body mass index was calculated as body weight (kg) divided by the squared height (m 2 ), and nurses were categorized on the basis of their body mass index as underweight (≤18.4), normal weight (18.5-23.9), and overweight (≥24) (23). Lifestyle variables included drinking caffeinated beverages at work (yes/no), sleep duration before and after the night shift, use of sleep-promoting drugs (yes/no), changes in diet during shift work, and the frequency and duration of naps.

Shift Work Characteristics
Shift work characteristics included shift schedule (two shifts/three shifts/other), direction of shift rotation (clockwise/counterclockwise), intervals between night shifts, number of years working in shift work, number of hours of working/overtime per week, number of months of shift work per year, number of night shifts per month, busyness during night shift (yes/no), psychological burden during shift work (yes/no), and physical discomfort during shift work (yes/no).

SWSD Assessment
SWSD was assessed with the following three questions: (a) Do you have a work schedule that sometimes overlaps with the time you usually sleep? (b) If yes, does this cause insomnia and/or excessive sleepiness due to reduced amount of sleep? (c) If yes, has this lasted for at least three months? Participants who responded "yes" to all three questions were classified as having SWSD (4).

Work-Family Balance Scale
The work-family balance was assessed using a 14-item workfamily balance scale, which is widely used to measure the workfamily inter-relation based on four domains, namely, workfamily conflict, family-work conflict, work-family enrichment, and family-work enrichment (24). This scale was adopted and modified by Chinese researchers and had good reliability and validity (25). Each item is rated on a 5-point Likert scale from 1 (definitely disagree) to 5 (definitely agree). Based on the total score, work-family balance was recategorized as high (>42) and low (≤42). In the present study, the Cronbach's alpha coefficient of the scale was 0.811.

Perceived Social Support Scale
Social support was assessed using the 12-item perceived social support scale, which is a reliable instrument for measuring perceived support from family, friends, and others (26). Responses for each item range from 1 (definitely disagree) to 7 (definitely agree). Based on the total score, shift nurses were categorized as having low , medium (37-60), and high (61-84) levels of perceived social support. The Chinese version has shown good psychometric properties with a Cronbach's alpha coefficient of 0.91 and a test-retest reliability of 0.85 (27). The Cronbach's alpha coefficient of the scale in this study was 0.980.

Fatigue Scale
The fatigue scale was used to measure the severity of fatigue (28). This 14-item questionnaire includes two dimensions: physical fatigue and mental fatigue. Nurses were categorized into different levels of fatigue on the basis of their total score on the scale: low (0-4), medium (5)(6)(7)(8)(9), and high (10)(11)(12)(13)(14). The Chinese version of the scale has been widely used with had good reliability and validity (29). The Cronbach's alpha coefficient of the scale in this study was 0.803.

Circadian Type Inventory
A 11-item circadian type inventory was previously developed to measure circadian rhythms (32). In the current study, we used the Chinese version of the scale to assess circadian flexibility and languidity (33). The two subscales used in this scale were divided into low, medium, and high levels based on the onethird and two-third cutoff points of the theoretical subscale score. A higher level on the flexibility subscale indicated greater adaptability to shifts, while a higher level on languidity indicated greater difficulty in overcoming sleepiness. The Cronbach's alpha coefficient of the scale in this study was 0.838.

International Physical Activity Questionnaire
The international physical activity questionnaire has been widely used to study the physical activity level (34). The questionnaire includes five parts: occupation, housework, transportation, leisure, and sitting. We used the Chinese version of the questionnaire, which has good reliability and validity (35). We divided the physical activity into three levels by calculating the total metabolic equivalent: low, medium, and high.

Data Collection
Baseline information was collected through the Wenjuanxing network platform between December 2020 to June 2021. The survey included eight electronic questionnaires, each of which would take about 10-15 min to complete. Nurses were sent these questionnaires via WeChat (a widely used instant messaging app used in China) in batches to avoid impatience. Subjects who did not participate in the survey within 1 week of sending the questionnaires received a WeChat reminder to complete the survey once per week until the questionnaires were completed. Liaisons were set in every hospital involved in this study, functioning as regulators for questionnaire completion and communicators with the researchers when facing unsolved problems. Only professionals who signed nondisclosure agreements had direct access to the data. To ensure data quality, we adopted the required questions and data logic control design for the electronic questionnaires.

Ethical Considerations
This study was carried out in accordance with the International Guidelines of Good Epidemiology Practice and the Declaration of Helsinki principles. The nurses' participation was voluntary, and written informed consents were obtained from all the respondents. The study was approved by the Ethics Committee of Scientific Research of Hospital with the registration number (KYLL-202011-085).

Data Analysis
All the variables in this study, including one dependent variable and 43 independent variables, were categorical variables. Descriptive statistics included frequency and percentage for them. The distribution of characteristics between shift nurses with and without SWSD was assessed using the chi-squared test or Fisher's exact test for categorical variables. Bonferroni correction was applied for correcting the multiple test and p < 0.0012 was considered statistically significant for the univariate analysis. Statistically significant variables in the univariate analyses were all considered for inclusion in a forward and backward stepwise multivariate logistic regression analysis to identify the main factors associated with SWSD and to estimate the adjusted odds ratios (ORs) and 95% confidence intervals (CIs). The final logistic regression model with Akaike Information Criterion (AIC) minimization was used, and variables in the final model with p < 0.05 denoted statistically significant differences. Logistic regression was used because this classic classifier has been widely used in disease prediction and it is easy to interpret the results.
In order to more accurately explore the importance of factors, we adopted a machine learning method, random forest. The random forest is a tree-based ensemble classification algorithm, which has been widely used to build disease prediction models with small variable sets and high generalization ability. We use grid search for hyperparameter adjustment when training the model. The data was divided into 7:3 training data and test data. The training set was used to train and obtain the final model, and the test set was used to draw receiver-operating characteristic (ROC) curves to determine the prediction performance of the model. The random forest was implemented using Python 3.7.4 (Python Software Foundation, USA) with Scikit-learn (Packt Publishing, UK). Other statistical analyses were carried out using R software version 4.0.5 (R Development Core Team, Vienna, Austria).

Respondent Characteristics and SWSD Prevalence
The participants (n = 4,655) of the present study were recruited from nurses participating in shift work in 14 hospitals across six cities in Shandong, China. Of these, nurses who did not complete the questionnaire or those whose questionnaires contained logical inconsistencies were eliminated, leading to the exclusion of 8.17% of respondents, leaving 4,275 participants. The study participants ranged in age from 21 to 55 years and 93.2% (3,985) were women. The average age of the respondents was 31.19 ± 5.48 years, and the average working duration was 8.62 ± 5.16 years. In this study, the prevalence of SWSD among participating nurses was as high as 48.5%. Table 1 shows all the respondents' features and univariate analysis of SWSD in detail.
Stepwise Multivariate Logistic Regression Analysis of Factors Associated With SWSD Table 2 shows the final forward and backward stepwise multivariate logistic regression model with adjusted OR and 95% CI of the main variables. In terms of physical and psychological factors, compared to nurses who reported a low level of fatigue, those who reported a moderate or high level of fatigue were significantly more likely to experience SWSD (OR = 1.35, 95% CI: 1.08-1.67; and OR = 2.50, 95% CI: 2.02-3.10, respectively). A high perceived stress level (OR = 4.63, 95% CI: 2.60-8.63) had a statistically significant association with SWSD. Additionally, shift nurses who experienced physical discomfort during their night shift (OR = 2.05, 95% CI: 1.74-2.43) or experienced psychological stress before or after the night shift (OR = 1.51, 95% CI: 1.14-2.01; and OR = 1.52, 95% CI: 1.30-1.77, respectively) were more likely to develop SWSD than shift nurses who did not have such experiences.
In terms of social and behavioral factors, nurses who had a good work-family balance (OR = 0.41, 95% CI: 0.29-0.56) or good social support (OR = 0.74, 95% CI: 0.60-0.91) were less likely to experience SWSD. Nurses who responded that they could get >8 h of sleep before night shift were statistically significantly more likely to have SWSD (OR = 1.95, 95% CI: 1.08-3.55) than those who got <3 h of sleep before night shift. In addition, napping 2-3 times per week was a protective factor for SWSD (OR = 0.77, 95% CI: 0.63-0.94). Shift nurses who often

Feature Importance Ranking of the Factors Associated With SWSD on the Random Forest Model
Before analyzing the predictors, we evaluated the performance of the two models. The results showed that the predictive efficiency  (Figure 1). Therefore, the model prediction performance were guaranteed in our study.
To illustrate the contribution of each factor to predictive capacity, we calculated the feature importance based on random forest. The feature importance of the top 50% of factors was visualized in Figure 2, including fatigue, physical discomfort, working overtime per week, using sleep medication, mental stress, diet, length of rest days, work-family balance, and departments. We found that the results of the two statistical methods were consistent, which can indirectly prove the robustness of the results of the logistic regression model.

Prevalence of SWSD Among Nurses
The incidence of SWSD among participating nurses was 48.5%, which is higher than the rates of 37.6% (11) and 35.2% (4) reported in two studies in Norwegian nurses, and the rate of 24.4% reported by The Nurses' Sleep Health Project in Japan (10). Previous studies on the prevalence of SWSD in different populations have reported a prevalence of 23.3% in Norwegian oil rig shift workers (36), 19.9% in Chinese textile shift workers (37), 12.7% in Australian shift workers (38), and 8.1% in American shift workers from the community (39).
This discrepancy in the prevalence among the different studies could be attributable to the different shift schedules of the target populations, different national conditions or different times. Although China and Japan have similar nurse shift schedules, we collected data from 14 hospitals of different levels in areas of different economic levels, so as to more comprehensively reflect the prevalence of SWSD in nurses. Moreover, the incidence of SWSD may be higher since the time of our survey was the COVID-19 pandemic. Our result was somewhat supported by a systematic review stating that 36.36% of healthcare workers experienced work-related sleep problems during the COVID-19 pandemic (40). In addition, differences in disease definitions or questionnaire content might also reflect this divergence. In 2009, some scholars developed a self-administered questionnaire based on the ICSD-2 (36). In 2014, ICSD-3 updated the criteria for SWSD, and Waage et al. updated the questionnaire content according to the updated criteria (4). In ICSD-2, SWSD is defined as insomnia and/or excessive sleepiness that overlaps with shift work hours and cannot be explained by other diseases (7). In ICSD-3, the criteria for SWSD included a reduction in total sleep time as a symptom and increased the minimum duration of symptoms associated with shift work scheduling from 1 to 3 months (8). ICSD-3 was more rigorous in terms of duration and total sleep time was added as an indicator to make the diagnosis more sensitive, leading to significant differences between the two criteria in estimating the incidence of SWSD. Unlike most previous studies based on ICSD-2 (10,11,36), we used the latest questionnaire based on the ICSD-3 criteria to obtain a more rigorous and sensitive diagnosis of SWSD.

Physical and Psychological Factors
Similar to the results of previous studies, characteristics of physical and psychological factors showed a statistically significant relationship with the prevalence of SWSD among nurses; however, there were some differences in specific variables. For instance, SWSD is essentially a circadian rhythm disorder that has been associated with polymorphisms in various molecular clock genes (41). Languidity and flexibility are two characteristics of the circadian rhythm, and several previous studies have shown that fatigue may influence SWSD, while flexibility is a protective factor (19,42,43). Contrary to our hypothesis, we found no relationship between SWSD and circadian rhythm characteristics. One possible reason for this unexpected finding could be that approximately 70% of the circadian rhythms of shift nurses measured by the questionnaire in this study were at a moderate level, suggesting that a more objective and sensitive measurement method could help explore the relationship between SWSD and circadian rhythm characteristics in greater detail. Furthermore, in our study, nurses who reported experiencing fatigue and physical discomfort (e.g., headache, heart discomfort, gastrointestinal symptoms) during their night shifts were more likely to develop SWSD than those who did not, supporting the results reported previously (15,23,42,44). Interestingly, high work stress during shift work and difficulty in falling asleep were previously found to be the typical symptoms of SWSD (45), consistent with our findings. In addition, the sensitivity of sleep to stress is a known predictor of SWSD (19,46,47), suggesting that reducing the mental burden of nursing work has potential benefits in reducing the incidence of SWSD. In addition, Afshari et al. reported that nurses have low levels of psychological resilience during the COVID-19 pandemic (48). Therefore, it is necessary to improve psychological resilience to relieve stress levels and indirectly reduce the impact of stress on SWSD.

Environmental and Occupational Factors
Haile et al. proposed that three shifts are statistically significantly related to SWSD (49), which is different from our finding. We found that nurses on two shifts had a higher prevalence of SWSD than nurses on three shifts, although the results of the stepwise multivariate logistic regression analysis showed no statistically significant correlation between SWSD and rotation schedule, which supported the results of Vanttola et al. (19). Thus, we will continue to investigate the effect of shift schedule on SWSD with prospective follow-up data to verify the current results. An interesting finding of our study was that the number of shift years had no effect on SWSD, but more than 6 months of shift scheduling per year was associated with a higher likelihood of having SWSD. In addition, SWSD was correlated with the average number of nights per month, and more than four night shifts were detrimental to nurses' sleep, which is similar to that reported in other studies (21,49). This suggests that in nursing management, the head nurse should schedule no more than four night shifts per month to prevent SWSD among nursing staff. Short rest periods (<11 h) between work shifts have been shown to increase the likelihood of SWSD in nurses (44,50). This supports our finding that more than 1 day of rest (>12 h) after a night shift was a protective factor for SWSD in our study. Moreover, our study reported a novel finding that a busy night shift with one nurse can easily lead to SWSD, but a night shift where two nurses take turns to rest reduces the odds of SWSD. An interval of 8 days or more between night shifts helps nurses adapt to shifts and reduces the odds of SWSD. Working overtime (>40 h/week) for more than a week per month affected sleep. Di Milia et al. also found a statistically significant association between SWSD and weekly work hours (43).

Social and Behavioral Factors
In our study, the marital status of nurses and whether or not they had children had no association with SWSD, but a higher workfamily balance or social support (i.e., family support and friend support) was associated with a lower likelihood of SWSD. Among nurses, the level of social support was statistically significantly and negatively associated with insomnia and daytime sleepiness (51), which supported our conclusions. Our results differed from those of previous studies (49), which only took into account marital status and children, possibly ignoring the more important variable of work-family balance. Afshari et al. reported that nurses have increased demands for social support and workfamily balance in the context of the COVID-19 pandemic (48). Therefore, future research can focus more on the influence of psychological factors related to family and society on SWSD.
Numerous studies have repeatedly identified the use of multiple hypnotic agents (such as exogenous melatonin and injectable sleep medications) as predictors of SWSD (21,49,52), which is consistent with the findings of our research. We also investigated the role of sleep time before and after the night shift, rest taken during the night shift, and weekly naps. Those with SWSD had greater sleep debt on shift days and slept less compensatory sleep on days off than those without SWSD (19). Hence, our results suggest that effective sleep supplementation should include sleeping for <8 h before the night shift, with 2-3 compensatory naps per week. It is vital to raise awareness and education on sleep hygiene practices among nurses (53). We innovatively explored the effects of diet on SWSD during shifts. In our study, we found that irregular eating during the shift was a predictor of SWSD, but the results should be interpreted with caution and studies should be carried out in the future to verify our findings. Our results show that highintensity physical activity was positively correlated with SWSD. However, a previous study contradicted our findings (54). One possible reason for the contradictory results is that we used a selfadministered questionnaire to assess physical activity whereas Hajo et al. used objective ActiGraph GT3X accelerometers for data collection (54,55). Therefore, future studies can combine self-reporting and objective monitoring to further explore the effect of physical activity on SWSD.

Limitations
One limitation is the cross-sectional nature of this analysis. Further studies should collect longitudinal data to better understand the causal relationship between factors and SWSD. In addition, this study only included hospitals in Shandong Province in China; therefore, caution should be exercised when generalizing the study's findings to other regions. Our study lacks objective indicators; thus, future studies with objective measurements will help to verify the accuracy of the results of the current study.

CONCLUSION
The prevalence of SWSD in the sampled shift nurses was 48.5% during the COVID-19 pandemic. To reduce the incidence of SWSD for nurses, head nurses and health care organizations need to develop more rational and scientific scheduling strategies and health education for nurses on the basis of the Integrated Framework for Population Health Risk Management. For example, fatigue, perceived stress, and body discomfort were among the physical and psychological factors associated with SWSD. Therefore, nursing managers should pay attention to the physical and mental health of nurses and offer timely interventions. In terms of environment and occupation, good scheduling strategies included placing two nurses on night shifts and allowing them to take turns to rest, arranging for shift nurses to rest for more than 1 day after night shifts, and controlling shift intervals for more than 8 days. However, arranging nurses in continuous shifts for more than 6 months, scheduling more than four night shifts per week, and working overtime were all bad scheduling strategies. Positive social and behavioral factors included 6-8 h of supplementary sleep after the night shift, 2-3 naps per week, eating regularly and ensuring adequate water intake per day, and improving family and social support. Therefore, health education for nurses should include good sleep strategies and maintenance of good family and social relations.

DATA AVAILABILITY STATEMENT
The datasets presented in this article are not readily available because the data from the Nurses' Health Cohort Study of Shandong needs time for data clearing and establishment of guidelines. We are planning on opening this data to the public in the future. Requests to access the datasets should be directed to caoyj@sdu.edu.cn.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Ethics Committee of Scientific Research of Shandong University Qilu Hospital. The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
YL and YW: methodology, formal analysis, data curation, software, writing-original draft, and visualization. RL and XL: writing-review and editing and project administration. XG, LL, and JL: investigation. FX, XJ, and YC: conceptualization, resources, supervision, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.

FUNDING
The study was funded by the National Key Research and Development Program of China (Grant number 2020YFC2003500). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.