Shyam Gangadharan1*
Thomas Lillicrap2
Ferdinand Miteff1,2Pablo Garcia-Bermejo1Thomas Wellings1Billy O'Brien3James Evans3Khaled Alanati1
Christopher Levi1,2
Mark W. Parsons2,4Andrew Bivard2,4
Carlos Garcia-Esperon1,2†
Neil J. Spratt1,2† for the Northern NSW Telestroke Investigators- 1Department of Neurology, John Hunter Hospital, University of Newcastle, Newcastle, NSW, Australia
- 2Hunter Medical Research Institute, University of Newcastle, Newcastle, NSW, Australia
- 3Department of Neurology, Gosford Hospital, Gosford, NSW, Australia
- 4Department of Neurology, Melbourne Brain Centre at the Royal Melbourne Hospital, University of Melbourne, Parkville, VIC, Australia
Background and Purpose: Telestroke aims to increase access to endovascular clot retrieval (ECR) for rural areas. There is limited information on transfer workflow for ECR in rural settings. We sought to describe the transfer metrics for ECR in a rural telestroke network with respect to decision making.
Methods: A retrospective cohort study was employed on consecutive patients transferred to the comprehensive stroke center (CSC) for ECR in a rural hub-and-spoke telestroke network between April 2013 and October 2019, by road or air. Key time-based metrics were analyzed.
Results: Sixty-two patients were included. Mean age was 66 years [standard deviation (SD), 14] and median National Institutes of Health Stroke Scale 13 [interquartile range (IQR), 8–18]. Median rural-hospital-door-to-CSC-door (D2D) was 308 min (IQR, 254–351), of which 68% was spent at rural hospitals [door-in-door-out (DIDO); 214 min; IQR, 171–247]. DIDO was longer for air transfers than road (P = 0.004), primarily because of a median 87 min greater decision-to-departure time (Decision-DO, P < 0.001). In multiple linear regression analysis, intubation but not thrombolysis was associated with significantly longer DIDO. The distance at which the extra speed of an aircraft made up for the delays involved in booking an aircraft was 299 km from the CSC.
Conclusions: DIDO is longer for air retrievals compared with road. Decision-DO represents the most important component of DIDO, being longer for air transfers. Systems for rapid transportation of rural ECR candidates need optimization for best patient outcomes, with decision support seen as a potential tool to achieve this.
Introduction
Endovascular clot retrieval (ECR) is the standard of care in the treatment of acute stroke patients with large vessel occlusion (LVO) (1). It has been shown to be major benefits in selected patients up to 24 h (2, 3), but earlier treatment leads to greater benefit (4–6). This presents logistical challenges in Australia, because patients from widely dispersed geographic regions are eligible for this treatment, but it is offered only in limited metropolitan centers (7, 8).
Telestroke is being increasingly used in Australia to help overcome geographical disparities in access to acute stroke care (9–11). Optimal implementation of ECR through telestroke requires efficient workflow from primary hospital to the comprehensive stroke center. Door-in-door-out (DIDO) time at the primary hospital has been thought to have the greatest impact on outcome for patients with LVO being transferred for ECR, among modifiable factors (12–14). A recent study in metropolitan Australia proposed that the target time for DIDO should be shortened to 45 min (15, 16).
Despite the fact that about 29% of Australia's population live in rural and remote areas, with people in very remote areas having a mortality rate almost 1.4 times as high as in major cities (17), there is a paucity of data on transfer workflow specific to rural Australia to guide further development of ECR in these areas. We aimed to describe the transfer metrics for ECR from rural hospitals to a regional comprehensive stroke center in a telestroke network in rural Australia, with respect to key points in clinical decision making.
Clinical decision making in stroke care is complex (18), with there being a correlation between decision delay in acute stroke and both pre-hospital and in-hospital delays (19–21). In addition to transfer workflow, we extended our study to look at the relationship between transfer metrics and clinical decision making so that the results might be more readily adapted to modify routine clinical practice and potentially identify areas for decision support. We sought to compare different transport modalities. Our main hypothesis was that DIDO for air transfers would be longer than for road transfers.
Methods
Setting
A “hub and spoke” telestroke network was developed to support delivery of reperfusion therapy (thrombolysis and ECR) as a 24/7 service to six rural (spoke) hospitals in the Hunter New England and Mid-North Coast local health districts of New South Wales, Australia. The comprehensive stroke center hub based at John Hunter Hospital (JHH) in Newcastle, Australia, covers a population of over 1.1 million people distributed over an area of 143,120 km2 (slightly larger than England). The average distance between the spoke sites and the hub is 227 km (range 167–423 km). All six spoke sites have advanced imaging (including brain CT angiography and CT perfusion) capability needed to select patients for reperfusion therapies, with all spoke sites capable of administering thrombolysis. The rural spoke hospitals were equipped with telehealth cameras, and the local physicians were trained in the face, arm, speech, time (FAST) scale. Only one of the spoke hospitals had neurologists or stroke physicians on staff who could guide local decisions about in-hours thrombolysis. All other decisions about reperfusion therapy in the acute phase were made by the telestroke neurologist at the hub (JHH) after reviewing the imaging and if required assessing the patient remotely through the telehealth cameras. Neurologists at Gosford District Hospital also participated in the telestroke roster, although only JHH accepted patients for ECR.
The decision to transfer for ECR was made by the telestroke neurologist on an individualized basis taking into account standard clinical, imaging, and patient factors. From November 2017, the time window for ECR was expanded to 24 h in image-selected patients following the release of pivotal ECR trial results (2, 3). When it was deemed necessary to give thrombolysis with ECR, thrombolysis was always initiated at the spoke hospital prior to transfer for ECR. Repeat imaging was not routinely performed at the hub because of the commitment to transfer patients across very large distances after the decision for transfer was made and the good baseline selection of patients using advanced imaging (22). Road transfers occurred by road ambulance directly from the spoke hospital to JHH. Air transfers by helicopter occurred from the rural airport servicing the respective spoke hospital to the helipad at JHH. Air transfers by fixed wing occurred from the rural airport servicing the respective spoke hospital to the regional airport based 26 km from JHH, after which the transfer to JHH took place by road ambulance. All connecting transfers to the rural airports from the respective spoke hospitals in order to facilitate air retrievals occurred by road ambulance. For more details of the telestroke setup, we direct the reader to our published experience (22, 23).
Study Design and Data Collection
We employed a retrospective cohort study design, where transfer metrics were collected on consecutive telestroke patients transferred to JHH for ECR between April 2013 and October 2019. Clinical data were collected prospectively from June 2016 and retrospectively prior to that. Patients were transferred by road or air (helicopter or fixed wing), depending on availability and hospital distance.
Data collected included baseline demographics, past medical history, National Institutes of Health Stroke Scale (NIHSS), advanced imaging characteristics, acute treatment decision, and mode of retrieval.
Outcome Measures and Analyses
Time metrics included DIDO, defined as the duration of time from arrival to departure at the rural hospital. DIDO was divided into two segments, door-in-to-decision (DI-Decision) and decision-to-door-out (Decision-DO). Decision was defined as the time that the telestroke neurologist at JHH contacted the rural hospital regarding the decision to transfer. Door-to-door (D2D) was the time from arrival at the rural hospital to arrival at the comprehensive stroke center. We also examined the total time from decision to arrival at JHH, and from arrival at the spoke site to arrival at JHH, correcting for the distance between spoke and hub.
Statistical Analysis
All statistical analyses were conducted in Stata Version 14.0 (StataCorp, USA). Univariate and multivariate linear regressions were performed to estimate the effect of each variable in the form of a coefficient with respect to DIDO, DI-Decision, and Decision-DO times. Adjusting each model for individual spoke site could not be done, as each spoke site had a small number of patients and almost exclusively used either air transport or road transport, making these effects difficult to separate reliably in this data set. Further, we needed to separate delays associated with air transfers from delays related to sites with reduced experience. Thus, spoke sites were grouped by experience into more experienced sites (those with more than 10 transfers) and less experienced sites (those with 10 or fewer transfers). This would also serve to correct for potentially faster triage and treatment at more experienced sites.
Ethics approval was gained from the Hunter New England Human Research Ethics Committee (HNEHREC Reference No: 13/02/20/5.06 and AU201712-15).
Results
Between April 2013 and October 2019, 1,087 patients were assessed by telestroke, of which 568 (52%) had confirmed ischemic strokes (when mimics and hemorrhages were excluded). Of all patients with confirmed strokes, 175 had LVO (31%). Of these patients with LVO, 75 (43%) were transferred with a view to receiving ECR. Twelve of these patients were transferred to other ECR centers owing to logistic factors relating to transport, weather, or unavailability of ECR at JHH. One patient was transferred to JHH who did not eventually receive ECR owing to established infarct found on repeat brain imaging performed at JHH due to severe unexpected delays during that transfer. A sample of 62 patients who underwent ECR at JHH was finally included in this analysis (Table 1). Mean age was 66 years [standard deviation (SD), 14 years], 34 (55%) were male, and median NIHSS was 13 [interquartile range (IQR), 8 to 18]. There were 42 air transfers (33 by helicopter and 9 by fixed wing) and 20 road transfers. The groups were well-matched overall.
Table 1. Baseline characteristics of different groups.
Transfer Metrics
Median D2D time was 308 min (IQR, 254 to 351). Median DIDO time was 214 min (IQR, 171 to 247). On average, 68% (SD, 10%) of D2D time was spent at the spoke hospital, and only 32% (SD, 10%) in transit. Median DI-Decision time was 70 min (IQR, 54 to 92). Median Decision-DO time was 130 min (IQR, 82 to 182). Median time from ED arrival at spoke sites to first telestroke neurologist contact was 26 min (IQR, 16 to 43).
Comparison of Air and Road Transfers
Spoke sites nearest to JHH transferred patients primarily by road, whereas more distant sites transferred patients primarily by air, with only two patients transferred by the method not normally used by their respective spoke site (Table 2).
Table 2. Distribution of transfers for rural spoke sites.
Air transfer was associated with a delay of 64 min in DIDO time (P = 0.004). This was entirely due to its association with longer Decision-DO times by 87 min (P < 0.001).
Separate univariate regressions of the time from decision to arrival at JHH against distance for air travel and road travel were performed and plotted (Figure 1). Extrapolating from these regressions suggested that the distance at which the extra speed of an aircraft made up for the delays involved in booking an aircraft was 299 km.
Figure 1. Comparison through univariate regression analysis of decision to arrival time at John Hunter Hospital for air and road transfers by distance. Extrapolating from these separate regressions suggested that the distance at which the extra speed of an aircraft made up for the delays involved in booking an aircraft was 299 km.
Other Factors Which Might Influence Workflow
Differences in the level of experience between individual hospitals, that is, between the more experienced and less experienced sites (Table 2), were not significant with respect to key transfer metrics.
Intubation was associated with a delay of 242 min in DIDO and a delay of 229 min in the DI-Decision time but was not significantly associated with the Decision-DO time (P = 0.61, Table 3).
Table 3. Effects of selected variables on components of DIDO in adjusted analysis.
Thrombolysis was associated with faster DI-Decision times by 26 min but slower Decision-DO times by 34 min, which canceled each other out with respect to DIDO times (meaning thrombolysis was not significantly associated with DIDO, P = 0.56).
Discussion
We found that air transfer was associated with a significant delay in DIDO time as compared with road. The main reason for this appeared to be the longer Decision-DO with air travel, which was probably attributable to multiple logistical factors that increased the complexity of organizing air retrieval. Our study indicates that Decision-DO appears to be a key component of DIDO, which has special relevance to the rural telestroke setup. While training and systems are already being used to minimize the DI-Decision time, which was also found to be too long in our study, Decision-DO in particular represents an important target for future optimization of transfer workflow for ECR in our area.
Our analysis suggests that the critical distance at which the speed of air travel makes up for the delays in booking aircraft is ~300 km in our telestroke network, although there is considerable uncertainty about the exact figure, owing to a lack of overlap in the distances covered by air transfers compared with road transfers. Nevertheless, this figure may serve in future as a guide to base decisions about allocation of retrieval resources. For instance, when applied to other health service networks with comparable demographics and infrastructure to ours, transfers below a similar threshold may be prioritized for road rather than air retrieval (24–27).
We discovered that DIDO in our study was 214 min, more than double that published in other metropolitan Australian networks (15, 16). Potential reasons for our relatively long DIDO time include the much larger distances needed to be traversed by retrieval services in our network (hence affecting transport arrival times to rural spoke sites) as well as the resource and workforce limitations inherent in the rural setting.
Intubation was associated with a significant delay in DIDO, whereas thrombolysis was not. Our results are in line with previous studies in which thrombolysis did not have a detrimental impact on workflow of ECR transfers (14–16). However, although thrombolysis had no clear effect on DIDO in our study, it was associated with longer Decision-DO, which was compensated by shorter DI-Decision. This indicates that patients who are identified as potentially suitable for thrombolysis are triaged and treated quickly in our rural spoke sites relative to those who present later. The fact that the Decision-DO time is extended suggests that the process of thrombolysis may actually delay transfers slightly. This effect might be potentially reduced by using tenecteplase instead of alteplase for thrombolysis, which is just a bolus dose and not an infusion.
Although there may be too many variables in the telestroke-guided retrieval process to fully control for, our study raises a few potential roles for decision support in this process, with the overarching goal of reducing delays. One such role might be to integrate the entire retrieval process, so there is a uniform approach for all transfers. For example, once the decision to transfer for ECR has been made in a particular case, a decision support system may work via a platform that helps to coordinate the retrieval process in real time and generate the most appropriate transfer strategy (taking into account factors such as distance, current weather conditions, availability of aircraft and landing stations, and need for intubation before transfer).
Another role for decision support might be in helping to activate the retrieval service. Ideally, the retrieval system could be activated for suspected LVO patients even before the formal decision by the telestroke neurologist. For example, decision support could facilitate ED and ambulance staff to activate the retrieval pathway even prior to consultation with the telestroke neurologist if certain clinical or imaging criteria are met, with the capacity to abort the transfer later if needed. This may serve to reduce Decision-DO in our population, although its implementation is likely to be hindered by transport limitations.
Besides the limited generalizability of our results due to small sample size and heterogeneity of our rural spoke sites, it was beyond the scope of our analysis to measure patient outcomes. Although we did not specifically look at patient outcomes, we instead studied key time metrics associated with DIDO, and by inference, the optimization of these metrics would likely serve to improve patient outcomes. The basis for this is that it is already well-established that delays in DIDO equate to delayed recanalization and hence are detrimental to patient outcomes (12). For similar reasons, we could not adequately compare helicopter and fixed wing retrievals. Some of our data were collected retrospectively. We were unable to account for less predictable factors, which could potentially affect transfers such as the effects of peak-hour traffic, aircraft availability, and refueling protocols.
In conclusion, DIDO in rural areas is longer than described in metropolitan areas. DIDO is longer for air retrievals compared with road. Decision-DO time represents the most important component of DIDO and hence a potential target for future interventions to improve transfer workflow. Decision-DO varies depending upon transport modality, being longer for air travel than road. Systems for rapid transportation of rural ECR candidates need optimization to minimize delays to treatment and ensure best patient outcomes, with decision support seen as a potential tool to achieve this.
Data Availability Statement
The datasets generated for this study will not be made publicly available. The full data-set contains potentially identifiable patient data. Request to access the data can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by Hunter New England Human Research Ethics Committee. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
SG: took the lead in writing the manuscript, in consultation with NS, CG-E, and TL. NS and CG-E: made an equal contribution to the article as joint senior authors and conceived the framework for the project with input from CL and together oversaw the overall operation of the telestroke network. In addition, TL, FM, PG-B, TW, BO'B, JE, KA, CG-E, and NS: collected and collated data for the article. TL conducted the statistical analyses. MP and AB: oversaw the imaging protocols. All authors contributed to designing and writing and the conception of the manuscript on behalf of the Northern NSW Telestroke investigators.
Funding
This study was based on registry information and did not receive any specific funding.
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.
Members (Besides Authors) of Northern NSW Telestroke Investigators
Ms Rachel Peake, Dr James Hughes, Dr Lisa Dark, Dr Nick Ryan, Dr Matt Shepherd, Dr Osama Ali, Dr Hugh Reid, Ms Fiona Minett, Ms Jaclyn Birnie, Ms Amanda Buzio, Dr Iain Bruce, Dr Alan Tankel, Ms Kim Parrey, Dr Matthew Kinchington, Dr Elizabeth Pepper, Dr Andre Loiselle, Dr James Thomas, Dr Jessica Stabler, Dr Mohammad Amin, Ms Michelle Russell, Ms Angela Royan, Mr Brett Roworth, Dr Mary Morgan.s.
References
1. Goyal M, Menon BK, van Zwam WH, Dippel DWJ, Mitchell PJ, Demchuk AM, et al. Endovascular thrombectomy after large-vessel ischaemic stroke: a meta-analysis of individual patient data from five randomised trials. Lancet. (2016) 387:1723–31. doi: 10.1016/S0140-6736(16)00163-X
2. 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
3. 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
4. Holodinsky JK, Patel AB, Thornton J, Kamal N, Jewett LR, Kelly PJ, et al. Drip and ship versus direct to endovascular thrombectomy: the impact of treatment times on transport decision-making. Eur Stroke J. (2018) 3:126–35. doi: 10.1177/2396987318759362
5. Sarraj A, Mlynash M, Savitz SI, Heit JJ, Lansberg MG, Marks MP, et al. Outcomes of thrombectomy in transferred patients with ischemic stroke in the late window: a subanalysis from the DEFUSE 3 trial. JAMA Neurol. (2019) 76:682–9. doi: 10.1001/jamaneurol.2019.0118
6. Saver JL, Goyal M, van der Lugt A, Menon BK, Majoie CBLM, Dippel DW, et al. Time to treatment with endovascular thrombectomy and outcomes from ischemic stroke: a meta-analysis. JAMA. (2016) 316:1279–88. doi: 10.1001/jama.2016.13647
8. Janssen PM, Venema E, Dippel DWJ. Effect of workflow improvements in endovascular stroke treatment: a systematic review and meta-analysis. Stroke. (2019) 50:665–74. doi: 10.1161/STROKEAHA.118.021633
9. Levine SR, Gorman M. Telestroke: the application of telemedicine for stroke. Stroke. (1999) 30:464–9. doi: 10.1161/01.STR.30.2.464
10. Bagot KL, Cadilhac DA, Hand PJ, Vu M, Bladin CF. Telemedicine expedites access to optimal acute stroke care. Lancet. (2016) 388:757–8. doi: 10.1016/S0140-6736(16)31360-5
11. Saler M, Switzer JA, Hess DC. Use of telemedicine and helicopter transport to improve stroke care in remote locations. Curr Treat Options Cardiovasc Med. (2011) 13:215–24. doi: 10.1007/s11936-011-0124-y
12. McTaggart RA, Moldovan K, Oliver LA, Dibiasio EL, Baird GL, Hemendinger ML, et al. Door-in-Door-out time at primary stroke centers may predict outcome for emergent large vessel occlusion patients. Stroke. (2018) 49:2969–74. doi: 10.1161/STROKEAHA.118.021936
13. Ranasinghe I, Turnbull F, Tonkin A, Clark RA, Coffee N, Brieger D. Comparative effectiveness of population interventions to improve access to reperfusion for ST-segment-elevation myocardial infarction in Australia. Circ Cardiovasc Qual Outcomes. (2012) 5:429–36. doi: 10.1161/CIRCOUTCOMES.112.965111
14. Regenhardt RW, Mecca AP, Flavin SA, Boulouis G, Lauer A, Zachrison KS, et al. Delays in the air or ground transfer of patients for endovascular thrombectomy. Stroke. (2018) 49:1419–25. doi: 10.1161/STROKEAHA.118.020618
15. Choi PMC, Tsoi AH, Pope AL, Leung S, Frost T, Loh PS, et al. Door-in-door-out time of 60 minutes for stroke with emergent large vessel occlusion at a primary stroke center. Stroke. (2019) 50:2829–34. doi: 10.1161/STROKEAHA.119.025838
16. Ng FC, Low E, Andrew E, Smith K, Campbell BCV, Hand PJ. Deconstruction of interhospital transfer workflow in large vessel occlusion: real-world data in the thrombectomy era. Stroke. (2017) 48:1976–9. doi: 10.1161/STROKEAHA.117.017235
17. Australian Institute of Health and Welfare. Australia's Health. (2018). Australia's health series no. 16. AUS 221.
18. Saposnik G, Johnston SC. Decision making in acute stroke care. Stroke. (2014) 45:2144–50. doi: 10.1161/STROKEAHA.114.005462
19. Mellor RM, Bailey S, Sheppard J, Carr P, Quinn T, Boyal A etal. Decisions and delays within stroke patients' route to the hospital: a qualitative study. Ann Emerg Med. (2015) 65:279–87. doi: 10.1016/j.annemergmed.2014.10.018
20. Faiz KW, Sundseth A, Thommessen B, Ronning OM. Factors related to decision delay in acute stroke. J Stroke Cerebrovasc Dis. (2014) 23:534–9. doi: 10.1016/j.jstrokecerebrovasdis.2013.05.007
21. Hsia AW, Castle A, Wing JJ, Edwards DF, Brown NC, Higgins TM etal. Understanding reasons for delay in seeking acute stroke care in an underserved urban population. Stroke. (2011) 42:1697–1701. doi: 10.1161/STROKEAHA.110.604736
22. 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. (2020) 26:367–73. doi: 10.1111/cns.13224
23. 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
24. Svenson JE, O'Connor JE, Lindsay MB. Is air transport faster? A comparison of air versus ground transport times for interfacility transfers in a regional referral system. Air Med J. (2006) 25:170–2. doi: 10.1016/j.amj.2006.04.003
25. Diaz MA, Hendey GW, Bivins HJ. When is the helicopter faster? A comparison of helicopter and ground ambulance transport times. J Trauma. (2005) 58:148–53. doi: 10.1097/01.TA.0000124264.43941.41
26. Munford BJ, Beehan SJ. Development of a model for air medical service utilisation. In: ADAC/International Society of Aeromedical Services AIRMED 96 Congress Report. Munich (1996). p. 597–602.
Keywords: telemedicine, decision, retrieval, thrombectomy, rural
Citation: Gangadharan S, Lillicrap T, Miteff F, Garcia-Bermejo P, Wellings T, O'Brien B, Evans J, Alanati K, Levi C, Parsons MW, Bivard A, Garcia-Esperon C and Spratt NJ (2020) Air vs. Road Decision for Endovascular Clot Retrieval in a Rural Telestroke Network. Front. Neurol. 11:628. doi: 10.3389/fneur.2020.00628
Received: 04 March 2020; Accepted: 28 May 2020;
Published: 17 July 2020.
Edited by:
Claus Ziegler Simonsen, Aarhus University Hospital, DenmarkReviewed by:
Christos Krogias, Ruhr University Bochum, GermanyMohammad Anadani, Washington University School of Medicine in St. Louis, United States
Copyright © 2020 Gangadharan, Lillicrap, Miteff, Garcia-Bermejo, Wellings, O'Brien, Evans, Alanati, Levi, Parsons, Bivard, Garcia-Esperon and Spratt. 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: Shyam Gangadharan, shyamg@outlook.com.au
†These authors have contributed equally to this work