This article has Open Peer Review reports available.
Quantifying the role of weather on seasonal influenza
- Marion Roussel^{1, 2}Email author,
- Dominique Pontier^{1, 2},
- Jean-Marie Cohen^{3},
- Bruno Lina^{4, 5} and
- David Fouchet^{1, 2}
https://doi.org/10.1186/s12889-016-3114-x
© Roussel et al. 2016
Received: 19 December 2015
Accepted: 12 May 2016
Published: 26 May 2016
Abstract
Background
Improving knowledge about influenza transmission is crucial to upgrade surveillance network and to develop accurate predicting models to enhance public health intervention strategies. Epidemics usually occur in winter in temperate countries and during the rainy season for tropical countries, suggesting a climate impact on influenza spread. Despite a lot of studies, the role of weather on influenza spread is not yet fully understood. In the present study, we investigated this issue at two different levels.
Methods
First, we evaluated how weekly (intra-annual) incidence variations of clinical diseases could be linked to those of climatic factors. We considered that only a fraction of the human population is susceptible at the beginning of a year due to immunity acquired from previous years. Second, we focused on epidemic sizes (cumulated number of clinical reported cases) and looked at how their inter-annual and regional variations could be related to differences in the winter climatic conditions of the epidemic years over the regions. We quantified the impact of fifteen climatic variables in France using the Réseau des GROG surveillance network incidence data over eleven regions and nine years.
Results
At the epidemic scale, no impact of climatic factors was highlighted. At the intra-annual scale, six climatic variables had a significant impact: average temperature (5.54 ± 1.09 %), absolute humidity (5.94 ± 1.08 %), daily variation of absolute humidity (3.02 ± 1.17 %), sunshine duration (3.46 ± 1.06 %), relative humidity (4.92 ± 1.20 %) and daily variation of relative humidity (4.46 ± 1.24 %). Since in practice the impact of two highly correlated variables is very hard to disentangle, we performed a principal component analysis that revealed two groups of three highly correlated climatic variables: one including the first three highlighted climatic variables on the one hand, the other including the last three ones on the other hand.
Conclusions
These results suggest that, among the six factors that appeared to be significant, only two (one per group) could in fact have a real effect on influenza spread, although it is not possible to determine which one based on a purely statistical argument. Our results support the idea of an important role of climate on the spread of influenza.
Keywords
Background
Influenza is one of the most significant diseases in humans, generating worldwide annual epidemics, which result in about three to five million cases of severe illness, and about 250,000 to 500,000 deaths [1]. Improving influenza knowledge about key epidemiological parameters such as survival, transmission and reproduction in hosts is essential to upgrade surveillance network and to develop more accurate predicting models. Better epidemic predictions would set up more appropriate public health prevention and intervention strategies.
Epidemics occur mainly during the winter season months in temperate countries [2–4] unlike in tropical and sub-tropical countries where they generally happen during the rainy season [5–8]. These differences suggest a climate impact on influenza spread. Climate might affect influenza diffusion (onset, duration, size) by impacting individuals’ contact rates (frequency and duration), population immunity and virus survival outside human body. The role of weather is however not fully understood [9] despite a lot of laboratory studies of host susceptibility according to environmental conditions [10–12] and mathematical modeling approaches analyzing the link between influenza morbidity or mortality and climatic factors [13–18].
Various climatic factors such as temperature, humidity, rainfalls, UV radiation, sunshine duration and wind speed might have an impact on influenza spread. In temperate countries, humidity and temperature might play an important role in influenza spread. Several laboratory works showed that a cold and dry weather promotes a higher virus survival outside human body and a better transmission [11, 19]. Cold air inhalation chills nasal epithelium leading to an inhibition of mechanical defenses of the respiratory mucosa and of the immune system [20]. Otherwise models explaining influenza epidemics (e.g., onset, peak, mortality) according to climatic factors reinforce the role of humidity and temperature in influenza spread in the United States [13, 15] as well as in Europe [16, 21]. Rainfalls might have an impact in tropical and sub-tropical countries such as in Central and South America [22–24] and in Asia [25, 26]. Another theory suggests a link between vitamin D secretion and influenza immunity, which is supported by experiments [27, 28]. As UV radiation is involved in vitamin D production, a lack of UV radiation in winter, for temperate countries, leads to a reduction of vitamin D production and might boost influenza epidemics [29, 30]. Dowell [31] also suggested a role of dark/light cycles and photoperiod on the immune systems caused by melatonin fluctuations. Thereby UV radiation and sunshine duration might have an indirect effect on influenza infections. Finally in China, Xiao et al. [32] proposed that a low wind speed contribute to influenza spread. In fact a strong wind speed may have a dispersive effect on influenza in the environment limiting its diffusion.
The aim of this study is to quantify the impact of several climatic factors such as temperature, humidity, and rainfalls, on influenza epidemics in France, a temperate country. The role of weather can be estimated based on the variation of influenza propagation in an area according to its climate variation. Usually studies compared observed to modeled epidemics taking into account climatic factors by comparing incidence or mortality within an epidemic year [13–18]. The impact of the climatic factors included in the model is supported if modeled and observed epidemics are similar. However little information is available about influenza transmission. Modeling approaches made a lot of hypotheses about the within host virus dynamic such as incubation and infectious periods R _{0} etc. Such hypotheses may have a strong impact on influenza propagation, which might lead to a misestimating of climatic effects. In order to reduce the set of model hypotheses, we built an autoregressive model based on the shape of the observed epidemics over time. We explained the intraseasonal variation of incidence of eleven French regions and for nine epidemic years (an epidemic year corresponds to October of a year until April of the year after) with the climatic variables listed before, to quantify their respective impact globally over all regions, then specifically in each region for significant climatic variables. The originality of our model is to consider that only a fraction of the human population is susceptible at the beginning of a year due to immunity acquired from previous years. Considering loss of immunity in modeling influenza epidemics might be important [33] even if almost no studies about influenza and climate take it into account to our knowledge. Here we called susceptible individuals people that could be infected and develop symptoms, as we only had data about infected people presenting symptoms. We then quantified potential effects of climatic factors on the interseasonal variation of influenza epidemics. To do that we built an autoregressive linear model that explains the epidemic size according to the average value of the climatic factors over an epidemic year for the nine epidemic years and the eleven French regions.
Methods
Data
Epidemiological data
Epidemiological data come from the Réseau des GROG (Regional Influenza Surveillance Group) sentinel network, which is a French surveillance network made up of general practitioners and pediatricians. These physician sentinels identify cases of respiratory pathogens including influenza. Each region has on average 25 sentinels (from 10 to 75 depending on regions and epidemic years) involved in the Réseau des GROG sentinel network. Every week from October to April, they describe in reports the intensity of their activity by giving the number of days they worked, the number of medical acts performed and the number of acute respiratory infection (ARI) defined as the sudden onset of at least one respiratory sign (cough, rhinitis, coryza, etc.) and at least one systemic sign suggesting an acute infectious context (fever, fatigue, headache, myalgia, malaise, etc.). In addition, sentinels randomly realize nasal/pharyngeal swab samples on patients with a less than 48 h ARI. Analysis of these samples allows virological confirmation of influenza infections. Using the weekly information reported by each physician sentinel (clinical reports and virological samples analysis), the Réseau des GROG sentinel network is able to provide an estimate of the number of influenza-infected individuals called the influenza incidence.
where ARI _{ GP }(t) and ARI _{ Ped }(t) stand for the number of ARI cases for week t, respectively, reported by general practitioners (GP) and pediatricians from the Réseau des GROG sentinel network. GP _{ region } and Ped _{ region } are, respectively, the number of GP and pediatricians of a region. GP _{ GROG participants }(t) and Ped _{ GROG participants }(t) represent the number of GP and pediatricians who participated in surveillance the week t, respectively. Age of infected individuals was not taken into account assuming that climatic factors have a uniform impact on influenza spread within the population.
Epidemiological data are available from the 2003–2004 epidemic year to the 2012–2013 epidemic year. However we excluded the 2009–2010 epidemic year where the H1N1 pandemic happened in order to study only seasonal epidemics.
Climatic data
Climatic data were provided by Météo-France, the French national meteorological service. We picked 65 meteorological stations (see Fig. 1) to collect data in order to estimate climatic factors that globally describe each region. We had information on temperature, relative humidity, absolute humidity, rainfalls, sunshine duration (very correlated to UV radiation), and wind speed (see Additional file 1). It is not necessarily easy to choose efficient climatic factors, as illustrated by Davis et al. [34] who highlighted the challenge of selecting an appropriate measure of the humidity covariate.
Definition of the climatic variables
Climatic variable | Formula |
---|---|
Average temperature | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations} average\ daily\ temperature\left(j,s\right)}} \) |
Daily variation of temperature | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations}\left( \max daily\ temperature\left(j,s\right)- \min daily temperature\left(j,s\right)\right)}} \) |
Relative weekly variation of temperature | Average temperature(week _{ t + 1}) − Average temperature(week _{ t }) |
Absolute weekly variation of temperature | |Average temperature(week _{ t + 1}) − Average temperature(week _{ t })| |
Average relative humidity | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations} average\ daily\ relative\ humidity\left(j,s\right)}} \) |
Daily variation of relative humidity | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations}\left( \max daily\ relative\ humidity\left(j,s\right)- \min daily\ relative\ humidity\left(j,s\right)\right)}} \) |
Relative weekly variation of relative humidity | Average relative humidity(week _{ t + 1}) − Average relative humidity(week _{ t }) |
Absolute weekly variation of relative humidity | |Average relative humidity(week _{ t + 1}) − Average relative humidity(week _{ t })| |
Average absolute humidity | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations} average\ daily\ absolute\ humidity\left(j,s\right)}} \) |
Daily variation of absolute humidity | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations}\left( \max daily\ absolute\ humidity\left(j,s\right)- \min daily\ absolute\ humidity\left(j,s\right)\right)}} \) |
Relative weekly variation of absolute humidity | Average absolute humidity(week _{ t + 1}) − Average absolute humidity(week _{ t }) |
Absolute weekly variation of absolute humidity | |Average absolute humidity(week _{ t + 1}) − Average absolute humidity(week _{ t })| |
Average wing speed | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days}{\displaystyle {\sum}_{s=0}^{stations} average\ daily\ wind\ speed\left(j,s\right)}} \) |
Rainfall height | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days} average\ daily\ rainfalls\ height\left(j,s\right)} \) |
Sunshine duration | \( \frac{1}{days}\cdot \frac{1}{stations}\cdot {\displaystyle {\sum}_{j=0}^{days} average\ daily\ sunshine\ duration\left(j,s\right)} \) |
Mathematical models
Climatic factors can impact influenza spread by both increasing the transmissibility of the virus and/or by increasing the susceptibility of its human host. One particularity of our data set is that the variability in influenza incidence is reported at different scales: the transmission scale (intraseasonal variation) and the epidemic scale (interseasonal variation). The impact of climatic factors may occur at the two scales in which it will be observed in a slightly different way.
At the transmission scale – during a seasonal epidemic of a given year in a given region – favorable climatic (for influenza diffusion) factors will lead to observe an increase in disease (apparent) transmission. At this scale we will search for significant associations between weekly variations of climatic factors and those of the disease apparent transmission rate (defined below). Different observed epidemics (in all regions and epidemic years) will be treated as independent replicates.
At the epidemic scale - the impact of a climatic factor (in a region over an entire epidemic year) may mainly be observed by the increase or decrease in the epidemic size (the total number of infected individuals). At this scale we will search for significant associations between the size of the epidemic and the average value of the different climatic factors (over an epidemic year in a region). Because both scales imply different response variables, they will be treated separately and independently.
Impact of climatic factors at the transmission scale
We built an auto-regressive statistical model with a lag of one week to explain variations in the weekly influenza incidence with climatic factors for eleven French regions over nine epidemic years.
where c quantifies the link between F _{ c } and β, d is a constant and ε _{2} is a random term independent of F _{ c } modeling the fluctuation in β independent of F _{ c }, i.e., due to other factors.
where Y(t) is the logarithm of the estimated number of infected individuals. ε = ε _{1} + ε _{2} is the total residual error and it is assumed to be distributed according to a Gaussian centered distribution with a standard deviation σ.
We defined \( \widehat{\alpha}=\frac{I_{Tmax}}{\widehat{N}} \), which provides an estimate of the proportion of individuals who developed the disease (with symptoms) in the pool of individuals that could have developed it. I _{ Tmax } denotes here the time at which the influenza surveillance ends (mid-April). α = 1 means that all individuals who could potentially become sick acquired the infection, and suggests that the disease has a sufficient transmission to reach the entire susceptible pool of the population. At the opposite α < 1 suggests that the virus spread was not sufficient to reach the entire susceptible pool.
where X, Z, U, V and W are formulas depending on R and Y. To take a few examples, be x a generic variable that can be a, b, c… x(0) means that x = 0 in the model; x(1) means that x is constant (intercept model); x(R) means that x depends on the region; x(R + Y) means that x depends on both the region and epidemic year in an additive way and x(R ⋅ Y) in a multiplicative way.
The most complicated model considered (i.e., the complete model) is not the model where all parameters depend multiplicatively on R and Y (R ⋅ Y), which would contain too many parameters to be tractable. Since a and b are shape parameters for the spread of the epidemic, it is reasonable to assume that they are characteristics of the region (a(R) and b(R)). d affects the average transmission rate of the virus. It can be different between regions (which show different demographic characteristics) and between epidemic years (because the circulating influenza strain is different from one epidemic year to the next), but it is reasonable to consider that it will only be slightly affected by the interaction between these two factors (d(R ⋅ Y)). That is why the most complicated model considered was a(R), b(R), c(R), d(R + Y), α(R ⋅ Y).
where sd _{ φ } stands for the standard deviation of the variable φ. over all epidemic weeks t, epidemic year Y and region R.
In total, fifteen climatic factors were tested, leading to potentially important problems of multiple testing. Since climatic factors are not independent, applying a simple Bonferroni correction would lead to a loss of statistical power [37]. Instead, we preferred a multiple testing correction based on permutation tests [38]. The idea of the permutation test we developed here is to keep the same values for all the climatic factors but to shuffle the week indexes, within a given region and a given epidemic year (in order to break the potential association between any climatic factor and the observed course of the epidemic). Mathematically, let us call F _{ t,Y,R } the value of the climatic factor F during the t ^{ th } week of region R and epidemic year Y. Let us call P a permutation of the week indexes t. The permuted climatic factors (F) associated to permutation P in region R and year Y will be defined by: F _{ P(t),Y,R }. The main advantage of this permutation procedure is that it conserves the within epidemic year and region correlation structure between the climatic factors. One permutation of the climatic factors is then defined as a set of permutations (one for each epidemic year in each region) leading to a set of permuted climatic factors in all regions and for all epidemic years. Note that these permuted factors have strictly no reason to be correlated with the apparent disease transmission rate (the permutation is purely random) and hence can be considered as realizations of the null hypothesis H0 “the apparent transmission rate is not linked to any climatic factor”.
We used the absolute value of the maximum estimated climatic factor coefficients (c _{ max }) as a test statistic for H0. We generated 10,000 permutations of climatic factors (see above) and for each one we calculated c _{ max }, leading to 10,000 independent realizations of c _{ max } under H0. The 95 % quantile of the distribution defines a significant threshold. Climatic factors are considered being significantly linked to the apparent transmission rate if the absolute value of their c estimate from data is above the defined threshold. Model parameters are estimated using maximum likelihood. Standard errors of the estimations of the model parameters were determined using the square roots of the diagonal elements of the covariance matrix (the inverse of the negative of the expected value of the Hessian matrix). Model implementation and permutation tests were performed in Python.
Impact of climatic factors at the epidemic scale
To evaluate the impact of climatic factors at the epidemic scale we considered the ratio of cumulated number of infected individuals across the entire epidemic period (from the first week of epidemic of the first region in epidemic to the last week of epidemic of the last region in epidemic) to the total population – an indicator of the epidemic size – as a response variable (ES).
where a _{0}, b and c are constant model parameters and a _{ Y } (respectively a _{ R }) models potential systematic variations in the epidemic size between epidemic years (respectively regions). These two terms account for the fact that some regions may be more prone to important epidemics (e.g., due to population demography) and the strains circulating some epidemic years can be more virulent or affect a larger set of the human population due to more important genetic differences with the strains of the previous epidemic years. \( \overline{F_{Y,R}} \) denotes the mean value of climatic factor F over the entire epidemic year.
Foremost we selected model parameters (a _{ Y }, a _{ R } and b) using an AIC criterion and then we assessed the impact of climatic factors.
Multiple hypothesis testing was corrected as in the previous section. Values of Y and R were shuffled together (pairs of values for Y and R were randomly re-attributed to all epidemics). For a permutation P, new climatic factors were built as \( \overline{F_{P\left(Y,R\right)}} \). The advantage of this permutation procedure is that, as above, it keeps the covariance structure between the climatic factors. As previously the permutation test is used to determine a significant threshold for the c coefficients using the maximum absolute estimated value of the c coefficients as a statistic.
Model parameters were estimated using the classical tools of linear models implemented in R3.1.2 [39].
Results
Impact of climatic factors at the transmission scale
AIC selection at the transmission scale
Model | Number of parameters | AIC criterion |
---|---|---|
a(1), b(1), c(0), d(1), α(1) | 5 | 1455.00 |
a(1), b(R), c(0), d(1), α(1) | 15 | 1460.48 |
a(1), b(1), c(0), d(R), α(1) | 15 | 1462.00 |
a(1), b(1), c(0), d(1), α(Y) | 13 | 1468.29 |
a(1), b(1), c(0), d(1), α(R) | 15 | 1471.00 |
a(1), b(1), c(0), d(Y), α(1) | 13 | 1801.52 |
a(R), b(1), c(0), d(1), α(1) | 15 | 1824.58 |
Global climatic impacts
Climatic factor | c | Standard deviation | Impact (%) |
---|---|---|---|
Average absolute humidity | −0.0612 | 0.0108 | 5.94 |
Average temperature | −0.0570 | 0.0109 | 5.54 |
Average relative humidity | −0.0505 | 0.0120 | 4.92 |
Daily variation of relative humidity | 0.0436 | 0.0124 | 4.46 |
Sunshine duration | 0.0340 | 0.0106 | 3.46 |
Daily variation of absolute humidity | −0.0307 | 0.0117 | 3.02 |
Besides the evaluation of impact of climatic factors at the transmission scale, the model built allowed the estimate of the susceptible population for an epidemic year \( \widehat{N} \) with the definition of \( \widehat{\alpha} \) that provides an estimate of the proportion of individuals who developed the disease in the pool of individuals that could have developed it. In the fifteen climatic models, estimates of α were included between 0.98 and 1 with a very low standard deviation (< 0.01).
Impact of climatic factors at the epidemic scale
AIC selection at the epidemic scale
Model | Number of parameters | AIC criterion |
---|---|---|
ES _{ Y,R } = a _{0} + a _{ Y } + a _{ R } + b ⋅ ES _{ Y − 1,R } | 19 | −241.64 |
ES _{ Y,R } = a _{0} + a _{ Y } + a _{ R } | 18 | −243.42 |
ES _{ Y,R } = a _{0} + a _{ Y } | 8 | −225.85 |
ES _{ Y,R } = a _{0} + b ⋅ ES _{ Y − 1,R } | 2 | −157.11 |
ES _{ Y,R } = a _{0} | 1 | −144.70 |
ES _{ Y,R } = a _{0} + a _{ Y } + b ⋅ ES _{ Y − 1,R } | 9 | −234.99 |
ES _{ Y,R } = a _{0} + a _{ R } | 11 | −135.74 |
Considering that variations in epidemic size could not be explained by our (measured) climatic variables, we then tried to decompose these variations into three sources. First variations in region characteristics (e.g., population size or non-measured climatic factors) can lead to systematic differences between regions. Second, temporal variations (e.g., in strain characteristics) can lead to systematic increase or decreased of epidemic sizes in all regions. Third, local conditions (in given epidemic years and regions) may also affect epidemic sizes. To quantify these three sources of variations, we built a model considering epidemic year and region as random variables: log(ES _{ Y,R }) = a _{0} + a _{ Y } + a _{ R } + ε, where a _{ Y } (respectively a _{ R }) is distributed according to a Gaussian centered distribution with a standard deviation σ _{ Y }. (respectively σ _{ R }). ε stands for the residual variations, taking into account the local variations of a given epidemic year and region, it is distributed according to a Gaussian centered distribution with a standard deviation σ _{ ε }. The homoscedasticity of the residuals is shown in Additional file 2: Figure S1.
Parameters were estimated with R.3.1.2 [39] using the package lme4 [43, 44]. We found \( {\widehat{\sigma}}_Y=0.036 \), \( {\widehat{\sigma}}_R=0.013 \) d \( {\widehat{\sigma}}_{\varepsilon }=0.0217 \) meaning that variations from one epidemic year to another one, from one region to another one and due to local conditions account for 50.9, 18.4 and 30.7 %, respectively.
Discussion
In the present paper, we presented the results of the analysis of the statistical link between influenza spread and fifteen climatic factors. Data were obtained from the French Réseau des GROG sentinel network. The network is based on voluntary practitioners who i) record acute respiratory infection and ii) randomly send nasal samples for an antigenic confirmation (or rejection) of influenza infection. Based on those two pieces of information, the Réseau des GROG sentinel network provides influenza incidence estimates of clinical cases. Two metrics were used for linking virus spread to climatic data: weekly incidence data of clinical cases and the epidemic size – measured as the total number of recorded clinical cases over the epidemic period.
Results of the analysis failed to isolate any correlation between epidemic size and climatic factors. Regarding weekly incidence data, we considered that incidence at time t was first affected by both the number of infected and susceptible individuals at time t − 1, as it is classically assumed in epidemic dynamic models of infectious diseases [36]. Six climatic factors were found to be significantly linked to influenza spread: average temperature, average absolute and relative humidity, daily variations of absolute and relative humidity as well as sunshine duration. However, a principal component analysis revealed that upon these six factors, two groups of three highly correlated factors could be separated. On a practical point of view, this implies that within each of the two groups, it is likely that only one factor has a biological link to influenza spread, the two remaining factors being linked to the disease spread because they are linked to the first factor (confounding effect).
The first group of factors is made up of average temperature and absolute humidity, and daily variations of absolute humidity. The role of a cold and dry weather on influenza spread has been highlighted from laboratory studies [19, 20] and modeling approaches in temperate countries [13, 15, 16, 21] including France [45]. Moreover models that included weekly variations of both temperature and absolute humidity in Israel [46] and in New York City [47] predicted reliable influenza epidemic estimations (better estimations with both factors than only one). That is why both the average temperature and absolute humidity seem to play an important role on the influenza spread.
The second group of factors is made up of average relative humidity, daily variations of relative humidity and sunshine duration. Both laboratory [11] and simulation [14] studies enhanced the impact of the relative humidity. About sunshine duration, a decrease of sunshine might favor influenza spread [31] but surprisingly our results showed a positive impact of sunshine duration on influenza epidemic spread. That is why the average relative humidity might impact influenza spread whereas sunshine duration might be a confounding factor.
Overall, the impact of the significant factors remained relatively low (a few percent). This is not surprising when we compare our finding with what is found in the literature (3 % impact of absolute humidity in the Netherlands - [21], less than 2 % impact of both absolute humidity and temperature on influenza mortality in the USA - [15]). However, it is important to raise reasonable hypotheses for explaining why the impact of climatic factors is found so low. First, low impacts can arise from the presence of important noise in data. The Réseau des GROG sentinel network is based on a limited number of voluntary practitioners, leading to noise in incidence estimates. Second, in order to obtain relatively reliable incidence estimates, we had to average incidence over entire regions. Climate and disease spread can be disparate within a region, leading to weaken the link between climatic factors and disease spread. Third, the model, which has a lag of one week (linking incidence at time t with the number of susceptible and infected individuals at time t − 1), can be a bit too simple. Actually simple compartmental models may not be sufficient to describe properly an influenza epidemic. Models are becoming more complex by, for example, taking into account more heterogeneous influenza transmission in the population (e.g., agent-based model) and including a contact network among people [48–50]. Finally, correlation between influenza spread and single climatic factors can be too simplistic. Climate can have a strong impact on disease spread, but on a more complex way involving several factors and potential interactions between these factors. Such combinations of factors were not considered in the model because it would have led to a huge number of hypotheses’ testing. Such an investigation of the most relevant combinations of climatic factors would be more relevantly achieved using descriptive statistics, but this was not the purpose of our study.
Another important question arising from our results is about of the disparity of the link of climatic factors with influenza spread using weekly incidence data and epidemic size data. The first obvious potential explanation is the lower statistical power associated to epidemic size data. Epidemic size is estimated only once per year while incidence is estimated every week. So epidemic size data contain less statistical information. An interesting alternative hypothesis could be that epidemic size and weekly incidence data capture different biological phenomena. Basically, incidence (corrected by the number of susceptible and infected individuals) may vary between weeks according to climatic factors for two reasons: i) because individuals are more susceptible to develop the clinical form of the infection and ii) because infection is more likely, i.e., the virus transmission rate increases. Epidemic size is schematically the result (product) of two phenomena: i) the proportion of individuals in the region that are susceptible to develop the clinical form of the disease upon infection and ii) the fraction of these individuals that will be reached by the virus, i.e., that will effectively become infected. If the latter phenomenon is linked to the virus transmission rate, the link is not linear. In particular, for large enough transmission rates, all susceptible individuals become infected during an epidemic and this term is poorly affected by the transmission rate. Interestingly, in that case, epidemic sizes are mainly an indicator of individuals’ susceptibility and hence contain information that differs from that of weekly incidence data.
The proportion of the susceptible (to the clinical disease) population that ultimately develops the disease is an important quantity for both data analysis interpretation and disease management. In data analysis, it will tell us how to interpret epidemic size data. When all susceptible individuals acquire the infection, then epidemic size is an indicator of the proportion of susceptible individuals in the population, i.e., the proportion of individuals that are in a healthy state (in terms of innate and acquire immunity) that does not permit them to control the disease upon infection. On a management point of view, if all individuals acquire the infection, this means that the virus transmission rate is high and reducing it will not necessarily lead to reduce its impact.
In our study, we introduced a term that we interpreted as the proportion of susceptible individuals who ultimately got infected. This is an interesting result, but which should be interpreted with great caution. First susceptibility is here defined as the ultimate development of the disease upon infection. It is hence not necessarily equivalent to susceptibility defined by antibody profiles. Second, it is important to recall that it is primarily a model parameter introduced for statistical convenience (i.e., a shape parameter). The fact that it equals one in our model only means that the decay in disease incidence at the end of the epidemic can be explained without having to assume any susceptible pool that would have escaped the infection. Since the study was not designed for estimating this biological quantity, we invite the reader not to interpret it as a formal estimation procedure of the proportion of susceptible individuals, but as a point raising interesting questions.
Several improvements could be brought to our analysis. First, it would be interesting to differentiate between the different subtypes of influenza. Influenza epidemics are often due to several subtypes that generate potentially shifted epidemics [51]. Practically, in our model this would imply that the number of susceptible individuals does not necessarily decreases with the cumulated number of influenza cases from all subtypes, but is subtype specific. Even though the use of permutation tests tends to reduce this problem, it would still be interesting to study the different subtypes separately because they might be differentially affected by climatic factors. Unfortunately, this information was not available in our data set.
The second interesting improvement that could be brought to our model is the consideration of different age-classes. Indeed, influenza is known to spread differentially within and between age-classes [52–54]. However, introducing age-classes in our model would tend to make it more complex. In the current paper we adopted a practical point of view by considering only the global spread of the epidemic without considering the heterogeneity of individuals that may exist within a population (age-classes, social classes, job-dependent degree of exposure, etc.).
Conclusion
Proper modeling of the relationship between climatic variables and infectious diseases spread and impact presents a challenging task. We presented a way to conciliate statistical and dynamical models of infectious diseases in a way that keeps the simplicity of statistical approach while introducing key knowledge about infectious dynamics (such as the decay of incidence after the epidemic peak).
We performed our study on two important influenza response variables at two levels: intra- and inter-annually. Linking variations of weekly incidence data with climatic factors is relevant because it allows anticipating the decay or increase in the number of cases of influenza in the weeks to come. The epidemic size is also a very important measure because it allows quantifying the impact of influenza according to climatic factors. This is especially valuable in the context of global climate changes to anticipate the future impact of influenza.
Declarations
Acknowledgements
We acknowledge the practitioners of Réseau des GROG sentinel network and the labs involved in the surveillance. We thank Isabelle Daviaud from Open Rome who sorted the epidemiological data from the Réseau des GROG sentinel network and Annick Auffray from Météo-France for her kindly help with meteorological data. This work was archived using the computing facilities of the CC LBBE/PRABI and of the CC IN2P3. It was performed within the framework of the LABEX ECOFECT (ANR‐ 11‐ LABX‐ 0048) of Université de Lyon, within the program “Investissements d’Avenir” (ANR‐ 11-IDEX‐ 0007) operated by the French National Research Agency (ANR).
Availability of data
Epidemiological data are available in the Additional file 3 and climatic data are available in the Météo-France website.
Authors’ contributions
All authors participated to the design of the study. MR conducted the analysis and prepared the initial draft of the manuscript. DF and DP supervised the analysis and writing of the manuscript. JMC is a general practioner specialist in influenza surveillance and BL is a virologist specialist in influenza virus. All authors contributed to the writing of and critically revised the manuscript. All authors approved the final version of the manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethical approval and consent to participate
Surveillance forms were routinely used in the influenza seasons, and oral informed consent was obtained from the ARI patient at the moment of swab taking in accordance with national regulations. All swab results and forms were anonymized by the laboratories before they were sent to the GROG network coordination, and only identified by the number given by each laboratory for virological tests. In accordance with the French applicable law n°2011–2012 of the 29th December, article 5, no clearance of an Ethics Committee is required in France for the retrospective analysis of anonymized data collected within routine influenza surveillance schemes.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
Authors’ Affiliations
References
- World Health Organization: Influenza (Seasonal). 2014 http://www.who.int/mediacentre/factsheets/fs211/en/.
- Viboud C, Boëlle P-Y, Pakdaman K, Carrat F, Valleron A-J, Flahault A. Influenza Epidemics in the United States, France, and Australia, 1972–1997. Emerg Infect Dis J. 2004;10:32–9.View ArticleGoogle Scholar
- Tamerius JD, Shaman J, Alonso WJ, Bloom-Feshbach K, Uejio CK, Comrie A, Viboud C. Environmental predictors of seasonal influenza epidemics across temperate and tropical climates. PLoS Pathog. 2013;9:e1003194.View ArticlePubMedPubMed CentralGoogle Scholar
- Finkelman BS, Viboud C, Koelle K, Ferrari MJ, Bharti N, Grenfell BT. Global Patterns in Seasonal Activity of Influenza A/H3N2, A/H1N1, and B from 1997 to 2005: Viral Coexistence and Latitudinal Gradients. PLoS One. 2007;2:e1296.View ArticlePubMedPubMed CentralGoogle Scholar
- Moura FEA, Perdigão ACB, Siqueira MM. Seasonality of influenza in the tropics: a distinct pattern in Northeastern Brazil. Am J Trop Med Hyg. 2009;81:180–3.PubMedGoogle Scholar
- Rao BL, Banerjee K. Influenza surveillance in Pune, India, 1978–90. Bull WHO. 1993;71:177–81.PubMedPubMed CentralGoogle Scholar
- Rao BL, Yeolekar LR, Kadam SS, Pawar MS, Kulkarni PB, More BA, Khude MR. Influenza surveillance in Pune, India, 2003. Southeast Asian J Trop Med Public Health. 2005;36:906–9.PubMedGoogle Scholar
- Dosseh A, Ndiaye K, Spiegel A, Sagna M, Mathiot C. Epidemiological and virological influenza survey in Dakar, Senegal: 1996-1998. Am J Trop Med Hyg. 2000;62:639–43.PubMedGoogle Scholar
- Fuhrmann C. The effects of weather and climate on the seasonality of influenza: what we know and what we need to know. Geography Compass. 2010;4:718–30.View ArticleGoogle Scholar
- Lowen AC, Steel J, Mubareka S, Palese P. High temperature (30 °C) blocks aerosol but not contact transmission of influenza virus. J Virol. 2008;82:5650–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowen AC, Mubareka S, Steel J, Palese P. Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathog. 2007;3:e151.View ArticlePubMed CentralGoogle Scholar
- McDevitt J, Rudnick S, First M, Spengler J. Role of absolute humidity in the inactivation of influenza viruses on stainless steel surfaces at elevated temperatures. Appl Environ Microbiol. 2010;76:3943–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaman J, Pitzer VE, Viboud C, Grenfell BT, Lipsitch M. Absolute humidity and the seasonal onset of influenza in the continental United States. PLoS Biol. 2010;8:e1000316.View ArticlePubMedPubMed CentralGoogle Scholar
- Żuk T, Rakowski F, Radomski JP. A model of influenza virus spread as a function of temperature and humidity. Comput Biol Chem. 2009;33:176–80.View ArticlePubMedGoogle Scholar
- Barreca AI, Shimshack JP. Absolute humidity, temperature, and influenza mortality: 30 years of county-level evidence from the United States. Am J Epidemiol. 2012;176:S114–S22.View ArticlePubMedGoogle Scholar
- van Noort SP, Águas R, Ballesteros S, Gabriela M, Gomes M. The role of weather on the relation between influenza and influenza-like illness. J Theor Biol. 2012;298:131–7.View ArticlePubMedGoogle Scholar
- Jaakkola K, Saukkoriipi A, Jokelainen J, Juvonen R, Kauppila J, Vainio O, Ziegler T, Rönkkö E, Jaakkola JJK, Ikäheimo TM, et al. Decline in temperature and humidity increases the occurrence of influenza in cold climate. Environ Health. 2014;13:1–8.View ArticleGoogle Scholar
- Chong KC, Goggins W, Zee BCY, Wang MH. Identifying meteorological drivers for the seasonal variations of influenza infections in a subtropical city - Hong Kong. Int J Environ Res Public Health. 2015;12:1560–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Lofgren E, Fefferman NH, Naumov YN, Gorski J, Naumova EN. Influenza seasonality: underlying causes and modeling theories. J Virol. 2007;81:5429–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Eccles R. An explanation for the seasonality of acute upper respiratory tract viral infections. Acta Otolaryngol. 2002;122:183–91.View ArticlePubMedGoogle Scholar
- te Beest DE, van Boven M, Hooiveld M, van den Dool C, Wallinga J. Driving factors of influenza transmission in the Netherlands. Am J Epidemiol. 2013;178:1469–77.View ArticleGoogle Scholar
- Soebiyanto RP, Clara W, Jara J, Castillo L, Sorto OR, Marinero S, de Antinori MEB, McCracken JP, Widdowson M-A, Azziz-Baumgartner E. The role of temperature and humidity on seasonal influenza in tropical areas: Guatemala, El Salvador and Panama, 2008–2013. PLoS One. 2014;9:e100659.View ArticlePubMedPubMed CentralGoogle Scholar
- Alonso WJ, Viboud C, Simonsen L, Hirano EW, Daufenbach LZ, Miller MA. Seasonality of influenza in Brazil: a traveling wave from the Amazon to the Subtropics. Am J Epidemiol. 2007;165:1434–42.View ArticlePubMedGoogle Scholar
- Mahamat A, Dussart P, Bouix A, Carvalho L, Eltges F, Matheus S, Miller MA, Quenel P, Viboud C. Climatic drivers of seasonal influenza epidemics in French Guiana, 2006–2010. J Infect. 2013;67:141–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Soebiyanto RP, Adimi F, Kiang RK. Modeling and predicting seasonal influenza transmission in warm regions using climatological parameters. PLoS One. 2010;5:e9450.View ArticlePubMedPubMed CentralGoogle Scholar
- Chumkiew S, Srisang W, Jaroensutasinee M, Jaroensutasinee K. Climatic factors affecting on influenza cases in Nakhon Si Thammarat. World Acad Sci Eng Technol. 2007;1:633–6.Google Scholar
- Helming L, Böse J, Ehrchen J, Schiebe S, Frahm T, Geffers R, Probst-Kepper M, Balling R, Lengeling A. 1α,25-dihydroxyvitamin D3 is a potent suppressor of interferon γ-mediated macrophage activation. Blood. 2005;106:4351–8.View ArticlePubMedGoogle Scholar
- Abu-Amer Y, Bar-Shavit Z. Impaired bone marrow-derived macrophage differentiation in vitamin D deficiency. Cell Immunol. 1993;151:356–68.View ArticlePubMedGoogle Scholar
- Cannell JJ, Vieth R, Umhau JC, Holick MF, Grant WB, Madronich S, Garland CF, Giovannucci E. Epidemic influenza and vitamin D. Epidemiol Infect. 2006;134:1129–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Urashima M, Segawa T, Okazaki M, Kurihara M, Wada Y, Ida H. Randomized trial of vitamin D supplementation to prevent seasonal influenza A in schoolchildren. Am J Clin Nutr. 2010;91:1255–60.View ArticlePubMedGoogle Scholar
- Dowell SF. Seasonal variation in host susceptibility and cycles of certain infectious diseases. Emerg Infect Dis. 2001;7:369–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiao H, Tian H, Lin X, Gao L, Dai X, Zhang X, Chen B, Zhao J, Xu J. Influence of extreme weather and meteorological anomalies on outbreaks of influenza A (H1N1). Chin Sci Bull. 2013;58:741–9.View ArticleGoogle Scholar
- Yaari R, Katriel G, Huppert A, Axelsen JB, Stone L. Modelling seasonal influenza: the role of weather and punctuated antigenic drift. J R Soc Interface. 2013;10:20130298.View ArticlePubMedPubMed CentralGoogle Scholar
- Davis RE, McGregor GR, Enfield KB. Humidity: a review and primer on atmospheric moisture and human health. Environ Res. 2016;144(Part A):106–16.View ArticlePubMedGoogle Scholar
- Roy M, Pascual M. On representing network heterogeneities in the incidence rate of simple epidemic models. Ecol Complex. 2006;3:80–90.View ArticleGoogle Scholar
- Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proc R Soc Lond A. 1927;115:700–21.View ArticleGoogle Scholar
- Nakagawa S. A farewell to Bonferroni: the problems of low statistical power and publication bias. Behav Ecol. 2004;15:1044–5.View ArticleGoogle Scholar
- Good P. Permutation tests: a practical guide to resampling methods for testing hypotheses. Verlag New York: Springer; 2000.View ArticleGoogle Scholar
- R Core team. R: A Language and Environment for Statistical Computing. 2014.Google Scholar
- Dray S, Dufour AB, Chessel D. The ade4 package-II: Two-table and K-table methods. R News. 2007;7:47–52.Google Scholar
- Dray S, Dufour A-B, et al. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22:1–20.View ArticleGoogle Scholar
- Chessel D, Dufour AB, Thioulouse J. The ade4 package-I-One-table methods. R News. 2004;4:5–10.Google Scholar
- Bates D, Maechler M, Bolker BM, Walker S. Fitting Linear Mixed-Effects Models using lme4. 2015.Google Scholar
- Bates D, Maechler M, Steven Walker BB: lme4: Linear mixed-effects models using Eigen and S4. 2015.Google Scholar
- Viboud C, Pakdaman K, Boëlle P-y, Wilson M, Myers M, Valleron A-J, Flahault A. Association of influenza epidemics with global climate variability. Eur J Epidemiol. 2004;19:1055–9.View ArticlePubMedGoogle Scholar
- Axelsen JB, Yaari R, Grenfell BT, Stone L. Multiannual forecasting of seasonal influenza dynamics reveals climatic and evolutionary drivers. Proc Natl Acad Sci U S A. 2014;111:9538–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaman J, Karspeck A. Forecasting seasonal outbreaks of influenza. Proc Natl Acad Sci U S A. 2012;109:20425–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Eubank S, Guclu H, Anil Kumar VS, Marathe MV, Srinivasan A, Toroczkai Z, Wang N. Modelling disease outbreaks in realistic urban social networks. Nature. 2004;429:180–4.View ArticlePubMedGoogle Scholar
- Lunelli A, Pugliese A, Rizzo C. Epidemic patch models applied to pandemic influenza: contact matrix, stochasticity, robustness of predictions. Math Biosci. 2009;220:24–33.View ArticlePubMedGoogle Scholar
- Balcan D, Gonçalves B, Hu H, Ramasco JJ, Colizza V, Vespignani A. Modeling the spatial spread of infectious diseases: the GLobal Epidemic and Mobility computational model. J Comput Sci. 2010;1:132–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Arkema JMS, Meijer A, Meerhoff TJ, Velden J, Paget WJ. Epidemiological and virological assessment of influenza activity in Europe, during the 2006-2007 winter. Euro Surveill. 2008;13:18958.PubMedGoogle Scholar
- Del Valle SY, Hyman JM, Hethcote HW, Eubank SG. Mixing patterns between age groups in social networks. Soc Netw. 2007;29:539–54.View ArticleGoogle Scholar
- Glass K, Mercer GN, Nishiura H, McBryde ES, Becker NG. Estimating reproduction numbers for adults and children from case data. J R Soc Interface. 2011;8:1248–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008;5:e74.View ArticlePubMedPubMed CentralGoogle Scholar