Article Text

other Versions

Download PDFPDF

International experience of mechanical thrombectomy during the COVID-19 pandemic: insights from STAR and ENRG
  1. Sami Al Kasab1,2,
  2. Eyad Almallouhi2,
  3. Ali Alawieh3,4,
  4. Michael R Levitt5,
  5. Pascal Jabbour6,
  6. Ahmad Sweid6,
  7. Robert M Starke7,8,
  8. Vasu Saini8,
  9. Stacey Q Wolfe9,
  10. Kyle M Fargen9,
  11. Adam S Arthur10,11,
  12. Nitin Goyal12,
  13. Abhi Pandhi12,
  14. Isabel Fragata13,
  15. Ilko Maier14,
  16. Charles Matouk15,
  17. Jonathan A Grossberg16,
  18. Brian M Howard17,18,
  19. Peter Kan19,
  20. Muhammad Hafeez20,
  21. Clemens M Schirmer21,
  22. R Webster Crowley22,
  23. Krishna C Joshi23,
  24. Stavropoula I Tjoumakaris6,
  25. Shakeel Chowdry24,
  26. William Ares25,
  27. Christopher Ogilvy26,
  28. Santiago Gomez-Paz26,
  29. Ansaar T. Rai27,
  30. Maxim Mokin28,
  31. Waldo Guerrero28,
  32. Min S Park29,
  33. Justin R Mascitelli30,
  34. Albert Yoo31,
  35. Richard Williamson32,
  36. Andrew Walker Grande33,
  37. Roberto Javier Crosa34,
  38. Sharon Webb35,
  39. Marios N Psychogios36,
  40. Andrew F Ducruet37,
  41. Christine A Holmstedt2,
  42. Andrew J Ringer38,
  43. Alejandro M Spiotta1
  44. On behalf of STAR collaborators
  1. 1 Neurosurgery, Medical University of South Carolina, Charleston, South Carolina, USA
  2. 2 Neurology, Medical University of South Carolina, Charleston, South Carolina, USA
  3. 3 Neurosurgery, Emory University, Atlanta, Georgia, USA
  4. 4 Microbiology and Immunology, Medical University of South Carolina, South Carolina, USA
  5. 5 Neurological Surgery, University of Washington School of Medicine, Seattle, Washington, USA
  6. 6 Neurological Surgery, Thomas Jefferson University Hospital, Philadelphia, Pennsylvania, USA
  7. 7 Neurological Surgery, University of Miami Miller School of Medicine, Miami Beach, Florida, USA
  8. 8 University of Miami Miller School of Medicine, Miami, Florida, USA
  9. 9 Neurosurgery, Wake Forest School of Medicine, Winston Salem, North Carolina, USA
  10. 10 Semmes-Murphey Neurologic and Spine Institute, Memphis, Tennessee, USA
  11. 11 Neurosurgery, University of Tennessee Health Science Center, Memphis, Tennessee, USA
  12. 12 Neurology, University of Tennessee Health Science Center, Memphis, Tennessee, USA
  13. 13 Neuroradiology, Centro Hospitalar de Lisboa Central, Lisboa, Portugal
  14. 14 Neurology, University Medicine Goettingen, Goettingen, NS, Germany
  15. 15 Neurosurgery, Yale University, New Haven, Connecticut, USA
  16. 16 Neurosurgery and Radiology, Emory University School of Medicine, Atlanta, Georgia, USA
  17. 17 Neurosurgery, Emory University School of Medicine, Atlanta, Georgia, USA
  18. 18 Radiology and Imaging Sciences, Emory University School of Medicine, Atlanta, Georgia, USA
  19. 19 Department of Neurosurgery, Baylor College of Medicine, Houston, Texas, USA
  20. 20 Department of Neurology, Baylor College of Medicine, Houston, Texas, USA
  21. 21 Neurosurgery, Geisinger Health System, Wilkes-Barre, Pennsylvania, USA
  22. 22 Rush University, Chicago, Illinois, USA
  23. 23 Neurological Surgery, Rush University Medical Center, Chicago, Illinois, USA
  24. 24 North Shore University Health System, Evanston, Illinois, USA
  25. 25 Neurosurgery, North Shore University Health System, Evanston, Illinois, USA
  26. 26 Neurosurgery, BIDMC, Boston, Massachusetts, USA
  27. 27 Department of Neurointerventional Radiology, West Virginia University, Morgantown, West Virginia, USA
  28. 28 Neurosurgery, University of South Florida, Tampa, Florida, USA
  29. 29 Neurosurgery, University of Virginia, Charlottesville, Virginia, USA
  30. 30 Neurosurgery, University of Texas Health Science Center at San Antonio, San Antonio, Texas, USA
  31. 31 Texas Stroke Institute, Fort Worth, Texas, USA
  32. 32 Stroke & Cerebrovascular Center, Allegheny Health Network, Pittsburgh, Pennsylvania, USA
  33. 33 Neurosurgery, Radiology and Neurology, University of Minnesota, Mendota Heights, Minnesota, USA
  34. 34 Endovascular Neurosurgery, Médica Uruguaya, Montevideo, Montevideo, Uruguay
  35. 35 Neurosurgery, Bon Secours St Francis Hospital, Greenville, South Carolina, USA
  36. 36 Department of Neuroradiology, Clinic of Radiology and Nuclear Medicine, University Hospital Basel, Basel, Switzerland
  37. 37 Department of Neurosurgery, Barrow Neurological Institute, Phoenix, Arizona, USA
  38. 38 Neurosurgery, Mayfield Clinic, University of Cincinnati, Cincinnati, Ohio, USA
  1. Correspondence to Dr Alejandro M Spiotta, Neurosurgery, Medical University of South Carolina, Charleston, SC 29425, USA; spiotta{at}


Background In response to the COVID-19 pandemic, many centers altered stroke triage protocols for the protection of their providers. However, the effect of workflow changes on stroke patients receiving mechanical thrombectomy (MT) has not been systematically studied.

Methods A prospective international study was launched at the initiation of the COVID-19 pandemic. All included centers participated in the Stroke Thrombectomy and Aneurysm Registry (STAR) and Endovascular Neurosurgery Research Group (ENRG). Data was collected during the peak months of the COVID-19 surge at each site. Collected data included patient and disease characteristics. A generalized linear model with logit link function was used to estimate the effect of general anesthesia (GA) on in-hospital mortality and discharge outcome controlling for confounders.

Results 458 patients and 28 centers were included from North America, South America, and Europe. Five centers were in high-COVID burden counties (HCC) in which 9/104 (8.7%) of patients were positive for COVID-19 compared with 4/354 (1.1%) in low-COVID burden counties (LCC) (P<0.001). 241 patients underwent pre-procedure GA. Compared with patients treated awake, GA patients had longer door to reperfusion time (138 vs 100 min, P=<0.001). On multivariate analysis, GA was associated with higher probability of in-hospital mortality (RR 1.871, P=0.029) and lower probability of functional independence at discharge (RR 0.53, P=0.015).

Conclusion We observed a low rate of COVID-19 infection among stroke patients undergoing MT in LCC. Overall, more than half of the patients underwent intubation prior to MT, leading to prolonged door to reperfusion time, higher in-hospital mortality, and lower likelihood of functional independence at discharge.

  • stroke
  • thrombectomy
  • complication

This article is made freely available for use in accordance with BMJ’s website terms and conditions for the duration of the covid-19 pandemic or until otherwise determined by BMJ. You may use, download and print the article for any lawful, non-commercial purpose (including text and data mining) provided that all copyright notices and trade marks are retained.

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.


The coronavirus disease 2019 (COVID-19) outbreak, which originated in Wuhan, China, has now spread to more than 200 countries worldwide and led to more than 300 000 deaths.1 The hallmark of infection with COVID-19 is a severe respiratory syndrome with symptoms ranging from mild infection with fever, cough, shortness of breath, and fatigue to severe pneumonia requiring hospitalization and assisted ventilation.2–4 Recent studies have reported neurologic manifestations in approximately 37% of patients with confirmed COVID-19 infection.5 6 A study from China reported ischemic stroke in about 5% of patients with COVID-19 infection, and most recently, a case series from New York reported five cases of large-vessel occlusion (LVO) in young patients, raising the concern for a possible association between COVID-19 infection and LVO in young patients.7 8 Additional studies have linked COVID-19 with a hypercoagulable state, including systemic venous and arterial thromboses.9–11 However, it is unclear whether there is a causal relationship between COVID-19 and LVO in younger patients.

Given the high contagiousness and rapid spread of COVID-19, along with reports of high rates of infection among healthcare workers, multiple national and international societies set forth guidelines based on expert recommendations to help guide stroke providers during the COVID-19 pandemic, in order to mitigate provider illness. Some of the recommendations included modifying the triage workflow and lowering intubation thresholds.12–17 As such, healthcare systems worldwide, in areas with high as well as low burden of COVID-19, adjusted protocols related to mechanical thrombectomy (MT), not necessarily adjusted to the actual local burden of COVID-19. The impact of such changes on MT time metrics and outcomes is unknown, and whether such changes should be implemented in all centers regardless of the COVID-19 burden is also unclear. It is possible that new protocol implementations,5 while aimed at protecting healthcare providers and patients, might not achieve this goal and possibly result in harmful delays and unnecessary use of resources. Furthermore, the prevalence of COVID-19 among patients undergoing MT remains unknown. In this study, we sought to determine the rate of COVID-19 infection among patients undergoing MT in a prospective, multicenter, international observational study. We also aimed to evaluate baseline features among patients undergoing MT during the pandemic and determining the impact of protocol changes on outcomes.


Centers and patient variables

We prospectively collected data from 28 thrombectomy-capable stroke centers from two working collaborations: Stroke Thrombectomy and Aneurysm Registry (STAR) and Endovascular Neurosurgery Research Group (ENRG), including centers in North America, South America, and Europe (figure 1). Consecutive patients who underwent MT for LVO during the peak 2 months of the COVID-19 surge (February–March or March–April 2020, depending on the individual site). We evaluated the rate of COVID-19 in MT patients during the COVID-19 surge period. Testing was performed for symptomatic patients only using real-time reverse transcriptase-polymerase chain reaction of nasopharyngeal swabs.18 Collected data included baseline demographics, admission National Institute Health Stroke Scale (NIHSS), whether or not patients received intravenous alteplase (tPA), time from last known normal (LKN) to hospital arrival, Alberta Stroke Program Early CT score (ASPECTS) for patients with anterior circulation strokes, location of occlusion, door to groin puncture time, thrombectomy procedure duration, number of attempts, complication rate, and rate of symptomatic intracranial hemorrhage. Using data from Johns Hopkins University Center for Systems Science and Engineering COVID-19 tracking project,19 we determined that centers that are located in counties with COVID-19 rate of ≥500 cases per 100 000 people (as of May 1 2020) were considered high COVID-19 burden counties (HCC), and centers with <500 cases per 100 000 people were considered low COVID-19 counties (LCC). Successful revascularization was defined as modified thrombolysis in cerebral infarction 2B-3. Symptomatic intracranial hemorrhage (sICH) was as worsening of 4 or more points in NIHSS attributed to hemorrhagic transformation,20 and favorable discharge outcome was defined as modified Rankin scale (mRS) of 0–2 on discharge. Procedural data was adjudicated by the investigators in each of the included centers.

Figure 1

Locations of the included centers in this study.

Statistical analysis

We used descriptive statistics to report patient demographic and clinical characteristics using median and IQR for continuous variables and percentages for categorical variables. Characteristics of groups were compared using the Wilcoxon rank-sum (Mann–Whitney U) test, Fisher’s exact test, and chi-square as appropriate. To evaluate the effect of COVID-19 infection on the discharge outcomes of MT patients, we used inverse probability of treatment weight (IPTW) using propensity scores to estimate a generalized linear model with a logit link function. The propensity scores were calculated using a binary logistic regression model for COVID-19 status using age, NIHSS, LKN to groin puncture, and location of LVO as predictors. Covariates in the propensity score model that had a less than 0.20 standardized difference before and after weighting were considered well-matched. Finally, to assess the relationship between intubation prior to MT and favorable discharge outcome (mRS 0–2), we estimated a generalized linear model with logit link to control for confounders including age, door to reperfusion time, tPA status, admission NIHSS, and location of LVO. All final models were tested for collinearity. An alpha level of 0.05 was used as the level of statistical significance. The analysis was conducted using SPSS v25 (IBM Corporation, New York, NY).


Patient population

During the COVID-19 surge period (February–April 2020), 458 patients underwent MT in 28 centers from five countries (USA, Germany, Switzerland, Portugal, and Uruguay) (figure 1). The median age was 71 (IQR 60–80) years' old, admission NIHSS 16 (10-21), LKN to hospital door duration 260 (IQR 119–506), and ASPECTS 8 (IQR 7–10). Our cohort included 210 (45.9%) females, 310 (67.7%) white patients, 46 (10%) patients with posterior circulation occlusion, and 184 (40.2%) patients who received IV tPA (table 1). Thirteen (2.8%) patients were confirmed COVID-19 positive out of 242 MT patients who were tested; two were confirmed positive on admission and 11 during the same hospitalization. Of the 242 patients that were tested, 50 were in HCC, of which nine (18%) were positive, and 192 were LCC of which 4 (2.1%) were positive (P<0.001).

Table 1

Patient, procedural, and outcome metrics of stroke patients treated with mechanical thrombectomy at the 28 international sites and divided by COVID-19 status

COVID-19 positive vs non-COVID-19 patients

In COVID-19 positive patients, the median age was 58 (IQR 50–71) with a range of 35–98, and three patients were younger than 50 years of age. COVID-19 patients had a higher NIHSS on admission (19 (IQR 16–24) vs 15 (IQR 10–20), P=0.047), otherwise there was no difference in the percentage of female patients (38.5% vs 46.1%, P=0.779), percentage of white patients (46.2% vs 68.3%, P=0.129), or LKN to hospital door duration (202 (IQR 92–606) vs 260 (IQR 119–506), P=0.557), in-hospital mortality rate (33.3% vs 24.1%, P=0.496), length of stay (8 (IQR 4–15) vs 5 (IQR 3–10), discharge mRS (4 (IQR 3–6) vs 4 (IQR 2–5), P=0.443), or rate of functional independence (mRS 0–2) at discharge (2 (16.7%) vs 94 (29.7%), P=0.52) between COVID-19 and non-COVID-19 patients. The full characteristics of MT patients with confirmed COVID-19 are summarized in table 1. On multivariate analysis using logistic regression model with IPTW and propensity score, COVID-19 status was not associated with in-hospital mortality (RR 2.19, 95% CI 0.585 to 8.199, P=0.244) or functional independence (mRS 0–2) on discharge (RR 0.493, 95% CI 0.101 to 2.412, P=0.383).

High vs low prevalence COVID-19 counties

Five centers with 104 patients treated with MT met the criteria for HCC, and the remaining 354 patients in 23 centers were considered LCC. There was a higher number of patients who tested positive for COVID-19 infection in HCC (9 (8.6%) vs 4 (1.1%), P≤0.001). Six of those positive COVID-19 cases inHCC came from a single institution where the incidence was 6/30 (20%). Of the nine patients from HCC, two were positive on admission and the remaining tested positive during admission. All four patients in LCC tested positive during admission. Table 2 summarizes baseline features and outcomes between HCC and LCC groups. Notably, there was a trend toward longer door to groin time in the LCC group (85 (45-127) min vs 68 (36-123) min, P=0.054). In addition, more patients received tPA in the LCC group (29 (27.9%) vs 155 (43.8%), P=0.004). There was no significant difference in any other variables between the two groups.

Table 2

Clinical characteristics and outcomes of patients in high and low COVID-19 counties

General anesthesia vs monitored anesthesia care

A total of 243 (53.1%) were intubated and underwent MT under general anesthesia (GA) (241 patients were intubated prior to the procedure and two patients required intubation during the procedure). The remaining 215 patients underwent MT with monitored anesthesia care (MAC) (n=215). Table 3 summarizes baseline demographics and outcomes in the GA and MAC groups. Compared with the MAC group, GA patients were younger (70 (IQR 58–79) vs 73 (IQR 63–82) years' old, P=0.02), more likely to have received tPA (110 (45.6%) vs 73 (34%), P=0.011), had higher rate of posterior circulation LVOs (31 (12.9%) vs 15 (7%), P=0.037), longer door to groin puncture time (92 (IQR 62–136) vs 67 (35-112) minutes, P<0.001), longer procedure duration (47 (IQR 92–188) vs 32 (19-56) minutes, P<0.001), and longer overall door to reperfusion time (138 (IQR 92–188) vs 100 (IQR 64–156) min, P≤0.001). Less patients in the GA group were functionally independent (mRS 0–2) at discharge compared with the MAC group (45 (24.5%) vs 51 (35.4%), P=0.03). There was no difference in any other variable between the two groups.

Table 3

Clinical characteristics and outcomes of patients who underwent mechanical thrombectomy under general anesthesia vs monitored anesthesia care

On multivariate analysis, intubation status was associated with higher probability of in-hospital mortality (RR 1.871, 95% CI 1.066 to 3.284, P=0.029) and lower probability of functional independence on discharge (0.53, 95% CI 0.318 to 0.884, P=0.015) after controlling for age, door to reperfusion duration, tPA status, admission NIHSS, and location of occlusion. Of note, discharge mRS was missing for 57 (23.6%) patients in the intubated group and 71 (33%) in the non-intubated group.


Several studies have evaluated the incidence of neurologic symptoms in patients with COVID-19. A study by Mao et al evaluated 214 patients with confirmed COVID-19 infection at three centers in China.5 In this study, 78 (36.4%) patients had neurologic manifestations. Acute cerebrovascular diseases were observed in five patients. Another Chinese study reported a similar incidence of acute cerebrovascular disease of 5.8% (13/221).21 Of these, 11 were arterial ischemic stroke, one was venous sinus thrombosis, and one was intracerebral hemorrhage. However, to date there are no reports regarding the incidence or prevalence of COVID-19 among patients presenting with emergent LVO. In this prospective, observational, multicenter international study we evaluated the rate of COVID-19 infection among patients undergoing MT: we also evaluated the impact of the COVID-19 pandemic on MT-related time metrics and clinical outcomes. We found that 2.8% of patients undergoing MT tested positive for COVID-19. In HCC, approximately 8.6% of patients undergoing MT tested positive for COVID-19 compared with 1.1% in LCC.

Cerebrovascular complications are reported in approximately 5% of patients admitted with COVID-19.7 Ischemic stroke in the setting of COVID-19 might occur concurrently either from common stroke risk factors such as hypertension, diabetes mellitus, hyperlipidemia, and atrial fibrillation, as a hypercoaguability- related complication due to COVID-19, a proinflammatory state, or due to decompensation in patients with underlying cerebrovascular disease such as extracranial or intra-cranial stenosis.7 8 22 23 As such, a recent report from New York raised the concern of a possible association between COVID-19 and LVO in young patients.8 Among the five reported patients however, three had vascular risk factors.8 In our study, while COVID-19 patients were relatively younger than non-COVID patients undergoing MT (58 vs 72), only 3/13 patients with confirmed COVID-19 infection were younger than 50 and the youngest was 35. While it is possible that COVID-19 could be associated with acute ischemic stroke in young, previously healthy patients, there is no evidence that this is a widespread phenomenon.

With the goal of protecting healthcare providers and optimizing neurologic outcomes among patients undergoing MT, and in light of potential and severe personal protective equipment shortages, various neuro-interventional societies have published recommendations based on expert consensus to guide interventionists performing MT. Some of the recommendations, especially related to the type of anesthetic use (MAC vs GA), required drastic changes in local protocols. The Society of Neurointerventional Surgery recommended a lower threshold for intubation prior to MT in patients who are either positive for COVID-19 or possibly could be positive (unknown or unconfirmed negative status).12 Similarly, the Society of Neuroscience in Anesthesiology & Critical Care recommended lower threshold for intubating patients with unknown COVID-19 status.16 While the intubation itself places the providers performing the intubation procedure at risk, the rationale for lowering thresholds for intubation involved protecting against unplanned aerosol generation, for example, a patient who decompensates during the MT procedure and requires intubation. On the other hand, the Society for Vascular and Interventional Neurology recommended considering MAC if possible to protect anesthesiologists from exposure and to reduce unnecessary intubation as well as conserve mechanical ventilator resources.13 Thus, there remains considerable uncertainty regarding ‘best practices’ as the pandemic progresses across the world.

Somewhat paradoxically, significantly higher number of patients were intubated in LCC prior to MT, which was associated with longer door to reperfusion times. Evidence from clinical trials indicates that GA is non-inferior to MAC in patients undergoing MT24–26 and in clinical practice there is considerable variability. In this study, 241/458 (52.6%) of the patients underwent planned pre-procedure GA, and only 2/194 (1%) required emergent, unplanned intra-procedural intubation. We found an approximately 25 min delay in door to groin times in patients undergoing GA compared with MAC. Similarly, door to reperfusion time was approximately 38 min longer in GA patients. In this study, despite relatively younger age for GA patients and higher rate of intravenous thrombolysis, independence at discharge was significantly lower in GA patients (24.5% vs 35.4%). By contrast, The Anesthesia During Stroke (AnStroke) and the General Or Local Anesthesia in Intra Arterial Therapy (GOLIATH) trials both showed no difference in long-term outcomes between patients undergoing MT with GA compared with MAC.27 28 However, in both trials there was no difference in door to groin puncture time or symptom onset to reperfusion duration between the GA and MAC groups. Similarly, there was no difference in time from symptom onset to reperfusion. A post hoc analysis of data collected in the MR CLEAN trial demonstrated an approximately 32 min delay with GA which was associated with worse outcome.29

In our study the majority of sites were not intubating most patients for MT prior to the pandemic, but rather did so following scientific society recommendations. It is therefore possible that the worsened outcomes in the GA group were as a direct result of the delay in door to reperfusion times compared with MAC patients. Changes in local protocols take time to become streamlined and when first implemented are often not optimized into the patient workflow. This consideration was likely exacerbated by other pandemic-related protocols affecting supporting services. Unlike centers that routinely perform MT under GA, centers that primarily perform MT under MAC are likely to encounter longer times to intubate and therefore delayed door to groin puncture times, as identified in this study. Beyond the actual intubation procedure itself, local protocols commonly dictated transporting patients to designated negative pressure rooms and mandated wait times prior to transporting up to the angiography suite. It is important to note that all these changes in patient care workflow likely contributed to the observed effect.

Strengths and limitations

This study’s strengths include the relatively large volume of patients and the prospective, multicenter nature of the data. The main limitations of the study include lack of long-term functional outcomes and missing discharge mRS for a considerable number of patients given that the study was undertaken between the months of January and April of this year and many of the included patients were still admitted when the analysis was conducted. Another limitation is that not all MT patients in our study were tested for COVID-19 which may have resulted in missing some of asymptomatic COVID-19 patients. Also, we cannot rule out false negativity in symptomatic patients either. In addition, this study cannot provide data related to whether COVID-19 is associated with LVO, particularly in the young especially given that HCC patients were under-represented in this cohort. This multicenter experience, however, provides preliminary demographic data on patients undergoing MT.


In this prospective, multicenter, international study, we found a low number of confirmed COVID-19 infections among stroke patients with LVO undergoing MT in LCC. Intubation prior to MT during the early stages adjustment to COVID-19 pandemic was associated with a higher in-hospital mortality rate and lower rate of functional independence on discharge.



  • SAK and EA are joint first authors.

  • Twitter @DrMichaelLevitt, @PascalJabbourMD, @AhmadSweidMD, @Starke_neurosurgery, @AdamArthurMD, @abhipandhi, @BrianHoward_MD, @brainvesseldoc, @SantiagoGP92, @alex_spiotta

  • Collaborators Dileep R Yavagal MD, Eric C Peterson MD, Daniel Raper MBBS, Patrick A Brown MD, M. Reid Gooch MD, Nabeel Herial MD, Ajith Thomas MD, Justin Moore MD, Felipe Albuquerque MD, Louis J. Kim MD, Melanie Walker MD, Michael Chen MD, Stephan Munich MD, Daniel Alan Hoit MD MPH, Violiza Inoa-Acosta MD, Christopher Nickele MD, Lucas Elijovich MD, Fernanda Rodriguez-Erazú PD Dr. med. Jan Liman, Michael Cawley, MD, Gustavo Pradilla MD, Brian Walcot, MD, Zeguang Re, MD PhD, Ryan Hebert MD, João Rei, MD, Jaime Pamplon, Rui Carvalho MD, Mariana Baptist, MD, Ana Nunes MD, Russell Cerejo MD, Ashis Tayal MD, Parita Bhuv, MD, Paul Hansen MD, Norman Ajiboye MD, Alex Brehm, MD

  • Contributors All authors have: provided a substantial contribution to the conception and design of the studies ‎and/or the acquisition and/or the analysis of the data and/or the interpretation of the data. ‎They have drafted the work or revised it for significant intellectual content and approved the ‎final version of the manuscript. They agree to be accountable for all aspects of the work, ‎including its accuracy and integrity.‎

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial, or not-for-profit sectors.

  • Map disclaimer The depiction of boundaries on this map does not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction, or area or of its authorities. This map is provided without any warranty of any kind, either express or implied.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement Data are available upon reasonable request. Additional data from this project can be acquired by contacting the ‎corresponding author.‎