- 1Department of Neurology, NO 984 Hospital of PLA, Beijing, China
- 2Department of Medical Psychology, NO 984 Hospital of PLA, Beijing, China
- 3Department of Medical Psychology, Army Medical University, Chongqing, China
- 4Department of Psychiatry, NO 984 Hospital of PLA, Beijing, China
Purpose: Rest-activity rhythms (RAR) is one of the most fundamental biological patterns that influence basic physical and mental health, as well as working performance. Nowadays, with the utilization of actigraphy device in detecting and analyzing RAR objectively, evidence has confirmed that RAR could be interrupted by shift work. Nurses need to participate shift work in their daily routine. The aim of the present study is to identify the RAR of nurses working rotating shifts and find out the relationship between RAR and executive function.
Methods: A total of 11 female nurses on day-shift (DS), 11 female nurses on rotating-shift (RS) were recruited. Demographic features, non-parametric variables of RAR as well as sleep variables according to actigraphy data, and a series of cognitive function evaluations during work time were collected.
Results: The Student’s t-test comparison revealed that RS group nurses showed significant shorter total sleep time (TST), lower level of interdaily stability (IS) and relative amplitude (RA), as well as poorer executive performance compared with DS nurses. The linear regression analysis demonstrated that several aspects of executive performance such as choice reaction test (CRT) (reflecting attention and concentration) and trail making test (TMT) -B (reflecting cognitive flexibility) were associated with IS in RS group nurses.
Conclusion: These results elucidated that RAR was disrupted for nurses working rotating shifts, and executive performance was associated with IS from day to day in nurses working rotating shifts.
Introduction
1A person’s physiological functions change periodically. These changes are known as biological rhythms and are induced by signals from internal molecular mechanisms. In human beings, sleep-wake activity follows a diurnal pattern which is called the rest–activity rhythm (RAR) (Li et al., 2021). In recent years, a variety of methods recording several days of continuous accelerometer data could be used to detect RAR features (Smagula et al., 2019). These methods made investigators measure and quantify RAR characteristics objectively possible. Furthermore, several studies implied that RAR characteristics might be associated with physical and mental health (Slyepchenko et al., 2019; Kume et al., 2020).
Shift work involves working irregular or unusual hours different from a daytime work schedule (Wang et al., 2011). Nurses, participating shift work have jobs which require them to be active during hours that constitute normal rest time for most people. It has been concluded that shift work could cause circadian disruption and sleep loss, resulting in cognitive impairments and even accidents (Kecklund and Axelsson, 2016).
Rest-activity rhythms is considered as a type of rhythm pattern parallel, but not exactly equal to sleep-wake cycle (Calogiuri et al., 2013). As has been mentioned by Calogiuri et al., “adding analysis of the rest-activity rhythm using rhythmometric procedures to traditional actigraphic studies could provide a deeper understanding of possible disorders of sleep in relation to the entrainment and vice versa.” It is reported that shift work and the associated modifications in daily routines affect the nurses’ circadian rhythm (Kang et al., 2015). Though amount of studies have elucidate the relationship between shift work and decline in work efficiency for nurses (Narciso et al., 2016; Merchaoui et al., 2017), whether disrupted RAR of nurses of shift work impact on their executive performance during work is not clear.
In this study, we sought to identify the RAR of nurses working in rotating shifts by analyzing actigraphy data through the lens of non-parametric variables, and find out the relationship between RAR and executive function in nurses working in rotating shifts.
Methods
Participants
From June 1, 2020 to December 31, 2020, 11 female nurses on day-shift (DS), 11 female nurses on rotating-shift (RS) were recruited. Our study was approved by the Academic Ethics Committee of the Biological Sciences Division of NO 984 Hospital of PLA in Beijing, China. Informed consent was obtained from each participant in the study.
Procedure
The nurses worked 2 × four slow, forward-rotating arranged RS (day shift-evening shift-night shift-off) or consecutive 8 days DS (08:00–16:00). Before beginning the actigraphic monitoring, demographic characteristics were collected. The exclusion criteria were: nurses who was pregnant or feeding baby younger than 3 years; in the habit of alcoholic drinking, strong tea or caffe; using of sedatives or hypnotic drugs; suffering from metabolic diseases, neurological diseases, renal diseases, or mental disorders.
Wrist Actigraphy
According to the procedure, each participant was instructed to wear an ActiGraph GT3X+ device (ActiGraph, Pensacola, FL, United States) on their non-dominant wrist for 24 h per day (except when bathing or swimming) for 8 days. Considering that longer wearing of actigraphy duration will impact on the compliance of nurses, we selected 8-day protocol. And 8-day (encompassing 2 day shift+2 evening shift+2 night shift+ 2 off) is the minimum possible to observe with some clarity the RAR in people with unconventional routines.
At the end of the wear period, data were downloaded using ActiLife software (ActiGraph, Pensacola, FL, United States). All data files were visually screened for sufficient wear time and then processed for analysis (Figure 1).
Figure 1. An example of a participant’s actogram, raw actigraphy data along with rest-activity rhythm (A) and sleep intervals (B) across the 24-h day.
Data Processing
We calculated the non-parametric features using actigraph count data (vector magnitude counts) as follows: (1) the interdaily stability (IS); (2) the intradaily variability (IV); (3) the time occurrence and corresponding activity counts of the most active 10 h period (M10) and of the least active 5 h period (L5); and (4) the relative amplitude (RA) (Van Someren et al., 1999).
The IS, which provides information about the RAR synchronization to supposedly stable environmental stimuli, is calculated from the mean 24 h profile as the following formula 1:
The IV, which provides information about the fragmentation in the RAR over the 24 h profile, calculated from hourly raw data as the following formula 2:
In equation of IS or IV, n is the total number of data, p is the number of data entries per day, −xh is the hourly average, −x is the average of all the data and−xi represents each hour of raw data. The IS ranges from 0 for Gaussian noise, which consisted of normal distribution noise, to 1 for a perfect synchronization. The IV measure ranges from 0 to 2, whose higher value indicates a more fragmented rhythm.
M10 is defined as the maximum sum of 10 consecutive hours of activity log. L5 is defined as the lowest sum of 5 consecutive hours of the activity log. RA is calculated as (M10-L5)/(M10+L5) (Witting et al., 1990).
Objective Sleep Variables Detection
Sleep variables were also collected using a triaxial accelerometer (GT3X+; ActiGraph, LLC, Pensacola, FL, United States); this is a small (4.6 × 3.3 × 1.5 cm, 19 g), easy-to-wear device that offers the possibility of recording and scoring sleep–wake patterns during the night (Kim et al., 2015). Data from the Actigraphy data were downloaded and analyzed using the software with a 60 s epoch. Bedtime and wake time from the sleep diary were used to define rest intervals for calculation of sleep-wake variables. Actigraphy variables shown here were sleep latency (SL), sleep efficiency (SE), total time in bed (TTB), total sleep time (TST), wake after sleep onset (WASO), times of awakenings (TA) and average duration of awakenings (ADA).
Cognitive Function Assessment
During the last DS, or night-shift work, similar to the procedure of Chang et al. (2011), all participants of each group completed a series of cognitive function tests, respectively. Considering that the main duties of nurses working night shifts at the study hospital include checking medical orders and prescriptions (Chang et al., 2011) correctly and quickly, the cognitive function assessments in the current study consisted of auditory verbal learning test-HuaShan version (AVLT-H), choice reaction test (CRT), digit symbol substitution test (DSST) and trail making test (TMT)-B (Zhào et al., 2019).
Auditory verbal learning test-HuaShan version is adapted by the Guo’s group (Yang et al., 2018) team based on the California Verbal Learning Test (CVLT) and Hong Kong auditory verbal learning test. This is a well-recognized measure for testing a person’s ability in different stages of working memory which includes encode, combine, store and recover verbal information (Zhang et al., 2019). Specifically, we informed the subjects to recall the words in advance. Then the neurologist read 12 Chinese words, including clothing, occupation, and flowers. We chose delayed memory mode, and the subject was asked to recall the words at the time point of 20 min after the words reading-out. CRT is a measurement of attention and concentration need to be completed on screen (Hou et al., 2020). During the test, stimulus (image of an arrowhead) was presented 40 times at random intervals on the screen. Participants were required to respond by pressing the button as quickly as possible with either hand. DSST is an evaluation of processing speed, derived from the battery of THINC-it (Hou et al., 2020). Participants are presented with six symbols, each representing one of six (1–6) digits. Subjects are asked to assign as many symbols as possible to the respective digits within 120 s. TMT-B is an assessment tool reflecting visual search speed and cognitive flexibility. It involves the participant connecting a trail through symbols (numbers 1, 2, 3……9 and corresponding Chinese characters 壹, 贰, 叁……玖) in the correct order (Wei et al., 2019). The time recorded began once the instruction to start was given and was stopped once the trail was completed.
Statistical Analyses
The differences between the groups demographic data, as well as the differences between the groups on sleep variables, Non-parameter variables of RAR, CRT, DSST, and TMT-B performance were analyzed by Student’s t-test. A linear regression analysis was selected to find out the relationship between RAR (IS, IV, and RA), sleep variables (SL, SE, TTB, TST, WASO, TA, and ADA) and each domain of executive function. The data were expressed as the mean ± standard deviation, and a P value <0.05 was considered statistically significant. All statistical analyses were carried out under the statistical software package SPSS, 22.0 (IBM Corp., Armonk, NY, United States).
Results
Nurses of both groups did not show significant difference in age, weight, height and body mass index (BMI). Most sleep variables except TST (357.34 ± 64.84 min vs. 414.59 ± 43.84 min, P = 0.025) did not reach statistical significance between RS group nurses and DS group nurses. On the contrary, non-parametric features of RAR such as IS (0.35 ± 0.10 vs. 0.59 ± 0.13, P = 0.000) and RA (0.80 ± 0.13 vs. 0.94 ± 0.02, P = 0.007), rather than IV (0.92 ± 0.18 vs. 0.91 ± 0.11, P = 0.933), differed statistically between group subjects. Details were shown in Table 1.
In addition, we evaluated the cognitive function of nurses in different groups. The RS group nurses spent significant longer time in completing CRT (0.38 ± 0.06 s vs. 0.32 ± 0.05 s, P = 0.045) and TMT (20.64 ± 3.01 vs. 17.55 ± 2.84, P = 0.022) relative to DS group nurses. Whereas, the scores on AVLT-H (9.72 ± 1.49 vs. 9.91 ± 1.30, P = 0.764) and DSST (23.00 ± 2.05 vs. 23.73 ± 1.95, P = 0.404) did not show obvious differences between groups. Details were shown in Table 1.
Lastly, the association between cognitive function and RAR variables were detected using linear regression analysis in RS group nurses. The CRT time was negatively associated with IS (P = 0.031, standardized β = −2.690). TMT-B time was negatively associated with IS (P = 0.004, standardized β = −4.214). The scatter graphs demonstrating the relationship between IS and CRT, as well as IS and TMT-B, were shown in Figures 2A,B, respectively. Neither AVLT-H score nor DSST score was associated with RAR variables. However, cognitive function did not correlate with sleep variables in the present study. Details were shown in Table 2.
Figure 2. Scatter graphs demonstrating the relationship between IS and CRT (A), as well as IS and TMT-B (B).
Table 2. Association between cognitive functional data and RAR non-parametric variables in RS group nurses.
Discussion
In the present study, RS group nurses showed disrupted RAR and poorer sleep quality, reflected by lower IS and RA, as well as shorter TST in comparison with DS group nurses, respectively. In the view of RAR findings, as far as we know, this is the first study revealing the reduced stability of activity patterns from day to day and decreased amplitude of RAR in nurses working rotating shifts. In terms of sleep quality findings, our findings are in consistent with the systemic review reported by Chang and Peng (2021), showing that the sleep quality of nurses working rotating shifts was poorer than that of nurses working fixed day shifts. The only interesting point was that we only found statistical difference in TST between groups. However, unlike other studies, it seemed that rotation shift did not obviously change the levels of other sleep variables except TST in the current study. There are several explanations for this phenomenon: (1) compared with fixed night shift, rotation shift was an alternative reducing the amount of night work as possible (Khan et al., 2021); (2) it has been evidenced clockwise rotation of shift work could lead to better sleep quality relative to other methods (Barton and Folkard, 1993); and (3) sleep quality was evaluated using subjective methods such as Pittsburgh Sleep Quality Index in previous researches (Chang and Peng, 2021), while objective sleep quality measurement based on actigraphy data was utilized in our study. Considering that rotating shift workers are susceptible to mental disorders, particularly anxiety and depression (Kalmbach et al., 2015), it could be inferred that subjective sleep quality evaluation might overestimate the impact of shift work on sleep quality in nurses.
Apart from RAR and sleep variables, we also evaluate cognitive function of nurses in both groups during their working period. In detail, RS group nurses show deficits in executive function (attention and cognitive flexibility), with working memory and processing speed relatively spared. Previous findings of Özdemir and colleagues (Özdemir et al., 2013) implied that shift work could influence cognitive function in more domains. This discrepancy might be due to the study design, since we chose to collect cognitive function data during working period, instead of a fixed time 8:00 am (within working period for day shift nurses, but out working period for night shift nurses). Furthermore, linear regression analysis in RS group nurses demonstrated that poorer executive function was associated with less stable day-to-day RAR, instead of sleep quality. This finding is similar to the results derived from early Parkinson’s disease (Wu et al., 2018).
Interestingly, the impact of shift work on nurses seemed to differ according to the characteristics of the cognitive tests in the current study. This is different from other previously published studies demonstrating obvious declines in memory tests (Thun et al., 2021) and DSST (Zion and Shochat, 2018). This disaccord can be explained by the linear regression results. For example, Thun et al. (2021) reported that memory problems were associated with quick return and short sleep duration, and Zion and Shochat. (2018) found that DSST score was associated with clock time and sleepiness. All these parameters were relevant to fast rotating shift, instead of slow rotating shift. Actually, there existed evidence that unstable daily routine was associated with both memory and executive function in patients with liver cirrhosis (Kim et al., 2021) and Alzheimer’s disease (Alfini et al., 2021). The reason why we only find less aspects of cognition associated with IS might be that the participants recruited in the present study were all young healthy adults. As has been reported previously, cognitive declines began from as early as 30–40 years old and continued into older age (Ferguson et al., 2021). Thus, we inferred young adults might be more tolerant than elderly/patients toward IS disrupted condition.
Combined with the decorrelation between IS and sleep variables, all these results pointed out that RAR (e.g., IS, IV) is not merely a proxy for sleep in its association with cognition.
Several limitations of this study warrant consideration. First, the sample size was small. Second, we did not included nurses working fixed shifts, because rotating working shift is the most commonly found type of shift work for clinical departments in China. Third, we only recruited female nurses, in future study, we would collect male participants to overcome this shortcoming.
In summary, nurses working rotating shifts showed disrupted RAR, low sleep quality and impaired executive function. Poor executive function was associated with less stability of activity patterns from day to day.
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 human participants were reviewed and approved by the Ethics Committee of NO 984 Hospital of PLA. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
HZ and XZ were responsible for data collection. QT was responsible for manuscript writing. ZF was responsible for data analysis. HZ was responsible for the study design. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Foundation (Grant No.: CLB18J026).
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.
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.
Acknowledgments
We thank Shuo Tian for the support of data collection.
References
Alfini, A., Albert, M., Faria, A. V., Soldan, A., Pettigrew, C., Wanigatunga, S., et al. (2021). Associations of actigraphic sleep and circadian rest/activity rhythms with cognition in the early phase of Alzheimer’s disease. Sleep Adv. 2:zab007. doi: 10.1093/sleepadvances/zpab007
Barton, J., and Folkard, S. (1993). Advancing versus delaying shift systems. Ergonomics 3, 59–64. doi: 10.1080/00140139308967855
Calogiuri, G., Weydahl, A., and Carandente, F. (2013). Methodological issues for studying the rest-activity cycle and sleep disturbances: a chronobiological approach using actigraphy data. Biol. Res. Nurs. 15, 5–12. doi: 10.1177/1099800411416224
Chang, W. P., and Peng, Y. X. (2021). Influence of rotating shifts and fixed night shifts on sleep quality of nurses of different ages: a systematic literature review and meta-analysis. Chronobiol. Int. 38, 1384–1396. doi: 10.1080/07420528.2021.1931273
Chang, Y. S., Wu, Y. H., Hsu, C. Y., Tang, S. H., Yang, L. L., and Su, S. F. (2011). Impairment of perceptual and motor abilities at the end of a night shift is greater in nurses working fast rotating shifts. Sleep Med. 12, 866–869. doi: 10.1016/j.sleep.2011.03.018
Ferguson, H. J., Brunsdon, V. E. A., and Bradford, E. E. F. (2021). The developmental trajectories of executive function from adolescence to old age. Sci. Rep. 11:1382. doi: 10.1038/s41598-020-80866-1
Hou, Y., Yao, S., Hu, S., Zhou, Q., Han, H., Yu, X., et al. (2020). Psychometric properties of the Chinese version of the THINC-it tool for cognitive symptoms in patients with major depressive disorder. J. Affect. Disord. 273, 586–591. doi: 10.1016/j.jad.2020.03.146
Kalmbach, D. A., Pillai, V., Cheng, P., Arnedt, J. T., and Drake, C. L. (2015). Shift work disorder, depression, and anxiety in the transition to rotating shifts: the role of sleep reactivity. Sleep Med. 16, 1532–1538. doi: 10.1016/j.sleep.2015.09.007
Kang, J. H., Miao, N. F., Tseng, I. J., Sithole, T., and Chung, M. H. (2015). Circadian activity rhythms and sleep in nurses working fixed 8-hr shifts. Biol. Res. Nurs. 17, 348–355. doi: 10.1177/1099800414549454
Kecklund, G., and Axelsson, J. (2016). Health consequences of shift work and insufficient sleep. BMJ 355:i5210. doi: 10.1136/bmj.i5210
Khan, W. A. A., Jackson, M. L., Kennedy, G. A., and Conduit, R. (2021). A field investigation of the relationship between rotating shifts, sleep, mental health and physical activity of Australian paramedics. Sci. Rep. 11:866. doi: 10.1038/s41598-020-79093-5
Kim, M., Liotta, E. M., Maas, M. B., Braun, R. I, Garcia-Canga, B., Ganger, D. R., et al. (2021). Rest-activity rhythm disturbance in liver cirrhosis and association with cognitive impairment. Sleep 44:zsaa288. doi: 10.1093/sleep/zsaa288
Kim, M., Sasai, H., Kojima, N., and Kim, H. (2015). Objectively measured night-to-night sleep variations are associated with body composition in very elderly women. J. Sleep Res. 24, 639–647. doi: 10.1111/jsr.12326
Kume, Y., Kodama, A., and Maekawa, H. (2020). Preliminary report; comparison of the circadian rest-activity rhythm of elderly Japanese community-dwellers according to sarcopenia status. Chronobiol. Int. 37, 1099–1105. doi: 10.1080/07420528.2020.1740725
Li, J., Somers, V. K., Lopez-Jimenez, F., Di, J., and Covassin, N. (2021). Demographic characteristics associated with circadian rest-activity rhythm patterns: a cross-sectional study. Int. J. Behav. Nutr. Phys. Act. 18:107. doi: 10.1186/s12966-021-01174-z
Merchaoui, I., Bouzgarrou, L., Mnasri, A., Mghanem, M., Akrout, M., Malchaire, J., et al. (2017). Influence of shift work on the physical work capacity of tunisian nurses: a cross-sectional study in two university hospitals. Pan. Afr. Med. J. 26:59. doi: 10.11604/pamj.2017.26.59.11279
Narciso, F. V., Barela, J. A., Aguiar, S. A., Carvalho, A. N., Tufik, S., and de Mello, M. T. (2016). Effects of shift work on the postural and psychomotor performance of night workers. PLoS One 11:e0151609. doi: 10.1371/journal.pone.0151609
Özdemir, P. G., Selvi, Y., Özkol, H., Aydın, A., Tülüce, Y., Boysan, M., et al. (2013). The influence of shift work on cognitive functions and oxidative stress. Psychiatry Res. 210, 1219–1225. doi: 10.1016/j.psychres.2013.09.022
Slyepchenko, A., Allega, O. R., Leng, X., Minuzzi, L., Eltayebani, M. M., Skelly, M., et al. (2019). Association of functioning and quality of life with objective and subjective measures of sleep and biological rhythms in major depressive and bipolar disorder. Aust. N. Z. J. Psychiatry 53, 683–696. doi: 10.1177/0004867419829228
Smagula, S. F., Gujral, S., Capps, C. S., and Krafty, R. T. (2019). A systematic review of evidence for a role of rest-activity rhythms in dementia. Front. Psychiatry 10:778. doi: 10.3389/fpsyt.2019.00778
Thun, E., Waage, S., Bjorvatn, B., Moen, B. E., Vedaa, Ø, Blytt, K. M., et al. (2021). Short sleep duration and high exposure to quick returns are associated with impaired everyday memory in shift workers. Nurs. Outlook 69, 293–301. doi: 10.1016/j.outlook.2020.09.008
Van Someren, E. J., Swaab, D. F., Colenda, C. C., Cohen, W., McCall, W. V., and Rosenquist, P. B. (1999). Bright light therapy: improved sensitivity to its effects on rest-activity rhythms in Alzheimer patients by application of nonparametric methods. Chronobiol. Int. 16, 505–518. doi: 10.3109/07420529908998724
Wang, X. S., Armstrong, M. E., Cairns, B. J., Key, T. J., and Travis, R. C. (2011). Shift work and chronic disease: the epidemiological evidence. Occup. Med. (Lond). 61, 78–89. doi: 10.1093/occmed/kqr001
Wei, W., Zhào, H., Liu, Y., and Huang, Y. (2019). Traditional trail making test modified into brand-new assessment tools: digital and walking trail making test. J. Vis. Exp. 153:e60456. doi: 10.3791/60456
Witting, W., Kwa, I. H., Eikelenboom, P., Mirmiran, M., and Swaab, D. F. (1990). Alterations in the circadian rest-activity rhythm in aging and Alzheimer’s disease. Biol. Psychiatry 27, 563–572. doi: 10.1016/0006-3223(90)90523-5
Wu, J. Q., Li, P., Stavitsky Gilbert, K., Hu, K., and Cronin-Golomb, A. (2018). Circadian rest-activity rhythms predict cognitive function in early parkinson’s disease independently of sleep. Mov. Disord. Clin. Pract. 5, 614–619. doi: 10.1002/mdc3.12692
Yang, K., Shen, B., Li, D. K., Wang, Y., Zhao, J., Zhao, J., et al. (2018). Cognitive characteristics in Chinese non-demented PD patients based on gender difference. Transl. Neurodegener. 7:16. doi: 10.1186/s40035-018-0120-1
Zhang, S., Luo, Y., Dong, Z., Chen, X., Cheng, X., Dong, Q., et al. (2019). Impact of periventricular hyperintensities and cystatin C on different cognitive domains in the population of non-demented elderly Chinese. J. Clin. Neurosci. 68, 201–210. doi: 10.1016/j.jocn.2019.05.053
Zhào, H., Wei, W., Do, E. Y., and Huang, Y. (2019). Assessing performance on digital clock drawing test in aged patients with cerebral small vessel disease. Front. Neurol. 10:1259. doi: 10.3389/fneur.2019.01259
Keywords: shift work, nurse, rest-activity rhythm, non-parametric analysis, executive functions
Citation: Zhao X, Tang Q, Feng Z and Zhào H (2022) Executive Performance Is Associated With Rest-Activity Rhythm in Nurses Working Rotating Shifts. Front. Neurosci. 16:805039. doi: 10.3389/fnins.2022.805039
Received: 29 October 2021; Accepted: 20 January 2022;
Published: 24 February 2022.
Edited by:
Laura Catherine Roden, Coventry University, United KingdomReviewed by:
Guilherme Umemura, University of São Paulo, BrazilBruno Da Silva Gonçalves, Federal University of São Paulo, Brazil
Copyright © 2022 Zhao, Tang, Feng and Zhào. 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: Zhengzhi Feng, fzz@tmmu.edu.cn; Hóngyi Zhào, zhaohongyi17@126.com
†These authors share first authorship