# Prevalence of Avian Influenza A(H5) and A(H9) Viruses in Live Bird Markets, Bangladesh

Younjung Kim, Paritosh K Biswas, Mohammad Giasuddin, Mahmudul Hasan, Rashed Mahmud, Yu-Mei Chang, Steve Essen, Mohammed A Samad, Nicola S Lewis, Ian H Brown, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Younjung Kim, Paritosh K Biswas, Mohammad Giasuddin, Mahmudul Hasan, Rashed Mahmud, et al. Prevalence of Avian Influenza A(H5) and A(H9) Viruses in Live Bird Markets, Bangladesh. Emerging Infectious Diseases, 2018, 24 (12), pp.2309-2316. 10.3201/eid2412.180879. hal-04598004

HAL Id: hal-04598004

https://hal.science/hal-04598004

Submitted on 3 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# Prevalence of Avian Influenza A(H5) and A(H9) Viruses in Live Bird Markets, Bangladesh

Younjung Kim,<fnref n="1" /> Paritosh K. Biswas,<fnref n="1" /> Mohammad Giasuddin,<fnref n="1" /> Mahmudul Hasan, Rashed Mahmud, Yu-Mei Chang, Steve Essen, Mohammed A. Samad, Nicola S. Lewis, Ian H. Brown, Natalie Moyen, Md. Ahasanul Hoque, Nitish C. Debnath, Dirk U. Pfeiffer, Guillaume Fournié

We conducted a cross-sectional study in live bird markets (LBMs) in Dhaka and Chittagong, Bangladesh, to estimate the prevalence of avian influenza A(H5) and A(H9) viruses in different types of poultry and environmental areas by using Bayesian hierarchical logistic regression models. We detected these viruses in nearly all LBMs. Prevalence of A(H5) virus was higher in waterfowl than in chickens, whereas prevalence of A(H9) virus was higher in chickens than in waterfowl and, among chicken types, in industrial broilers than in cross-breeds and in digenous breeds. LBMs with ≥1 wholesaler were more frequently contaminated by A(H5) virus than retail-only LBMs. Prevalence of A(H9) virus in poultry and level of environmental contamination were also higher in LBMs with ≥1 wholesaler. We found a high level of circulation of both avian influenza viruses in surveyed LBMs. Prevalence was influenced by type of poultry, environmental site, and trading patterns.

Low pathogenicity avian influenza A(H9N2) virus and highly pathogenic avian influenza A(H5N1) virus are endemic in poultry populations in Bangladesh (1–4). In addition to their adverse effect on poultry production, these viruses have resulted in sporadic influenza cases in humans (2,3). Because there is potential for generating novel reassortant variants between them or with other virus subtypes, their persistent circulation in poultry poses a serious threat to animal and human health globally (5–9).

Author affiliations: The Royal Veterinary College, Hatfield, UK (Y. Kim, Y.-M. Chang, N. Moyen, D.U. Pfeiffer, G. Fournié); City University of Hong Kong, Hong Kong, China (Y. Kim, D.U. Pfeiffer); Chittagong Veterinary and Animal Sciences University, Chittagong, Bangladesh (P.K. Biswas, R. Mahmud, M.A. Hoque); Bangladesh Livestock Research Institute, Dhaka, Bangladesh (M. Giasuddin, M. Hasan, M.A. Samad); Animal Health and Plant Agency, Weybridge, UK (S. Essen, N.S. Lewis, I.H. Brown); Food and Agriculture Organization of the United Nations, Dhaka (N.C. Debnath)

DOI: https://doi.org/10.3201/eid2412.180879

Live bird markets (LBMs) form the backbone of poultry trade in many countries in Asia. Birds of different types and from different geographic areas are introduced daily into LBMs and might be caged together, promoting local transmission of multiple virus subtypes and generating opportunities for reassortment (10–12). Surveys and routine surveillance have described the abundance and diversity of avian influenza A viruses (AIVs) in LBMs in AIV-endemic countries, including Bangladesh (1,4,11,13–21). However, only the proportion of positive samples is usually reported, without accounting for the hierarchical data structure, especially the clustering of sampled poultry per LBM. Therefore, a robust assessment of AIV prevalence in LBMs is lacking, although this knowledge is essential to understand AIV epidemiology and optimize surveillance design.

Multiple poultry species and, for each poultry species, multiple breeds are offered for sale in LBMs in Bangladesh. Desi, Sonali, and broiler are the most commonly traded chicken types. Desi, which means “local” in Ben-gali, are indigenous chicken breeds raised in backyard farms. Sonali is a cross-breed of the Rhode Island Red cocks and Fayoumi hens. Broilers are industrial white-feathered breeds. In addition to varying levels of susceptibility, different poultry types might be raised in different farming systems and traded through different value chains (i.e., the range of activities that businesses perform to deliver products to customers), therefore being exposed to different pathogens and pathogen loads (22,23). However, the proportion of AIV-positive samples is generally reported as an overall estimate or stratified only by poultry species. Likewise, the relative position of individual LBMs in a regional or national live poultry trading network might also affect AIV prevalence; the sources from which poultry are supplied to traders and the time they spend in LBMs influence the likelihood of virus introduction and amplification in LBMs and vary depending on traders being wholesalers or retailers (12,24).

<footnote n="1">¹These authors contributed equally to this article.</footnote>
# Prevalence of Avian Influenza A(H5) and A(H9) Viruses in Live Bird Markets, Bangladesh

Younjung Kim,<fnref n="1" /> Paritosh K. Biswas,<fnref n="1" /> Mohammad Giasuddin,<fnref n="1" /> Mahmudul Hasan, Rashed Mahmud, Yu-Mei Chang, Steve Essen, Mohammed A. Samad, Nicola S. Lewis, Ian H. Brown, Natalie Moyen, Md. Ahasanul Hoque, Nitish C. Debnath, Dirk U. Pfeiffer, Guillaume Fournié

We conducted a cross-sectional study in live bird markets (LBMs) in Dhaka and Chittagong, Bangladesh, to estimate the prevalence of avian influenza A(H5) and A(H9) viruses in different types of poultry and environmental areas by using Bayesian hierarchical logistic regression models. We detected these viruses in nearly all LBMs. Prevalence of A(H5) virus was higher in waterfowl than in chickens, whereas prevalence of A(H9) virus was higher in chickens than in waterfowl and, among chicken types, in industrial broilers than in cross-breeds and in digenous breeds. LBMs with ≥1 wholesaler were more frequently contaminated by A(H5) virus than retail-only LBMs. Prevalence of A(H9) virus in poultry and level of environmental contamination were also higher in LBMs with ≥1 wholesaler. We found a high level of circulation of both avian influenza viruses in surveyed LBMs. Prevalence was influenced by type of poultry, environmental site, and trading patterns.

Low pathogenicity avian influenza A(H9N2) virus and highly pathogenic avian influenza A(H5N1) virus are endemic in poultry populations in Bangladesh (1–4). In addition to their adverse effect on poultry production, these viruses have resulted in sporadic influenza cases in humans (2,3). Because there is potential for generating novel reassortant variants between them or with other virus subtypes, their persistent circulation in poultry poses a serious threat to animal and human health globally (5–9).

Author affiliations: The Royal Veterinary College, Hatfield, UK (Y. Kim, Y.-M. Chang, N. Moyen, D.U. Pfeiffer, G. Fournié); City University of Hong Kong, Hong Kong, China (Y. Kim, D.U. Pfeiffer); Chittagong Veterinary and Animal Sciences University, Chittagong, Bangladesh (P.K. Biswas, R. Mahmud, M.A. Hoque); Bangladesh Livestock Research Institute, Dhaka, Bangladesh (M. Giasuddin, M. Hasan, M.A. Samad); Animal Health and Plant Agency, Weybridge, UK (S. Essen, N.S. Lewis, I.H. Brown); Food and Agriculture Organization of the United Nations, Dhaka (N.C. Debnath)

DOI: https://doi.org/10.3201/eid2412.180879

Live bird markets (LBMs) form the backbone of poultry trade in many countries in Asia. Birds of different types and from different geographic areas are introduced daily into LBMs and might be caged together, promoting local transmission of multiple virus subtypes and generating opportunities for reassortment (10–12). Surveys and routine surveillance have described the abundance and diversity of avian influenza A viruses (AIVs) in LBMs in AIV-endemic countries, including Bangladesh (1,4,11,13–21). However, only the proportion of positive samples is usually reported, without accounting for the hierarchical data structure, especially the clustering of sampled poultry per LBM. Therefore, a robust assessment of AIV prevalence in LBMs is lacking, although this knowledge is essential to understand AIV epidemiology and optimize surveillance design.

Multiple poultry species and, for each poultry species, multiple breeds are offered for sale in LBMs in Bangladesh. Desi, Sonali, and broiler are the most commonly traded chicken types. Desi, which means “local” in Ben-gali, are indigenous chicken breeds raised in backyard farms. Sonali is a cross-breed of the Rhode Island Red cocks and Fayoumi hens. Broilers are industrial white-feathered breeds. In addition to varying levels of susceptibility, different poultry types might be raised in different farming systems and traded through different value chains (i.e., the range of activities that businesses perform to deliver products to customers), therefore being exposed to different pathogens and pathogen loads (22,23). However, the proportion of AIV-positive samples is generally reported as an overall estimate or stratified only by poultry species. Likewise, the relative position of individual LBMs in a regional or national live poultry trading network might also affect AIV prevalence; the sources from which poultry are supplied to traders and the time they spend in LBMs influence the likelihood of virus introduction and amplification in LBMs and vary depending on traders being wholesalers or retailers (12,24).

<footnote n="1">¹These authors contributed equally to this article.</footnote>
However, such information is generally poorly document-
ed, or even ignored.

To address these issues, we conducted a cross-section-
al survey in the 2 largest cities in Bangladesh, Dhaka and
Chittagong, during February–March 2016. First, we esti-
mated prevalence of influenza A(H5) and A(H9) viruses
in marketed poultry and the LBM environment. We also
accounted for the clustering effect at LBM level by using
Bayesian hierarchical logistic regression models. Second,
we assessed the effect of type of poultry and environmental
site, and the position of LBMs in the poultry value chain on
AIV prevalence.

# Materials and Methods

Sample Collection

An LBM was defined as an open space in which ≥2 poultry stalls sell live poultry at least once a week, and only those selling >400 poultry/day were considered eligible for this study. We aimed to sample 40 LBMs, and from each of these LBMs, 60 birds and 50 environmental sites (sample size calculations in online Technical Appendix 1, https://wwwnc.cdc.gov/EID/article/24/12/18-0879-Techapp1.pdf). We used a stratified cluster sampling design. For poultry, LBMs, stalls within selected LBMs, and birds within selected stalls constituted the primary, secondary, and tertiary sampling units, respectively. For environmental sites, LBMs constituted primary sampling units and environmental sites within selected LBMs constituted secondary sampling units.

We stratified LBMs by city for Dhaka and Chittagong and, within each city, by poultry sales into large and small LBMs, hypothesizing that the risk for AIV infection varies between geographic locations and the number of poultry traded. Also, simple random sampling with too small a sample size of LBMs was not likely to capture diversity of LBM types because the distribution of LBMs as a function of their size tended to be right-skewed; the largest LBMs were often wholesale markets (24,25). We hypothesized that samples of different origins have different AIV prevalences and thus stratified birds and environmental sites into 5 types of poultry and 10 types of environmental sites commonly found to be contaminated with AIV (26) (online Technical Appendix 2 Tables 1, 2, https://wwwnc.cdc.gov/EID/article/24/12/18-0879-Techapp2.pdf).

At the first sampling stage, we sampled 40 LBMs.
The number of LBMs selected in Dhaka (n = 26) and Chit-
tagong (n = 14) was proportional to the number of LBMs
eligible in each city (n = 80 for Dhaka and n = 36 for Chit-
tagong). In each city, we further stratified LBMs by size:
50% of the selected LBMs were large, trading the highest
number of poultry (13 largest LBMs in Dhaka and 7 largest
LBMs in Chittagong); 50% were small, randomly selected

from the bottom 50% of eligible LBMs in terms of number
of poultry traded.

At the second sampling stage, we randomly selected stalls and environmental sites in each LBM independently for each type of poultry and environmental site. We created a list of stalls selling each poultry type for each LBM. We then selected stalls from these lists by using a random number generator. Likewise, for each type of environmental site, we selected sites from a list of sites identified in each LBM by using a random number generator. We collected 1 swab specimen from each environmental site. We pooled 5 swab specimens collected from the same LBM and site type.

At the third sampling stage, for each poultry type, we randomly selected 5 birds from each of the stalls selected for that type and collected cloacal and oropharyngeal swab specimens from each of the selected birds. We pooled 5 swab specimens collected from the same stall and poultry type separately for cloacal and oropharyngeal swab specimens. We transported samples collected in Chittagong on the day of sampling to the Chittagong Veterinary and Animal Sciences University (Chittagong) and samples collected in Dhaka on the day of sampling to the Bangladesh Livestock Research Institute (Dhaka). Samples were stored at -80°C until diagnostic laboratory processing.

Sample Screening

We screened pools for AIVs by using a real-time reverse
transcription PCR (RT-PCR) and specific primers and
probes (27,28). We extracted virus RNA by using the Mag-
MAX RNA Isolation Kit (QIAGEN, Hilden, Germany) and
reverse transcribed and amplified virus RNA by using the
AgPath-ID One-Step RT-PCR (ThermoFisher Scientific,
Waltham, MA, USA). We then screened a pool with a cy-
cle threshold (C<sub>t</sub>) <40 for the AIV matrix gene for the H5
and H9 genes. Results were considered positive for the H5
subtype if C<sub>t</sub> <38 and positive for the H9 subtype if C<sub>t</sub>
(27,28). A pool was considered positive for AIV if its C<sub>t</sub>
for the AIV matrix gene <38 or if it was positive for any
of the H5 and H9 subtypes. A given group of 5 birds was
considered positive if any of its cloacal and oropharyngeal
pools showed a positive result.

Bayesian Hierarchical Logistic Regression Models

We developed 2-level Bayesian hierarchical logistic regres-
sion models to estimate LBM-level, bird-level, and environ-
mental swab specimen-level prevalence from pooled swab
samples, accounting for lower-level (swab specimens) and
higher-level (LBMs) risk factors. We developed separate
models for poultry and environmental samples to avoid pa-
rameters related to different sampling units interfering with
each other (29). We used LBM type (retail or mixed), city
(Chittagong or Dhaka), and size (small or large) as LBM-
However, such information is generally poorly document-
ed, or even ignored.

To address these issues, we conducted a cross-section-
al survey in the 2 largest cities in Bangladesh, Dhaka and
Chittagong, during February–March 2016. First, we esti-
mated prevalence of influenza A(H5) and A(H9) viruses
in marketed poultry and the LBM environment. We also
accounted for the clustering effect at LBM level by using
Bayesian hierarchical logistic regression models. Second,
we assessed the effect of type of poultry and environmental
site, and the position of LBMs in the poultry value chain on
AIV prevalence.

Materials and Methods

Sample Collection

An LBM was defined as an open space in which ≥2 poultry stalls sell live poultry at least once a week, and only those selling >400 poultry/day were considered eligible for this study. We aimed to sample 40 LBMs, and from each of these LBMs, 60 birds and 50 environmental sites (sample size calculations in online Technical Appendix 1, https://wwwnc.cdc.gov/EID/article/24/12/18-0879-Techapp1.pdf). We used a stratified cluster sampling design. For poultry, LBMs, stalls within selected LBMs, and birds within selected stalls constituted the primary, secondary, and tertiary sampling units, respectively. For environmental sites, LBMs constituted primary sampling units and environmental sites within selected LBMs constituted secondary sampling units.

We stratified LBMs by city for Dhaka and Chittagong and, within each city, by poultry sales into large and small LBMs, hypothesizing that the risk for AIV infection varies between geographic locations and the number of poultry traded. Also, simple random sampling with too small a sample size of LBMs was not likely to capture diversity of LBM types because the distribution of LBMs as a function of their size tended to be right-skewed; the largest LBMs were often wholesale markets (24,25). We hypothesized that samples of different origins have different AIV prevalences and thus stratified birds and environmental sites into 5 types of poultry and 10 types of environmental sites commonly found to be contaminated with AIV (26) (online Technical Appendix 2 Tables 1, 2, https://wwwnc.cdc.gov/EID/article/24/12/18-0879-Techapp2.pdf).

At the first sampling stage, we sampled 40 LBMs.
The number of LBMs selected in Dhaka (n = 26) and Chit-
tagong (n = 14) was proportional to the number of LBMs
eligible in each city (n = 80 for Dhaka and n = 36 for Chit-
tagong). In each city, we further stratified LBMs by size:
50% of the selected LBMs were large, trading the highest
number of poultry (13 largest LBMs in Dhaka and 7 largest
LBMs in Chittagong); 50% were small, randomly selected

from the bottom 50% of eligible LBMs in terms of number
of poultry traded.

At the second sampling stage, we randomly selected stalls and environmental sites in each LBM independently for each type of poultry and environmental site. We created a list of stalls selling each poultry type for each LBM. We then selected stalls from these lists by using a random number generator. Likewise, for each type of environmental site, we selected sites from a list of sites identified in each LBM by using a random number generator. We collected 1 swab specimen from each environmental site. We pooled 5 swab specimens collected from the same LBM and site type.

At the third sampling stage, for each poultry type, we randomly selected 5 birds from each of the stalls selected for that type and collected cloacal and oropharyngeal swab specimens from each of the selected birds. We pooled 5 swab specimens collected from the same stall and poultry type separately for cloacal and oropharyngeal swab specimens. We transported samples collected in Chittagong on the day of sampling to the Chittagong Veterinary and Animal Sciences University (Chittagong) and samples collected in Dhaka on the day of sampling to the Bangladesh Livestock Research Institute (Dhaka). Samples were stored at -80°C until diagnostic laboratory processing.

Sample Screening

We screened pools for AIVs by using a real-time reverse
transcription PCR (RT-PCR) and specific primers and
probes (27,28). We extracted virus RNA by using the Mag-
MAX RNA Isolation Kit (QIAGEN, Hilden, Germany) and
reverse transcribed and amplified virus RNA by using the
AgPath-ID One-Step RT-PCR (ThermoFisher Scientific,
Waltham, MA, USA). We then screened a pool with a cy-
cle threshold (C<sub>t</sub>) <40 for the AIV matrix gene for the H5
and H9 genes. Results were considered positive for the H5
subtype if C<sub>t</sub> <38 and positive for the H9 subtype if C<sub>t</sub>
(27,28). A pool was considered positive for AIV if its C<sub>t</sub>
for the AIV matrix gene <38 or if it was positive for any
of the H5 and H9 subtypes. A given group of 5 birds was
considered positive if any of its cloacal and oropharyngeal
pools showed a positive result.

Bayesian Hierarchical Logistic Regression Models

We developed 2-level Bayesian hierarchical logistic regres-
sion models to estimate LBM-level, bird-level, and environ-
mental swab specimen-level prevalence from pooled swab
samples, accounting for lower-level (swab specimens) and
higher-level (LBMs) risk factors. We developed separate
models for poultry and environmental samples to avoid pa-
rameters related to different sampling units interfering with
each other (29). We used LBM type (retail or mixed), city
(Chittagong or Dhaka), and size (small or large) as LBM-
level risk factors. We defined an LBM with only retailers (i.e., trader selling poultry to end-users only) as retail and an LBM with ≥1 wholesaler (i.e., trader selling poultry to other traders) as mixed. Each LBM-level risk factor (LBM type, city, and size) was assessed separately because we could not include them simultaneously in any given model due to the small number of sampled LBMs. In models for poultry samples, we differentiated birds into 1) chicken and waterfowl or 2) broiler, Desi, Sonali (i.e., chicken types), and waterfowl. Waterfowl were not differentiated further because of the small number of pools collected from ducks and geese. In models for environmental samples, we differentiated environmental sites into stall and slaughter areas or classified as environmental area without differentiation. We ran models (online Technical Appendix 1) by using a Markov Chain Monte Carlo simulation in JAGS (30) and R.3.4.2 (31).

# Results

## Descriptive Results for Pooled Swab Samples

We collected 477 pairs of cloacal and oropharyngeal pooled samples from 2,384 birds, and 400 environmental pooled samples from 2,000 environmental sites in 40 LBMs in Chittagong and Dhaka. Each pool was composed of 5 swab specimens, except for 1 pair of cloacal and oropharyngeal pools made from 4 swab specimens collected from geese. We collected 12 pairs of cloacal and oropharyngeal pooled samples from all LBMs, except for 11 pairs from 3 LBMs. We sampled chickens in all LBMs (8–12 pairs/LBM), and waterfowl in 25 LBMs (0–4 pairs/LBM). Broilers accounted for most samples (32.1%), followed by Desi (26.6%) and Sonali (25.6%). Ducks accounted for 76% of 75 pool pairs collected from waterfowl and geese accounted for 24%. We collected 10 environmental pools in each LBM (stall areas: 4–8 pools, slaughter areas, 2–6 pools).

Of 47.4% (416/877) pools considered positive for AIV, 6.5% pools were negative for the AIV matrix gene but positive for any of the H9 and H5 subtypes. The H9 subtype (63.2% positive pools) was detected more frequently than the H5 subtype (21.6%), and 12.3% of pools were positive for both subtypes and 27.4% of pools were negative for both subtypes. Although 80.0% of the LBMs had >1 A(H5) virus-positive poultry or environmental pool, 97.5% had ≥1 A(H9) virus-positive poultry or environmental pool. We determined the prevalence of pools that were positive for A(H5) and A(H9) viruses according to sample and LBM type (Table 1).

Approximately 33.3% of pools collected from waterfowl were positive for A(H5) virus, whereas only 5.5% of those collected from chickens were positive. In contrast, the prevalence of A(H9) virus-positive pools was higher in chickens (36.3%) than in waterfowl (18.7%). Among

waterfowl, ducks (19.3%) and geese (16.7%) had a similar prevalence of A(H9) virus-positive pools, but the prevalence of A(H5) virus-positive pools was higher in ducks (36.8%) than in geese (22.2%). For both H5 and H9 subtypes, the prevalence of positive pools was higher for oropharyngeal samples (8.6% for H5 and 31.9% for H9) than for cloacal samples (3.6% and 9.9%) in all surveyed poultry types (online Technical Appendix 2 Table 1).

Approximately 25% of environmental pools were positive for A(H9) virus, and the prevalence of positive pools was higher in slaughter areas (31.5%), especially knives and boards used for slaughter and processing, than stall areas (20.2%). The prevalence of A(H5) virus-positive environmental pools was lower (10.8%) and did not vary between slaughter and stall areas (online Technical Appendix 2 Table 2).

## Bayesian Model Results

Convergence was achieved for all models; the Gelman and Rubin statistic was <1.001 and the effective sample size was >10,000 for all parameters. For each AIV subtype, the best models reasonably predicted the number of positive pools (online Technical Appendix 2 Figure 1). In the best H5 models (i.e., lowest deviance information criterion), A(H5) virus prevalence differed according to poultry species (chicken, waterfowl), but not according to the type of environmental site. In contrast, in the best H9 models, A(H9) virus prevalence differed according to type of poultry (broiler, Desi, Sonali, waterfowl) and environmental site (slaughter and stall area). For both subtypes, LBM size and city did not improve model fit when compared with LBM type. For ease of comparison between the 2 AIV subtypes, we report LBM-level, bird-level, and environmental swab specimen-level prevalences of A(H5) and A(H9) viruses on the basis of the best H9 models with LBM type (Tables 1, 2). This reporting did not affect interpretation of results, and we provide estimates obtained with more parsimonious models (online Technical Appendix 2 Tables 3–6).

LBM-level A(H5) virus prevalence was lower in retail LBMs than in mixed LBMs, and the posterior median estimate was ≈100% for mixed LBMs. However, among contaminated LBMs, levels of virus detection in birds and environmental areas did not vary between LBM types, but A(H5) virus prevalence in waterfowl was ≈6 times higher than in chickens (Figure). The prevalence did not vary between chicken breeds or environmental areas.

In contrast to that for A(H5) virus, we found that the posterior median estimate of the LBM-level A(H9) virus prevalence was ≈100% for retail and mixed LBM groups, but the level of virus detection in birds and environmental areas was higher for mixed LBMs than for retail LBMs. A(H9) virus prevalence was highest in broilers and lowest in waterfowl. The prevalence in broilers was 3.8 times as
level risk factors. We defined an LBM with only retailers (i.e., trader selling poultry to end-users only) as retail and an LBM with ≥1 wholesaler (i.e., trader selling poultry to other traders) as mixed. Each LBM-level risk factor (LBM type, city, and size) was assessed separately because we could not include them simultaneously in any given model due to the small number of sampled LBMs. In models for poultry samples, we differentiated birds into 1) chicken and waterfowl or 2) broiler, Desi, Sonali (i.e., chicken types), and waterfowl. Waterfowl were not differentiated further because of the small number of pools collected from ducks and geese. In models for environmental samples, we differentiated environmental sites into stall and slaughter areas or classified as environmental area without differentiation. We ran models (online Technical Appendix 1) by using a Markov Chain Monte Carlo simulation in JAGS (30) and R.3.4.2 (31).

## Results

### Descriptive Results for Pooled Swab Samples

We collected 477 pairs of cloacal and oropharyngeal pooled samples from 2,384 birds, and 400 environmental pooled samples from 2,000 environmental sites in 40 LBMs in Chittagong and Dhaka. Each pool was composed of 5 swab specimens, except for 1 pair of cloacal and oropharyngeal pools made from 4 swab specimens collected from geese. We collected 12 pairs of cloacal and oropharyngeal pooled samples from all LBMs, except for 11 pairs from 3 LBMs. We sampled chickens in all LBMs (8–12 pairs/LBM), and waterfowl in 25 LBMs (0–4 pairs/LBM). Broilers accounted for most samples (32.1%), followed by Desi (26.6%) and Sonali (25.6%). Ducks accounted for 76% of 75 pool pairs collected from waterfowl and geese accounted for 24%. We collected 10 environmental pools in each LBM (stall areas: 4–8 pools, slaughter areas, 2–6 pools).

Of 47.4% (416/877) pools considered positive for AIV, 6.5% pools were negative for the AIV matrix gene but positive for any of the H9 and H5 subtypes. The H9 subtype (63.2% positive pools) was detected more frequently than the H5 subtype (21.6%), and 12.3% of pools were positive for both subtypes and 27.4% of pools were negative for both subtypes. Although 80.0% of the LBMs had >1 A(H5) virus-positive poultry or environmental pool, 97.5% had ≥1 A(H9) virus-positive poultry or environmental pool. We determined the prevalence of pools that were positive for A(H5) and A(H9) viruses according to sample and LBM type (Table 1).

Approximately 33.3% of pools collected from waterfowl were positive for A(H5) virus, whereas only 5.5% of those collected from chickens were positive. In contrast, the prevalence of A(H9) virus-positive pools was higher in chickens (36.3%) than in waterfowl (18.7%). Among

waterfowl, ducks (19.3%) and geese (16.7%) had a similar prevalence of A(H9) virus-positive pools, but the prevalence of A(H5) virus-positive pools was higher in ducks (36.8%) than in geese (22.2%). For both H5 and H9 subtypes, the prevalence of positive pools was higher for oropharyngeal samples (8.6% for H5 and 31.9% for H9) than for cloacal samples (3.6% and 9.9%) in all surveyed poultry types (online Technical Appendix 2 Table 1).

Approximately 25% of environmental pools were positive for A(H9) virus, and the prevalence of positive pools was higher in slaughter areas (31.5%), especially knives and boards used for slaughter and processing, than stall areas (20.2%). The prevalence of A(H5) virus-positive environmental pools was lower (10.8%) and did not vary between slaughter and stall areas (online Technical Appendix 2 Table 2).

### Bayesian Model Results

Convergence was achieved for all models; the Gelman and Rubin statistic was <1.001 and the effective sample size was >10,000 for all parameters. For each AIV subtype, the best models reasonably predicted the number of positive pools (online Technical Appendix 2 Figure 1). In the best H5 models (i.e., lowest deviance information criterion), A(H5) virus prevalence differed according to poultry species (chicken, waterfowl), but not according to the type of environmental site. In contrast, in the best H9 models, A(H9) virus prevalence differed according to type of poultry (broiler, Desi, Sonali, waterfowl) and environmental site (slaughter and stall area). For both subtypes, LBM size and city did not improve model fit when compared with LBM type. For ease of comparison between the 2 AIV subtypes, we report LBM-level, bird-level, and environmental swab specimen-level prevalences of A(H5) and A(H9) viruses on the basis of the best H9 models with LBM type (Tables 1, 2). This reporting did not affect interpretation of results, and we provide estimates obtained with more parsimonious models (online Technical Appendix 2 Tables 3–6).

LBM-level A(H5) virus prevalence was lower in retail LBMs than in mixed LBMs, and the posterior median estimate was ≈100% for mixed LBMs. However, among contaminated LBMs, levels of virus detection in birds and environmental areas did not vary between LBM types, but A(H5) virus prevalence in waterfowl was ≈6 times higher than in chickens (Figure). The prevalence did not vary between chicken breeds or environmental areas.

In contrast to that for A(H5) virus, we found that the posterior median estimate of the LBM-level A(H9) virus prevalence was ≈100% for retail and mixed LBM groups, but the level of virus detection in birds and environmental areas was higher for mixed LBMs than for retail LBMs. A(H9) virus prevalence was highest in broilers and lowest in waterfowl. The prevalence in broilers was 3.8 times as
<table><caption>**Table 1.** Prevalence of avian influenza A(H5) and A(H9) viruses in pooled poultry and environmental samples, Chittagong and Dhaka, Bangladesh*</caption><thead><tr><th rowspan="2">Sample type</th><th rowspan="2">No. pools</th><th colspan="2">A(H5) virus prevalence, %</th><th colspan="2">A(H9) virus prevalence, %</th></tr><tr><th>Observed pool level</th><th>Estimated bird level (95% HDI)†</th><th>Observed pool level</th><th>Estimated bird level (95% HDI)†</th></tr></thead><tbody><tr><td colspan="6"><strong>Poultry</strong></td></tr><tr><td>Retail LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Broiler</td><td>96</td><td>5.2</td><td>0.9 (0–9.4)</td><td>37.5</td><td>10.8 (2.3–22.5)</td></tr><tr><td>Sonali</td><td>62</td><td>3.2</td><td>1.4 (0–13.2)</td><td>32.3</td><td>6.6 (1.2–14.5)</td></tr><tr><td>Desi</td><td>61</td><td>3.4</td><td>1.3 (0–12.2)</td><td>29.5</td><td>6.8 (1.3–14.9)</td></tr><tr><td>Waterfowl</td><td>20</td><td>50.0</td><td>8.1 (0–46.8)</td><td>25.0</td><td>2.8 (0.3–7.0)</td></tr><tr><td>Mixed LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Broiler</td><td>57</td><td>1.8</td><td>0.9 (0–4.0)</td><td>47.4</td><td>13.1 (1.2–30.1)</td></tr><tr><td>Sonali</td><td>60</td><td>10.0</td><td>1.4 (0–5.7)</td><td>31.7</td><td>8.0 (0.4–19.8)</td></tr><tr><td>Desi</td><td>66</td><td>9.1</td><td>1.3 (0–5.2)</td><td>39.4</td><td>8.3 (0.5–20.4)</td></tr><tr><td>Waterfowl</td><td>55</td><td>27.3</td><td>7.6 (0–24.6)</td><td>16.4</td><td>3.4 (0.1–9.7)</td></tr><tr><td colspan="6"><strong>Environmental site</strong></td></tr><tr><td>Retail LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Stall area</td><td>101</td><td>5.9</td><td>1.5 (0–10.4)</td><td>16.8</td><td>3.2 (0.1–9.1)</td></tr><tr><td>Slaughter area</td><td>99</td><td>7.1</td><td>1.4 (0–10.1)</td><td>25.3</td><td>6.2 (0.2–16.6)</td></tr><tr><td>Mixed LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Stall area</td><td>102</td><td>15.7</td><td>3.1 (0–11.3)</td><td>23.5</td><td>5.2 (0.1–14.1)</td></tr><tr><td>Slaughter area</td><td>98</td><td>14.3</td><td>3.0 (0–11.0)</td><td>37.8</td><td>9.9 (0.4–25.0)</td></tr></tbody></table>

high as that in waterfowl and 1.6 times as high as that in Desi and Sonali (Figure). The environmental swab specimen-level prevalence was ≈2 times as high for slaughter areas than for stall areas (Figure).

Discussion

We detected A(H5) and A(H9) viruses in marketed poultry and environmental sites in nearly all LBMs sampled in Chittagong and Dhaka. The prevalence of A(H5) virus was higher in waterfowl than in chickens, whereas the prevalence of A(H9) virus was higher in chickens than waterfowl and also varied among chicken types, being more prevalent in broilers than in Desi and Sonali breeds. Slaughter areas were more frequently contaminated by A(H9) virus than stall areas. Whereas mixed LBMs were more frequently contaminated by A(H5) virus than were retail LBMs, prevalence of A(H9) virus was higher in mixed LBMs than in retail LBMs for birds and environmental areas.

AIVs were ubiquitous in surveyed LBMs. The LBM-level prevalence of A(H5) virus in Bangladesh was higher than in other AIV-endemic countries, including Egypt (32) and Vietnam (16). For both AIV subtypes, LBM-level prevalence was also higher than in another study conducted

in Chittagong (21), which found that 17.5% of LBMs had ≥1 environmental sample pool contaminated by A(H5) virus and 12.5% of LBMs had ≥1 environmental sample pool contaminated by A(H9) virus. This difference might have been caused by different sampling schemes; in our study, we collected a larger number of pools per LBM.

Bird-level prevalence was also higher than that reported in other AIV-endemic countries, including Bangladesh (1,4,16,19). However, care must be taken when comparing these results because studies used different study designs and sample screening protocols over different periods. Bird-level prevalence for contaminated LBMs was much higher than for virologic surveys conducted in backyard and commercial farms in Bangladesh (1,4,33,34). This finding suggests that virus transmission was amplified along the value chain from farms to LBMs. Overcrowding and continuous supply of susceptible birds of different species and breeds might have created conditions promoting the silent transmission of AIVs within these markets (10).

Our results suggest that birds in LBMs with a mixture of wholesalers and retailers were at higher risk for infection than birds in LBMs with primarily retail poultry businesses. Poultry value chains supplying different business types might differ structurally, thereby affecting the risk for

<table><caption>**Table 2.** Prevalence of avian influenza A(H5) and A(H9) viruses in LBMs, Chittagong and Dhaka, Bangladesh*</caption><thead><tr><th>Sample type</th><th>No. LBMs</th><th>H5 virus median prevalence, % (95% HDI)</th><th>H9 virus median prevalence, % (95% HDI)</th></tr></thead><tbody><tr><td colspan="4"><strong>Poultry</strong></td></tr><tr><td>Retail LBM</td><td>20</td><td>69.9 (40.2–100.0)</td><td>96.4 (85.5–100.0)</td></tr><tr><td>Mixed LBM</td><td>20</td><td>92.0 (72.3–100.0)</td><td>96.0 (84.0–100.0)</td></tr><tr><td colspan="4"><strong>Environmental sites</strong></td></tr><tr><td>Retail LBM</td><td>20</td><td>76.5 (47.2–100.0)</td><td>94.9 (80.5–100.0)</td></tr><tr><td>Mixed LBM</td><td>20</td><td>93.2 (75.5–100.0)</td><td>96.0 (84.0–100.0)</td></tr></tbody></table>

<footnote>*Desi, which means "local" in Bengali, are indigenous chicken breeds raised in backyard farms. Sonali is a cross-breed of the Rhode Island Red cocks and Fayoumi hens. HDI, high-density interval; LBM, live bird market.</footnote>

<footnote>†Bird and environmental swab specimen-level prevalence in contaminated live bird markets. Median values are reported.</footnote>

<footnote>*Prevalence estimates were made by using the best H9 models with LBM type. HDI, high-density interval; LBM, live bird market.</footnote>
<table><caption>**Table 1.** Prevalence of avian influenza A(H5) and A(H9) viruses in pooled poultry and environmental samples, Chittagong and Dhaka, Bangladesh*</caption><thead><tr><th rowspan="2">Sample type</th><th rowspan="2">No. pools</th><th colspan="2">A(H5) virus prevalence, %</th><th colspan="2">A(H9) virus prevalence, %</th></tr><tr><th>Observed pool level</th><th>Estimated bird level (95% HDI)†</th><th>Observed pool level</th><th>Estimated bird level (95% HDI)†</th></tr></thead><tbody><tr><td colspan="6"><strong>Poultry</strong></td></tr><tr><td>Retail LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Broiler</td><td>96</td><td>5.2</td><td>0.9 (0–9.4)</td><td>37.5</td><td>10.8 (2.3–22.5)</td></tr><tr><td>Sonali</td><td>62</td><td>3.2</td><td>1.4 (0–13.2)</td><td>32.3</td><td>6.6 (1.2–14.5)</td></tr><tr><td>Desi</td><td>61</td><td>3.4</td><td>1.3 (0–12.2)</td><td>29.5</td><td>6.8 (1.3–14.9)</td></tr><tr><td>Waterfowl</td><td>20</td><td>50.0</td><td>8.1 (0–46.8)</td><td>25.0</td><td>2.8 (0.3–7.0)</td></tr><tr><td>Mixed LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Broiler</td><td>57</td><td>1.8</td><td>0.9 (0–4.0)</td><td>47.4</td><td>13.1 (1.2–30.1)</td></tr><tr><td>Sonali</td><td>60</td><td>10.0</td><td>1.4 (0–5.7)</td><td>31.7</td><td>8.0 (0.4–19.8)</td></tr><tr><td>Desi</td><td>66</td><td>9.1</td><td>1.3 (0–5.2)</td><td>39.4</td><td>8.3 (0.5–20.4)</td></tr><tr><td>Waterfowl</td><td>55</td><td>27.3</td><td>7.6 (0–24.6)</td><td>16.4</td><td>3.4 (0.1–9.7)</td></tr><tr><td colspan="6"><strong>Environmental site</strong></td></tr><tr><td>Retail LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Stall area</td><td>101</td><td>5.9</td><td>1.5 (0–10.4)</td><td>16.8</td><td>3.2 (0.1–9.1)</td></tr><tr><td>Slaughter area</td><td>99</td><td>7.1</td><td>1.4 (0–10.1)</td><td>25.3</td><td>6.2 (0.2–16.6)</td></tr><tr><td>Mixed LBM</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Stall area</td><td>102</td><td>15.7</td><td>3.1 (0–11.3)</td><td>23.5</td><td>5.2 (0.1–14.1)</td></tr><tr><td>Slaughter area</td><td>98</td><td>14.3</td><td>3.0 (0–11.0)</td><td>37.8</td><td>9.9 (0.4–25.0)</td></tr></tbody></table>

high as that in waterfowl and 1.6 times as high as that in Desi and Sonali (Figure). The environmental swab specimen-level prevalence was ≈2 times as high for slaughter areas than for stall areas (Figure).

Discussion

We detected A(H5) and A(H9) viruses in marketed poultry and environmental sites in nearly all LBMs sampled in Chittagong and Dhaka. The prevalence of A(H5) virus was higher in waterfowl than in chickens, whereas the prevalence of A(H9) virus was higher in chickens than waterfowl and also varied among chicken types, being more prevalent in broilers than in Desi and Sonali breeds. Slaughter areas were more frequently contaminated by A(H9) virus than stall areas. Whereas mixed LBMs were more frequently contaminated by A(H5) virus than were retail LBMs, prevalence of A(H9) virus was higher in mixed LBMs than in retail LBMs for birds and environmental areas.

AIVs were ubiquitous in surveyed LBMs. The LBM-level prevalence of A(H5) virus in Bangladesh was higher than in other AIV-endemic countries, including Egypt (32) and Vietnam (16). For both AIV subtypes, LBM-level prevalence was also higher than in another study conducted

in Chittagong (21), which found that 17.5% of LBMs had ≥1 environmental sample pool contaminated by A(H5) virus and 12.5% of LBMs had ≥1 environmental sample pool contaminated by A(H9) virus. This difference might have been caused by different sampling schemes; in our study, we collected a larger number of pools per LBM.

Bird-level prevalence was also higher than that reported in other AIV-endemic countries, including Bangladesh (1,4,16,19). However, care must be taken when comparing these results because studies used different study designs and sample screening protocols over different periods. Bird-level prevalence for contaminated LBMs was much higher than for virologic surveys conducted in backyard and commercial farms in Bangladesh (1,4,33,34). This finding suggests that virus transmission was amplified along the value chain from farms to LBMs. Overcrowding and continuous supply of susceptible birds of different species and breeds might have created conditions promoting the silent transmission of AIVs within these markets (10).

Our results suggest that birds in LBMs with a mixture of wholesalers and retailers were at higher risk for infection than birds in LBMs with primarily retail poultry businesses. Poultry value chains supplying different business types might differ structurally, thereby affecting the risk for

<table><caption>**Table 2.** Prevalence of avian influenza A(H5) and A(H9) viruses in LBMs, Chittagong and Dhaka, Bangladesh*</caption><thead><tr><th>Sample type</th><th>No. LBMs</th><th>H5 virus median prevalence, % (95% HDI)</th><th>H9 virus median prevalence, % (95% HDI)</th></tr></thead><tbody><tr><td colspan="4"><strong>Poultry</strong></td></tr><tr><td>Retail LBM</td><td>20</td><td>69.9 (40.2–100.0)</td><td>96.4 (85.5–100.0)</td></tr><tr><td>Mixed LBM</td><td>20</td><td>92.0 (72.3–100.0)</td><td>96.0 (84.0–100.0)</td></tr><tr><td colspan="4"><strong>Environmental sites</strong></td></tr><tr><td>Retail LBM</td><td>20</td><td>76.5 (47.2–100.0)</td><td>94.9 (80.5–100.0)</td></tr><tr><td>Mixed LBM</td><td>20</td><td>93.2 (75.5–100.0)</td><td>96.0 (84.0–100.0)</td></tr></tbody></table>

<footnote>*Desi, which means "local" in Bengali, are indigenous chicken breeds raised in backyard farms. Sonali is a cross-breed of the Rhode Island Red cocks and Fayoumi hens. HDI, high-density interval; LBM, live bird market.</footnote>

<footnote>†Bird and environmental swab specimen-level prevalence in contaminated live bird markets. Median values are reported.</footnote>

<footnote>*Prevalence estimates were made by using the best H9 models with LBM type. HDI, high-density interval; LBM, live bird market.</footnote>
<figure><img src="image_1.png" /></figure>

**Figure.** Bird and environmental swab specimen-level avian influenza A(H5) and A(H9) virus prevalence ratios, Bangladesh. Dotted lines indicate H5 subtypes, and solid lines indicate H9 subtypes. Diamonds indicate median values, and horizontal bars indicate 95% high-density interval of a given prevalence ratio. Asterisks (*) indicate reference groups for each comparison. Desi, which means "local" in Bengali, are indigenous chicken breeds raised in backyard farms. Sonali is a cross-breed of the Rhode Island Red cocks and Fayoumi hens.

introduced birds being already infected. Wholesalers gen-
erally trade a larger number of birds from more diverse geo-
graphic origins than do retailers (23) and therefore might
have increased likelihood of virus introduction into mixed
LBMs. Moreover, because wholesalers might sell birds to
retailers in the same LBM (23), virus amplification might
be increased through the presence of wholesalers.

The higher prevalence of A(H9) virus in broilers than
in Sonali and Desi might result from differences in the struc-
ture of their respective value chains (23). Depending on the
chicken type, different value chain actors might be involved
and their trading practices might differ (23). The amount
of time chickens spend with traders, the density at which
chickens are kept in flocks of traders, and the frequency of
contact with chickens from other flocks might vary with
chicken type. The greater number of broilers marketed in
surveyed LBMs might mean that broilers are more likely
than Desi and Sonali to be sourced from large numbers of
flocks, which are then mixed in densely populated trucks
during transport to LBMs, promoting AIV transmission.

However, these prevalence patterns might also be caused
by varying levels of genetic susceptibility to AIV infection
(35,36). Further investigations are needed to disentangle
the possible influences of trade-related and genetic factors
on AIV transmission in these chicken types. The higher
level of contamination with A(H9) virus in slaughter areas
than in stall areas suggests that, in the absence of appropri-
ate biosecurity measures, slaughtering is likely to expose
humans to AIVs by fomite transmission (37).

Co-circulation of A(H5) and A(H9) viruses arouses
concerns over evolution of novel reassortant variants (5–
8). Detection of both subtypes in some poultry pools sug-
gests that these subtypes co-circulated near each other or
in the same host during the study period. Although A(H5)
viruses have considerable variability in their ability to in-
fect, cause disease, and be transmitted among waterfowl
(38), waterfowl are generally known to be less suscep-
tible to highly pathogenic avian influenza A(H5N1) vi-
ruses (39). Therefore, waterfowl could harbor this virus
but remain asymptomatic and serve as a potential host
<figure><img src="image_3.png" /></figure>

**Figure.** Bird and environmental swab specimen-level avian influenza A(H5) and A(H9) virus prevalence ratios, Bangladesh. Dotted lines indicate H5 subtypes, and solid lines indicate H9 subtypes. Diamonds indicate median values, and horizontal bars indicate 95% high-density interval of a given prevalence ratio. Asterisks (*) indicate reference groups for each comparison. Desi, which means "local" in Bengali, are indigenous chicken breeds raised in backyard farms. Sonali is a cross-breed of the Rhode Island Red cocks and Fayoumi hens.

introduced birds being already infected. Wholesalers gen-
erally trade a larger number of birds from more diverse geo-
graphic origins than do retailers (23) and therefore might
have increased likelihood of virus introduction into mixed
LBMs. Moreover, because wholesalers might sell birds to
retailers in the same LBM (23), virus amplification might
be increased through the presence of wholesalers.

The higher prevalence of A(H9) virus in broilers than
in Sonali and Desi might result from differences in the struc-
ture of their respective value chains (23). Depending on the
chicken type, different value chain actors might be involved
and their trading practices might differ (23). The amount
of time chickens spend with traders, the density at which
chickens are kept in flocks of traders, and the frequency of
contact with chickens from other flocks might vary with
chicken type. The greater number of broilers marketed in
surveyed LBMs might mean that broilers are more likely
than Desi and Sonali to be sourced from large numbers of
flocks, which are then mixed in densely populated trucks
during transport to LBMs, promoting AIV transmission.

However, these prevalence patterns might also be caused
by varying levels of genetic susceptibility to AIV infection
(35,36). Further investigations are needed to disentangle
the possible influences of trade-related and genetic factors
on AIV transmission in these chicken types. The higher
level of contamination with A(H9) virus in slaughter areas
than in stall areas suggests that, in the absence of appropri-
ate biosecurity measures, slaughtering is likely to expose
humans to AIVs by fomite transmission (37).

Co-circulation of A(H5) and A(H9) viruses arouses
concerns over evolution of novel reassortant variants (5–
8). Detection of both subtypes in some poultry pools sug-
gests that these subtypes co-circulated near each other or
in the same host during the study period. Although A(H5)
viruses have considerable variability in their ability to in-
fect, cause disease, and be transmitted among waterfowl
(38), waterfowl are generally known to be less suscep-
tible to highly pathogenic avian influenza A(H5N1) vi-
ruses (39). Therefore, waterfowl could harbor this virus
but remain asymptomatic and serve as a potential host
for genesis of novel AIVs in the presence of other virus subtypes. Also, the high level of A(H9) virus circulation among chickens could provide an ideal environment for virus diversification and selection in the LBM system. The different prevalence patterns in chickens and water-fowl observed suggest that these poultry species should be separated in LBMs and that active surveillance of novel reassortant variants should be implemented.

This study had some limitations. First, our models only accounted for clustering of sampled birds at the LBM level, but not at stall level. It is plausible that clustering of samples at stall level has less influence on AIV infection probability across the study population than clustering at the LBM level because stallholders in a given LBM in Bangladesh are likely to be supplied by the same traders and trade between each another (23). However, potential risk factors at stall level, such as ducks and hygiene level (21), might cause heterogeneous levels of AIV infection across stalls.

Second, our models did not account for the fact that sampling units in each stratum were selected with unequal probabilities. Although we selected different numbers of birds for each poultry type to account for variations in poultry populations, birds were still selected with different probabilities because their populations varied between clusters and strata. This selection might have resulted in larger SEs and thus less precise estimates compared with what could have been obtained with proportional sample sizes. Moreover, the overall prevalence might have been biased toward prevalence in samples selected with higher probabilities.

Third, our models assumed perfect sensitivity and specificity of real-time RT-PCR for pooled samples. The assays used in this study are considered highly sensitive and specific (27,28), and previous studies did not report any differences in virus detection for pooled and individual samples (40–42). Furthermore, virus detection in our study was based on parallel interpretation of cloacal or oropharyngeal sample test results (i.e., positive if $≥1$ was positive). However, pools that were negative for the AIV matrix gene but positive for any of the H5 and H9 subtypes indicate that accounting for actual test sensitivity and specificity would enable more robust prevalence estimation. Virus isolation might be attempted for RT-PCR-positive pools to assess the viability of virus material. However, this testing was not attempted in our study. Each pool consisted of swab specimens from different birds or environmental sites. Thus, multiple AIV subtypes and virus species, including Newcastle disease viruses, could be present in the same pool and interfere with growth of each virus in chicken eggs (43). Should such studies be replicated, the collection of individual swab specimens and their pooling at the laboratory is recommended to enable analysis of individual swab specimens that formed a virus-positive pool.

Fourth, we collected samples over a short period to reduce variability that could arise from seasonal variations in AIV prevalence. We focused on winter months, which are often reported to be periods of higher risk for AIV infection (44). Therefore, our estimates only represented AIV prevalence during that period and did not capture seasonal changes.

Contrary to previous cross-sectional studies, our approach enabled us to estimate AIV prevalence not only by poultry species but also by chicken type and account for the type of LBMs in which sampled poultry were marketed. Despite most AIV surveys and surveillance activities being based on multistage sampling, single-level analytic methods are generally used to analyze their results, while ignoring within-market correlation in poultry infection status. Accounting for this effect by incorporating LBM-specific random effects in a hierarchical model, and enabling mutual influence between bird-level, environmental swab specimens-level, and LBM-level parameters, improved the reliability of prevalence estimates (29). When applied to other settings, this approach needs to be adapted on the basis of an understanding of the variety of poultry value chains. Information about LBM locations and about trading practices and numbers and types of poultry sold within these LBMs is rarely readily available and would need to be collected to inform the study design.

In conclusion, LBMs surveyed in Bangladesh were highly contaminated by A(H5) and A(H9) viruses. The level of virus detection was associated with the type of poultry and environmental area and the presence of wholesalers in LBMs. These findings need to be included in the design of risk-based surveillance and control interventions aimed at reducing AIV prevalence, human exposure, and the risk for emergence of novel virus reassortant variants.

## Acknowledgments

We thank Eric Brum for his support during study implementation and the participants involved in the study.

This study was supported by the BALZAC research program "Behavioural adaptations in live poultry trading and farming systems and zoonoses control in Bangladesh" (BB/L018993/1) and is 1 of 11 programs supported by the Zoonoses and Emerging Livestock Systems, a joint research initiative between the Biotechnology and Biological Sciences Research Council, the Defence Science and Technology Laboratory, the Department for International Development, the Economic and Social Sciences Research Council, the Medical Research Council, and the Natural Environment Research Council.

## About the Author

Mr. Kim is a doctoral student at the College of Veterinary Medicine and Life Sciences, City University of Hong Kong,
for genesis of novel AIVs in the presence of other virus subtypes. Also, the high level of A(H9) virus circulation among chickens could provide an ideal environment for virus diversification and selection in the LBM system. The different prevalence patterns in chickens and water-fowl observed suggest that these poultry species should be separated in LBMs and that active surveillance of novel reassortant variants should be implemented.

This study had some limitations. First, our models only accounted for clustering of sampled birds at the LBM level, but not at stall level. It is plausible that clustering of samples at stall level has less influence on AIV infection probability across the study population than clustering at the LBM level because stallholders in a given LBM in Bangladesh are likely to be supplied by the same traders and trade between each another (23). However, potential risk factors at stall level, such as ducks and hygiene level (21), might cause heterogeneous levels of AIV infection across stalls.

Second, our models did not account for the fact that sampling units in each stratum were selected with unequal probabilities. Although we selected different numbers of birds for each poultry type to account for variations in poultry populations, birds were still selected with different probabilities because their populations varied between clusters and strata. This selection might have resulted in larger SEs and thus less precise estimates compared with what could have been obtained with proportional sample sizes. Moreover, the overall prevalence might have been biased toward prevalence in samples selected with higher probabilities.

Third, our models assumed perfect sensitivity and specificity of real-time RT-PCR for pooled samples. The assays used in this study are considered highly sensitive and specific (27,28), and previous studies did not report any differences in virus detection for pooled and individual samples (40–42). Furthermore, virus detection in our study was based on parallel interpretation of cloacal or oropharyngeal sample test results (i.e., positive if $≥1$ was positive). However, pools that were negative for the AIV matrix gene but positive for any of the H5 and H9 subtypes indicate that accounting for actual test sensitivity and specificity would enable more robust prevalence estimation. Virus isolation might be attempted for RT-PCR-positive pools to assess the viability of virus material. However, this testing was not attempted in our study. Each pool consisted of swab specimens from different birds or environmental sites. Thus, multiple AIV subtypes and virus species, including Newcastle disease viruses, could be present in the same pool and interfere with growth of each virus in chicken eggs (43). Should such studies be replicated, the collection of individual swab specimens and their pooling at the laboratory is recommended to enable analysis of individual swab specimens that formed a virus-positive pool.

Fourth, we collected samples over a short period to reduce variability that could arise from seasonal variations in AIV prevalence. We focused on winter months, which are often reported to be periods of higher risk for AIV infection (44). Therefore, our estimates only represented AIV prevalence during that period and did not capture seasonal changes.

Contrary to previous cross-sectional studies, our approach enabled us to estimate AIV prevalence not only by poultry species but also by chicken type and account for the type of LBMs in which sampled poultry were marketed. Despite most AIV surveys and surveillance activities being based on multistage sampling, single-level analytic methods are generally used to analyze their results, while ignoring within-market correlation in poultry infection status. Accounting for this effect by incorporating LBM-specific random effects in a hierarchical model, and enabling mutual influence between bird-level, environmental swab specimens-level, and LBM-level parameters, improved the reliability of prevalence estimates (29). When applied to other settings, this approach needs to be adapted on the basis of an understanding of the variety of poultry value chains. Information about LBM locations and about trading practices and numbers and types of poultry sold within these LBMs is rarely readily available and would need to be collected to inform the study design.

In conclusion, LBMs surveyed in Bangladesh were highly contaminated by A(H5) and A(H9) viruses. The level of virus detection was associated with the type of poultry and environmental area and the presence of wholesalers in LBMs. These findings need to be included in the design of risk-based surveillance and control interventions aimed at reducing AIV prevalence, human exposure, and the risk for emergence of novel virus reassortant variants.

## Acknowledgments

We thank Eric Brum for his support during study implementation and the participants involved in the study.

This study was supported by the BALZAC research program "Behavioural adaptations in live poultry trading and farming systems and zoonoses control in Bangladesh" (BB/L018993/1) and is 1 of 11 programs supported by the Zoonoses and Emerging Livestock Systems, a joint research initiative between the Biotechnology and Biological Sciences Research Council, the Defence Science and Technology Laboratory, the Department for International Development, the Economic and Social Sciences Research Council, the Medical Research Council, and the Natural Environment Research Council.

## About the Author

Mr. Kim is a doctoral student at the College of Veterinary Medicine and Life Sciences, City University of Hong Kong,
Hong Kong, China. His primary research interests are the socioeconomic, cultural, and epidemiologic factors that shape zoonotic and animal infectious disease transmission, including avian influenza.

## References

1. Negovetich NJ, Feeroz MM, Jones-Engel L, Walker D, Alam SM, Hasan K, et al. Live bird markets of Bangladesh: H9N2 viruses and the near absence of highly pathogenic H5N1 influenza. PLoS One. 2011;6:e19311. http://dx.doi.org/10.1371/journal.pone.0019311

2. Chakraborty A. Outbreak of mild respiratory disease caused by H5N1 and H9N2 infections among young children in Dhaka, Bangladesh, 2011. Health Science Bulletin. 2011;9:5-12.

3. World Health Organization. Cumulative number of confirmed human cases for avian influenza A(H5N1) reported to WHO, 2003–2017; 2017 [cited 2017 Oct 9]. http://www.who.int/influenza/human_animal_interface/H5N1_cumulative_table_archives/en/

4. Turner JC, Feeroz MM, Hasan MK, Akhtar S, Walker D, Seiler P, et al. Insight into live bird markets of Bangladesh: an overview of the dynamics of transmission of H5N1 and H9N2 avian influenza viruses. Emer Microbes Infect. 2017;6:e12. http://dx.doi.org/10.1038/emi.2016.142

5. Lin YP, Shaw M, Gregory V, Cameron K, Lim W, Klimov A, et al. Avian-to-human transmission of H9N2 subtype influenza A viruses: relationship between H9N2 and H5N1 human isolates. Proc Natl Acad Sci U S A. 2000;97:9654–8. http://dx.doi.org/10.1073/pnas.160270697

6. Monne I, Yamage M, Dauphin G, Claes F, Ahmed G, Giasuddin M, et al. Reassortant avian influenza A(H5N1) viruses with H9N2-PB1 gene in poultry, Bangladesh. Emer Infect Dis. 2013;19:1630–4. http://dx.doi.org/10.3201/eid1910.130534

7. Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, et al. Human infection with a novel avian-origin influenza A (H7N9) virus. N Engl J Med. 2013;368:1888–97. http://dx.doi.org/10.1056/NEJMoa1304459

8. Chen H, Yuan H, Gao R, Zhang J, Wang D, Xiong Y, et al. Clinical and epidemiological characteristics of a fatal case of avian influenza A H10N8 virus infection: a descriptive study. Lancet. 2014;383:714–21. http://dx.doi.org/10.1016/S0140-6736 (14)60111-2

9. Lee DH, Bertran K, Kwon JH, Swayne DE. Evolution, global spread, and pathogenicity of highly pathogenic avian influenza H5Nx clade 2.3.4.4. J Vet Sci. 2017;18(S1):269–80. http://dx.doi.org/10.4142/jvs.2017.18.S1.269

10. Webster RG. Wet markets—a continuing source of severe acute respiratory syndrome and influenza? Lancet. 2004;363:234–6. http://dx.doi.org/10.1016/S0140-6736(03)15329-9

11. Nguyen DC, Uyeki TM, Jadhao S, Maines T, Shaw M, Matsuoka Y, 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–12. http://dx.doi.org/10.1128/JVI.79.7.4201-4212.2005

12. Fournié G, Guitian FJ, Mangtani P, Ghani AC. Impact of the implementation of rest days in live bird markets on the dynamics of H5N1 highly pathogenic avian influenza. J R Soc Interface. 2011;8:1079–89. http://dx.doi.org/10.1098/rsif.2010.0510

13. Liu M, He S, Walker D, Zhou N, Perez DR, Mo B, et al. The influenza virus gene pool in a poultry market in South central China. Virology. 2003;305:267–75. http://dx.doi.org/10.1006/viro.2002.1762

14. Okamatsu M, Nishi T, Nomura N, Yamamoto N, Sakoda Y, Sakurai K, et al. The genetic and antigenic diversity of avian influenza viruses isolated from domestic ducks, muscovy ducks, and chickens in northern and southern Vietnam, 2010–2012.

Virus Genes. 2013;47:317–29. http://dx.doi.org/10.1007/s11262-013-0954-7

15. Phan MQ, Henry W, Bui CB, Do DH, Hoang NV, Thu NT, et al. Detection of HPAI H5N1 viruses in ducks sampled from live bird markets in Vietnam. Epidemiol Infect. 2013;141:601–11. http://dx.doi.org/10.1017/S0950268812001112

16. Nguyen DT, Bryant JE, Davis CT, Nguyen LV, Pham LT, Loth L, et al. Prevalence and distribution of avian influenza A(H5N1) virus clade variants in live bird markets of Vietnam, 2011–2013. Avian Dis. 2014;58:599–608. http://dx.doi.org/10.1637/10814-030814-Reg

17. Huang Y, Zhang H, Li X, Hu S, Cai L, Sun Q, et al. Detection and genetic characteristics of H9N2 avian influenza viruses from live poultry markets in Hunan Province, China. PLoS One. 2015;10:e0142584. http://dx.doi.org/10.1371/journal.pone.0142584

18. Chen LJ, Lin XD, Guo WP, Tian JH, Wang W, Ying XH, et al. Diversity and evolution of avian influenza viruses in live poultry markets, free-range poultry and wild wetland birds in China. J Gen Virol. 2016;97:844–54. http://dx.doi.org/10.1099/jgv.0.000399

19. Thuy DM, Peacock TP, Bich VTN, Fabrizio T, Hoang DN, Tho ND, et al. Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014. Infect Genet Evol. 2016;44:530–40. http://dx.doi.org/10.1016/j.meegid.2016.06.038

20. Wang X, Wang Q, Cheng W, Yu Z, Ling F, Mao H, et al. Risk factors for avian influenza virus contamination of live poultry markets in Zhejiang, China during the 2015–2016 human influenza season. Sci Rep. 2017;7:42722. http://dx.doi.org/10.1038/srep42722

21. Sayeed MA, Smallwood C, Imam T, Mahmud R, Hasan RB, Hasan M, et al. Assessment of hygienic conditions of live bird markets on avian influenza in Chittagong metro, Bangladesh. Prev Vet Med. 2017;142:7–15. http://dx.doi.org/10.1016/j.prevetmed.2017.04.009

22. Biswas PK, Biswas D, Ahmed S, Rahman A, Debnath NC. A longitudinal study of the incidence of major endemic and epidemic diseases affecting semi-scavenging chickens reared under the Participatory Livestock Development Project areas in Bangladesh. Avian Pathol. 2005;34:303–12. http://dx.doi.org/10.1080/03079450500178972

23. Moyen N, Ahmed G, Gupta S, Tenzin T, Khan R, Khan T, et al. A large-scale study of a poultry trading network in Bangladesh: implications for control and surveillance of avian influenza viruses. BMC Vet Res. 2018;14:12. http://dx.doi.org/10.1186/s12917-018-1331-5

24. Fournié G, Tripodi A, Nguyen TT, Nguyen VT, Tran TT, Bisson A, et al. Investigating poultry trade patterns to guide avian influenza surveillance and control: a case study in Vietnam. Sci Rep. 2016;6:29463. http://dx.doi.org/10.1038/srep29463

25. Molia S, Boly IA, Duboz R, Coulibaly B, Guitian J, Grosbois V, et al. Live bird markets characterization and trading network analysis in Mali: implications for the surveillance and control of avian influenza and Newcastle disease. Acta Trop. 2016;155:77–88. http://dx.doi.org/10.1016/j.acnatropica.2015.12.003

26. Indriani R, Samaan G, Gultom A, Loth L, Irianti S, Adjid R, et al. Environmental sampling for avian influenza virus A (H5N1) in live-bird markets, Indonesia. Emer Infect Dis. 2010;16:1889–95. http://dx.doi.org/10.3201/eid1612.100402

27. Monne I, Ormelli S, Salviato A, De Battisti C, Bettini F, Salomoni A, et al. Development and validation of a one-step real-time PCR assay for simultaneous detection of subtype H5, H7, and H9 avian influenza viruses. J Clin Microbiol. 2008;46:1769–73. http://dx.doi.org/10.1128/JCM.02204-07

28. Heine HG, Foord AJ, Wang J, Valdeter S, Walker S, Morrissy C, et al. Detection of highly pathogenic zoonotic influenza virus H5N6 by reverse-transcriptase quantitative polymerase chain reaction. Virol J. 2015;12:18. http://dx.doi.org/10.1186/s12985-015-0250-3
Hong Kong, China. His primary research interests are the socioeconomic, cultural, and epidemiologic factors that shape zoonotic and animal infectious disease transmission, including avian influenza.

## References

1. Negovetich NJ, Feeroz MM, Jones-Engel L, Walker D, Alam SM, Hasan K, et al. Live bird markets of Bangladesh: H9N2 viruses and the near absence of highly pathogenic H5N1 influenza. PLoS One. 2011;6:e19311. http://dx.doi.org/10.1371/journal.pone.0019311

2. Chakraborty A. Outbreak of mild respiratory disease caused by H5N1 and H9N2 infections among young children in Dhaka, Bangladesh, 2011. Health Science Bulletin. 2011;9:5-12.

3. World Health Organization. Cumulative number of confirmed human cases for avian influenza A(H5N1) reported to WHO, 2003–2017; 2017 [cited 2017 Oct 9]. http://www.who.int/influenza/human_animal_interface/H5N1_cumulative_table_archives/en/

4. Turner JC, Feeroz MM, Hasan MK, Akhtar S, Walker D, Seiler P, et al. Insight into live bird markets of Bangladesh: an overview of the dynamics of transmission of H5N1 and H9N2 avian influenza viruses. Emer Microbes Infect. 2017;6:e12. http://dx.doi.org/10.1038/emi.2016.142

5. Lin YP, Shaw M, Gregory V, Cameron K, Lim W, Klimov A, et al. Avian-to-human transmission of H9N2 subtype influenza A viruses: relationship between H9N2 and H5N1 human isolates. Proc Natl Acad Sci U S A. 2000;97:9654–8. http://dx.doi.org/10.1073/pnas.160270697

6. Monne I, Yamage M, Dauphin G, Claes F, Ahmed G, Giasuddin M, et al. Reassortant avian influenza A(H5N1) viruses with H9N2-PB1 gene in poultry, Bangladesh. Emer Infect Dis. 2013;19:1630–4. http://dx.doi.org/10.3201/eid1910.130534

7. Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, et al. Human infection with a novel avian-origin influenza A (H7N9) virus. N Engl J Med. 2013;368:1888–97. http://dx.doi.org/10.1056/NEJMoa1304459

8. Chen H, Yuan H, Gao R, Zhang J, Wang D, Xiong Y, et al. Clinical and epidemiological characteristics of a fatal case of avian influenza A H10N8 virus infection: a descriptive study. Lancet. 2014;383:714–21. http://dx.doi.org/10.1016/S0140-6736 (14)60111-2

9. Lee DH, Bertran K, Kwon JH, Swayne DE. Evolution, global spread, and pathogenicity of highly pathogenic avian influenza H5Nx clade 2.3.4.4. J Vet Sci. 2017;18(S1):269–80. http://dx.doi.org/10.4142/jvs.2017.18.S1.269

10. Webster RG. Wet markets—a continuing source of severe acute respiratory syndrome and influenza? Lancet. 2004;363:234–6. http://dx.doi.org/10.1016/S0140-6736(03)15329-9

11. Nguyen DC, Uyeki TM, Jadhao S, Maines T, Shaw M, Matsuoka Y, 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–12. http://dx.doi.org/10.1128/JVI.79.7.4201-4212.2005

12. Fournié G, Guitian FJ, Mangtani P, Ghani AC. Impact of the implementation of rest days in live bird markets on the dynamics of H5N1 highly pathogenic avian influenza. J R Soc Interface. 2011;8:1079–89. http://dx.doi.org/10.1098/rsif.2010.0510

13. Liu M, He S, Walker D, Zhou N, Perez DR, Mo B, et al. The influenza virus gene pool in a poultry market in South central China. Virology. 2003;305:267–75. http://dx.doi.org/10.1006/viro.2002.1762

14. Okamatsu M, Nishi T, Nomura N, Yamamoto N, Sakoda Y, Sakurai K, et al. The genetic and antigenic diversity of avian influenza viruses isolated from domestic ducks, muscovy ducks, and chickens in northern and southern Vietnam, 2010–2012.

Virus Genes. 2013;47:317–29. http://dx.doi.org/10.1007/s11262-013-0954-7

15. Phan MQ, Henry W, Bui CB, Do DH, Hoang NV, Thu NT, et al. Detection of HPAI H5N1 viruses in ducks sampled from live bird markets in Vietnam. Epidemiol Infect. 2013;141:601–11. http://dx.doi.org/10.1017/S0950268812001112

16. Nguyen DT, Bryant JE, Davis CT, Nguyen LV, Pham LT, Loth L, et al. Prevalence and distribution of avian influenza A(H5N1) virus clade variants in live bird markets of Vietnam, 2011–2013. Avian Dis. 2014;58:599–608. http://dx.doi.org/10.1637/10814-030814-Reg

17. Huang Y, Zhang H, Li X, Hu S, Cai L, Sun Q, et al. Detection and genetic characteristics of H9N2 avian influenza viruses from live poultry markets in Hunan Province, China. PLoS One. 2015;10:e0142584. http://dx.doi.org/10.1371/journal.pone.0142584

18. Chen LJ, Lin XD, Guo WP, Tian JH, Wang W, Ying XH, et al. Diversity and evolution of avian influenza viruses in live poultry markets, free-range poultry and wild wetland birds in China. J Gen Virol. 2016;97:844–54. http://dx.doi.org/10.1099/jgv.0.000399

19. Thuy DM, Peacock TP, Bich VTN, Fabrizio T, Hoang DN, Tho ND, et al. Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014. Infect Genet Evol. 2016;44:530–40. http://dx.doi.org/10.1016/j.meegid.2016.06.038

20. Wang X, Wang Q, Cheng W, Yu Z, Ling F, Mao H, et al. Risk factors for avian influenza virus contamination of live poultry markets in Zhejiang, China during the 2015–2016 human influenza season. Sci Rep. 2017;7:42722. http://dx.doi.org/10.1038/srep42722

21. Sayeed MA, Smallwood C, Imam T, Mahmud R, Hasan RB, Hasan M, et al. Assessment of hygienic conditions of live bird markets on avian influenza in Chittagong metro, Bangladesh. Prev Vet Med. 2017;142:7–15. http://dx.doi.org/10.1016/j.prevetmed.2017.04.009

22. Biswas PK, Biswas D, Ahmed S, Rahman A, Debnath NC. A longitudinal study of the incidence of major endemic and epidemic diseases affecting semi-scavenging chickens reared under the Participatory Livestock Development Project areas in Bangladesh. Avian Pathol. 2005;34:303–12. http://dx.doi.org/10.1080/03079450500178972

23. Moyen N, Ahmed G, Gupta S, Tenzin T, Khan R, Khan T, et al. A large-scale study of a poultry trading network in Bangladesh: implications for control and surveillance of avian influenza viruses. BMC Vet Res. 2018;14:12. http://dx.doi.org/10.1186/s12917-018-1331-5

24. Fournié G, Tripodi A, Nguyen TT, Nguyen VT, Tran TT, Bisson A, et al. Investigating poultry trade patterns to guide avian influenza surveillance and control: a case study in Vietnam. Sci Rep. 2016;6:29463. http://dx.doi.org/10.1038/srep29463

25. Molia S, Boly IA, Duboz R, Coulibaly B, Guitian J, Grosbois V, et al. Live bird markets characterization and trading network analysis in Mali: implications for the surveillance and control of avian influenza and Newcastle disease. Acta Trop. 2016;155:77–88. http://dx.doi.org/10.1016/j.acnatropica.2015.12.003

26. Indriani R, Samaan G, Gultom A, Loth L, Irianti S, Adjid R, et al. Environmental sampling for avian influenza virus A (H5N1) in live-bird markets, Indonesia. Emer Infect Dis. 2010;16:1889–95. http://dx.doi.org/10.3201/eid1612.100402

27. Monne I, Ormelli S, Salviato A, De Battisti C, Bettini F, Salomoni A, et al. Development and validation of a one-step real-time PCR assay for simultaneous detection of subtype H5, H7, and H9 avian influenza viruses. J Clin Microbiol. 2008;46:1769–73. http://dx.doi.org/10.1128/JCM.02204-07

28. Heine HG, Foord AJ, Wang J, Valdeter S, Walker S, Morrissy C, et al. Detection of highly pathogenic zoonotic influenza virus H5N6 by reverse-transcriptase quantitative polymerase chain reaction. Virol J. 2015;12:18. http://dx.doi.org/10.1186/s12985-015-0250-3
29. Kruschke JK. Doing Bayesian data analysis: a tutorial with R, JAGS, and Stan. 2nd ed. New York: Academic Press; 2014.

30. Plummer M. JAGS version 3.4.0 user manual, 2013 [cited 2018 Aug 23]. http://www.stats.ox.ac.uk/~nicholls/MScMCMC15/ jags_user_manual.pdf

31. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2016 [cited 2018 Aug 23]. https://www.r-project.org/

32. Abdelwhab EM, Selim AA, Arafa A, Galal S, Kilany WH, Hassan MK, et al. Circulation of avian influenza H5N1 in live bird markets in Egypt. Avian Dis. 2010;54:911-4. http://dx.doi.org/10.1637/9099-100809-RESNOTE.1

33. Khatun A, Giasuddin M, Islam KM, Khanom S, Samad MA, Islam MR, et al. Surveillance of avian influenza virus type A in semi-scavenging ducks in Bangladesh. BMC Vet Res. 2013;9:196. http://dx.doi.org/10.1186/1746-6148-9-196

34. Sarkar S, Khan SU, Mikolon A, Rahman MZ, Abedin J, Zeidner N, et al. An epidemiological study of avian influenza A (H5) virus in nomadic ducks and their raising practices in northeastern Bangladesh, 2011-2012. Influenza Other Respi Viruses. 2017;11:275-82. http://dx.doi.org/10.1111/irv.12438

35. Blohm U, Weigend S, Preisinger R, Beer M, Hoffmann D. Immunological competence of different domestic chicken breeds against avian influenza infection. Avian Dis. 2016;60(Suppl): 262-8. http://dx.doi.org/10.1637/11159-051615-RegR

36. Ruiz-Hernandez R, Mwangi W, Peroval M, Sadeyen JR, Ascough S, Balkissoon D, et al. Host genetics determine susceptibility to avian influenza infection and transmission dynamics. Sci Rep. 2016;6:26787. http://dx.doi.org/10.1038/srep26787

37. Fournié G, Høg E, Barnett T, Pfeiffer DU, Mangtani P. A systematic review and meta-analysis of practices exposing humans to avian influenza viruses, their prevalence, and rationale. Am J Trop Med Hyg. 2017;97:376-88. http://dx.doi.org/10.4269/ajtmh.17-0014

38. Hulse-Post DJ, Sturm-Ramirez KM, Humberd J, Seiler P, Govorkova EA, Krauss S, et al. Role of domestic ducks in the propagation and biological evolution of highly pathogenic H5N1 influenza viruses in Asia. Proc Natl Acad Sci U S A. 2005;102:10682-7. http://dx.doi.org/10.1073/pnas.0504662102

39. Kim JK, Negovetich NJ, Forrest HL, Webster RG. Ducks: the “Trojan horses” of H5N1 influenza. Influenza Other Respi Viruses. 2009;3:121-8. http://dx.doi.org/10.1111/j.1750-2659.2009.00084.x

40. Ladman BS, Spackman E, Gelb J Jr. Comparison of pooling 11 or 5 oropharyngeal swabbings for detecting avian influenza virus by real-time reverse transcription-PCR in broiler chickens. Avian Dis. 2012;56:227-9. http://dx.doi.org/10.1637/9839-062011-ResNote.1

41. Arnold ME, Slomka MJ, Coward VJ, Mahmood S, Raleigh PJ, Brown IH. Evaluation of the pooling of swabs for real-time PCR detection of low titre shedding of low pathogenicity avian influenza in turkeys. Epidemiol Infect. 2013;141:1286-97. http://dx.doi.org/10.1017/S0950268812001811

42. Spackman E, Pedersen JC, McKinley ET, Gelb J Jr. Optimal specimen collection and transport methods for the detection of avian influenza virus and Newcastle disease virus. BMC Vet Res. 2013;9:35. http://dx.doi.org/10.1186/1746-6148-9-35

43. Ge S, Zheng D, Zhao Y, Liu H, Liu W, Sun Q, et al. Evaluating viral interference between Influenza virus and Newcastle disease virus using real-time reverse transcription-polymerase chain reaction in chicken eggs. Virol J. 2012;9:128. http://dx.doi.org/10.1186/1743-422X-9-128

44. Park AW, Glass K. Dynamic patterns of avian and human influenza in east and southeast Asia. Lancet Infect Dis. 2007;7:543-8. http://dx.doi.org/10.1016/S1473-3099(07)70186-X

Address for correspondence: Younjung Kim, City University of Hong Kong, Rm 1B-410, 4/F, Block 1, To Yuen Bldg, 31 To Yuen St, Kowloon, Hong Kong, China; email: younjung.kim@my.cityu.edu.hk

## EID Podcast:
Emerging Infectious Diseases Cover Art

Byron Breedlove, managing editor of the journal, elaborates on aesthetic considerations and historical factors, as well as the complexities of obtaining artwork for Emerging Infectious Diseases.

EMERGING
INFECTIOUS DISEASES

Modeling

CDC

January 2017

<figure><img src="image_2.png" /></figure>

EMERGING
INFECTIOUS DISEASES

CDC

November 2017

<figure><img src="image_3.png" /></figure>

EMERGING
INFECTIOUS DISEASES

National Bacteriology Unit

CDC

August 2017

<figure><img src="image_4.png" /></figure>

EMERGING
INFECTIOUS DISEASES

Disease Patterns

CDC

October 2017

<figure><img src="image_5.png" /></figure>

EMERGING
INFECTIOUS DISEASES

Control Navigation System Infections

CDC

July 2017

<figure><img src="image_6.png" /></figure>

EMERGING
INFECTIOUS DISEASES

Vaccine-Neutral Dissemination

CDC

August 2014

<figure><img src="image_7.png" /></figure>

Visit our website to listen:

https://www2c.cdc.gov/podcasts/player.asp?f=8646224

EMERGING
INFECTIOUS DISEASES
29. Kruschke JK. Doing Bayesian data analysis: a tutorial with R, JAGS, and Stan. 2nd ed. New York: Academic Press; 2014.

30. Plummer M. JAGS version 3.4.0 user manual, 2013 [cited 2018 Aug 23]. http://www.stats.ox.ac.uk/~nicholls/MScMCMC15/ jags_user_manual.pdf

31. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2016 [cited 2018 Aug 23]. https://www.r-project.org/

32. Abdelwhab EM, Selim AA, Arafa A, Galal S, Kilany WH, Hassan MK, et al. Circulation of avian influenza H5N1 in live bird markets in Egypt. Avian Dis. 2010;54:911-4. http://dx.doi.org/10.1637/9099-100809-RESNOTE.1

33. Khatun A, Giasuddin M, Islam KM, Khanom S, Samad MA, Islam MR, et al. Surveillance of avian influenza virus type A in semi-scavenging ducks in Bangladesh. BMC Vet Res. 2013;9:196. http://dx.doi.org/10.1186/1746-6148-9-196

34. Sarkar S, Khan SU, Mikolon A, Rahman MZ, Abedin J, Zeidner N, et al. An epidemiological study of avian influenza A (H5) virus in nomadic ducks and their raising practices in northeastern Bangladesh, 2011-2012. Influenza Other Respi Viruses. 2017;11:275-82. http://dx.doi.org/10.1111/irv.12438

35. Blohm U, Weigend S, Preisinger R, Beer M, Hoffmann D. Immunological competence of different domestic chicken breeds against avian influenza infection. Avian Dis. 2016;60(Suppl): 262-8. http://dx.doi.org/10.1637/11159-051615-RegR

36. Ruiz-Hernandez R, Mwangi W, Peroval M, Sadeyen JR, Ascough S, Balkissoon D, et al. Host genetics determine susceptibility to avian influenza infection and transmission dynamics. Sci Rep. 2016;6:26787. http://dx.doi.org/10.1038/srep26787

37. Fournié G, Høg E, Barnett T, Pfeiffer DU, Mangtani P. A systematic review and meta-analysis of practices exposing humans to avian influenza viruses, their prevalence, and rationale. Am J Trop Med Hyg. 2017;97:376-88. http://dx.doi.org/10.4269/ajtmh.17-0014

38. Hulse-Post DJ, Sturm-Ramirez KM, Humberd J, Seiler P, Govorkova EA, Krauss S, et al. Role of domestic ducks in the propagation and biological evolution of highly pathogenic H5N1 influenza viruses in Asia. Proc Natl Acad Sci U S A. 2005;102:10682-7. http://dx.doi.org/10.1073/pnas.0504662102

39. Kim JK, Negovetich NJ, Forrest HL, Webster RG. Ducks: the “Trojan horses” of H5N1 influenza. Influenza Other Respi Viruses. 2009;3:121-8. http://dx.doi.org/10.1111/j.1750-2659.2009.00084.x

40. Ladman BS, Spackman E, Gelb J Jr. Comparison of pooling 11 or 5 oropharyngeal swabbings for detecting avian influenza virus by real-time reverse transcription-PCR in broiler chickens. Avian Dis. 2012;56:227-9. http://dx.doi.org/10.1637/9839-062011-ResNote.1

41. Arnold ME, Slomka MJ, Coward VJ, Mahmood S, Raleigh PJ, Brown IH. Evaluation of the pooling of swabs for real-time PCR detection of low titre shedding of low pathogenicity avian influenza in turkeys. Epidemiol Infect. 2013;141:1286-97. http://dx.doi.org/10.1017/S0950268812001811

42. Spackman E, Pedersen JC, McKinley ET, Gelb J Jr. Optimal specimen collection and transport methods for the detection of avian influenza virus and Newcastle disease virus. BMC Vet Res. 2013;9:35. http://dx.doi.org/10.1186/1746-6148-9-35

43. Ge S, Zheng D, Zhao Y, Liu H, Liu W, Sun Q, et al. Evaluating viral interference between Influenza virus and Newcastle disease virus using real-time reverse transcription-polymerase chain reaction in chicken eggs. Virol J. 2012;9:128. http://dx.doi.org/10.1186/1743-422X-9-128

44. Park AW, Glass K. Dynamic patterns of avian and human influenza in east and southeast Asia. Lancet Infect Dis. 2007;7:543-8. http://dx.doi.org/10.1016/S1473-3099(07)70186-X

Address for correspondence: Younjung Kim, City University of Hong Kong, Rm 1B-410, 4/F, Block 1, To Yuen Bldg, 31 To Yuen St, Kowloon, Hong Kong, China; email: younjung.kim@my.cityu.edu.hk

<figure><img src="image_4.png" /></figure>