Estimating Influenza-associated Mortality in Korea: The 2009-2016 Seasons
Article information
Abstract
Objectives
Estimating influenza-associated mortality is important since seasonal influenza affects persons of all ages, causing severe illness or death. This study aimed to estimate influenza-associated mortality, considering both periodic changes and age-specific mortality by influenza subtypes.
Methods
Using the Microdata Integrated Service from Statistics Korea, we collected weekly mortality data including cause of death. Laboratory surveillance data of respiratory viruses from 2009 to 2016 were obtained from the Korea Centers for Disease Control and Prevention. After adjusting for the annual age-specific population size, we used a negative binomial regression model by age group and influenza subtype.
Results
Overall, 1 859 890 deaths were observed and the average rate of influenza virus positivity was 14.7% (standard deviation [SD], 5.8), with the following subtype distribution: A(H1N1), 5.0% (SD, 5.8); A(H3N2), 4.4% (SD, 3.4); and B, 5.3% (SD, 3.7). As a result, among individuals under 65 years old, 6774 (0.51%) all-cause deaths, 2521 (3.05%) respiratory or circulatory deaths, and 1048 (18.23%) influenza or pneumonia deaths were estimated. Among those 65 years of age or older, 30 414 (2.27%) all-cause deaths, 16 411 (3.42%) respiratory or circulatory deaths, and 4906 (6.87%) influenza or pneumonia deaths were estimated. Influenza A(H3N2) virus was the major contributor to influenza-associated all-cause and respiratory or circulatory deaths in both age groups. However, influenza A(H1N1) virus–associated influenza or pneumonia deaths were more common in those under 65 years old.
Conclusions
Influenza-associated mortality was substantial during this period, especially in the elderly. By subtype, influenza A(H3N2) virus made the largest contribution to influenza-associated mortality.
INTRODUCTION
Influenza is a common acute respiratory disease caused by influenza virus. Although influenza is frequently confused with the common cold, some individuals with influenza can become severely ill, to the point of hospitalization or even death. According to estimates of the World Health Organization, 5-10% of adults and 20-30% of children are infected by the influenza virus every year [1-3]. Therefore, measuring the disease burden of seasonal influenza is of major interest among infection epidemiologists [4-6].
A simple way to measure the disease burden is to estimate the mortality of influenza by comparing the mortality rate to the rate of influenza virus positivity. However, because of the seasonality of influenza virus and the annual variability of subtypes that cause epidemics [7], it is necessary to include complicated calculations on time and the positivity rate of each virus subtype to obtain accurate estimates of influenza-associated mortality. Furthermore, when we count only people who die directly from influenza, we cannot estimate the substantial impact that influenza virus may exert by inducing severe complications that lead to death [7,8]. Additionally, we need to explore various causes of death when estimating influenza-associated deaths because records for mortality may not directly include information about viral infection. Therefore, few studies have considered all these factors when measuring the disease burden of influenza. This study aimed to estimate influenza-associated deaths considering influenza type, patients’ age, and various causes of death, and to explore differences in influenza severity by age group through these results. We used a negative binomial regression model for the comprehensive estimation, following previous studies [9-12].
METHODS
Data Sources and Measures
Mortality data were obtained from 2009 through 2016 for all deaths in Korea. We extracted age-specific mortality from the Microdata Integrated System (MDIS) [13], generated by Statistics Korea. The causes of death were coded by the International Statistical Classification of Diseases and Related Health Problems 10th revision in this dataset and included demographic features such as age, sex, and residence. We constructed three different models by cause of death: all-cause deaths, respiratory (diseases of the respiratory system, J00-J99) or circulatory (diseases of the circulatory system, I00-I99) deaths, and influenza or pneumonia (influenza and pneumonia, J09-J18) deaths. We also selected other chronic causes of deaths such as cancer (neoplasms, C00-D48) to compare the results. Age-specific weekly deaths were divided by age-specific mid-year population estimates, which were also obtained from the MDIS. The weekly age-specific mortality rate was used as a predictive variable.
The Korea Centers for Disease Control and Prevention (KCDC) provided influenza virus surveillance data through the Korea Influenza and Respiratory Viruses Surveillance System (KINRESS) [14]. Thirty-six to 90 representative institutes reported tested and positive samples for the following subtypes of influenza virus: influenza A(H1N1), influenza A(H3N2), and influenza B virus. Since the KCDC started separate monitoring of each influenza virus type in 2009, we used data from 2009 to 2016 in our models. Since the influenza surveillance system is organized by season globally, we structured each season to extend from the 36th week of one year to the 35th week of the next year.
Statistical Analysis
The model included a total of 7 seasons, from 2009 to 2016, including the 2009 H1N1 pandemic. We used negative binomial regression models for the weekly age-specific mortality data. Various statistical methods have been developed for estimating influenza-associated deaths [15-30], but there is no single way to estimate the exact scale of deaths. Poisson regression analysis is broadly accepted as the predominant estimation method. However, negative binomial regression is also commonly used for modeling count data, especially when there is a suspicion of overdispersion. Since the Poisson distribution does not allow greater variability, the distribution of age-specific deaths and the positivity rate of each influenza virus in this study needed additional free parameters for model fitness. Therefore, to estimate influenza-associated deaths in 2009-2016, we used a negative binomial regression model [31] considering influenza subtypes, age, and various causes of death. An additional process of checking the model assumptions and the goodness-of-fit was performed. Mortality data were transformed into weekly deaths and divided into 2 age groups: under 65 years old and 65 years or older. To consider linear and non-linear time trends, we included these factors as covariates and the age-specific population as the offset. The prediction model can be expressed as:
where Y represents the predicted number of deaths, α is the offset term for the age-specific population, and β6 through β8 represent the estimated coefficient for each viral subtype. After predicting age-specific mortality by cause of death, we calculated viral subtype–specific deaths as the count data using this model. The total number of influenza-associated deaths was expressed as the sum of the count for each viral subtype. To deal with weekly variance and autocorrelation, the Durbin-Watson test was performed and the Akaike information criterion was compared. Finally, we inserted a 1-week lag for mortality in each model. All statistical analyses were conducted in R version 3.4.0 (https://cran.r-project.org/bin/windows/base/old/3.4.0/).
Ethics Statement
All data analyzed in this study are available from each organization with permission. Since no personal data were collected or stored, there was no personal information. Code for re-analysis is available from the first author upon request.
RESULTS
From the 2009-2010 season through the 2015-2016 season, 1 859 890 all-cause deaths, 566 475 underlying respiratory or circulatory deaths, and 77 850 underlying influenza or pneumonia deaths were observed (Table 1). Deaths among those at least 65 years of age accounted for 72.3% (all-cause deaths), 85.3% (respiratory or circulatory deaths), and 92.5% (influenza or pneumonia deaths) of deaths.
In the respiratory virus surveillance data, an average of 14.6% of specimens tested positive for influenza virus, with the following subtype distribution: influenza A(H1N1) virus, 5.0%; influenza A(H3N2) virus, 4.3%; and influenza B virus, 5.3% (Table 2). Influenza A(H1N1) virus was the major cause of epidemics in the 2009-2010 and 2010-2011 seasons, while influenza A(H3N2) virus was the major cause of epidemics in the 2011-2012, 2012- 2013 and 2014-2015 seasons. Influenza B virus was the main cause of the epidemic in the 2013-2014 season. In the 2015-2016 season, both influenza A(H1N1) and influenza B virus were epidemic.
The overall estimated influenza-associated all-cause deaths were 5313 per year (2.02% of all-cause deaths, 10.59 per 100 000 population). In particular, the population 65 years of age or older made the largest contribution to this mortality (2.28% of all-cause deaths among those at least 65 years old, 74.12 per 100 000 population). The estimated influenza-associated mortality figures by age group are summarized in Table 3 and Figure 1. Of the all-cause deaths, 6774 (0.51%) deaths in those under 65 years old and 30 414 deaths (2.27%) in those at least 65 years old were estimated to be influenza-associated deaths. While influenza A(H3N2) and influenza B virus contributed to deaths in the younger population (under 65 years old) at a similar rate (0.50% and 0.48% of all-cause deaths each), influenza A(H3N2) virus made the largest contribution of the 3 virus subtypes to deaths in the older population (at least 65 years old, 1.50%). In total, 2521 (3.05%) and 16 411 (3.42%) respiratory or circulatory deaths were estimated as influenza-associated deaths in the age groups of under 65 years and 65 years or older, respectively. The overall influenza-associated respiratory or circulatory deaths were not significantly different by age group. Similarly to all-cause deaths, influenza A(H3N2) virus made the largest contribution to respiratory or circulatory deaths among those at least 65 years old (2.27%). Furthermore, 1048 (18.23%) and 4906 (6.87%) influenza or pneumonia deaths were estimated to have occurred in the age groups of under 65 years and 65 years or older, respectively, and 10.05% of the influenza or pneumonia deaths in individuals under 65 years old were due to influenza A(H1N1) virus. In addition to season, age group and the cause of death contributed to differences in estimated deaths by influenza subtype. We only considered 3 causes of death and 3 influenza subtypes since no other causes of deaths or respiratory viruses were fitted to our model.
Estimated influenza-associated deaths per 100 000 people by influenza virus type and age. A, C and E shows estimated all-cause, respiratory or circulatory, and influenza or pneumonia deaths by influenza virus in under 65 years old group, respectively. B, D and F shows estimated all-cause, respiratory or circulatory, and influenza or pneumonia deaths by influenza virus in over 65 years old group, respectively. All estimated deaths are visualized by weekly estimated data.
DISCUSSION
There were approximately 5313 influenza-associated deaths per year (10.59 per 100 000 people, 74.12 per 100 000 people above 65 years old). This is the first study, to our knowledge, to estimate age-specific influenza-associated deaths in Korea using a negative binomial regression model and considering lag time and model fit. By adding the estimated mortality of each influenza subtype, we were able to apply the model more accurately. We also considered other respiratory viruses as covariates in the model, but they were excluded from the final model due to statistical non-significance. Moreover, we classified the population by age group and evaluated the suitability of the model, making it possible to estimate excess mortality by age group.
The results of this study can be compared to those of a prior study that estimated 2900 influenza-associated deaths per year in Korea [23]. Park et al. [23] used multiple linear models to estimate the excess mortality rate caused by influenza from 2003 through 2013. However, since it is inappropriate to assume fixed fatality among viral subtypes and age groups, and considering that death is a rare and non-negative event, a count data model might be more persuasive than a linear model for estimating influenza-associated deaths. Furthermore, the KCDC reorganized the respiratory virus surveillance system, KINRESS, in 2006, changing the institutions to be monitored by specifying the percentage of target organizations in each region. They also changed the types of respiratory viruses to be monitored, and the laboratory method for detecting positivity. Therefore, influenza positivity rates before and after 2006 are non-comparable. Additionally, the results of this study are consistent with those reported in studies from other countries [5,16,24,26,28,29], most of which reported 10 or more estimated deaths per 100 000 people, with recent studies showing much higher numbers of estimated deaths among the elderly.
In our analysis of estimated influenza-associated deaths by type of influenza, influenza A(H3N2) virus caused the most excess mortality, accounting for 65.8% of influenza-associated all-cause deaths. This is a meaningful point to focus on because this study included the 2009-2010 season, when influenza A(H1N1) virus became pandemic worldwide. Based on this fact, it can be inferred that the infectivity of influenza A(H1N1) virus was high, but the mortality rate was low, or that the later influenza viruses’ mortality rate was higher than that of influenza A(H1N1) virus. Regardless of the study period, a higher mortality rate of influenza A(H3N2) virus has also been proven in other countries, perhaps because of this subtype’s severity in vulnerable populations [8,23,31-34].
Another noteworthy finding of this study is the difference in the proportion of deaths attributable to different influenza subtypes between age groups. According to estimates of influenza-associated mortality, older individuals (65 years or older) were more vulnerable to influenza A(H3N2) virus. Considering that the deaths of individuals 65 years of age and older account for 85.3% of respiratory or circulatory deaths, the proportion of estimated influenza-associated respiratory or circulatory deaths in the older population is not negligible. Since the mortality data were extracted from death certificates, influenza or pneumonia deaths may have been over-estimated, especially in the 2009-2010 pandemic season. Furthermore, it cannot be guaranteed that other chronic causes of deaths made a constant contribution to weekly all-cause deaths. Therefore, as supported by previous studies [8,15,35,36], estimated all-cause deaths or influenza or pneumonia deaths could have been overestimated or underestimated. Therefore, it can be assumed that respiratory or circulatory deaths provide the most accurate estimate, compared to the other two types of estimates, in the older population.
There are several limitations of this study. The younger age group could not be subdivided into more age groups since we used weekly mortality data for the estimation, so that there were many 0 cells in the under-15 age group. This could lead to inaccurate estimation when using a negative binomial regression model as the estimation method. Therefore, we estimated excess mortality in children and adults under 65 years of age within a single group. Fortunately, this methodological choice could only lead to underestimation of influenza-associated mortality in young adults because influenza or pneumonia deaths and respiratory or circulatory deaths occur more frequently in adults than children. Secondly, since we only used respiratory virus surveillance data for the analysis, it was not possible to consider the effects of other viral or bacterial infections in a particular week. If another infectious disease were epidemic during these seasons, it could have caused an overestimation, but there were no other life-threatening pandemic infectious diseases in this period. Finally, mortality data can have inherent uncertainty because if cause-of-death certificates are filled out at a busy moment, doctors may ignore less important or indirect causes. For this reason, respiratory or circulatory deaths may be more accurate for estimating actual influenza-associated deaths.
The high overall rate of influenza-associated deaths indicates that the burden of influenza is never light, and as noticed in the USA [37] and Europe [38], we may use surveillance data for calculating influenza-associated deaths in real time for the management of high-risk groups in each epidemic. The difference in the severity of each virus type in specific age groups also has important implications. Among the 3 types of estimated deaths in this study, the older (65 years and above) population had higher mortality in all-cause deaths or respiratory or circulatory deaths, particularly with influenza A(H3N2) virus. However, the younger (under 65 years) population had a higher mortality due to influenza A(H1N1) virus, except for all-cause deaths. Additionally, influenza B virus should be considered more cautiously in them than in the elderly. Further studies are needed to explain this variance in influenza epidemics.
Notes
CONFLICT OF INTEREST
The authors have no conflicts of interest associated with the material presented in this paper.
AUTHOR CONTRIBUTIONS
Conceptualization: KH, BCC. Data curation: KH, BCC, SS. Formal analysis: KH. Funding acquisition: None. Methodology: KH, BCC, SS. Project administration: BCC. Visualization: KH. Writing - original draft: KH, BCC, SS. Writing - review & editing: KH, BCC, SS.
ACKNOWLEDGEMENTS
None.