ORIGINAL RESEARCH article
The Need for Structured Strategies to Improve Stroke Care in a Rural Telestroke Network in Northern New South Wales, Australia: An Observational Study
- 1School of Medicine and Public Health, The University of Newcastle, Newcastle, NSW, Australia
- 2Department of Neurology, John Hunter Hospital, Newcastle, NSW, Australia
- 3Department of Neuroscience, Gosford Hospital, Gosford, NSW, Australia
- 4Department of Neurology, Melbourne Brain Center at Royal Melbourne Hospital, Parkville, VIC, Australia
- 5Department of Neurology, Liverpool hospital, Liverpool, NSW, Australia
- 6Department of Neurology, The University of Utah, Salt Lake City, UT, United States
- 7Hunter New England Health, Wallsend, NSW, Australia
Introduction: A telestroke network in Northern New South Wales, Australia has been developed since 2017. We theorized that the telestroke network development would drive a progressive improvement in stroke care metrics over time.
Aim: This study aimed to describe changes in acute stroke workflow metrics over time to determine whether they improved with network experience.
Methods: We prospectively collected data of patients assessed by telestroke who received multimodal computed tomography (mCT) and were diagnosed with ischemic stroke or transient ischemic attack from January 2017 to July 2019. The period was divided into two phases (phase 1: January 2017 – October 2018 and phase 2: November 2018 – July 2019). We compared median door-to-call, door-to-image, and door-to-decision time between the two phases.
Results: We included 433 patients (243 in phase 1 and 190 in phase 2). Each spoke site treated 1.5–5.2 patients per month. There were Door-to-call time (median 39 in phase 1, 35 min in phase 2, p = 0.18), and door-to-decision time (median 81.5 vs. 83 min, p = 0.31) were not improved significantly. Similarly, in the reperfusion therapy subgroup, door-to-call time (median 29 vs. 24.5 min, p = 0.12) and door-to-decision time (median 70.5 vs. 67.5 min, p = 0.75) remained substantially unchanged. Regression analysis showed no association between time in the network and door-to-decision time (coefficient 1.5, p = 0.32).
Conclusion: In our telestroke network, acute stroke timing metrics did not improve over time. There is the need for targeted education and training focusing on both stroke reperfusion competencies and the technical aspects of telestroke in areas with limited workforce and high turnover.
The increased availability of telestroke in Australia means that physicians from rural areas can consult expert stroke neurologists at comprehensive stroke centers (1). This support leads to an increasing number of patients receiving reperfusion therapies (2–5), being particularly relevant after the recent trials expanding the window for endovascular thrombectomy (EVT) (6, 7) and thrombolysis (8).
Since April 2013, a telestroke network using multimodal computed tomography (mCT - including brain non-contrast CT, CT angiography, and CT perfusion) has been developed in Northern New South Wales, Australia. This network had a progressive expansion, accelerated since 2017, increasing the number of patients assessed via telestroke from 48 in 2016 to over 600 in 2019 (9, 10).
In August 2019, as part of a support system for regional and rural stroke care, the clinical trial TACTICS (Trial of Advanced CT Imaging and Combined Education Support for Drip and Ship) was launched. TACTICS aimed to boost access to reperfusion therapies using a multifaceted education and training program targeting opportunities for workflow improvement in regional and rural acute stroke care (ACTRN12619000750189). In this study, we aim to describe the acute stroke treatment timing metrics in our telestroke network following full implementation of telestroke care across the network in 2017 but prior to commencement of the TACTICS trial. We hypothesized there would be trends to improvement in workflows with growth in network experience over time (11, 12).
Materials and Methods
The network commenced with a single spoke hospital in April 2013, with a second site added in 2014 (13). In 2017, the telestroke network incorporated three new spoke hospitals, finishing its expansion in 2018 with the inclusion of the last hospital. The rural population that the telestroke network covers has increased from 110,000 in 2014 to 300,000 people (13) in 2018, distributed over an area of 143,120 km (slightly larger than England) (2). The John Hunter Hospital (Newcastle) is the comprehensive stroke center (the hub hospital), and sole hospital providing EVT to the area. In our network, all stroke patients except those who require EVT are treated in spoke hospitals.
Regular visits to the spoke sites were performed by the hub team (average of 1–2 per year), and teaching sessions focused on acute stroke pathway were organized (2–3 per year). Immediate feedback after complex cases was regularly given via local stroke coordinators.
All six spoke hospitals were equipped with cameras in the Emergency Department and the local radiographers were trained in multimodal CT acquisition. We refer the readers to previous publications, in which our telestroke workflow has been described in detail (9, 10).
Study Design and Data Collection
Patients assessed by telestroke between January 2017 and July 2019 were included.
Patients from the 2013–2016 period were not included in the analysis because the network expansion and infrastructure deployment were still in progress. During this period just 112 patients were assessed.
We divided the analysis in two phases: Phase 1, from January 2017 to October 2018 and phase 2, from November 2018 to July 2019. Phase 1 was selected to begin at the formalization of the central telestroke network and end once the majority of sites had participated in the network for at least 12 months. We excluded one site that was not initiated telestroke until September 2018 and this site contributed fewer than 20 consults per year. Clinical data was prospectively collected including baseline demographics, National Institute of Health Stroke Scale (NIHSS), mCT characteristics (vessel occlusion status, ischemic core and penumbra volumes), and modified Rankin Scale (mRS) before, and 3 months after, stroke. Large vessel occlusion (LVO) was defined as an occlusion in the proximal segment of the middle cerebral artery (M1 or M2), terminal internal carotid artery, basilar artery or combined intra- and extracranial occlusions. Date and time of stroke onset (or last time seen well), telestroke consult initiation, treatment decision time, and thrombolysis time were also collected. Onset-to-door time was calculated as the difference between stroke onset/last time seen well and Emergency Department (ED) arrival. Door-to-call time was calculated as the difference between patient arrival at the rural hospital and the call to the telestroke neurologist. Door-to-image time was calculated as the difference between patient arrival and mCT initiation. Door-to-decision time was calculated as difference between patient arrival and time when the neurologist notified the spoke team of the plan for the patient (candidate/not candidate for reperfusion therapies). Door-to-needle time was calculated as the difference between patient arrival time and thrombolysis bolus initiation. Patients who suffered a stroke in hospital were excluded from door-to-call, door-to-imaging and door-to-decision analyses. Patients were diagnosed with a stroke mimic by telestroke neurologists if the clinical presentation was atypical for stroke and follow up brain CT/MRI was normal.
The mCT imaging protocol included brain non-contrast CT, CT angiography, and CT perfusion at baseline and either brain CT or MRI follow-up as per routine clinical practice. All perfusion images were analyzed using MIStar (Apollo Medical Imaging Technology, Melbourne, Australia). Tissue with a delay time >3 s and a relative cerebral blood flow >30% of the contralateral hemisphere was defined as penumbra, and the tissue with a delay time >3 s and a relative cerebral blood flow <30% was defined as ischemic core (14).
Groups were compared with the independent samples Mann–Whitney U-test. For categorical variables, Pearson's Chi squared test or Fisher's exact test were used as appropriate. Statistical significance was set at 0.05. Time metrics were also assessed using a regression model. To examine the effect of telestroke experience, we set treatment timings as a dependent variable and set the time elapsed since January 2017, NIHSS, ED arrival within 4.5 h from the onset, and reperfusion therapy as independent variables, since our previous research showed symptom severity, fast ED arrival and reperfusion therapy related to faster treatment timings (15). Partial missing data from the database are random and the proportion of cases missing data was relatively small, thus we used the case deletion method. All statistical analysis was performed on Stata version 14 (Statacorp, USA). Results are presented as mean ± standard deviation (SD) or median [Interquartile range-IQR].
From January 2017 to July 2019, 827 patients were assessed by the telestroke network, 434 in phase 1 and 393 in phase 2. mCT was performed in 591 patients (71.5%) (323 and 268 patients in phase 1 and 2, respectively). Of those patients receiving mCT, 433 (73.3%) were confirmed to have ischemic stroke or transient ischemic attack (TIA) and 145 (24.5%) were stroke mimics. 236 patients (28.5%) did not undergo urgent mCT for the following reasons: minor symptoms at presentation (73 cases), initial presentation unlikely to be a stroke based on neuro-assessment (59 cases) or presence of hemorrhage in brain NCCT (38 cases) (Figure 1).
Figure 1. Diagnoses in each phase and the reasons why mCT was not used. mCT, multimodal Computed Tomography (CT); TIA, Transient Ischemic Attack; ICH, Intracranial Hemorrhage; NCCT, non-contrast CT; mRS, modified Rankin Scale.
The median age of the 433 patients with confirmed ischemic stroke/TIA was 74 years [IQR 64–81] and 277 were male (64%). Most of the patients were independent before their stroke (mRS 0–2, 391 patients, 90.3%). The median baseline NIHSS was 4 [2–9], and 121 patients (28.0%) had a large vessel occlusion. The median penumbra volume was 5 [0–41] ml and the median core volume was 1 [0–7] ml. A total of 109 patients (25.2%) received reperfusion therapies; 74 patients (68.0%) underwent thrombolysis, 58 received EVT (53.2%), and 23 (21.1%) patients underwent combined therapy.
As a whole, the median door-to-call time was 36 [IQR 21–63] min, the median door-to-image was 50 [33–77] min, and the median door-to-decision was 82 [64–118] min. There was no significant difference between phases for these metrics. Phase 1 had a median door-to-call of 39 min compared to 35 min in phase 2 (p = 0.18), median door-to-image was 49 min in phase 1 and 54 min in phase 2 (p = 0.36) and median door-to-decision was 81.5 min in phase 1 and 83 min in phase 2 (p = 0.31) (Table 1, Figure 2A). The results did not change even when experience was analyzed based on the date of each individual site's telestroke initiation (data not shown).
Table 1. Clinico-radiological characteristics and workflow metrics of ischemic stroke and transient ischemic attack patients assessed with mCT.
Figure 2. Superimposed bar graphs of Door-to-Call, Door-to-Image, Door-to-Decision, and Door-to-Needle time. Numbers in each bar are median time and [IQR] (min). (A) Confirmed ischemic stroke patients and TIA patients in phase 1 and 2. (B) Patients who underwent reperfusion therapy. (C) Patients who underwent thrombolysis.
Regression analysis showed no significant association between the time from telestroke network initiation and any stroke metrics (door-to-call, door-to-image, and door-to decision times) except call-to-decision time, which showed 2.0 min increase in delay per 3 months (95% CI 0.6–3.4, p < 0.001) (Table 2).
Table 2. Estimated effects on treatment timings of the experience, NIHSS, ED arrival within 4.5 h and the use of reperfusion therapy.
Reperfusion Therapy Subgroup
Of the 109 patients who received reperfusion therapy, the median age was 73 [64–80] years and 66 were male (60.6%), and most of them (98.2%) were independent at baseline. The baseline NIHSS was 9 [6–15], and 84 (77.1%) had a large vessel occlusion. The baseline characteristics between phase 1 and 2 for these patients were similar, but there was a higher proportion of male sex in phase 1 (69.1% in phase 1 vs. 46.3% in phase 2, p = 0.02) and a higher proportion of large vessel occlusion strokes in phase 2 (70.6 vs. 87.8%, p = 0.04). There was a significant decrease in the rate of thrombolysis in phase 2 compared to phase 1 (76.5 vs. 53.7%, p = 0.01) but an increase in use of EVT (45.6 vs. 65.9%, p = 0.04) (Table 3).
Table 3. Clinico-radiological characteristics and workflow metrics of patients who underwent reperfusion therapy.
The median door-to-call time was 26.5 min [17–43], median door-to-image was 35 min [21–52], and median door-to-decision was 70 min [55–95]. All these metrics were slightly shorter in phase 2, but none of the differences was significant. Door-to-call was 29 min in phase 1 vs. 24.5 in phase 2 (p = 0.12), door-to-image 37.5 min vs. 34.5 min (p = 0.87) and door-to-decision was 70.5 min in phase 1 and 67.5 min in phase 2 (p = 0.82); see Table 3 and Figure 2B.
In the thrombolysed subgroup, door-to-needle time were slightly shorter in phase 2 (93.5 min in phase 1 vs. 89 min in phase 2, p = 0.5) but this difference was not significant. No difference was found in any of the other metrics (Figure 2C).
Regression analysis demonstrated no significant association between time since the telestroke network initiation and door-to-needle (Adjusted β = 1.1 min, p = 0.62) or decision-to-needle times (Adjusted β =1.2 min, p = 0.15).
In this study, we described the temporal trends in acute stroke workflow metrics in our telestroke network. We noted no significant improvements in process of care metrics over a period of 2.5 years. Although, some publications suggest that participation in telestroke networks improves treatment timings in acute stroke (11, 12), this was not reflected in the data from our network. Door-to-call time slightly improved in phase 2 compared to phase 1 (from 39 to 35 min), however, door-to-image and door-to-decision times were slightly extended.
There are several possible explanations for this slight increase in time metrics in our telestroke network. Firstly, there is a limited number of acute stroke patients that each spoke hospitals' physicians could see (1.5–5.2 stroke patients per month in the study period). In addition, the rural hospitals suffer from workforce shortages and a high turnover across all career stages, from specialists to trainees (16). Therefore, clinical experience does not concentrate within a specific team but disseminates. Staff experience is an essential factor to improve acute stroke quality, and there is a report that showed center volume had a strong effect on shortened Door-to-Needle time (17). The hub stroke team has provided regular visits and feedback through stroke coordinators at the spoke sites, nevertheless, our result showed that this support was not enough to improve acute stroke timing metrics for spoke hospitals in our telestroke network. Therefore, there is a need for structured strategies including stroke-specific education and training, and a system that encourages and motivates health professionals in spoke hospitals. Rural settings often face well-recognized barriers to system change, including facility leadership support, competing demands, stakeholder resistance, resource, and technical support constraints (18–21). We speculate that a lack of improvement in acute stroke time metrics might be ascribed to these disadvantages.
Previous studies have indicated that telestroke networks can successfully improve their performance level with appropriate initiatives to improve the quality of acute stroke care (22, 23). Aiming to improve acute stroke workflows and times to reperfusion, a trial across the network with a focus in multifaceted education and support for workflow re-design (the TACTICS trial), is now underway. The intervention includes teleconferences within the network and reviews of current acute stroke workflow to identify barriers and provide prioritized solutions and strategies as well as stroke-specific education to health professionals (ACTRN12619000750189). The comprehensive stroke center team collaborate with the network hospitals by simulating action planning and ongoing clinical practice improvement activities.
Besides reperfusion therapy, all patients who suffer stroke can benefit from specialized stroke-unit care, and telestroke networks can play a role in this (24). Each spoke site has a stroke unit that provides care to all stroke patients admitted to that site, but the function of the telestroke network is currently limited to the provision of acute reperfusion therapy. Establishing a system that supports and educates health professionals in spoke hospitals regarding general acute stroke care through telestroke remains a long-term goal.
Time metrics in our telestroke network are comparatively long, considering the guideline recommendation (25). Workforce shortages (emergency teams must treat multiple patients simultaneously), inefficient workflows (absence of pre-notification from paramedics, telestroke activation call after initial imaging), and the infrequent use of acute reperfusion therapies by health professionals in spoke hospitals all contribute to the slow time metrics (Supplementary Table 1). In addition, our telestroke network uses mCT as the standard image for acute stroke care. We acknowledge that mCT adds 10–13 min to the door-to-decision/needle times in the acute phase, but CT perfusion is required to select candidates for EVT in the extended window and provides invaluable information regarding tissue prognosis in the standard window. In a geographically disparate network like ours, where transfers can take over 5 h and transport assets are extremely limited, careful triage to avoid unnecessary transfers is essential (9, 10).
There are some limitations in this study, most importantly the small number of patients treated. The requirement for patients to undergo mCT to be included in this analysis introduces a potential bias, as stroke patients who did not undergo mCT due to technical failures were excluded. However, the occurrence rate of the technical failures was low (3% of the whole sample in this study) and random in nature, so those patients are unlikely to affect the result in either a positive or negative direction. Lastly, the generalizability of the study is limited, and it may only be applicable to other rural telestroke networks with similar structure and limited workforce. Conversely, strengths include the fact that all the data were collected prospectively, and various analyses of the timing data consistently show no significant difference between groups, including the regression analysis which did not require any arbitrary time division. While this study may not be generalizable to all telestroke networks, the problems of limited experience in spoke hospitals, high staff turnover and work-force shortages are common problems in rural areas both in Australia and overseas (26, 27).
In conclusion, our study shows that in our telestroke network, acute stroke timing metrics did not improve over time. There is the need for targeted education and training focusing on both stroke reperfusion competencies, and the technical aspects of telestroke in areas with limited workforce and high turnover.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The studies involving human participants were reviewed and approved by the Hunter New England Health Ethics Committee. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
CG-E, NS, and CL conceived the present study. CG-E, FM, PG-B, SG, BC, WO'B, JE, and KA acquired the data, interpreted imaging data, and reviewed the manuscript. TL analyzed and interpreted the data and edited the manuscript. YK and CG-E analyzed the data and drafted and edited the manuscript. AB, MP, JM, NS, and CL supervised the study and critically reviewed the manuscript. All authors read and approved the final manuscript.
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.
The authors would like to thank all contributing collaborators of the Northern NSW Telestroke investigators for this project (See Supplementary Material).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.645088/full#supplementary-material
2. Wechsler LR, Demaerschalk BM, Schwamm LH, Adeoye OM, Audebert HJ, Fanale CV, et al. Telemedicine quality and outcomes in stroke: a scientific statement for healthcare professionals from the American Heart Association/American Stroke Association. Stroke. (2017) 48:e3–25. doi: 10.1161/STR.0000000000000114
3. Kepplinger J, Barlinn K, Deckert S, Scheibe M, Bodechtel U, Schmitt J. Safety and efficacy of thrombolysis in telestroke: a systematic review and meta-analysis. Neurology. (2016) 87:1344–51. doi: 10.1212/WNL.0000000000003148
4. Sørensen T, Dyb K, Rygh E, Salvesen R, Thomassen L. A qualitative description of telemedicine for acute stroke care in Norway: technology is not the issue. BMC Health Serv Res. (2014) 14:643. doi: 10.1186/s12913-014-0643-9
6. Nogueira RG, Jadhav AP, Haussen DC, Bonafe A, Budzik RF, Bhuva P, et al. Thrombectomy 6 to 24 hours after stroke with a mismatch between deficit and infarct. N Engl J Med. (2018) 378:11–21. doi: 10.1056/NEJMoa1706442
7. Albers GW, Marks MP, Kemp S, Christensen S, Tsai JP, Ortega-Gutierrez S, et al. Thrombectomy for stroke at 6 to 16 hours with selection by perfusion imaging. N Engl J Med. (2018) 378:708–18. doi: 10.1056/NEJMoa1713973
8. Ma H, Campbell BCV, Parsons MW, Churilov L, Levi CR, Hsu C, et al. Thrombolysis guided by perfusion imaging up to 9 hours after onset of stroke. N Engl J Med. (2019) 380:1795–803. doi: 10.1056/NEJMoa1813046
9. Demeestere J, Sewell C, Rudd J, Ang T, Jordan L, Wills J, et al. The establishment of a telestroke service using multimodal CT imaging decision assistance: “Turning on the fog lights.” J Clin Neurosci. (2017) 37:1–5. doi: 10.1016/j.jocn.2016.10.018
10. Garcia-Esperon C, Soderhjelm Dinkelspiel F, Miteff F, Gangadharan S, Wellings T, O'Brien B, et al. Implementation of multimodal computed tomography in a telestroke network: five-year experience. CNS Neurosci Ther. (2019) 26:367–73. doi: 10.1111/cns.13224
11. Sanders KA, Patel R, Kiely JM, Gwynn MW, Johnston LH. Improving telestroke treatment times in an expanding network of hospitals. J Stroke Cerebrovasc Dis. (2016) 25:288–91. doi: 10.1016/j.jstrokecerebrovasdis.2015.09.030
12. Moreno A, Schwamm LH, Siddiqui KA, Viswanathan A, Whitney C, Rost N, et al. Frequent hub-spoke contact is associated with improved spoke hospital performance: results from the Massachusetts General Hospital telestroke network. Telemed e-Health. (2018) 24:678–83. doi: 10.1089/tmj.2017.0252
13. Australian Bureau of Statistics. Census 2016. (2016). Available online at: https://www.abs.gov.au/websitedbs/D3310114.nsf/Home/Census?OpenDocument&ref=topBar (accessed June 16, 2020).
15. Lillicrap T, Pinheiro A, Miteff F, Garcia-Bermejo P, Gangadharan S, Wellings T, et al. No evidence of the “Weekend Effect” in the Northern New South Wales telestroke network. Front Neurol. (2020) 11:130. doi: 10.3389/fneur.2020.00130
16. Australian Institute of Health and Welfare. Medical Practitioners Workforce 2015. (2016). Available online at: https://www.aihw.gov.au/reports/workforce/medical-practitioners-workforce-2015 (accessed March 27, 2020).
17. Strbian D, Ahmed N, Wahlgren N, Lees KR, Toni D, Roffe C, et al. Trends in door-to-thrombolysis time in the safe implementation of stroke thrombolysis registry: effect of center volume and duration of registry membership. Stroke. (2015) 46:1275–80. doi: 10.1161/STROKEAHA.114.007170
18. Levi CR, Attia JA, D'Este C, Ryan AE, Henskens F, Kerr E, et al. Cluster-randomized trial of thrombolysis implementation support in Metropolitan and Regional Australian Stroke Centers: lessons for individual and systems behavior change. J Am Heart Assoc. (2020) 9:e012732. doi: 10.1161/JAHA.119.012732
19. Tappen RM, Wolf DG, Rahemi Z, Engstrom G, Rojido C, Shutes JM, et al. Barriers and facilitators to implementing a change initiative in long-term care using the INTERACT® quality improvement program. Health Care Manag. (2017) 36:219–30. doi: 10.1097/HCM.0000000000000168
20. Geerligs L, Rankin NM, Shepherd HL, Butow P. Hospital-based interventions: a systematic review of staff-reported barriers and facilitators to implementation processes. Implement Sci. (2018) 13:36. doi: 10.1186/s13012-018-0726-9
22. Nguyen-Huynh MN, Klingman JG, Avins AL, Rao VA, Eaton A, Bhopale S, et al. Novel telestroke program improves thrombolysis for acute stroke across 21 hospitals of an integrated healthcare system. Stroke. (2018) 49:133–9. doi: 10.1161/STROKEAHA.117.018413
23. Mong R, Tiah L, Wong M, Tan C. Improving telestroke treatment times through a quality improvement initiative in a Singapore emergency department. Singapore Med J. (2019) 60:69–74. doi: 10.11622/smedj.2018083
25. Jauch EC, Saver JL, Adams HP, Bruno A, Connors JJB, Demaerschalk BM, et al. Guidelines for the early management of patients with acute ischemic stroke: a guideline for healthcare professionals from the American Heart Association/American Stroke Association. Stroke. (2013) 44:870–947. doi: 10.1161/STR.0b013e318284056a
Keywords: telestroke, acute stroke care, multimodal computed tomography, door-to-needle time, thrombolysis, thrombectomy
Citation: Kashida Y, Garcia-Esperon C, Lillicrap T, Miteff F, Garcia-Bermejo P, Gangadharan S, Chew BLA, O'Brien W, Evans J, Alanati K, Bivard A, Parsons M, Majersik JJ, Spratt NJ, Levi C and the members of Northern NSW Telestroke investigators for this project (2021) The Need for Structured Strategies to Improve Stroke Care in a Rural Telestroke Network in Northern New South Wales, Australia: An Observational Study. Front. Neurol. 12:645088. doi: 10.3389/fneur.2021.645088
Received: 22 December 2020; Accepted: 11 March 2021;
Published: 09 April 2021.
Edited by:Tae-Hee Cho, Hospices Civils de Lyon, France
Reviewed by:Klaus Fassbender, Saarland University Hospital, Germany
Christoph Gumbinger, Heidelberg University, Germany
Copyright © 2021 Kashida, Garcia-Esperon, Lillicrap, Miteff, Garcia-Bermejo, Gangadharan, Chew, O'Brien, Evans, Alanati, Bivard, Parsons, Majersik, Spratt, Levi and the members of Northern NSW Telestroke investigators for this project. 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: Yumi Tomari Kashida, firstname.lastname@example.org
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship