Influenza A(H7N9) virus emerged in eastern China in February 2013 and continues to circulate in this region, but its ecology is poorly understood. In April 2013, the Guangdong Provincial Center for Disease Control and Prevention (CDC) implemented environmental and human syndromic surveillance for the virus. Environmental samples from poultry markets in 21 city CDCs (n = 8,942) and respiratory samples from persons with influenza-like illness or pneumonia (n = 32,342) were tested; viruses isolated from 6 environmental samples and 16 patients were sequenced. Sequence analysis showed co-circulation of 4 influenza A(H7N9) virus strains that evolved by reassortment with avian influenza A(H9N2) viruses circulating in this region. In addition, an increase in human cases starting in late 2013 coincided with an increase in influenza A H7 virus isolates detected by environmental surveillance. Co-circulation of multiple avian influenza viruses that can infect humans highlights the need for increased surveillance of poultry and potential environmental sources.
Human infection with a novel avian-origin influenza A(H7N9) virus was first identified in eastern China in February 2013, and a major outbreak occurred from the end of March through the beginning of May (
After H7N9 virus was initially detected in eastern China, the Guangdong Provincial Center for Disease Control and Prevention (CDC) implemented environmental surveillance and syndromic surveillance in humans. After H7N9 virus infection in humans was reported, enhanced sentinel hospital surveillance and environmental sampling programs were implemented to identify possible cases and analyze the evolution of the virus. The first case of human infection with H7N9 virus in Guangdong was confirmed in a patient with severe pneumonia on August 10, 2013. This case-patient was a 51-year-old female poultry worker at a local LPM in Huizhou city in Guangdong Province. No other cases were reported until October 2013, when a wave of cases began in several provinces, including Guangdong. We present results of ongoing environmental and human syndromic surveillance for H7N9 virus circulation in Guangdong Province, as well as virologic analyses related to the natural history of this outbreak.
Since April 16, 2013, enhanced surveillance for influenza A(H7N9) virus has been conducted in 28 sentinel hospitals and 23 collaborating laboratories in 21 cities in Guangdong Province, China. All specimens are first tested for avian influenza A virus; specimens with positive results are then further tested in the laboratories of the 21 city CDCs for subtypes H5, H7, and H9 by using commercial real-time reverse transcription PCR (RT-PCR) (Liferiver, Shanghai ZJ Bio-Tech, Shanghai, China). Results are further verified at the Guangdong Provincial CDC. H7-positive specimens are tested for N9 gene segments by using real-time RT-PCR (
In April 2013, the Guangdong Provincial CDC launched an environmental surveillance program to monitor for avian influenza viruses in LPMs. Environmental samples were collected from LPMs in Guangdong Province during April 15, 2013–February 28, 2014. Each week during April 15–May 31, 2013, twenty environmental samples per market were collected from selected markets in 21 cities at the prefecture level. The environmental sampling was done by collecting wet-swab specimens from poultry feces, chicken epilator surfaces, chopping board surfaces, cage surfaces, and sewage. If a human H7N9 case was confirmed and the case-patient had exposure in an LPM, >20 environmental samples would be collected from that market.
Samples were collected from poultry feces by swabbing the surfaces of the chicken epilator, chopping boards, and cages 4–8 times with separate cotton-tipped swabs (Copan Italia, Brescia, Italy). The swabs were then inserted into a tube containing 3 mL of virus transport medium (Copan Italia) and stored at 4°C before shipping to the CDC laboratory. At arrival in the laboratory, samples were vortexed and swabs were discarded. Total RNA and DNA were extracted simultaneously from 200 μL supernatant of the sample by using the QIAamp minElute Virus Spin Kit (QIAGEN, Crawley, UK), according to the manufacturer’s instructions. Then, 5 μL of eluate was tested for influenza A; thereafter, testing for subtypes H5, H7, and H9 was done by commercial real-time reverse transcription PCR (RT-PCR) (Liferiver, Shanghai ZJ Bio-Tech). H7-positive specimens were further subjected to N9 gene-specific RT-PCR. Dual positive samples were blindly passaged for 2–3 generations in 9- to 10-day-old embryonated chicken eggs for virus isolation. Hemagglutinin (HA)–positive isolates were collected and further subtyped by HI and NI assays using a panel of reference antiserum, as described by Huang et al. (
For all selected isolates, all 8 gene segments (HA, neuraminidase [NA], nucleoprotein [NP], polymerase basic 1 [PB1] and 2 [PB2], polymerase acidic [PA], matrix [M], and nonstructural [NS]) were sequenced by using a next-generation sequencing strategy for influenza A virus sequencing with the Ion PGM System and PathAmp FluA Reagents (Life Technologies, Carlsbad, CA, USA). Viral RNA extraction was performed by using the QIAamp Viral RNA Mini Kit (QIAGEN, Hilden, Germany). Reverse transcription PCR amplification of all 8 gene segments were performed by using PathAmp Flu A Preamplification Reagents (Life Technologies). Amplicons were purified and quantitated by using the Ampure XP purification kit (Beckman, Brea, CA, USA). The genomic libraries were prepared with the Ion Xpress Plus Fragment Library Kit (Life Technologies). Enzymatic fragmentation was used for the 200-bp read protocol with a 10-min incubation time. Samples were assigned barcodes by using the Ion Xpress Barcode Adapters 1–32 Kit (Life Technologies). Automated template preparation was performed by using the Ion OneTouch 2 System (Ion PGM Template OT2 200 Kit; Life Technologies). Final sequencing was performed with the Ion PGM Sequencing 200 Kit v2 (Life Technologies) using the Ion 316 Chip V2.
Genome sequence assembly was done with the pathogen analysis program of the Ion PGM server (Life Technologies). Multiple sequence alignment against previously published complete genome sequences of H7N9 and H9N2 viruses was performed by using ClustalW (
Maximum-likelihood trees were estimated for all 8 gene segments by using MEGA version 5.2 (
The initial trees were constructed in PhyML by using general time reversible + I + Γ4 as the evolutionary model (
A total of 34,342 samples from patients with pneumonia or influenza-like illness from 21 cities in Guangdong Province were tested during April 2013–February 2014. Of these, samples from 81 patients from 10 cities tested positive for influenza A(H7N9) virus. With the exception of 1 case confirmed in August 2013, all other cases were found starting in October 2013. In total, 6 cases were confirmed in 2013 and 75 in early 2014.
Patient characteristics are summarized in
| Characteristic | No. (%) case-patients |
|---|---|
| Sex* | |
| M | 48 (59.3) |
| F | 33 (40.7) |
| Age group, y† | |
| ≤20 | 12 (14.8) |
| 21–40 | 13 (16.0) |
| 41–60 | 25 (30.9) |
| ≥61 | 31 (38.3) |
| Exposure history | |
| Live poultry market visit | 47 (58.0) |
| Household poultry breeding | 10 (12.3) |
| Poultry processing | 5 (6.2) |
| Mixed | 11 (13.6) |
| None known | 8 (9.9) |
| Clinical use of oseltamivir‡ | |
| Used | 65 (80.2) |
| Unknown | 16 (19.8) |
| Clinical condition | |
| Influenza-like illness | 11 (13.6) |
| Pneumonia | 14 (17.3) |
| Mild influenza | 2 (2.5) |
| Severe influenza | 54 (66.7) |
| Outcome | |
| Death | 22 (27.2) |
| Recovery | 45 (55.6) |
| Hospitalization | 14 (17.3) |
*M:F sex ratio 1.5:1. †Median patient age 52.0 y. ‡Mean time from symptom onset to diagnosis 8.5 d; mean time from symptom onset to use of oseltamivir 5.2 d.
A total of 8,942 swab samples from LPMs were collected and screened in laboratories of city CDCs at the prefecture level and at Guangdong Provincial CDC (
Distribution of influenza A H9 (A) and H7 (B) viruses, Guandong Province, China, 2013–2014. Shading indicates percentage of environmental swab specimens from live poultry markets in each region that were positive for each influenza subtype by reverse transcription PCR. Circles indicate locations of human cases; larger circles indicate higher numbers of cases.
Whole-genome sequences were determined for the H7N9 virus isolates from 16 patients and 6 environmental samples and for the 9 H9 viruses collected during environmental surveillance. Phylogenetic trees were constructed for each internal gene segment against all currently available H7N9 and H9N2 virus sequences from the National Center for Biotechnology Information and GISAID. Phylogenetic analysis of the whole-genome sequences showed the clear separation of PB2, PB1, and NS gene sequences of the Guangdong patient isolates from those detected from patients from eastern China (
Timed phylogenies of internal gene segments of influenza A(H7N9) viruses detected in 16 patients from Guangdong Province, China, 2008–2014. A) Polymerase basic 2; B) polymerase basic 1; C) polymerase acidic; D) nucleoprotein; E) matrix; F) nonstructural. Avian H9N2 viruses are shown in green, avian H9N2 viruses from Guangdong Province in blue, human H7N9 sequences from Guangdong Province in pink, and human H7N9 sequences from the main case cluster in 2013 in eastern China in red. Blue boxes indicate clusters of patient sequences on separate branches; numbers correspond with the reassortant numbers listed in
Gene sequences on the same branch were reviewed to identify region of origin; these sequences matched internal genes from H9N2 viruses detected in southern China (
| Sample identification | PB2 variant | PB1 variant | PA variant | NP variant | M variant | NS variant |
|---|---|---|---|---|---|---|
| Patient samples | ||||||
| HZ-Y/01/2013 | 1 | 1 | 0 | 1 | 0 | 1 |
| FS-Y/019/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| FS2-Y/033/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| FS2-T/034/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| SZ-T/035/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| SZ-Y/036/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| A/Hong Kong/734/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| DG-Y/02/2013 | 1 | 1 | 1 | 1 | 0 | 1 |
| DG-T/03/2013 | 1 | 1 | 1 | 1 | 0 | 1 |
| YJ-T/04/2013 | 1 | 1 | 1 | 1 | 0 | 1 |
| YJ-T/05/2013 | 1 | 1 | 1 | 1 | 0 | 1 |
| YJ2-T/012/2014 | 1 | 1 | 1 | 1 | 0 | 1 |
| YJ2-Y/013/2014 | 1 | 1 | 1 | 1 | 0 | 1 |
| FS-Y/029/2014 | 1 | 1 | 1 | 1 | 0 | 1 |
| Hong Kong 2212982 | 1 | 1 | 1 | 1 | 0 | 1 |
| GZ-T/010/2014 | 1 | 1 | 2 | 1 | 0 | 1 |
| FS-Y/031/2014 | 1 | 1 | 2 | 1 | 0 | 1 |
| SZ-T/026/2014 | 1 | 1 | 0 | 2 | 0 | 1 |
| Environmental samples | ||||||
| Guangdong/02620/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| Guangdong/24997/2014 | 1 | 1 | 0 | 1 | 0 | 1 |
| Guangdong/0092/2014 | 1 | 1 | 1 | 1 | 0 | 1 |
| Guangdong/02124/2014 | 1 | 1 | 1 | 1 | 0 | 1 |
| Guangdong/25003/2014 | 1 | 1 | 1 | 1 | x | 1 |
| Guangdong/02125/2014 | 1 | 1 | 2 | 1 | 0 | 1 |
| Variant origin | SC | SC | PA-1: SC | NP-1: SC | SC | |
| PA-2: EC | NP-2:EC |
*Lineages shown correlated with those shown
The HA and NA genes of all H7N9 isolates from Guangdong Province we analyzed had high sequence identity (98.2%–99.7%) to viruses from other regions of China detected throughout 2013. However, the Guangdong Province isolates branched off and clustered with strains detected in Hong Kong during 2014 (
The outbreak of influenza A(H7N9) virus infection in humans in Guangdong Province that started in late 2013 coincided with the emergence of influenza A H7 viruses in LPMs as determined by environmental surveillance. The strong association between detection of human H7N9 cases and presence of H7 viruses in LPMs confirms that these market environments are sources of human exposure (
The data from our study suggest sources of human infection in Guangdong Province are local. All sequenced viruses had evolved from the original strains through genetic reassortment with internal genes originating from influenza A(H9N2) viruses commonly found in southern China. The HA sequences of these viruses formed a distinct branch in phylogenetic trees, which indicates that reassortment occurred some time ago and suggests that local sources of virus transmission in Guangdong Province have not been detected by surveillance activities. However, available data are insufficient to determine how long influenza A(H7N9) viruses have been circulating in this region.
The results from environmental surveillance suggest a higher prevalence of H9N2 virus in Guangdong Province poultry environments than in poultry environments in other regions of China (
In summary, our findings indicate that multiple strains of H7N9 and H9N2 influenza viruses are circulating in poultry in Guangdong Province, creating an environment that is rich for reassortment of these viruses and that poses an ongoing risk for human infection. Continued and uncontrolled co-circulation of multiple avian and other influenza viruses that can infect humans pose a potential pandemic threat; increased animal surveillance and further study of the ecology of influenza viruses are essential.
Phylogenetic trees of all gene segments of currently available influenza A(H7N9) viruses, showing locations of H7N9 and H9N2 viruses from Guangdong Province, China.
This work was supported by 12th five-year-major projects of China’s Ministry of Public Health (grant no. 2012zx10004-213); the PREDICT Surveillance Animal Human Interface Project (grant no. Gvf: 06-09-057-02, to G.V.F.); and by the CASTELLUM grant, Dutch Ministry of Economic Affairs, Agriculture and Innovation (to M.K.).
Dr Ke is a senior scientist at Guangdong Provincial CDC, Guangdong, China, and is involved in the diagnosis and molecular epidemiology research of zoonoses and emerging infectious diseases.