{'fraction_non_alphanumeric': 0.023340627279358133, 'fraction_numerical': 0.017852801222604287, 'mean_word_length': 5.198493003229279, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
National Influenza Surveillance in the Philippines from 2006 to 2012: seasonality and circulating strains

Background The results of routine influenza surveillance in 13 regions in the Philippines from 2006 to 2012 are presented, describing the annual seasonal epidemics of confirmed influenza virus infection, seasonal and alert thresholds, epidemic curve, and circulating influenza strains. Methods Retrospective analysis of Philippine influenza surveillance data from 2006 to 2012 was conducted to determine seasonality with the use of weekly influenza positivity rates and calculating epidemic curves and seasonal and alert thresholds using the World Health Organization (WHO) global epidemiological surveillance standards for influenza. Results Increased weekly influenza positive rates were observed from June to November, coinciding with the rainy season and school opening. Two or more peaks of influenza activity were observed with different dominant influenza types associated with each peak. A-H1N1, A-H3N2, and two types of B viruses circulated during the influenza season in varying proportions every year. Increased influenza activity for 2012 occurred 8 weeks late in week 29, rather than the expected week of rise of cases in week 21 as depicted in the established average epidemic curve and seasonal threshold. The intensity was severe going above the alert threshold but of short duration. Southern Hemisphere vaccine strains matched circulating influenza virus for more surveillance years than Northern Hemisphere vaccine strains. Conclusions Influenza seasonality in the Philippines is from June to November. The ideal time to administer Southern Hemisphere influenza vaccine should be from April to May. With two lineages of influenza B circulating annually, quadrivalent vaccine might have more impact on influenza control than trivalent vaccine. Establishment of thresholds and average epidemic curve provide a tool for policy-makers to assess the intensity or severity of the current influenza epidemic even early in its course, to help plan more precisely resources necessary to control the outbreak. Influenza surveillance activities should be continued in the Philippines and funding for such activities should already be incorporated into the Philippine health budget. Electronic supplementary material The online version of this article (doi:10.1186/s12879-016-2087-9) contains supplementary material, which is available to authorized users.


Background
The World Health Organization recommends that seasonal and alert thresholds and average epidemic curves be established in countries as tools for early detection of influenza outbreaks to help control annual influenza epidemics [1,2]. Seasonal influenza epidemics occur all over the world and cause substantial economic burden through health care costs and absenteeism [3]. Although surveillance is now routinely conducted in many countries providing seasonality data and recommendations for timely vaccination [4,5], not all countries have established influenza thresholds or epidemic curves. Developed countries like the United States (US), United Kingdom (UK), most European countries, Australia, and New Zealand [6][7][8][9][10] are already using thresholds for influenza surveillance calculated through various methods [11,12].
In contrast, few countries in Asia have established their own thresholds and epidemic curves to determine if an influenza outbreak has started while concurrently monitoring influenza activity [2]. The purpose of this article is to report the results of Philippine National Influenza Surveillance from 2006 to 2012, particularly on the establishment of seasonal and alert thresholds, and average epidemic curve according to methods in the WHO manual [1], define the seasonality of influenza, and describe the circulating influenza strains.

Study design
We retrospectively analyzed weekly trends in influenza virus activity using Philippine National Influenza Surveillance data from 2006 to 2012.

Surveillance site
The Philippines is located within the latitude and longitude of 13°00′ N, 122°00′ E [13]. The Philippines includes 300,000 km 2 of total land area and is divided into three island groups: Luzon, Visayas, and Mindanao. It is administratively divided into 17 regions. It has a tropical climate with the rainy season from June to November and the dry season from December to May [14].

Establishment of influenza sentinel sites
In 2004, the Research Institute for Tropical Medicine (RITM) became a National Influenza Center (NIC) and received a grant from the Centers for Disease Control and Prevention (CDC) in the United States (US), to establish a National Influenza Surveillance Network. Surveillance activities started in ten sentinel sites in five regions in June 2005. From 2006 to October 2008, the sentinel sites had increased to 36 distributed in 13 out of 17 regions in the country. There were 18 health centers and 18 Outpatient Departments (OPD) of tertiary hospitals. Eight of the health centers were located in the National Capital Region (NCR) (Fig. 1). In November 2008, a CDC-funded influenza Burden of Disease study (BOD) [15] was set up with 16 health centers in the Cordillera Autonomous Region (CAR) (Fig. 1). The NIC asked the investigators to contribute influenza virus infection data to the national influenza surveillance. There were now 52 sentinel sites by 2009. Sentinel sites were chosen based on the high number of clinic consultations among the health centers in the city, a high population density, and the region being a migratory pathway of birds mixing with local population of fowls and poultry.

Data collection Clinical and epidemiologic data on influenza-like-illness
Clinical and epidemiologic data were collected in health centers and outpatient departments or admissions in hospitals designated as influenza surveillance sites, from patients with influenza-like-illness (ILI). The health center is the first line of primary health care in the Philippines, and is located in the community it is serving with a typical population of 15,000 to 250,000. One physician, one nurse, one midwife, one medical technologist and one sanitary inspector usually man the health center. It is responsible for the delivery of the health programs of the Department of Health. These health services include but are not limited to morbidity consultations and management, expanded program of immunization, prenatal care, family planning services and such disease control programs like tuberculosis, diarrhea, and acute respiratory infections. As the first line of primary health care, sick individuals who patronize the services at the health centers are either mildly to severely sick. Severely ill patients consulting the center are immediately referred to the hospitals.
Data collected were written into case report forms and included signs and symptoms of cough, fever, sore throat, rhinitis or colds, and history of travel and treatment received for illness before consulting the health center. ILI was defined as any patient with history of fever of sudden onset or measured temperature of >37.7°C and cough or sore throat [16]. In children ≤3 years of age, ILI was defined only as fever and colds. "Colds" in these children was synonymous with runny nose. Sentinel sites were visited twice a week (Monday and Tuesday) to collect nasopharyngeal (NP) or oropharyngeal (OP) swabs from patients with ILI. Only clinical information and epidemiologic data were collected from Wednesday to Friday.

Specimen collection and transport of specimens
A standard operating procedure was created for the collection and transport of appropriate clinical specimens (NP or OP swabs) from patients with ILI throughout each year and especially during the influenza season or suspected outbreaks. Influenza surveillance officers (ISO) employed by the project or trained personnel like midwives or nurses in the sentinel sites collected the NP/OP swabs. Swabs were transferred to virus transport medium (VTM) and stored in refrigerators before transporting to RITM in cold boxes with ice or frozen gel packs through private couriers or by ISO within a day or two of specimen collection.

Identification of influenza virus
For the ILI sentinel surveillance clinical specimens (NP or OP swabs), we used viral culture to detect influenza virus and other respiratory viruses (Parainfluenza 1,2,3, Respiratory Syncytial Virus (RSV), Herpes Simplex Virus (HSV), adenovirus, enterovirus and rhinovirus) using standard procedures recommended by the WHO Global Influenza Surveillance Network (WHO-GISN) [17]. From viral culture, direct Immunofluorescence (IF) test of the Imagen kit was used to identify influenza virus types A and B. We performed Hemagglutination Inhibition (HAI) test using the WHO Influenza Reagent kit provided by the WHO Collaborating Centre for Reference and Research on Influenza (WHO-CCRI), Melbourne, Australia for further influenza virus strain characterization. For the confirmation of Parainfluenza 1, 2, 3, RSV, HSV and Adenovirus in viral culture, we used Immunofluorescence test while Acid stability testing was used to identify enterovirus and rhinovirus.
The influenza BOD study used real-time reversetranscriptase polymerase chain reaction (RT-PCR) to detect Influenza A from clinical specimens (NP or OP swabs) using the US CDC method or the Applied Biosystems Pandemic H1N1/09 Assay Set v. 2.0 (Life Technologies, CA, USA) [18]. Further subtyping was done to identify seasonal H1 and H3 subtypes using the method of the Centre for Health Protection, Department of Health, Hong Kong SAR [19] and Pandemic H1N1 2009 subtype [18] for Influenza A positive samples. Conventional RT-PCR methodology was used to detect Influenza B and RSV [20].

Data management and statistical analysis
The National Influenza Surveillance database containing information from CRFs that were initially encoded into EPI-INFO software version 6.04 (public domain statistical software from CDC) was used to obtain data on weekly counts of ILI cases, ILI cases with and without specimens, and those that were positive for influenza and other viruses. Analysis was conducted using MS-Excel, STATA 10 (STATA Corporation, College Station, TX, USA), and R 2.15.0 software [21]. We calculated seasonal and alert thresholds, and epidemic curves using both counts (number of influenza positive cases) and proportions (percent positivity rates). We chose to present the results according to percent positivity rates per week.

Determination of influenza season and seasonal threshold
Seasonal threshold was defined as the level of influenza activity that signalled the start and end of the annual influenza season [1]. We defined influenza season as the period in which weekly influenza positive rates were above the average influenza positivity rate for a particular surveillance year for at least three consecutive weeks. Weekly influenza positivity rate was the per cent of ILI with specimens per week, which were positive for influenza virus. We noted the onset, end, and duration of the influenza season for each surveillance year. The onset of influenza season was the week in which positivity rate was above the weekly average positivity rate for that year and which continued for three consecutive weeks. The end of the influenza season was the first week in which the positivity rate was below the average weekly positivity rate for that year and which continued for three consecutive weeks. Once the onset and end of the influenza season was known, weekly data from this period were deleted creating a pre-and post-influenza period. This was done for surveillance years 2006,2007,2008,2010, and 2011 and the average weekly influenza positivity rate for all these years was determined and considered as the seasonal threshold. We excluded 2009 because a pandemic year would not reflect findings for a regular influenza season.

Determination of average epidemic curve and alert threshold
We used the weekly influenza positivity rates for the calculation of average epidemic curve and alert threshold utilizing the 2006, 2007, 2008, 2010 and 2011 data. The average epidemic curve was the usual level of influenza activity in and out of season and alert threshold was the level above which influenza activity was higher than most years [1]. Using Excel, the transmission peaks of 5 years of data were aligned around the median week of peak occurrence during the influenza season for each year studied. Then we calculated the average weekly influenza positivity rate for each week centered on the median peak week of transmission. A 4-week running average was used to smooth the curve that represented the average epidemic curve. The standard deviation (SD) of the mean for each week was calculated to define extreme values, after which a curve was created for those values based on 1.65 SD above the mean representing the upper 90% confidence interval of the mean. This value was used as the alert threshold for severe seasons. After this procedure, weekly data for 2012 was plotted to the curve to determine how 2012 compared with the average epidemic curve, seasonal and alert thresholds calculated from the previous 5 years of observation.

Comparison of circulating influenza viruses with Northern and Southern Hemisphere vaccine antigens
Comparison was done by simply matching the influenza viruses by type, subtype and/or lineage (in the case of influenza B) and by circulating strain for each surveillance year with the Northern (NH) and Southern Hemisphere (SH) vaccine antigens recommended for the influenza seasons for the same surveillance year. Influenza seasonality and duration Influenza virus activity was observed all-year round but with distinct weekly positivity rates above the average weekly influenza positivity rate for a year starting between weeks 24 to 31 and ending during weeks 41 to 52 ( Table 2). Table 2  We observed that within the influenza season described above, there was more than one peak of increased proportion of influenza positive cases (Fig. 2). In 2006, using the 4-week moving average, three peaks were seen. The first peak was higher than the other peaks. Dual peaks were clearly distinct for 2007, 2008 and 2011. 2010 also had a minor second peak, which could only be detected by looking at the bar graph. Starting 2008 up to 2011, we observed minor peaks during the first 11 weeks of the year. However these were just single or two consecutive weeks of increased weekly influenza positive rates.

Number
Seasonal and alert threshold, and epidemic curve Figure 3 shows the non-aligned and aligned transmission peaks of influenza activity from 5 years of surveillance data. Figure 4 shows the plot of percent of cases for 2012 compared to the seasonal threshold (5%), average epidemic curve and alert threshold. The first 2 weeks of 2012 showed weekly influenza positivity rates which were higher than the alert threshold followed by weekly influenza positivity rates that were lower than the seasonal threshold but only up to week 29 when influenza season started, 8 weeks after the expected week of rise of cases in week 21. The intensity was severe going above the average epidemic curve but of short duration. There was a sharp drop of weekly influenza positivity rates after week 30 lower than the alert threshold but still above the seasonal threshold. Thereafter, a minor peak was observed in week 35 conforming to the dual peak pattern observed for the other years. The decrease in the percentage of weekly influenza positivity rates below the seasonal threshold was sustained until week 40 with another small peak occurring in week 46.

Comparison of circulating influenza strains and vaccine strains
Ninety-two percent of SH vaccine strains matched that of the circulating strains in 2006, compared to only 69% for NH strains. This was also observed in 2008 (69% versus 15.12%), and 2010 (66% versus 37%). Percent matched strains were lowest in 2009 (the pandemic year) with both SH and NH at 10.18% (Table 3).

Discussion
In this study we have established the seasonal and alert thresholds, and average epidemic curve based on 5 years of influenza surveillance data using the method described in the WHO surveillance manual [1], described the seasonality of influenza in the Philippines, and presented the influenza strains circulating in the country from 2006 to 2011. While influenza viruses were detected year-round, increased activity was seen from June to November. We used the established thresholds and epidemic curve to assess the influenza season for 2012, its onset and end, and its severity in relation to these thresholds. We also showed that intensity of influenza activity was different for each year, and that different strains of influenza viruses circulated every year. Analysis of influenza activity using weekly detections showed more than one peak of activity within the influenza season from June to November of each year. Different dominant strains were associated with each peak. Two influenza B lineages circulated for five out of six surveillance years. Circulating strains matched the Southern Hemisphere vaccine more often than Northern Hemisphere vaccine strains.
The WHO surveillance manual has emphasized the need for standardized tools to establish thresholds and epidemic curves that could be used to compare influenza season between countries, and to assess the severity of the current season in a country to previous seasons of that same country. The method recommended in the manual [1] is relatively simple and uses software, which is readily available, and can be implemented by NICs as long as five or more years of influenza surveillance data are available. Knowing the usual baseline level of disease and the seasonal pattern as a point of reference aids in determining whether the current season is atypical both in timing and relative severity compared to previous ones. This information can help improve the accuracy of clinical diagnosis, appropriate use of antiviral medication, and the uptake and timeliness of seasonal influenza vaccines. As of this writing, we have not seen studies in Southeast Asia reporting established alert and seasonal thresholds or epidemic curves using the WHO surveillance manual. Although there is one publication from Thailand wherein Early Aberration Reporting System (EARS) and Cumulative sum (CUSUM) were used, calculations were only for ILI [22]. Routine influenza weekly monitoring reports from surveillance activities in Asia included ILI and influenza strains but compared these to historical influenza activity [23]. We have not been able to obtain information from countries in Asia on the use of the WHO method or moving epidemic method (MEM) to establish influenza seasonal thresholds or epidemic curves. At this point in time, there are already many countries in Asia, which have more than 5 years of surveillance data. The WHO method for calculating thresholds and epidemic curves could very well be applied in these countries.
The Philippines, Bangladesh, Cambodia, India, the Lao People's Democratic Republic, Thailand and Vietnam have similar influenza seasonality [5]. Seasonality has been associated not only with climatological but also latitudinal variations [5,[24][25][26]. Influenza virus infection was found throughout the year in these countries, but more than 60% of influenza positivity rates were observed during the months of June to November [5]. The capital cities of these countries lie north of the equator (between 11.6 to 28.7°N). Thus, April to May is the most appropriate month for influenza vaccination for the Philippines and countries with similar seasonality [5]. In areas in India north of 30°N however, there was also increased influenza activity during winter, so that a separate vaccination timing has been recommended [5,24,25] for this country. In contrast, tropical countries in Asia below 11°N like Indonesia (Jakarta, 6.2°N), Malaysia (Kuala Lumpur, 3.2°N), and Singapore (1.3°N), influenza was found all throughout the year with no distinct influenza season [5]. Specific analysis on the association between influenza activity and climate variables has not been done in this study but June to November coincides with the rainy season in the Philippines similar to that of Cambodia [27], Myanmar [28], North Thailand [29], Bangladesh [5,30], Lao [31], and some parts of India South of Srinagar [24,25]. Factors other than the rainy season might also play a role for the increased influenza virus activity starting the month of June in the Philippines. Schools open in June for all levels in the Philippines, suggesting the effect of crowding on influenza transmission.
More than one peak of influenza activity was observed during the surveillance years but these peaks were confined to the defined influenza season from June to November. The presence of the different dominant influenza strains for each of the peaks may be through the following mechanisms: virus interaction leading to interference between viruses, break in transmission due to school vacation, development of immunity to the dominant virus, or a combination of these factors. Goldstein et al. [32] have noted that there was a negative association between strains' incidences indicating that high infection rates with one strain can interfere with the transmission of other strains. Raoult [33] has postulated that the dynamics of influenza virus subtypes against those of other subtypes and even other respiratory viruses is complex and interference between these viruses might impact on their transmissibility in humans. School closure has been associated with reduced transmission of influenza virus because of decreased contact rate [34,35]. School vacation or semestral break occurs at around the end of October to the first week of November. We have observed that influenza activity decreases during this time of school disruption. A break in the transmission of the first peak dominant virus causes co-circulating influenza strains to increase in activity leading to the development of a second peak. Increased transmission of the first dominant influenza strain may later lead to development of sufficient immunity among the population, which in turn will stop or decrease the spread of the virus. Consequently increased activity of co-circulating strains [36,37] occurs. The influenza strains circulating in the Philippines from 2006 to 2011 frequently matched the Southern rather than the Northern Hemisphere vaccine strains. The vaccine currently used in the Philippines is the Southern Hemisphere vaccine. Influenza B strains both of the Yamagata and the Victoria lineages were found during the years of surveillance. With these findings, the appropriate timing of vaccination should be from April to May, the current choice of Southern Hemisphere influenza vaccine should be continued, and most importantly, the vaccine should contain two influenza B lineages. Therefore, quadrivalent rather than trivalent vaccines would be a better choice. There are several limitations to this study. The first limitation is that we used collated data from all sentinel sites to analyze seasonality. Sixty-four percent of influenza positive cases came from Luzon where Manila is located, 15% from the Visayas islands, and 21% from Mindanao. Since majority of the cases came from the Luzon islands, the results may be representative for only the north of the Philippines. The actual latitudinal location of the country extends from 5 to 20°N. Mindanao islands are situated below 9.7°N up to 5°N. The second limitation is that we were not able to compare the WHO method to other analytic methods. The third limitation is that data from this study mostly came from the pediatric age group (<15 years of age). There were very few adult ILI consultations probably because adults are able to handle ILI much better than children. The  fourth limitation is that we could only match vaccine and circulating strains through names of the strains detected and not by genomic sequencing or antigenic testing. Although the Philippine National Surveillance followed the CDC definition for ILI, the definition of ILI was different from other countries and this may be considered a fifth limitation. Despite these limitations, this study has given a valuable insight into the measurement of seasonal and alert thresholds, and epidemic curves that will help DOH influenza program managers in the management of annual influenza epidemics in the country.

Conclusions
Analysis of surveillance data from 2006 to 2012 has provided information about seasonal thresholds, epidemic curves, and circulating strains to guide health personnel in the control of influenza virus infection. Although influenza virus circulated throughout the year, a distinct seasonality was observed from June to November. The ideal time therefore to administer Southern Hemisphere influenza vaccines is from April to May. Quadrivalent vaccine might have more impact on influenza control than trivalent vaccine. Policy-makers are now capable of assessing the intensity or severity of the current influenza epidemic even early in its course, and plan more precisely resources necessary to control the outbreak. Influenza surveillance activities should be continued in the Philippines and funding for such activities should already be incorporated into the Philippine health budget. The results of seasonality patterns will be of immense help in the analysis of future studies on excess mortality during influenza seasons in the Philippines. Continuation of surveillance activities is needed not only to obtain data on circulating strains but also to recalculate the thresholds and epidemic curves as influenza surveillance continues. More sentinel sites should be added in the Mindanao islands. If possible, a separate analysis for seasonality is conducted for the Mindanao area. Further research should be conducted on the association of meteorological and other factors on seasonality of influenza virus infection. Genomic sequencing and antigenic testing of Philippine influenza virus isolates should be conducted to determine more accurately matching of vaccine and circulating strains. Funding Funding for the national influenza surveillance as well as the BOD study conducted by the RITM-Tohoku University collaboration came from the Centers for Disease Control and Prevention, Atlanta, Georgia, USA. The National Influenza Surveillance of the DOH conducted the surveillance according to a protocol submitted to the funding agency. The funding agency (CDC) did not participate in the design of the study and collection, analysis, and interpretation of data except for help by one of the authors in the comparison or matching of the vaccine strains circulating in the country and the strains available in the Southern and Northern Hemisphere vaccine.

Availability of data and materials
The dataset supporting the conclusions of this article comes from the Philippine National Influenza Surveillance dataset (2006 to 2012) which included data from the Influenza BOD study conducted in Baguio City by the RITM-Tohoku University starting November 2008. As such, the authors do not own the dataset and are unable to share the data to the public. However, we are sharing the excel file which we used to calculate the thresholds and epidemic curve (Additional file 1).
Authors' contributions MGL helped in the revisions of the concept and design of the surveillance starting 2007 to 2012 leading to changes in the acquisition of data; wrote the manuscript, analyzed and interpreted the data; gave final approval of the version to be published and agreed to be accountable for all aspects of the work. MTI, AGT, CCD, VLT, AVB, RMO, ESM, VFA, and ANB helped in the revisions of the concept and design of the surveillance starting 2007 to 2012 leading to changes in the acquisition of data; analyzed data; critically appraised the manuscript; and gave final approval of the version to be published and agreed to be accountable for all aspects of the work. LTN helped in the analysis and interpretation of results particularly in the calculation of thresholds and epidemic curves, and analysis of circulating strains; wrote some parts of the manuscript and critically appraised and proof-read the draft manuscript; gave final approval of the version to be published and agreed to be accountable for all aspects of the work. TR helped in the interpretation of the matching of vaccine and circulating strains; critically appraised the manuscript; and gave final approval of the version to be published and agreed to be accountable for all aspects of the work.