key: cord-0003468-vojbn0p4 authors: Hsieh, Ying‐Hen title: Pandemic influenza A (H1N1) during winter influenza season in the southern hemisphere date: 2010-07-08 journal: Influenza Other Respir Viruses DOI: 10.1111/j.1750-2659.2010.00147.x sha: 27274e02683a0348194933d6f5669741346b7113 doc_id: 3468 cord_uid: vojbn0p4 Please cite this paper as: Hsieh Ying‐Hen. (2010) Pandemic influenza A (H1N1) during winter influenza season in the southern hemisphere. Influenza and Other Respiratory Viruses 4(4), 187–197. Background Countries in the southern hemisphere experienced sizable epidemics of pandemic influenza H1N1 in their winter season during May–August, 2009. Methods We make use of the Richards model to fit the publicly available epidemic data (confirmed cases, hospitalizations, and deaths) of six southern hemisphere countries (Argentina, Brazil, Chile, Australia, New Zealand, and South Africa) to draw useful conclusions, in terms of its reproduction numbers and outbreak turning points, regarding the new pH1N1 virus in a typical winter influenza season. Results The estimates for the reproduction numbers of these six countries range from a high of 1·53 (95% CI: 1·22, 1·84) for confirmed case data of Brazil to a low of 1·16 (1·09, 1·22) for pH1N1 hospitalizations in Australia. For each country, model fits using confirmed cases, hospitalizations, or deaths data always yield similar estimates for the reproduction number. Moreover, the turning points for these closely related outbreak indicators always follow the correct chronological order, i.e., case–hospitalization–death, whenever two or more of these three indicators are available. Conclusions The results suggest that the winter pH1N1 outbreaks in the southern hemisphere were similar to the earlier spring and later winter outbreaks in North America in its severity and transmissibility, as indicated by the reproduction numbers. Therefore, the current strain has not become more severe or transmissible while circulating around the globe in 2009 as some experts had cautioned. The results will be useful for global preparedness planning of possible tertiary waves of pH1N1 infections in the fall/winter of 2010. The pandemic influenza H1N1 2009 (pH1N1) was first identified by the Centers for Disease Control in two children in April. 1 As of 14 February 2010, worldwide more than 212 countries and overseas territories or communities have reported laboratory-confirmed cases of pandemic influenza H1N1 2009, including at least 15921 deaths. 2 As predicted earlier by many experts (e.g., 3) , the southern hemisphere also witnessed a sizable epidemic of pH1N1 in its winter season during May-August, with attack rates exceeding those in a typical influenza season in these countries. Moreover, it is important to monitor changes in antigenicity, severity, transmissibility, and antiviral resistance in the southern hemisphere to compare with those of the spring and winter outbreaks in the northern hemisphere. We construct simple mathematical models using the epidemiologic data for these southern hemisphere countries to draw timely and useful conclusions regarding the role played by the new pH1N1 strain under the setting of a typical winter influenza season. Data We obtain the 2009 pH1N1 data for Argentina, Chile, and Brazil used in this study from the respective Ministries of Health of Argentina, 4 Chile, 5 and Brazil 6 websites. The New Zealand data were accessed from the Institute of Environmental Science and Research of New Zealand website. 7 The Australian data were obtained from Department of Health and Ageing website. 8 The South Africa data were accessed from the National Institute for Communicable Diseases (NICD) website. 9 Methods Richards 10 proposed the following model to study the growth of biologic populations: Here, the prime '¢' denotes the derivative or time rate of change and the time unit is in days or epidemiological weeks (or e-weeks that start on Sunday and end on Saturdays), depending on the time unit of the data used. C(t) is the cumulative number of cases (confirmed, hospitalization, or death) at time unit (day or week) t, K is the 'maximum' case number or the final outbreak size over a single wave of outbreak, r is the per capita growth rate of the cumulative number of cases, while a is an exponent of deviation of the epidemiologic curve. The solution of Richards model can be explicitly given in terms of its model parameters and the initial value C(0) as CðtÞ ¼ K½1 þ e ÀrðtÀt m Þ À1=a . The parameter t m is related to the turning point t i (defined as the time when the rate of case accumulation changes from increasing to decreasing or vice versa) of the epidemic by the simple formula t m = t i + (ln a) ⁄ r, where ln denotes the natural logarithm function. By fitting cumulative case data of a particular outbreak, these four quantities, namely, the growth rate r, carrying capacity K, inflection point t i , and exponent of deviation a of the model, can be obtained simultaneously using standard software such as MATLAB or SAS. In this model formulation, the basic reproduction number R 0 is given by the formula R 0 ¼ exp ðrTÞ, where T is the disease generation time (or generation interval) defined as the average time interval from one person being infected to the time when an infection by this infected individual occur. The estimation of mean generation time could be difficult because it is often impractical to report observed generation time, and sophisticated statistical procedures have been developed to attain this goal (see e.g., 11, 12 ) . Moreover, it has been shown mathematically 11 that, given the growth rate r, the expression R 0 ¼ exp ðrTÞ provides the upper bound of basic reproduction number regardless of the distribution of the generation interval that is being used. In the case where there is little prior human immunity for the virus under consideration, the estimate for R 0 approximates the number of infections caused by an infectious individual entering an immunologically naïve population. Readers are referred to [13] [14] [15] for more technical details regarding the Richards model. The basic premise of the Richards model is that the incidence curve of a single wave of infections consists of a single peak of high incidence, resulting in an S-shaped cumulative epidemic curve with a single turning point for the outbreak. The turning point t i can be easily pinpointed by locating the inflection point of the cumulative case curve, i.e., the moment at which the trajectory begins to decline. This quantity in time has obvious epidemiologic importance, indicating either the beginning of a new wave of infections (i.e., moment of acceleration after deceleration) or the end of the current wave of infections (i.e., moment of deceleration after acceleration). One of the advantages in using the Richards model for mathematical modeling is to fit the accumulative case number, which helps to smooth out the stochastic variations in epidemic curve owing to variations in data collection. When there are indeed two or more waves of infections, a variation of the S-shaped Richards model has been proposed, 14 which distinguishes between two types of turning points: one in the initial S curve that signifies the first turning point that ends initial exponential growth (or a downturn in case number) and a second type of turning point in the epidemic curve where the growth rate of the number of cumulative cases begin to increase again (or a upturn in case number), signifying the beginning of the next wave of infections. This variation of the Richards model provides a systematic method of determining whether an outbreak is 1-wave or multi-wave in nature, and is designed to isolate the main turning points from those peaks and valleys resulting from the data stochasticity. More details on the multi-wave Richards model can be found in, [14] [15] [16] where 2003 SARS incidence curves for the Great Toronto Area and Singapore are shown to contain two peaks (local maximum or turning point of the first type) and one valley (local minimum or turning point of second type). The readers are also referred to 14 for its application to 2-wave dengue outbreak in Taiwan in 2007. We note that the estimates of R 0 should be valid for the new H1N1 strain, although discrepancy might exist because of prior immunity. 17 Moreover, no mass intervention measures (e.g., community-level isolation ⁄ quarantine ⁄ closure) were implemented at any time during the course of the outbreaks, although some measures were taken in each country with varying degrees of impact. To be conservative in our conclusions, we will henceforth denote the reproduction number obtained here as the effective reproduction number R of the outbreak in question. We fit the cumulative pandemic influenza A (H1N1) 2009 case data from various southern hemisphere countries described earlier to the Richards model and its multi-wave variants. The country-wise results of model parameter estimation, using least squares techniques, are detailed in the following text. For the purpose of computing R 0 , we make use of the mean estimated generation interval (and its 95% CI) of T = 1AE91 days (95% CI: 1AE30-2AE71) as given in, 18 which was estimated from early Mexico novel H1N1 data before April 30. The daily 2009 pH1N1 confirmed case, hospitalization [of all severe respiratory infections or sever respiratory infection (SRI)], and pH1N1 death data of Argentina 4 during May 19 to August 23 are fitted to the Richards model. The results of model parameter estimation for Argentina are given in Tables 1 and 2. Table 1 shows that the daily confirmed pH1N1 case data fit both 1-wave and 2-wave Richards models. Using the daily confirmed case data, the 1-wave model fit gives a turning point of June 26, while the 2-wave yields 3 turning points, one 'downturn' tuning point for each wave around June 1 and June 24, and an 'upturn' turning point on June 4 signifying the start of the second wave. For comparison, both the daily SRI hospitalization data (6 ⁄ 1-7 ⁄ 25) and daily pH1N1 death data (6 ⁄ 12-8 ⁄ 20) fit 1-wave Richards models ( Table 2 ). The turning point using daily SRI hospitalization data is June 28, while the turning point for the daily H1N1 death data is July 1. For the purpose of comparison with the estimates for the final case number K for each wave of outbreak, we also provide the actual case numbers in the table legends. The estimates for the reproduction number, using the daily confirmed pH1N1 case data, are R = 1AE41 (95% CI: 1AE23, 1AE60), and 1AE84 (95% CI: 0AE86, 2AE82) and 1AE34 (95% CI: 1AE19, 1AE48), respectively, for the 2-wave model. The corresponding estimates are R = 1AE27 (95% CI: 1AE16, 1AE38) and 1AE28 (95% CI: 1AE16, 1AE39), respectively, using SRI hospitalization and pH1N1 death data. The graphs of the theoretical curves are shown in Figure 1A -D. The weekly pH1N1 hospitalization data by onset week of eweek 20-33 (May 10 to August 15) and daily death data of June 1 to August 18 in Chile 6 fit 1-wave Richards model (Table 3 and Figure 1E -F). The turning point for weekly pH1N1 hospitalization data is e-week 26, whereas the turning point for the daily death data is July 11. The estimated reproduction numbers are 1AE24 (95% CI: 1AE14, 1AE35) for weekly pH1N1 hospitalization data and 1AE29 (95% CI: 1AE17, 1AE41) for the daily death data. Note that the confirmed pH1N1 case data for Chile did not converge for earlier Ministry of Health of Chile Pandemic Virus Report reported in July, most likely attributable to the transient state of data collection during earlier stages of the outbreak (e.g., changing case definitions and level of surveillance), and was not given in the Pandemic Influenza Report in August 5 and onward. The weekly pH1N1 confirmed case data by onset week of e-week 21-33 (May 17 to August 15) in Brazil 6 also fits both 1-wave and 2-wave Richards models (Table 4 ). The turning point for 1-wave model is e-week 31, while the two respective 'downturn' turning points for the two waves are e-week 26 and e-week 31. The estimated reproduction numbers are 1AE53 (95% CI: 1AE22, 1AE84) for 1-wave model, and 1AE66 (95% CI: 1AE23, 2AE09) and 1AE39 (95% CI: 0AE99, 1AE79) for the respective two waves of the 2-wave model. The graph for model fit is given in Figure 1G . We first note that the New Zealand data are given in calendar weeks ending on Sundays, 7 as opposed to e-weeks that end on Saturdays. The weekly pH1N1 confirmed case data by onset week of weeks ending on April 12-September 20 and weekly pH1N1 hospitalization data by hospitalization, onset, or report week of weeks ending on June 7-September 13 in New Zealand 7 fit 1-wave Richards models (Table 5 and Figure 2A-D) . Moreover, the sentinel and non-sentinel surveillance data by reporting week of weeks ending on June 14-September 20 and May 3-September 20, respectively, also fit 1-wave Richards models. The turning point for the weekly pH1N1 confirmed case data is the week ending July 5 and for the hospitalization data is the succeeding week ending July 12, giving strong evidence that the turning point for the better had occurred by early July. The turning point for the sentinel and non-sentinel surveillance data are the weeks ending July 19 and July 12, respectively. The estimated reproduction numbers are 1AE22 (95% CI: 1AE13, 1AE31) for the weekly confirmed case data; 1AE29 (95% CI: 1AE16, 1AE41) for weekly pH1N1 hospitalization data; 1AE18 (95% CI: 1AE10, 1AE25) for the sentinel; and 1AE22 (95% CI: 1AE13, 1AE31) for the non-sentinel surveillance data. The weekly pH1N1 confirmed case and hospitalization data by onset week in Australia 8 fit the 1-wave Richards model (Table 6 and Figure 2E -F). The turning point for both the weekly pH1N1 confirmed case data and the hospitalization data is e-week 30, showing that the down turning point of case number had occurred by late July. The estimated reproduction numbers are 1AE19 (95% CI: 1AE11, 1AE26) for the weekly confirmed case data; 1AE16 (95% CI: 1AE09, 1AE22) for weekly pH1N1 hospitalization data. The weekly pH1N1 confirmed case data of weeks ending on July 5 to September 20 in South Africa 9 fit the 1-wave Richards model (Table 7 and Figure 2G ). The turning point for the weekly pH1N1 confirmed case data is the week ending August 16. The estimated reproduction number is 1AE22 (95% CI: 1AE11, 1AE33). Several alternate estimates of the generation time T have been reported. Hence, we perform sensitivity analysis on R for Argentina and New Zealand data using estimates of T obtained from pH1N1 data in US: T = (2AE6, 3AE2) 19 and T = 2AE5 20 ; in Netherlands: T = 2AE7 (SD = 1AE1); 21 and in Japan: T = 2AE83 (Variance = 1AE26) [Nishiura, personal communications]. The resulting estimates for R are given in Tables 8 and 9 . Reproduction number 23 as well as numerous estimates for spring ⁄ summer outbreaks of other countries in literature (see Table 2 , 23). The estimates for R in Argentina and New Zealand obtained by using different estimates of generation interval (T) for 2009 pH1N1 in the literature are also in close agreement (Tables 8 and 9 ). The results indicate that higher values of T will raise all the estimates for R, as can be expected from the formula for R given earlier. For the Argentina hospitalization data with T = 1AE91 18 , R = 1AE27 is obtained; whereas using a T value of around 2AE7 21 we have R = 1AE40 (see Table 8B ). Both estimates are within reasonable range. However, the discrepancy becomes larger for larger values of R (see next to last row, Table 8A ), thus indicating the dependence of R on accurate estimate of T. The corresponding estimates for other countries are similar and hence are omitted for brevity. The estimates for New Zealand are substantially lower than the range of 1AE49-2AE55 for reproduction number obtained in, 24 using the case data of 585 confirmed cases and 38 probable cases from 28 May to 28 June. We note, however, that 63 imported cases were first removed for the data in. 24 This omission most likely resulted in an overestimation of the initial growth rate, by ignoring these imported cases that must have collectively played a significant role in the early local transmissions in the island nation. We note that the role of imported cases is problematic to quantitate when estimating transmissibility, as there are no 'local' source of infection for these imported cases. However, by discounting the imported cases, which in the case of New Zealand data accounts for more than 10% of total cases, one ignores the infections caused by the imported infective and hence inevitably overestimates the transmissibility of the remaining local cases. There are several factors relating to imported cases that might impact the estimation of R, such as the timing and the size of the imported cases, the nature of the resulting local infections (e.g., household, community, nosocomial), and not the least of all the epidemiological characteristics of the disease under consideration. For example, a few imported cases at the early stages that set off a fast-spreading disease locally (e.g., pandemic influenza) would have a decidedly different effect when compared to the situation with imported cases scattered over a long period of time for a less transmissible disease with no apparent asymptomatic infections such as SARS, unless there are very fastdeveloping local clusters of infections or nosocomial spread, as in the 2003 SARS outbreak. A study on pH1N1 epidemic in Victoria, Australia, indicates that, by using methods designed to account for the undetected transmission caused by unidentified imported cases, the mean estimate for R reduces dramatically from 2AE4 to 1AE6. 26 The reproduction number obtained for Brazil is only slightly higher than other southern hemisphere countries studied (see Figure 3 ), but with a much larger confidence interval. In the cases of Argentina and Brazil where both 1-wave and 2-wave Richards models fit the cumulative case data, R is always larger in the first phase in May to June, perhaps indicative of the initial severity or lack of immunity ⁄ intervention during early stages of a novel virus outbreak. We note that higher R for earlier data was also observed in, 22 after correction for ascertainment bias. It is particularly interesting to note that for countries such as Argentina, Chile, New Zealand, and Australia, where multiple estimates were obtained from distinct epidemic data (case, hospitalization, and deaths), the estimates of R are highly comparable. Therefore, in events where the disease epidemic curve is not readily available or does not converge (as in the case of Chile), the hospitalization and deaths data, while not necessarily accurate for detecting the Among the South American countries, the downturn turning point of the outbreak occurred earlier in Argentina and Chile in June, while Brazil did not experience its downtown until late July. New Zealand (early July) had its turning point slightly ahead of Australia (late July). For South Africa the outbreak did not turn for the better until mid-August. Moreover, in countries where multiple epidemic data (cases, hospitalizations, and deaths) are available for model fit, the resulting turning points for reporting of cases, hospitalizations, and deaths always follow the proper chrono- logical order of disease progression (Table 10 and Figure 4) , including that of Australia where the turning points for reporting of cases and hospitalizations fall on the same e-week. This gives further credence to the detection of turning points of an outbreak, not only for disease incidence but also for the temporal changes in disease progression and treatment of the infected individuals. We note that in general the turning points of an outbreak correspond to peaks and valleys in the epidemic curve. However, as a result of stochastic variations, multiple spikes often appear around the peaks. Moreover, some temporary deceleration or acceleration in the cumulative rate might not be related to a declining epidemic. The Richards model has the advantage of isolating a small range of time period during which the turning points (up turns and down turns) may have occurred. [13] [14] [15] [16] 25 We also note that the estimates for the outbreak sizes, K, also closely approximate the actual case numbers for the periods of Table 9 . Estimates of effective reproduction number R for New Zealand data using estimates of T obtained from pH1N1 data of: US with (a) T = (2AE6, 3AE2) 19 14, 16 and hence are not repeated here. In addition to existing pre-immunity, estimates of R may also be influenced by factors such as social mixing patterns, access to medical care, and level of intervention measures in each country. To address this issue, model specific for each country is required to account for the social mixing and healthcare infrastructure of each country, as well as any specific interventions that were implemented, but is beyond the scope of this work. Numerous sophisticated models have been proposed (e.g., 27-30) that give time-dependent reproduction number R, which could change with time owing to interventions and ⁄ or susceptible depletion but requires more detailed and high-quality data. Nevertheless, the Richards model provides a swift and viable method to quick estimate of R as well as identifying the important turning points of the outbreak, with a minimal requirement of epidemic data that are typically available to the public. It is interesting to note that, for Argentina, the hospitalization data for confirmed pH1N1 cases are not available on the web; therefore, we used the hospitalization data for severe respiratory infections, which include but are not limited to the confirmed pH1N1 cases. However, the resulting estimate for reproduction number is surprisingly similar to the estimates using the case data or the death data. Moreover, the estimated turning point for hospitalizations is also consistent with those of cases and deaths. Most of the southern hemisphere countries studied are of substantial size, and thus the epidemic had different timings in different regions. For example, in Australia, the epidemic in Victoria was a full month ahead of some other states might well have impacted the epidemic in other states. In such cases, instead of using national totals, the detailed data for each region are needed to account for the regional heterogeneity. However, other factors such as the level of mobility between regions must then be considered, which would call for a complicated model with spatial elements, which is again beyond the scope of this study. In summary, our results suggest that the winter pandemic A (H1N1) outbreak in the southern hemisphere, as indicated by its reproduction number, was not dissimilar to the earlier spring epidemics 23 or the fall ⁄ winter outbreaks in the northern hemisphere based on preliminary results of an ongoing study of the fall ⁄ winter outbreak in Taiwan using the same modeling approach. Therefore, the current strain, while circulating around the globe, has not become more severe or transmissible either in the southern or in the northern hemisphere, as some experts had cautioned. The results, which corroborate earlier findings that the 2009 pandemic had not been severe in the southern hemisphere, 31 will be useful for global preparedness planning of possible tertiary waves of pH1N1 infections in the fall ⁄ winter of both southern and northern hemispheres in 2010. Swine Influenza A (H1N1) Infection in Two Children -Southern California Pandemic (H1N1) 2009 -update 88 Managing and Reducing Uncertainty in an Emerging Influenza Pandemic INFLUENZA PANDÉ MICA (H1N1) Brazilian Ministry of Health Influenza Surveillance Report No. 19, for the week ending 18 The National Institute for Communicable Diseases (NICD) of South Africa. NICD Weekly Situation Report A flexible growth function for empirical use How generation intervals shape the relationship between growth rates and reproductive numbers The ideal reporting interval for an epidemic to objectively interpret the epidemiological time course Real-time forecast of multi-wave epidemic outbreaks Richards Model: a simple procedure for real-time prediction of outbreak severity Turning points, reproduction number, and impact of climatological events on multi-wave dengue outbreaks Understanding influenza transmission, immunity and pandemic threats Pandemic potential of a strain of Influenza A (H1N1): early findings The transmissibility and control of pandemic influenza A (H1N1) virus Estimation of the reproductive number and the serial interval in early phase of the 2009 influenza A ⁄ H1N1 pandemic in the USA Epidemiology and control of influenza A(H1N1)v in the Netherlands: the first 115 cases Initial human transmission dynamics of the pandemic (H1N1) 2009 virus in North America Transmission dynamics and impact of pandemic infl uenza A (H1N1) 2009 virus Estimating the reproduction number of the novel influenza A virus (H1N1) in a Southern Hemisphere setting: preliminary estimate in New Zealand Intervention measures, turning point, and reproduction number for dengue Early transmission characteristics of influenza A(H1N1)v in Australia: Victorian state Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures Estimating in real time the efficacy of measures to control emerging communicable diseases Real Time Bayesian estimation of the epidemic potential of emerging infectious diseases A likelihood-based method for real-time estimation of the serial interval and reproductive number of an epidemic Pandemic H1N1 influenza lessons from the southern hemisphere YHH is supported by grants (NSC 97-2314-B-039-013-MY3, NSC 97-2118-M-039-004, and CMU 97 323). The author is grateful to the reviewers for their insightful comments and suggestions, which significantly improved this manuscript.