key: cord-005129-2e6b5tm2
authors: Chaudhry, Mamoona; Ahmad, Maqbool; Rashid, Hamad Bin; Sultan, Bakhat; Chaudhry, Haroon Rashid; Riaz, Aayesha; Shaheen, Muhammad Shabir
title: Prospective study of avian influenza H9 infection in commercial poultry farms of Punjab Province and Islamabad Capital Territory, Pakistan
date: 2016-10-20
journal: Trop Anim Health Prod
DOI: 10.1007/s11250-016-1159-6
sha: 
doc_id: 5129
cord_uid: 2e6b5tm2

A prospective study was conducted from November 2013 to February 2014 to estimate the spatial clustering; cumulative incidence and risk factors associated with avian influenza (AI) subtype H9 infection on commercial poultry farms of Pakistan. A total of 400 farms were enrolled and followed during the study period. Among these, 109 farms submitted samples suspected for AI to the laboratory, and only 47 farms were confirmed positive by hemagglutinin inhibition (HI) test. Data was collected from these 109 farms about their demography, management, and biosecurity practices. The cumulative incidence of H9N2 was 11.75 % (95 % confidence interval (CI) 8.76–15.23). The highest number of cases (40.42 %) was reported in January. One most likely cluster (p = 0.009, radius = 4.61 km) occurred in the Kasur district. Multivariable logistic regression analyses showed that the presence of wild birds on the farms (odds ratio (OR) = 16.18; 95 % CI 3.94–66.45) was independently associated with H9N2 infection. Cleaning of cages before delivery on farm (OR = 0.16; 95 % CI = 0.06–0.47), presence of a footbath at the entrance of farm (OR = 0.24; 95 % CI 0.08–0.79), and changing of gloves (OR = 0.33; 95 % CI 0.11–0.99) were protective factors against H9N2 infection. Reducing the exposure to risk factors and adapting biosecurity measures may reduce the risk of AI H9N2 infection on commercial poultry farms in Pakistan.

Influenza is a continuing threat to human and animal health. Every year, thousands of people are infected with seasonal influenza and may be exposed to subtypes of avian (H5, H6, H7, H9, and H10) and swine (H1 and H3) origin (García-Sastre and Schmolke 2014).

Avian influenza viruses (AIVs) of subtype H9N2 have spread widely since their first identification in turkeys in Wisconsin, USA, in 1966 (Homme and Easterday 1970) .

H9N2 viruses were isolated from pigs in 1998 and were subsequently isolated from humans with an influenza-like illness in both Hong Kong and Mainland China (Peiris et al. 1999) . H9N2 are significantly important due to their extensive circulation in domestic poultry in different regions of world from the Far East to the Middle East (Fusaro et al. 2011) . Genetic analysis of H9N2 viruses has showed extensive re-assortment of these viruses with many subtypes of AIVs including HPAI H5N1 and H7N3 viruses (Chaudhry et al. 2015; Fusaro et al. 2011) .

In Pakistan, commercial poultry production has attained the shape of an industry in recent years with investment of billions of rupees. Since 1995, AIV subtypes H9, H7, and H5 are responsible for five massive epidemics in Pakistan affecting poultry and poultry products across the country (Naeem et al. 2007) . Although H9N2 viruses are of low pathogenicity, the frequent heavy losses caused by them have raised serious concerns for the poultry industry in many countries.

Advancement has been made in disease investigations with new tools like geographic information system (GIS), which is used for spatiotemporal analysis of important emerging infections, e.g., severe respiratory syndrome (SARS), AIV H5N1, and influenza A (H1N1) (Tiensin et al. 2009; Martin et al. 2011; Lai et al. 2013) . Disease clustering can be detected by using space-time scan statistics (Kulldorff et al. 2005) .

Few studies have examined the association of risk factors with AI on poultry farms in Pakistan (Abbas et al. 2012; Chaudhry et al. 2015) . To date, very little information is available on spatial clustering of H9 infection in this region. Awareness about risk factors responsible for disease introduction and spatial clustering is critically important in developing risk-based surveillance strategies, policies, and timely recommendation for control. The primary objectives of this study were to calculate attack rate (AR) of H9N2 infection and to identify risk factors associated with this infection among poultry farms of Pakistan. The other objective was to identify any clustering of unusually high number of H9 cases than expected for early detection of any emerging outbreak of this disease in different areas of Pakistan when only the number of cases is available.

A prospective study was conducted from November 2013 to February 2014. All commercial poultry farms of Pakistan raising domesticated poultry for sale were considered as the target population of study. The final study population was commercial poultry farms submitting samples for laboratory analysis to the collaborating private poultry laboratory for routine screening and suspected infections. Each commercial farm was taken as a sampling unit. All poultry farms, which were included in the study, were considered negative for H9 at the start of study due to the absence of any influenza or influenzalike illness in the flock. None of broiler flock was vaccinated against H9, H5, or H7. Breeder and layers were vaccinated against H9.

A total of 400 commercial poultry farms of different production categories (breeders, broiler, and layer farms) located in Punjab Province and Islamabad Capital Territory of Pakistan were enrolled in the study. Out of these 400 farms, only 109 farms submitted samples to laboratory for suspected infection with AIV, and a pretested questionnaire was filled from the owner/supervisor of these 109 farms in a face-to-face interview after explaining the objectives of study to the farmers. Prior to interview, written consent of the owner/ attendant was obtained. The questionnaire contained 26 questions about risk factors, which were known to influence the disease occurrence and were selected after reviewing literature about AI (Nishiguchi et al. 1999; Ward et al. 2008; McQuiston et al. 2005; Fang et al. 2008; Woo and Park 2008; Abbas et al. 2012; Chaudhry et al. 2015; Nishiguchi et al. 2007 ) and from the observations of technical staff working on these farms. The farmers were requested to provide five to ten dead birds from total mortality on farm, which were carefully examined by conducting postmortem examination for specific disease lesions. Typical pathological lesions in respiratory system, i.e., rhinitis, sinusitis, congestion, and inflammation in the trachea (Swayne 2008) , were suspected for AIV. Confirmation was done by Anigen Rapid AIV Ag Detection Kit (BIONOTE Inc., Korea). The outcome of interest was H9 status, i.e., infected and non-infected farms. Samples confirmed by rapid test were further tested by virus isolation in embryonated chicken eggs, and subtyping was done by hemagglutinin inhibition (HI) test.

AR of H9N2 was calculated (Thrusfield 2007) . All biologically plausible and relevant variables were screened in univariable analysis by using glm function of the epicalc package (version, 2.15.1.0) in R statistical software (available at http://www.R-project.org). A multivariable model was derived by forward stepwise selection procedure (Dohoo et al. 2003) . Variables with significant univariable relationship at P < 0.25 were selected for inclusion in the final model. Odds ratios (ORs) and corresponding 95 % confidence intervals (CIs) were calculated (Hosmer and Lemeshow 2000) .

All laboratory-confirmed cases of H9N2 between December 2013 to February 2014 were geocoded to street addresses. SatScan software version 9.1.1 developed by Martin Kulldorff, Havard Medical School (Boston, USA) prospective space-time permutation scan statistic module was used to detect a local excess of events and to test if this excess could have occurred by chance (Hyder et al. 2011 ). This method consists of thousands of cylinders that move across space and/or time. Each cylinder has a base, which represents geographical area (in this study, a commercial farm), and height, which is time (in this study, a day). The base of each cylinder comprised a maximum of 50 % of the population, while height was a maximum of 50 % of the study time (60 days). The cylinder with more observed cases than expected, with respect to cases reported outside the cylinder, is called Bmost likely cluster.^For each location and size of the cylinder, the number of observed and expected cases is counted. Among these, the most Bunusual^excess of observed cases is noted. The statistical significance of this cluster is then evaluated taking into account the multiple testing stemming from the many potential cluster locations and sizes evaluated (Kulldorff et al. 2005 ). ArcGIS version 10 was used for the map display.

During the study, we followed 400 commercial poultry farms of which 109 submitted samples. Among these 109 farms, 47 ). Among the infected farms, majority (66 %) reported 5-9 % mortality due to H9 (Table 1) . The AR was highest in Lahore district (11/28) followed by Kasur (24/38) and Sheikhupura (6/13) districts (Fig. 1) . No sample was positive for Newcastle disease virus, H5, and H7 AIVs.

The study showed a high AR in the month of January (21/ 34) followed by December (11/29), November (9/23), and lowest incidence in February (6/23) (Fig. 2) .

Out of 26, 11 factors were selected for inclusion in final model of multivariable analysis (Table 2) . Factors with p > 0.25 were excluded from further analysis.

In the final multivariable model, four factors were identified as significant (Table 3) . Among those four factors, one factor was identified as risk factor (OR > 1), i.e., wild birds on the farm, and three factors were proved to be protective factor (OR < 1), namely cleaning of cages before entering the farm area, having foot bath/dipping area at the entrance of farm, and workers change gloves while entry into bird area. 

Total examined farms in Punjab Province were 107, while 2 farms were examined from Islamabad Capital Territory (Fig. 3) .

From November 1, 2013 to February 28, 2014, one most likely cluster (p = 0.009, radius = 4.61 km) occurred in the Kasur district of Punjab, Pakistan (Fig. 4) . This signal had four cases observed over 25 days when 0.52 cases were expected [relative risk (RR) = 7.67], with a null occurrence rate of once every 111 days.

Attack rate was highest in Lahore district (17/28) followed by Kasur (24/38) and Sheikhupura (6/13) districts. The reason for this high AR could be the high density of commercial poultry farms in these districts. Association of HPAI H5N1 and high and medium density of poultry farms has been studied and was identified as a risk factor in different countries (Henning et al. 2009; Abbas et al. 2012) . Poultry farm densities in these regions ranged from 0.05 to 4 farms per square kilometer (Hamilton et al. 2009 ). As distance between the farms plays a significant role in the transmission of infection, poultry farms in these areas are exposed almost every year during the endemic (Abbas et al. 2012; Chaudhry et al. 2015) . Furthermore, in densely populated poultry areas, movement of vehicles and people from farm to farm is considerably high, subsequently facilitating the spread of this virus through fomites. In the current study, the two districts (Lahore and Kasur) are on the main highway road (Grand Trunk road), which have a significant number of commercial poultry farms alongside. The vehicles over loaded with infected birds or mortality move on this road routinely. Previously, many studies have demonstrated that proximity to major roads was associated with avian influenza (Ward et al. 2008; Chaudhry et al. 2015) . Movement of veterinarian and para-veterinary staff between different farms to implement control measures or to investigate mortalities can also contribute to the spread of virus among farms. In Pakistan, poultry farmers who apply strict biosecurity measures sometimes relax these rules for visitors, who enter the bird areas (Chaudhry et al. 2015) .

Spatiotemporal analysis detected a most likely cluster in Kasur district (p = 0.009, radius = 4.61 km). Poultry farm in this cluster has a significantly higher risk of being infected as compared to poultry farms outside the cluster (RR = 7.67). These results represent the most important hotspot of expected outbreak and are valuable for improving knowledge and understanding of spatial pattern of H9N2 in specific areas. Targeted surveillance of these districts is needed for early detection of any future AI outbreak and their possible re-assortment. Spatial analyses have been used to study different outbreaks of human and avian influenza infections (Tiensin et al. 2009; Leveau et al. 2015) .

Risk factor analyses showed that the presence of wild birds on farm could enhance the probability of infection (OR = 16.18; 95 % CI 3.94-66.45). Wild birds could serve as a potential source of propagation of AI virus especially when biosecurity measures are poorly implemented on farm. They can act both as mechanical and biological vectors (shedding the virus in droppings) and are important source of introducing virus to new areas (McQuiston et al. 2005; Henning et al. 2009; Chaudhry et al. 2015) .

Cleaning of cages before delivery was strongly associated with decrease in risk of H9 (OR = 0.16, 95 % CI 0.06-0.47). Though statistically non-significant, cleaning of cages has been studied as a protective factor in lowering risk of H9 infection (Chaudhry et al. 2015) . Poultry traders are well aware about the importance of cleaning cages and vehicles as an effective biosecurity measure (Kurscheid et al. 2015) .

Presence of a footbath/dipping area at the entrance of farm and changing of gloves are also important as part of biosecurity measures on farm and have proved effective in decreasing risk of AI previously (Biswas et al. 2009; Chaudhry et al. 2015) .

This study found evidence of clustering, in space and time, and identified some well-known factors mainly responsible for increasing risk of AIV infection. Enhancing good management practices and strict biosecurity can lower the risk of infection among poultry farms. Spatial clustering of disease provides information to health authorities to more effectively target and improve their surveillance and control strategies in affected areas.

Acknowledgments The authors are highly indebted to respondents of the commercial farms, who participated in study for data collection. We also acknowledge the support of staff of commercial poultry laboratory for providing us data from logbook about case farms, which were included in the study. Without their great cooperation and help, this study would not have been possible.

Ethical approval This article does not contain any studies with animals performed by any of the authors.

Informed consent The manuscript does not contain clinical studies or patient data. The owners of the commercial farms were briefed about the objective of study, and informed consent was obtained from all individual participants included in the study to collect data.

The authors declare that they have no conflict of interest.

Contact structure and potential risk factors for avian influenza transmission among open-sided chicken farms in Kamalia, an important poultry rearing area of Pakistan

Risk factors for infection with highly pathogenic influenza A virus (H5N1) in commercial chickens in Bangladesh

A case-control study to identify risk factors associated with avian influenza subtype H9N2 on commercial poultry farms in Pakistan

Veterinary epidemiologic research

Environmental factors contributing to the spread of H5N1 avian influenza in Mainland China

Phylogeography and evolutionary history of reassortant H9N2 viruses with potential human health implications

Avian influenza A H10N8?a virus on the verge?

Are the Australian poultry industries vulnerable to large outbreaks of highly pathogenic avian influenza?

Farm-and flock-level risk factors associated with highly pathogenic avian influenza outbreaks on small-holder duck and chicken farms in the Mekong Delta of Viet Nam

Avian influenza virus infections: I. Characteristics of influenza A-turkey-Wisconsin-1966 virus

Applied logistic regression

Use of spatiotemporal analysis of laboratory submission data to identify potential outbreaks of new or emerging diseases in cattle in Great Britain

A space-time permutation scan statistics for disease outbreak detection

Knowledge and perceptions of highly pathogenic avian influenza (HPAI) among poultry traders in live bird markets in Bali and Lombok

Spatio-temporal and stochastic modeling of the severe acute respiratory syndrome (SARS)

Spatiotemporal trends of cases of pandemic influenza A (H1N1)pdm09 in Argentina

Spatial distribution and risk factors of highly pathogenic avian influenza (HPAI) H5N1 in China

Evaluation of risk factors for the spread of low pathogenicity H7N2 avian influenza virus among commercial poultry farms

Avian influenza in Pakistan: outbreaks of low-and high-pathogenicity avian influenza in Pakistan during

Risk factors for the introduction of avian influenza virus into commercial layer chicken farms during the outbreaks caused by a low-pathogenic H5N2 virus in Japan in 2005

Human infection with influenza H9N2

Avian influenza

Veterinary Epidemiology

Ecologic risk factor investigation of clusters of avian influenza A (H5N1) virus infection in Thailand

Environmental and anthropogenic risk factors for highly pathogenic avian influenza subtype H5N1 outbreaks in Romania

Seroprevalence of low pathogenic avian influenza (H9N2) and associated risk factors in the Gyeonggi-do of Korea during


