Molecular epidemiology of influenza A (H5N1) viruses, Bangladesh, 2007-2011

Epidemiology and phylodynamics of multiple clades of H5N1 circulating in domestic duck farms in different production systems in Bangladesh

Waterfowl are considered to be natural reservoirs of the avian influenza virus (AIV). However, the dynamics of transmission and evolutionary patterns of AIV and its subtypes within duck farms in Bangladesh remain poorly documented. Hence, a cross-sectional study was conducted in nine districts of Bangladesh between 2019 and 2021, to determine the prevalence of AIV and its subtypes H5 and H9, as well as to identify risk factors and the phylodynamics of H5N1 clades circulating in domestic duck farms. The oropharyngeal and cloacal swab samples were tested for the AIV Matrix gene (M-gene) followed by H5, H7, and H9 subtypes using rRT-PCR. The exploratory analysis was performed to estimate AIV and its subtype prevalence in different production systems, and multivariable logistic regression model was used to identify the risk factors that influence AIV infection in ducks. Bayesian phylogenetic analysis was conducted to generate a maximum clade credibility (MCC) tree and the maximum likelihood method to determine the phylogenetic relationships of the H5N1 viruses isolated from ducks. AIV was detected in 40% (95% CI: 33.0-48.1) of the duck farms. The prevalence of AIV was highest in nomadic ducks (39.8%; 95% CI: 32.9-47.1), followed by commercial ducks (24.6%; 95% CI: 14.5-37.3) and backyard ducks (14.4%; 95% CI: 10.5-19.2). The H5 prevalence was also highest in nomadic ducks (19.4%; 95% CI: 14.0-25.7). The multivariable logistic regression model revealed that ducks from nomadic farms (AOR: 2.4; 95% CI: 1.45-3.93), juvenile (AOR: 2.2; 95% CI: 1.37-3.61), and sick ducks (AOR: 11.59; 95% CI: 4.82-32.44) had a higher risk of AIV. Similarly, the likelihood of H5 detection was higher in sick ducks (AOR: 40.8; 95% CI: 16.3-115.3). Bayesian phylogenetic analysis revealed that H5N1 viruses in ducks belong to two distinct clades, 2.3.2.1a, and 2.3.4.4b. The clade 2.3.2.1a (reassorted) has been evolving silently since 2015 and forming at least nine subgroups based on >90% posterior probability. Notably, clade 2.3.4.4b was introduced in ducks in Bangladesh by the end of the year 2020, which was genetically similar to viruses detected in wild birds in Japan, China, and Africa, indicating migration-associated transmission of an emerging panzootic clade. We recommend continuing AIV surveillance in the duck production system and preventing the intermingling of domestic ducks with migratory waterfowl in wetlands.

Annual trading patterns and risk factors of avian influenza A/H5 and A/H9 virus circulation in turkey birds (Meleagris gallopavo) at live bird markets in Dhaka city, Bangladesh

The impacts of the avian influenza virus (AIV) on farmed poultry and wild birds affect human health, livelihoods, food security, and international trade. The movement patterns of turkey birds from farms to live bird markets (LBMs) and infection of AIV are poorly understood in Bangladesh. Thus, we conducted weekly longitudinal surveillance in LBMs to understand the trading patterns, temporal trends, and risk factors of AIV circulation in turkey birds. We sampled a total of 423 turkeys from two LBMs in Dhaka between May 2018 and September 2019. We tested the swab samples for the AIV matrix gene (M-gene) followed by H5, H7, and H9 subtypes using real-time reverse transcriptase-polymerase chain reaction (rRT-PCR). We used exploratory analysis to investigate trading patterns, annual cyclic trends of AIV and its subtypes, and a generalized estimating equation (GEE) logistic model to determine the factors that influence the infection of H5 and H9 in turkeys. Furthermore, we conducted an observational study and informal interviews with traders and vendors to record turkey trading patterns, demand, and supply and turkey handling practices in LBM. We found that all trade routes of turkey birds to northern Dhaka are unidirectional and originate from the northwestern and southern regions of Bangladesh. The number of trades from the source district to Dhaka depends on the turkey density. The median distance that turkey was traded from its source district to Dhaka was 188 km (Q1 = 165, Q3 = 210, IQR = 45.5). We observed seasonal variation in the median and average distance of turkey. The qualitative findings revealed that turkey farming initially became reasonably profitable in 2018 and at the beginning of 2019. However, the fall in demand and production in the middle of 2019 may be related to unstable market pricing, high feed costs, a shortfall of adequate marketing facilities, poor consumer knowledge, and a lack of advertising. The overall prevalence of AIV, H5, and H9 subtypes in turkeys was 31% (95% CI: 26.6-35.4), 16.3% (95% CI: 12.8-19.8), and 10.2% (95% CI: 7.3-13.1) respectively. None of the samples were positive for H7. The circulation of AIV and H9 across the annual cycle showed no seasonality, whereas the circulation of H5 showed significant seasonality. The GEE revealed that detection of AIV increases in retail vendor business (OR: 1.71; 95% CI: 1.12-2.62) and the bird's health status is sick (OR: 10.77; 95% CI: 4.31-26.94) or dead (OR: 11.33; 95% CI: 4.30-29.89). We also observed that winter season (OR: 5.83; 95% CI: 2.80-12.14) than summer season, dead birds (OR: 61.71; 95% CI: 25.78-147.75) and sick birds (OR 8.33; 95% CI: 3.36-20.64) compared to healthy birds has a higher risk of H5 infection in turkeys. This study revealed that the turkeys movements vary by time and season from the farm to the LBM. This surveillance indicated year-round circulation of AIV with H5 and H9 subtypes in turkey birds in LBMs. The seasonality and health condition of birds influence H5 infection in birds. The trading pattern of turkey may play a role in the transmission of AIV viruses in the birds. The selling of sick turkeys infected with H5 and H9 highlights the possibility of virus transmission to other species of birds sold at LBMs and to people.

Epidemiology and molecular characterization of avian influenza virus in backyard poultry of Chattogram, Bangladesh

Ducks, the natural reservoir of avian influenza virus (AIV), act as reassortment vessels for HPAI and low pathogenic avian influenza (LPAI) virus for domestic and wild bird species. In Bangladesh, earlier research was mainly focused on AIV in commercial poultry and live bird markets, where there is scanty literature reported on AIV in apparently healthy backyard poultry at the household level. The present cross-sectional study was carried out to reveal the genomic epidemiology of AIV of backyard poultry in coastal (Anowara) and plain land (Rangunia) areas of Bangladesh. We randomly selected a total of 292 households' poultry (having both chicken and duck) for sampling. We administered structured pre-tested questionnaires to farmers through direct interviews. We tested cloacal samples from birds for the matrix gene (M gene) followed by H5 and H9 subtypes using real-time reverse transcriptase-polymerase chain reaction (rRT-PCR). All AIV-positive samples were subjected to four-gene segment sequencing (M, PB1, HA, and NA gene). We found that the prevalence of AIV RNA at the household level was 6.2% (n = 18; N = 292), whereas duck and chicken prevalence was 3.6% and 3.2%, respectively. Prevalence varied with season, ranging from 3.1% in the summer to 8.2% in the winter. The prevalence of subtypes H5 and H9 in backyard poultry was 2.7% and 3.3%, respectively. The phylogenetic analysis of M, HA, NA, and PB1 genes revealed intra-genomic similarity, and they are closely related to previously reported AIV strains in Bangladesh and Southeast Asia. The findings indicate that H5 and H9 subtypes of AIV are circulating in the backyard poultry with or without clinical symptoms. Moreover, we revealed the circulation of 2.3.2.1a (new) clade among the chicken and duck population without occurring outbreak which might be due to vaccination. In addition to routine surveillance, molecular epidemiology of AIV will assist to gain a clear understanding of the genomic evolution of the AIV virus in the backyard poultry rearing system, thereby facilitating the implementation of effective preventive measures to control infection and prevent the potential spillover to humans.

Prevalence and Distribution of Avian Influenza Viruses in Domestic Ducks at the Waterfowl-Chicken Interface in Wetlands

Ducks are a natural reservoir of influenza A viruses (IAVs) and can act as a reassortment vessel. Wetlands, such as Hakaluki and Tanguar haor in Bangladesh, have unique ecosystems including domestic duck (Anas platyrhynchos domesticus) rearing, especially household and free-range ducks. A cross-sectional study was, therefore, conducted to explore avian influenza status and its distribution and risk factors in the wetland areas. During the three consecutive winters of 2015-2017, specifically in December of these years, we collected a total of 947 samples including blood, oropharyngeal and cloacal swabs from domestic ducks (free-range duck (n = 312 samples) and household ducks (n = 635 samples) in wetlands. We screened serum samples using a nucleoprotein competitive enzyme-linked immunosorbent assay (c-ELISA) to estimate seroprevalence of IAV antibodies and swab samples by real-time reverse transcriptase polymerase chain reaction (rRT-PCR) to detect IA viral M gene. Eleven (11) M gene positive samples were subjected to sequencing and phylogenetic analysis. Serological and viral prevalence rates of IAVs were 63.8% (95% CI: 60.6-66.8) and 10.7% (8.8-12.8), respectively. Serological and viral RNA prevalence rates were 51.8% (95% CI: 47.2-56.4) and 10.2% (7.6-13.3) in Hakaluki haor, 75.6% (71.5-79.4) and 11.1% (8.5-14.3) in Tanguar haor, 66.3% (62.5-69.9) and 11.2% (8.8-13.9) in household ducks and 58.7% (52.9-64.2) and 9.6% (6.5-13.4) in free-range ducks, respectively. The risk factors identified for higher odds of AI seropositive ducks were location (OR = 2.9, 95% CI: 2.2-3.8, p < 0.001; Tanguar haor vs. Hakaluki haor), duck-rearing system (OR = 1.4, 1.1-1.8, household vs. free-range), farmer's education status (OR = 1.5, 1.2-2.0, p < 0.05 illiterate vs. literate) and contact type (OR = 3.0, 2.1-4.3, p < 0.001; contact with chicken vs. no contact with chicken). The risk factors identified for higher odds of AI RNA positive ducks were farmer's education status (OR = 1.5, 1.0-2.3, p < 0.05 for illiterate vs literate), contact type (OR = 2.7, 1.7-4.2, p < 0.001; ducks having contact with chicken vs. ducks having contact with waterfowl). The phylogenetic analysis of 11 partial M gene sequences suggested that the M gene sequences detected in free-range duck were very similar to each other and were closely related to the M gene sequences of previously reported highly pathogenic avian influenza (HPAI) and low pathogenic avian influenza (LPAI) subtypes in waterfowl in Bangladesh and Southeast Asian countries. Results of the current study will help provide significant information for future surveillance programs and model IAV infection to predict the spread of the viruses among migratory waterfowl, free-range ducks and domestic poultry in Bangladesh.

Changing risk awareness and personal protection measures for low to high pathogenic avian influenza in live-poultry markets in Taiwan, 2007 to 2012

Background

Outbreaks of low and high pathogenic avian influenza (LPAI, HPAI) H5N2 in chickens have occurred in Taiwan since 2003 and 2012, respectively. Fully understanding the different awareness, attitudes and protective behaviors adopted by workers in live-poultry markets (LPMWs) and local community residents (CRs) to face the challenges of LPAI and HPAI is very important to minimize viral adaptations to human populations.

Methods

A structural questionnaire containing information on respondents’ occupation, personal risk awareness, attitudes toward different policies, and preventative measures was administered. The two-stage survey (before and after HPAI H5N2 outbreaks) was conducted from 2007 to 2012, including: (1) 430 LPMWs and 418 CRs at LPMs from different geographical areas of Taiwan after the government announced outbreaks of LPAI H5N2 during 2007–2009, and (2) 73 LPMWs and 152 CRs at two LPMs in central Taiwan after the HPAI H5N2 outbreaks in 2012. The chi-squared test and logistic regression were applied for univariate and multivariate analyses, respectively.

Results

Before HPAI-H5N2 outbreaks, higher educated respondents demonstrated greater risk awareness and concerns regarding AI. However, LPM-workers protected themselves less from AI viruses (AIVs) and had lower acceptance of human or avian influenza vaccines. Most importantly, the participants who opposed (versus agreed with) the policy on banning live-poultry slaughtering at LPMs reported lower awareness of government prevention and control policies [Odds Ratio (OR): 0.76, 95 % Confidence Interval (CI): 0.56–1.01] or practiced preventive measures (OR: 0.42, 95 % CI: 0.25–0.70).

After HPAI-H5N2 outbreaks, the risk awareness about AI in central Taiwan significantly increased [LPAI to HPAI LPMWs: 34.6 to 65.6 %, p < 0.05; CRs: 44.0 to 76.5 %, p < 0.05] and LPMWs’ belief in the effectiveness of vaccination to prevent human or avian influenza virus infection strikingly decreased (92.3 to 68.5 %, p < 0.05).

Conclusions

Risk awareness depends on high or low pathogenicity of AIVs, working in LPMs, levels of education, age, and proximity to the sites of severe AI outbreaks. Regardless of novel LPAI or HPAI virus reassortants that pose public health risks, prompt and clear risk communication focusing on both correct information about AIVs and the most appropriate preventive measures are important for effective prevention of human infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-015-0987-8) contains supplementary material, which is available to authorized users.

Molecular evolution of H5N1 highly pathogenic avian influenza viruses in Bangladesh between 2007 and 2012

In Bangladesh, highly pathogenic avian influenza (HPAI) virus subtype H5N1 was first detected in February 2007. Since then the virus has become entrenched in poultry farms of Bangladesh. There have so far been seven human cases of H5N1 HPAI infection in Bangladesh with one death. The objective of the present study was to investigate the molecular evolution of H5N1 HPAI viruses during 2007 to 2012. Partial or complete nucleotide sequences of all eight gene segments of two chicken isolates, five gene segments of a duck isolate and the haemagglutinin gene segment of 18 isolates from Bangladesh were established in the present study and subjected to molecular analysis. In addition, full-length sequences of different gene segments of other Bangladeshi H5N1 isolates available in GenBank were included in the analysis. The analysis revealed that the first introduction of clade 2.2 virus in Bangladesh in 2007 was followed by the introduction of clade 2.3.2.1 and 2.3.4 viruses in 2011. However, only clade 2.3.2.1 viruses could be isolated in 2012, indicating progressive replacement of clade 2.2 and 2.3.4 viruses. There has been an event of segment re-assortment between H5N1 and H9N2 viruses in Bangladesh, where H5N1 virus acquired the PB1 gene from a H9N2 virus. Point mutations have accumulated in Bangladeshi isolates over the last 5 years with potential modification of receptor binding site and antigenic sites. Extensive and continuous molecular epidemiological studies are necessary to monitor the evolution of circulating avian influenza viruses in Bangladesh.

The spread of highly pathogenic avian influenza (subtype H5N1) clades in Bangladesh, 2010 and 2011

Since the global spread of highly pathogenic avian influenza H5N1 during 2005-2006, control programs have been successfully implemented in most affected countries. HPAI H5N1 was first reported in Bangladesh in 2007, and since then 546 outbreaks have been reported to the OIE. The disease has apparently become endemic in Bangladesh. Spatio-temporal information on 177 outbreaks of HPAI H5N1 occurring between February 2010 and April 2011 in Bangladesh, and 37 of these outbreaks in which isolated H5N1 viruses were phylogenetically characterized to clade, were analyzed. Three clades were identified, 2.2 (21 cases), 2.3.4 (2 cases) and 2.3.2.1 (14 cases). Clade 2.2 was identified throughout the time period and was widely distributed in a southeast-northwest orientation. Clade 2.3.2.1 appeared later and was generally confined to central Bangladesh in a north-south orientation. Based on a direction test, clade 2.2 viruses spread in a southeast-to-northwest direction, whereas clade 2.3.2.1 spread west-to-east. The magnitude of spread of clade 2.3.2.1 was greater relative to clade 2.2 (angular concentration 0.2765 versus 0.1860). In both cases, the first outbreak(s) were identified as early outliers, but in addition, early outbreaks (one each) of clade 2.2 were also identified in central Bangladesh and in northwest Bangladesh, a considerable distance apart. The spread of highly pathogenic avian influenza H5N1 in Bangladesh is characterized by reported long-distance translocation events. This poses a challenge to disease control efforts. Increased enforcement of biosecurity and stronger control of movements between affected farms and susceptible farms, and better surveillance and reporting, is needed. Although the movement of poultry and equipment appears to be a more likely explanation for the patterns identified, the relative contribution of trade and the market chain versus wild birds in spreading the disease needs further investigation.