Abstract

Background. Coronavirus disease 2019 (COVID-19) had a significant impact on the National Health Service in the United Kingdom (UK), with over 35 000 cases reported in London by July 30, 2020. Detailed hospital-level information on patient characteristics, outcomes, and capacity strain is currently scarce but would guide clinical decision-making and inform prioritisation and planning. Methods. We aimed to determine factors associated with hospital mortality and describe hospital and ICU strain by conducting a prospective cohort study at a tertiary academic centre in London, UK. We included adult patients admitted to the hospital with laboratory-confirmed COVID-19 and followed them up until hospital discharge or 30 days. Baseline factors that are associated with hospital mortality were identified via semiparametric and parametric survival analyses. Results. Our study included 429 patients: 18% of them were admitted to the ICU, 52% met criteria for ICU outreach team activation, and 61% had treatment limitations placed during their admission. Hospital mortality was 26% and ICU mortality was 34%. Hospital mortality was independently associated with increasing age, male sex, history of chronic kidney disease, increasing baseline C-reactive protein level, and dyspnoea at presentation. COVID-19 resulted in substantial ICU and hospital strain, with up to 9 daily ICU admissions and 41 daily hospital admissions, to a peak census of 80 infected patients admitted in the ICU and 250 in the hospital. Management of such a surge required extensive reorganisation of critical care services with expansion of ICU capacity from 69 to 129 beds, redeployment of staff from other hospital areas, and coordinated hospital-level effort. Conclusions. COVID-19 is associated with a high burden of mortality for patients treated on the ward and the ICU and required substantial reconfiguration of critical care services. This has significant implications for planning and resource utilisation.

1. Introduction

Coronavirus disease 2019 (COVID-19), an infectious syndrome caused by SARS-CoV-2, appeared in December 2019 and evolved into a pandemic that caused more than eleven million cases and 530 000 deaths worldwide by July 2020 [1]. The high number and acuity of patients resulted in an unprecedented demand for hospitalisation and critical care services in many affected countries. Developed areas such as Wuhan (China), Lombardy (Italy), and New York (United States of America) reported a surge in critically ill patients, which quickly led to significant strain on healthcare systems through shortages in the Intensive Care Unit (ICU) beds, equipment, and trained personnel [24].

There was concern that the United Kingdom (UK) would face similar challenges, particularly in densely populated areas like London. The first infection was reported on January 30, 2020, and by late July, the country had recorded a large number of cases, with over 35 000 confirmed infections in the greater area of London alone [5].

Despite the high number of hospital admissions with COVID-19, the existing peer-reviewed literature in the UK remains restricted to large population-level studies and small retrospective cohorts with few details on clinical management and hospital-level strain [68]. Major reports published so far focus either on large but incomplete hospital cohorts with missing data and short follow-up time [7] or exclusively on ICU admissions [8] with no connection between the two populations. Extrapolation from international settings is difficult due to significant differences in population characteristics [2], health system organisation [2, 4, 9], and perceived health system strain [10].

Detailed local patient-level information on characteristics and outcomes, as well as institution-level information on service pressures, would guide clinical decision-making and inform effective prioritisation and resource allocation in the future, particularly in circumstances when a surge in demand places the National Health Service again at risk of being overwhelmed. The aim of our study was to describe the clinical characteristics and course of hospitalised patients with COVID-19, as well as its broader resource implications for a tertiary academic hospital in London, UK.

2. Methods

2.1. Study Design

We conducted a prospective cohort study at King’s College Hospital (KCH), a tertiary academic centre in London, UK. Institutional (IRAS-256619; April 15, 2020) and regional (Health Research Authority; IRAS-256619; April 22, 2020) review board approvals waived the need for ethics committee review and the need for informed patient consent.

All consecutive patients tested for SARS-CoV-2 infection using reverse-transcriptase polymerase chain reaction (RT-PCR) assays of respiratory tract samples between February 25 and March 31, 2020, were considered eligible for study participation. Inclusion criteria were age of 18 years or above and laboratory-confirmed SARS-CoV-2 infection, which was defined as at least one positive RT-PCR nasopharyngeal swab [11]. Patients with missing identifiers and missing SARS-CoV-2 test results and those transferred from other hospitals were excluded. Included patients were followed up until death, hospital discharge, or 30 days after hospital admission. Follow-up was concluded on April 30, 2020. For patients with multiple hospital or ICU admissions, only the first admission was recorded. Transfers between different ICUs within the hospital were considered part of the same admission. The primary outcome was mortality at hospital discharge or at 30 days. Secondary outcomes were ICU mortality, as well as hospital and ICU capacity strain, measured in bed occupancy. The report of our findings is based on the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement [12]. More details regarding the institutional setting and pandemic surge are available in Supplementary File 1.

2.2. Data Collection and Management

Trained members of the clinical team extracted anonymised data from electronic health records (EHRs) [13]. We recorded the following for all included patients: age, sex, ethnicity, area-level socioeconomic deprivation (Index of Multiple Deprivation [14]) clinical frailty, medical comorbidities, Age-adjusted Charlson Comorbidity Index (ACCI), body mass index (BMI), prior residence, self-reported presenting symptoms, reason for hospital admission, laboratory, microbiological, and imaging tests, COVID-19-specific drug treatments, treatment limitations such as Treatment Escalation Plans (TEPs) or do-not-resuscitate (DNR) orders, length of hospital stay, and outcome at hospital discharge or at 30 days.

For patients who were initially admitted to the ward, we also recorded peak temperature and National Early Warning Score version 2 (NEWS2) score [15], maximum level of respiratory support, and development of complications such as hypoxia, hypotension, tachycardia, and depressed consciousness. For patients who were admitted to the ICU, we additionally recorded reason for ICU admission, daily Sequential Organ Failure Assessment (SOFA) score [16], organ support (mechanical ventilation, renal replacement therapy, and circulatory support), arterial partial pressure of oxygen (PaO2), fraction of inspired oxygen (FiO2), PaO2/FiO2 ratio (PFR), use of medications such as vasopressors, pulmonary vasodilators, and neuromuscular blocking drugs, use of prone positioning or extracorporeal membrane oxygenation (ECMO), insertion of tracheostomy, length of ICU stay, and outcome at ICU discharge or at 30 days.

We recorded comorbidities and symptoms based on electronic case note review and categorised variables according to clinical relevance and current literature [4, 10, 17]. We used only validated laboratory results and official reports of imaging studies. Initial tests for newly admitted COVID-19 patients refer to those performed within 24 hours of hospital admission; for patients already admitted, they refer to tests performed within 24 hours of COVID-19 diagnosis. Clinical frailty was assessed on a 9-category scale [18], using information available in the EHR; patients with a score above 4 were considered frail [19]. Treatment escalation plans (TEPs) refer to structured assessments of patients’ suitability regarding specific aspects of treatment such as organ support or ICU admission [20]. Acute kidney injury (AKI) was defined according to Kidney-Disease Improving Global Outcomes (KDIGO) criteria [21]. More details regarding definitions of collected data are provided in Supplementary File 2.

2.3. Statistical Methods

Descriptive analyses are presented as median (IQR) or number (%), and we avoided univariate comparisons between groups. To identify factors associated with the time to death at hospital discharge or at 30 days, we performed multivariable Cox proportional hazards and parametric survival analyses. Covariate inclusion followed a structured approach. We tested the proportional hazards assumption and investigated interactions of the included covariates with sex and age. We also assessed whether the baseline survival experience differed by categories of age, sex, frailty, and ethnicity, using stratified Cox regression. The effect of each included covariate was quantified by calculating adjusted hazard ratios (HRs) with 95% confidence intervals (95% CI). We followed a similar approach for the parametric analysis and plotted the hazard function over time for six hypothetical patients, in order to show the effect of each included covariate on the hazard of death. Statistical tests were 2-sided, with an α-level of 0.05 for statistical significance. We did not impute any missing data. Analyses were performed using Stata/MP version 15.1 (StataCorp). Details regarding the statistical approach are provided in Supplementary File 3.

3. Results

Between February 25 and March 31, 2020, 2728 patients were tested for SARS-CoV-2 infection at KCH. After excluding 170 patients (6%), 2558 patients (94%) were screened for inclusion; 2129 (83%) of them did not meet the inclusion criteria, and 429 (17%) were included in the study. Among them, 353 patients (82%) were treated only on the ward and 76 (18%) were treated in the ICU. The study flow diagram is shown in Figure 1.

Baseline patient characteristics are shown in Table 1. Patients had a median age of 65 years (IQR 52–81), and 61% were of ethnic minority background. The majority (61%) were overweight and 90% of them had at least one serious comorbidity; the most common were hypertension (52%), diabetes (37%), and chronic kidney disease (CKD) (14%). The most common presenting complaints included cough (62%), fever (62%), and dyspnoea (43%), but only 42% of patients were hospitalised for hypoxaemic respiratory failure. Initial chest X-rays were normal in 23% of patients and showed diffuse, bilateral infiltrates in more than half.

Compared to ward patients, those admitted to the ICU were younger and less commonly frail and had a lower burden of chronic comorbidities, as described by the ACCI. They were, however, more likely to be diabetic and to be hospitalised for respiratory failure, with significantly more deranged initial laboratory tests and abnormal chest imaging.

Admission laboratory findings for all patients are shown in Table 2. Lymphopenia was common, as were elevations in C-reactive protein (CRP) (median 75 mg/L, IQR 28–143) and creatinine (median 92 mg/L, IQR 68–130). Raised lactate dehydrogenase, ferritin, creatine kinase, and D-dimers were also common, albeit measured in a minority of patients.

3.1. Patients Treated on the Ward

Most patients developed significant morbidity during their hospital admission. Fifty-two percent met the criteria for critical care outreach team activation (NEWS2 score above six), predominantly for respiratory failure (Table 3). On chest computed tomography (CT), more than 70% of scans showed bilateral diffuse infiltrates and 25% revealed pulmonary embolism (PE). PE was more frequently identified in ICU patients (33%). Subsequent laboratory tests revealed worsening lymphopenia, raised CRP (median 159 mg/L, IQR 81–299) and intracellular enzymes, and impaired hepatic, renal, and haemostatic function (Table 2). Overall, 21% of patients developed KDIGO stage three AKI, but the incidence among ICU patients (67%) was much higher than that among ward patients (11%). We recorded only mild degrees of myocardial involvement, evidenced by small increases in values of high-sensitivity troponin T and unremarkable echocardiographic findings. Secondary infections were common in the ICU-treated group: 50% had a positive respiratory tract sample and 41% had a positive blood culture. Samples from the respiratory tract were frequently positive for Gram-negative organisms and fungi.

Treatment limitations were placed in 61% of patients overall, most commonly on hospital admission (median 0 days, IQR 0–1). The majority (>80%) of these limitations involved TEPs and DNR orders. Two-thirds of patients with a TEP were considered frail. Among patients with treatment limitations that were not considered frail, many had significant medical comorbidities such as active malignancy and stroke.

3.2. Patients Treated in ICU

The clinical trajectory of 76 patients who were admitted to the ICU is described in Table 4. Many patients developed precipitous respiratory failure and required ICU admission directly from the emergency department (57%). Multiple organ failure was common with 97% requiring mechanical ventilation for hypoxic respiratory failure, 80% requiring pharmacologic circulatory support, and 57% requiring renal replacement therapy (RRT) during their ICU stay. The severity of illness was reflected in the high SOFA scores, which remained elevated even after two weeks of ICU stay. The median duration of mechanical ventilation was 12 days (IQR 6–23) and that of renal replacement was 11 days (IQR 4–17). Patients frequently required rescue oxygenation strategies (54%), which included neuromuscular blocking drugs, inhaled prostacyclin, and prone positioning. Extracorporeal membrane oxygenation (ECMO) was used only in 5% of ICU patients. A significant proportion (38%) of patients required a tracheostomy, which was performed at a median of 16.5 days (IQR 14–20) after ICU admission. Almost all laboratory parameters were more deranged among ICU patients, compared to those treated on the ward (Table 2).

3.3. Outcomes

Unadjusted mortality at hospital discharge or at 30 days was 26% overall, 23% for ward patients, and 38% for those treated in the ICU. In unadjusted Cox survival analysis, hospital mortality was associated with age, sex, ACCI score, frailty category, CRP, creatinine, CKD, diabetes, and dyspnoea or fever as presenting complaints. After adjustment, it was independently associated with increasing age (HR 1.07 per decade above 40 years; 95% CI 1.04–1.09, ), male sex (HR 2.31; 95% CI 1.52–3.50, ), raised admission CRP level (HR 1.03 per 10 mg·l−1 increments above the upper normal limit of 5 mg·l−1; 95% CI 1.01–1.04, ), history of CKD (HR 1.87; 95% CI 1.21–2.89, ), and dyspnoea as a presenting symptom (HR 1.88; 95% CI 1.24–1.86, ). Ethnicity or level of deprivation was not associated with mortality in unadjusted or adjusted analyses. Stratification of the model by categories of age (≤60, >60 years), sex, frailty, or ethnicity did not provide evidence of different baseline hazards. The effect of each included covariate, based on parametric modelling, is shown in Figure 2, with six examples of hypothetical patients. In the parametric model, male sex and dyspnoea on presentation had the largest impact on the hazard function. More details regarding the results of the parametric analysis are available in Supplementary File 3.

At the end of the 30-day follow-up, 15 of 353 ward patients (4%) remained admitted in hospital and 12 of 270 (4%) were readmitted after hospital discharge. Among all ICU patients, 13 (17%) remained admitted in the ICU at the end of follow-up. Among the 50 patients discharged from the ICU, three (6%) were readmitted to the ICU, 25 (50%) were discharged from the hospital, three (6%) died on the ward, and 22 (44%) remained hospitalised on the ward.

During the study period, the hospital experienced a sudden surge in critically ill patients, and details regarding the number of hospital and ICU admissions are shown in Figure 3. Over a period of approximately six weeks, 20–40 patients were admitted daily to the hospital with COVID-19 and between five and ten of them required admission to the ICU for organ support. At the peak of the pandemic, both the hospital and ICUs experienced significant capacity strain, with 28% of the hospital capacity of 900 beds and 70% of the ICU surge capacity of 129 beds taken up by patients with COVID-19. This increased demand was supported by changes to the service described in Table 5. Major changes included anaesthetic cover for newly opened ICU beds, re-deployment of the entire hospital workforce to support ICU, delivery of critical care in nonconventional areas like operating rooms, and the introduction of teams dedicated to tasks like prone positioning, communication with families, and tracheostomies. More details are available in Supplementary File 4.

4. Discussion

We describe data from 429 patients hospitalised with COVID-19 in an academic hospital in London UK, of which 76 (18%) were treated in the ICU. We observed a disproportionate burden of COVID-19 hospitalisation and ICU admission in patients from an ethnic minority background, consistent with reports from similarly diverse areas [9] and the rest of the UK [7, 8, 22]. The range of clinical presentations was broad, similar to reports of the severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) [23, 24]. Respiratory involvement was very frequent, but not universally present. Hence, the reliance on the initial presence of respiratory failure or abnormal chest imaging results to make decisions regarding patient infectiousness and the requirement for isolation may be problematic.

Overall hospital mortality was identical to that reported in a recent large UK observational study [7], but patients admitted to the ward had higher morbidity and mortality than previously reported [4, 17, 25]. A significant proportion experienced clinical deterioration with cardiorespiratory compromise and AKI, meeting the criteria for critical care outreach activation. Patients treated on the ward required considerable input from critical care services, with the ICU outreach team seeing almost four out of ten admitted patients. Similar patients may have been treated in the ICU in other described cohorts, and this may explain the difference in unadjusted mortality.

The contribution of the critical care outreach team, as evidenced by the number of patients seen, was crucial in our cohort; however, no comparison data exist. As part of their role, outreach teams commonly bridge the gap between ward and ICU by providing expertise, monitoring, and interventions outside the ICU environment and, importantly, by engaging in philosophy-of-care discussions with patients and their surrogate decision-makers [26]. In the context of many fixed constraints such as the number of ICU beds, the outreach service is a delivery model that greatly expands the availability of expertise and advanced treatments on the ward. Depending on the institutional context, however, outreach services require appropriate equipment and trained personnel such as doctors, nurses, and allied health professionals [26]. This has significant implications on staff planning and resource allocation during a pandemic.

The majority of ward-treated patients in our study had a decision to limit the escalation of therapy recorded during their admission, which represents a departure from previous practice. Most patients had treatment limitations placed upon admission. In addition to being clinically frail, serious comorbidities such as active malignancy and stroke were more common in these patients. Within the framework of generic national guidance regarding treatment escalation and ICU admission [27], we believe that the COVID-19 pandemic promoted more proactive communication about goals of care between physicians and their patients. In fact, Palliative Care Service in our institution faced a dramatic increase in ward referrals over the study period [28], but comparative data from other UK settings are lacking. The resource implications of delivering a comprehensive palliative care response during a pandemic are unknown but likely to be substantial [2931].

Patients treated in the ICU had high morbidity and mortality. Few of them had trials of CPAP or NIV, and almost all required invasive mechanical ventilation, which differs from the experience in the UK [7] and other countries [10]. Pulmonary involvement was significant, with longer periods of mechanical ventilation and more frequent recourse to rescue oxygenation therapies than previously described [9]. Extrapulmonary organ involvement was also common and led to higher utilisation of pharmacological circulatory support and RRT than that reported across the UK [8] and in recent COVID-19 cohorts from Italy (27.8%) [32], the USA (31%) [9], and China (25%) [33]. The high proportion of patients need organ support in the ICU could be explained by the fact that many patients with less severe organ failure were treated on the ward, with the aid of the outreach service. Similarly, the high incidence of significant extrapulmonary involvement supports the impression that critically ill COVID-19 patients are more similar to those with MERS [24], whereas SARS was a predominantly respiratory disease with single organ failure [34].

Hospital mortality in mechanically ventilated COVID-19 patients was similar to that of unselected cohorts of severe adult respiratory distress syndrome (ARDS) patients who required ventilation (46.1%), as well as that of mechanically ventilated patients with SARS (45%) and MERS (52.4%) [24, 34, 35]. In our survival analysis, we identified age, male sex, CKD, elevated CRP, and dyspnoea at presentation as important baseline prognostic factors associated with 30-day mortality after hospitalisation. These are broadly similar to patient characteristics described in reports from other countries [9, 36, 37] and the UK [7]. We did not, however, identify an independent association with ethnicity, which is also in line with findings from a recent large US cohort study [38]. Early identification of patient characteristics associated with poor clinical outcomes and higher resource utilisation is important when faced with the need to prioritise patients in the face of high demand for critical care services, both in the ICU and on the ward.

The COVID-19 pandemic led to a sudden surge in the overall number of critically ill patients across the hospital, as well as in the intensity and duration of treatments they required. This resulted in significant strain for the ICU in terms of capacity and service delivery. ICU bed capacity was effectively doubled; critically ill patients on the ward required frequent input by the outreach team, while patients in the ICU required prolonged mechanical ventilation and RRT, frequent rescue oxygenation therapies, and tracheostomy. To support this increased demand for critical care services, several operational changes were instituted, with support from Anaesthetic Department. These were based on local adaptation of lessons from previous infectious disease outbreaks and information from countries that had been already affected by COVID-19 [2, 3, 39]. The implementation of this service delivery model required the provision of a considerable amount of education and training. Also, these changes were implemented in a short time window and required collaboration across the entire hospital. Our lesson from this process is that clinicians and health administrators need to consider a rapidly scalable model, as health systems could easily become overwhelmed, and our experience is directly comparable to that from other, similarly affected, metropolitan areas [40].

Our study has several strengths. First, it includes a clearly defined cohort with near-complete 30-day outcomes, thus minimising selection bias. Second, it provides a granular clinical description of patient trajectories with previously unknown information regarding ICU outreach involvement and palliative care practice. Third, the inclusion of patients treated both on the ward and in the ICU offers a more complete picture of the organ dysfunctions experienced and the corresponding resource implications. Finally, we provide unique details regarding hospital and ICU capacity strain and service reorganisation.

Our study results must, however, be viewed in light of its methodological limitations. The prospective design relied on abstracting data from the EHR and is susceptible to missing data and recording bias. We identified a number of factors associated with mortality but do not imply any causal relationships. Finally, this study was performed in a single centre in the UK and this may limit its generalisability to other institutional contexts.

5. Conclusions

Our study identified age, sex, CKD, baseline CRP, and respiratory involvement as factors associated with hospital mortality in COVID-19 patients admitted to a large academic hospital in London, UK. Patients treated on the ward and in the ICU had a significant extrapulmonary disease, frequent treatment limitations, and a high burden of mortality. Medical wards caring for COVID-19 patients experienced a substantially increased workload and required frequent ICU outreach input. In order to provide effective care under pandemic surge conditions, significant reorganisation took place within the hospital. As a result, when planning an effective response to support patients with COVID-19, clinicians and health administrators should consider the need for additional critical care resources in the ICU, as well as the wider hospital.

Abbreviations

ACCI:Age-adjusted Charlson Comorbidity Index
AHP:Allied health professions
AKI:Acute kidney injury
ARDS:Adult respiratory distress syndrome
AST:Aspartate aminotransferase
AVPU:“Alert, verbal, pain, unresponsive” scale
BMI:Body mass index
CKD:Chronic kidney disease
COVID-19:Coronavirus disease-19
CPAP:Continuous positive airway pressure
CRP:C-reactive protein
CT:Computed tomography
DNR:Do-not-resuscitate
ECMO:Extracorporeal membrane oxygenation
EHR:Electronic health record
FiO2:Fraction of inspired oxygen
GGO:Ground-glass opacities
GGT:Gamma-glutamyl transferase
HR:Hazard ratio
HIV:Human immunodeficiency virus
ICU:Intensive Care Unit
IMD:Index of Multiple Deprivation
IQR:Interquartile range
LDH:Lactate dehydrogenase
LV:Left ventricle
MERS:Middle East respiratory syndrome
NEWS2:National Early Warning Score version 2
NIV:Noninvasive ventilation
PaO2:Arterial partial pressure of oxygen
PE:Pulmonary embolism
PFR:PaO2/FiO2 ratio
PPE:Personal protective equipment
RRT:Renal replacement therapy
RT-PCR:Reverse-transcriptase polymerase chain reaction
SARS-CoV-2:Severe acute respiratory syndrome coronavirus-2
SARS:Severe acute respiratory syndrome
SBP:Systolic blood pressure
SOFA:Sequential Organ Failure Assessment
SOP:Standard operating procedure
TEP:Treatment escalation plan
UK:United Kingdom
USA:United States of America
WBC:White blood cells.

Data Availability

The dataset generated and analysed during the current study is not publicly available due to Institutional Review Board approval restrictions.

Disclosure

This paper has been published as a preprint and is available at the following site: https://www.medrxiv.org/content/10.1101/2020.07.16.20155069v1. This research did not receive any funding and was performed as part of the authors’ employment with King’s College Hospital NHS Foundation Trust.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

SV and RM conceived the study and its design, had full access to the data, and take responsibility for the integrity of the data and accuracy of the analysis. SV, AW, VM, SC, CLS, JP, KL, TP, CS, KA, BA, RB, FJ, SAH, WB, and RM extracted and entered the data. SV and RM contributed to data analyses. SV, AW, VM, SC, CLS, JP, KL, WB, and RM contributed to data interpretation. SV and RM drafted the manuscript. All authors critically revised the drafted manuscript and approved the submitted manuscript.

Acknowledgments

We thank Professor Alistair McGuire (Department of Health Policy, London School of Economics) for his valuable help in the preparation of the manuscript. We also thank King’s College Hospital Business Intelligence Unit for providing us with heatmap images of the pandemic evolution in the area. We thank “Easy Mobile Forms Software, The DataDyne Group, LLC” for providing 2 months of free access to their data collection tool Magpi. Finally, we would like to thank our patients and fellow health-care workers for providing outstanding patient care at considerable personal risk.

Supplementary Materials

Supplementary File 1. File name: Supplementary_file_1.pdf. Title: institutional setting. Description: details regarding the institutional setting of King’s College Hospital and pandemic surge in the catchment area, including its temporal evolution with the use of heat maps. Supplementary File 2. File name: Supplementary_file_2.pdf. Title: definition of variables. Description: details regarding the definitions of variables used in the study. Supplementary File 3. File name: Supplementary_file_3.pdf. Title: statistical methods. Description: details regarding the statistical approach including the selection of variables and models. Description of the parametric regression equation. Additional figure describes the effect of each variable on the primary outcome of hospital mortality. Supplementary File 4. File name: Supplementary_file_4.pdf. Title: institutional surge planning and resource implications. Description: details regarding the reconfiguration of services implemented during the pandemic surge and the resources required. (Supplementary Materials)