Next Article in Journal
Evaluation of Unconventional Supplements to the Diet of Intensively Reared Agouti (Dasyprocta leporina) in Trinidad, West Indies
Previous Article in Journal
Improvement of the Enterotoxigenic Escherichia coli Infection Model for Post-Weaning Diarrhea by Controlling for Bacterial Adhesion, Pig Breed and MUC4 Genotype
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of High-Risk Areas for the Spread of Highly Pathogenic Avian Influenza in Central Luzon, Philippines

by
Roderick Salvador
1,2,
Neil Tanquilut
3,4,
Kannika Na Lampang
5,
Warangkhana Chaisowwong
5,
Dirk Pfeiffer
6 and
Veerasak Punyapornwithaya
5,*
1
College of Veterinary Science and Medicine, Central Luzon State University, Science City of Muñoz 3120, Philippines
2
Faculty of Veterinary Medicine, Chiang Mai University, Chiang Mai 50200, Thailand
3
College of Veterinary Medicine, Pampanga State Agricultural University, Pampanga 2011, Philippines
4
School of Environmental Science and Management, University of the Philippines, Los Banos 4030, Philippines
5
Veterinary Public Health and Food Safety Centre for Asia Pacific (VPHCAP), Faculty of Veterinary Medicine Chiang Mai University, Chiang Mai 50200, Thailand
6
Veterinary Epidemiology, Economics and Public Health Group, Royal Veterinary College, London AL9 7TA, UK
*
Author to whom correspondence should be addressed.
Vet. Sci. 2020, 7(3), 107; https://doi.org/10.3390/vetsci7030107
Submission received: 6 July 2020 / Revised: 2 August 2020 / Accepted: 4 August 2020 / Published: 8 August 2020

Abstract

:
Highly pathogenic avian influenza virus (HPAIV) is a major problem in the poultry industry. It is highly contagious and is associated with a high mortality rate. The Philippines experienced an outbreak of avian influenza (AI) in 2017. As there is always a risk of re-emergence, efforts to manage disease outbreaks should be optimal. Linked to this is the need for an effective surveillance procedure to capture disease outbreaks at their early stage. Risk-based surveillance is the most effective and economical approach to outbreak management. This study evaluated the potential of commercial poultry farms in Central Luzon to transmit HPAI by calculating their respective reproductive ratios ( R 0 ). The reproductive number for each farm is based on the spatial kernel and the infectious period. A risk map has been created based on the calculated R0. There were 882 (76.63%) farms with R 0 < 1. Farms with R 0 ≥ 1 were all located in Pampanga Province. These farms were concentrated in the towns of San Luis (n = 12) and Candaba (n = 257). This study demonstrates the utility of mapping farm-level R 0 estimates for informing HPAI risk management activities.

1. Introduction

Highly pathogenic avian influenza virus (HPAIV) has caused disease outbreaks among domestic poultry and wild bird populations, resulting in significant economic losses in many countries around the world. The potential of the infection to cause illnesses in humans is also a major concern. The Philippines confirmed an outbreak of HPAI H5N6 subtype, based on follow-up tests done in an Australian laboratory in August of 2017. Commercial poultry farms in barangays San Carlos and Santa Rita of San Luis, Pampanga were affected. After a week, cases were confirmed in Nueva Ecija, a province adjacent to Pampanga. A second wave of the outbreak was reported on 12 November 2017 in Cabiao, Nueva Ecija. A total of 23 farms were identified to have bird flu infection in the two provinces from July to November 2017. Sixteen farms were in Pampanga and seven in Nueva Ecija. Eight towns in Pampanga were affected (Apalit, Bacolor, Candaba, Lubao, Mabalacat, San Luis, San Simon, SantoTomas) and four towns in Nueva Ecija were included (Cabiao, Jaen, San Isidro, Zaragosa). During the outbreak, 15 layer chicken farms, 6 duck farms, and 2 quail farms were infected. An outbreak investigation identified key factors affecting the spread of the disease among farms including low biosecurity, having multiple species reared on the farm, no outbreak monitoring system, inability to identify the disease, and uncontrolled trade in poultry commodities. There were 45,682 birds that had died from the disease, while 393,515 birds were killed during the implementation of relevant outbreak management programs. The estimated daily loss as a consequence of the epidemic was about US$3.6 million [1]. Although increased mortality figures among poultry could not be definitively attributed to HPAI, there were accounts that many birds had died from symptoms consistent with HPAI several months prior to the confirmation of the outbreak.
Given the widespread under-reporting of this epidemic that has been associated with inadequate compensation and fear of trade restrictions, passive surveillance measures will be of limited effectiveness. Effective surveillance protocols enable authorities to detect disease outbreaks at an early stage, which should then trigger an immediate control response. As a result, the affected livestock industry will be prevented from incurring severe economic losses. While this may be the ideal situation, in most countries around the world the government’s veterinary service may not be able to cope with the necessary surveillance demands. This is due to the need for the simultaneous monitoring of many different diseases among various species and the need for a single agency to cover a large territory. Consequently, present surveillance strategies mainly rely on reporting by livestock keepers.
Implementation of an effective active surveillance protocol requires increased veterinary manpower, but most municipalities and cities are not prepared to invest in strengthening the capacity of their veterinary services. A more cost-effective alternative to active surveillance is a risk-based surveillance approach. By focusing on the high-risk strata of the population at risk, it is possible to conserve valuable resources. The identification of high-risk strata requires a good understanding of the epidemiological and socio-economic drivers of the spread of the relevant infectious disease. Risk factor studies have been used to identify spatial drivers of infection risk, and these have been used to generate risk maps [2,3,4,5]. Such outputs can be used to identify geographical areas where surveillance efforts should be implemented. Furthermore, several authors have described the successful use of a transmission kernel method for generating risk maps [6,7,8,9,10,11,12,13,14]. In brief, a transmission kernel represents the relative risk of transmission in terms of the distance to an infected farm. It is used as a proxy for all potential routes of infection that are spread between farms [7,12,15,16]. The method allows for the calculation of farm-level basic reproductive ratio values ( R 0 ). The R0 is the expected number of secondary infections that occur among susceptible farms, as well as those that are caused by the first infectious farm throughout the initial infectious period [8,17]. A farm has the potential to cause an epidemic if its R 0 > 1, while it cannot if its R 0 < 1 [14]. The identification of farms that can serve as sources of epidemics can be set as the focus of the surveillance initiatives.
Highly pathogenic avian influenza is included in the priority list of diseases being monitored by the Philippine government. Therefore, this study evaluated the potential of commercial poultry farms in Central Luzon for transmitting HPAI by calculating their basic reproductive ratio numbers. This information was then used to produce a risk map based on the calculated R0.

2. Materials and Methods

2.1. Study Area

The study was conducted in Central Luzon, Philippines. The Central Luzon region is composed of 7 provinces including Aurora, Bataan, Bulacan, Nueva Ecija, Pampanga, Tarlac, and Zambales (Figure 1). According to the data of the Philippine Statistics Authority, Central Luzon is the region with the highest density of ducks (32%) as of January 1, 2019 [18]. Moreover, the region is home to 68.8% of the total ducks kept in commercial farms. Likewise, the region is also home to the highest number of chickens with a 35.9% share of the country’s total chicken population over the same period of time [19]. A commercial farm is defined as a poultry operation with a minimum of 500 layers or 1000 broilers. For combined operation, a farm with at least 100 layers and 100 broilers is considered commercial.

2.2. Data Sources

HPAI outbreak data were provided by the Bureau of Animal Industry (BAI), Philippines. The poultry farms included in this study were based on the registration records for commercial poultry operations of the Provincial Veterinary Offices (PVO). The farm profiles (poultry population by species, farm coordinates) were compiled from farm visits that had been conducted throughout the region. Each farm’s maximum poultry capacity, and not the actual bird population at the time of the visits, were used in the analysis.

2.3. Model Formulation

The calculation of the farm level reproductive number ( R 0 ) as described by Boender et al. [6,11] requires the computation of the transmission kernel (h) and the stochastic infectious period (T). The transmission kernel (h) was used to represent all forms of disease transmission between susceptible farms [15]. This parameter has been used to HPAI [7,20,21,22] and other infectious animal disease models like foot and mouth disease [15,23,24]. It refers to the infection hazard posed by an infectious animal from an infected farm to a susceptible animal of a susceptible farm. It is a function of the Euclidean distance (r) between farms. We adopted the kernel formulation (Equation (1)) from the studies of Boender et al. [6,11] and Dorigatti et al. [22]:
h ( r i j ) = h 0 ( 1 + ( r i j / r 0 ) α )
The r i j refers to the Euclidean distance between an infectious farm j and a susceptible farm i. The h 0 is the maximum hazard rate, occurring when r 0 = 0. The r 0 influences how far the hazard rate extends over distance. The α parameter influences the rate of decay in the hazard rate from the maximum. A dataset for the pairwise Euclidean distances (r) between farms was created. The distance between farms was calculated using the distGeo function of the geosphere in R statistical software version 3.5.3 (R Core Team, Vienna, Austria) [25].
To calculate the likelihood function, the force of infection on a susceptible farm i at time t was given by Equation (2) as described in [6]:
λ i ( t ) = j i h ( r i j )   1   [   j   is   infectious   ]
With λ i ( t ) as the force of infection or the cumulative hazard rate experienced by farm i on day t, the probability that farm i is infected on day t is:
q i ( t ) = 1 e λ i ( t )
The probability that farm i remains uninfected up to day t is:
r i ( t ) = e s = 1 t 1 λ i ( s )
Given the risk of infection above, the log-likelihood function is given by:
= k ϵ K t = 1 t m a x 1 λ k ( t ) l ϵ Λ t = 1 t c u l , l 1 λ l ( t ) m ϵ M t = 1 t i n f , m 1 λ m ( t ) + m ϵ M l o g ( 1 e λ m ( t i n f , m ) )
where the set K contains all farms that remained uninfected and that were not culled; Λ contains the farms that were not infected but that were culled (at times t c u l ,   l ); M contains the farms that were infected (at times t i n f ,   m ). With Equation (5), the estimates for the transmission kernel parameter values ( h 0 , r 0 , and α) were calculated through maximum likelihood estimate (MLE) as described in [6] using the 2017 Philippines HPAI epidemic data. The bbmle functions of the R program [25] was used for this procedure.
Following the procedure described by Boender et al. [6], the basic reproduction number was estimated for each poultry farm in Central Luzon using the spatial transmission kernel ((h(rij)) and the stochastic infectious period (T) of farms i. The reproduction number of farms i, R 0 , is given by:
R 0 = j i ( 1 E   [ e h ( r i j ) T i ]
The computed R0 after 1000 iterations was used as the reference in developing the risk map. The doParallel functions of R [25] was used to calculate the reproductive ratio of each farm. Farms with at least 1.0 R 0 were classified as high-risk farms. The risk map was created using QGIS (Open Source Geospatial Foundation Project, Zurich, Switzerland) [26].

3. Assumptions

Upon infection, a farm passes through a latency period of 1–2 days [22,27,28,29]. In this stage, the avian influenza virus develops within the host until the host becomes infectious by shedding the virus. Viral inoculation experiments of HPAIV showed that birds can start shedding the virus after 0.76 day [30], or 1.1 day [31]. Though it is possible that birds can start to become infectious in less than 24 h after infection, we adopted the results of Vander Goot [32]. Their results showed that a latent period of 1–2 days fits the data better than the model with a latent period of less than one day. This is likewise adopted by Boender et al. [6], and Dorigatti et al. [22] in their HPAI simulation studies. The detection of HPAI in the farm is assumed to be on days 5–7. This is based on mortality thresholds [33]. Several studies have shown mortality thresholds to be the most reliable indicator of AI infection in poultry farms [34,35]. The assumptions used here are not specific for the H5N6 subtype. However, the range of transmission rate of the HPAI H5N6 virus in domestic poultry is close to or within the range of transmission rates found for other HPAI subtypes. This is generally applicable in cases where transmission experiments are limited. The infection experiment by Jiao (2016) showed that inoculated chicken died at 4–7 days post inoculation (dpi) [36]. Chicken succumbed by 4–5 dpi based on the results of Kim (2018) [37]. We assumed that suspicion of HPAIV infection due to increased mortalities will happen between 5 and 7 days after infection. Based on the response of farm managers to a survey, they will be notifying authorities about the presence of the disease 5–7 days post detection. Using the outbreak data as reference, culling of an infected farm is initiated 1–3 days after receipt of the report. The days required for completion of culling of infected farms after reporting were 1–3 days, depending on the size of the farm. Specifically, for farms with a population ≤50,000, the culling needs 1 day; for farms with 50,000 < population ≤100,000, the culling needs 2 days; for farms with population >100,000, the culling needs 3 days. The infectious period lasts until the completion of culling.

4. Results

4.1. Farm Location

The locations of the poultry farms included in this study are shown in Figure 2. There was a total of 1151 poultry farms that kept chickens (n = 880), ducks (n = 260), or quail (n = 11). The duck farms identified in the study were all located in Pampanga Province. The poultry farm distribution data in the study area, and their average flock size are presented in Table 1.
Poultry production levels are considered high in the provinces of Nueva Ecija, Pampanga, and Bulacan. These provinces represent 78.8% of the total poultry farms in the region. Likewise, they share 70.3% of the total poultry population. Large scale operations were observed in the provinces of Bulacan, Tarlac, Nueva Ecija, and Zambales considering their high average flock sizes. In contrast, Pampanga was able to significantly contribute to the regional poultry production levels because of the high number of poultry farms (n = 485) despite its lower average flock size. Aurora had the smallest number of farms. This may be due to the logistical challenges of poultry farming that exist in the area. This would likely be a consequence of its remoteness and it being in a mountainous area.

4.2. Farm Practices

At this study site, ducks were raised mainly for the purposes of egg production. Incorporated into the management protocol of duck farms is the concept of foraging. This provides ducks the opportunity to consume snails which are believed to improve the shell quality of the eggs produced by these ducks. Foraging sites include swamps and empty rice paddies. This practice requires frequent movement of the flock in search of foraging areas. Ducks were moved via trucks with customized battery cages. Foraging in a particular geographical area may range from 3 to 7 days before the ducks are moved to a different area. Fresh and processed eggs were then sold to traders or directly to consumers through bakeries and restaurants.
In contrast, chickens are always kept within the boundaries of the farm. Chicken raising operations can be dedicated either to layer or broiler chickens. At the end of a production cycle, poultry manure can be sold to vegetable farms as fertilizer. Broilers were either bought by traders or by the company partner if the operation was under contract. For the layers, the spent hens were sold to poultry traders.

4.3. Transmission Parameters

Maximum likelihood estimation based on the 23 infected farms during the 2017 HPAI outbreak resulted in the following parameter estimates: h 0 = 0.0012 day−1 (CI: 0.0001–0.1), r 0 = 3.4 km (CI: 1.001–10.0), and α = 1.4 (CI: 1.001–5.0), with = 206.99 and AIC = 419.98.

4.4. Farm-Level Basic Reproductive Ratio Values

The summary statistics for the farm-level farm reproductive ratio are presented in Table 2. The median of the farm-level basic reproductive ratio values for the 1151 commercial farms in Central Luzon was 0.424. There were 882 (76.63%) farms with a R0 of less than 1. Chickens are the predominant species (96.2%, 849/882) raised in these farms. Meanwhile, ducks are the major species (86.98%, 234/269) in farms with a R 0 ≥ 1.
Notably, around 77% (882/1151) of the total poultry farms in Central Luzon are outside the high-risk areas. These farms include 26 duck farms, 7 quail farms, and 849 chicken farms (Figure 3). All farms with a R 0 ≥ 1 are located in Pampanga Province. These farms are concentrated in the towns of San Luis (n = 12) and Candaba (n = 257). This area is composed of 234 duck farms, 31 chicken farms, and 4 quail farms.

5. Discussion

In this study, farm-level basic reproductive ratio values were calculated and used to generate a risk map for the transmission of HPAI in Central Luzon. The transmission kernel was derived from the profiles of the 23 farms infected during the 2017 HPAI epidemic in the Philippines. These farms were identified through reports from farm owners to the BAI. It is acknowledged that the results from this modeling exercise need to be used cautiously, since only this single outbreak could be used for generating the transmission kernel and there may have been underreporting of farm locations and disease events. The findings are discussed below in the context of the wider literature to allow an assessment of their epidemiological validity.
To construct the model, the factors used to account for the species-specific transmissibility were developed by Hayama [8], Keeling [10], and Dorigatti [22]. The number of animals in the farm is likewise added to the equation used by Keeling [10] and Hayama [23]. The inclusion of a species-specific transmission coefficient is considered important if there are variations in susceptibility and the occurrence of infections among different species, such as is the case with avian influenza. Unfortunately, no further parameters could be derived from this small dataset. The approach used here has also been applied in order to generate risk maps for FMD [11] and HPAI [6] in the Netherlands.
Moreover, with the available dataset, kernel density estimation may be the best approach. Other spatial techniques include the cluster detection method, the relative space method, and risk factor identification [38]. The research studies of Ahmed et al. [39], Loth et al. [2], and Tiensin et al. [40] likewise were aimed at identifying high-risk areas for the spread of HPAI. Both Ahmed and Loth used data acquired from the Bangladesh outbreaks, while Tiensin et al. [40] used data obtained from the Thailand outbreaks. Bangladesh had two epidemic waves of HPAI in two consecutive years, while Thailand generated outbreak data from 2004 to 2007. In these cases, the availability of large and more detailed datasets allowed for the application of the spatial cluster detection method and risk factor identification.
The calculated R 0 can be used to inform risk-based surveillance. Currently, the Philippines does not have laws providing veterinary authorities with the legal powers to investigate the health status of animal farms in the absence of a disease outbreak. But if high risk farms can be identified based on scientific criteria, such as the R0, it might be possible to convince farmers of the need for risk-based surveillance protocols. Large poultry industry companies could use the spatial patterns of R0 to identify low-risk geographical areas to target in their production expansion initiatives.
The results indicate that duck farms represent the majority (87%) of farms with an R 0 ≥ 1 (Table 2). Wild and domestic ducks have been associated with the spread and maintenance of avian influenza viruses (AIV) [40,41,42]. Free grazing ducks played a significant role in the maintenance and transmission of infections during the second wave of HPAI epidemic of Thailand in 2004 [43] The AIV of the H5N6, H5N2, and H9N3 subtypes were isolated from healthy ducks for sale in live markets in Hanoi, Vietnam [44]. The study of Sturm-Ramirez et al. (2005) showed that the HPAI H5N1 virus causes minimal signs of diseases in ducks. Consequently, they can spread HPAI silently and efficiently among domestic and wild birds in Asia [45].
Part of the production management procedure for duck farmers involves grazing in post-harvest rice paddy fields, and thereby, requires frequent movement between habitat sites which the ducks share with other wild birds. Analysis of the HPAI-H5N1 epidemic in Thailand revealed that outbreaks are commonly reported in the central and northern regions where duck population densities are high [46]. Similarly, low pathogenic avian influenza (LPAI) subtypes H4N6 and H3N8 were isolated from sentinel flocks of ducks in Thailand several years after the epidemic [47]. The 2017 HPAI H5N6 outbreak in the Philippines confirmed that HPAI is virulent among domestic ducks. This was also observed among H5N1 infected ducks in 2002. The virus killed many aquatic bird species. Earlier, ducks infected with highly pathogenic H5 or H7 avian influenza virus are consistently asymptomatic or had very mild disease. The outbreaks in Thailand, Vietnam, and Japan from 2002 to 2004 were pathogenic to ducks. Hulse-Post et al. [48] reported that there is a possibility for AIV to revert to its non-pathogenic form among the duck population. Their study showed a trend toward decreased pathogenicity among ducks. The virus remained highly pathogenic to chickens and humans. This finding demonstrates the limitation of using only clinical signs to determine infection with HPAI among the duck population. Consequently, the reliance on passive surveillance will not be effective. Thus, regular monitoring through testing among the duck population should be included in the strategy against HPAIV.
The vast majority of farms with R 0 ≥ 1 were located in Candaba, Pampanga (95.5%, n = 257). The Candaba swamp that is in this town is a recognized stopover for migratory birds from the Northern Hemisphere, especially for birds from China and Russia. During the dry season, part of the swamp is utilized for rice production [49]. Epidemiological studies of HPAI outbreaks have identified the relevance of the proximity to natural wetlands and water bodies as an important risk factor [50,51,52,53]. These areas provide an interface for contact between domestic and wild ducks.
The identification of the risk factors associated with the spread of HPAI is outside the objectives of this study. The spatial risks presented above are brought out as the observed features of the farms that are at high risk for outbreaks (high farm density, high duck farm density, and proximity to the Candaba swamp). The risks that were identified agree with the findings of prior studies. Programs to address these spatial risks should be integrated into the plans being put into place to address the avian influenza concern. The potential indirect contact between domestic ducks and infected migratory birds through sharing the same foraging areas cannot be avoided. Feeding in post-harvest rice paddy fields is integral to the production system of duck farms. Surveillance efforts should be comprehensive among the duck population, and both low and highly pathogenic AIV should be covered. To supplement active surveillance practices, sentinel flocks can be set up to monitor the presence of AIV.
The risk of spread is more likely to be observed in areas with high levels of farm density. The affected area in the 2010 FMD epidemic in Japan was the Miyazaki Prefecture. This prefecture is home to the highest density of cattle farms (165.0 farms/100 km2) and the second highest density of pig farms (9.8 farms/100 km2). Proximity (≤1.5 km) to an infected farm has been associated with the spread of HPAI in Italy [53]. The HPAI outbreaks in Thailand occurred in areas with high levels of poultry farm-density [46].
In part, the density of farms can be regulated through the development of policies that limit the number of farms allowed in an area. This can be implemented in current low-risk areas. Regulations, however, can only cover commercial farms as these are required to obtain government permits for their operation. The establishment of small-scale operations cannot be controlled unless a mandatory registration of farms with at least 50 birds is implemented. To address the risk of outbreaks spreading into high-density areas, farms must be required to maintain high levels of biosecurity. There should be information campaigns on improved biosecurity practices among poultry raisers in both high- and low-risk areas. Another aspect of the programs to be instituted would include the identification and description of the trading dynamics that exist between the farms located in high-risk areas and their customers. This information will facilitate contact-tracing in the event of an outbreak. With the risk map as a guide, management strategies in the event of an HPAI outbreak are dependent upon the location of the index farm. Quarantine procedures and the culling of infected farms located outside high-risk areas would be sufficient in controlling the potential spread. In contrast, farms with a high probability of disease transmission require additional pre-emptive culling to mitigate the possibility of an outbreak. One limitation of this study is the unknown number of unregistered quail and duck farms in contrast to poultry farms, which by law have been officially registered. While 261 duck farms were registered, we believe that there were other commercial duck farms outside Pampanga. Poultry farms were identified based on the records of Provincial Veterinary Offices. Apparently, a number of small-scale duck and quail producers were not compliant in registering their operations. Compared to the operation of other species, a commercial quail farm does not require much space, and the transportation of quail eggs to the market can be done in concealed ways such as in vans or in the boots of cars. Quails are known to be highly susceptible to HPAI, and the clinical disease period can last longer than in chickens. Thus, quail farms are known to have high potential for transmitting the disease [54,55]. It is therefore likely that the current map is biased due to missing observations that are associated with specific farm species that are known to be important risk factors for avian influenza.
Based on the findings from this study, it is recommended to expand the coverage of the risk map to a national scale. Thus, mandatory registration of all poultry holdings will be required. This will then potentially allow to identify further high-risk areas where more intensive surveillance should be implemented. Furthermore, a national database of poultry farm profiles would allow epidemic modeling projects related to other avian diseases.

6. Conclusions

In conclusion, the HPAI risk map was generated by using the transmission kernel estimated from the 2017 epidemic data. This method allowed us to identify high-risk areas by calculating the farm-level basic reproductive ratio. The resulting risk map facilitates visual communication of the findings to decision-makers, and it thereby serves as an effective tool in support of the development of HPAI control and prevention strategies for specific geographic areas. Importantly, risk maps are dynamic, and need to be critically reviewed and updated once new information becomes available. Extending the risk mapping process to the national level with inclusion of small-scale poultry farms (of at least 50 birds) may allow identification of other high-risk areas in the country.

Author Contributions

R.S., V.P., K.N.L., and W.C. contributed to conception and design. N.T. assisted in the collection of data. R.S. performed the computer simulations and drafted the article. V.P. supervised and reviewed the analysis. V.P. and D.P. improved the material for publication. All authors have read and agreed to the published version of the manuscript.

Funding

The publication of this article was funded by Veterinary Public Health and Food Safety Centre for Asia Pacific (VPHCAP), Faculty of Veterinary Medicine, Chiang Mai University. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

The authors would like to thank Ronnie Domingo, Janice Sabagay-Garcia, Anthony Bucad of the Bureau of Animal Industry, the managers of the participating farms and everyone who extended their support in realizing this study.

Conflicts of Interest

The authors declare no conflict of interest.

Data Availability

The data used in this study is available at figshare. DOI: https://doi.org/10.6084/m9.figshare.11987604.v1. The geographical coordinates of farms used in the analysis were slightly modified due to privacy issue.

References

  1. Lee, H.; Lao, A. Transmission dynamics and control strategies assessment of avian influenza A (H5N6) in the Philippines. Infect. Dis. Model. 2018, 3, 35–59. [Google Scholar]
  2. Loth, L.; Gilbert, M.; Osmani, M.G.; Kalam, A.M.; Xiao, X. Risk factors and clusters of highly pathogenic avian influenza H5N1 outbreaks in Bangladesh. Prev. Vet. Med. 2010, 96, 104–113. [Google Scholar] [CrossRef] [Green Version]
  3. Loth, L.; Gilbert, M.; Wu, J.; Czarnecki, C.; Hidayat, M.; Xiao, X. Identifying risk factors of highly pathogenic avian influenza (H5N1 subtype) in Indonesia. Prev. Vet. Med. 2011, 102, 50–58. [Google Scholar] [CrossRef] [Green Version]
  4. Pfeiffer, D.U.; Minh, P.Q.; Martin, V.; Epprecht, M.; Otte, M.J. An analysis of the spatial and temporal patterns of highly pathogenic avian influenza occurrence in Vietnam using national surveillance data. Vet. J. 2007, 174, 302–309. [Google Scholar] [CrossRef] [PubMed]
  5. Dhingra, M.S.; Dissanayake, R.; Negi, A.B.; Oberoi, M.; Castellan, D.; Thrusfield, M.; Linard, C.; Gilbert, M. Spatio-temporal epidemiology of highly pathogenic avian influenza (subtype H5N1) in poultry in eastern India. Spat. Spatiotemporal. Epidemiol. 2014, 11, 45–57. [Google Scholar] [CrossRef] [PubMed]
  6. Boender, G.J.; Hagenaars, T.J.; Bouma, A.; Nodelijk, G.; Elbers, A.R.W.; De Jong, M.C.M.; Van Boven, M. Risk maps for the spread of highly pathogenic avian influenza in poultry. PLoS Comput. Biol. 2007, 3, e71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Boender, G.J.; Elbers, A.R.W.; De Jong, M.C.M. Spread of avian influenza in The Netherlands: Identifying areas at high risk. Vet. Ital. 2007, 43, 605–609. [Google Scholar] [PubMed]
  8. Hayama, Y.; Yamamoto, T.; Kobayashi, S.; Muroga, N.; Tsutsui, T. Evaluation of the transmission risk of foot-and-mouth disease in Japan. Virology 2015, 9, 1167–1170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. McCullagh, M. Detecting hotspots in time and space. Int. Symp. Exhib. Geoinf. 2006, 349, 1–18. [Google Scholar]
  10. Keeling, M.J.; Woolhouse, E.J.; Shaw, D.J.; Matthews, L.; Chase-Topping, M.; Haydon, D.T.; Cornell, S.J.; Kappey, J.; Wilesmith, J.; Grenfell, B.T. Dynamics of the 2001 UK foot and mouth epidemic: Stochastic dispersal in a heterogeneous landscape. Science 2001, 294, 813–817. [Google Scholar] [CrossRef] [Green Version]
  11. Boender, G.J.; Van Roermund, H.J.; De Jong, M.C.M.; Hagenaars, T.J. Transmission risks and control of foot-and-mouth disease in The Netherlands: Spatial patterns. Epidemics 2010, 2, 36–47. [Google Scholar] [CrossRef] [PubMed]
  12. Tildesley, M.J.; Keeling, M.J. Is R 0 a good predictor of final epidemic size: Foot-and-mouth disease in the UK. J. Theor. Biol. 2009, 258, 623–629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. LeMenach, A.; Legrand, J.; Grais, R.F.; Viboud, C.; Valleron, A.-J.; Flahault, A. Modeling spatial and temporal transmission of foot-and-mouth disease in France: Identification of high-risk areas. Vet. Res. 2005, 36, 699–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hartemink, N.A.; Purse, B.V.; Meiswinkel, R.; Brown, H.E.; De Koeijer, A.; Elbers, A.R.W.; Boender, G.J.; Rogers, D.J.; Heesterbeek, J.A.P. Mapping the basic reproduction number (R0) for vector-borne diseases: A case study on bluetongue virus. Epidemics 2009, 1, 153–161. [Google Scholar] [CrossRef]
  15. Ster, I.C.; Ferguson, N.M. Transmission parameters of the 2001 foot and mouth epidemic in Great Britain. PLoS ONE 2007, 2, e502. [Google Scholar]
  16. Boender, G.J.; Meester, R.; Gies, E.; De Jong, M.C.M. The local threshold for geographical spread of infectious diseases between farms. Prev. Vet. Med. 2007, 82, 90–101. [Google Scholar] [CrossRef]
  17. Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A.J. On the definition and the computation of the basic reproduction ratio R o in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef] [Green Version]
  18. Bersales, L.G.S. Duck situation report January–December 2018. Philipp. Stat. Auth. 2019, 1–37, ISSN 2546-0625. [Google Scholar]
  19. Bersales, L.G.S. Chicken situation report January–December 2018. Philipp. Stat. Auth. 2019, 1–45, ISSN 2546-0625. [Google Scholar]
  20. Truscott, J.; Garske, T.; Chis-ster, I.; Guitian, J.; Pfeiffer, D.; Snow, L.; Wilesmith, J.; Ferguson, N.M.; Ghani, A.C. Control of a highly pathogenic H5N1 avian influenza outbreak in the GB poultry flock. Proc. R. Soc. B Biol. Sci. 2007, 1623, 2287–2295. [Google Scholar] [CrossRef] [Green Version]
  21. Bonney, P.J.; Malladi, S.; Boender, G.J.; Weaver, J.T.; Ssematimba, A.; Halvorson, D.A.; Cardona, C.J. Spatial transmission of H5N2 highly pathogenic avian influenza between Minnesota poultry premises during the 2015 outbreak. PLoS ONE 2018, 13, e0204262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Dorigatti, I.; Mulatti, P.; Rosà, R.; Pugliese, A.; Busani, L. Modelling the spatial spread of H7N1 avian in fl uenza virus among poultry farms in Italy. Epidemics 2010, 2, 29–35. [Google Scholar] [CrossRef] [PubMed]
  23. Hayama, Y.; Yamamoto, T.; Kobayashi, S.; Muroga, N.; Tsutsui, T. Mathematical model of the 2010 foot-and-mouth disease epidemic in Japan and evaluation of control measures. Prev. Vet. Med. 2013, 112, 183–193. [Google Scholar] [CrossRef] [PubMed]
  24. Keeling, M.J. Models of foot-and-mouth disease. Proc. R. Soc. B Biol. Sci. 2005, 272, 1195–1202. [Google Scholar] [CrossRef] [Green Version]
  25. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2008; ISBN 3-900051-07-0. [Google Scholar]
  26. QGIS Development Team QGIS Geographic Information System. Open Source Geospatial Foundation Project 2019. Available online: http://www.qgis.org (accessed on 5 January 2020).
  27. Swayne, D.E.; Beck, J.R. Experimental study to determine if low-pathogenicity and high-pathogenicity avian influenza viruses can be present in chicken breast and thigh meat following intranasal virus inoculation. Avian Dis. 2005, 49, 81–85. [Google Scholar] [CrossRef]
  28. Tian, G.; Zhang, S.; Li, Y.; Bu, Z.; Liu, P.; Zhou, J.; Li, C.; Shi, J.; Yu, K.; Chen, H. Protective efficacy in chickens, geese and ducks of an H5N1-inactivated vaccine developed by reverse genetics. Virology 2005, 341, 153–162. [Google Scholar] [CrossRef] [Green Version]
  29. Van Der Goot, J.A.; Koch, G.; De Jong, M.C.M.; Van Boven, M. Quantification of the effect of vaccination on transmission of avian influenza (H7N7) in chickens. Proc. Natl. Acad. Sci. USA 2005, 102, 18141–18146. [Google Scholar] [CrossRef] [Green Version]
  30. Hénaux, V.; Samuel, M.D.; Bunck, C.M. Model-based evaluation of highly and low pathogenic avian influenza dynamics in wild birds. PLoS ONE 2010, 5, e10997. [Google Scholar] [CrossRef]
  31. Spekreijse, D.; Bouma, A.; Koch, G.; Stegeman, J.A. Airborne transmission of a highly pathogenic avian influenza virus strain H5N1 between groups of chickens quantified in an experimental setting. Vet. Microbiol. 2011, 152, 88–95. [Google Scholar] [CrossRef] [Green Version]
  32. Van der Goot, J.A.; De Jong, M.C.M.; Koch, G.; Boven, V.M. Comparison of the transmission characteristics of low and high pathogenicity avian influenza A virus (H5N2). Epidemiol. Infect. 2003, 131, 1003–1013. [Google Scholar] [CrossRef]
  33. Dorea, F.C.; Vieira, A.R.; Hofacre, C.; Waldrip, D.; Cole, D.J. Stochastic model of the potential spread of highly pathogenic avian influenza from an infected commercial broiler operation in Georgia. Avian Dis. 2010, 5, e172–e173. [Google Scholar] [CrossRef]
  34. Elbers, A.A.R.W.; Fabri, T.H.F.; De Vries, T.S.; De Wit, J.J.; Pijpers, A.; Koch, G. The highly pathogenic Avian influenza A (H7N7) virus epidemic in The Netherlands in 2003—Lessons learned from the first five outbreaks. Avian Dis. 2004, 48, 691–705. [Google Scholar] [CrossRef] [PubMed]
  35. Savill, N.J.; Rose, S.G.S.; Woolhouse, M.E.J.; Interface, J.R.S.; Savill, N.J.; Rose, S.G.S.; Woolhouse, M.E.J. Detection of mortality clusters associated with highly pathogenic avian influenza in poultry: A theoretical analysis. J. R. Soc. Interf. 2008, 5, 1409–1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Jiao, P.; Cui, J.; Song, Y.; Song, H.; Zhao, Z.; Wu, S.; Qu, N.; Wang, N.; Ouyang, G.; Liao, M. New reassortant H5N6 highly pathogenic avian influenza viruses in Southern China, 2014. Front. Microbiol. 2016, 7, 1–12. [Google Scholar] [CrossRef]
  37. Lee, E.K.; Lee, Y.N.; Kye, S.J.; Lewis, N.S.; Brown, I.H.; Sagong, M.; Heo, G.B.; Kang, Y.M.; Cho, H.K.; Kang, H.M.; et al. Characterization of a novel reassortant H5N6 highly pathogenic avian influenza virus clade 2.3.4.4 in Korea, 2017 correspondence. Emerg. Microbes Infect. 2018, 7, 1–8. [Google Scholar]
  38. Moore, D.A.; Carpenter, T.E. Spatial analytical methods and GIS: Use in health research and epidemiology. Epidemiol. Rev. 1999, 21, 143–161. [Google Scholar] [CrossRef] [Green Version]
  39. Ahmed, S.S.U.; Ersbøll, A.K.; Biswas, P.K.; Christensen, J.P. The space-time clustering of highly pathogenic avian influenza (HPAI) H5N1 outbreaks in Bangladesh. Epidemiol. Infect. 2010, 138, 843–852. [Google Scholar] [CrossRef] [Green Version]
  40. Tiensin, T.; Ahmed, S.S.U.; Rojanasthien, S.; Songserm, T.; Ratanakorn, P.; Chaichoun, K.; Kalpravidh, W.; Wongkasemjit, S.; Patchimasiri, T.; Chanachai, K.; et al. Ecologic risk factor investigation of clusters of Avian influenza A (H5N1) virus infection in Thailand. J. Infect. Dis. 2009, 199, 1735–1743. [Google Scholar] [CrossRef] [Green Version]
  41. Gilbert, M.; Chaitaweesub, P.; Parakamawongsa, T.; Premashthira, S.; Tiensin, T.; Kalpravidh, W.; Wagner, H.; Slingenbergh, J. Free-grazing ducks and highly pathogenic Avian Influenza, Thailand. Emerg. Infect. Dis. J. 2006, 12, 227–234. [Google Scholar] [CrossRef] [Green Version]
  42. Wibawa, H.; Bingham, J.; Nuradji, H.; Lowther, S.; Payne, J.; Harper, J.; Junaidi, A.; Middleton, D.; Meers, J. Experimentally infected domestic ducks show efficient transmission of Indonesian H5N1 highly pathogenic Avian influenza virus, but lack persistent viral shedding. PLoS ONE 2014, 9, e83417. [Google Scholar] [CrossRef]
  43. Tiensin, T.; Nielen, M.; Vernooij, H.; Songserm, T.; Kalpravidh, W.; Chotiprasatintara, S.; Chaisingh, A.; Wongkasemjit, S.; Chanachai, K.; Thanapongtham, W.; et al. Transmission of the highly pathogenic Avian influenza virus H5N1 within flocks during the 2004 epidemic in Thailand. J. Infect. Dis. 2007, 196, 1679–1684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Nguyen, D.C.; Uyeki, T.M.; Jadhao, S.; Maines, T.; Shaw, M.; Matsuoka, Y.; Smith, C.; Rowe, T.; Lu, X.; Hall, H.; et al. Isolation and characterization of Avian influenza viruses, including highly pathogenic H5N1, from poultry in live bird markets in Hanoi, Vietnam, in 2001. J. Virol. 2005, 79, 4201–4212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sturm-Ramirez, K.M.; Hulse-Post, D.J.; Govorkova, E.A.; Humberd, J.; Seiler, P.; Puthavathana, P.; Buranathai, C.; Nguyen, T.D.; Chaisingh, A.; Long, H.T.; et al. Are ducks contributing to the endemicity of highly pathogenic H5N1 influenza virus in Asia? J. Virol. 2005, 79, 11269–11279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Thanawat, T.; Prasit, C.; Thaweesak, S.; Arunee, C.; Wirongrong, H.; Chantanee, B.; Tippawon, P.; Sith, P.; Alongkorn, A.; Marius, G.; et al. Highly pathogenic Avian influenza H5N1, Thailand, 2004. Emerg. Infect. Dis. J. 2005, 11, 1664. [Google Scholar]
  47. Boonyapisitsopa, S.; Chaiyawong, S.; Nonthabenjawan, N. Sentinel model for influenza A virus monitoring in free-grazing ducks in Thailand. Vet. Microbiol. 2016, 182, 35–43. [Google Scholar] [CrossRef] [PubMed]
  48. Hulse-Post, D.J.; Sturm-Ramirez, K.M.; Humberd, J.; Seiler, P.; Govorkova, E.A.; Krauss, S.; Scholtissek, C.; Puthavathana, P.; Buranathai, C.; Nguyen, T.D.; et al. Role of domestic ducks in the propagation and biological evolution of highly pathogenic H5N1 influenza viruses in Asia. Proc. Natl. Acad. Sci. USA 2005, 102, 10682–10687. [Google Scholar] [CrossRef] [Green Version]
  49. Melendres, R.G. The utilization of Candaba swamp from prehistoric to present time: Evidences from archaeology, history and ethnography. Bhatter Coll. J. Multidiscip. Stud. 2014, 4, 81–93. [Google Scholar]
  50. Biswas, P.K.; Christensen, J.P.; Ahmed, S.S.U.; Barua, H.; Das, A.; Rahman, M.H.; Giasuddin, M.; Hannan, A.S.M.A.; Habib, A.M.; Debnath, N.C. Risk factors for infection with highly pathogenic influenza A virus (H5N1) in commercial chickens in Bangladesh. Vet. Rec. 2009, 164, 743–746. [Google Scholar] [CrossRef]
  51. Ward, M.P.; Maftei, D.; Apostu, C.; Suru, A. Environmental and anthropogenic risk factors for highly pathogenic avian influenza subtype H5N1 outbreaks in Romania, 2005—2006. Vet. Res. Commun. 2008, 32, 627–634. [Google Scholar] [CrossRef]
  52. Chunxiang, C.A.O.; Min, X.U.; Chaoyi, C.; Yong, X.U.E.; Shaobo, Z.; Liqun, F.; Wuchun, C.A.O.; Hao, Z.; Mengxu, G.A.O.; Qisheng, H.E.; et al. Risk analysis for the highly pathogenic avian influenza in Mainland China using meta-modeling. Chin. Sci. Bull. 2010, 55, 4168–4178. [Google Scholar]
  53. Fang, L.; Vlas, S.J.D.; Liang, S.; Looman, C.W.N.; Gong, P.; Xu, B.; Yang, H.; Richardus, J.H.; Cao, W. Environmental factors contributing to the spread of H5N1 Avian influenza in Mainland China. PLoS ONE 2008, 3, e2268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Perez, D.R.; Lim, W.; Seiler, J.P.; Yi, G.; Peiris, M.; Shortridge, K.F.; Webster, R.G. Role of quail in the interspecies transmission of H9 influenza A viruses: Molecular changes on HA that correspond to adaptation from ducks to chickens. J. Virol. 2003, 77, 3148–3156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Webster, R.G.; Guan, Y.; Peiris, M.; Walker, D.; Krauss, S.; Zhou, N.N.; Govorkova, E.A.; Ellis, T.M.; Dyrting, K.C.; Sit, T.; et al. Characterization of H5N1 influenza viruses that continue to circulate in geese in Southeastern China. J. Virol. 2002, 76, 118–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Map of the Philippines showing the study provinces.
Figure 1. Map of the Philippines showing the study provinces.
Vetsci 07 00107 g001
Figure 2. Map of Central Luzon, Philippines showing the physical location of all 1151 commercial poultry farms. Chicken farms are represented by blue dots. Quail farms are represented by green squares. Duck farms are represented by red triangles.
Figure 2. Map of Central Luzon, Philippines showing the physical location of all 1151 commercial poultry farms. Chicken farms are represented by blue dots. Quail farms are represented by green squares. Duck farms are represented by red triangles.
Vetsci 07 00107 g002
Figure 3. Risk map of highly pathogenic avian influenza transmission in Central Luzon, Philippines.
Figure 3. Risk map of highly pathogenic avian influenza transmission in Central Luzon, Philippines.
Vetsci 07 00107 g003
Table 1. Distribution of poultry farms in the provinces of Central Luzon.
Table 1. Distribution of poultry farms in the provinces of Central Luzon.
ProvincePoultry FarmsAverage Population
Aurora98222
Bataan9054,147
Bulacan15492,858
Nueva Ecija26857,823
Pampanga48527,257
Tarlac11287,231
Zambales33104,164
Table 2. Summary of farm-level basic reproductive ratio values.
Table 2. Summary of farm-level basic reproductive ratio values.
R0DucksQuailChickensTotal
0.01–0.5086705719
0.51–0.99181144163
1.00–1.50831223
1.51–2.00135015150
2.01–2.5911496

Share and Cite

MDPI and ACS Style

Salvador, R.; Tanquilut, N.; Na Lampang, K.; Chaisowwong, W.; Pfeiffer, D.; Punyapornwithaya, V. Identification of High-Risk Areas for the Spread of Highly Pathogenic Avian Influenza in Central Luzon, Philippines. Vet. Sci. 2020, 7, 107. https://doi.org/10.3390/vetsci7030107

AMA Style

Salvador R, Tanquilut N, Na Lampang K, Chaisowwong W, Pfeiffer D, Punyapornwithaya V. Identification of High-Risk Areas for the Spread of Highly Pathogenic Avian Influenza in Central Luzon, Philippines. Veterinary Sciences. 2020; 7(3):107. https://doi.org/10.3390/vetsci7030107

Chicago/Turabian Style

Salvador, Roderick, Neil Tanquilut, Kannika Na Lampang, Warangkhana Chaisowwong, Dirk Pfeiffer, and Veerasak Punyapornwithaya. 2020. "Identification of High-Risk Areas for the Spread of Highly Pathogenic Avian Influenza in Central Luzon, Philippines" Veterinary Sciences 7, no. 3: 107. https://doi.org/10.3390/vetsci7030107

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

Article Metrics

Back to TopTop