Genetic drift evolution under vaccination pressure among H5N1 Egyptian isolates
Virology Journalvolume 8, Article number: 283 (2011)
The highly pathogenic H5N1 is a major avian pathogen that intensively affects the poultry industry in Egypt even in spite of the adoption of vaccination strategy. Antigenic drift is among the strategies the influenza virus uses to escape the immune system that might develop due to the pressure of extensive vaccination. H5N1 mutates in an intensified manner and is considered a potential candidate for the possible next pandemic with all the catastrophic consequences such an eventuality will entail.
H5N1 was isolated from the pooled organ samples of four different affected flocks in specific pathogen free embryonated chicken eggs (SPF-ECE). A reverse transcriptase polymerase chain reaction (RT-PCR) was performed to the haemagglutingin and neuraminidase. Sequencing of the full length haemagglutingin was performed. Sequence analyses of the isolated strains were performed and compared to all available H5N1 from Egyptian human and avian strains in the flu database. Changes in the different amino acid that may be related to virus virulence, receptor affinity and epitope configuration were assigned and matched with all available Egyptian strains in the flu database.
One out of the four strains was found to be related to the B2 Egyptian lineage, 2 were related to A1 lineage and the 4th was related to A2 lineage. Comparing data obtained from the current study by other available Egyptian H5N1 sequences remarkably demonstrates that amino acid changes in the immune escape variants are remarkably restricted to a limited number of locations on the HA molecule during antigenic drift. Molecular diversity in the HA gene, in relevance to different epitopes, were not found to follow a regular trend, suggesting abrupt cumulative sequence mutations. However a number of amino acids were found to be subjected to high mutation pressure.
The current data provides a comprehensive view of HA gene evolution among H5N1 subtype viruses in Egypt. Egyptian H5N1-AIVs are constantly undergoing genetic changes and reveal a complex pattern of drifts. These findings raise the concerns about the value of using influenza vaccines in correlation with the development of antigenic drift in influenza epidemics.
The highly pathogenic H5N1 influenza virus is currently panzootic in Egyptian poultry populations and crosses species barriers to humans and animals. H5N1 viruses are classified into 10 distinct initial clades (numbered 0-9) based on the HA gene . H5N1 virus strains from European-Middle Eastern-African (EMA) origin were assigned to three clades (EMA-1-3)  and referred to as sub-clades 2.2.1.-2.2.3 . Egyptian H5N1 viruses were exclusively classified as clade 2.2.1, which were further subdivided into two major groups A (A1-A5) and B (B1-B2) . Following the widespread of HPAI H5N1 in Egypt, authorities began culling and vaccination to control the spread of the disease in poultry. H5N1 and several H5N2 inactivated vaccines were imported in a trial to control the H5N1 in Egypt. Nevertheless, it still endemic in the country and continues to constitute a major challenge to poultry population. The long-term pattern of endemic status increases the opportunity for the emergence of potential pandemic strains through further adaptation by genetic mutation or reassortment . The HA gene is a target of neutralizing antibodies and a classic example of an antigenically drifting protein where a large number of point mutations accumulate in the epitope regions or antibody combining regions . Antigenic drift at the epitope regions represents the most important hindrance to vaccine development since an effective vaccination is possible only when the epidemic strain matches with the vaccine strain . In additions, suboptimal vaccination strategies could result in enhanced antigenic drift accelerated by the immune pressure due to prior immunization exerted on the replicating viruses and continuous interspecies and intraspecies transmission . However, it is generally accepted that the highly error-prone replication of influenza viruses and viral genome reassortment enhance the robustness of an efficient evolutionary capacity of the virus [10, 11]
In the present paper, we isolated and characterized four H5N1 strains from Egyptian commercial chicken farms to get a deep insight into possible intra- and inter-subtype sequence variations. The genetic and phylogenetic properties of these viruses were determined and compared to each other and to all available Egyptian published sequences in the flu database to determine the vigor of the mutational events of HA. We tried to address the mutations change over time in HA in correlation to their functional impact.
Materials and methods
Liver, spleen and proventriculus organ samples were collected from succumbed birds from four infected poultry farms (F7-F10) in three different Governorates (Fayoum, Qaluobiya and Giza). Birds in all farms received the H5N2 or H5N1 vaccine 2-3 weeks before developing the characteristic signs of the disease. Organs were routinely processed individually, infected materials were pooled for each farm, centrifuged at 500 × g for 10 min. and gentamicin sulfate solution (50 mg/ml) was added prior to inoculation into the allantoic cavity of five, 10-day-SPF-ECE (100 μl/egg). Inoculated embryos were incubated at 37°C for 24-48 hrs.
The allantoic fluid of pooled infected embryos from each sample (F7-F10) was subjected to haemagglutination in 25-μl volume in 96-well HI plates. Two fold serial dilution of the infected allantoic fluid was performed for each sample. Equal volumes of 1% suspension of chicken erythrocytes were dispensed to each well . Isolates possessed HA titers ≥ 3 log2 were suspected to be positive. The microtiter plate was left at room temperature until the negative control well exhibited a tight well circumscribed button of unagglutinated sedimented cells. The result of each sample was numerically reported.
Avian influenza virus antigen detection
The allantoic fluid that showed haemagglutination was screened for AIV group and H5 antigens using the rapid chromatographic strip test (Animal Genetic Inc. Korea) according to manufacturer instructions.
Viral RNA extraction and RT PCR
Viral RNA was extracted from virus containing the chorioallantoic membrane (CAM) homogenate of each sample using a SV Total RNA Isolation System (Promega Corporation, Madison, Wis.). RT-PCR amplification for the full length of both NA and HA genes of H5N1 were performed using the Verso™ 1 step RT PCR (Thermo Fisher Scientific Inc.). While a single set of primers flanking 1345 bp was used for NA, for HA, four sets of primers, flanking overlapping regions of the full length gene, were used (Table 1).
RT-PCR amplicons were subjected to electrophoresis in 1% agarose gel. HA specific bands were excised and purified with the QIAquick gel extraction kit (Qiagen). Each purified RT-PCR products was sequenced directly in both forward and reverse directions (Macrogen, Korea). Sequences were trimmed to remove amplicon primer-linker sequence and assembled. All sequence data used in this study are available in the GenBank database (Accession No: JF357720-JF357723).
Multiple sequence analysis and phylogenetic tree
Comparative analyses were performed using the CLUSTAL W multiple sequence alignment program, Mega 4.1 . AIV representative sequences used for the alignments were obtained from the GenBank and EMBL databases. The phylogenetic tree was constructed by using the neighbour-joining method with Kimura two-parameter distances by using the Mega 4.1. The reliability of internal branches was assessed by 1000 bootstrap replications and the p-distance substitution model. N-linked potential glycosylation sites were analyzed with the NetNglyc server with predicted threshold values above 0.5.
Deduced amino acid sequence analysis
We analyzed the HA deduced amino acid sequences of 4 isolated strains (F7-F10) and compared them with Egyptian H5N1 isolates available in the flu database. We used a complete collection of 315 H5N1 avian strains isolated from 2006 to 2010 in Egypt (2006 [n = 37], 2007 [n = 78], 2008 [n = 82], 2009 [n = 74], 2010[n = 44]), as well as 97 Egyptian human isolates, in order to screen the sites that had undergone amino acid intensive changes in recently determined different epitopes . Amino acid residues associated with mammalian virulence were also screened. Amino acid residues were examined in the flu database .
Building 3D model
The SWISS MODEL service was utilized to build the 3D models for the HA protein of F7-F10 isolates, by finding the exact templates. The PDB files were exposed by using the module of the Molsoft Internal Coordinate Mechanics (ICM-Pro) software.
Results and discussion
In the current study, four H5N1 strains were isolated from infected poultry broiler farms suffering from respiratory distress and variable mortalities; 2 from northern Egypt (Ismalia [F9] and Qaluobiya [F10]) and 2 from Middle Egypt (Fayoum Governorate [F7] and [F8]). The strains were isolated from pool organ samples that were inoculated in the SPF-ECE. All embryos died 24-48 hrs post inoculation. The infected allantoic fluid showed 3-4 haemagglutination log2 titers when tested with chicken RBC. The allantoic fluids reacted positively with chromatographic group and H5 antigen detecting strips. The RT-PCR for both HA and NA genes of H5N1 were found positive (data not shown). The present phylogenetic analysis, based on the analysis of the HA gene, showed that the cumulative genetic drifts in the HA gene clustered Egyptian isolates into two main lineages (A[A1-A5] and B[B1-B5]) (Figure 1). F7, A/chicken/Egypt/F7/2009 and F8, A/chicken/Egypt/F8/2009 were found to be related to lineage A1; F9 (A/chicken/Egypt/F9/2009) to lineage A2; and F10 (A/chicken/Egypt/F10/2009) to lineage B2 (Figure 1). Eight amino acid substitutions were found in the variant F10 (Figure 2) and other Egyptian variants in lineage B at the amino acid positions P74S, D 97N, H110R (epitope A), S123P, R140G (epitope B) and antigenic site), F144Y (epitope B), N165H (glycosylation site) and A184E.
The variability of amino acid in the 5 epitopes (A-E)  for the isolates under investigation and other available Egyptian sequences (2006-2010) in different years were screened. Residues, 43 (E epitope), 71 (E epitope), 110 (A epitope), 129 (A epitope and receptor binding site), 140 (B epitope), 141 (B epitope), 144 (B epitope), 151 (D epitope), 154 (D epitope, and glycosylation site), 185 (D epitope), 192 (D epitope) and 226 (D epitope) were subjected to high mutational pressure (Table 2).
F7 and F8 viruses and other viruses in the A1 group showed a conspicuous deletion of amino acid 129 in the HA protein (H5 numbering) (Table 2). The amino acid residue number 129 (S in F9 and L in F10) is located in the epitope A (Table 2) and part of the receptor binding site and also belongs to an antigenic site . The 129 deletion was not present in the ancestor A/Goose/Guangdong/1/96 strain  or in the H5N1 viruses originally introduced into Egypt in 2006. S129L substitution and 129 deletion were detected in 2007 (L in 6/78 isolates and deletion in 6/78 isolates), 2008 (L in 16/82, deletion in 0/82), 2009(deletion in 36/74, L 22/74) and 2010 (deletion in 26/44 and L in 18/44) (Table 2). Serine at site 129 is a receptor-binding site for α-2,3-linked sialic acids. Previously, it was reported that a mutation at this site (S129L) decreased virulence of HPAIV H5N1 in mice . However, a structural study on the H5N1 HA has reported that change from serine to leucine facilitates access to α-2,6-linked sialic acids . The 129S deletion was not detected in the 2006 human isolates (0/16), however, in 2007-2010, the 129 deletion was detected in (5/24) of 2007, (2/6) of 2008, (31/34) of 2009 and (12/12) 2010 human cases. Isolates from nonhuman mammals (equine, n = 1) like those in the majority of human in 2009 and 2010 possessed the 129-deletion. Interestingly, there was a similar deletion in the HA protein at the corresponding position of all human seasonal H1N1 and H3N2 viruses. It should be noted that this position is close to a domain modulating receptor interaction. Likewise, it is also interesting that strains with this deletion seem to evolve towards a receptor usage that is similar to human H1N1.
The deduced amino acid exchanges, as with most H5N1 Egyptian strains, are distributed across the HA protein and are found in the polybasic cleavage motif, whose consensus sequence is PGERRRKKR/GLF for clade 2.2 viruses. All the isolated strains possessed a multibasic HA0 amino acid cleavage sequence characteristic for the highly pathogenic H5N1 strains clade 2.2. While the sequence diversity of the haemagglutinin open reading frame of F7 and F8 had cleavage site of the consensus of 2.2, F10 PQGEGRRKKR/GLF showed (R325G) substitution and F9 showed two amino acid substitutions in the cleavage site PQGKSRRKKR/GLF, (E324K and R325S). All variations within the polybasic cleavage motif probably do not affect cleavage of the HA precursor molecule since F7-F10 have the RX(K/R)R consensus motif, for the cleavage by furin or the subtilisin-like protein convertases  which remains conserved. The significance of the change in the cleavage site, however, remains to be understood.
The HA sequences contained different N-linked potential glycosylation sites, F10 had eight potential sites at the following positions: 11, 23, 72, 236, 273, 286, 484 and 543. F7 and F8 possess 6 sites, at 11, 23, 164, 285, 483, 542 whereas F9 possessed only 5 sites at 23, 165, 286, 484, 543. F10 acquired two additional glycosylation sites (72 [epitope E] and 236 [polybasic cleavage site]) but lost a glycosylation site at165 that was located near the epitope D. F9 lost the glycosylation site at 11 position (H5 Numbering).
Lycett et al. specified 6 amino acid residues (86V, 124S, L/N138, T/S156, E/R212, T263) that are linked to the virulence of H5N1 in mammals . T156 and T263 were also present in F10. V86 was also present in 1/78 (2007), 1/82 (2008),5/74 (2009) isolates but none of the human isolates. T156 was present in 16/82, 23/74 and 17/44 of avian isolates in 2008, 2009 and 2010 isolates respectively, however it was present in 1/94 human isolates (A/Egypt/3300-NAMRU3/2008). T263 was present in F7-F10 isolates in addition to all human and 314/315 avian isolates (Table 3). On the other hand, while S129L and A156T substitutions were found to be associated with viral adaptation to poultry , both substitutions were found in F10 isolate.
The K189R amino acid substitution at the receptor-binding site was reported as an escape mutant change  and was further clearly shown to be significantly antigenic . Only, 4/315 isolates of the Egyptian GenBank sequences showed K189 while K189R was present in the F7-F10 isolates as well as 92/94 human and 296/315 avian Egyptian strains. Accordingly, it does not seem to be linked to escape mutants in the Egyptian strains. Duvvuri et al., analyzed potential mutational changes in L138Q, R140K and A156T . The amino acids at sites 138 and 140 were involved in the carbohydrate specificity of HA molecule towards host receptor . L138Q was present in the F7-F10, 80/82 of 2008 isolates, 72/74 of 2009 isolates and 42/44 of 2010 isolates and all human strains. On the other hand, R140K substation was present in only 8/78 of 2007 avian isolates but none of the human strains.
HA2 29E (H5 numbering 359) was detected in the F7-F10 isolates (Figure 2) and in all Egyptian human and avian isolates other than the A/chicken/Egypt/34-2/2008 and the A/chicken/Egypt/3046NAMRU3-CLEVB62/2007. All viruses in clade 1 and 2 possessed E29 in the HA2 subunit. The K29E substitution resulted in loss of a proteasomal cleavage site. The Glu-29 in the HA2 subunit may have been helpful to H5N1 viruses since such change induced the loss of one proteasomal cleavage site, which means that the virus can provoke a lower cytotoxic T cell response . All human and avian strains possessed E359 with the exception of two isolates; A/chicken/Egypt/3046NAMRU3-CLEVB62/2007 and A/chicken/Egypt/34-2/2008, which possessed K359 and N359 respectively.
Egypt is a good candidate of the study of AIV genetic evolution under high vaccination pressure as there is a high avian population density, significant inter-species and intra-species transmission and as the vaccination has been a regular veterinary practice since 2006. We identified several important differences in HA gene from the 4 H5N1 strains under investigation and the available Egyptian H5N1 sequences in the flu database and specified amino acids substitutions that: distinguished variant from classical parent strains, those linked to the virulence of H5N1 in mammals, and those present in different HA epitopes. The HA of the Egyptian H5N1 strains was found to undergo intensive genetic changes in a complex pattern. It has been speculated that the avian influenza vaccination has resulted in accelerated genetic drifts in the HA gene sequence . Nevertheless, other researchers have suggested that viral replication, per se, allows for the expression of drifts among subsequent viral populations [29–31]. The presence of mutant strain, that could replicate and shed in spite of vaccination and possesses amino acid substitutions linked to virulence in mammals, necessitates the importance of employing more rigorous requisites than those being utilized by the World Organization of Animal Health (OIE) which only considers protection afforded by the vaccine in terms of bird survival rates.
World Health Organization, World Organization for Animal Health, Food and Agriculture Organization: H5N1 Evolution Working Group. Toward a unified nomenclature system for highly pathogenic avian influenza virus (H5N1). Emerg Infect Dis 2008,14(7):e1. [http://www.ncbi.nlm.nih.gov/pubmed/18598616]
Salzberg SL, Kingsford C, Cattoli G, Spiro DJ, Janies DA, Aly MM, Brown IH, Couacy-Hymann E, De Mia GM, Dung do H, Guercio A, Joannis T, Maken AAS, Osmani A, Padalino I, Saad MD, Savić V, Sengamalay NA, Yingst S, Zaborsky J, Zorman-Rojs O, Ghedin E, Capua I: Genome analysis linking recent European and African influenza (H5N1) viruses. Emerg Infect Dis 2007, 13: 713-718.
Writing Committee of the Second World Health Organization: Consultation on Clinical Aspects of Human Infection with Avian Influenza A (H5N1) Virus: Update on Avian Influenza A (H5N1) Virus Infection in Humans. N Engl J Med 2008, 358: 261-273.
Abdel-Moneim AS, Shany SA, Fereidouni SR, Eid BT, El-Kady MF, Starick E, Harder T, Keil GM: Sequence diversity of the haemagglutinin open reading frame of recent highly pathogenic avian influenza H5N1 isolates from Egypt. Arch Virol 2009, 154: 1559-1562. 10.1007/s00705-009-0461-2
Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y: Evolution and ecology of influenza-A viruses. Microbiol Rev 1992, 56: 152-179.
Webster RG, Laver WG, Air GM, Schild GC: Molecular mechanisms of variation in influenza viruses. Nature 1982, 296: 115-121. 10.1038/296115a0
Wilson IA, Cox NJ: Structural basis of immune recognition of influenza virus hemagglutinin. Annu Rev Immunol 1990, 8: 737-771. 10.1146/annurev.iy.08.040190.003513
Stӧhr K: Influenza-WHO cares. Lancet Infect Dis 2002, 2: 517. 10.1016/S1473-3099(02)00366-3
Abdelwhab EM, Hafez HM: An overview of the epidemic of highly pathogenic H5N1 avian influenza virus in Egypt: epidemiology and control challenges. Epidemiol Infect 2011,139(5):647-657. [http://www.ncbi.nlm.nih.gov/pubmed/21281550] 10.1017/S0950268810003122
Ito T, Goto H, Yamamoto E, Tanaka H, Takeuchi M, Kuwayama M, Kawaoka Y, Otsuki K: Generation of a highly pathogenic avian influenza A virus from an avirulent field isolate by passaging in chickens. J Virol 2001, 75: 4439-4443. 10.1128/JVI.75.9.4439-4443.2001
Bragstad K, Nielsen L, Fomsgaard A: The evolution of human influenza A viruses from 1999 to 2006: a complete genome study. Virol J 2008, 5: 40. 10.1186/1743-422X-5-40
Beard CW: Serological procedures. In In American Association of Avian Pathologists. Iowa: Kendall/Hunt Publishing Company; 1989.
Abdel-Moneim AS, Abdel-Ghany AE, Shany SA: Isolation and characterization of highly pathogenic avian influenza virus subtype H5N1 from donkeys. J Biomed Sci 2010, 14: 25.
Kumar S, Tamura K, Jakobsen IB, Nei M: Molecular evolutionary genetics analysis software. Bioinformatics 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244
Gupta R, Jung E, Brunak S: Prediction of N-glycosylation sites in human proteins.2004. [http://www.cbs.dtu.dk/services/NetNGlyc/]
Duvvuri VRSK, Duvvuri B, Cuff WR, Wu GE, Wu J: Role of positive selection pressure on the evolution of H5N1 hemagglutinin. Genomics Proteomics and Bioinformatics 2009, 7: 47-56. 10.1016/S1672-0229(08)60032-7
Influenza Virus Resource2004. [http://www.ncbi.nlm.nih.gov/genomes/FLU/]
Xu X, Subbarao K, Cox NJ, Guo Y: Genetic characterization of the pathogenic influenza A/Goose/Guangdong/1/96 (H5N1) virus: similarity of its hemagglutinin gene to those of H5N1 viruses from the 1997 outbreaks in Hong Kong. Virology 1999, 261: 15-19. 10.1006/viro.1999.9820
Kaverin NV, Rudneva IA, Rudneva IA: Epitope mapping of the hemagglutinin molecule of a highly pathogenic H5N1 influenza virus by using monoclonal antibodies. J Virol 2007, 81: 12911-12917. 10.1128/JVI.01522-07
Li G, Tao S, Wang X: Sequence and epitope analysis of surface proteins of avian influenza H5N1 viruses from Asian patients. Chin Sci Bull 2006, 51: 2472-2481. 10.1007/s11434-006-2140-y
Veljkovic V, Veljkovic N, Muller CP, Müller S, Glisic S, Perovic V, Kӧhler H: Characterization of conserved properties of hemagglutinin of H5N1 and human influenza viruses: possible consequences for therapy and infection control. BMC Struct Biol 2009, 9: 21. 10.1186/1472-6807-9-21
Zhou A, Webb G, Zhu X, Steiner DF: Proteolytic processing in the secretory pathway. J Biol Chem 1999, 274: 20745-20748. 10.1074/jbc.274.30.20745
Lycett SJ, Ward MJ, Lewis FI, Poon AF, Kosakovsky PSL, Brown AJ: Detection of mammalian virulence determinants in highly pathogenic avian influenza H5N1 viruses: multivariate analysis of published data. J Virol 2009, 83: 9901-9910. 10.1128/JVI.00608-09
Perdue ML, Suarez DL: Structural features of the avian influenza virus hemagglutinin that influence virulence. Vet Microbiol 2000, 74: 77-86. 10.1016/S0378-1135(00)00168-1
Philpott M, Hioe C, Sheerar M, Hinshaw V: Hemagglutinin mutations related to attenuation altered cell tropism of a virulent avian influenza A virus. J Virol 1990, 64: 2941-2947.
Ilyushina NA, Rudneva IA, Gambaryan AS, Tuzikov AB, Bovin NV: Receptor specificity of H5 influenza virus escape mutants. Virus Res 2004, 100: 237-241. 10.1016/j.virusres.2003.12.032
Nielsen M, Lundegaard C, Lund O, Kesmir C: The role of the proteasome in generating cytotoxic T cell epitopes: insights obtained from improved predictions of proteasomal cleavage. Immunogenetics 2005, 57: 33-41. 10.1007/s00251-005-0781-7
Lee CW, Senne DA, Suarez DL: Effect of vaccine use in the evolution of Mexican lineage H5N2 avian influenza virus. J Virol 2004, 78: 8372-8381. 10.1128/JVI.78.15.8372-8381.2004
Ellis JS, Zambon MC: Combined PCR-heteroduplex mobility assay for detection and differentiation of influenza A viruses from different animal species. J Clin Microbiol 2001, 39: 4097-4102. 10.1128/JCM.39.11.4097-4102.2001
Gambaryan A, Tuzikov A, Pazynina G, Bovin N, Balish A, Klimov A: Evolution of the receptor binding phenotype of influenza A (H5) viruses. Virology 2006, 344: 432-438. 10.1016/j.virol.2005.08.035
Widjaja L, Ilyushina N, Webster RG, Webby RJ: Molecular changes associated with adaptation of human influenza A virus in embryonated chicken eggs. Virology 2006, 350: 137-145. 10.1016/j.virol.2006.02.020
Wu WL, Chen Y, Wang P, Song W, Lau SY, Rayner JM, Smith GJD, Webster RG, Peiris JSM, Lin T, Xia N, Guan Y, Chen H: Antigenic profile of avian H5N1 viruses in Asia from 2002 to 2007. J Virol 2008, 82: 1798-1807. 10.1128/JVI.02256-07
Bright RA, Ross TM, Subbarao K, Robinson HL, Katz JM: Impact of glycosylation on the immunogenicity of a DNA-based influenza H5 HA vaccine. Virology 2003, 308: 270-278. 10.1016/S0042-6822(03)00008-4
Chen H, Bright RA, Subbarao K, Smith C, Cox NJ, Katz JM, Matsuoka Y: Polygenic virulence factors involved in pathogenesis of 1997 Hong Kong H5N1 influenza viruses in mice. Virus Res 2007, 128: 159-163. 10.1016/j.virusres.2007.04.017
The authors declare that they have no competing interests.
ASA carried out all the experiments including designing the experiments, acquisition of data, analysis and interpretation of data. and drafted the manuscript. MAA and MFE and have helped in acquisition of samples, data analysis and revising the manuscript. All authors have read and approved the final manuscript.