lunes, 17 de septiembre de 2018

Molecular Evolution, Diversity, and Adaptation of Influenza A(H7N9) Viruses in China - Volume 24, Number 10—October 2018 - Emerging Infectious Diseases journal - CDC

Molecular Evolution, Diversity, and Adaptation of Influenza A(H7N9) Viruses in China - Volume 24, Number 10—October 2018 - Emerging Infectious Diseases journal - CDC





Volume 24, Number 10—October 2018

Research

Molecular Evolution, Diversity, and Adaptation of Influenza A(H7N9) Viruses in China

Changwen Ke1Comments to Author , Jayna Raghwani1, Rhys Pryce, Thomas A. Bowden, Julien Thézé, Shanqian Huang, Yingchao Song, Lirong Zou, Lijun Liang, Ru Bai, Yi Jing, Pingping Zhou, Min Kang, Lina Yi, Jie Wu2, and Oliver G. Pybus2
Author affiliations: Guangdong Provincial Center for Disease Control and Prevention, Guangzhou, China (J. Lu, Y. Song, L. Zou, L. Liang, R. Bai, Y. Jing, P. Zhou, M. Kang, L. Yi, J. Wu, C. Ke)Guangdong Provincial Institution of Public Health, Guangzhou (J. Lu, P. Zhou, L. Yi)University of Oxford, Oxford, UK (J. Raghwani, R. Pryce, T.A. Bowden, J. Thézé, O.G. Pybus)Beijing Normal University, Beijing, China (S. Huang)

Abstract

The substantial increase in prevalence and emergence of antigenically divergent or highly pathogenic influenza A(H7N9) viruses during 2016–17 raises concerns about the epizootic potential of these viruses. We investigated the evolution and adaptation of H7N9 viruses by analyzing available data and newly generated virus sequences isolated in Guangdong Province, China, during 2015–2017. Phylogenetic analyses showed that circulating H7N9 viruses belong to distinct lineages with differing spatial distributions. Hemagglutination inhibition assays performed on serum samples from patients infected with these viruses identified 3 antigenic clusters for 16 strains of different virus lineages. We used ancestral sequence reconstruction to identify parallel amino acid changes on multiple separate lineages. We inferred that mutations in hemagglutinin occur primarily at sites involved in receptor recognition or antigenicity. Our results indicate that highly pathogenic strains likely emerged from viruses circulating in eastern Guangdong Province during March 2016 and are associated with a high rate of adaptive molecular evolution.
Since its first detection in March 2013, avian influenza A(H7N9) virus has caused 1,534 human infections that, as of November 30, 2017, had resulted in 608 deaths. Recurrent waves of human cases have been reported in 27 provinces in China, indicating sustained transmission of H7N9 viruses (1). Moreover, since its emergence, H7N9 virus has reassorted with influenza A(H9N2) viruses that co-circulate in China, resulting in an increasingly diverse array of virus genomes (24). The fifth influenza epidemic wave (2016–17) was marked by a notable increase in the number of human cases (677 during September 2016–May 2017), making it the largest outbreak of influenza A(H7N9) since 2013. Moreover, geographic distribution of human cases suggests that H7N9 virus is now more widespread and that residences of patients have shifted gradually from urban to semiurban and rural areas (1,57). These epidemiologic observations have raised public health concerns.
Previous molecular surveillance studies suggested that H7N9 virus lineages originate in 2 densely populated areas, the Yangtze River Delta region in eastern China and the Pearl River Delta region in central Guangdong Province (8). Preliminary epidemiologic data suggested that most human infections in the current fifth epidemic wave were caused by viruses from the Yangtze River Delta region (5) (previously named lineage C viruses) (3). These viruses, in contrast to viruses from the Pearl River Delta region (previously named lineage B viruses) (3), appear to exhibit reduced cross-reactivity with existing candidate vaccine virus strains (9). Furthermore, a subset of lineage C isolates has also acquired a highly pathogenic (HP) phenotype (5,10,11).
These observations suggest that the increased epidemic activity of H7N9 viruses in China might be driven, at least in part, by ongoing virus evolution and adaptation. Decreased cross-reactivity and increased pathogenicity of some H7N9 viruses was discovered only recently (9), and the genetic diversity and evolution of the current fifth epidemic wave of these viruses are not yet well understood. Information necessary to clarify this situation includes geographic distribution of currently circulating H7N9 virus lineages, origin and genetic composition of newly emerged HP H7N9 viruses, and evolutionary and structural characterization of mutations associated with the fifth epidemic wave of these viruses.
We report 47 hemaglutinnin (HA) and 43 neuraminidase (NA) gene sequences of human-derived and poultry-derived H7N9 viruses that were isolated during 2015–2017 in Guangdong Province, China. We conducted structural and evolutionary analyses of these strains and characterized the evolution and emergence of currently circulating H7N9 viruses in China.

Materials and Methods

Ethics
This study was approved by the institutional ethics committee of the Center for Disease Control and Prevention of Guangdong Province. Written consent was obtained from patients or their guardian(s) when samples were collected. Patients were informed about the study before providing written consent, and data were anonymized for analysis.
Sample Collection
Samples from persons with suspected cases of influenza A(H7N9) were initially tested for avian influenza A virus in provincial clinics in Guangdong Province. Specimens with positive results were subsequently analyzed (12,13). For poultry-related samples, we obtained samples from locations where poultry were housed and processed (e.g., cages, feeding troughs, defeathering machines) (12). Respiratory specimens were collected from persons with suspected cases of influenza A(H7N9) by the Ministry of Health of China.
Sequence Alignment
For phylogenetic studies, we sequenced 47 HA and 41 NA sequences from 20 human samples and 28 poultry-related samples; all belonged to the fourth and fifth epidemic waves of influenza A(H7N9) (GISAID [https://www.gisaid.org/] accession nos. EPI866538–77, 972231–6, 972238–303, 974029, 974523, 974539–42, 997159–60, and 1171786–93). These new H7N9 sequences were combined with all available H7N9 gene sequences whose sampling dates and locations were known. Two gene sequence datasets were generated: H7, HA (n = 737) and N9, NA (n = 610). We constructed multiple sequence alignments by using ClustalW (14) and edited these sequences manually by using AliView (15).
Molecular Clock Phylogenetic Analysis
We estimated molecular clock phylogenies by using the Bayesian Markov Chain Monte Carlo approach implemented in BEAST version 1.8 (16) as described (4). We computed 4 independent Markov Chain Monte Carlo runs of 1.5 × 108 steps for each alignment and extracted a subset of 2,000 phylogenies from the posterior tree distribution, subsequently used as an empirical tree distribution for phylogeographic analyses (17). We computed maximum clade credibility trees for each dataset by using TreeAnnotator (16).
Phylogeographic Analysis of Influenza A(H7N9) Epidemic
We used the discrete phylogeographic method (18) implemented in BEAST to investigate spatial dynamics of H7N9 virus lineages from 6 regions in China as classified in a previous study (4). The 6 locations were eastern China (Anhui, Shanghai, Zhejiang, Jiangsu, and Shandong); central China (Jiangxi and Hunan); northern China (Beijing, Henan, Hebei, and Xinjiang); southeastern China (Fujian); central Guangdong Province (Guangzhou, Huizhou, Foshan, Dongguan, Zhongshan, Shenzhen, Jiangmen, Zhaoqing Yangjiang, Maoming, and Yunfu); and eastern Guangdong Province (Meizhou, Heyuan, Chaozhou, Jieyang, Shantou, Shanwei, and Shaoguan).
Because sporadic human cases detected in Malaysia and Taiwan were believed to have originated in China, we used available epidemiologic information to assign their location to the most likely source in China. Hong Kong and central Guangdong Province were treated as a single location because of their proximity to each other. We analyzed reported H7N9 virus infections and virus sequences (Technical Appendix Table 1). To estimate directionality of virus lineage movement, we used asymmetric continuous-time Markov chain phylogeographic model (19) and a Bayesian stochastic search variable selection procedure (18).
Inferring Phylogenetic Distribution of Amino Acid Changes
We investigated phylogenetic positions of amino acid changes among H7N9 virus isolates by using HA and NA maximum clade credibility trees. We estimated maximum posterior probability amino acid sequences for each internal node by using BEAST with a Jones–Taylor–Thornton amino acid substitution model (20), gamma-distributed among-site rate heterogeneity (21), and a strict molecular clock model. To infer amino acid substitutions along the trunk branches of the H7N9 phylogeny, we mapped amino acid changes onto internal branches by using a Java script (available on request). Trunk branches corresponded to internal branches that subtended >5 terminal nodes in the fifth influenza epidemic wave.
Structure-Based Mapping Analysis
We used the crystal structure of the HA (Protein Data Bank no. 4BSE) (22) and NA (Protein Data Bank no. 2C4L) glycoproteins from an influenza A(H7N9) virus to map amino acid changes identified by evolutionary analysis. We performed residue mapping onto the H7 and N9 structures by using PyMol (23). We calculated solvent accessibility for trimeric hemagglutinin with the ligands removed by using ESPript (24) and identified receptor-binding residues by using CONTACT in CCP4 (25).
Positive Selection Analyses
To identify sites under positive selection, we used methods implemented in HyPhy (26) to estimate the dN/dS ratio of codons in HA. These methods included single-likelihood ancestor counting (27), fixed effects likelihood (27), mixed effects model of evolution (28), and the fast unconstrained Bayesian approximation approach (29).
Estimating Rates of Virus Molecular Adaptation
We estimated rates of adaptive substitution in H7N9 virus HA and NA genes by using an established population genetic method related to the McDonald-Kreitman test (30,31). We used a consensus of H7N9 first-wave sequences as an outgroup to estimate derived and ancestral mutational site frequencies in each subsequent wave. Specifically, we classified polymorphisms into 3 categories according to their frequency in the population (low, 0%–15%; intermediate, 15%–75%; and high, 75%–100%). We calculated the number of adaptive substitutions from the number of synonymous and nonsynonymous sites in each category and assessed statistical uncertainty by using a bootstrap approach (1,000 replicates) (30,31).
Serologic Analysis
We obtained serum samples from 4 patients with influenza A(H7N9) 2–3 weeks after clinical symptoms were observed. We performed hemagglutination inhibition assays by using different lineages of H7N9 viruses as antigens (online Technical Appendix). Three lineage C1 strains, 4 lineage C2 strains, 5 lineage B strains, and 4 HP strains were used as antigens (Table). We calculated serum titer for each H7N9 strain as the highest reciprocal serum dilution providing complete hemagglutination inhibition.
 Top

Results

Molecular Epidemiology of Viruses Isolated during 2013–2017
Thumbnail of Regression of root-tip divergence estimated from hemagglutinin gene of influenza A(H7N9) viruses, China. Arrow indicates the time of the most recent common ancestor of the epidemic lineage.
Figure 1. Regression of root-tip divergence estimated from hemagglutinin gene of influenza A(H7N9) viruses, China. Arrow indicates the time of the most recent common ancestor of the epidemic lineage.
Thumbnail of Genetic evolution and spatial spread of epidemic lineage of influenza A(H7N9) viruses, China, 2013–2017. Bayesian maximum clade credibility tree of the hemagglutinin gene is shown. Black bars to the right of the tree indicate sequences (from waves 4 and 5) from other studies (1,5), and red bars indicate sequences reported in this study from Guangdong Province. Branch colors indicate most probable ancestral locations of each branch. Three major lineages (A, B, and C) of H7N9 viruses
Figure 2. Genetic evolution and spatial spread of epidemic lineage of influenza A(H7N9) viruses, China, 2013–2017. Bayesian maximum clade credibility tree of the hemagglutinin gene is shown. Black bars to the right of...
During 2013–2017, the influenza A(H7N9) virus epidemic lineage was geographically structured and classified into 3 major lineages, A, B, and C, in accordance with the lineage naming scheme used in a previous study (3). H7N9 virus has evolved in a clock-like manner (i.e., there is a strong linear relationship between genetic divergence and sampling time; correlation coefficient 0.95) (Figure 1). The estimated time to the most recent common ancestor (TMRCA) of H7N9 virus HA sequences is November 2012 (95% credible region October–December 2012). The corresponding molecular clock phylogeny for NA (Technical Appendix Figure 1) also shows A–C lineages and has a similar estimated TMRCA of September 2012 (95% credible region July–October 2012). However, the topology of the NA phylogeny differs from that of HA, suggesting reassortment between HA and NA during emergence of the H7N9 virus epidemic lineage (Figure 2Technical Appendix Figure 1).
Thumbnail of Geographic location and lineage classification of 374 influenza A(H7N9) human viruses, China. Values in parentheses indicate number of sequenced viruses from each region. Pie charts indicate approximate percentages of each virus lineage (A, B, C, or unclustered). Sequences from Xinjiang Province in northern China are not shown.
Figure 3. Geographic location and lineage classification of 374 influenza A(H7N9) human viruses, China. Values in parentheses indicate number of sequenced viruses from each region. Pie charts indicate approximate percentages of each virus...
Different H7N9 virus lineages are associated with different epidemiologic patterns (Figures 23). Specifically, most (86%, 32/37) lineage B viruses that were isolated during the fourth and fifth influenza epidemic waves descended from viruses circulating in central Guangdong Province during earlier epidemic seasons (Figure 2). In addition, lineage B viruses isolated from the fourth and fifth influenza waves were almost exclusively restricted to central (rather than eastern) Guangdong Province (Figures 23). In contrast, viruses in eastern China, composed of 2 lineages (A and C) have been exported to and become dominant in multiple regions as the epidemic has progressed (3). These findings indicate a comparatively broader geographic dissemination (Figure 3Technical Appendix Figure 1).
Thumbnail of Reconstruction of amino acid changes along trunk of lineage B phylogenies of influenza A(H7N9) viruses, China. Maximum clade credibility tree of hemagglutinin gene sequences from lineage B is shown. Branches are colored according to geographic locations, as in Figure 3. Thicker lines indicate the trunk lineage leading up to the current fifth influenza epidemic wave. Amino acid changes along the trunk are indicated. Red branches indicate sites undergoing parallel amino acid changes a
Figure 4. Reconstruction of amino acid changes along trunk of lineage B phylogenies of influenza A(H7N9) viruses, China. Maximum clade credibility tree of hemagglutinin gene sequences from lineage B is shown. Branches are...
The new isolates from eastern Guangdong Province, combined with isolates from eastern China (1,5), suggest that recent H7N9 virus activity is driven primarily by lineage C viruses (Figure 2). The estimated TMRCA of lineage C is December 2013 (95% highest posterior density October 2013–January 2014). For lineage C, we observed 2 clades (Figure 4). The larger of these clades (C1) circulates mainly in central and eastern China, and the smaller clade (C2) is found predominantly in eastern Guangdong Province. Clade C2 also includes recently identified HP viruses (Figures 124).
To investigate these HP viruses, we undertook retrospective screening of poultry-related samples collected in Guangdong Province during January 2016–February 2017 and identified 7 HP influenza virus isolates that belong to the HP cluster within C2 (Figure 2). These HP viruses also form a distinct cluster within lineage C viruses in the NA phylogeny Technical Appendix Figure 1). Our analyses indicated that the HP clade likely emerged from clade C2 viruses circulating in eastern Guangdong Province in 2016.
Adaptive Evolution in Virus C Lineage
Thumbnail of Reconstruction of amino acid changes along trunk of lineage C phylogenies of influenza A(H7N9) viruses, China. Maximum clade credibility tree of hemagglutinin gene sequences from lineage C is shown. Branches are colored according to geographic locations, as in Figure 3. Thicker lines indicate the trunk lineage leading up to the current fifth influenza epidemic wave. Amino acid changes along the trunk are indicated. Red branches indicate sites undergoing parallel amino acid changes a
Figure 5. Reconstruction of amino acid changes along trunk of lineage C phylogenies of influenza A(H7N9) viruses, China. Maximum clade credibility tree of hemagglutinin gene sequences from lineage C is shown. Branches are...
We then investigated whether the increasing prevalence of lineage C viruses might be associated with virus adaptation. We combined ancestral sequence reconstruction of lineage B and C HA gene sequences (Figures 45) by mapping residues that have undergone changes onto the crystal structure of the trimeric hemagglutinin. Our analysis identified several notable amino acid substitutions that occurred along the internal branches of lineage C viruses (Figure 4).
Thumbnail of Structural analysis of amino acid changes in hemagglutinin in lineages B and C of influenza A(H7N9) viruses, China. Crystal structure of the homotrimeric H7 hemagglutinin bound to a human receptor analog (Protein Data Bank no. 4BSE) (27) (A) and rotated 90° counterclockwise (B) are shown. Two of the 3 protomers are displayed with high transparency to aid visualization. The carbon Cα positions of salient features are shown as spheres. Blue indicates receptor-binding residues, red ind
Figure 6. Structural analysis of amino acid changes in hemagglutinin in lineages B and C of influenza A(H7N9) viruses, China. Crystal structure of the homotrimeric H7 hemagglutinin bound to a human receptor analog...
Around the time of the second influenza epidemic wave, ancestral viruses of lineage C acquired several amino acid changes in HA, specifically L177I, G386A, S489R, and S128N (based on H3 sequence numbering). Three of these mutations (G386A, S489R, and S128N) are located in solvent-accessible regions of HA (Figure 6Technical Appendix Table 2). Furthermore, S128N was found within the 130 loop and is proximal to the receptor surface (distance ≈20 Å) (Figure 6).
We found by evolutionary analysis that several HA sites acquired amino acid mutations independently in different phylogenetic clades. First, 4 mutations (A135V, L177I, M236I, and S489N) occurred independently along the trunk lineages that gave rise to the current lineage B and C viruses (Figure 5). Three of these mutations (A135V, M236I, and S489N) were observed only in the fourth and fifth influenza epidemic waves of lineage B (Figure 5). Second, comparison of the C1 and C2 clades also identified parallel amino acid changes within lineage C (A122T/P and M236I) (Figure 4).
The observation of parallel amino acid changes along those particular lineages (Technical Appendix Tables 2, 3) that have persisted until the fifth influenza epidemic wave (i.e., parallel changes between lineages B and C and between the C1 and C2 clades) is suggestive of convergent, adaptive molecular evolution. The parallel changes in lineage C (A122T/P and M236I) are estimated to be fully or partially solvent accessible and the A135V mutation is located at the receptor-binding site (Figure 6). One subclade of lineage B viruses appears to have acquired mutations A135V and S489N within the last 12 months (Figure 5). Therefore, we suggest that this subclade should be closely monitored in the future.
Within the C2 clade, we found that HA acquired 7 amino acid changes (I48T, A122P, K173E, L226Q, M236I, I326V, and E393K) on the internal branch immediately ancestral to the HP cluster. This internal branch represents a period of approximately 1 year (Figure 4). Although all of these changes appeared in residues with partial or full solvent accessibility, mutations K173E, L226Q, and I326V are particularly noteworthy because they have arisen at or near known antigenic, receptor-binding, and proteolytic cleavage sites, respectively (Figure 6). Furthermore, these mutations in the HP cluster also coincide with appearance of a 4-amino acid insertion (KRTA) near the HA1-HA2 proteolytic cleavage site (Figure 4). A subset of HA substitutions (at sites 57, 114, 140, 226, and 276) that occurred on the trunk branches of lineages B and C viruses was also found to be under positive selection on the basis of dN/dS ratios we estimated by using methods implemented in HyPhy (Technical Appendix Table 4).
We also investigated whether amino acid changes in the HA gene during emergence of influenza A(H7N9) virus have been driven by adaptive evolution similar to that observed for seasonal human influenza (30). We found evidence for adaptive evolution in HA genes of B and C virus lineages. We estimated that lineage B adapted at a rate of 0.80 (interquartile range [IQR] 0.21–0.95) adaptive substitutions across the whole HA gene per year and lineage A at a rate of 0.60 (IQR 0.10–1.18) adaptive substitutions per year. Within lineage C, the estimated adaptation rate of the C1 clade is ≈1.84 (IQR 1.09–2.14) adaptive substitutions per year and that for the C2 clade (which includes the HP cluster) is 3.12 (IQR 2.45–3.79) adaptive substitutions per year. These results indicate molecular adaptation across the whole H7N9 lineage and suggest that adaptation is faster in the 2 C clades than in the A and B lineages. Previous estimates of rates of adaptive substitution were 1.02 fixations per year in the whole HA gene for seasonal human influenza A(H1N1) virus and 1.52 fixations per year in the whole HA gene for influenza A(H3N2) virus (30). In this context, the rate of adaptive evolution observed for lineage C here is notable and raises concern for ongoing evolution of these viruses.
Antigenic Properties
We collected serum samples from 4 patients infected with H7N9 virus during 2015 and 2017 (Table). For patients 3 and 4, the corresponding virus strains were isolated and sequenced. Phylogenetic analysis indicated that patient 3 was infected with clade C1 virus and that patient 4 was infected with HP virus. Hemagglutination inhibition results suggested the presence of 3 antigenic clusters among the 16 H7N9 virus strains selected. Serum samples from patients 1, 2, and 3 showed similar patterns, reacting robustly to clade C1 viruses and moderately to clade C2 and lineage B viruses but poorly to HP viruses. A serum sample from a patient infected with an HP H7N9 virus appeared to react strongly to all H7N9 virus strains.
 Top

Discussion

Our results show that H7N9 viruses of lineage C, which were prevalent in the recent fifth influenza epidemic wave in China, comprise 2 geographically distinct clades (C1 and C2) that have undergone adaptive evolution. Clade C1 is found primarily in eastern and central China and clade 2 in Guangdong Province, and both clades appear to have circulated in bird populations for ≈3 years. Our ancestral state reconstruction analysis provides crucial evidence that 2 successful lineages of H7N9 viruses (lineages B and C) have experienced multiple parallel amino acid changes (Figures 45), suggesting the possible action of convergent molecular evolution.
We also observed a higher rate of virus adaptation in eastern Guangdong Province (C2 clade compared with C1). Although clades C1 and C2 are phylogenetically closely related, serum from a clade C1 virus-infected patient has moderate reactivity with C2 strains from 2015–2016 and poor reactivity to the HP virus from 2016–2017. The higher adaptation rate and antigenic changes in clade C2 are of concern from a public health perspective. Introduction of HP avian influenza into domestic poultry might constitute a serious risk, as demonstrated by emergence of goose–Guangdong lineage HP H5N1 viruses, which spilled back into wild birds and caused the longest global outbreak of HP avian influenza to date (33).
Parallel amino acid changes in clades C1 and C2 occurred at 2 sites in HA (122 and 236) (Figure 6). Furthermore, we observed 4 mutations that emerged independently in lineages B and C (sites 135, 177, 236, and 489). These results suggest adaptive convergent molecular evolution. Site 135 is located in the receptor-binding region and is near antigenic site A, as defined by Wiley et al. (34). Thus, the observed A135V mutation might modulate receptor affinity and contribute to immune escape (Figure 6), as observed in influenza A(H7N1) and A(H7N7) viruses (35,36).
Specifically, experimental studies indicate that threonine at position 135 in the HP H7N7 virus A/Netherlands/219/2003 confers broad-scale resistance to neutralizing monoclonal antibodies against the earliest strain of H7N9 virus (A/Shanghai/02/2013) (37). Furthermore, the World Health Organization has reported that recent clade C1 viruses (but not those of lineage B) react less to postinfection ferret antiserum raised against the A/Anhui/1/2013 and A/Shanghai/2/2013-derived candidate vaccine strains (9). Consistent with this finding, we found that most clade C1 viruses isolated in 2015 have the A135V mutation. However, this mutation was only detected in a small proportion of recent lineage B viruses (Figure 5).
In this study, we performed a preliminary evaluation of the antigenicity of H7N9 viruses by using patient serum samples collected in 2015 and 2017. Without serum raised in response to early strains from 2013, we cannot discriminate antigenic change between strains from 2013 and those from 2017. However, the limited data we have indicate the presence of 3 antigenic clusters among the 4 phylogenetic clusters circulating during the fifth influenza epidemic wave (Table). HA1 positions 109–301 (H3 numbering) include the A–E antigenic epitopes, which are known to determine antigenicity of influenza A viruses (34). The amino acid changes responsible for the antigenic differences between clade C1 and other clades were located in antigenic site A (position 140, H3 numbering; Technical Appendix Figure 2).
The mutation R140K has been observed in viruses isolated from ferrets experimentally infected with avian influenza A(H7N9) viruses and has been linked to antigenic drift of influenza A(H5N1) viruses (3840). By comparing the sequences of the HP H7N9 virus cluster and other clade C2 viruses, we found a substitution in antigenic site E (position 173, H3 numbering) that could underlie antigenic change in HP H7N9 viruses (Technical Appendix Figure 2). In future work, we aim to explore the roles of these substitutions in determining viral antigenicity in the context of H7N9 virus genomes by using reverse genetics.
HP H7N9 viruses belonging to lineage C2 were first detected in late 2016, but spread greatly in geographic range in early 2017 (41). Several mechanisms for the genesis of an HP virus from a low pathogenicity virus have been proposed, including transcription errors (42), stepwise amino acid substitutions (43), or recombination (44). For H7, emergence of HP viruses is attributed to nonhomologous recombination resulting in simultaneous insertion of several amino acids at the HA cleavage site. These insertions might be derived from host 28S rRNA sequence (45) or from other influenza gene segments, such as the matrix (46) and nucleoprotein genes (44). The 12-nt insert in the HP H7N9 virus strains is 100% identical to a region in the polymerase basic 1 gene in multiple avian influenza A viruses, including subtypes H3N2, H6N2, and H9N2, but is not present in the polymerase basic 1 gene of HP H7N9 virus. H9N2 virus is the most frequently detected avian influenza virus in chickens in China, and the detection rate of this subtype in environmental samples from live poultry markets is ≈20% during the influenza epidemic season (13). Therefore, co-infection with H7N9 and other avian influenza viruses, such as influenza A(H9N2) viruses, could, in theory, lead to insertion of a polybasic cleavage site by nonhomologous recombination.
Recent studies have shown that the HP H7N9 virus is more pathogenic in mice, and more thermally stable, than low pathogenicity A/Anhui/1/2013 virus (47,48). Current surveillance indicates that HP H7N9 viruses have spread to several provinces in China and are responsible for large influenza outbreaks in poultry in central and northern China that show high mortality rates (http://www.fao.org/ag/againfo/programmes/en/empres/H7N9/situation_update.html). This finding raises the possibility of global dissemination of H7N9 viruses through migration of wild birds, in a manner similar to that observed for HP H5N1 viruses first identified in Guangdong Province (32). Although vaccination of poultry against H7N9 virus has now been implemented in some regions of China, virus adaptation and spatial distribution should be more closely monitored.
 Top
Dr. Lu is a virologist at the Guangdong Provincial Center for Disease Control and Prevention, Guangzhou, China. His primary research interests are epidemiology, evolution, and transmission of viruses associated with emerging infectious diseases.
 Top

Acknowledgments

We thank Tommy Lam for performing sequence alignments.
This study was supported by the National Natural Science Foundation of China (grant 81501754) and the National Key Research and Development Program of China (grant 2016YFC1200201). T.A.B. is supported by the Medical Research Council of the United Kingdom (grant MR/L009528/1). The Wellcome Trust Centre for Human Genetics is supported by the Wellcome Trust Centre (grant 203141/Z/16/Z).
C.K., J.W., and J.L. designed the study; J.L., Y.S., L.Z., L.L., R.B., Y.J., P.Z., M.K., and L.Y. collected samples and performed genome sequencing; J.L., J.R., R.P., T.A.B., J.T., S.H., and O.G.P. analyzed data; J.L., J.R., R.P., T.A.B., and O.G.P. interpreted data; J.L., J.R., R.P., T.A.B., and O.G.P. prepared the figures; and J.L., J.R., R.P., T.A.B., and O.G.P. wrote the article.
 Top

References

  1. Xiang  NLi  XRen  RWang  DZhou  SGreene  CMet al. Assessing change in avian influenza A(H7N9) virus infections during the fourth epidemic—China, September 2015–August 2016.MMWR Morb Mortal Wkly Rep2016;65:13904DOIPubMed
  2. Cui  LLiu  DShi  WPan  JQi  XLi  Xet al. Dynamic reassortments and genetic heterogeneity of the human-infecting influenza A (H7N9) virus. Nat Commun2014;5:3142DOIPubMed
  3. Lam  TTZhou  BWang  JChai  YShen  YChen  Xet al. Dissemination, divergence and establishment of H7N9 influenza viruses in China. Nature2015;522:1025DOIPubMed
  4. Wu  JLu  JFaria  NRZeng  XSong  YZou  Let al. Effect of live poultry market interventions on influenza A(H7N9) virus, Guangdong, China. Emerg Infect Dis2016;22:210412DOIPubMed
  5. Iuliano  ADJang  YJones  JDavis  CTWentworth  DEUyeki  TMet al. Increase in human infections with avian influenza A(H7N9) virus during the fifth epidemic—China, October 2016–February 2017. MMWR Morb Mortal Wkly Rep2017;66:2545DOIPubMed
  6. Wang  XJiang  HWu  PUyeki  TMFeng  LLai  Set al. Epidemiology of avian influenza A H7N9 virus in human beings across five epidemics in mainland China, 2013-17: an epidemiological study of laboratory-confirmed case series. Lancet Infect Dis2017;17:82232DOIPubMed
  7. Artois  JJiang  HWang  XQin  YPearcy  MLai  Set al. Changing geographic patterns and risk factors for avian influenza A(H7N9) infections in humans, China. Emerg Infect Dis2018;24:8794DOIPubMed
  8. Wang  DYang  LZhu  WZhang  YZou  SBo  Het al. Two outbreak sources of influenza A (H7N9) viruses have been established in China. J Virol2016;90:556173DOIPubMed
  9. Zhu  WZhou  JLi  ZYang  LLi  XHuang  Wet al. Biological characterisation of the emerged highly pathogenic avian influenza (HPAI) A(H7N9) viruses in humans, in mainland China, 2016 to 2017. Euro Surveill2017;22:30533DOIPubMed
  10. Ke  CMok  CKPZhu  WZhou  HHe  JGuan  Wet al. Human infection with highly pathogenic avian influenza A(H7N9) virus, China. Emerg Infect Dis2017;23:133240DOIPubMed
  11. Zhang  FBi  YWang  JWong  GShi  WHu  Fet al. Human infections with recently-emerging highly pathogenic H7N9 avian influenza virus in China. J Infect2017;75:715DOIPubMed
  12. Ke  CLu  JWu  JGuan  DZou  LSong  Tet al. Circulation of reassortant influenza A(H7N9) viruses in poultry and humans, Guangdong Province, China, 2013. Emerg Infect Dis2014;20:203440DOIPubMed
  13. Lu  JWu  JZeng  XGuan  DZou  LYi  Let al. Continuing reassortment leads to the genetic diversity of influenza virus H7N9 in Guangdong, China. J Virol2014;88:8297306DOIPubMed
  14. Larkin  MABlackshields  GBrown  NPChenna  RMcGettigan  PAMcWilliam  Het al. Clustal W and Clustal X version 2.0. Bioinformatics2007;23:29478DOIPubMed
  15. Larsson  AAliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics2014;30:32768DOIPubMed
  16. Drummond  AJSuchard  MAXie  DRambaut  ABayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol2012;29:196973DOIPubMed
  17. Lemey  PRambaut  ABedford  TFaria  NBielejec  FBaele  Get al. Unifying viral genetics and human transportation data to predict the global transmission dynamics of human influenza H3N2. PLoS Pathog2014;10:e1003932DOIPubMed
  18. Lemey  PRambaut  ADrummond  AJSuchard  MABayesian phylogeography finds its roots. PLOS Comput Biol2009;5:e1000520DOIPubMed
  19. Edwards  CJSuchard  MALemey  PWelch  JJBarnes  IFulton  TLet al. Ancient hybridization and an Irish origin for the modern polar bear matriline. Curr Biol2011;21:12518DOIPubMed
  20. Jones  DTTaylor  WRThornton  JMThe rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci1992;8:27582.PubMed
  21. Yang  ZAmong-site rate variation and its impact on phylogenetic analyses. Trends Ecol Evol1996;11:36772DOIPubMed
  22. Xiong  XMartin  SRHaire  LFWharton  SADaniels  RSBennett  MSet al. Receptor binding by an H7N9 influenza virus from humans. Nature2013;499:4969DOIPubMed
  23. Schrödinger  L. The PyMOL molecular graphics system. Version 1.8. New York: Schrödinger, LLC; 2015.
  24. Gouet  PCourcelle  EStuart  DIMétoz  FESPript: analysis of multiple sequence alignments in PostScript. Bioinformatics1999;15:3058DOIPubMed
  25. Winn  MDBallard  CCCowtan  KDDodson  EJEmsley  PEvans  PRet al. Overview of the CCP4 suite and current developments. Acta Crystallogr D Biol Crystallogr2011;67:23542DOIPubMed
  26. Pond  SLFrost  SDMuse  SVHyPhy: hypothesis testing using phylogenies. Bioinformatics2005;21:6769DOIPubMed
  27. Kosakovsky Pond  SLFrost  SDNot so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol2005;22:120822DOIPubMed
  28. Murrell  BWertheim  JOMoola  SWeighill  TScheffler  KKosakovsky Pond  SLDetecting individual sites subject to episodic diversifying selection. PLoS Genet2012;8:e1002764DOIPubMed
  29. Murrell  BMoola  SMabona  AWeighill  TSheward  DKosakovsky Pond  SLet al. FUBAR: a fast, unconstrained bayesian approximation for inferring selection. Mol Biol Evol2013;30:1196205DOIPubMed
  30. Bhatt  SHolmes  ECPybus  OGThe genomic rate of molecular adaptation of the human influenza A virus. Mol Biol Evol2011;28:244351DOIPubMed
  31. Raghwani  JBhatt  SPybus  OGFaster adaptation in smaller populations: counterintuitive evolution of HIV during childhood infection. PLOS Comput Biol2016;12:e1004694DOIPubMed
  32. Kapczynski  DRPantin-Jackwood  MGuzman  SGRicardez  YSpackman  EBertran  Ket al. Characterization of the 2012 highly pathogenic avian influenza H7N3 virus isolated from poultry in an outbreak in Mexico: pathobiology and vaccine protection. J Virol2013;87:908696DOIPubMed
  33. Subbarao  KKlimov  AKatz  JRegnery  HLim  WHall  Het al. Characterization of an avian influenza A (H5N1) virus isolated from a child with a fatal respiratory illness. Science1998;279:3936DOIPubMed
  34. Wiley  DCWilson  IASkehel  JJStructural identification of the antibody-binding sites of Hong Kong influenza haemagglutinin and their involvement in antigenic variation. Nature1981;289:3738DOIPubMed
  35. Monne  IFusaro  ANelson  MIBonfanti  LMulatti  PHughes  Jet al. Emergence of a highly pathogenic avian influenza virus from a low-pathogenic progenitor. J Virol2014;88:437588DOIPubMed
  36. de Wit  EMunster  VJvan Riel  DBeyer  WERimmelzwaan  GFKuiken  Tet al. Molecular determinants of adaptation of highly pathogenic avian influenza H7N7 viruses to efficient replication in the human host. J Virol2010;84:1597606DOIPubMed
  37. Schmeisser  FVasudevan  AVerma  SWang  WAlvarado  EWeiss  Cet al. Antibodies to antigenic site A of influenza H7 hemagglutinin provide protection against H7N9 challenge. PLoS One2015;10:e0117108DOIPubMed
  38. Belser  JAGustin  KMPearce  MBMaines  TRZeng  HPappas  Cet al. Pathogenesis and transmission of avian influenza A (H7N9) virus in ferrets and mice. Nature2013;501:5569DOIPubMed
  39. Cattoli  GMilani  ATemperton  NZecchin  BBuratin  AMolesti  Eet al. Antigenic drift in H5N1 avian influenza virus in poultry is driven by mutations in major antigenic sites of the hemagglutinin molecule analogous to those for human influenza virus. J Virol2011;85:871824DOIPubMed
  40. Xu  LBao  LDeng  WDong  LZhu  HChen  Tet al. Novel avian-origin human influenza A(H7N9) can be transmitted between ferrets via respiratory droplets. J Infect Dis2014;209:5516DOIPubMed
  41. Yang  LZhu  WLi  XChen  MWu  JYu  Pet al. Genesis and spread of newly emerged highly pathogenic H7N9 avian viruses in mainland China. J Virol2017;91:e0127717.PubMed
  42. García  MCrawford  JMLatimer  JWRivera-Cruz  EPerdue  MLHeterogeneity in the haemagglutinin gene and emergence of the highly pathogenic phenotype among recent H5N2 avian influenza viruses from Mexico. J Gen Virol1996;77:1493504DOIPubMed
  43. Horimoto  TRivera  EPearson  JSenne  DKrauss  SKawaoka  Yet al. Origin and molecular changes associated with emergence of a highly pathogenic H5N2 influenza virus in Mexico.Virology1995;213:22330DOIPubMed
  44. Suarez  DLSenne  DABanks  JBrown  IHEssen  SCLee  CWet al. Recombination resulting in virulence shift in avian influenza outbreak, Chile. Emerg Infect Dis2004;10:6939DOIPubMed
  45. Kapczynski  DRPantin-Jackwood  MGuzman  SGRicardez  YSpackman  EBertran  Ket al. Characterization of the 2012 highly pathogenic avian influenza H7N3 virus isolated from poultry in an outbreak in Mexico: pathobiology and vaccine protection. J Virol2013;87:908696DOIPubMed
  46. Pasick  JHandel  KRobinson  JCopps  JRidd  DHills  Ket al. Intersegmental recombination between the haemagglutinin and matrix genes was responsible for the emergence of a highly pathogenic H7N3 avian influenza virus in British Columbia. J Gen Virol2005;86:72731DOIPubMed
  47. Shi  JDeng  GKong  HGu  CMa  SYin  Xet al. H7N9 virulent mutants detected in chickens in China pose an increased threat to humans. Cell Res2017;27:140921DOIPubMed
  48. Imai  MWatanabe  TKiso  MNakajima  NYamayoshi  SIwatsuki-Horimoto  Ket al. A highly pathogenic avian H7N9 influenza virus isolated from a human is lethal in some ferrets infected via respiratory droplets. Cell Host Microbe2017;22:615626.e8DOIPubMed
 Top

Figures

Table

 Top
Cite This Article

DOI: 10.3201/eid2410.171063
1These authors contributed equally to this article.
2These senior authors jointly supervised this study.

No hay comentarios:

Publicar un comentario