Skip to main content

Excess hospitalizations and mortality associated with seasonal influenza in Spain, 2008–2018

Abstract

Background

Influenza may trigger complications, particularly in at-risk groups, potentially leading to hospitalization or death. However, due to lack of routine testing, influenza cases are infrequently coded with influenza-specific diagnosis. Statistical models using influenza activity as an explanatory variable can be used to estimate annual hospitalizations and deaths associated with influenza. Our study aimed to estimate the clinical and economic burden of severe influenza in Spain, considering such models.

Methods

The study comprised ten epidemic seasons (2008/2009–2017/2018) and used two approaches: (i) a direct method of estimating the seasonal influenza hospitalization, based on the number of National Health Service hospitalizations with influenza-specific International Classification of Diseases (ICD) codes (ICD-9: 487–488; ICD-10: J09-J11), as primary or secondary diagnosis; (ii) an indirect method of estimating excess hospitalizations and deaths using broader groups of ICD codes in time-series models, computed for six age groups and four groups of diagnoses: pneumonia or influenza (ICD-9: 480–488, 517.1; ICD-10: J09–J18), respiratory (ICD-9: 460–519; ICD-10: J00–J99), respiratory or cardiovascular (C&R, ICD-9: 390–459, 460–519; ICD-10: I00–I99, J00–J99), and all-cause. Means, excluding the H1N1pdm09 pandemic (2009/2010), are reported in this study.

Results

The mean number of hospitalizations with a diagnosis of influenza per season was 13,063, corresponding to 28.1 cases per 100,000 people. The mean direct annual cost of these hospitalizations was €45.7 million, of which 65.7% was generated by patients with comorbidities. Mean annual influenza-associated C&R hospitalizations were estimated at 34,894 (min: 16,546; max: 52,861), corresponding to 75.0 cases per 100,000 (95% confidence interval [CI]: 63.3–86.3) for all ages and 335.3 (95% CI: 293.2–377.5) in patients aged ≥ 65 years. We estimate 3.8 influenza-associated excess C&R hospitalizations for each hospitalization coded with an influenza-specific diagnosis in patients aged ≥ 65 years. The mean direct annual cost of the estimated excess C&R hospitalizations was €142.9 million for all ages and €115.9 million for patients aged ≥ 65 years. Mean annual influenza-associated all-cause mortality per 100,000 people was estimated at 27.7 for all ages.

Conclusions

Results suggest a relevant under-detected burden of influenza mostly in the elderly population, but not neglectable in younger people.

Peer Review reports

Introduction

Influenza is a viral infection accountable for seasonal epidemics, leading to a substantial disease burden throughout the world [1]. The infection not only triggers mild symptoms but also causes severe complications, which may result in hospitalization or even death of the patient [1,2,3].

Some patients are at a higher risk of developing severe complications due to their age or medical conditions. The World Health Organization (WHO) recommends that those countries initiating or expanding their programs for seasonal influenza vaccination should consider the following targets for vaccination (not in the order of priority): health workers, individuals with comorbidities and underlying conditions, older adults, and pregnant women [4]. The WHO also states that countries may consider additional (sub)populations for vaccination, such as children, and other groups at a high risk of severe influenza in congregate-living settings [4].

Increasing seasonal vaccination within priority groups is a key strategy to lessen the burden of influenza [5]. Prior to the COVID-19 pandemic, Spanish official estimates indicated vaccination coverage rates (VCRs) in those aged 65 or above below the 75% VCR goal [6, 7]. During 2008/2009, VCR in this at-risk group was 65.4%. It has gradually declined over the following nine years, reaching 55.7% in 2017/2018 [8].

The intensity and severity of influenza epidemics are variable and often unpredictable, depending on the VCRs, the circulating strains of influenza viruses, and vaccine effectiveness against the circulating virus, among other variables [3, 9, 10]. Therefore, it is important to continuously monitor the morbidity and mortality, and specific impacts of influenza by age, risk groups, and groups of potential complications [3, 11].

Health administrative data can be used to efficiently assess the burden of diseases in broad populations [12]. These algorithms may be applied to hospitalized patients even when laboratory data are not available and tests are not routinely performed, to help improve the epidemiologic studies of influenza. However, there may be variations in terms of the national coding practices [12]. In Canada, a study by Hamilton et al. reports that even the best performing influenza algorithm had a specificity of 99% but a sensitivity of 73% [12]. The challenges of estimating the disease burden of influenza infection in the context of insufficient influenza laboratory testing and coding practices have led to an increasing use of regression techniques to estimate excess morbidity and mortality due to influenza [13].

In Spain, available studies are mostly focused on hospitalized cases coded with an International Classification of Diseases (ICD) influenza code and 2014/15 is the most recently analyzed season [14,15,16,17,18,19,20]. There are no national estimates available for morbidity and, those available for mortality during this study, were not updated [21, 22]. León-Gómez et al. (2015) published estimates on Pneumonia and Influenza (P&I) excess mortality for population aged ≥ 65 years, but these include only data until season 2011/2012 [21]. The annual estimates published by the surveillance system include a characterization of confirmed severe influenza cases that are hospitalized and indicates excess mortality periods as well. However, influenza-associated excess deaths were not quantified at the time of this study neither in total nor by age.

There is an undeniable need for updated national estimates on the influenza-associated excess hospitalization and mortality in Spain, considering the impact of influenza on other respiratory and non-respiratory complications. The contrast of these estimates with those obtained through the direct analysis of influenza-coded cases can be of added value, to understand the influenza under-detection factor, and its evolution over the course of time.

This study aims to address this knowledge gap by estimating the clinical and economic burden of severe influenza in Spain from season 2008/09 to season 2017/18 across distinct age groups, by: (i) describing and estimating the cumulative seasonal incidence of Spanish National Health Service (NHS) hospitalizations with a diagnosis of influenza by age group and comorbidities; (ii) estimating the number of excess influenza-associated hospitalizations and deaths during influenza epidemics, considering broader groups of diagnoses.

Methods

Study design

The Burden of Acute Respiratory Infections (BARI) is a real-world evidence study assessing the clinical and economic burden of acute respiratory infections (influenza and respiratory syncytial virus) in Spain and Portugal [23,24,25,26]. We are reporting the results for the burden of severe influenza in Spain, measured through hospitalizations and deaths. Results from a similar analysis performed for Portugal are reported by Froes et al. [24].

This publication uses the following two approaches: (i) a direct method of estimating seasonal influenza incidence, based on the number of NHS hospitalizations with influenza-specific ICD codes; (ii) an indirect method of estimating excess hospitalizations and deaths using broader groups of ICD codes in time-series ecological models.

Hospitalizations with a diagnosis of influenza

Influenza episodes were defined as those coded with ICD-9-MC 487 or 488, or ICD-10-ES J09, J10 or J11, in any primary or secondary diagnosis field. These codes were found to be highly accurate ICD algorithms to identify influenza hospitalizations as per Hamilton et al. [12]. Total hospitalizations excluded admissions related to routine birth, programmed activity, and those admissions in which the primary diagnosis was related to musculoskeletal, alcohol, or mental disease. Diagnostic information from the influenza episode was also used to identify individuals who had at least one medical condition regarded as a risk factor for severe influenza (Additional file 1: Table S1), in any primary or secondary diagnosis field, considering the following conditions: pregnancy, diabetes mellitus, respiratory or lung disease, cardiovascular, immunocompromised, chronic liver disease, and chronic kidney disease.

Influenza-associated excess hospitalization and mortality

According to the primary diagnosis, influenza-associated excess hospitalization was computed for four groups of diagnoses [27, 28]: pneumonia or influenza (P&I, ICD-9-MC: 480–488, 517.1; or ICD-10-ES: J09-J18), respiratory (R, ICD-9-MC: 460–519; or ICD-10-ES: J00-J99), respiratory or cardiovascular (C&R, ICD-9-MC: 390–459, 460–519; or ICD-10-MC: I00-I99, J00-J99), and all-cause (any ICD-9/10 diagnosis). The aforementioned groups of diagnoses and the primary cause of death in the death certificate were used to compute influenza-associated excess mortality.

Data sources

Hospital discharge data

Anonymized administrative data on hospitalizations (January 2008-December 2018) were retrieved from the IASIST Projected Hospitalizations Database (PHDB) for Spain, including public hospital episodes. The Minimum Basic Datasets (MBDS) collect anonymous information from inpatient episodes, including demographic, administrative, and clinical information data, including any ICD diagnosis and procedures coded for every inpatient episode. The PHDB receives MBDS from 189 hospitals of the NHS. The PHDB represents approximately 55% to 65% of the total NHS episodes per year. Inferential methodology is used to estimate the final universe of patients who have been hospitalized for each specific pathology [29]. This methodology estimates data on hospitalization for the total Spanish NHS universe, which were used in this study. The PHDB database was used to identify the hospitalizations with a diagnosis of influenza and the hospitalizations associated with other groups of diagnoses used for the influenza-associated excess hospitalization analysis.

These data were also used to describe the in-hospital case fatality risk based on discharge status, corresponding to the in-hospital death during influenza episode at hospitalization.

Death certificate data

The monthly reports for death per cause per age group, as in the death certificate, were obtained from the Instituto Nacional de Estadística (INE, National Statistics Institute) in October 2020 [30].

Influenza activity data

The primary predictor of influenza excess hospitalizations and deaths was the overall weekly incidence rate of influenza-like-illness (ILI), which was obtained from ScVGE—Sistema centinela de Vigilancia de la Gripe en España (National Epidemiology Surveillance Network) in November 2020 [31]. This rate is estimated by ScVGE based on data collected by each sentinel network belonging to the ScVGE. The ScVGE uses the the clinical definition of influenza—as proposed by the European Union [32]—to compute this rate [33]. Being close to the WHO definition of ILI, it will be treated as such throughout this manuscript [34].

Demographic data

Data on age-specific annual resident population estimates were downloaded from INE’s website to calculate rates the incidence rates of cases per 100,000 people [35].

Statistical analysis

An ecological approach was used to estimate the hospitalization costs, number of influenza-associated excess hospitalizations, and deaths; weekly morbidity and monthly mortality data were explicitly modeled against indicators of influenza activity (weekly moving average, broken down by season with a 1-week lag for hospitalizations and monthly indicator broken down by season for death) using cyclic regression models.

Different Poisson cyclic models (time series) were used, where age- and cause-specific hospitalization and mortality data were explained by the ILI incidence [13, 36, 37], as well as time trends and seasonal terms, using a log link [38, 39]. The model used for hospitalization is given below.

$$\mathbf{l}\mathbf{o}\mathbf{g}\left({\varvec{E}}\left(\frac{{{\varvec{H}}}_{{\varvec{t}}}^{\left({\varvec{c}},{\varvec{a}}\right)}}{{\varvec{P}}{\varvec{o}}{\varvec{p}}}\right)\right)={\boldsymbol{\alpha }}_{0}+{\boldsymbol{\alpha }}_{1}{\varvec{t}}+{\boldsymbol{\alpha }}_{2}{{\varvec{t}}}^{2}+\boldsymbol{ }{\sum }_{{\varvec{i}}=1}^{12}{\left[{{\varvec{\beta}}}_{{\varvec{i}}}{\varvec{c}}{\varvec{o}}{\varvec{s}}\left(\frac{{\varvec{i}}{\varvec{t}}{\varvec{\pi}}}{\mathrm{52,14}}\right)+{{\varvec{\gamma}}}_{{\varvec{i}}}{\varvec{s}}{\varvec{i}}{\varvec{n}}\left(\frac{{\varvec{i}}{\varvec{t}}{\varvec{\pi}}}{\mathrm{52,14}}\right)\right]}_{{\varvec{i}}}+{\varvec{\delta}}{{\varvec{I}}{\varvec{N}}{\varvec{F}}}_{{\varvec{t}}}$$

where \({{\varvec{H}}}_{{\varvec{t}}}^{\left({\varvec{c}},{\varvec{a}}\right)}\) the weekly number of hospitalizations for cause c in a specific age group a, and \({\varvec{P}}{\varvec{o}}{\varvec{p}}\) the population offset. The terms t, t2, cos, and sin are terms representing the time trend and seasonality. \({{\varvec{I}}{\varvec{N}}{\varvec{F}}}_{{\varvec{t}}}\) represents the weekly influenza activity data). A backward selection method was used to identify significant time trends and seasonal terms in the age- and cause-specific models.

Baseline hospitalization and mortality data were calculated using the Poisson models as the model-expected values when the variable ILI incidence was set to zero. The number of influenza-associated excess hospitalizations and mortality was the difference between the expected hospitalizations and deaths, as estimated by the Poisson model, which incorporated the ILI incidence rate as the indicator of influenza activity; the number of hospitalizations and deaths were estimated by the baseline model and the indicator was not considered in this regard. The number of influenza-associated hospitalizations was defined as the sum of the weekly excesses for each epidemic season. The number of influenza-associated deaths was defined as the sum of the monthly excesses. Every year, influenza season was defined as the period from September to June.

The correlation between the values predicted by the model and those observed indicated the performance of the model, using Pearson’s correlation. The mean absolute percentage error (MAPE), defined as the average of the percentage difference between predicted and observed values for all weeks (or months in case of mortality), was also computed. All analyses were performed using SAS 9.4.

For the excess hospitalization and mortality analysis, the hospitalizations and deaths were categorized into four groups of diagnoses (P&I, respiratory, respiratory, or cardiovascular, and all-cause). The population was stratified in the following age groups: 0–4 years, 5–18, 19–49, 50–64, 65–74, ≥ 75, and ≥ 65 years.

Confidence intervals (CIs) of 95% were calculated for the estimated influenza-associated excess hospitalizations and deaths by causes and by age groups. A significance level of 5% was considered for the association of excess estimates with the ILI incidence variable (p-value < 0.05).

The estimated influenza cases were divided by the Spanish population in the respective age group to calculate cases per 100,000 people.

The indicator of cost per unit of hospital production (UHP) was used to compute costs of those hospitalizations with a diagnosis of influenza. This indicator relates the operating costs incurred by the hospital to carry out its activity to the production carried out by the hospital. The UHP value was established for each type and level of hospital in the Spanish National Health System [40, 41]. The average cost per UHP was updated annually through the Hospital TOP 20 Program; the degree of complexity of hospitalization episode, considering variables related to the patient and specific episode, was calculated using the 3 M™ All Patient Refined Diagnosis Related Groups (APR DRGs) system (version 32). The cost of influenza-associated excess hospitalization was computed by multiplying the number of estimated excess hospitalizations by the mean cost per hospitalization, by cause and by age group.

The estimates of influenza-associated excess hospitalization and mortality were compared to the number of primary or secondary influenza‐specific diagnoses captured from NHS hospital discharge records and fatalities observed during these episodes.

The following characteristics of hospitalizations with a diagnosis of influenza—approach (i) – were analyzed: the contribution (in percentage) of each age group or comorbidity status to total hospitalizations with a diagnosis of influenza; the mean length of stay (LoS), calculated as the time from hospitalization until discharge; the in-hospital case fatality risk (corresponding to the hospitalizations with a diagnosis of influenza that resulted in death) during the respective inpatient episode; and the mean hospitalization cost per patient, computed by dividing the total cost of hospitalizations with a diagnosis of influenza by the number of patients hospitalized with a diagnosis of influenza. The analysis was stratified by age group and comorbidity status (with or without comorbidities).

Means are reported only for nine seasons, excluding the H1N1pdm09 pandemic (2009/2010). Three-year means, including the most recent seasons (2015/2016 to 2017/2018), are also included.

Results

Descriptive statistics for unmodeled hospitalizations and deaths

In 2018, Spain had a resident population of 46.7 million people, of which 19.2% were ≥ 65 years old [35]. During the study period, a total of 37.0 million all-cause hospitalizations were estimated in Spanish public hospitals. C&R, R, and P&I diagnoses accounted for 26.4%, 10.8%, and 2.9% of hospitalizations, respectively. Patients aged ≥ 65 years contributed to 44.2% of all-cause, 65.6% of C&R, 55.5% of R, and 65.6% of P&I hospitalizations. A total of 4.0 million people died in Spain over the same period, of whom 3.4 million (84.9%) were ≥ 65 years old.

Influenza epidemiologic burden

Hospitalizations with a diagnosis of influenza

A total of 121,393 hospitalizations with a diagnosis of influenza were registered between 2008/2009 and 2017/2018, of which 83,515 (68.8%) were during the last three analyzed seasons. Excluding 2009/2010, which was affected by the H1N1 pandemic, the mean annual number of hospitalizations was 13,063 per season and 28.1 cases per 100,000 peopleFootnote 1 (Table 1). The maximum number of annual hospitalizations was observed in 2017/18 (37,619; 80.6 per 100,000).

Table 1 Hospitalizations with a diagnosis of influenza by age group and epidemic season, in absolute number of cases and incidence rates per 100,000 people, in Spanish Public Hospitals between 2008/2009 and 2017/2018

Patients aged ≥ 65 years contributed to 56.7% of hospitalizations with a diagnosis of influenza1, with a mean of 88.2 cases per 100,000 over nine seasons, reaching 289.3 cases per 100,000 in 2017/18 (Table 1). Patients with comorbidities accounted for 59.0% of hospitalizations, 32.0% of whom were < 65 years old (Table 2)1.

Table 2 Summary of the Characteristics of the Hospitalizations With a Diagnosis of Influenza by Age Group and Presence of Comorbidities in Spanish Public Hospitals Between 2008/2009 and 2017/2018

The mean LoS stood at 9.4 days1. The highest mean LoS was observed in patients aged < 65 with comorbidities (12.0 days)1. Those patients between 50 and 64 years of age presented the highest mean LoS (11.7 days)1. The mean in-hospital case fatality risk was 5.2%, which increased with age—from 0.3% in those aged < 5 years to 7.8% in those aged ≥ 75 years1. Patients with comorbidities, who were < 65 years old, presented 2.8 times higher mean in-hospital case fatality risk than those patients of the same age and without comorbidities (4.5% vs. 1.6%, respectively)1. These indicators were detailed by age group and presence of comorbidities in Table 2.

Influenza-associated excess hospitalization

The estimated all-age influenza-associated excess hospitalizations over the 10-year study period was 162,691 in the P&I, 282,246 in the R, 326,890 in the C&R, and 280,072 in the all-cause model (Table 3). The estimated excess hospitalizations per group of diagnoses for patients aged ≥ 65 years was 96,517 in the P&I, 205,845 in the R, 253,816 in the C&R, and 235,921 in the all-cause model. The maximum number of influenza-associated excess hospitalizations was observed in 2017/18 for P&I, R, and C&R diagnoses, whereas the maximum excess hospitalizations was estimated to have occurred in 2015/2016 for all-cause diagnoses, followed by 2017/2018. The models using C&R hospitalizations showed the best performance for the all ages and the ≥ 65 years age groups (Additional file 1: Table S2).

Table 3 Estimated Influenza-Associated Excess Hospitalizations for All Ages and for Population Aged 65 Years Old or Over, in Absolute Number of Cases and per 100,000 People, by Group of Diagnoses and Epidemic Season in Spanish Public Hospitals Between 2008/2009 and 2017/2018

The mean annual influenza-associated C&R hospitalizations over nine seasons were estimated at 34,8941 (Table 4). The mean annual excess C&R hospitalization rate per 100,000 for the Spanish population was estimated at 75.0 (95% CI: 63.3–86.3) for all ages and 335.3 (95% CI: 293.2–377.5) in the population aged ≥ 65 years1. In the population aged < 65 years, the highest mean annual excess C&R hospitalization rate per 100,000 was estimated in those aged between 50 to 64 years (54.8 [95% CI: 44.8–64.9])1.

Table 4 Estimated Influenza-Associated Excess Respiratory or Cardiovascular Hospitalizations, in Absolute Number of Cases and per 100,000 People, by Age Group and Epidemic Season in Spanish Public Hospitals Between 2008/2009 and 2017/2018

The absolute and relative influenza-associated excess hospitalizations by age groups and by diagnoses (P&I, respiratory cause alone and all-causes) are displayed in Supplementary Materials (Additional file 1: Table S4).

Influenza-associated excess mortality

The total number of all-age influenza-associated excess deaths in the ten-year study period was 12,063 in the P&I, 38,075 in the R, 74,694 in the C&R, and 117,389 in the all-cause model (Table 5). The estimated excess deaths per group of diagnoses in the population aged ≥ 65 years was 11,010 in the P&I, 35,364 in the R, 68,484 in the C&R, and 108,939 in the all-cause model. The highest estimated influenza-associated deaths in P&I and R groups of diagnoses was seen during 2017/2018. The C&R and all-cause groups showed maximum influenza-associated deaths in 2014/2015. The models using all-cause deaths yielded the best model performance (Additional file 1: Table S3).

Table 5 Estimated Influenza-Associated Excess Deaths for All Ages and for Population Aged 65 Years Old or Over, in Absolute Number of Cases and per 100,000 People, by Group of Diagnoses and Epidemic Season in Spain Between 2008/2009 and 2017/2018

The mean annual influenza-associated all-cause deaths per 100,000 people – over nine seasons –was estimated at 27.7 for the all ages group (Table 6). This ratio (9.7 per 100,000) was 10.8 times higher in those aged 50–64 years compared with the population aged 19–49 years (0.9 per 100,000)1. The population aged ≥ 65 years showed a high annual excess all-cause influenza-associated deaths in five seasons; the specifics are as follows: 2008/09 (14,193 deaths; 187.2 deaths per 100,000 people); 2011/12 (18,744; 232.7), 2014/15 (24,282; 285.1), 2016/17 (19,074; 217.8), and 2017/18 (16,685; 187.7). All these seasons showed a dominance or co-predominance of A(H3N2) and/or B/Yamagata strains.

Table 6 Estimated Influenza-Associated Excess All-Cause Deaths, in Absolute Number of Cases and per 100,000 People, by Age Group and Epidemic Season in Spain Between 2008/2009 and 2017/2018

The absolute and relative influenza-associated excess all-cause deaths stratified by age groups and by diagnoses groups are displayed in Additional file 1: Table S5.

Influenza economic burden

Hospitalizations with a diagnosis of influenza

Over ten years, the 121,393 hospitalizations with a diagnosis of influenza have costed NHS €422.8 million among the Spanish population. Costs varied with the epidemic season, from €5.9 million in 2008/09 to €120.5 million in 2017/18. The mean annual cost was €45.7 million (excluding the 2009/10 pandemic season), of which 52.8% was generated by patients aged ≥ 65 years1 (Fig. 1). Over the study period, patients with comorbidities, regardless of age, accounted for 65.7% of the costs of hospitalizations with a diagnosis of influenza1. At-risk patients, either because of their age (≥ 65 years) or underlying medical conditions, accounted for 79.9% of influenza hospitalization costs1. Table 2 presents the mean influenza hospitalization cost per patient, according to their age and comorbidities profile.

Fig. 1
figure 1

Direct Cost of Hospitalizations With a Diagnosis of Influenza by Age Group and Epidemic Season in Spanish Public Hospitals Between 2008/2009 and 2017/2018 (Million €)

Influenza-associated excess hospitalization

The direct cost of all-age influenza-associated excess hospitalizations over the ten-year study period was estimated at €484.9 million in the P&I, €826.1 million in the R, €1339.1 million in the C&R, and €1015.2 million in the all-cause model.

The mean annual cost of influenza-associated excess C&R hospitalizations over 9 seasons was estimated at 142.9 M€ (min: 68.6 M€; max: 216.6 M€) for the all ages group, €115.9 million for the population aged ≥ 65 years, €24.2 million for the population aged 50 to 64 years, and €8.8 million for the population aged 19 to 49 years. The population aged ≥ 65 years accounted for 76.4% of the cost, and younger adults accounted for 22.1% of the cost (Fig. 2).

Fig. 2
figure 2

Direct Cost of Influenza-Associated Excess Respiratory or Cardiovascular Hospitalizations by Age Group and Epidemic Season in Spanish Public Hospitals Between 2008/2009 and 2017/2018 (Million €). The total cost corresponds to the cost of all-age influenza-associated excess C&R hospitalization. The percentage for each age group was computed as a proportion of the sum of influenza-associated excess C&R hospitalization estimated for each individual age group

Discussion

This is the first study to estimate the global epidemiologic and economic morbidity and mortality burden of influenza and to capture complications associated with influenza in Spain by integrating two approaches: cases actually coded with influenza-specific diagnosis and cases associated with influenza based on statistical excess modeling techniques. The results were assessed across ten epidemic seasons, six age groups, and four groups of diagnoses. The mortality data represent all deaths registered in Spain. The hospitalization data represent all estimated NHS hospitalizations, exclusive of the private sector.

Findings from this study confirm the high burden of severe influenza in Spain. Furthermore, regardless of the selected group of diagnoses, the excess modeling analyses show that the burden of severe influenza will be underestimated if only administrative data with influenza-specific codes is used (both for hospitalization and mortality).

Mean NHS hospitalizations with a diagnosis of influenza per 100,000 per season was 28.1 for all ages and 88.2 in the population aged ≥ 65 years. These are higher than those reported by San-Román-Montero et al. [14] for seasons 2009/10 to 2014/15 due to the much higher admission rates observed in seasons 2015/16 to 2017/18. If the mean annual influenza-associated excess NHS hospitalizations are considered instead, this ratio would rise to 35.3 cases per 100,000 (95% CI: 32.6–38.0) for all ages and 123.8 (95% CI: 113.4–133.4) for the population aged ≥ 65 years in the P&I; 63.4 (95% CI: 56.0–70.3) and 268.1 (95% CI: 240.4–294.2) in the R; 75.0 (95% CI: 63.3–86.3) and 335.3 (95% CI: 293.2–377.5) in the C&R; and 66.1 (95% CI: 31.3–99.8) and 309.7 (95% CI: 223.7–393.3) in the all-cause model, respectively. Hospitalizations with a diagnosis of influenza during the study period represented only 37.4% of the estimated influenza-associated excess C&R hospitalizations in all-age group and 26.3% in the population aged ≥ 65 years. Czaja et al. reported similar ratios for Colorado, USA, estimating four influenza-associated excess C&R hospitalizations for each influenza-coded hospitalization from hospital discharge records, in the population aged ≥ 65 years (vs. our 3.8 ratio over 9 seasons) [27]. However, the data also suggest an improvement in testing and/or coding for influenza in Spain, as the mean over the last three analyzed seasons increased to 64.8% in all-age and 48.0% in population aged ≥ 65.

The mean annual influenza-associated deaths per 100,000 people were estimated at 2.8 (95% CI: 2.3–3.4) for all ages and 14.4 (95% CI: 11.5–17.3) for population aged ≥ 65 years for P&I; 9.0 (95% CI: 7.2–10.6) and 46.3 (95% CI: 37.1–55.3) for R; 17.8 (95% CI: 14.3–21.4) and 90.9 (95% CI: 72.4–109.8) for C&R; and, 27.7 (95% CI: 21.9–33.5) and 144.4 (95% CI: 112.4–176.9) for all-cause deaths, respectively. Deaths observed during NHS hospitalizations with a diagnosis of influenza represented only 5.3% and 4.3% of the estimated influenza-associated excess all-cause deaths (regardless of place of death, within Spain) in the all ages group and in the population aged ≥ 65 years, respectively.

In both approaches, most of the hospitalizations and deaths were observed in the elderly patients. However, relevant excess hospitalization was also observed in younger age groups, particularly in those aged between 50 and 64 years old. This age group also reported the longest hospital stays and the highest mean hospitalization cost per patient. Beyond age, underlying medical conditions play an important role and must be adequately addressed. However, the population with comorbidities had longer hospital stays, higher hospitalization costs, and higher in-hospital case fatality risk irrespective of the age. The highest mean LoS was observed in patients with comorbidities, aged < 65 years (12.0 days), which had an impact on hospitalization costs but also on lost productivity; the latter was not quantified in this study.

As expected, excess morbidity and mortality estimates varied considerably with the groups of diagnoses, age, and seasons used in the model. Although the performance of each model was computed, the analysis alone cannot determine the importance of all‐cause or cause‐specific data in estimating morbidity and mortality, which was also an issue raised by other authors [42]. Nonetheless, it is extensively studied that the influenza infection can trigger various complications—beyond respiratory complications -, leading to hospitalizations and deaths [3]. Some authors who have computed influenza-associated respiratory deaths for Spain and other geographies have stated that their estimates were conservative, as evidence suggested that excess mortality estimates could potentially double when cardiovascular deaths were included [43, 44]. This duplication was indeed observed in our study between the R and C&R models, both for all ages and for those aged ≥ 65 years.

Other authors have explored the influenza-associated excess mortality in Spain using different groups of diagnoses. In the GLaMOR project, Paget et al. (2022) estimated the average influenza-associated respiratory deaths in Spain per 100,000 people between 2002 and 2011 (excluding 2009) to be 5.4 (2.9–7.3) for all ages and 28.3 (14.4–41.4) for population aged ≥ 65 years, which is lower than our estimates for the same group of diagnoses [43]. In another study, Schmidt et al. estimated the influenza‐associated excess mortality for (i) all‐cause, (ii) respiratory and circulatory, (iii) circulatory, (iv) respiratory, and (v) pneumonia and influenza, for the 2015/2016 and 2016/2017 influenza seasons, in three countries, including Spain [42]. In Spain, the mean excess mortality for the two seasons based on all‐cause, C&R, R, and P&I was 37.9, 19.6, 10.1, and 3.3 per 100,000 people, respectively (all ages) [42]. In all groups of diagnoses, Schmidt et al. estimates are slightly higher but closer to the ones from our study.

These differences reflect the variability of results and difficulty in comparing different studies of this nature [45]. The reasons for the variability in results may come mainly from differences in seasons considered in each study (different circulating strains, level of match between vaccine and circulating strains, vaccine coverage rates) and possibly from differences in the used methods (e.g., statistical models, causes and age groups definitions, data source); however, the accurate reason is yet to be explored [45]. A possible explanation for the lower estimates reported by Paget et al. may be the different viral characteristics across the analyzed seasons. In our study, A(H3N2) was predominant or co-predominant in 66.7% of the seasons [46] versus 44.4% as observed in the GLaMOR study [43, 47]. In another article, the authors report that the dominant subtype circulating in each season was a predictor of excess respiratory mortality [47]. The most severe seasons coincided with influenza A(H3N2) in older adults and it was A(H1N1)pdm09 in the younger group [47]. Methods and data sources also differed with those from Paget et al.

Results of this study should be interpreted in light of several limitations. First, the study is based on extrapolated administrative data, and subject to coding errors or missing information. No laboratory data or information on the vaccination status of subjects were available. Data on deaths per cause per age group from death certificate could be retrieved only on a monthly basis. The study period has included two ICD systems, but ICD10 brought a greater granularity and it may have impacted coding practices. The excess modeling analysis helped overcome some of these limitations. As for the excess hospitalization and mortality model, we used ILI incidence as an indicator of influenza activity, without further control variables, based on published literature [24, 27, 38]. Other variables could have also been used, such as laboratory data on the circulation of other respiratory viruses [48, 49]. In France, a comparison of models using ILI and percent influenza positive demonstrated that ILI was the most statistically relevant indicator of mortality and that, most importantly, both indicators produced similar influenza burden estimates [38]. Nonetheless, as indicated by Froes et al., these findings may not necessarily hold in other countries [24]. The choice to use ILI indicator for influenza activity in the model may be a limitation of our work, as ILI may capture diseases caused by other respiratory pathogens. Hence excess mortality modelled with ILI cannot strictly be attributed to influenza alone. Still, the models employed in the BARI study presented global high performance in terms of correlation and MAPE, as presented in the Supplementary Materials. In the excess hospitalization and mortality models, we used the “all-ages” ILI incidence as an indicator of influenza activity, as reported by the national authorities [31] and without further control variables. The use of age-stratified incidence rates throughout the study period, if made available, should improve the performance of the models per age group.

Cost estimates are based on mean estimated DRG costs at a national level, although distinct DRG costs may be observed across autonomous communities. Also, despite providing a broader perspective on the cost of treating severe influenza, the study still shows only a fraction of the problem. Hospitalizations in the private setting, the burden of milder influenza cases with self-management, or those cases managed in primary health care setting were not included. Indirect costs of lost productivity (from patients or their parents) and the medium and long-term costs related to the posterior treatment of complications caused or aggravated by influenza were not computed. Published studies for Spain report that hospitalization costs represent less than half of the total direct costs of influenza (with its weight varying according to patients’ age and comorbidities status), reflecting the importance of considering the burden in other settings of care [20, 23].

Overall, the present study supports the need for increased testing for influenza in Spain to tackle the current underestimation of influenza burden [27]. The use of a multiplier approach when computing influenza’s burden, as performed in other geographies could be helpful [50]. Further studies analyzing the impact of influenza on specific complications could enable more targeted interventions. Regional assessments of the burden of influenza could also provide further insights for policy evaluation and decision-making.

Conclusion

The BARI study contributed to a better understanding of the burden of influenza in Spain, showing the impact of influenza and its complications on hospitalization and mortality. Also, administrative data without a systematic testing for influenza would lead to under-detection of influenza, thus underestimating the impact of influenza virus. We conclude that there is still a substantial burden of influenza on the elderly population and younger population in Spain, which needs to be addressed.

Availability of data and materials

The data that support the findings of this study are available from IQVIA, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of IQVIA. Those wishing to request the data from this study should contact the author Mafalda Carmo.

Notes

  1. Mean for nine seasons (2008/2009 to 2017/2018), ie, excluding the H1N1pdm09 pandemic (2009/10).

References

  1. Troeger CE, Blacker BF, Khalil IA, Zimsen SRM, Albertson SB, Abate D, Abdela J, Adhikari TB, Aghayan SA, Agrawal S, et al. Mortality, morbidity, and hospitalisations due to influenza lower respiratory tract infections, 2017: an analysis for the Global Burden of Disease Study 2017. Lancet Respir Med. 2019;7(1):69–89.

    Article  Google Scholar 

  2. Paules C, Subbarao K. Influenza. The Lancet. 2017;390(10095):697–708.

    Article  Google Scholar 

  3. Macias AE, McElhaney JE, Chaves SS, Nealon J, Nunes MC, Samson SI, Seet BT, Weinke T, Yu H. The disease burden of influenza beyond respiratory illness. Vaccine. 2021;39(Suppl 1):A6–14.

    Article  CAS  Google Scholar 

  4. Mondiale de la Santé O, Organization WH. Vaccines against influenza: WHO position paper–May 2022–Vaccins antigrippaux: note de synthèse de l’OMS–mai 2022. Weekly Epidemiological Record= Relevé épidémiologique hebdomadaire 2022, 97(19):185–208.

  5. https://www.euro.who.int/en/health-topics/communicable-diseases/influenza/vaccination.

  6. MCSBS: Coberturas de vacunación frente a gripe en ≥65 años, personas de 60–64 años, embarazadas y personal sanitario. Comunidades autónomas. Campaña 2018–2019 (actualización). In. Edited by Dirección General de Salud Pública CeI, Subdirección General de Promoción de la Salud y Vigilancia en Salud Pública; 2019.

  7. MCSBS: Coberturas de vacunación frente a gripe en ≥65 años, personas de 60–64 años, embarazadas y personal sanitario. Comunidades autónomas. Campaña 2020–2021. In. Edited by Dirección General de Salud Pública CeI SGdPdlSyP; 2021.

  8. MCSBS: Evolución de cobertura de vacunación antigripal en población ≥ 65 años. España, temporadas 2008–2009 a 2017–2018 (actualización). Edited by Dirección General de Salud Pública CeI, Subdirección General de Promoción de la Salud y Vigilancia en Salud Pública.; 2018

  9. Bresee J, Hayden F. Epidemic Influenza - Responding to the Expected but Unpredictable. N Engl J Med. 2013;368:89.

    Article  Google Scholar 

  10. Vaccine Effectiveness: How Well Do Flu Vaccines Work? https://www.cdc.gov/flu/vaccines-work/vaccineeffect.htm

  11. Rodrigues E, Machado A, Silva S, Nunes B. Excess pneumonia and influenza hospitalizations associated with influenza epidemics in Portugal from season 1998/1999 to 2014/2015. Influenza Other Respir Viruses. 2018;12(1):153–60.

    Article  Google Scholar 

  12. Hamilton MA, Calzavara A, Emerson SD, Djebli M, Sundaram ME, Chan AK, Kustra R, Baral SD, Mishra S, Kwong JC. Validating International Classification of Disease 10th Revision algorithms for identifying influenza and respiratory syncytial virus hospitalizations. PLoS ONE. 2021;16(1): e0244746.

    Article  CAS  Google Scholar 

  13. Thompson WW, Moore MR, Weintraub E, Cheng P-Y, Jin X, Bridges CB, Bresee JS, Shay DK. Estimating Influenza-Associated Deaths in the United States. Am J Public Health. 2009;99(S2):S225–30.

    Article  Google Scholar 

  14. San-Román-Montero JM, Gil Prieto R, Gallardo Pino C, Hinojosa Mena J, Zapatero Gaviria A. Gil de Miguel A: Inpatient hospital fatality related to coding (ICD-9-CM) of the influenza diagnosis in Spain (2009–2015). BMC Infect Dis. 2019;19(1):700–700.

    Article  Google Scholar 

  15. Pérez-Rubio A, Platero L, Eiros Bouza JM. Gripe estacional en España: carga clínica y económica y programas de vacunación. Med Clin. 2019;153(1):16–27.

    Article  Google Scholar 

  16. Galante M, Garin O, Sicuri E, Cots F, García-Altés A, Ferrer M, Dominguez À, Alonso J. Health services utilization, work absenteeism and costs of pandemic influenza A (H1N1) 2009 in Spain: a multicenter-longitudinal study. PLoS ONE. 2012;7(2):e31696–e31696.

    Article  CAS  Google Scholar 

  17. Badia Llach X, Roset Gamisans M, Francés Tudel JM, Alvarez Sanz C, Rubio Terrés C. Study of flu costs. Aten Primaria. 2006;38(5):260–7.

    Article  Google Scholar 

  18. Albertí C, Orriols R, Manzanera R, Jardí J: Flu and other Acute Respiratory Infections in the Working Population. The Impact of Influenza A (H1N1) Epidemic. Archivos de Bronconeumología ((English Edition)) 2010, 46(12):634–639.

  19. Soto M, Sampietro-Colom L, Vilella A, Pantoja E, Asenjo M, Arjona R, Hurtado JC, Trilla A, Alvarez-Martínez MJ, Mira A, et al. Economic Impact of a New Rapid PCR Assay for Detecting Influenza Virus in an Emergency Department and Hospitalized Patients. PLoS ONE. 2016;11(1):e0146620–e0146620.

    Article  Google Scholar 

  20. Gil de Miguel Á, Eiros Bouza JM, Martínez Alcorta LI, Callejo D, Miñarro C, Vallejo-Aparicio LA, García A, Tafalla M, Cambronero MDR, Rodríguez R, et al. Direct Medical Costs of Four Vaccine-Preventable Infectious Diseases in Older Adults in Spain. Pharmacoeconomics. 2022;6(4):509–18.

    Article  Google Scholar 

  21. León-Gómez I, Delgado-Sanz C, Jiménez-Jorge S, Flores V, Simón F, Gómez-Barroso D, Larrauri A, de Mateo OS. Excess mortality associated with influenza in Spain in winter 2012. Gac Sanit. 2015;29(4):258–65.

    Article  Google Scholar 

  22. López-Cuadrado T, de Mateo S, Jiménez-Jorge S, Savulescu C, Larrauri A. Influenza-related mortality in Spain, 1999–2005. Gac Sanit. 2012;26(4):325–9.

    Article  Google Scholar 

  23. Gil-de-Miguel Á, Martinón-Torres F, Díez-Domingo J, de Lejarazu Leonardo RO, Pumarola T, Carmo M, Drago G, López-Belmonte JL, Bricout H, de Courville C, et al. Clinical and economic burden of physician-diagnosed influenza in adults during the 2017/2018 epidemic season in Spain. BMC Public Health. 2022;22(1):2369.

    Article  Google Scholar 

  24. Froes F, Carmo M, Lopes H, Bizouard G, Gomes C, Martins M, Bricout H, de Courville C, de Sousa JC, Rabaçal C, et al. Excess hospitalizations and mortality associated with seasonal influenza in Portugal, 2008–2018. BMC Infect Dis. 2022;22(1):726.

    Article  CAS  Google Scholar 

  25. Martinón-Torres F, Carmo M, Platero L, Drago G, López-Belmonte J, Bangert M, Díez-Domingo J, Garcés-Sánchez M. Clinical and economic burden of respiratory syncytial virus in Spanish children: the BARI study. BMC Infect Dis. 2022;22(1):1–15.

    Article  Google Scholar 

  26. Bandeira T, Carmo M, Lopes H, Gomes C, Martins M, Guzman C, Bangert M, Rodrigues F, Januário G, Tomé T. Burden and severity of children's hospitalizations by respiratory syncytial virus in Portugal, 2015–2018. Influenza and Other Respiratory Viruses 2022.

  27. Czaja CA, Miller L, Colborn K, Cockburn MG, Alden N, Herlihy RK, Simões EAF. State-level estimates of excess hospitalizations and deaths associated with influenza. Influenza Other Respir Viruses. 2020;14(2):111–21.

    Article  Google Scholar 

  28. https://www.cdc.gov/flu/about/burden/how-cdc-estimates.htm#References.

  29. Lindner L, García-Sánchez R, Alvarez C, Betegón L, Badia X. Hospitalizaciones por hipoglucemia grave en pacientes con diabetes mellitus en España. Rev Clin Esp. 2013;213(8):370–6.

    Article  CAS  Google Scholar 

  30. INE: Estadísticas de Defunciones según la Causa de Muerte, 2008 - 2018. In: Edited by Estadística INd; 2020.

  31. ScVGE. Tasas de incidencia semanal de gripe, temporadas 2007–08 a 2018–19. In: Edited by Epidemiológica. ScdVdlGeESRNdV; 2020.

  32. Decision A. 253/EC Laying Down Case Definitions for Reporting Communicable Diseases to the Community Network Under Decision No 2119/98/EC of the European Parliament and of the Council.In: Official Journal of the European Union 2012.

  33. ISCIII. Guía de procedimientos para la vigilancia de gripe en España, Octubre 2019. In: Edited by Microbiología GdVdGdCNdELdVRyGdCNd; 2019.

  34. Fitzner J, Qasmieh S, Mounts AW, Alexander B, Besselaar T, Briand S, Brown C, Clark S, Dueger E, Gross D. Revision of clinical case definitions: influenza-like illness and severe acute respiratory infection. Bull World Health Organ. 2018;96(2):122.

    Article  Google Scholar 

  35. INE. Población residente por grupo de edad, a 1 enero de cada año, 2008–2019. In: Instituto Nacional de Estadística; 2020.

  36. Thompson WW, Shay DK, Weintraub E, Brammer L, Cox N, Anderson LJ, Fukuda K. Mortality Associated With Influenza and Respiratory Syncytial Virus in the United States. JAMA. 2003;289(2):179.

    Article  Google Scholar 

  37. Smith DJ. Mapping the Antigenic and Genetic Evolution of Influenza Virus. Science. 2004;305(5682):371–6.

    Article  CAS  Google Scholar 

  38. Lemaitre M, Carrat F, Rey G, Miller M, Simonsen L, Viboud C. Mortality burden of the 2009 A/H1N1 influenza pandemic in France: comparison to seasonal influenza and the A/H3N2 pandemic. PLoS ONE. 2012;7(9):e45051–e45051.

    Article  CAS  Google Scholar 

  39. Yang L, Chan KP, Cowling BJ, Chiu SS, Chan KH, Peiris JSM, Wong CM. Excess mortality associated with the 2009 pandemic of influenza A(H1N1) in Hong Kong. Epidemiol Infect. 2011;140(9):1542–50.

    Article  Google Scholar 

  40. García-Eroles L, Illa C, Arias A, Casas M. Los Top 20 2000: objetivos, ventajas y limitaciones del método. Rev Calid Asist. 2001;16(2):107–16.

    Article  Google Scholar 

  41. Sánchez-Martínez F, Abellán-Perpiñán JM, Martínez-Pérez JE, Puig-Junoy J. Cost accounting and public reimbursement schemes in Spanish hospitals. Health Care Manag Sci. 2006;9(3):225–32.

    Article  Google Scholar 

  42. Schmidt SSS, Iuliano AD, Vestergaard LS, Mazagatos-Ateca C, Larrauri A, Brauner JM, Olsen SJ, Nielsen J, Salomon JA, Krause TG. All-cause versus cause-specific excess deaths for estimating influenza-associated mortality in Denmark, Spain, and the United States. Influenza Other Respir Viruses. 2022;16(4):707–16.

    Article  CAS  Google Scholar 

  43. Paget J, Danielle Iuliano A, Taylor RJ, Simonsen L, Viboud C, Spreeuwenberg P. Estimates of mortality associated with seasonal influenza for the European Union from the GLaMOR project. Vaccine. 2022;40(9):1361–9.

    Article  Google Scholar 

  44. Iuliano AD, Roguski KM, Chang HH, Muscatello DJ, Palekar R, Tempia S, Cohen C, Gran JM, Schanzer D, Cowling BJ, et al. Estimates of global seasonal influenza-associated respiratory mortality: a modelling study. Lancet (London, England). 2018;391(10127):1285–300.

    Article  Google Scholar 

  45. Roguski KM, Rolfes MA, Reich JS, Owens Z, Patel N, Fitzner J, Cozza V, Lafond KE, Azziz-Baumgartner E, Iuliano AD. Variability in published rates of influenza-associated hospitalizations: A systematic review, 2007–2018. J Glob Health. 2020;10(2): 020430.

    Article  Google Scholar 

  46. III IdSC. Informes de Vigilancia de la Gripe en España, temporadas 2008–2009 hasta la temporada 2017–2018. In: Edited by España SdVdlGe; 2008–2018.

  47. Paget J, Spreeuwenberg P, Charu V, Taylor RJ, Iuliano AD, Bresee J, Simonsen L, Viboud C. Global mortality associated with seasonal influenza epidemics: New burden estimates and predictors from the GLaMOR Project. J Glob Health. 2019;9(2): 020421.

    Article  Google Scholar 

  48. Leung VK, Wong JY, Barnes R, Kelso J, Milne GJ, Blyth CC, Cowling BJ, Moore HC, Sullivan SG. Excess respiratory mortality and hospitalizations associated with influenza in Australia, 2007–2015. Int J Epidemiol. 2022;51(2):458–67.

    Article  Google Scholar 

  49. Nielsen J, Krause TG, Mølbak K. Influenza-associated mortality determined from all-cause mortality, Denmark 2010/11-2016/17: The FluMOMO model. Influenza Other Respir Viruses. 2018;12(5):591–604.

    Article  Google Scholar 

  50. Influenza Hospitalization Surveillance Network (FluSurv-NET). https://www.cdc.gov/flu/weekly/influenza-hospitalization-surveillance.htm.

Download references

Acknowledgements

The authors would like to thank all investigators who participated in the BARI study, and José María Guillén for assistance during the project design and results assessment.

Funding

The BARI study was funded by Sanofi. Martinón-Torres F has received support for the present work from the Instituto de Salud Carlos III (Proyecto de Investigación en Salud, Acción Estratégica en Salud): Fondo de Investigación Sanitaria (FIS; PI070069/PI1000540/PI1601569/PI1901090) del plan nacional de I + D + I and ‘fondos FEDER’ and Proyectos GAIN Rescata-Covid_IN845D 2020/23 (GAIN, Xunta de Galicia).

Author information

Authors and Affiliations

Authors

Contributions

All authors have contributed to the protocol revision, critical analysis, and discussion of the study results, and collaborated for the redaction and revision of the manuscript. MC and GB have conducted the technical analysis of the BARI study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to T. Pumarola.

Ethics declarations

Ethics approval and consent to participate

The study was conducted following the ethical principles of the Declaration of Helsinki and the local regulation, including privacy laws. The protocol of the BARI study was validated by a panel of clinical experts, classified by the Agency of Medicines and Medical Devices (AEMPS) as an observational study and approved by the Ethics Committee of Hospital Clinic de Barcelona (HCB/2020/1132). Informed consent to participate was waived by the same ethics committee that approved the study (Ethics Committee of Hospital Clinic de Barcelona).

Consent for publication

Not applicable.

Competing interests

Drago G, López-Belmonte JL, Bricout H, and de Courville C are the employees of Sanofi and may hold shares and/or stock options in the company. Carmo M, Bizouard G are the employees of IQVIA. Gil-de-Miguel A and Díez-Domingo J have received fees from Sanofi. Martinón-Torres F received honoraria from GSK group of companies, Pfizer Inc, Sanofi, MSD, Seqirus, and Janssen for taking part in advisory boards and expert meetings and for acting as a speaker in congresses outside the scope of the submitted work. Martinón-Torres F has also acted as the principal investigator in randomized controlled trials of the above-mentioned companies as well as Ablynx, Gilead, Regeneron, Roche, Abbott, Novavax, and MedImmune, with honoraria paid to his institution. Ortiz de Lejarazu R received honoraria from AstraZeneca, GSK, MSD, Roche, Sanofi, and Seqirus for taking part in advisory boards and expert meetings and for acting as a speaker in congresses outside the scope of the submitted work. Ortiz de Lejarazu R has been the Director of the National Influenza Centre, Valladolid, Spain, until 2019, without any salary for that position and has participated in regular commitments of influenza and other respiratory virus surveillance. Pumarola T received honoraria from Roche, Sanofi, Seegene, and Seqirus for taking part in advisory boards and expert meetings and for acting as a speaker in congresses outside the scope of the submitted work. Redondo Margüello E received honoraria from GSK group of companies, Pfizer Inc, Sanofi, and MSD, for taking part in advisory boards and expert meetings and for acting as a speaker in congresses outside the scope of the submitted work.

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. Table S1.

Diagnostic codes used to identify comorbidities/risk factors for influenza: a) in people ≥ 5 years old; b) in children < 5 years old. Table S2. Performance of the excess hospitalization model per age group and cause. Table S3. Performance of the excess deaths model per age group and cause. Table S4. Estimated influenza-associated excess pneumonia or influenza, respiratory, and all-cause excess hospitalizations, in absolute and per 100,000 people, by age group and epidemic season in Spanish public hospitals between 2008/2009 and 2017/2018. Table S5. Estimated influenza-associated excess pneumonia or influenza, respiratory, and respiratory or cardiovascular deaths, in absolute and per 100,000 people, by age group and epidemic season in Spain between 2008/2009 and 2017/2018. Table S6. P-value associated to the ILI variables of estimated influenza-associated excess hospitalization. Table S7. P-value associated to the ILI variables of estimated influenza-associated excess deaths.

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

Pumarola, T., Díez-Domingo, J., Martinón-Torres, F. et al. Excess hospitalizations and mortality associated with seasonal influenza in Spain, 2008–2018. BMC Infect Dis 23, 86 (2023). https://doi.org/10.1186/s12879-023-08015-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-023-08015-3

Keywords