Skip to main content

The direct disease burden of COVID-19 in Belgium in 2020 and 2021

Abstract

Background

Burden of disease estimates have become important population health metrics over the past decade to measure losses in health. In Belgium, the disease burden caused by COVID-19 has not yet been estimated, although COVID-19 has emerged as one of the most important diseases. Therefore, the current study aims to estimate the direct COVID-19 burden in Belgium, observed despite policy interventions, during 2020 and 2021, and compare it to the burden from other causes.

Methods

Disability-adjusted life years (DALYs) are the sum of Years Lived with Disability (YLDs) and Years of Life Lost (YLLs) due to disease. DALYs allow comparing the burden of disease between countries, diseases, and over time. We used the European Burden of Disease Network consensus disease model for COVID-19 to estimate DALYs related to COVID-19. Estimates of person-years for (a) acute non-fatal disease states were calculated from a compartmental model, using Belgian seroprevalence, social contact, hospital, and intensive care admission data, (b) deaths were sourced from the national COVID-19 mortality surveillance, and (c) chronic post-acute disease states were derived from a Belgian cohort study.

Results

In 2020, the total number of COVID-19 related DALYs was estimated at 253,577 [252,541 − 254,739], which is higher than in 2021, when it was 139,281 [136,704 − 142,306]. The observed COVID-19 burden was largely borne by the elderly, and over 90% of the burden was attributable to premature mortality (i.e., YLLs). In younger people, morbidity (i.e., YLD) contributed relatively more to the DALYs, especially in 2021, when vaccination was rolled out. Morbidity was mainly attributable to long-lasting post-acute symptoms.

Conclusion

COVID-19 had a substantial impact on population health in Belgium, especially in 2020, when COVID-19 would have been the main cause of disease burden if all other causes had maintained their 2019 level.

Peer Review reports

Introduction

On the 11th of March 2020, the World Health Organization (WHO) declared the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak as a pandemic [1]. In Belgium, the first infections with SARS-CoV-2 were reported by the end of February 2020. A rapid increase in cases resulted in a lockdown of the country on the 18th of March 2020 and the peak of the first wave was reached at the beginning of April 2020 [2]. After a national lockdown, virus circulation increased again, causing a high number of infections, hospitalizations, and casualties over different waves in 2020 and 2021 [3]. To prevent the further transmission of SARS-CoV-2 in 2020 and 2021, people were encouraged to maintain protective public health measures such as physical distancing, extensive testing and tracing followed by isolation of cases and quarantine of close contacts, wearing facemasks in public, and from 2021 onwards when vaccines became available, COVID-19 vaccination.

Monitoring population health is essential to identify where health improvements are required, to develop targeted prevention and intervention strategies, and to guide policies for addressing unmet health needs and identifying health priorities. To this end, burden of disease assessments allow estimates of the occurrence of morbidity and mortality to be translated into a single measure, the Disability-Adjusted Life Years (DALYs), thereby facilitating comparisons with other causes of disease and injury [4]. This composite DALY measure is calculated by standardizing the effects of morbidity and mortality as health loss over time [5]. For the morbidity component, severity distributions and disability weights (DWs) indicate how burdensome living with a disease is; whereas age-conditional life expectancies reveal the disease burden caused by premature mortality, implying that death at a younger age constitutes a greater burden, ceteris paribus [6]. The direct impact of COVID-19 on individual health occurs via three main pathways, (1) the acute infectious phase with a morbidity component (e.g., feeling sick, being hospitalized, etc.), (2) the acute infectious phase with a mortality component (i.e., COVID-19 related fatalities), and (3) a chronic phase with long-lasting morbidity beyond the period of the acute symptomatic infection. During the acute symptomatic stage of the infection, COVID-19 can cause a wide spectrum of symptoms, such as nausea, coughing, and loss of taste and smell that could reduce the quality of life (expressed by DWs between 0 and 1, with 0 reflecting “perfect health” and 1 reflecting “equivalent to death”), and, in the most extreme case, results in death [7]. If symptoms persist beyond the acute infectious period, i.e., post-COVID-19 condition, symptoms could contribute to a chronic morbidity component, whereby some symptoms persist for longer periods [8]. Additionally, there is mounting evidence that COVID-19 leads to long-term elevated risks of relatively common causes of disease burden, such as cardiovascular and neurological complications, including ischemic and hemorrhagic stroke, cognition and memory disorders, peripheral nervous system disorders, episodic disorders, extrapyramidal and movement disorders, mental health disorders, musculoskeletal disorders, sensory disorders, Guillain–Barré syndrome, and encephalitis or encephalopathy [9, 10]. However, quantifying the indirect disease burden resulting from COVID-19 is a challenge due to its complex and multi-faceted nature. The estimation of DALYs requires the existence of and accessibility to a wide variety of data sources to match the different disease states of the underlying YLLs and YLDs model (e.g., cases, hospitalizations, and deaths).

Clearly, pandemics, such as COVID-19, also exert indirect impacts beyond COVID-19 patients and their personal family, friends and employers [11]. For example, perceived risks of acquiring COVID-19 and imposed non-pharmaceutical interventions influenced the frequency and severity of mental health issues [12]. Furthermore, both risk perceptions and capacity problems due to COVID-19, caused postponements and cancellations of unrelated preventive and curative healthcare [13,14,15]. Even if these indirect impacts are beyond the scope of the current study, it represents an essential first step, by attempting to estimate the extent of the direct impact on population health in Belgium, while focusing on the first two calendar years of the pandemic (2020 and 2021).

This study aimed to estimate the COVID-19 burden including direct acute and direct post-acute consequences using the DALY approach at the national level in Belgium anno 2020 and 2021 by relying on a consensus burden of disease model, developed for COVID-19 by the European Burden of Disease Network [16].

Methods

The current study has adopted the European Burden of Disease Network consensus disease model for COVID-19 using a prevalence-based approach to estimate COVID-19 related DALYs [16]. DALYs reflect the healthy life years lost due to diseases and constitute a morbidity component, i.e., Years Lived with Disability (YLDs), and a mortality component, i.e., Years of Life Lost (YLLs). Estimates of the burden of disease were calculated at the Belgian national level. Direct inputs for the calculation of the disease burden were (1) the output from an age-structured stochastic compartmental model, which has been calibrated on age-specific seroprevalence data, confirmed COVID-19 hospitalizations and deaths among others, while accounting for changing contact patterns, vaccination uptake and different SARS-CoV-2 variants of concern [17], (2) a cohort in the Belgian population called COVIMPACT including 6,913 participants who reported to experience COVID-19 symptoms after a positive test, for which at least two follow-up points were available (December 2022) to assess the proportion of acute COVID-19 patients that develop post-acute long-lasting symptoms [18], and (3) COVID-19 mortality surveillance systems in Belgium. The total disease burden in DALYs is calculated as the sum of the morbidity in YLDs and mortality in YLLs components.

Morbidity (Years Lived with Disability)

Data sources

YLDs were estimated using the person-time estimates from the stochastic compartmental model developed for Belgium. This stochastic transmission model as described by Abrams et al. (2021) has previously been adapted to include vaccination and the emergence of several variants of concern. [17] The model is calibrated on early seroprevalence data, genomic surveillance data, hospital admission data, mortality data, and social contact data from the Belgian CoMix survey [19]. Gradually accumulating and waning of naturally-acquired immunity in the population is accounted for, as well as immunity induced by vaccination. The vaccine-type and age-specific uptake over time is based on the reported data by Sciensano, derived from Epistat. The model results contain stochastic variation in the transmission process and parameter uncertainty based on 40 model parameter configurations. The calibration process relies on likelihood-based MCMC sampling until convergence of a single MCMC chain after which 40 posterior samples of the joint parameter distribution are generated. Each of these posterior samples is then used to simulate 10 stochastic realizations of the model output, leading to a total of 400 simulated values for, for example, the total number of (a)symptomatic infections per day between the 1st of March 2020 (the beginning of the COVID-19 crisis) and the 31st of December 2021 stratified by 10-year age groups. Symptomatic infections were characterized in three compartments: (1) not requiring hospitalization; (2) requiring hospitalization but not intensive care unit admission; and, (3) requiring intensive care unit admission. Based on these estimates, we calculated the total number of person-years by compartment (i.e., health state), whereby the number of people that were in each of the aforementioned compartments and the time spent in the compartment (i.e., duration of illness) was considered. As the proportion of people that continue to experience symptoms after the acute symptomatic phase is currently unknown and debated [8, 20, 21], we have estimated the proportion of patients that develop post-acute consequences and the duration of these post-acute consequences based on the longitudinal Belgian COVIMPACT cohort study, which, as described by Smith et al. (2022), was designed for that purpose [18, 22]. This population-based cohort included a total of 6,913 participants in December 2022 with an age between 18 and 102 years. 35.1% of participants reported to be men, 64.8% to be women, and 0.1% described themselves as the category other. The COVIMPACT study assesses self-reported symptoms at the time of infection, and at regular 3-month intervals after infection between the 29th of April 2021 and the 31st of December 2022. In the calculations, YLDs caused by COVID-19 were only considered during the respective calendar years, i.e., the prevalence approach [23].

Health states

The different health states as defined by the European Burden of Disease Network consensus model [16] that were considered in this study are listed in Table 1. A total of five different health states are described, of which the health states “asymptomatic”, “mild/moderate”, “severe”, and “critical” were approximated by the output from the stochastic compartmental model for Belgium. The health state “post-acute consequences” was approximated based on the results from the COVIMPACT study, where information on “tiredness” and “pain” was used to align the study outcomes with the definition of the “post-acute consequences” health state. In total 180 of 5,866 participants (3.1%, 95% uncertainty interval: 2.6 − 3.6%) reported being “tired” and “feeling pain over the body”, as specified in the health state definition, three months after confirmation of the infection. The distribution of the duration of symptoms was inferred by considering different parametric distributions for which the parameters were estimated while accounting for left-truncation and the interval-censored nature of the data. More specifically, a truncated log-normal distribution with mean and standard deviation parameters µ and σ equal to 5.10 (95% CI: 5.07–5.13) and 0.26 (95% CI: 0.24–0.27), respectively, was selected. Considering a minimum duration of 3 months (90 days), the average duration of post-acute symptoms in patients that developed post-acute consequences within this cohort equalled 260 days with a standard deviation of 44.9 days. More detailed information on the applied methodology to model the “post-acute consequences” health state can be found in Appendix 1. The total number of person-years in this health state was calculated by applying a microsimulation decision tree model with 400 iterations, whereby each case that would transit from a “symptomatic” compartment (i.e., community-based cases, cases requiring hospitalization, and cases requiring intensive care unit admission) to the “recovered” compartment in the compartmental model had a probability of 3.1% (95% CI: 2.6 − 3.6%) to develop long-term symptoms of “tiredness” and “bodily pain” with a duration generated from the aforementioned truncated log-normal distribution.

Table 1 Health states and disability weights for estimating the COVID-19 burden in Belgium

Disability weights

DWs vary between 0 (no impact, perfect health) to 1 (equivalent to death) and reflect how incapacitating a health state is for a panel. For the different health states considered in this study, DWs were based on the Global Burden of Disease (GBD) 2019 study for infectious diseases (i.e., lower respiratory tract infections) [24]. For critical care, the DW was based on the European Disability Weights Study (Table 1) [25].

For each age category, YLDs were calculated by summing across all health states i, the product of the number of person-years (Ni), and the corresponding disability weight (DWi) [16]:

$$YLDs=\sum _{i}{N}_{i}{DW}_{i}$$

The total YLDs are calculated by summing age-specific YLDs.

Mortality (Years of Life Lost)

YLLs were calculated using an extracted dataset (extracted in December 2022) of the epidemiological mortality surveillance of COVID-19 in Belgium [26, 27]. This database is managed by Sciensano, the Belgian institute for health. Sciensano adopted an inclusive definition for COVID-19 death notifications, by also including possible COVID-19 related deaths, whereby patients met the following criteria: “at least one of the following symptoms that appear with no other obvious cause: cough, dyspnoea, thoracic pain, acute anosmia or dysgeusia” or “at least two of the following symptoms with no other obvious cause: fever, muscle pain, fatigue, rhinitis, sore throat, headache, anorexia, watery diarrhoea, acute confusion, sudden fall” or “deterioration if the patient shows chronic respiratory symptoms”. This resulted in a robust correlation between COVID-19 and all-cause mortality [2, 26]. A death was classified as a COVID-19 related death if it occurred in “a confirmed, radiologically-confirmed or possible case that occurred in any setting (hospitals, at home, nursing homes, and other settings) unless a clear alternative cause of death unrelated to COVID-19 is identified”. The database includes all COVID-related deaths during the first two calendar years of the pandemic (2020–2021). More information on the COVID-19-related mortality has been published elsewhere (see Renard et al., 2021 [26]). Missing data on age and sex were imputed by its marginal or conditional probability as described in Appendix 2.

For each age group, YLLs were estimated by multiplying the number of deaths due to COVID-19, denoted by D, by the specific remaining life expectancy RLE from the GBD reference life table 2019 [28]:

$$YLLs=D\times RLE$$

The total YLLs because of COVID-19 was obtained by summing all age-specific YLL estimates. The calculated YLLs in each respective calendar year considered all expected years of life lost (also beyond the respective calendar year) [23]. Future life years were not discounted in this analysis.

Burden of disease (DALYs)

The total disease burden of COVID-19 was estimated in absolute DALYs, in DALYs per 100,000 inhabitants, and in age-standardized DALYs per 100,000 individuals using the Belgian age distribution obtained from the Belgian statistical office (Statbel), which reflects the Belgian population on the 1stof January 2020 and 2021 [29]. After calculating the DALYs, we compared the estimates for 2020 and 2021 with the latest GBD DALY results for the top ten health conditions ranked at GBD level 3 for Belgium in 2019. We applied no discounting on the DALYs, YLLs, or YLDs components in the estimation.

Model uncertainty

For the calculations of the YLLs we relied on the Belgian COVID-19 mortality register, whereby we did not redistribute causes of deaths. To account for the uncertainty in the remaining life expectancy tables, we based the YLLs calculation on the GBD 2019 life table, and performed a sensitivity analysis by calculating the YLLs based on the Belgian national life table. However, the GBD 2019 life table and Belgian national lifetable refer to different constructs. More specifically, the GBD life table was constructed based on the lowest observed age-specific mortality rates by location and sex across all estimation years from all locations with populations over 5 million in 2016. The Belgian national life tables on the other hand indicate how a generation would die out throughout its existence, if it was exposed to the same mortality conditions per age as in the specified year. The results of this sensitivity analysis can be found in Appendix 3. For the YLDs, the model applied a Monte Carlo approach including 1000 replications to account for the uncertainty with respect to (1) the disease model (i.e., the DWs in the different health states as described in Appendix 4), (2) the estimates of the health-state specific person-years (i.e., based on the stochastic compartmental model) [30], (3) the estimated probability of having long-term symptoms, and (4) the estimated distribution of the duration of long-term symptoms. In the first step, the 400 iterations from the stochastic compartmental model and the 400 iterations from the microsimulation model for long-COVID symptoms were resampled to obtain 1000 iterations. In a second step, the uncertainty of the different disability weights was approximated by a Beta distribution that matches the mean and 95% uncertainty interval (betaExpert-function as available from the “prevalence” package in R [31, 32]) yielding a random sample of 1000 DW estimates per health state. In a third step, the uncertainty for post-acute symptoms was modelled considering the uncertainty around the proportion of patients that would develop post-acute symptoms and their duration. To approximate the proportion of patients that could develop post-acute COVID-19 related symptoms, the same Beta distribution approach as described earlier was applied. Uncertainty related to the duration of post-acute symptoms, which was described by a truncated log-normal distribution, was accounted for by (1) sampling from a normal distribution with mean equal to the µ and σ parameter estimates and standard deviation the standard error around these estimates to account for the uncertainty around the parameter estimates of the log-normal distribution, and by (2) sampling a random duration based on a log-normal distribution characterized by the sampled parameters in step (1). Uncertainty intervals were defined as the 2.5th and 97.5th percentile of the specific disease estimated distributions.

Results

In Belgium, the total observed person-time of COVID-19 cases in all health states was estimated at 27,499 person-years (95% uncertainty interval: [24,740 − 30,086]) in 2020, of which hospitalizations in intensive care accounted for 362 person-years [338–385], community-based cases with severe symptomatology accounted for 1,684 person-years [1,598-1,765], community-based cases with mild/moderate symptomatology accounted for 12,880 person-years [10,445 − 14,909], and post-acute symptoms accounted for 12,574 person-years [11,101 − 14,002]. In 2021, the observed person-time of all COVID cases was estimated at 49,732 person-years [45,032–54,911], of which hospitalizations in intensive care accounted for 432 person-years [407–457], community-based cases with severe symptomatology accounted for 1,610 person-years [1,531-1,695], community-based cases with moderate symptomatology accounted for 15,620 person-years [12,815 − 18,452], and post-acute symptoms accounted for 32,070 person-years [28,106 − 36,030]. The number of reported COVID-related deaths was 19,846 in 2020, and 8,557 in 2021. A visual representation of the time spent in each health state by age and time can be found in Appendix 5 A and 5B.

As shown in Table 2 and depicted in Fig. 1, COVID-19 has led to 253,577 DALYs [252,541 − 254,739] in 2020 and 139,281 DALYs [136,704 − 142,306] in 2021. Most of the disease burden was associated with premature mortality (proportion of YLLs equal to 98.5% in 2020 and 94.0% in 2021). The disease burden was unequally distributed among age categories with 45.9% and 29.6% of the total disease burden occurring in the oldest (80+) age category in 2020 and 2021, respectively, and 41.9% and 49.7% in people with an age between 60 and 79 years in 2020 and 2021, respectively. The total disease burden was the lowest in the population less than 20 years old (0.3% and 0.4% in 2020 and 2021, respectively).

Table 2 Estimated disease burden (mean and 95% uncertainty interval) of COVID-19 in Disability-Adjusted Life Years (DALYs), Years Lived with Disability (YLDs) and Years of Life Lost (YLLs) by age category and total for the calendar years 2020 and 2021
Fig. 1
figure 1

Total estimated disease burden of COVID-19 expressed as Disability-Adjusted Life Years (DALYs) per age category for the period 2020–2021. YLLs = Years of Life Lost; YLDs = Years Lived with Disability

Morbidity (years lived with disability)

In 2020, a total of 3,863 YLDs [2,826-5,025] were observed due to acute symptomatic COVID-19 infections and post-acute symptoms after infection including mild/moderate cases (654 YLDs; 16.9%), severe cases (223 YLDs; 5.8%), critical cases (236 YLDs; 6.1%), and post-acute cases (2,750 YLDs; 71.2%). In 2021, a total of 8,303 YLDs [5,726 − 11,327] were estimated due to acute symptomatic COVID-19 infections including mild/moderate cases (791 YLDs; 9.5%), severe cases (213 YLDs; 2.6%), critical cases (213 YLDs; 3.4%), and post-acute cases (7,016 YLDs; 84.5%). The percentage of observed DALYs attributable to morbidity is highest among people younger than 20, whereby it was higher in 2021 (72.4%) as compared to 2020 (23.3%). In the group aged between 20 and 39, 15.1% of the disease burden was attributable to YLDs in 2020 and 33.1% in 2021. In the older age groups, YLDs accounted for less than 15.0% of the total disease burden in 2020 and 2021.

Mortality (years of Life Lost)

During the first year of the pandemic, 249,714 years of life were lost due to premature mortality in Belgium. In 2021, there was a significant decrease in premature mortality, with an estimated 130,979 YLLs. The YLLs per death was estimated at 12.6 in 2020 compared to 15.3 in 2021. The percentage of DALYs attributable to mortality was the highest among people aged between 60 and 79 years, and older than 80 years. In 2020, it was estimated at 98.9% and 99.4% in the age groups of 60 and 79 years, and older than 80 years, respectively. In 2021 it was estimated at 96.5% and 97.2% in the age groups of 60 and 79 years, and older than 80 years, respectively. Over 80.0% of the YLLs occurred in the group of people with an age over 60 years in 2020 and 2021. YLLs estimates based on the Belgian national life expectancy table in comparison with the GBD 2019 table can be found in Appendix 3.

Discussion

This study aimed to estimate the direct COVID-19 disease burden in Belgium during the first two calendar years of the pandemic. We estimated a total COVID-19 disease burden of 253,577 DALYs [252,541 − 254,739] in 2020, and a total of 139,281 DALYs [136,704 − 142,306] in 2021. More than 80% of the observed disease burden occurred in the population of 60 years of age and older. Of the total observed disease burden, more than 90% could be attributed to premature mortality.

Compared to the disease burden estimates for 2019 by the GBD, the disease burden of COVID-19 in 2020 was 40,254 DALYs higher than ischemic heart disease, the disease with the largest contribution to the overall disease burden, and 183,277 DALYs more than lower respiratory infections [33]. Note that COVID-19 only contributed to the disease burden in the final 10 months of the year 2020, whereas the estimates for these two other causes covered the full year in 2019. If disease burden estimates for all other causes in 2019 could be transferred to 2020, COVID-19 would be the most important cause of disease burden in 2020. The COVID-19 disease burden substantially dropped from 253,577 DALYs in 2020 to 139,281 DALYs in 2021, and COVID-19 would then rank as the 6th most important cause of illness in 2021, in between chronic obstructive pulmonary disorder (142,937 DALYs) and stroke (134,777 DALYs) [33]. In addition, COVID-19 would supersede the impact of lower respiratory infections, which is the next group of vaccine-preventable diseases in the top 20, and ranked at place 13th in 2019. When limiting to the morbidity component (YLDs), the disease burden of COVID-19 in 2020 could be ranked as 57th between gout and urticaria, while in 2021, the observed disease burden of COVID-19 could be ranked as 44th between other sense organ diseases and viral skin diseases. These comparisons should however be interpreted with caution, as it is so far not yet clear how the COVID crisis affected the disease burden of non-COVID conditions.

Our results show that, although public health measures such as telework, mask-wearing, isolation, test strategies, social distancing, and travel restrictions were adopted in Belgium to limit the impact of COVID-19 on population health (see Appendix 6 for more details), the estimated burden in 2020 remained high [34]. Although our study does not allow to infer causality for the large drop in DALYs in 2021 compared to 2020, the biomedical plausibility that the vaccination campaign, which started in Belgium on the 28th of December 2020, caused the change in remaining disease burden is high. Indeed, 2021 saw increasing vaccine uptake, reduced social distancing and other public health measures, more severe SARS-CoV-2 variants and a shift in burden towards less lethal, less severe and more prevalent COVID-19 cases, as evidenced by the randomised controlled vaccine trials [35]. In addition, the build-up of natural immunity against SARS-CoV-2 could also play an important role in the observed drop. Even though social distancing and other public health measures were reduced in 2021 and more severe SARS-CoV-2 variants emerged, we observed a shift in burden towards less lethal, less severe, and more prevalent COVID-19 cases, as one would expect from the evidence provided by the randomized controlled vaccine trials. This could also partly explain the increased YLD in 2021 compared to 2020.

In Europe, DALY estimates have been calculated using the European Burden of Disease Network consensus disease model for COVID-19 during the first calendar year of the pandemic in Denmark, Ireland, Scotland, Germany, Malta, the Netherlands, and France [6, 36,37,38,39,40,41]. Caution is required when comparing burden of disease estimates across countries, as differences between countries in terms of COVID-19 disease burden are partly explained by differences in data that is available in those countries. Although some countries applied models that account for possible over- or underreporting of symptomatic cases or deaths, the data that is sourced in the consensus disease model remains vulnerable to the country-specific situation. For example, using a broad definition for COVID-19 related death or differences in testing strategy in the first months of the pandemic will affect the disease burden estimates. The estimated COVID-19 burden in Belgium in 2020 (2,206.4; 95% uncertainty interval: 2,197.4-2,216.5) was higher compared to the reported burden in other European countries in 2020, our estimated total burden per 100,000 inhabitants was 1.18 times higher than the burden estimated for Scotland [37], 1.35 times higher compared to the Netherlands [40], 2.03 times higher than Malta [42], and 1.50 times higher than France [41]. Other countries (i.e., Ireland, Germany, and Denmark) either used a different standardized life table or did not consider post-acute symptoms. Although the share of the disease burden attributable to morbidity is rather small, the study for Scotland also applied an SEIR compartmental model to estimate the number of community-based cases, which might partly explain the relatively higher burden estimates for Scotland and Belgium during the first year of the pandemic. In 2020, the test and tracing strategy might not have been optimized, resulting in many missed reported cases during the first months of the pandemic. In addition, the death count in Belgium is based on a more inclusive definition for a COVID-19-related death compared to other European countries [2, 26, 27]. To date, no information on COVID-19 mortality is available based on death certificates as the compilation of causes of death from these certificates has a delay of 2 to 3 years. Once available, a validity check could be performed against these certificates.

Although the incidence of observed COVID-19 infections remained high in 2021, the overall age-standardized disease burden decreased from 2,206.4 [2,197.4-2,216.5] per 100,000 inhabitants in 2020 to 1,204.4 [1,182.1-1,230.7] per 100,000 inhabitants. This drop is mainly attributable to changes in mortality as only a relatively small share (< 10.0%) of the total burden of COVID-19 was attributable to morbidity. Although the share of morbidity is lower compared to mortality in the estimated DALYs, there is a danger of over-interpreting YLDs health estimates in isolation. For example, preventing cases with severe symptomatology would probably not only decrease the YLDs related to severe symptomatology but also prevent YLLs related to severe symptomatology, as patients with severe symptoms were more likely to decease. The share of the disease burden attributable to morbidity was substantially larger in the younger group. Studies have shown that there is a strong age-mortality relation in the population of patients diagnosed with COVID-19, with a higher risk to decease in older ages [43]. However, the group of adolescents and young adults are more vulnerable in terms of mental health after an acute COVID-19 infection compared to, for example, older individuals [44, 45]. However, this was not accounted for in this study. The total disease burden of COVID-19 attributable to its post-acute consequences has increased from 71.2% to 2020 to 84.5% in 2021, a finding similar to other European reports. Scotland, for example, reported the highest contribution of post-acute consequences to YLDs (76%) [37], followed by Malta (60%) [42].

Strengths and limitations

This is the first study that estimated the disease burden of COVID-19 in Belgium up to 2021. We successfully combined national data sources and models on COVID-19 infections, hospitalisations, seroprevalence, variants of concern, vaccination and mortality to capture the direct impact of COVID-19 on population health. Uncertainties around the disease model and the person-year estimates were considered, and we have highlighted their impact on YLDs and DALYs estimates. To estimate the YLLs, the standard GBD life table was used to facilitate the comparison of these results with other studies, achieve consistency with other cause-specific disease burden estimates. However, we acknowledge the ongoing disagreements on the use of standard life table [46]. Therefore, we have performed a sensitivity analysis (which can be found in Appendix 3) by comparing the estimated YLLs based on the GBD with the estimated YLLs based on the Belgian national life table. Lastly, we adopted a standardized approach based on a consensus model that allowed us to produce comparable estimates with other countries.

Besides its strengths, the current study had some limitations. First, the current Belgian compartmental model made, for example, no distinction between men and women, which made between-sex comparisons of the burden of disease impossible. The stochastic model is calibrated on hospitalizations, which is not prone to testing strategies, but assumes a direct link between hospitalizations and cases. Next, the YLLs were estimated based on an ad hoc COVID-19 mortality surveillance set up during a crisis situation including a broad definition of COVID-19, and not using the standard practice of ICD-10 nomenclature. These results may therefore over- or underestimate the YLLs, however, until causes of death estimates are available, this assessment cannot be made. Nevertheless, this approach allowed to estimate the COVID-19 mortality during the pandemic based on the latest available data. Third, YLDs of post-acute consequences of COVID-19 were calculated using transition probabilities based on a new cohort with only a limited time of follow-up. The estimated cases of post-acute symptoms of COVID-19 should be interpreted with caution, as there is an ongoing debate on the definition of post-acute consequences of COVID-19 [47, 48], and if these post-acute consequences are causally related to the acute infection. These discussions contribute to the widely varying estimates that have been reported so far [49]. In the methodology to estimate the disease progression, we fitted a parametric survival model with a log-normal distribution, which is only an approximation of the duration in symptoms. As these underlying assumptions might not hold, their corresponding estimates should be interpreted with caution. Lastly, we do not account for comorbidities but we estimate the disease burden of COVID-19 by age, which we assume as proxy for the frailty of patients (e.g., by suffering from comorbidities).

Implications for public health

Our estimates of DALYs due to COVID-19 in 2020 and 2021, and its YLLs and YLDs composites, highlight the direct and immediate impact in terms of mortality and morbidity in acute infections and post-acute consequences of COVID-19 among the Belgian population. This study highlights the importance of available data collection practices especially in the context of more specific information on symptomatology, symptom duration, and post-acute consequences of COVID-19. These estimates are useful for informing public health decision-makers regarding the COVID-19 burden among subgroups of the population who were the most affected by mortality and to target the high-risk sub-groups to limit the impact of COVID-19 on mortality estimates. The current study identified that COVID-19 likely caused more disease burden during the first calendar year (2020) compared to the second year (2021), which emphasizes the importance of preparedness for potential future pandemics in an attempt to reduce the initial burden of a pandemic during its initial months, when there is no natural immunity against the virus, or no vaccine available. Furthermore, it is important to not only evaluate the direct impact of the COVID-19 pandemic on population health, but also study the long-term indirect impact such as the disruption of the utilisation of vital health care services [50], increased mental health problems [51], risk of increase in health inequalities [52], and economic losses. Therefore, future studies should include measures at regional levels, and attempt to measure the indirect impact of COVID-19 as a risk factor on the collateral damage to health care services, mental health and increase of incidences of non-communicable diseases associated with COVID-19.

Conclusion

COVID-19 had an important effect on the population health in Belgium during its first two calendar years. Most of the direct population health loss was attributable to mortality, and mainly carried by people aged 60 and above. Our results highlight that the post-acute consequences of COVID-19 might strongly contribute to COVID-19 morbidity and can have a considerable impact on population health with, consequently, a higher burden on the health care system and health-related quality of life. Further research is required to evaluate the direct longer term and indirect impact of COVID-19 and study potential inequalities associated with the health burden of COVID-19.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the COVID-19 repository from Sciensano, https://epistat.sciensano.be/covid/.

References

  1. WHO/Europe | Coronavirus Disease (COVID-19) outbreak—about the virus. https://www.euro.who.int/en/health-topics/health-emergencies/coronavirus-covid-19/novel-coronavirus-2019-ncov. Accessed 20 May 2022.

  2. Bustos Sierra N, Bossuyt N, Braeye T, Leroy M, Moyersoen I, Peeters I, et al. All-cause mortality supports the COVID-19 mortality in Belgium and comparison with major fatal events of the last century. Arch Public Health. 2020;78:1–8.

    Article  Google Scholar 

  3. Epistat—COVID-19 monitoring. https://epistat.wiv-isp.be/covid/. Accessed 20 May 2022.

  4. Murray CJ, Lopez AD, Organization WH, et al. The global burden of disease: a comprehensive assessment of mortality and disability from diseases, injuries, and risk factors in 1990 and projected to 2020: summary. World Health Organization; 1996.

  5. Lopez AD, Mathers CD, Ezzati M, et al. Global Burden of Disease and Risk Factors. Washington (DC): The International Bank for Reconstruction and Development / The World Bank; 2006. Available from: https://www.ncbi.nlm.nih.gov/books/NBK11812/. Co-published by Oxford University Press, New York.

  6. Haneef R, Schmidt J, Gallay A, Devleesschauwer B, Grant I, Rommel A, et al. Recommendations to plan a national burden of disease study. Arch Public Health. 2021;79:1–8.

    Article  Google Scholar 

  7. Guan W, Ni Z, Hu Y, Liang W, Ou C, He J, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. 2020;382:1708–20.

    Article  CAS  PubMed  Google Scholar 

  8. Soriano JB, Murthy S, Marshall JC, Relan P, Diaz JV; WHO Clinical Case Definition Working Group on Post-COVID-19 Condition. A clinical case definition of post-COVID-19 condition by a Delphi consensus. Lancet Infect Dis. 2022;22(4):e102-7. https://doi.org/10.1016/S1473-3099(21)00703-9.

  9. Xu E, Xie Y, Al-Aly Z. Long-term neurologic outcomes of COVID-19. Nat Med. 2022;28(11):2406-15. https://doi.org/10.1038/s41591-022-02001-z.

  10. Taquet M, Sillett R, Zhu L, Mendel J, Camplisson I, Dercon Q, et al. Neurological and psychiatric risk trajectories after SARS-CoV-2 infection: an analysis of 2-year retrospective cohort studies including 1 284 437 patients. Lancet Psychiatry. 2022;9:815–27.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Beutels P, Edmunds WJ, Smith R. Partially wrong? Partial equilibrium and the economic analysis of public health emergencies of international concern. Health Econ. 2008;17:1317–22.

    Article  CAS  PubMed  Google Scholar 

  12. Rens E, Smith P, Nicaise P, Lorant V, Van den Broeck K. Mental distress and its contributing factors among young people during the first wave of COVID-19: a belgian survey study. Front Psychiatry. 2021;12: 575553.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Douglas M, Katikireddi SV, Taulbut M, McKee M, McCartney G. Mitigating the wider health effects of covid-19 pandemic response. BMJ. 2020;369:m1557.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Danhieux K, Buffel V, Pairon A, Benkheil A, Remmen R, Wouters E, et al. The impact of COVID-19 on chronic care according to providers: a qualitative study among primary care practices in Belgium. BMC Fam Pract. 2020;21:1–6.

    Article  Google Scholar 

  15. De Pelsemaeker M-C, Guiot Y, Vanderveken J, Galant C, Van Bockstal MR. The impact of the COVID-19 pandemic and the associated Belgian governmental measures on cancer screening, surgical pathology and cytopathology. Pathobiology. 2021;88:46–55.

    Article  PubMed  Google Scholar 

  16. Wyper GM, Assunção RM, Colzani E, Grant I, Haagsma JA, Lagerweij G, et al. Burden of disease methods: a guide to calculate COVID-19 disability-adjusted life years. Int J Public Health. 2021;66:619011.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Abrams S, Wambua J, Santermans E, Willem L, Kuylen E, Coletti P, et al. Modelling the early phase of the Belgian COVID-19 epidemic using a stochastic compartmental model and studying its implied future trajectories. Epidemics. 2021;35: 100449.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Smith P, Proesmans K, Van Cauteren D, Demarest S, Drieskens S, De Pauw R, et al. Post COVID-19 condition and its physical, mental and social implications: protocol of a 2-year longitudinal cohort study in the Belgian adult population. Arch Public Health. 2022;80:1–12.

    Article  Google Scholar 

  19. Coletti P, Wambua J, Gimma A, Willem L, Vercruysse S, Vanhoutte B, et al. CoMix: comparing mixing patterns in the Belgian population during and after lockdown. Sci Rep. 2020;10:1–10.

    Article  Google Scholar 

  20. Chen C, Haupert SR, Zimmermann L, Shi X, Fritsche LG, Mukherjee B. Global prevalence of post COVID-19 condition or Long COVID: a meta-analysis and systematic review. J Infect Dis. 2022;226(9):1593–607.

    Article  CAS  PubMed  Google Scholar 

  21. Davis HE, McCorkell L, Vogel JM, Topol EJ. Long COVID: major findings, mechanisms and recommendations. Nat Rev Microbiol. 2023:1–14.

  22. Smith P, Charafeddine R, Drieskens S, De Pauw R, De Ridder K, Demarest S, et al. COVIMPACT studie: long COVID en de lichamelijke, psychische en sociale gevolgen – resultaten van 3 maanden opvolging. Brussels: Sciensano; 2021.

    Google Scholar 

  23. Schroeder SA. Incidence, prevalence, and hybrid approaches to calculating disability-adjusted life years. Popul Health Metr. 2012;10:1–7.

    Article  Google Scholar 

  24. Salomon JA, Haagsma JA, Davis A, de Noordhout CM, Polinder S, Havelaar AH, et al. Disability weights for the global burden of disease 2013 study. Lancet Glob Health. 2015;3:e712–723.

    Article  PubMed  Google Scholar 

  25. Haagsma JA, Maertens de Noordhout C, Polinder S, Vos T, Havelaar AH, Cassini A, et al. Assessing disability weights based on the responses of 30,660 people from four European countries. Popul Health Metr. 2015;13:1–15.

    Article  Google Scholar 

  26. Renard F, Scohy A, Van der Heyden J, Peeters I, Dequeker S, Vandael E, et al. Establishing an ad hoc COVID-19 mortality surveillance during the first epidemic wave in Belgium, 1 March to 21 June 2020. Euro Surveill. 2021;26:2001402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Molenberghs G, Faes C, Verbeeck J, Deboosere P, Abrams S, Willem L, et al. COVID-19 mortality, excess mortality, deaths per million and infection fatality ratio, Belgium, 9 March 2020 to 28 June 2020. Euro Surveill. 2022;27:2002060.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Global Burden of Disease Collaborative Network. Global Burden of Disease Study 2019 (GBD 2019) Reference Life Table. 2021.

  29. Statistics Belgium. Bevolking naar woonplaats, nationaliteit, burgerlijke staat, leeftijd en geslacht. 2022. https://statbel.fgov.be/nl/open-data/bevolking-naar-woonplaats-nationaliteit-burgerlijke-staat-leeftijd-en-geslacht-10.

  30. Bilcke J, Beutels P, Brisson M, Jit M. Accounting for methodological, structural, and parameter uncertainty in decision-analytic models: a practical guide. Med Decis Making. 2011;31:675–92.

    Article  PubMed  Google Scholar 

  31. Powell M, Ebel E, Schlosser W. Considering uncertainty in comparing the burden of illness due to foodborne microbial pathogens. Int J Food Microbiol. 2001;69:209–15.

    Article  CAS  PubMed  Google Scholar 

  32. Devleesschauwer B, Torgerson P, Charlier J, Levecke B, Praet N, Roelandt S, et al. Prevalence: tools for prevalence assessment studies. 2022.

  33. Vos T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, et al. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the global burden of disease study 2019. Lancet. 2020;396:1204–22.

    Article  Google Scholar 

  34. Haug N, Geyrhofer L, Londei A, Dervic E, Desvars-Larrive A, Loreto V, et al. Ranking the effectiveness of worldwide COVID-19 government interventions. Nat Hum Behav. 2020;4:1303–12.

    Article  PubMed  Google Scholar 

  35. Rotshild V, Hirsh-Raccah B, Miskin I, Muszkat M, Matok I. Comparing the clinical efficacy of COVID-19 vaccines: a systematic review and network meta-analysis. Sci Rep. 2021;11:22777.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Rommel A, von der Lippe E, Plass D, Ziese T, Diercke M, an der Heiden M, et al. The COVID-19 disease burden in Germany in 2020: years of life lost to death and disease over the course of the pandemic. Dtsch Ärztebl Int. 2021;118:145.

    PubMed  PubMed Central  Google Scholar 

  37. Wyper G, Fletcher E, Grant I, McCartney G, Fischbacher C, Harding O, et al. Measuring disability-adjusted life years (DALYs) due to COVID-19 in Scotland, 2020. Arch Public Health. 2022;80:1–10.

    Article  Google Scholar 

  38. Pires SM, Redondo HG, Espenhain L, Jakobsen LS, Legarth R, Meaidi M, et al. Disability adjusted life years associated with COVID-19 in Denmark in the first year of the pandemic. BMC Public Health. 2022;22:1–9.

    Article  Google Scholar 

  39. Moran DP, Pires SM, Wyper G, Devleesschauwer B, Cuschieri S, Kabir Z. Estimating the direct disability-adjusted life years associated with SARS-CoV-2 (COVID-19) in the Republic of Ireland: the first full year. Int J Public Health. 2022;83.

  40. McDonald SA, Lagerweij GR, de Boer P, de Melker HE, Pijnacker R, Gras LM, et al. The estimated disease burden of acute COVID-19 in the Netherlands in 2020, in disability-adjusted life-years. 2021.

  41. Haneef R, Fayad M, Fouillet A, Sommen C, Bonaldi C, Wyper GM, et al. Direct impact of COVID-19 by estimating disability-adjusted life years at national level in France in 2020. PLoS One. 2023;18:e0280990.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Cuschieri S, Calleja N, Devleesschauwer B, Wyper G. Estimating the direct Covid-19 disability-adjusted life years impact on the Malta population for the first full year. BMC Public Health. 2021;21:1–7.

    Article  Google Scholar 

  43. Dessie ZG, Zewotir T. Mortality-related risk factors of COVID-19: a systematic review and meta-analysis of 42 studies and 423,117 patients. BMC Infect Dis. 2021;21:855.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Singh S, Roy D, Sinha K, Parveen S, Sharma G, Joshi G. Impact of COVID-19 and lockdown on mental health of children and adolescents: a narrative review with recommendations. Psychiatry Res. 2020;293: 113429.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Penninx BW, Benros ME, Klein RS, Vinkers CH. How COVID-19 shaped mental health: from infection to pandemic effects. Nat Med. 2022;22:1–11.

    Google Scholar 

  46. Devleesschauwer B, McDonald SA, Speybroeck N, Wyper G. Valuing the years of life lost due to COVID-19: the differences and pitfalls. Int J Public Health. 2020;65:719–20.

    Article  PubMed  Google Scholar 

  47. Fernández-de-Las-Peñas C, Palacios-Ceña D, Gómez-Mayordomo V, Cuadrado ML, Florencio LL. Defining post-COVID symptoms (post-acute COVID, long COVID, persistent post-COVID): an integrative classification. Int J Environ Res Public Health. 2021;18: 2621.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Fernández-de-Las-Peñas C. Long COVID: current definition. Infection. 2021:1–2.

  49. Fernández-de-Las-Peñas C, Palacios-Ceña D, Gómez-Mayordomo V, Florencio LL, Cuadrado ML, Plaza-Manzano G, et al. Prevalence of post-COVID-19 symptoms in hospitalized and non-hospitalized COVID-19 survivors: a systematic review and meta-analysis. Eur J Intern Med. 2021;92:55–70.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Moynihan R, Sanders S, Michaleff ZA, Scott AM, Clark J, To EJ, et al. Impact of COVID-19 pandemic on utilisation of healthcare services: a systematic review. BMJ Open. 2021;11: e045343.

    Article  PubMed  Google Scholar 

  51. Bruggeman H, Smith P, Berete F, Demarest S, Hermans L, Braekman E, et al. Anxiety and depression in Belgium during the first 15 months of the COVID-19 pandemic: a longitudinal study. Behav Sci. 2022;12: 141.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Bambra C, Riordan R, Ford J, Matthews F. The COVID-19 pandemic and health inequalities. J Epidemiol Community Health. 2020;74:964–8.

    PubMed  Google Scholar 

Download references

Acknowledgements

This study was conducted within the framework of the Belgian National Burden of Disease Study (BeBOD), coordinated by Sciensano, the Belgian institute for health.

Funding

This study received funding from the European Union’s Horizon 2020 research and innovation programme - project EpiPose (Grant Agreement Number 101003688). This work reflects only the authors’ view. The European Commission is not responsible for any use that may be made of the information it contains.

Author information

Authors and Affiliations

Authors

Contributions

BDV, LW, and RDP conceptualized the paper, performed the main analysis, and took the lead on writing the original draft. JJ contributed to the methodology of the mortality data. PS was responsible for the data and analysis on the post-acute phase. BDV, GW and SMP contributed to the burden of disease protocol and validation of the results. NVG contributed to the delivery and validation of the data on hospitalizations. PB, SA, NF, NH contributed to the analysis of cases by the implementation of the extended compartmental model. SA contributed to reviewing and validating the analysis. DVC contributed to the validation of the data on the COVID-19 cases. JJ, PS, AS, GW, SMP, NVG, PB, NF, SA, DVC, NS, NH reviewed and edited the final draft of the paper. Funding for this analysis was indirectly acquired by BDV, LW, RDP, PB, SA, and NH.

Corresponding author

Correspondence to Robby De Pauw.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1. Appendix 1.

Modeling long-lasting symptoms after a confirmed SARS-CoV-2 infection based on the longitudinal follow-up data from the Belgian COVIMPACT study. Appendix 2. Imputation of the COVID-19 mortality dataset. Appendix 3. Comparison of YLL estimates of the COVID-19 burden in Belgium using the Global Burden of Disease Remaining Life Expectancy Table (2019) and the Belgian National Life Expectancy Table (2019). Appendix 4. Overview of distributions of disability weights approximated by a beta distribution using the 95% uncertainty interval around the DW. Appendix 5A. Figure with an overview of prevalent cases and deaths in each health state for the year 2020 by day and age category. Appendix 5B. Figure with an overview of prevalent cases and deaths in each health state for the year 2021 by day and age category. Appendix 6. Table with an overview of important public health measures taken in response to the COVID-19 crisis during the period 2020-2021.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Devleesschauwer, B., Willem, L., Jurčević, J. et al. The direct disease burden of COVID-19 in Belgium in 2020 and 2021. BMC Public Health 23, 1707 (2023). https://doi.org/10.1186/s12889-023-16572-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-16572-0

Keywords