Previous Article | Next Article ![]()
Journal of Virology, August 2004, p. 8372-8381, Vol. 78, No. 15
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.15.8372-8381.2004
Southeast Poultry Research Laboratory, Agricultural Research Service, U.S. Department of Agriculture, Athens, Georgia 30605,1 National Veterinary Services Laboratories, Animal and Plant Health Inspection Service, U.S. Department of Agriculture, Ames, Iowa 500102
Received 3 January 2004/ Accepted 16 March 2004
|
|
|---|
|
|
|---|
It is well documented that human influenza virus undergoes frequent antigenic drift, which is the accumulation of point mutations in the antigenic domain of the HA protein (3). As a result, viruses with a slightly changed antigenic structure emerge and can escape the host's acquired immunity, whether this immunity is acquired by natural infection or vaccination. Therefore, to maintain optimal protection by vaccination, the presently prevailing strains of influenza virus need to be included in each year's influenza vaccine, requiring yearly reevaluation and frequent changes to the vaccine formulation (1, 6). In contrast, broad subtype-specific immunity to HPAI has been observed following parenteral AI vaccination in chickens (30, 31). In the United States, H7N2 LPAI viruses have been circulating in live-bird markets since 1994 and progressive genetic drift in the HA gene has also been observed (23, 25). However, an existing commercial AI vaccine prepared from a 1997 seed stock H7N2 virus was able to provide protection against an H7N2 virus isolated in 2002 from a turkey in Virginia (35). In another study, which used an HPAI challenge, vaccine strains with HA protein sequences showing as little as 86% similarity to that of the challenge strain still provided good protection in chickens (29). One possible difference between influenza viruses infecting humans and poultry might relate to the different life spans of the hosts. Humans have multiple exposures to influenza viruses, either through natural infection or vaccination, and have a level of population immunity, but because of their short productive lives, poultry are usually naïve to influenza virus, particularly since vaccination for influenza remains uncommon. This likely results in different types of selection pressure on the virus in each host and may partially explain the broader protection provided by AI vaccine in chickens. In this context, long-term use of AI vaccine in Mexico and adjacent countries presents us with a unique opportunity to examine the evolution of AI virus in the presence of vaccine pressure.
The first suspected cases of AI in poultry in Mexico were detected in late 1993, and LPAI (H5N2) virus was first identified in May 1994. By the time influenza was confirmed in Mexican poultry, it had already spread widely, and eventually the LPAI virus mutated to the HPAI (H5N2) form of the virus in two different regions in Mexico, Puebla and Queretaro (36). The Mexican government started using vaccination in 1995 to aid in the control of both HPAI and LPAI viruses. Over 1 billion doses of inactivated emulsified vaccine have been authorized for use from 1995 to 2001, and 459 million doses of recombinant fowl pox-vectored AI vaccine have been authorized from 1998 to 2001 (37). Though HPAI viruses have not been reported since 1996, LPAI virus continues to circulate in Mexico, and genetically related viruses have also circulated in Guatemala and El Salvador since at least 2000 and 2001, respectively (27). Though some of the earlier Mexican H5N2 isolates have been characterized previously (7, 8, 14), there has been no further characterization of field isolates after the introduction of the vaccine. Thus, it was speculated that antigenic drift might have occurred in this rapidly changing RNA virus. In the present study, we gathered 52 Mexican lineage H5N2 viruses isolated from 1993 through 2002 and conducted molecular and antigenic analysis. Our study demonstrates that the Mexican lineage H5N2 viruses have undergone antigenic drift away from the vaccine strain and highlights the effect of vaccine use on the evolution of AI virus.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Mexican lineage H5N2 AI virus isolates analyzed in this study
|
RNA extraction and RT-PCR. Using an RNeasy mini kit (QIAGEN, Valencia, Calif.), viral RNA from infectious allantoic fluid from ECE was extracted according to a modified protocol as previously described (22). Standard reverse transcription (RT) PCR was carried out with a QIAGEN one-step RT-PCR kit with primers specific for influenza virus. The primer sequences and amplification conditions are available from the authors upon request. The PCR product was separated on an agarose gel by electrophoresis, and amplicons of the appropriate size were subsequently excised from the gel and extracted with a QIAGEN gel extraction kit.
Sequencing and phylogenetic analysis. Direct sequencing was performed with a PRISM Ready Reaction DyeDeoxy Terminator Cycle Sequencing kit (Perkin-Elmer, Foster City, Calif.) run on a 3700 or 3730 automated sequencer (Perkin-Elmer). The nucleotide sequences were compared initially by using the Megalign program (DNASTAR, Madison, Wis.) with the Clustal V alignment algorithm. Pairwise sequence alignments were also performed in the Megalign program to determine nucleotide and amino acid sequence similarities. Phylogenetic comparisons of the aligned sequences for each gene segment were generated using the maximum parsimony method in a heuristic search using PAUP 4.0b10 software (Sinauer Associates, Inc., Sunderland, Mass.) (33).
Virus titration. The infectious doses of viruses used for the challenge experiment and for vaccine preparation were determined using ECE as previously described (38). Tracheal swab samples obtained from vaccine experiments were titrated by quantitative real-time RT-PCR assay. In brief, swabs from individual birds were suspended in 1.5 ml of brain heart infusion medium (Difco, Detroit, Mich.). RNA was extracted from fluid containing tracheal swabs by using an RNeasy mini kit (QIAGEN), and quantitative real-time RT-PCR was performed with primers for type A influenza virus as described previously (16a). The 50% egg infective doses (EID50) of virus from the clinical samples were interpolated from the cycle thresholds by using standard curves generated from known amounts of control RNA (101.0 to 106.0 EID50/ml). The details of the procedure and the validation of the assay have been described previously (22).
Production of hyperimmune antisera. Seven field isolates (Table 2) were selected based on the topology of the HA1 phylogenetic tree (Fig. 1), and antisera were prepared against these isolates and also from a commercial vaccine. The seven field isolates were grown in 10-day-old ECE, and the infectious allantoic fluid was pooled for each virus. Infective titers were determined prior to inactivation with 0.1% betapropiolactone (Sigma, St. Louis, Mo.), and the inactivated virus was used to produce an oil emulsion vaccine as previously described (24). Three-week-old specific-pathogen-free (SPF) chickens were inoculated subcutaneously with either the commercial vaccine or the inactivated oil emulsion vaccines made at SEPRL. The vaccine dose for the seven field viruses contained 107.0 to 108.0 EID50 per 0.5 ml before viral inactivation. Three weeks later, the chickens were booster vaccinated by the same route with the same amount of vaccine. Sera were harvested and inactivated at 56°C for 30 min before being used in the HI and VN procedures.
|
View this table: [in a new window] |
TABLE 2. Genetic similarity of the Mexican lineage isolates
|
![]() View larger version (52K): [in a new window] |
FIG. 1. Phylogenetic tree based on nucleotide sequences of the HA1 and NS genes from Mexican lineage isolates. The tree, generated by the maximum parsimony method with PAUP 4.0b10, is the result of heuristic search and midpoint rooting.
|
The VN test was performed in ECE by using the diluted-serum constant-virus (beta) procedure (9). Briefly, antiserum was mixed with virus for 1 h at room temperature. This mixture was then inoculated into ECE. Three days postinoculation (dpi), allantoic fluid was examined for hemagglutinating activity to determine the presence of the virus. The Reed and Muench formula was used to calculate endpoint titers for homologous and heterologous neutralization (21). Antigenic relatedness values were then calculated using the Archetti and Horsfall formula (2). The isolates and the antisera used in the HI and VN tests are listed in Table 3.
|
View this table: [in a new window] |
TABLE 3. Cross-HI and cross-VN test results expressed as percentages of antigenic relationshipa
|
|
View this table: [in a new window] |
TABLE 4. Vaccine protection against LPAI challenge
|
|
View this table: [in a new window] |
TABLE 5. Cross-protection against HPAI virus (CK/Puebla/8623/607/94) challenge
|
GenBank accession numbers. The nucleotide sequences presented in this article have been deposited in the GenBank database under accession numbers AY497063 to AY497191.
|
|
|---|
Genetic analysis of H5 HA1. The coding sequences for the HA1 gene segments from 34 isolates were newly determined. These include six isolates from Guatemala and two isolates from El Salvador. A total of 52 Mexican lineage viruses that were isolated between 1993 and 2002 were compared. The percentages of nucleotide and amino acid similarities among all HA1 sequences ranged from 86.0 to 99.8% and 84.6 to 99.7%, respectively. According to a phylogenetic tree based on the HA1 sequence, at least four distinct clusters of viruses were observed within the Mexican lineage isolates (Fig. 1). Viruses isolated before 1996 were mostly divided into the previously described Jalisco and Puebla sublineages (7), except one isolate, A/CK/Puebla/28159-474/95. The majority of the early isolates belonged to the Jalisco sublineage. However, several clusters of isolates were found within this Jalisco sublineage. Nine viruses isolated in Mexico after 1995 formed a separate clade from the Puebla and Jalisco sublineages and were designated as sublineage A. Viruses isolated in Guatemala and El Salvador were closely related to each other and were part of a fourth Mexican clade of viruses, including the early isolate A/CK/Puebla/28159-474/95, and were designated as sublineage B.
Since phylogenetic analysis demonstrated heterogeneity in the HA gene in the early isolates, we made the consensus sequence of the HA1 gene with 18 viruses isolated before 1995 (Table 1). The sequence of the vaccine strain showed only two amino acid differences from the consensus sequence (Fig. 2). To determine the mutations that may have occurred after the introduction of the vaccine, we compared the sequence of the vaccine strain with sequences from isolates belonging to two recent sublineages (Fig. 2). A total of 20 and 27 amino acid substitution sites, which indicated that mutations had become fixed in a progressive manner from earlier isolates to recent isolates, were detected in sublineages A and B, respectively. Fifteen substitution sites overlapped between the two sublineages. Six substitutions (at positions 136, 137, 154, 181, 188, and 275) occurred in regions previously proposed as antigenic sites A, B, and C in H1 and H3 structures (5, 40). Two amino acid changes (at positions 126 and 136) occurred in the site that was previously demonstrated in an H5 escape mutant (15, 19, 20). The amino acid change from D to N at position 126 resulted in the formation of a glycosylation site that was also observed in an H5 escape mutant. Though many substitutions occurred in regions or sites other than the previously proposed antigenic sites, the majority of changes were located in the distal globular head of HA1. Furthermore, most of the overlapping substitution sites between A and B sublineage viruses were located in the outer surface of the globular head of HA1. The locations of substitution sites are illustrated in the three-dimensional structure of the H5 HA (11) (Fig. 3).
![]() View larger version (60K): [in a new window] |
FIG. 2. Alignment of HA1 amino acid sequences of Mexican lineage isolates with that of the consensus sequence made with 18 early isolates. The underlined residues are potential glycosylation sites, and the residues in the open boxes are previously identified antigenic sites with H1 and H3 molecules. , substitution sites demonstrated in an H5 escape mutant; , amino acid substitution sites in sublineage A and B viruses analyzed in this study. Overlapping sites are circled.
|
![]() View larger version (65K): [in a new window] |
FIG. 3. Ribbon diagram of the monomer of H5 (A/Duck/Singapore/3/97) HA. Front (a) and back (b) views are shown. Location of amino acid changes in HA1 and potential glycosylation sites of Mexican lineage isolates are labeled. The color scheme is as follows: green, amino acid changes in sublineage A virus; blue, amino acid changes in sublineage B virus; red, overlapping amino acid changes in sublineage A and B viruses; pink, potential glycosylation sites; purple, termini of HA1 (N1 and C1) and HA2 (N2 and C2).
|
By comparing the consensus sequence with sequences of consecutive isolates in both sublineages, the rates of fixation of mutations were calculated to be 12 x 103 and 10 x 103 substitutions per nucleotide site per year for sublineages A and B, respectively. Evolutionary rates based on amino acid changes were calculated to be 4.2 x 103 and 4.0 x 103 substitutions per nucleotide site per year, respectively.
Genetic analysis of NS and remaining genes. The coding sequence for the NS gene segments from 47 isolates were newly determined and a total of 52 Mexican lineage viruses were phylogenetically analyzed. All of the Mexican lineage isolates belonged to subtype B (group B) and assorted into a similar topology, as we observed in the HA tree. The substitution rates for sublineages A and B were calculated to be 6.0 x 103 and 12 x 103 substitutions per nucleotide site per year, respectively.
The remaining genes (N2, M, NP, PA, PB1, and PB2) of selected isolates from different sublineages were also sequenced and compared with other sequences available in GenBank. All Mexican lineage isolates sequenced were closely related (Table 2) to each other and no evidence of reassortment was observed. All of the isolates examined had stalk deletions in the N2 gene, which is thought to be a characteristic of chicken-adapted AI viruses (17).
One-way serologic test of vaccine antisera to the field isolates. To determine the cross-reactivity of vaccine antisera to the viruses from different sublineages, the HI test was conducted with sera obtained from chickens vaccinated with commercial vaccine (Table 1). The results of the HI assay demonstrated that there is high cross-reactivity of vaccine antisera to the viruses belonging to the Jalisco sublineage, which was the source of the seed strain for the vaccine. Compared to the HI titer of the sera to the homologous vaccine strain, the HI titers against most of the Jalisco sublineage isolates were similar or within a fourfold difference. The HPAI viruses belonging to the Puebla sublineage demonstrated a fourfold loss in HI activity. In contrast, more than a fourfold difference was observed with several sublineage A viruses, and the sublineage B viruses had as much as a 16-fold difference in HI activity compared to the vaccine strain.
Cross-HI and -VN tests. To further determine the antigenic relatedness among isolates belonging to different sublineages, eight isolates were selected on the basis of phylogenetic and one-way serologic data. We then conducted cross-HI tests with hyperimmune sera produced in chickens against representative isolates of the different sublineages. The percentage of antigenic relatedness (R value) was then determined with the Archetti and Horsfall method (2) by calculation of homologous and heterologous titer ratios (Table 3). In general, serologically related isolates have R values of more than 50%. The vaccine strain shared more than 50% relatedness only with isolates belonging to the Jalisco and Puebla sublineages. Isolate A/CK/Puebla/28159-474/95, which showed high cross-reactivity with vaccine sera in the one-way HI test, demonstrated only 32% relatedness to the vaccine strain in the cross-HI test. Isolates from Guatemala and El Salvador were antigenically close to each other (71%) but shared very low antigenic relatedness to any of the other isolates examined (less than 32%). To verify these results, we conducted VN tests with five isolates used in the HI test and found results similar to those of the HI tests (Table 3). In general, the VN test was more discriminative than the HI test and hyperimmune sera from vaccine-immunized chickens did not effectively neutralize the sublineage A and B viruses. Therefore, it can be concluded that Jalisco and Puebla sublineage viruses are antigenically similar to the vaccine strain, whereas sublineage A and B viruses are antigenically different from the vaccine strain and each other.
Efficacy of the commercial vaccine against LPAI virus challenge. To compare the efficacies of the present vaccine against antigenically different viruses, birds were challenged with two viruses from recent sublineages after a single vaccination with commercial vaccine. The vaccine strain was also used for a challenge comparison. Three groups of unvaccinated birds were also challenged with each virus for comparison. The efficacies of the vaccine were compared in terms of the ability to prevent virus replication in the trachea, and the results are summarized in Table 4. In birds challenged with the vaccine strain, large amounts of virus were detected only in the birds that showed low HI titers (1:16 in twofold serial dilution). In most of the birds that had HI titers greater than 32, the vaccine prevented detectable virus shedding, with only a few birds shedding small amounts of virus that were about 2 logs lower than virus detected from unvaccinated and challenged control birds. In contrast, large amounts of virus, comparable to the amount of virus detected in unvaccinated challenged birds, were detected in birds challenged with two recent viruses and the amount of virus shedding was independent of the HI antibody titer of the birds.
Cross-protection against HPAI virus challenge. To evaluate the effect of the antigenic difference between the vaccine and challenge strains in the protection of chickens against HPAI challenge, three groups of birds were vaccinated with strains that show various degrees of antigenic relatedness to the challenge virus (Table 5). Among birds vaccinated with a homologous strain as the challenge virus, only 5 of 16 shed low amounts of virus (mean titer of 101.45 EID50/ml) at 3 dpi. In contrast, all birds vaccinated with commercial vaccine and most birds (14 of 16) vaccinated with a Guatemala isolate shed higher amounts of virus (mean titers of 103.03 and 103.32 EID50/ml, respectively) from the trachea. The amount of viral shed was approximately 1 log lower than that observed in the unvaccinated and challenged control birds. However, no clinical signs were observed in any of the vaccinated birds and the birds remained healthy during the 10-day observation period. All unvaccinated and challenged control birds demonstrated typical HPAI clinical signs, such as edema of the head and face, cyanosis of unfeathered skin, etc., and seven of eight birds died, with a mean time to death of 6.6 days.
|
|
|---|
In one-way HI tests, we observed high cross-reactivity of vaccine antisera to viruses belonging to the Jalisco sublineage, which includes the vaccine seed strain (Table 1). In contrast, differences of more than fourfold were observed between sublineage A and B viruses. This antigenic difference was further confirmed in cross-HI and cross-VN tests, in which sublineage A and B viruses showed less than 50% antigenic relatedness to the vaccine strain (Table 3). In contrast, HP virus from the Puebla lineage was about 50% related to the vaccine strain and this may partially explain the early disappearance of Puebla lineage HP virus from Mexico after vaccination was introduced. From the above mentioned genetic and serologic findings, it is likely that sublineage A viruses are undergoing antigenic drift away from the vaccine strain, though additional analysis with more recent isolates should be conducted to confirm this conclusion. The sublineage B viruses, possibly originated from A/CK/Puebla/28159-474/95, may have a selective advantage due to their antigenic differences from the vaccine strain and are also drifting further away from the vaccine strain.
In general, in vitro serologic tests assessing similarity are more discriminating than in vivo cross-protection studies (10, 13). Thus, we conducted a vaccine protection test to further determine the antigenic relationship in vivo as well as to assess the efficacies of the presently available commercial vaccine against two recent viruses. The in vivo test further confirmed the serologic data indicating that the commercial vaccine tested was not able to prevent virus shedding when chickens were challenged with antigenically different field isolates (Table 4). Though clinical signs were not observed with these LP viruses in experimental conditions, the shedding of H5 subtype LP viruses is of concern because of the potential for transmission to naïve flocks and the possibility of mutation to the HP form of the virus (7, 14, 16, 18).
Previous studies have demonstrated broad cross-protection against HPAI virus in terms of morbidity and mortality (30, 31). Our study confirms the previous finding that chickens can be protected from clinical disease resulting from HPAI challenge even with a high degree of antigenic difference (Table 5). However, the level of virus shedding in the trachea correlated with the antigenic differences of vaccine and challenge strains. This discrepancy likely reflects the fact that vaccines for AI do not prevent infection but that antibodies in the bloodstream may effectively prevent the systemic phase of disease caused by HPAI viruses.
The H1 and H3 subtype viruses from humans comprise a variety of HA antigenic drift variants, and the antigenic epitopes have been well defined within the H3 molecule's three-dimensional structure (5, 40). In contrast, avian HA subtypes retain a relatively stable antigenic structure, and selection and characterization of monoclonal antibody escape mutants has been used as the only means of identifying antigenic epitopes (15, 19, 20). In this study, we mapped the amino acid changes on the three-dimensional structure of the H5 HA molecule (11) and most of the amino acid substitutions were found in the membrane distal globular head of HA1 (Fig. 3). Between sublineages A and B, we found 15 overlapping substitution sites that were located in the outer surface of the globular head of HA1. Compared to what has been reported previously, six substitutions were found to occur in regions previously proposed as antigenic sites A, B, and C in H1 and H3 structures and two amino acid changes occurred in sites that were previously demonstrated with H5 monoclonal antibody escape mutants (Fig. 2). Though several substitutions occurred in sites other than that previously identified with H5 escape mutants, it is possible that the choice of monoclonal antibodies could have limited the diversity of the selected escape mutants in previous studies. The amino acid change of aspartic acid to asparagine at position 126 that resulted in the formation of a glycosylation site is of interest because glycosylation of viral antigens can mask antigenic epitopes and therefore is an important process in the generation of new viruses. Interestingly, this acquisition of glycosylation was also observed in an H5 escape mutant and in H2 subtype escape mutants at equivalent positions (15, 34). It is plausible that the substitution sites described here may play an important role in drift variation in the field with constant vaccine pressure, and characterization of other H5 subtype antigenic variants will further elucidate the importance of individual substitution sites.
In this study, we demonstrate that presently prevailing Mexican lineage AI viruses, which belong to two different sublineages, are antigenically different from the vaccine strain. Although circumstantial evidence supports the idea that vaccine pressure is driving this phenomenon, other possibilities for explaining the antigenic drift cannot be ruled out. Whatever the cause, the antigenic drift of H5N2 Mexican lineage viruses away from the vaccine strain is of concern for the effectiveness of vaccination strategies within the poultry industry. Vaccination is used not only to prevent clinical disease but also to reduce viral shedding. A reduction in viral shedding reduces the chance of virus spreading from infected vaccinated flocks to uninfected flocks. This is the principle behind barrier or "ring" vaccination around a quarantine zone during an outbreak, which can be a valuable control tool. Vaccination, however, cannot be used alone for the control of AI and must be accompanied by other control measures, including quarantines, controlled depopulation (taking birds to slaughter facilities under tight controls), and increased surveillance. The almost total failure of the present vaccine strain to prevent viral shedding when birds were challenged with recent field isolates was expected because of the extreme difference in cross-HI and -VN activities and greatly reduces the value of vaccination as part of a control and eradication program. Our in vitro and in vivo findings indicate that to maintain optimal protection with our AI vaccines, we should evaluate how well the prevailing isolates match the present AI vaccine formulation and update the vaccine by using criteria similar to those used for human influenza vaccines (1, 6)
We thank Suzanne Deblois, Joan Beck, and the SAA sequencing facility for technical assistance and Roger Brock for animal care assistance. We also thank Erica Spackman, David Swayne, and Terrence Tumpey for critical reviews of the manuscript.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»