Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 03 December 2024
Sec. Epilepsy
This article is part of the Research Topic Prevention of and early intervention for behavioral health disorders in people with epilepsy View all 3 articles

Circadian rhythm and epilepsy: a nationally representative cross-sectional study based on actigraphy data

Tianyou Tang\nTianyou Tang1 YuDong Zhou\n\n
\nYuDong Zhou1 * Xuan Zhai\n,\n\n
\nXuan Zhai1 ,2 *
  • 1Department of Neurosurgery Children’s Hospital of Chongqing Medical University, Chongqing, China
  • 2Chongqing Key Laboratory of Child Neurodevelopment and Cognitive Disorders, Chongqing, China

Objective: The study aims to assess the relationship between epilepsy and circadian rhythms.

Method: This study included a cohort of 7,410 participants sourced from the 2013–2014 National Health and Nutrition Examination Survey (NHANES) database. The investigation focused on the comparative analysis of seven nonparametric indices associated with circadian rhythms (Interdaily Stability (IS), Intradaily Variability (IV), Relative Amplitude (RA), L5, M10, L5 start time, and M10 start time) between the overall population and patients with epilepsy. Logistic regression analysis was utilized to assess the potential correlation between the rest-activity circadian rhythm patterns and the presence of epilepsy within the cohort.

Results: Compared to the general population, individuals with epilepsy exhibited lower values of IS and M10. Multivariable logistic regression analysis, when IS, RA, and M10 were categorized into four groups based on quartiles, revealed that the odds ratio (IS: OR = 0.36, 95% CI: 0.13, 0.89; RA: OR = 0.25, 95% CI: 0.06, 0.77; M10: OR = 0.24, 95% CI: 0.06, 0.73) for the highest quartile was lower than that for the lowest quartile. Furthermore, after adjustment for confounding factors, participants in the highest quartile compared to those in the lowest quartile of IV and M10 start time demonstrated a higher prevalence of epilepsy.

Conclusion: Individuals with epilepsy demonstrate significant alterations in circadian rhythms.

Highlights

• We transformed the objective accelerometer data into circadian rhythm information.

• We conducted a cross-sectional study using the NHANES database to explore the association between circadian rhythms and epilepsy.

• Individuals with epilepsy exhibit significant changes in circadian rhythms.

Introduction

Epilepsy is a severe neurological disease characterized by recurrent seizures, with a lifetime prevalence ranging from 1 to 3% (1, 2). The circadian rhythm is an intrinsic part of nearly all biological functions, exhibiting a 24-h daily cycle (3). The circadian rhythm of rest-activity is considered a reliable marker of the circadian system (4).

It has long been observed that epilepsy tends to manifest at specific times of the day. Consequently, researchers have endeavored to explore the relationship between the 24-h circadian rhythm and epilepsy. With advancements in technologies such as long-term electroencephalographic recordings, our understanding of the role of the circadian rhythm in epilepsy has been significantly enhanced (5, 6). However, existing studies have not clearly delineated the relationship between epilepsy and the circadian rhythm, and there is also a scarcity of clinical research based on objective data in this regard (7, 8). Liguori et al. included 22 patients with epilepsy and 17 controls, finding that patients with epilepsy exhibited lower synchrony and higher fragmentation of the rest-activity rhythm (9). This is currently the only study that has explored the association between the two using activity data. To address this gap in the literature, we leverage the robust data from the National Health and Nutrition Examination Survey (NHANES), a nationally representative cross-sectional survey database that contains detailed information on disease conditions and various lifestyle factors, thus demonstrating a high degree of generalizability. Specifically, this study quantified accelerometer data from the NHANES 2013–2014 dataset to measure the 24-h rest-activity rhythm in persons with epilepsy, aiming to yield novel insights into this critical area of research.

Method

Participants

Participants were sourced from the NHANES dataset. NHANES is a multi-stage, stratified, nationally representative cross-sectional study (10). The flowchart depicting the participant selection process is presented in Figure 1. Ultimately, a total of 7,410 participants were included. The NHANES protocol has been approved by the Research Ethics Review Board of the National Center for Health Statistics in the United States, and all participants provided informed consent.

Figure 1
www.frontiersin.org

Figure 1. The inclusion of study participants.

Data on rest–activity rhythm

Rest-activity rhythm data were collected from the Physical Activity Monitor (PAM) data. Participants in NHANES aged 3 years and older were instructed to wear ActiGraph GT3X+ devices continuously for seven consecutive days to collect objective data on 24-h movement patterns during both wakefulness and sleep. The PAM utilized in this study captured acceleration data on the x, y, and z axes at intervals of 1/80 s (80 Hz). Minute summary files and daily summary files were extracted as aggregated data. The triaxial values of the Microelectromechanical System Inertial Measurement Units (MIMS) at the minute level, denoted by the variable name PAXMTSM, were employed in the computation of the rest-activity rhythm. Accelerometer count values pertaining to a patient were designated as missing when PAXMTSM was recorded as “−0.01.” Additionally, if the device failed to continuously capture triaxial values, the accelerometer count values were also treated as missing. The procedures for accelerometer recording and data preprocessing have been outlined in previous studies (11, 12).

The R package “nparACT” was utilized to calculate the resultant nonparametric variables that characterize rest-activity rhythms, a methodology extensively elucidated in previous literature (13): (a) interdaily stability (IS), which gauges the fidelity of the rest-activity pattern in adherence to the light–dark cycle (IS ≃ 0 denotes resemblance to Gaussian noise, while IS ≃ 1 signifies optimal stability); (b) intradaily variability (IV), quantifying the extent of fragmentation (IV ≃ 0 represents a flawless sine wave, and IV ≃ 2 reflects Gaussian noise); (c) relative amplitude (RA), delineating the disparity between the peak activity period encompassing 10 consecutive hours (M10) and the trough activity period spanning 5 consecutive hours (L5) over an average 24-h period. (d) M10 start time denoting the initiation time of the peak activity phase; and (e) L5 start time offering insights into the commencement time of the nadir activity phase.

The definition of epilepsy

The epilepsy-related information was derived from questionnaire dataset known as “prescription medications.” Participants self-reported their use of medications to treat “epilepsy and recurrent seizures” (International Classification of Diseases: G40). Furthermore, we manually screened the medications used by these participants, and those who were not taking standard antiseizure medications (ASM) were excluded from the analysis. Details are presented in Datasheet 1.

Confounding variables

All confounding variables were discretized into categorical data. The demographic information was based on self-reported data, including age (3–19, 20–39, 40–59, > = 60 years), gender (female, male), race/ethnicity (Mexican, other Hispanic, non-Hispanic Black, non-Hispanic White, Other/multiracial), annual family income (<= $44,999, > = $45,000), and marriage status (never, married, separated). Body mass index (BMI) data was from examination data and classified into 0–18.5, 18.5–25, 25–30, >30 groups.

Statistical analysis

We computed the unweighted frequencies and weighted proportions of categorical variables for both the epilepsy and non-epilepsy cohorts. Continuous variables were expressed as weighted means along with their corresponding 95% confidence intervals (CI). Subsequently, we employed multivariable logistic regression analysis using survey weights to ascertain the odds ratio (OR) of epilepsy with respect to rest-activity rhythm parameters. The strength of the association was evaluated using regression coefficients or ORs, accompanied by their respective 95% confidence intervals. We assessed three models with different sets of covariates. Model 1 did not involve any adjustments for additional covariates. Model 2 incorporated age, gender, and race as covariates. Model 3 further adjusted for annual family income and BMI. Statistical significance was defined as a p-value below 0.05. All statistical analyses were performed using the R software package.

Result

Baseline characteristics

All means (SD) and percentages were weighted based on the NHANES sampling frame. A total of 7,410 participants were included in the combined 2013–2014 NHANES sample. The mean age was 39.62 years (SD: 40.48), 52% were male, and 63% were non-Hispanic white. Among them, 53 individuals (0.75%) were identified as epilepsy patients. Detailed information can be found in Table 1.

Table 1
www.frontiersin.org

Table 1. Characteristics of participants.

Different rest-activity circadian rhythm

Differences in the circadian rest-activity rhythm indices between epilepsy patients and healthy controls are presented in Table 2, with weighted means and standard deviations. Compared to healthy participants, epilepsy patients exhibited a lower IS value (p < 0.05), indicating lower daytime stability in epilepsy patients compared to healthy individuals. Additionally, the M10 of epilepsy patients was significantly lower than that of normal individuals.

Table 2
www.frontiersin.org

Table 2. Differences in circadian rhythms between individuals with epilepsy and the normal population.

Relationship between rest-activity circadian rhythm and epilepsy

We classified the circadian rhythm parameters based on quartiles and used logistic regression analysis to examine the relationship between different circadian rhythm parameters and epilepsy. In all models, individuals categorized in the highest quartile of IS, RA, and M10 exhibited a reduced prevalence of epilepsy in contrast to those classified in the lowest quartile. Additionally, following adjustments for all covariates, individuals falling into the highest quartile of IV and M10 start time demonstrated an elevated prevalence of epilepsy in comparison to those in the lowest quartile (IV: OR = 2.52, 95% CI: 1.04–7.09; M10 start time: OR = 3.13, 95% CI: 1.39–7.58) (see Table 3).

Table 3
www.frontiersin.org

Table 3. Relationship between rest-activity circadian rhythm and epilepsy.

We further employed a restricted cubic spline (RCS) model to capture the nonlinear relationships between IS, RA, M10 and epilepsy. We observed that these relationships remained statistically significant (Figure 2). Although the outcomes of the nonlinear examination of the RCS exhibited minor discrepancies, the overarching patterns of the reliant and autonomous variables remained largely coherent throughout the diagrams.

Figure 2
www.frontiersin.org

Figure 2. RCS models for epilepsy and rest-activity circadian rhythm. The RCS assessment evaluates the correlation of IS, RA, and M10 with epilepsy. The red solid line corresponds to the central estimate, while the purple shaded area represents the 95% confidence interval.

Discussion

Despite advances in EEG technology deepening our understanding of the association between circadian rhythms and the risk of epileptic seizures, there is a clear lack of clinical studies based on objective data. Previously, only one study utilized actigraphy to assess the relationship between rest-activity cycles and epilepsy (9). In contrast to this study, our cross-sectional study has a larger sample size, enabling more robust statistical analyses. Additionally, the NHANES database represents a wide spectrum of the U.S. population, while the study by Claudio Liguori et al. was based solely on patients visiting a healthcare center in Rome, potentially introducing selection bias.

Participants in NHANES aged 3 years and older were instructed to wear ActiGraph GT3X+ devices continuously for 7 days to collect information on 24-h movement during wakefulness and sleep. Actigraphy is considered the optimal tool for studying circadian rhythms and sleep (14). The current study represents the first comparison of circadian rhythms in individuals with epilepsy and the general population, providing the most comprehensive objective analysis of epilepsy and circadian rhythms to date. Through rigorous preprocessing and data transformation, we were able to directly identify differences in seven circadian rhythm parameters between individuals with epilepsy and the normal population. The significant differences identified were in IS and M10: individuals with epilepsy exhibited a more synchronized rest-activity rhythm with the light–dark cycle compared to the normal population, with increased daytime activity levels.

After multivariable logistic regression analysis, we found that Participants with greater IS were linked with a decreased prevalence of epilepsy, whereas individuals with elevated levels of IV displayed a converse relationship. These are all evidence of a close association between disruptions in the 24-h rest-activity circadian rhythm and epilepsy. Although the directionality of the associations in this cross-sectional study cannot be determined, existing evidence suggests a complex bidirectional relationship between circadian rhythms and epilepsy (15). Circadian rhythms are an inherent part of the biological physiological 24-h daily rhythm, influencing seizure occurrence and activity. This association may be related to the regulation of membrane excitability and neurotransmitter imbalance in neuronal circuits by circadian rhythms (16), as well as the involvement of core clock genes and circadian regulatory factors in influencing the expression of genes associated with epilepsy (17). During sleep, synchronous discharge of thalamocortical networks promotes the generation of non-REM sleep oscillations (18), while the circadian release of cortisol and melatonin may influence seizure occurrence (19). Studies indicate that the circadian rhythm system influences epilepsy susceptibility and seizure generation by regulating the timing of wakefulness and sleep stages and brain function (20). Sleep and wakefulness are considered the most reliable predictors of seizure occurrence, especially in temporal lobe epilepsy, where seizures are more frequent during sleep or exclusively occur during sleep. Furthermore, with increasing time awake, cortical excitability also increases, with a greater reduction in nighttime cortical excitability compared to daytime, potentially increasing the risk of seizures in the morning after sleep deprivation. Animal models also support the regulatory role of circadian rhythms in seizure occurrence, such as rodent models demonstrating endogenous circadian distribution of seizures (21, 22). These findings emphasize the significant impact of sleep–wake behaviors and circadian rhythms on epilepsy susceptibility and seizures, particularly in shaping the temporal patterns of seizures by regulating the timing of wakefulness and sleep stages and brain function. Therefore, understanding the circadian patterns of seizures and trends influenced by various factors, including epilepsy type, seizure onset region, and age, is crucial for predicting seizures, improving treatment, and ensuring patient safety.

Furthermore, our analysis revealed that participants in the highest quartile of M10, a variable indicative of the level of physical activity, exhibited a 76% reduced prevalence of epilepsy in contrast to those categorized in the lowest quartile. We believe this may be related to discouragement of physical exercise in individuals with epilepsy (23). In the past, exercise was thought to potentially trigger seizures by increasing stress, hyperventilation, liver enzyme metabolism, and exercise-related injuries (24). However, in reality, physical activity has no adverse effects for most patients, and animal models and cohort studies have confirmed that physical activity can improve epilepsy and its comorbidities (25, 26). Therefore, healthcare professionals should provide health education and encourage individuals with epilepsy to engage in safe physical activities.

Although this study has yielded exciting results, there are several limitations. Firstly, the type of epilepsy (generalized or focal) and the site of seizure onset (frontal or temporal lobe) play a role in the circadian epilepsy pattern. We identified individuals with epilepsy as those self-reporting seizures treated with standard ASM medication for epilepsy, and this information is not sufficient to determine the classification of epilepsy. Secondly, we utilized actigraphy data to generate indices of circadian rhythms. We were unable to verify the quality of the transformed data. However, rigorous preprocessing was conducted to ensure the quality of the seven indices. Thirdly, our raw data did not indicate specific timing of epileptic seizures, which hindered further elucidation of seizure patterns. Fourthly, the study did not specify whether participants were seizure-free while taking antiepileptic medications (ASM), which is crucial for understanding the relationship between circadian rhythms and seizure occurrence. Therefore, using ‘the intake of at least one antiepileptic medication (ASM)’ or the ICD code (G40) as proxy variables for identifying epilepsy may inherently introduce a potential risk of bias and limit the ability to draw conclusions about seizure control in this population. We advocate for more research based on objective data approaching from the perspective of circadian rhythms to offer new insights into predicting epilepsy.

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.

Author contributions

TT: Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft. YZ: Conceptualization, Methodology, Supervision, Writing – review & editing. XZ: Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by Chongqing Technology Innovation and Application Development Special Project (CSTB2024TIAD-LUX0001).

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.

Abbreviations

NHANES, National Health and Nutrition Examination Survey; IS, Interdaily Stability; IV, Intradaily Variability; RA, Relative Amplitude; ASM, Antiseizure medications; RCS, Restricted cubic spline.

References

1. Berg, AT, Berkovic, SF, Brodie, MJ, Buchhalter, J, Cross, JH, van Emde Boas, W, et al. Revised terminology and concepts for organization of seizures and epilepsies: report of the ILAE commission on classification and terminology, 2005-2009. Epilepsies. (2010) 51:676–85. doi: 10.1111/j.1528-1167.2010.02522.x

PubMed Abstract | Crossref Full Text | Google Scholar

2. Thijs, RD, Surges, R, O'Brien, TJ, and Sander, JW. Epilepsy in adults. Lancet. (2019) 393:689–701. doi: 10.1016/S0140-6736(18)32596-0

PubMed Abstract | Crossref Full Text | Google Scholar

3. Hofstra, WA, and de Weerd, AW. The circadian rhythm and its interaction with human epilepsy: a review of literature. Sleep Med Rev. (2009) 13:413–20. doi: 10.1016/j.smrv.2009.01.002

PubMed Abstract | Crossref Full Text | Google Scholar

4. Mulè, A, Bruno, E, Pasanisi, P, Galasso, L, Castelli, L, Caumo, A, et al. Sex differences in rest-activity circadian rhythm in patients with metabolic syndrome. Front Physiol. (2021) 12:641461. doi: 10.3389/fphys.2021.641461

PubMed Abstract | Crossref Full Text | Google Scholar

5. Spencer, DC, Sun, FT, Brown, SN, Jobst, BC, Fountain, NB, Wong, VS, et al. Circadian and ultradian patterns of epileptiform discharges differ by seizure-onset location during long-term ambulatory intracranial monitoring. Epilepsia. (2016) 57:1495–502. doi: 10.1111/epi.13455

PubMed Abstract | Crossref Full Text | Google Scholar

6. Anderson, CT, Tcheng, T, Sun, FT, Fau, SF, Morrell, MJ, and Morrell, MJ. Day-night patterns of Epileptiform activity in 65 patients with long-term ambulatory Electrocorticography. J Clin Neurophysiol. (2015) 32:406–12. doi: 10.1097/WNP.0000000000000183

PubMed Abstract | Crossref Full Text | Google Scholar

7. Bijlenga, D, Fronczek, R, Gorter, EJ, and Thijs, RD. Vigilance and circadian function in daytime and nocturnal epilepsy compared to controls. Epilepsy Res. (2023) 197:107238. doi: 10.1016/j.eplepsyres.2023.107238

PubMed Abstract | Crossref Full Text | Google Scholar

8. Unterberger, I, Gabelia, D, Prieschl, M, Chea, K, Hofer, M, Högl, B, et al. Sleep disorders and circadian rhythm in epilepsy revisited: a prospective controlled study. Sleep Med. (2015) 16:237–42. doi: 10.1016/j.sleep.2014.09.021

PubMed Abstract | Crossref Full Text | Google Scholar

9. Liguori, C, Spanetta, M, Fernandes, M, Izzi, F, Placidi, F, and Mercuri, NB. More than sleep and wake disturbances: an actigraphic study showing the sleep-wake pattern dysregulation in epilepsy. Seizure. (2022) 94:95–9. doi: 10.1016/j.seizure.2021.11.024

Crossref Full Text | Google Scholar

10. Chiappa, M, and Porter, KS. National health and nutrition examination survey: plan and operations, 1999-2010. Vital Health Stat. (2013) 1:1–37.

PubMed Abstract | Google Scholar

11. Xu, Y, Su, S, McCall, WV, Isales, C, Snieder, H, and Wang, XA-O. Rest-activity circadian rhythm and impaired glucose tolerance in adults: an analysis of NHANES 2011-2014. BMJ Open Diabetes Res Care. (2022) 10:e002632. doi: 10.1136/bmjdrc-2021-002632

PubMed Abstract | Crossref Full Text | Google Scholar

12. Cespedes Feliciano, EM, Quante, M, Weng, J, Mitchell, JA, James, P, Marinac, CR, et al. Actigraphy-derived daily rest-activity patterns and body mass index in community-dwelling adults. Sleep. (2017) 40:zsx168. doi: 10.1093/sleep/zsx168

PubMed Abstract | Crossref Full Text | Google Scholar

13. Gonçalves, BS, Adamowicz, T, Louzada, FM, Moreno, CR, and Araujo, JF. A fresh look at the use of nonparametric analysis in actimetry. Sleep Med Rev. (2015) 20:84–91. doi: 10.1016/j.smrv.2014.06.002

Crossref Full Text | Google Scholar

14. Ancoli-Israel, S, Cole, R, Fau, AC, Chambers, M, Fau, CM, Moorcroft, W, et al. The role of actigraphy in the study of sleep and circadian rhythms. Sleep. (2003) 26:342–92. doi: 10.1093/sleep/26.3.342

PubMed Abstract | Crossref Full Text | Google Scholar

15. Khan, S, Nobili, L, Khatami, R, Loddenkemper, T, Cajochen, C, Dijk, DJ, et al. Circadian rhythm and epilepsy. Lancet Neurol. (2018) 17:1098–108. doi: 10.1016/S1474-4422(18)30335-1

Crossref Full Text | Google Scholar

16. Zaitsev, AA-O . Molecular and cellular mechanisms of epilepsy 2.0. Int J Mol Sci. (2023) 24:17464. doi: 10.3390/ijms242417464

Crossref Full Text | Google Scholar

17. Partch, CL, Green, CB, and Takahashi, JS. Molecular architecture of the mammalian circadian clock. Trends Cell Biol. (2014) 24:90–9. doi: 10.1016/j.tcb.2013.07.002

Crossref Full Text | Google Scholar

18. Saper, CB, Te Fau, S, Lu, J, and Lu, J. Hypothalamic regulation of sleep and circadian rhythms. Nature. (2005) 437:1257–63. doi: 10.1038/nature04284

Crossref Full Text | Google Scholar

19. Gamble, KL, Berry, R, Frank, SJ, and Young, ME. Circadian clock control of endocrine factors. Nat Rev Endocrinol. (2014) 10:466–75. doi: 10.1038/nrendo.2014.78

Crossref Full Text | Google Scholar

20. Stirling, RA-O, Hidajat, CM, Grayden, DA-O, D'Souza, WJ, Naim-Feil, J, Dell, KL, et al. Sleep and seizure risk in epilepsy: bed and wake times are more important than sleep duration. Brain. (2023) 146:2803–13. doi: 10.1093/brain/awac476

Crossref Full Text | Google Scholar

21. Pitsch, J, Becker, AJ, Schoch, S, Müller, JA, de Curtis, M, and Gnatkovsky, V. Circadian clustering of spontaneous epileptic seizures emerges after pilocarpine-induced status epilepticus. Epilepsia. (2017) 58:1159–71. doi: 10.1111/epi.13795

Crossref Full Text | Google Scholar

22. Quigg, M, Straume, M, Menaker, M, and Bertram, EH. Temporal distribution of partial seizures: comparison of an animal model with human partial epilepsy. Ann Neurol. (1998) 43:748–55. doi: 10.1002/ana.410430609

Crossref Full Text | Google Scholar

23. Arida, RM, Fau, CE, da Silva, AC, Fau, d SA, Scorza, FA, and Scorza, FA. Physical activity and epilepsy: proven and predicted benefits. Sports Med. (2008) 38:607–15. doi: 10.2165/00007256-200838070-00006

Crossref Full Text | Google Scholar

24. Pimentel, J, Tojal, R, and Morgado, J. Epilepsy and physical exercise. Seizure. (2015) 25:87–94. doi: 10.1016/j.seizure.2014.09.015

Crossref Full Text | Google Scholar

25. Kanner, AM . The treatment of depressive disorders in epilepsy: what all neurologists should know. Epilepsia. (2013) 1:3–12. doi: 10.1111/epi.12100

Crossref Full Text | Google Scholar

26. Tutkun, E, Ayyildiz, M, Fau Agar, E, and Agar, E. Short-duration swimming exercise decreases penicillin-induced epileptiform ECoG activity in rats. Acta Neurobiol Exp (Wars). (2010) 70:382–9. doi: 10.55782/ane-2010-1810

Crossref Full Text | Google Scholar

Keywords: epilepsy, circadian rhythm, NHANES, cross-sectional study, actigraphy data

Citation: Tang T, Zhou Y and Zhai X (2024) Circadian rhythm and epilepsy: a nationally representative cross-sectional study based on actigraphy data. Front. Neurol. 15:1496507. doi: 10.3389/fneur.2024.1496507

Received: 14 September 2024; Accepted: 31 October 2024;
Published: 03 December 2024.

Edited by:

Janelle Wagner, Medical University of South Carolina, United States

Reviewed by:

Arturo Garay, Norberto Quirno Medical Education and Clinical Research Center (CEMIC), Argentina
Elisa Evangelista, Centre Hospitalier Universitaire de Nîmes, France

Copyright © 2024 Tang, Zhou and Zhai. 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: Xuan Zhai, emhhaXh1YW5AY3FtdS5lZHUuY24=; YuDong Zhou, MzgyNzY2MjcxQHFxLmNvbQ==

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.