key: cord-0684298-2lk6w4va
authors: Dongelmans, Dave A.; Termorshuizen, Fabian; Brinkman, Sylvia; Bakhshi-Raiez, Ferishta; Arbous, M. Sesmu; de Lange, Dylan W.; van Bussel, Bas C. T.; de Keizer, Nicolette F.
title: Characteristics and outcome of COVID-19 patients admitted to the ICU: a nationwide cohort study on the comparison between the first and the consecutive upsurges of the second wave of the COVID-19 pandemic in the Netherlands
date: 2022-01-13
journal: Ann Intensive Care
DOI: 10.1186/s13613-021-00978-3
sha: 42d5936bd6be804e5ba95fee7a1b36e5dfb98b08
doc_id: 684298
cord_uid: 2lk6w4va

BACKGROUND: To assess trends in the quality of care for COVID-19 patients at the ICU over the course of time in the Netherlands. METHODS: Data from the National Intensive Care Evaluation (NICE)-registry of all COVID-19 patients admitted to an ICU in the Netherlands were used. Patient characteristics and indicators of quality of care during the first two upsurges (N = 4215: October 5, 2020–January 31, 2021) and the final upsurge of the second wave, called the ‘third wave’ (N = 4602: February 1, 2021–June 30, 2021) were compared with those during the first wave (N = 2733, February–May 24, 2020). RESULTS: During the second and third wave, there were less patients treated with mechanical ventilation (58.1 and 58.2%) and vasoactive drugs (48.0 and 44.7%) compared to the first wave (79.1% and 67.2%, respectively). The occupancy rates as fraction of occupancy in 2019 (1.68 and 1.55 vs. 1.83), the numbers of ICU relocations (23.8 and 27.6 vs. 32.3%) and the mean length of stay at the ICU (HRs of ICU discharge = 1.26 and 1.42) were lower during the second and third wave. No difference in adjusted hospital mortality between the second wave and the first wave was found, whereas the mortality during the third wave was considerably lower (OR = 0.80, 95% CI [0.71–0.90]). CONCLUSIONS: These data show favorable shifts in the treatment of COVID-19 patients at the ICU over time. The adjusted mortality decreased in the third wave. The high ICU occupancy rate early in the pandemic does probably not explain the high mortality associated with COVID-19. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13613-021-00978-3.

Various studies reported high mortality rates among COVID-19 patients admitted to the Intensive Care Unit (ICU). These rates appear to be higher than reported Open Access *Correspondence: F.Termorshuizen@amsterdamumc.nl 1 National Intensive Care Evaluation (NICE) Foundation, PO Box 23640, 1100 EC Amsterdam, The Netherlands Full list of author information is available at the end of the article Dongelmans et al. Annals of Intensive Care (2022) 12:5 among ICU patients with other types of viral pneumonia [1] . During the COVID-19 pandemic in the Netherlands, many ICUs were scaled up above their maximum capacity. Of note, the Netherlands healthcare system has less ICU beds per capita than other European countries, i.e. 6.4 per 100,000. For instance, Belgium with 15.9 or Germany with 29.2 have significantly more ICU beds per 100,000 inhabitants [2] . This combined with both the high number of patients and their longer length of stay (when compared to average ICU patients) caused a limitation in the numbers of available ICU beds. This potentially has had adverse effects on the quality of care, and may have increased the mortality risk. To utilize ICU beds most effectively, preventing overload due to new cases of high emergency (related or not related to , ICU patients were transferred between regional and (inter)national ICUs [3] . A relocation is not risk free, especially for mechanically ventilated patients [4] .

In the present study, we describe the clinical characteristics, the length of ICU stay, relocation rate, and mortality of ICU patients infected with SARS-CoV-2 during the first wave, the period in-between and the consecutive upsurges thereafter of the COVID-19 pandemic in the Netherlands. The aim is to examine whether there is evidence for an improvement in the quality of care for COVID-19 patients at the ICU, as indicated by a shorter duration of stay, a lower relocation rate, and a lower hospital mortality rate during the epidemic upsurges after the first wave. This will also be helpful to understand whether the high mortality risk associated with COVID-19 is an inevitable effect of the virus infection itself or was also due to the (imminent) inadequate response of the health care system in the early phases of the pandemic.

The National Intensive Care Evaluation (NICE) registry is a voluntary quality registry in which all Dutch ICUs participate. This registry includes prospectively collected clinical data of all patients admitted to an ICU and includes demographics, major comorbidities, physiology, clinical course, ICU length of stay, in-hospital mortality and mechanical ventilation. Of note, mechanical ventilation in our data set includes invasive as well as non-invasive ventilation, there is no differentiation between the two. High flow nasal oxygen (HFNOT) is excluded as a mechanical ventilation mode. These routinely collected data are extracted from the electronic health record and after validation uploaded on a monthly basis. According to clear definitions and using rigorous data quality checks maximum data quality was ensured [5] . The purpose of NICE is to provide feedback on performance indicators to ICUs, thus enabling ICUs to monitor and improve their quality of care. From the start of the COVID-19 outbreak in the Netherlands, the Dutch government requested all ICUs to record all suspected and confirmed COVID-19 patients admitted to the ICU. Therefore, the existing NICE data infrastructure was expanded with a module allowing for daily recording of admission-and discharge dates, and survival status at ICU-and hospital discharge of COVID-19 patients. Clinical data of the COVID-19 patients were linked afterwards when uploaded according to the regular NICE processes. NICE data infrastructure allowed that COVID-19 patients could be accurately tracked throughout subsequent hospital admissions. This made it possible to combine the length of ICU stay between hospitals and final hospital survival status, thus taking transfers into account.

A confirmed COVID-19 patient was defined as follows: a positive SARS-CoV-2 Reverse transcription polymerase chain reaction (RT-PCR) on a nasopharyngeal swab or a CT-scan consistent with COVID-19 (i.e., a CO-RADS score of ≥ 4 in combination with the clinical diagnosis viral pneumonia) [6] .

In 2020 and 2021, two big waves of ICU admissions due to COVID-19 infections were observed in the Netherlands with a distinct period in between ( Fig. 1 and Additional file 1: Fig. S1 ). We decided to categorize calendar time into the following time periods under investigation: the first wave (February 1-May 24, 2020), the episode between the waves (May 25-October 4, 2020), the second wave-first and second upsurge (October 5, 2020-January 31, 2021), and the second wave-final upsurge (February 1, 2021-June 30, 2021). Albeit the final upsurge of the second wave does not appear to be a separated wave, that is, did not follow after a substantial decline in cases towards zero, for ease of reference we call the final upsurge of the second wave the 'third wave' . The period of this third wave is chosen based on visual inspection and because of the fact that the vaccination program in the Netherlands commenced on January 6th 2021 with a focus on vulnerable and older people at first. About half a million people had received their first vaccination on February the first. From then on vaccination was sped up and about a million people received a vaccination per month. Patients categorized in a time period under investigation were analyzed until death or hospital discharge, that is, beyond the time period's limits. Thus, a death event in period two of a patient admitted in period one was assigned to period one.

The Medical Ethics Review Committee of the Academic Medical Center waived the need for informed consent [reference number W21_091 # 21.102].

The number of patients, patient-and treatment characteristics, duration of total ICU-and hospital stay (combining subsequent hospital and ICU stays within a treatment trajectory of an unique patient), occupancy rate calculated as the number of occupied beds at the day of ICU admission of the patient at the pertinent ICU as percentage of the daily average in 2019 specific for that ICU, and relocation rate were described for the four periods separately. The in-hospital mortality, that is death at the ICU or death at the hospital ward after ICU discharge, was estimated as percentage and analyzed in a logistic regression model with wave (first, period inbetween, second and third) as the main covariate. The Odds ratios (ORs) of hospital mortality were adjusted for age, sex, BMI, the APACHE-IV mortality probability, and the ICU occupancy rate at the day of ICU admission. The length of ICU stay was analyzed in a multivariable Cox regression model with ICU discharge as outcome event. Again, wave was the main variable of interest and the Hazard ratio (HR) of ICU discharge was adjusted for the above-mentioned co-variates.

A higher or lower mortality rate may lead to a shorter or longer mean length of ICU stay, respectively. To examine whether shifts in ICU length of stay are independent of possible shifts in mortality, the analysis was performed in two different ways: first with time of death included as an ICU discharge event, and second with time of death included as reason for censoring. The last analysis is probably more valid, as it reduces the influence of short stays caused by early death events, and thus reflects trends in length of stay assuming that all patients survived the ICU admission. If no substantial differences between both methods are found, this indicates that shifts in mortality did not influence shifts in mean length of stay.

All analyses were performed using the R statistical environment (version 3.6.1) (R Foundation for Statistical Computing, Vienna, Austria). We report P values and effect estimates with 95% confidence intervals.

From February 2020 until the first of July 2021 there were 12,723 COVID-19 patients admitted to 78 Dutch ICUs. Of these, we excluded 5.4%, since their clinical data were not available, leaving 12,030 patients (94.6%) for the final analyses. Of note 2.6% of the patients were not discharged from the hospital at the time of the final analysis.

In Table 1 , the demographic and clinical characteristics of included patients are given, for the full cohort, and stratified by period (i.e. Wave 1, Period in-between, Wave 2, and Wave 3). During Wave 1 compared to Wave 2, the proportion of patients aged 75 years or older was 13.6 vs. 17.0% (Additional file 1: Table S1 ). In Wave 3 the proportion of this group decreased to 10.7%. The In accordance with this finding, the proportion of patients with 2 or more comorbidities was higher during Wave 2 compared to Wave 1 (12.9% vs. 7.5%) and decreased in Wave 3 (9.0%). The prevalence figures for each of the comorbidities were higher during Wave 2, compared to either Wave 1 or the period in-between, whereas in Wave 3 these decreased again (Additional file 1: Table S1 ). In particular, renal insufficiently increased from 2.7% in Wave 1 to 5.9% in Wave 2 and decreased to 3.4 in Wave 3. The APACHE-IV estimated risk of in-hospital mortality remained relatively stable with a slight increase in Wave 2 (Table 1) . Over time there was an increase in the percentage of obese patients BMI (> 30 kg/m2).

As presented in Table 2 , the mean (sd) time spend in hospital before admission to the ICU increased almost a day, from 1.7 (3.0) days in the Wave 1 to 2.4 (9.6) days in the Wave 3. The percentage of patients who were mechanically ventilated at the moment of admission decreased between Wave 1 from 48.1% to 25.1% in Wave 2 and this remained in the same range in Wave 3 (24.7%). The percentage of patients who were mechanically ventilated in the first 24 h of admission was 79.1% vs. 58.1% for Wave 1 and Wave 2, respectively, and this remained stable in Wave 3 (58.2%). The use of vasoactive drugs decreased as well (67.2% vs. 48.0%) and in Wave 3 this was even lower (44.7%). The ICU occupancy rates were on average lower both during Wave 3 and Wave 2 and the period in-between than the occupancy rate during Wave 1 (168%, 155% and 109% vs. 183%, respectively). In accordance with this finding, the relocation rates during the last periods were lower than the relocation rate during Wave 1 (23.8%, 27.6%, and 22.7% vs. 32.3%, respectively). 

The crude percentage of patients who died in the hospital during Wave 1 was 29.9%. This crude hospital mortality decreased to 27.1% in the in-between period, increased to 32.0% during Wave 2 and then decreased to 23.4% in Wave 3 (Table 2 ). After adjustment for age, gender, BMI, and the APACHE-IV mortality probability, the odds of hospital death in Wave 3 was significantly lower than the odds of hospital death in Wave 1: OR = 0.79 [0.70-0.89]. Further adjustment for the ICU occupancy rate at the day of ICU admission did not modify these results: OR = 0.80 [0.71-0.90] (Fig. 2 , Additional file 1: Table S2 ).

The average length of ICU stay was 20.6 days during Wave 1 and this decreased to 16.3 days during the in-between period and to 17.2 days during Wave 2 and decreased even further to 16.0 days in Wave 3 ( Table 2 ). The significantly shorter duration until ICU discharge during Wave 2 and Wave 3 compared to Wave 1 was also shown in a Cox regression model with death as censoring event and after adjustment for age, gender, BMI, and the APACHE-IV mortality probability (HR = 1.31 [1.23-1.38] and 1.45 [1.37-1.54]) (Fig. 3 , Additional file 1: Table S3 ). Further adjustment for the ICU occupancy rate did not substantially modify this result. Similar results were found when death was included as discharge event instead of censoring event (Additional file 1: Table S3 ).

In this study, we compared the second wave, split into two parts, with the first wave of ICU COVID-19 admissions in the Netherlands. We found a shift to patients with a less favorable prognosis at ICU admission during the first and second upsurge of Wave 2, i.e. towards older patients with more comorbidities and a higher APACHE-IV mortality probability. In addition, we found a shift to favorable outcome and treatment characteristics, that is a shorter mean length of stay at the ICU, less relocations to other hospitals, and a lower number of patients treated with mechanical ventilation and vasopressor drugs. Furthermore, during the final upsurge of Wave 2 ('Wave 3') Table S3 , death as censoring event). A Hazard ratio of ICU discharge higher than 1.00 implies a comparatively high rate of discharge and, thus, a shorter length of stay at the ICU. Wave 1, February, 2020-May 24, 2020. Period in-between, May 25, 2020-October 4, 2020. Wave 2: first and second upsurge of Wave 2, October 5, 2020-January 31, 2021. Wave 3: final upsurge of Wave 2, February 1-June 30, 2021 the age declined. It is unfortunate that no improvement in the patient survival could be established up until Wave 3. However, we cannot exclude that selection of patients with the best chance of ICU survival in Wave 1, due to health care driven scarce in ICU beds in the Netherlands, might have led to underestimation of the differences in survival.

We found that case-mix adjusted hospital mortality of ICU patients decreased in Wave 3 of the pandemic in the Netherlands. This suggests that insights in the pathophysiology of COVID-19 with increasingly appropriate treatments and vaccination and improved national healthcare system response by effective logistic-and organizational preparations (such as patient spreading) entailed more efficient care for COVID-19 patients.

Comparing our cohort of COVID-19 patients with respect to age to other cohorts of ICU patients shows that our patients are comparable to patients in other countries described in a review by Serafin et al. [7] . An important observation is that the age of COVID-19 patients in our study was higher during Wave 2. This was also observed in a study comparing the first and second wave in France [8] . In a large cohort of ICU patients in Brazil pre-vaccination, the age of patients admitted to the ICU declined over time [9] . In our view, the age of patients admitted to the ICU is influenced by selection of patients on one hand and the vaccination program on the other. For instance, in the Netherlands general practitioners played an important role in making it possible for patients to stay at home and receive tailored (palliative) care [10] . Scarcity of resources, i.e. availability of ICU beds, could also influence patient selection. The ICU occupancy rate was, as we showed, on average much higher during Wave 1 than during the in-between period and Wave 2 and Wave 3 ( Table 2) . Whether or not this truly has played a role in the selection of patients remains an important question. During Wave 3 the vaccination program started targeting vulnerable and older subjects at first.

In our cohort the number of comorbidities per patient was higher in Wave 2 and declined in Wave 3. In the cohort described by ICNARC from the UK there was also a slight increase in comorbidities [11] , in contrast with the French cohort, where there was a non-significant decrease of comorbidities. These trends in number of comorbidities could again be explained by changing selection criteria over time. In our cohort we feel that a more liberal attitude towards admission of patients with comorbidities was the case in Wave 2. The decline in comorbidities in Wave 3 may be associated with the decrease in age due to the vaccination program. Moreover, the differences between countries could also be explained by the number of available ICU beds per country [2] .

With respect to the need for mechanical ventilation, our study showed a decrease between Wave 1 and Waves 2 and 3. The same pattern was observed in Germany, where the percentage of ICU patients on mechanical ventilation decreased from 64 to 54% [12] . In France the need for invasive ventilation (this does thus not include non-invasive mechanical ventilation) in ICU dropped from 88 to 64% [8] . The reason for this decrease in our cohort is clearly not a decline in severity of illness. It could, however, be due to the fact that alternative means of providing oxygen support were used more often in ICU patients. Although recognized early as a possible respiratory therapy in hypoxemic failure, there was reluctance to use of HFNO due to presumed risk of transmission of the virus in the aerosols produced while using HFNO at the start of the pandemic [13] . Unfortunately, the use of HFNO is not in our database to further analyze this. The variation in the use of mechanical ventilation is seen in cohorts worldwide [14] .

The length of ICU stay of patients was lower in Waves 2 and 3. Since patients in our cohort showed a stable disease severity and an increase in the number of comorbidities, the opposite, i.e. an increase of length of stay, would have been expected in Wave 2. In our view the decrease in length of ICU stay could have had several reasons. The first possibility is a change in the treatment for the most characteristic symptom of COVID-19 patients, hypoxia. As we showed, there was a decrease in use of mechanical ventilation. On the other hand, since invasive mechanical ventilation has inherent negative effects and may cause damage to lung, it could be that COVID-19 patients benefited from this change [15] .

The second possibility is that because commencing mechanical ventilation requires the use of muscle relaxants and sedatives and use of the latter is known to lead to an increase in the length of stay [16] . Both mechanical ventilation and sedation are only applied when necessary and are lifesaving treatments. Yet another possibility would be that the ubiquitous use of dexamethasone led to a change in disease pattern, i.e. patients were further in their disease process when admitted to the ICU [17] . This is supported by the fact that patients stayed longer in the ward before being admitted to the ICU. Our results show that between Wave 1 and the Wave 3 the patients stayed on average almost a day longer on the ward. This may be explained by increased testing capacity over time leading to earlier diagnosis within the disease trajectory of individual patients affecting earlier hospital admission. Another possibility might have been that extra resources in the Wave 1 were more difficult to organize than in Waves 2 and 3. Not only bed utilization rates were higher but also extra personnel was not yet sufficiently trained during Wave 1, which might have affected quality of care [18] .

The crude in-hospital mortality of COVID-19 ICU patients in our study remained at the same level of around 30% and declined in Wave 3 to below 25%. The fact that medication such as corticosteroids and tociluzimab had become standard of care may explain a part of this reduction. Second and more importantly, the vaccination programme commenced focusing on old and or fragile groups first. After adjustment for case-mix and other confounders this difference in mortality between the different periods was still apparent. Some studies show markedly different mortality levels, for instance the French cohort showed an ICU mortality of 50% [8] . In a Brazilian cohort, a striking different course of mortality over time was shown, in-hospital mortality decreased from a low 18.0% to an even lower 9.8% [9] . However, in a systematic review describing mainly the first wave and including 32 papers, the mean ICU mortality was comparable with our study, i.e., 30% [7] .

A major strength of our study was that we were able to report on the first and second pandemic waves in the Netherlands covering all COVID-19 patients admitted to the ICU nationwide, serving a population of around 17.5 million people. We used an existing quality registry which for the COVID-19 pandemic was serving as one of the nation's ways to track and model the pandemic in the Netherlands. Notably, transferred patients' stays in multiple ICUs were combined, thus including each patient only once in our cohort, which is a major advantage over other cohorts that did not combine multiple stays at different ICUs. A limitation of our approach is that we did not have information on treatment protocols applied to the patients in our cohort. In the Netherlands, the use of dexamethasone and toziluzimab, as well as anticoagulant prophylaxis and other treatment, in these patients was directed by the scientific communities and institutions [19] . We also have no information on the use of noninvasive mechanical ventilation in relation to the use of invasive mechanical ventilation. We only have data on the use of mechanical ventilation as a whole. In addition, important reported comorbidities such as hypertension and hypercholesterolemia were not available in our data set [20] . Yet, another limitation is that we do not have data on the way COVID-19 was confirmed per patient. The fact is that in the second wave there was sufficient capacity to confirm the diagnosis with a PCR, whereas in the first wave, we had to rely more on radiologic confirmation. The change in confirmation method might have slightly affected the number of false positives and/ or negatives. However, we do not expect that this has influenced the observed mortality rate in COVID-19 patients in our cohort. Since COVID-19 which leads to an ICU admission is a clinical diagnosis, the confirmation is done as one needs to be as certain of the diagnosis as possible. We are confident that the vast majority of patients was correctly diagnosed with COVID-19 notwithstanding the confirmation method.

In this study, the subsequent pandemic waves of ICU admissions due to SARS-CoV-2 infection in the Netherlands were compared. We found a shift to patients who were slightly older and had more comorbidities in the first two upsurges of Wave 2 and a decline in age and number of comorbidities in the final upsurge of Wave 2 ('Wave 3'). In Waves 2 and 3 a favorable shift in treatment characteristics with shorter ICU stay, less relocations, and less mechanical ventilation was seen. This indicates more efficient care for COVID-19 patients admitted to the ICU. Moreover, the high ICU occupancy rate early in the pandemic does probably not explain the high mortality associated with COVID-19. The adjusted in-hospital mortality decreased in Wave 3. 

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s13613-021-00978-3.

Additional file 1: Figure S1 . Mean occupancy rate at the ICU (as fraction of the average number of patients in 2019) per week during the COVID-19 waves. Table S1 . Patient characteristics. Table S2 . Crude and adjusted logistic regression showing Odds ratios of hospital death during Wave 2, Wave 3, and the period in-between compared to Wave 1 (see Fig. 2 ). Table S3 . Crude and adjusted Cox regression showing higher Hazard ratios of ICU discharge during Wave 2, Wave 3, and the period in-between compared to Wave 1 (see Fig. 3 ). A Hazard ratio of ICU discharge higher than 1.00 implies a comparatively high rate of discharge and, thus, a shorter length of stay at the ICU. Shown are an analysis with death as ICU discharge event and an analysis with death as censoring event.

ICNARC COVID-19 Team. COVID-19 in critical care: epidemiology of the first pandemic wave across England, Wales and Northern Ireland

The variability of critical care bed numbers in Europe

Transport of COVID-19 and other highly contagious patients by helicopter and fixed-wing air ambulance: a narrative review and experience of the Swiss air rescue Rega

ICU MOVE Investigators and the Pragmatic Critical Care Research Group. Effect of interhospital ICU relocation on patient physiology and clinical outcomes

Data Resource Profile: the Dutch National Intensive Care Evaluation (NICE) registry of admissions to adult intensive care units

COVID-19 Standardized Reporting Working Group of the Dutch Radiological Society. CO-RADS: a categorical CT assessment scheme for patients suspected of having COVID-19-definition and evaluation

Clinical course and outcomes of critically ill patients with COVID-19 infection: a systematic review

Comparison between first and second wave among critically ill COVID-19 patients admitted to a French ICU: no prognostic improvement during the second wave? Crit Care

Evolving changes in mortality of 13,301 critically ill adult patients with COVID-19 over 8 months

Intensieve en palliatieve covid-19-zorg door huisarts [Intensive covid-19 care in general practice; nationwide registration study

A new risk prediction model for critical care: the Intensive Care National Audit & Research Centre (ICNARC) model

Case characteristics, resource use, and outcomes of 10 021 patients with COVID-19 admitted to 920 German hospitals: an observational study

High flow nasal oxygen therapy to avoid invasive mechanical ventilation in SARS-CoV-2 pneumonia: a retrospective study

International variation in the management of severe COVID-19 patients

Ventilator-induced lung injury

Sedation in the intensive care unit: a systematic review

Corticosteroid use in COVID-19 patients a systematic review and metaanalysis on clinical outcomes

Organizational characteristics: effect on outcome of ICU COVID-19 patients in Belgium-authors' reply

Incidence of thrombotic complications and overall survival in hospitalized patients with COVID-19 in the second and first wave

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

We would like to acknowledge the effort of all Dutch hospitals to record their COVID-19 patients three times a day and additionally all the ICUs that also collected the comprehensive clinical data of the COVID-19 patients to NICE-DB. Without their effort this study would not have been possible. We also like to thank all employees of NICE Research & Support for the technical development of the COVID-19 Online Module, processing the data, and support under high work pressure. 

Authors' contributions DD: contributed to the acquisition of data, the conception and design of the study, drafted and wrote the manuscript, and contributed to the interpretation of the results; FT: contributed to the design of the study, performed the statistical analysis, wrote parts of the manuscript and contributed to the interpretation of the results, and approved the version to be submitted; SB: contributed to the design of the study, the performing of the statistical analysis and the interpretation of the results, and approved the version to be submitted; FB-R: contributed to the design of the study, the interpretation of the results and the writing of the manuscript, and approved the version to be submitted; SA: contributed to the acquisition of data, the conception of the study, the interpretation of the results and the writing of the manuscript, and approved the version to be submitted; DdL: contributed to the acquisition of data, the conception of the study, the interpretation of the results and the writing of the manuscript, and approved the version to be submitted; BvB: contributed to the acquisition of data, the conception of the study, the interpretation of the results and the writing of the manuscript, and approved the version to be submitted; NK: contributed to the conception and design of the study, the interpretation of the results and the writing of the manuscript, and approved the version to be submitted. The Collaborators of the Dutch COVID-19 Research Consortium contributed to the acquisition of data, read and reviewed the manuscript, and approved the version to be submitted. All authors have read and approved the final manuscript. 

De-identified aggregated data used for the study will be made available upon reasonable request.

The Medical Ethics Review Committee of the Academic Medical Center had no objection to this study and waived the need for informed consent [reference number W21_091 # 21.102].

Not applicable.

The authors declare that they have no competing interests.