Reference scenario: 2009 H1N1 influenza pandemic in the absence of interventions
In the absence of interventions, the simulated epidemics show a mean peak weekly incidence of 8.9 symptomatic cases (95 % CI: 1.6–13.0) per 10,000 individuals and an overall infection attack rate of 19.2 % (95 % CI: 17.9–20.6 %), see Fig.1a and c. In model simulations, the number of weekly cases is correlated with student absenteeism, reaching a peak of 9.6 % absent students on average on the same week when the epidemic shows its maximum incidence (Fig. 1b). Simulated infection attack rate (Fig. 1d) is highly dependent on age and compares well with serosurvey data on the 2009 H1N1 influenza pandemic in Italy [17, 22]: school-age individuals represent clearly the most affected portion of society (with 56.5 % of individuals aged 6–18 years affected), while individuals aged 65 years or older represent the less affected age group (with an infection attack rate lower than 10 %).
Effectiveness of school closure strategies
Our results highlight that, independently of the number of theoretical weeks lost, the identification of the optimal excess absenteeism values triggering closures is a major determinant of policy effectiveness (see Fig. 2). Specifically, we found that closing the entire school (strategy S) when the optimal excess absenteeism threshold is adopted leads to about 4.4, 9.9, 13.8, and 19.4 % attack rate reduction for one, two, three, and four theoretical weeks lost, respectively (Fig. 2a) – an assessment of the variability of these estimates is provided in the Additional file 1: SM. Higher reductions are observed for gradual strategies (strategies CS and CGS), which are both capable of leading to reductions of about 8, 17, 23, and 25 % for one, two, three, and four theoretical weeks lost, respectively (Fig. 2a). Such a difference between the strategies is imputable to the effective number of days lost by students; in fact, as in gradual strategies also classes, and possibly grades, can remain closed without the need of closing the entire school, the actual amount of school days lost by students due to the application of strategies can be higher than the maximum number of weeks that a school can remain closed (see Fig. 2b). In particular, we found that attack rate reduction is strongly correlated with the actual number of school days lost per student (correlation 0.99, p-value < 0.0001, for all strategies), and a linear model with regression coefficient 1.12 is capable to well explain the variation (coefficient of determination 0.99, p-value < 0.0001), which means that an increase of one (effective) school day lost per student corresponds to a decrease in the attack rate of about 1.1 %. Such a pattern is consistent for the three strategies (see Additional file 1: SM).
All evaluated strategies are based on monitoring student absenteeism irrespective to the cause; however, an important distinction has to be made: strategy S is based on the monitoring of student absenteeism in the whole school, while strategies CS and CGS are based on the monitoring of absenteeism in classes, which are at least one order of magnitude smaller than schools. Clearly, the closure of single classes as a consequence of an increase in absenteeism not related to influenza yields a more acceptable burden than the closure of the entire school. On the other hand, this implies that the same percentage of absent students triggering closure corresponds to a much higher number of absent students for strategy S than for gradual strategies, suggesting that small variations in the excess absenteeism threshold would have a higher impact in strategy S than in the others. This introduces a second important aspect that has to be taken into account when considering school closure policies as intervention options: the stability of the effectiveness of the strategies with respect to the choice of the absenteeism threshold. In fact, a highly effective strategy that however requires a very precise level of absenteeism to be effective would be hardly applicable in real world situations. We consider a strategy to be stable if the range of excess absenteeism thresholds yielding at least 80 % (or 90 %) of the maximum reduction possible is wide. As clearly highlighted in Fig. 2a, strategy S has a low level of stability, while strategies CS and CGS are remarkably more stable: an excess absenteeism in the range 3.3–12.5, 0–9.5, 0–8.2, and 0–7.8 % (strategy CS) and 2.7–12.4, 0–9.9, 0–8.5, and 0–7.8 % (strategy CGS) leads to an attack rate reduction close to that obtained for the optimal value for one, two, three, and four theoretical weeks lost, respectively. These values also suggest that, in case decision makers were inclined to keep the social cost low, it would be more effective to close schools (and classes) when there is less uncertainty on the cause of student absenteeism.
All evaluated strategies are highly effective in reducing the peak weekly incidence: an increase of one (effective) school day lost per student corresponds to a decrease of the peak week incidence of about 2.6 % (see Additional file 1: SM). In particular, for one theoretical week lost, we estimate a maximum peak week incidence reduction of about 19.6 % for strategy S, and of about 34.5 % for both gradual strategies; for two or more theoretical weeks lost, the estimated maximum reduction reaches about 39–48 % for strategy S, and about 51–56 % for both gradual strategies (Fig. 2c). Again, gradual strategies show a much higher level of stability of the results with respect to the threshold of excess absenteeism than strategy S. Finally, we also found that, except for two theoretical weeks lost, the optimal threshold levels of excess absenteeism for reducing attack rate are similar to those for reducing peak incidence, highlighting that the analyzed school closure policies are highly beneficial for mitigating influenza spread in several ways. In the case of two theoretical weeks lost, the highest reduction in peak week incidence is estimated to occur when the policies are applied in the early phase of the epidemic (with all schools closing almost simultaneously), while for reducing the attack rate it would be preferable to have less uncertainty on the causes of student absenteeism by waiting to reach a higher excess absenteeism.
Finally, results suggest that the mitigation effects of strategies tend to vanish for increasingly high values of the excess absenteeism; this is due to the growing difficulty in reaching such high values of absenteeism, which results in poor (or failure in the) application of strategies.
Sensitivity of school closure strategies to intervention parameters
So far, we have shown that the effectiveness of the three strategies strongly depends on the number of theoretical weeks lost and on the threshold in student absenteeism. However, there are other parameters whose effect on the regulation of strategies needs to be investigated. In particular, in order to limit the effort required to each school (e.g., in terms of logistic organization), decision makers may be inclined to consider 2 or 4 weeks lost per student as an affordable cost, but only in form of a single closure episode. We show that, assuming a fixed number of theoretical weeks lost, results are quite stable for all strategies with respect to changes in the number of closure events (see Fig. 3A and Additional file 1: SM), provided that the proper threshold for the excess absenteeism of students is chosen (Fig. 3B). In general, we found that single closure episodes achieve better results for higher threshold levels (i.e. when there is less uncertainty on the cause of such high absenteeism rates), since longer closures are able to disrupt disease transmission more effectively.
For all strategies, results remain very stable also with respect to changes in the start of student absenteeism monitoring in schools and, consequently, in the application of school closure policies (see Fig. 3A and Additional file 1: SM). In particular, we found slightly better results when starting to monitor student absenteeism after a higher incidence (i.e., more than 75 weekly symptomatic cases per 100,000 individuals) in the total population has been observed, suggesting again that the three strategies are more effective if applied when there is low uncertainty that an epidemic is under way. Such finding is confirmed by the optimal value for the threshold on student absenteeism, which generally increases as the trigger incidence decreases (see Fig. 3B). In sum, these findings suggest that monitoring influenza-like illness (ILI) incidence - which collects information that is routinely reported to surveillance systems of health authorities - to ascertain the presence of an ongoing influenza epidemic would be sufficient as trigger for starting the monitoring of student absenteeism in schools.
Since background student absenteeism is variable by school type (i.e., larger in secondary than in primary schools), by considering excess absenteeism to determine school (and class) closure, we are assuming different absolute levels of absenteeism for each school type. Nonetheless, our findings show that this choice does not undermine the mitigation potential of the three strategies; indeed, results remains consistent if we consider a unique threshold on student absenteeism in all types of schools instead of excess absenteeism (Fig. 3A and B).
Sensitivity of school closure strategies to biological and behavioral characteristics
Epidemic spread and effectiveness of intervention options largely depend on several factors (such as reproduction number and severity) that, unlike parameters regulating intervention strategies, are not under control of decision makers. Our findings highlight the predominant role of reproduction number and age-specific susceptibility to infection in determining the effectiveness of school closure policies. In particular, as shown in Fig. 4A, we found that, although school closure is highly effective in a 2009 pandemic scenario, for remarkably higher reproduction numbers (effective reproduction number of 1.8) none of the three strategies is capable to effectively reduce the attack rate (reductions result to be lower than 5 %). This suggests that in the case of highly transmissible influenza viruses, combined strategies relying on several interventions (e.g., vaccination as well as household quarantine) would be required to mitigate the spread of influenza pandemic/epidemic [32–34]. On the other hand, by fixing the value of the reproduction number and varying the age-specific susceptibility to infection, we found that the more susceptible to infection children are with respect to adults, the more effective school closure policies become. In particular, in a scenario where children are twice as susceptible to infection as in the 2009 pandemic, at three theoretical weeks lost strategies CS and CGS entail an attack rate reduction of about 50 % (Fig. 4A) and a peak week incidence reduction larger than 50 % (see Additional file 1: SM). Again, gradual strategies CS and CGS result to be more effective than strategy S and remarkably less sensitive to the excess absenteeism threshold, consequently allowing for a more manageable application in real-world situations (see Fig. 4A and B).
The effectiveness of school closure policies appears to be much less affected by changes in the symptomatic probability and by variations in the transmission rate of symptomatic infectious individuals with respect to asymptomatic ones: comparable attack rate reductions are observed for all three strategies (Fig. 4A). Finally, we found that gradual strategies are more effective also if contacts between students at school are homogeneous, suggesting that the trigger on classes works better than that on the whole school (see Fig. 4A). In fact, more homogenous contacts within school entail a larger dispersion of the infection inside the school and thus a larger number of classes to be closed simultaneously by the strategies and, in turn, a more immediate closure of the entire school as a reaction to influenza spread (see Additional file 1: SM); on the other hand, closures of single classes (possibly asynchronous) prevail if transmission is related only to contacts within the class, allowing the infection to circulate in non-closed classes and hence propagate to the general population. For the two intermediate scenarios, one considering the default within-school contact patterns resulting from data analysis [18, 27] and one where classes, grades and the overall school have the same weight in infection transmission, we found that the different school closure policies yield very similar effectiveness. Strategy S is generally not sensitive to variations in within-school contact patterns (Fig. 4A).
Comparison with targeted school closure strategies for the UK
We compare our results with those presented in Fumanelli et al. [16] relative to the UK. The two studies differ in three substantial modeling aspects: first, our model accounts for data on time use and on face-to-face contacts in a primary school (see Section Methods), while in [16] a simpler model based on a static network of contacts between individuals was employed; second, we investigate the effectiveness of school closure strategies at the level of a single municipality, while [16] was focused on the national scale; third, the analyzed country is different. Clearly, considering two different countries requires resorting to nation-specific data, for instance on the socio-demographic characteristics of the population and thus on contact pattern, background absenteeism, school size and structure. In addition, although both works consider the 2009 H1N1 pandemic as the reference scenario, the different characteristics of the two countries give origin to two different epidemics even in the absence of interventions (e.g., R0 = 1.4 here, R0 = 1.5 in [16]). However, we are considering the same pandemic scenario and strategies S and CGS are implemented the same way, thus a comparison between the two is possible. We find that attack rate reductions are comparable in case of strategy S and even more similar for strategy CGS (Fig. 5a); moreover, both models show that strategy CGS is much more stable than strategy S with respect to the choice of the absenteeism threshold (Fig. 5b), suggesting that the gradual closure approach may be preferable in real-world situations. In both models the optimal excess absenteeism threshold decreases as the number of weeks of closure increases (Fig. 5b); its value differs between the two models probably as a consequence of the different average level of background absenteeism, average school and class size of the two countries.