Next Article in Journal
Nitrogen Accumulation in Oyster (Crassostrea gigas) Slurry Exposed to Virucidal Cold Atmospheric Plasma Treatment
Previous Article in Journal
Imaging of Intracellular and Plasma Membrane Pools of PI(4,5)P2 and PI4P in Human Platelets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relative Risk Prediction of Norovirus Incidence under Climate Change in Korea

1
Department of Mathematics and Statistics, Chonnam National University, Gwangju 61186, Korea
2
Microbial Safety Division, National Institute of Agricultural Sciences, Rural Development Administration, Wanju 55365, Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Life 2021, 11(12), 1332; https://doi.org/10.3390/life11121332
Submission received: 28 October 2021 / Revised: 19 November 2021 / Accepted: 30 November 2021 / Published: 2 December 2021
(This article belongs to the Topic Climate Change and Infectious Disease)

Abstract

:
As incidences of food poisoning, especially norovirus-induced diarrhea, are associated with climate change, there is a need for an approach that can be used to predict the risks of such illnesses with high accuracy. In this paper, we predict the winter norovirus incidence rate in Korea compared to that of other diarrhea-causing viruses using a model based on B-spline added to logistic regression to estimate the long-term pattern of illness. We also develop a risk index based on the estimated probability of occurrence. Our probabilistic analysis shows that the risk of norovirus-related food poisoning in winter will remain stable or increase in Korea based on various Representative Concentration Pathway (RCP) scenarios. Our approach can be used to obtain an overview of the changes occurring in regional and seasonal norovirus patterns that can help assist in making appropriate policy decisions.

1. Introduction

As the reproduction of microorganisms is substantially affected by weather-related factors such as temperature and humidity, climate change is extremely likely to cause changes in the seasonal patterns of food poisoning incidents caused by microorganisms. Therefore, as climate conditions continue to change, it is expected that adjustments to food safety policies related to food will become necessary, and scientific methodologies that can accurately predict long-term fluctuations in food poisoning patterns must be developed to assist administrators in proactively devising reasonable policies to respond to the effects of climate change. Most advanced countries already have access to vast amounts of weather and food poisoning-related clinical data; the key challenge is deriving useful information from these big data. According to a recent 2016 study, more than 40 million disability-adjusted life years were lost in children under the age of five due to cases of diarrhea and related deaths [1]. Although the incidence of diarrhea may be reduced by economic development, some reports have found that climate change damages urban infrastructure and reduces overall water availability [2]. The spread of diarrhea is complex and depends on many factors. Infectious diarrhea can be caused by a variety of pathogens, and it is affected by both host susceptibility and environmental components. Therefore, there is an urgent need for research that quantifies the impact of diarrheal illness as well as improved predictive estimates [3,4,5,6]. This paper chooses to consider norovirus specifically for two main reasons. First, variations in the incidence pattern of norovirus are clinically significant. As diagnostic technology continues advancing, norovirus has been increasingly identified as a major cause of food poisoning [7]. The World Health Organization (WHO) reports that norovirus is one of the top five causes of death from food poisoning [8]. More than 200,000 people die from norovirus every year, including 70,000 children in developing countries [9]. In the United States, norovirus causes 23 million cases of acute gastroenteritis (AGE) every year [10]. According to the Centers for Disease Control and Prevention (CDC), norovirus causes annual economic losses of more than $2 billion in the United States alone. In other words, norovirus continues to have a negative impact on both human life and the economy worldwide. Second, norovirus has a distinct seasonal pattern, with the frequency of cases rising in the winter [11]. Because of this, it is also called “the winter vomiting disease” [12,13]. New norovirus variants often emerge on cruise ships in the winter and early spring. In the United States, there were particularly high numbers of deaths from norovirus in the winter months in both 2002–2003 and 2006–2007. These seasonal patterns in disease prevalence make it easier to analyze the relationship between weather-related factors and food poisoning using probabilistic methods. As described above, the main causes of norovirus are contaminated food and human contact; however, many papers examining the relationship between norovirus incidence and climate support the assertion of this paper. Previous studies have shown that norovirus is negatively correlated with mean temperature [14] and that it is a wintertime phenomenon, at least in the temperate northern hemisphere [15]. Another study reviewing norovirus incidence over the previous decade showed differences in case numbers between influenza and norovirus infections, with norovirus still showing strong incidence in certain high population density areas [16]. A recent study examined temperature and relaxation rates among patients with diarrhea given the climate change scenarios in Japan. That study estimated the future probability of the disease based on a simple comparison by region [17]. The above study explains the importance of temperature-related future changes in diarrheal patients. To simplify the model and actively identify temperature-related effects, among various weather-related factors, the current study focused on average daily temperature among various weather-related factors. In this paper, we also attempted to predict the variations in norovirus incidence that will occur as the temperature continues to rise in winter due to climate change.
We introduce an index that we developed to quantify the incidence rate of norovirus-induced diarrhea according to Representative Concentration Pathway (RCP) scenarios that explain climate change, and we used this index to calculate the future risk of disease due to climate change. In this analysis, a generalized additive linear model (GALM) method using B-spline was used to compensate for the fact that the probability of dependence on the explanatory variable increases or decreases in simple logistic regression. In addition, the risk assessment was made more objective by using the relative risk index (RRI) to represent risk.

2. RCP Scenarios

In order to establish an effective response and adaptation plan for risk factors that will appear due to climate change, the Intergovernmental Panel on Climate Change (IPCC) introduced a set of climate change scenarios. The RCP scenarios consider both recent observations along with greenhouse gas reduction technologies. The IPCC has recommended that researchers use the RCP scenarios as basic data for studies that involve climate change impact evaluations; in particular, four scenarios have been introduced: RCP2.6, RCP4.5, RCP6.0, and RCP8.5. The numbers in the scenario names indicate the increasing amounts of net radiation energy experienced by the year 2100 [18].
The Korea Meteorological Administration (KMA) has contributed to the RCP scenarios (http://www.climate.go.kr (accessed on 5 August 2021)). The KMA developed RCP scenarios using five regional climate models (HadGEM3-RA, RegCM4, SNURCM, GRIMs, WRF) [19,20,21,22,23]. As climate models contain uncertainty associated with several sources, e.g., imprecise initial conditions, some statistical refinement techniques and bias correction methods are applied to the output of the model, before providing data to the end users [24]. From these five models they created two model ensembles. The first model is MME4s, based on data from HadGEM3-RA, RegCM4, SNURCM, and WRF. The second one is MME5s, based on data from HadGEM3-RA, RegCM4, SNURCM, GRIMs, and WRF. Climate projections until the year 2100 are evaluated under RCP2.6 and 6.0 for MME4s, and under RCP4.5 and 8.5 for MME5s. The latitude and longitude of the RCP scenario data range from 33° N to 39° N and 124.5° E to 132.0° E, respectively, which includes the entirety of the southern Korean Peninsula.

3. Data Descriptions

From 2005 to 2018, the total number of patients in Korea reported to have diarrhea was 283,651. Of these cases, 24,642 were caused by norovirus, representing a rate of 8.69%, thus making norovirus a very significant cause of diarrhea. These clinical data come from the Korea Centers for Disease Control and Prevention Portal (http://www.kdca.go.kr/ (accessed on 3 September 2020)). Table 1 summarizes the known norovirus outbreaks by gender and age and the rate of norovirus diarrhea cases compared to diarrhea cases from all-causes. The incidence rates of norovirus do not significantly differ by gender, with the rates in males and females being 8.82% and 8.53%, respectively. In contrast, the effect of age is clearly significant. In particular, the incidence rates of norovirus in the age groups of 0–5, 6–15, 16–59, and over 59 are 13%, 3.37%, 4.79%, and 8.74%, respectively. This shows that very young people are significantly more vulnerable to norovirus than those in other age groups. When dividing Korea into northern, central, and southern regions, there were slight differences in the rates of norovirus outbreaks between regions, with rates of 10.14%, 6.05%, and 9.26%, respectively. In the southern part of the country along the coast, diet and temperature are expected to affect the risk of norovirus outbreaks. The northern part of the country also has a somewhat high risk as it includes Seoul, which is particularly densely populated, leading to higher disease risk. The results of the logistic regression analysis used to calculate the odds ratio by gender, age, and region for norovirus occurrence are shown in Table 2. The odds ratio for a female (vs. male) is 0.101, which is not significant. Age had a significant effect: compared to those aged 0–5, the odds ratios were 0.646, 0.350, and 0.245 for individuals aged 6–15, 16–59, and 60 and older, respectively. This indicates that individuals under five years of age are substantially more vulnerable to norovirus infection than those in other age groups. In addition, the odds ratio for the central and southern regions of the country relative to the northern region were significantly different at 0.777 and 1.062, respectively. Figure 1 shows the regions described in Table 2 on a map for better understanding. We used R 4.2.1 for all data analysis [25,26,27,28,29,30].

4. Methods

First, the relationship between weather-related factors and the occurrence of food poisoning caused by norovirus was analyzed stochastically. Based on the results of that analysis, we calculated how the seasonal pattern of norovirus incidence is expected to change in the future due to climate change. In Korea, incidences of diarrhea are reported daily by region and cause. We used average daily temperature data from the Meteorological Data Open Portal (https://data.kma.go.kr/ (accessed on 5 August 2021)). The number of diarrhea patients with norovirus may vary as winter temperatures rise rapidly in accordance with the climate projections under the considered RCP scenarios.
In addition, a generalized additional model (GAM) was employed to express the additive model of the functions including these bases. The logistic model using GAM helps us to predict the probability of the occurrence of norovirus according to the average temperature. In the GAMs that were used, the linear components β j x i j   of the model were replaced with β j f j ( x i j ) [31,32,33,34]. π ( y i | x i 1 , x i 2 , , x i p ) is used to denote the probability that norovirus cases appear along with the explanatory variables x i 1 , x i 2 , , x i p . The GALM assumes that
l o g [ π ( y i | x i 1 , x i 2 , , x i p ) 1 π ( y i | x i 1 , x i 2 , , x i p ) ] = β 0 + β 1 f 1 ( x i 1 ) + β 1 f 1 ( x i 1 ) + + β p f p ( x i p )  
where we used the functions f 1 , f 2 , , f p to smooth bases estimated in various forms, such as with B-spline, cubic splines, and natural cubic splines. We used recursive 4th B-spline bases as supports of the mean temperature as follows [34,35,36,37,38]. Let U be a set of m + 1 non-decreasing numbers, u 0 u 1 u 2 u m .
N i , 0 ( u ) = { 1   i f   u i u u i + 1 0   otherwise  
N i , p ( u ) = u u i u i + p u i N i , p 1 ( u ) + u i + p + 1 u u i + p + 1 u i + 1 N i + 1 , p 1 ( u )
We create a logistic regression model using the previously generated N i , p as f p ( x i p ) and the presence or absence of norovirus as the dependent variable. A backfitting algorithm is typically used for estimations when the splines are complex. However, in this paper, the splines are calculated and then used to estimate the logistic regression because they can be calculated with relative ease. We calculate the binary logistic regression model using iteratively reweighted least squares (IRLS), which is equivalent to maximizing the log-likelihood of a Bernoulli distributed process using Newton’s method as follows.
parameters :   w T = [ β 0 , β 1 , β 2 , ]
explanatory   variables :   x ( i ) = [ β 0 , β 1 , β 2 , ] T
X = [ 1 x 1 ( 1 ) x 2 ( 1 ) 1 x 1 ( 2 ) x 2 ( 2 )   ]
response   variables :   y ( i ) = [ y ( 1 ) , y ( 2 ) , ] T
expected   value   of   the   Bernoulli   distribution :   μ ( i ) = 1 1 + e w T x ( i )
μ = [ μ ( 1 ) , μ ( 2 ) , ]
using iterative algorithm
w k + 1 = ( X T S K X ) 1 X T ( S k X w k + y μ k )
The average temperature that maximizes the probability of GALM is defined as the maximum rate temperature, and the intervals of height corresponding to 0.8 and 0.9 times the maximum rate are respectively defined as the risk interval and the high-risk interval.
We employed the RRI to understand the effect of temperature-dependent probabilities in the RCP scenarios. RRI measures how the probability of occurrence is driven by temperature in GALM models.
Relative   risk   index   ( RRI ) = y i A D e v i a t i o n i × χ ( D e v i a t i o n i )
where A is defined as the focus set, D e v i a t i o n i is defined as π ( y i | · ) —the critical value, and χ represents the indicator function.
Larger RRIs indicate an increased risk of infection by norovirus. That is, the probability in a particular group is greater than the mean probability of being infected by norovirus among diarrhea patients. For example, the incidence of norovirus among all diarrhea patients aged 0–5 is 13.00%, so the critical value is 0.13. We can compute the RRI for age, region, and X-year RCP scenarios.
We calculated the RRI in each RCP scenario for the winter months in Korea (December to February). We organized the 451,351 grids in the RCP scenarios into 96,172 grids based on their location in the Korean peninsula to calculate the RRIs, and then we used the GALM model to predict the probability of the incidence of norovirus on a daily basis. Figure 2 shows the grid, made up of 96,172 points, we used. In a logistic regression model the outcome of y i is 0 or 1, with 1 indicating an event occurring (such as becoming infected with a disease) and 0 indicating no event occurring. In a typical logistic model, the probability of a response increases or decreases depending on the value of the explanatory variable. Therefore, the quadratic form, which increases and decreases in order to obtain an optimum value, is not applicable here. To compensate for this, a number of orthogonal bases were used in this paper to make the quadratic form feasible.

5. Results

Figure 3 shows the GALM analysis using B-spline. The GALM model is more useful for predicting the probability of occurrence than the generalized linear model. The findings indicate that the rate of occurrence in patients in different age groups differs significantly depending on temperature. From the analysis of Figure 3 it is evident that people under the age of 15 are more than twice as likely to contract norovirus compared to those over the age of 15. Especially, the two age groups containing patients under the age of 15 show the highest probability of contracting norovirus when the temperature is near 0 degrees, while the other two age groups, 16–59 years and over 60 years, have the highest probability of being infected by a norovirus when the temperature is between −10 and 0. In other words, the temperature interval in which the highest probability of norovirus occurs was found to be lower for those over 15 years of age compared to those under 15 years of age.
Table 3 lists the maximum incidence rates, maximum rate temperatures, risk intervals, and high-risk intervals of the four age groups. Table 3 shows the temperatures at which patients in each group are vulnerable to norovirus; this information could prove very helpful in preventing norovirus. As noted in the previous section, the proportion of all diarrhea patients aged 0–5 who contract norovirus is 13%, which is significantly higher than the corresponding proportions in all other age groups. The predicted norovirus infection rate among patients aged 0–5 years reaches a maximum of 23.8% at −2.0 °C, which is significantly greater than the overall observed rate of 13.00%. A risk interval above 19.0% corresponds to an average daily temperature between −10.4 °C and 6.4 °C. In addition, a high-risk interval over 21.4% corresponds to an average daily temperature between −7.8 °C and 3.8 °C. The predicted norovirus infection rate among all diarrhea patients aged 6–15 years reaches a maximum of 21.2% at −0.5 °C; this is significantly greater than the overall observed rate of 8.74% for the same group. A risk interval above 17.0% corresponds to an average daily temperature of between −5.7 °C and 4.6 °C. A high-risk interval exceeding 19.1% corresponds to an average daily temperature between −4.1 °C and 3.0 °C. The predicted norovirus infection rate among all diarrhea patients aged 16–59 reaches a maximum of 11.9% at −5.8 °C, which is more than double the overall observed rate of 4.79% for the same group. A risk interval exceeding 9.5% corresponds to an average daily temperature between −14.3 °C and 2.7 °C. The high-risk interval exceeding 10.7% corresponds to an average daily temperature between −11.7 °C and 0.1°C. The predicted norovirus infection rate among all diarrhea patients aged 60 years or older reaches a maximum of 7.7% at −4.6 °C, which is significantly greater than the overall observed rate of 3.37% for the same group. The risk interval exceeding 6.2% corresponds to an average daily temperature between −12.9 °C and 3.6 °C. A high-risk interval, which exceeds 7.0%, corresponds to an average daily temperature between −10.3 °C and 1.1 °C. In particular, the incidence rates in patients aged 0–5 and 6–15 are much greater than the average rates. In this paper, we constructed scenarios with rates greater than the average (13% for 0–5 years old and 8.74% for 6–15 years old) where norovirus occurs as a priority for each age group.
For each age group, the RRI was calculated using the critical value of the incidence of norovirus among diarrhea patients. Table 4, Table 5, Table 6 and Table 7 and Figure 4 show the resulting values and 25 and 97.5 percentile of these RRIs by age relative to the critical values in the years 2030, 2050, 2070, and 2100 according to RCP scenarios 2.6, 4.5, 6.0, and 8.5. These results were obtained using percentile bootstrap methods over 100 trials. In Figure 4, the RRI shows a similar pattern in the same age groups and RCP scenarios. In the RCP 2.6 scenario, the RRI for the 0–5 years, 16–59 years, and over 60 years age groups will increase by 2050 and 2070, but in 2100 it is similar or only slightly higher than in 2030. In RCP 4.5, RRI tends to decrease, but it rises briefly in 2070. In RCP 6.0, the RRI increases until 2070 and then decreases slowly, and by 2100 the RRI is decreasing rapidly. In RCP 8.5, the RRI is maintained until 2050, then decreases significantly until 2100. In Figure 5, Figure 6, Figure 7 and Figure 8, we present the probability of being infected with norovirus according to the temperature calculated by GALM for each grid on the map of the Korean Peninsula. Overall, it can be seen that the more rapid the climate change, the lower the probability of infection. In addition, in terms of the regional analysis, the probability of infection decreases from the southern and eastern coasts. In Figure 6, it can be seen that in all RCPs, the probability of infection in the northeast area is low in 2030 and in 2050. It is lower in the coastal areas than in other areas for two age groups, 16 to 60 and over 60 years groups.

6. Discussion and Conclusions

In the previous section, we calculated the probability of norovirus occurrence for different age groups by applying the GALM model using average daily temperature as an explanatory variable, and we calculated the risk index according to RCP scenarios based on the RRI index presented in this paper. These results are consistent with the general observation that food poisoning caused by norovirus occurs most often in winter. Considering that the average winter temperature in Korea is 1.5 °C in December, −1.0 °C in January, and 1.1 °C in February, Korea’s winter season is favorable for norovirus to occur. The results show that the risk of norovirus occurrence increases in 2050 and 2070 but decreases by the year 2100 under the RCP 2.6 scenario for 0–5 years, 15–60 years, and over 60 years age groups. For all age groups, the risk of norovirus occurrence decreases over time under the RCP 4.5 scenario, despite a brief increase in 2070. Under the RCP 6.0 scenario, for all age groups, the risk increases until 2070 and then decreases to 2100. Under the RCP 8.5 scenario, for all age groups, the risk remains stable until 2050 and then decreases to 2100. In the prediction probability plots, we can see that the probability of occurrence decreases from near the coast, and that the probability of occurrence in the northeast area is low in 2030 and 2050 in all RCP scenarios. In general, the incidence of norovirus is predicted to decrease when the temperature rises due to climate change. However, the findings of this study indicate that, even if the temperature rises, the temperature in Korea will remain within the risk interval of temperature for norovirus occurrence, so the risk of norovirus will either remain stable or increase. Despite the predicted increase in temperature due to climate change, it is necessary to continue to prepare for future norovirus outbreaks in Korea, among those under the age of 15. Based on these simulation results, if the winter temperature rises, the possibility of food poisoning due to domestic norovirus will increase, and effective disease control measures will be needed. Of course, our ability to predict changes in the pattern of norovirus outbreaks simply by looking at the average daily temperature alone may be somewhat limited. A number of factors should be considered to improve the accuracy of these predictions, such as humidity and precipitation in particular, as these are closely related to microbial reproduction. In fact, there have been reports that lower water temperatures increase the viability of noroviruses in their natural environments [39,40], and, statistically, lower water temperatures in the Ontario River have been associated with many cases of food poisoning caused by noroviruses [41]. Norovirus can be transmitted not only through fecal and oral routes but also through aerosols caused by vomiting, so an increase in humidity is likely to increase pathogen viability [40]. In a similar vein, there have been reports that high humidity increases the viability of rotavirus [42]. In order to reasonably predict changes in the seasonal patterns of norovirus outbreaks, it is necessary to consider not only climate-related factors but also human behavior and socioeconomic changes that may occur with climate change. The cooler temperatures in winter cause people to remain inside longer, thus increasing the incidence of food poisoning from norovirus, meaning that changes in human behavior can affect the incidence of norovirus.
This paper calculated the risk index according to the temperature changes suggested in various RCP scenarios, but this approach has limitations, as the occurrence of norovirus is affected by other latent factors. For example, outbreaks of food poisoning caused by norovirus are often related to the consumption of fresh fruits and vegetables, so a significant increase in the consumption of fresh fruits and vegetables in winter could lead to a marked increase in norovirus food poisoning regardless of the climate conditions. In addition, 42.5% of cases of food poisoning caused by norovirus in Korea from 2000 to 2007 are known to have been caused by careless food handling by cooks, so improved hygiene practices among individuals working in large centers that handle fresh agricultural products could reduce incidences of norovirus. Nevertheless, based on the results of this study, we predict that the risk of norovirus in winter will either remain stable or increase in Korea despite climate change leading to warmer winters. This approach provides an overview of the likely changes that will occur in regional and seasonal norovirus patterns, which are expected to help appropriate policy decisions.

Author Contributions

Formal analysis—I.C.; Methodology Visualization—I.C., T.-K.K., J.P. and H.-Y.K.; Writing original draft—I.C., T.-K.K., J.P. and H.-Y.K.; Writing review and editing—I.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Project No.PJ015845)” Rural Development Administration, Republic of Korea.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors thank Yire Shin for his assistance about the RCP scenarios.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Troeger, C.; Colombara, D.V.; Rao, P.C.; Khalil, I.A.; Brown, A.; Brewer, T.G.; Guerrant, R.L.; Houpt, E.R.; Kotloff, K.L.; Misra, K.; et al. Global disability-adjusted life-year estimates of long-term health burden and undernutrition attributable to diarrhoeal diseases in children younger than 5 years. Lancet Glob. Health 2018, 6, e255–e269. [Google Scholar] [CrossRef] [Green Version]
  2. WHO (World Health Organization). Quantitative Risk Assessment of the Effects of Climate Change on Selected Causes of Death, 2030s and 2050s; Hales, S., Kovats, K., Lloyd, S., Campbell-Lendrum, D., Eds.; World Health Organization: Geneva, Switzerland, 2014. [Google Scholar]
  3. Kotloff, K.L.; Nataro, J.P.; Blackwelder, W.C.; Nasrin, D.; Farag, T.H.; Panchalingam, S.; Wu, Y.; Sow, S.O.; Sur, D.; Breiman, R.F.; et al. Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): A prospective, case-control study. Lancet 2013, 382, 209–222. [Google Scholar] [CrossRef]
  4. Walker, C.L.F.; Rudan, I.; Liu, L.; Nair, H.; Theodoratou, E.; Bhutta, Z.A.; O’Brien, K.; Campbell, H.; Black, R.E. Global burden of childhood pneumonia and diarrhoea. Lancet 2013, 381, 1405–1416. [Google Scholar] [CrossRef]
  5. Guerrant, R.L.; DeBoer, M.D.; Moore, S.; Scharf, R.; Lima, A. The impoverished gut—a triple burden of diarrhoea, stunting and chronic disease. Nat. Rev. Gastroenterol. Hepatol. 2013, 10, 220–229. [Google Scholar] [CrossRef] [PubMed]
  6. Prüss, A.; Kay, D.; Fewtrell, L.; Bartram, J. Estimating the burden of disease from water, sanitation, and hygiene at a global level. Environ. Health Perspect. 2002, 110, 537–542. [Google Scholar]
  7. Hall, A.J.; Glass, R.I.; Parashar, U.D. New insights into the global burden of noroviruses and opportunities for prevention. Expert Rev. Vaccines 2016, 15, 949–951. [Google Scholar] [CrossRef] [Green Version]
  8. Havelaar, A.H.; Kirk, M.D.; Torgerson, P.R.; Gibb, H.J.; Hald, T.; Lake, R.J.; Praet, N.; Bellinger, D.C.; de Silva, N.R.; Gargouri, N.; et al. World Health Organization global estimates and regional comparisons of the burden of foodborne disease in 2010. PLoS Med. 2015, 12, e1001923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Lopman, B.A. Global burden of norovirus and prospects for vaccine development. In CDC Foundation Global Burden Report; CDC: Atlanta, GA, USA, 2015; pp. 1–44. [Google Scholar]
  10. Mead, P.S.; Slutsker, L.; Griffin, P.M.; Tauxe, R.V. Food-related illness and death in the United States reply to Dr. Hedberg. Emerg. Infect. Dis. 1999, 5, 607–625. [Google Scholar] [CrossRef] [PubMed]
  11. Lopman, B.A.; Reacher, M.; Gallimore, C.; Adak, G.K.; Gray, J.J.; Brown, D.W. A summertime peak of “winter vomiting disease”: Surveillance of noroviruses in England and Wales, 1995 to 2002. BMC Public Health 2003, 3, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Verhoef, L.; Depoortere, E.; Boxman, I.; Duizer, E.; Van Duynhoven, Y.; Harris, J.; Johnsen, C.; Kroneman, A.; Le Guyader, S.; Lim, W.; et al. Emergence of new norovirus variants on spring cruise ships and prediction of winter epidemics. Emerg. Infect. Dis. 2008, 14, 238–243. [Google Scholar] [CrossRef] [PubMed]
  13. McSwiggan, D.A.; Cubitt, D.; Moore, W. Calicivirus associated with winter vomiting disease. Lancet 1978, 311, 1215. [Google Scholar] [CrossRef]
  14. Kim, J.G.; Kim, J.S.; Kim, J.G. Characteristics of Norovirus Food Poisoning Outbreaks in Korea in the 2000s. J. Food Prot. 2021, 84, 472–480. [Google Scholar] [CrossRef] [PubMed]
  15. Ahmed, S.M.; Lopman, B.A.; Levy, K. A systematic review and meta-analysis of the global seasonality of norovirus. PLoS ONE 2013, 8, e75922. [Google Scholar]
  16. Stegmaier, T.; Oellingrath, E.; Himmel, M.; Fraas, S. Differences in epidemic spread patterns of norovirus and influenza seasons of Germany: An application of optical flow analysis in epidemiology. Sci. Rep. 2020, 10, 1–14. [Google Scholar] [CrossRef]
  17. Onozuka, D.; Gasparrini, A.; Sera, F.; Hashizume, M.; Honda, Y. Modeling future projections of temperature-related excess morbidity due to infectious gastroenteritis under climate change conditions in Japan. Environ. Health Perspect. 2019, 127, 077006. [Google Scholar] [CrossRef] [PubMed]
  18. Moss, R.H.; Babiker, M.; Brinkman, S.; Calvo, E.; Carter, T.; Edmonds, J.A.; Elgizouli, I.; Emori, S.; Lin, E.; Hibbard, K.; et al. Towards New Scenarios for Analysis of Emissions, Climate Change, Impacts, and Response Strategies; Intergovernmental Panel on Climate Change Secretariat (IPCC): Geneva, Switzerland, 2008; ISBN 978-92-9169-125-8. [Google Scholar]
  19. Juang, H.M.H.; Hong, S.Y.; Kanamitsu, M. The NCEP Regional Spectral Model: An update. Bull. Am. Meteorol. Soc. 1997, 78, 2125–2143. [Google Scholar] [CrossRef] [Green Version]
  20. Kang, H.S.; Cha, D.H.; Lee, D.K. Evaluation of the mesoscale model/land surface model (MM5 = LSM) coupled model for East Asian summer monsoon simulations. J. Geophys. Res. 2005, 110, D10. [Google Scholar] [CrossRef] [Green Version]
  21. Skamarock, W.C.; Klemp, J.B.; Dudhia, J.; Gill, D.O.; Barker, D.M.; Wang, W.; Powers, J.G. A Description of the Advanced Research WRF Version 2; National Center For Atmospheric Research Boulder Co Mesoscale and Microscale Meteorology Div: Boulder, CO, USA, 2005; 100p. [Google Scholar]
  22. Giorgi, F.; Coppola, E.; Solmon, F.; Mariotti, L.; Sylla, M.B.; Bi, X.; Elguindi, N.; Diro, G.T.; Nair, V.; Giuliani, G.; et al. RegCM4: Model description and preliminary tests over multiple CORDEX domains. Clim. Res. 2012, 52, 7–29. [Google Scholar] [CrossRef] [Green Version]
  23. Kim, S.; Kim, H.; Jung, Y.; Heo, J.H. Assessment of Frequency Analysis using Daily Rainfall Data of HadGEM3-RA Climate Model. J. Wetl. Res. 2019, 21, 51–60. [Google Scholar]
  24. Jo, W.S.; Kim, H.Y.; Kim, B.J. Climate change alters diffusion of forest pest: A model study. J. Korean Phys. Soc. 2017, 70, 108–115. [Google Scholar] [CrossRef] [Green Version]
  25. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 31 August 2021).
  26. Wickham, H. The Split-Apply-Combine Strategy for Data Analysis. J. Stat. Softw. 2011, 40, 1–29. Available online: https://www.jstatsoft.org/v40/i01/ (accessed on 31 August 2021). [CrossRef] [Green Version]
  27. Canty, A.; Ripley, B. Boot: Bootstrap R (S-Plus) Functions. R Package Version 1.3-20. 2017. [Google Scholar]
  28. Garnier, S.; Ross, N.; Rudis, B.; Sciaini, M.; Camargo, A.P.; Scherer, C.; Garnier, S. Viridis: Default Color Maps from ‘Matplotlib’. R Package Version 0.5.1. Available online: https://CRAN.R-project.org/package=viridis (accessed on 31 August 2021).
  29. Dowle, M.; Srinivasan, A.; Gorecki, J.; Chirico, M.; Stetsenko, P.; Short, T.; Lianoglou, S.; Lianoglou, E.; Antonyan, E.; Bonsch, M.; et al. Data. Table: Extension of ‘Data.Frame’. R Package Version 1.11.8. 2018. Available online: https://CRAN.R-project.org/package=data.table (accessed on 31 August 2021).
  30. Bivand, R.; Keitt, T.; Rowlingson, B.; Pebesma, E.; Sumner, M.; Hijmans, R.; Baston, D.; Warmerdam, F.; Ooms, J. Rgdal: Bindings for the ‘Geospatial’ Data Abstraction Library. R Package Version 1.3-9. 2019. Available online: https://CRAN.R-project.org/package=rgdal (accessed on 31 August 2021).
  31. Hastie, T.; Tibshirani, R. Generalized additive models: Some applications. J. Am. Stat. Assoc. 1987, 82, 371–386. [Google Scholar] [CrossRef]
  32. Hastie, T.; Tibshirani, R. Generalized additive models for medical research. Stat. Methods Med Res. 1995, 4, 187–196. [Google Scholar] [CrossRef] [PubMed]
  33. Hastie, T.J.; Tibshirani, R.J. Generalized Additive Models, 1st ed.; CRC Press: Boca Raton, FL, USA, 1990; ISBN 0-412-34390-8. [Google Scholar]
  34. Wold, S. Spline functions in data analysis. Technometrics 1974, 16, 1–11. [Google Scholar] [CrossRef]
  35. Craven, P.; Wahba, G. Smoothing noisy data with spline functions. Numer. Math. 1978, 31, 377–403. [Google Scholar] [CrossRef]
  36. Hastie, T.J.; Tibshirani, R.J.; Friedman, J. The Elements of Statistical Learning: Data Mining, Inference, and Prediction; Springer: New York, NY, USA, 2009; ISBN 978-0-387-84857-0. [Google Scholar]
  37. Perperoglou, A.; Sauerbrei, W.; Abrahamowicz, M.; Schmid, M. A review of spline function procedures in R. BMC Med. Res. Methodol. 2019, 19, 1–16. [Google Scholar] [CrossRef] [Green Version]
  38. Wahba, G. Spline Models for Observational Data; Society for Industrial and Applied Mathematics, University City Science Center: Philadelphia, PA, USA, 1990; ISBN 978-0-898712-44-5. [Google Scholar]
  39. Rohayem, J. Norovirus seasonality and the potential impact of climate change. Clin. Microbiol. Infect. 2009, 15, 524–527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Xiao, S.; Tang, J.W.; Li, Y. Airborne or fomite transmission for norovirus? A case study revisited. Int. J. Environ. Res. Public Health 2017, 14, 1571. [Google Scholar] [CrossRef] [Green Version]
  41. Greer, A.L.; Drews, S.J.; Fisman, D.N. Why “winter” vomiting disease? Seasonality, hydrology, and norovirus epidemiology in Toronto, Canada. Ecohealth 2009, 6, 192–199. [Google Scholar] [CrossRef] [PubMed]
  42. Nathanson, N.; Martin, J.R. The epidemiology of poliomyelitis: Enigmas surrounding its appearance, epidemicity, and disappearance. Am. J. Epidemiol. 1979, 110, 672–692. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Three study areas in the South Korean Peninsula separated based on latitude and administrative divisions.
Figure 1. Three study areas in the South Korean Peninsula separated based on latitude and administrative divisions.
Life 11 01332 g001
Figure 2. The 96,172 grid points based on the land of the southern Korean peninsula to calculate the RRIs and the GALM model.
Figure 2. The 96,172 grid points based on the land of the southern Korean peninsula to calculate the RRIs and the GALM model.
Life 11 01332 g002
Figure 3. Predicted norovirus incidence rates for four age groups.
Figure 3. Predicted norovirus incidence rates for four age groups.
Life 11 01332 g003
Figure 4. Plots of RRI in 2030, 2050, 2070, and 2100 under four RCP scenarios (RCP2.6, RCP4.5, RCP6.0, and RCP8.5).
Figure 4. Plots of RRI in 2030, 2050, 2070, and 2100 under four RCP scenarios (RCP2.6, RCP4.5, RCP6.0, and RCP8.5).
Life 11 01332 g004aLife 11 01332 g004b
Figure 5. The prediction probability plots of norovirus incidence by the GALM model for Age 0–5 years.
Figure 5. The prediction probability plots of norovirus incidence by the GALM model for Age 0–5 years.
Life 11 01332 g005aLife 11 01332 g005b
Figure 6. The prediction probability plots of norovirus incidence by the GALM model for Age 6–15 years.
Figure 6. The prediction probability plots of norovirus incidence by the GALM model for Age 6–15 years.
Life 11 01332 g006
Figure 7. The prediction probability plots of norovirus incidence by the GALM model for Age 16–59 years.
Figure 7. The prediction probability plots of norovirus incidence by the GALM model for Age 16–59 years.
Life 11 01332 g007
Figure 8. The prediction probability plots of norovirus incidence by the GALM model for Age over 60 years.
Figure 8. The prediction probability plots of norovirus incidence by the GALM model for Age over 60 years.
Life 11 01332 g008aLife 11 01332 g008b
Table 1. Clinical characteristics of Korean patients with diarrhea (2005–2018).
Table 1. Clinical characteristics of Korean patients with diarrhea (2005–2018).
CharacteristicNumber of Patients with Norovirus(%)Total Diarrhea Patients
Total24,642 (8.69)283,651
GenderMale13,787 (8.82)156,397
Female10,855 (8.53)127,254
Age0–517,390 (13.00)133,737
6–152062 (8.74)23,605
16–593139 (4.79)65,524
over 602051 (3.37)60,785
RegionSouth Area11,706 (9.26)126,467
Central Area4452 (6.05)73,533
North Area8484 (10.14)83,651
Table 2. Results of logistic regression analysis.
Table 2. Results of logistic regression analysis.
Covariatesβexp (β)sep-Value95% CI
Intercept−1.8840.1520.012<0.0010.148–0.156
GenderMale Reference
Female0.011.010.0140.4610.983–1.037
Age0–5 Reference
6–15−0.4370.6460.024<0.0010.616–0.678
16–59−1.0510.3500.02<0.0010.336–0.364
over 60−1.4050.2450.024<0.0010.234–0.257
RegionSouth Area Reference
Central Area−0.2520.7770.019<0.0010.749–0.806
North Area0.061.0620.015<0.0011.031–1.094
Table 3. Calculated maximum norovirus incidence rate and temperature by age group.
Table 3. Calculated maximum norovirus incidence rate and temperature by age group.
AgeMaximum
Rate
Maximum Rate
Temperature
Risk IntervalHigh-Risk Interval
RateTemperature RangeRateTemperature Range
0–523.8%−2 °C19.0%−10.4 °C6.4 °C21.4%−7.8 °C3.8 °C
6–1521.2%−0.5 °C17.0%−5.7 °C4.6 °C19.1%−4.1 °C3 °C
16–5911.9%−5.8 °C9.5%−14.3 °C2.7 °C10.7%−11.7 °C0.1 °C
Over 607.7%−4.6 °C6.2%−12.9 °C3.6 °C7.0%−10.3 °C1.1 °C
Table 4. RRI for the 0–5 years age group in the four RCP scenarios.
Table 4. RRI for the 0–5 years age group in the four RCP scenarios.
AgeRCPYearRRIse2.5%97.50%
0–5RCP2.62030744,363.19156.249744,238.988744,475.598
2050798,557.16558.093798,448.003798,682.025
2070799,461.84358.131799,364.242799,581.800
2100769,219.31364.379769,104.059769,342.982
RCP4.52030797,518.05950.979797,417.997797,620.212
2050746,507.32254.491746,404.296746,603.098
2070765,325.67866.426765,193.106765,443.671
2100719,243.67072.033719,098.98719,370.906
RCP6.02030773,885.89859.472773,778.616774,007.359
2050819,199.54145.961819,120.154819,287.547
2070796,606.39248.595796,516.908796,685.359
2100680,299.53178.89680,118.010680,425.951
RCP8.52030793,144.66052.257793,056.644793,233.400
2050789,608.88458.676789,479.547789,718.438
2070696,860.63577.764696,731.282697,019.472
2100637,696.71476.576637,551.772637,827.478
Table 5. RRI for the 6–15 years age group in the four RCP scenarios.
Table 5. RRI for the 6–15 years age group in the four RCP scenarios.
AgeRCPYearRRIse2.5%97.50%
6–15RCP2.62030759,573.41587.567759,392.791759,752.307
2050830,048.75687.292829,897.173830,171.827
2070847,114.17185.711846,989.130847,286.762
2100824,541.47685.232824,408.661824,737.405
RCP4.52030797,031.13992.947796,868.583797,198.477
2050780,405.36588.667780,244.899780,579.088
2070818,259.56091.263818,110.020818,439.662
2100744,777.347109.334744,567.330744,988.180
RCP6.02030814,398.39476.988814,255.317814,547.288
2050873,433.72878.062873,283.645873,577.149
2070850,347.59585.542850,199.707850,538.008
2100721,090.103102.792720,891.629721,264.156
RCP8.52030839,317.62180.411839,185.498839,501.208
2050842,927.70082.334842,793.436843,091.020
2070726,317.494119.128726,139.886726,580.614
2100675,713.129112.328675,513.487675,906.920
Table 6. RRI for the 16–59 years age group in the four RCP scenarios.
Table 6. RRI for the 16–59 years age group in the four RCP scenarios.
AgeRCPYearRRIse2.5%97.50%
16–59RCP2.62030443,918.39956.697443,813.478444,020.317
2050487,006.04742.343486,927.932487,083.852
2070478,341.41636.410478,277.169478,404.044
2100441,578.84046.270441,503.240441,667.933
RCP4.52030514,711.08838.449514,641.502514,797.870
2050434,542.39249.690434,429.309434,637.481
2070444,944.11850.705444,840.469445,052.900
2100419,674.71055.697419,577.711419,794.594
RCP6.02030470,768.44946.386470,685.615470,851.935
2050493,738.98138.454493,661.984493,798.576
2070469,803.79844.892469,712.005469,881.028
2100374,173.54560.717374,068.956374,298.359
RCP8.52030472,046.92447.130471,954.943472,142.872
2050464,492.95247.240464,407.449464,577.836
2070397,364.82963.693397,242.417397,481.481
2100343,059.16556.321342,955.761343,160.810
Table 7. RRI for those aged more than 60 years in the four RCP scenarios.
Table 7. RRI for those aged more than 60 years in the four RCP scenarios.
AgeRCPYearRRIse2.5%97.50%
over 60RCP2.62030282,203.50426.677282,156.802282,257.290
2050308,661.95526.281308,604.289308,703.138
2070304,635.01725.072304,589.212304,683.773
2100283,677.65728.946283,622.792283,730.317
RCP4.52030321,812.21320.200321,776.156321,848.276
2050277,891.56830.253277,824.773277,943.525
2070284,957.85633.514284,888.233285,012.275
2100267,956.16331.890267,900.299268,015.764
RCP6.02030298,005.49126.261297,947.565298,047.385
2050314,090.96021.365314,052.541314,135.816
2070300,210.74029.942300,160.387300,264.328
2100241,975.56836.335241,903.235242,032.298
RCP8.52030300,852.81524.281300,802.806300,888.892
2050296,927.58126.964296,875.699296,971.395
2070254,937.22935.084254,858.145254,992.024
2100222,659.98135.031222,595.774222,712.245
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, T.-K.; Paek, J.; Kim, H.-Y.; Choi, I. Relative Risk Prediction of Norovirus Incidence under Climate Change in Korea. Life 2021, 11, 1332. https://doi.org/10.3390/life11121332

AMA Style

Kim T-K, Paek J, Kim H-Y, Choi I. Relative Risk Prediction of Norovirus Incidence under Climate Change in Korea. Life. 2021; 11(12):1332. https://doi.org/10.3390/life11121332

Chicago/Turabian Style

Kim, Tae-Kyoung, Jayeong Paek, Hwang-Yong Kim, and Ilsu Choi. 2021. "Relative Risk Prediction of Norovirus Incidence under Climate Change in Korea" Life 11, no. 12: 1332. https://doi.org/10.3390/life11121332

APA Style

Kim, T. -K., Paek, J., Kim, H. -Y., & Choi, I. (2021). Relative Risk Prediction of Norovirus Incidence under Climate Change in Korea. Life, 11(12), 1332. https://doi.org/10.3390/life11121332

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop