Previous Article | Next Article ![]()
Journal of Virology, March 2008, p. 2295-2304, Vol. 82, No. 5
0022-538X/08/$08.00+0 doi:10.1128/JVI.02267-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

,#
Raúl Rabadán,3,
Arnold J. Levine,3 and
Peter Palese1,2*
Departments of Microbiology,1 Medicine, Mount Sinai School of Medicine, 1 Gustave L. Levy Place, New York, New York 10029,2 Institute for Advanced Study, Einstein Dr., Princeton, New Jersey 085403
Received 18 October 2007/ Accepted 7 December 2007
|
|
|---|
|
|
|---|
Two simple models have been hypothesized for the packaging of influenza vRNA: random incorporation and selective incorporation (20). The random-incorporation model assumes that a common structural feature which enables them to be randomly incorporated into budding virions is present on all vRNAs. Assuming equal concentrations of vRNA segments in the cytoplasm, the probability that eight different segments are packed in one virion is very small (P value = 8!/88 = 0.0024), indicating that, if this is the case, only a few virus particles will be viable. If instead of eight segments we consider that more than eight are incorporated, the probability that eight of them are different goes to more-reasonable numbers (1, 7). If 12 vRNAs were packaged randomly, the mathematical models suggest that infectivity increases to approximately 10%, which is comparable with the percentage for the experimental data (4). As only 1 to 2% of the weight of the influenza virus particle is vRNA, it is difficult to accurately quantify the exact number of vRNA segments packaged.
The selective-incorporation model suggests that each vRNA segment contains a unique "packaging signal" allowing it to act independently and be selectively packaged. The use of green fluorescent protein (GFP) packaging construct incorporation signals has been described for all eight segments (9, 10, 12, 15, 17, 25). However, it is still not understood how these signals direct the packaging of vRNA into budding virions. Liang et al. (12) demonstrated that both PB1 and PA vRNA needed at least 40 5' and 66 3' coding nucleotides in order to package efficiently and that PB2 vRNA required only 80 5' coding nucleotides, with packaging being independent of coding nucleotides at the 3' end of the vRNA. It was also demonstrated that the packaging signals from individual vRNA segments work in conjunction. Pairing of a 3' signal from one segment and a 5' signal from a different segment did not direct the packaging of a GFP-based vRNA into budding virions. The result for the PB2 vRNA is consistent with previously reported data for PB2 defective interfering vRNA (6). Dos Santos Afonso et al. (5) also demonstrated the importance of the 5' end of the PB2 vRNA and the necessity of it being adjacent to the 5' untranslated region. This was demonstrated by the insertion of a flag epitope tag fused to the PB2 open reading frame. Placement of the tag at the 5' end of the vRNA required duplication of the last 109 coding nucleotides of PB2, whereas the tag sequence could be inserted before the start codon at the 3' end of the vRNA with no effect on vRNA packaging.
Using an approach similar to that of Liang et al. (12) with GFP packaging constructs, Muramoto et al. (15) also analyzed the requirements of the 3' and 5' coding regions as packaging signals for the polymerase vRNAs. Although their results differ slightly from those of the previous study, this was explained by the different approaches taken in the two studies, with the study by Liang et al. competing the packaging of GFP constructs with the packaging of the parental vRNA. In contrast to that of Liang et al. (12), the approach of Muramoto et al. (15) required 120 nucleotides at the 3' end of the PB2 vRNA for efficient packaging, which was inconsistent with previous data suggesting that this end is not important for packaging. The Muramoto et al. study also suggested a hierarchy of vRNA incorporation, with data showing that the loss of incorporation of one segment leads to reduced incorporation of other segments. This effect was most prominent when the PB2 vRNA was omitted (15). This hierarchy and the intersegment interactions are supported by thin-section electron micrographs showing viral ribonucleoproteins organized in a distinct pattern within virions (16).
To summarize, increasing evidence from many different laboratories that supports the selective-incorporation model has been accumulated. First, there are vRNA regions within the coding sequence at both the 5' and 3' ends of almost every segment that are sufficient for packaging foreign RNA sequences. Second, loss of incorporation of one segment alters the incorporation of the others, and finally, electron micrographs have shown that there is a clear segmental pattern within a virion.
If the packaging regions are crucial for the viability of the virus, we expect them to be conserved. Factorizing out population structure and amino acid conservation, we have been able to find some conserved regions at the ends of each of the different segments. Possible explanations for the presence of these conserved regions in the vRNA have to do with RNA structure, RNA-RNA interactions, and RNA-protein interactions. We are able to show that these regions of highly conserved sequence are important for packaging of vRNA. Mutations in key regions of these highly conserved sequences resulted in 5- to 20-fold reductions in the levels of packaging of the mutated vRNAs compared with the level of wild-type (wt) vRNA packaging. We are also able to confirm the results of Muramoto et al. (15) demonstrating that the reduction in packaging of one vRNA results in reductions in other vRNAs.
While this study was being completed, Gog et al. (11) reported a result for the PB2 segment similar to ours. In that study, synonymous changes were introduced into a GFP packaging construct and the effect on packaging of this construct was analyzed. Changes in individual conserved residues resulted in a 20- to 50-fold reduction in packaging. In contrast to the experimental results reported herein, the packaging of the GFP construct was in competition with that of the parental PB2 vRNA, which was provided by infection with wt virus.
|
|
|---|
Analysis of influenza virus sequences. To analyze the genomes of influenza A viruses for conserved regions within the polymerase vRNA, approximately 600 avian sequences for each segment were aligned (667, 591, and 621 sequences for segments 1, 2 and 3, respectively). Sequences were obtained from the Influenza Virus Resource, National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/genomes/FLU/). Avian sequences were chosen because there is more sequence variation within the published sequences and to avoid the bias within the human sequences in the public database. The bias arises in part from the fact that a large percentage of sequences in the database have been isolated in New York and New Zealand within the past 10 years and from the fact that the human viruses present a more severe bottleneck structure than the avian viruses.
Constructs and cloning. The plasmids used for the rescue of recombinant influenza A/WSN/33 (WSN) virus have been described previously (8). To generate recombinant WSN viruses with mutated PB1, PB2, or PA vRNA segments, pPOLI plasmids encoding the WSN gene of interest were subjected to site-directed mutagenesis using a Stratagene QuikChange XL site-directed mutagenesis kit (Stratagene, La Jolla, CA). For the generation of recombinant A/Puerto Rico/8/34 (PR/8) viruses with mutated PA vRNA segments, we utilized a Stratagene QuikChange XL site-directed mutagenesis kit (Stratagene) for mutagenesis of a pGEM plasmid containing the PA gene of PR/8. Sequences of each mutated construct were confirmed by automated sequencing. Correctly mutated PR/8 PA genes were then subcloned into pDZ (21). All primer sequences and vector maps are available upon request.
Reverse genetics for recombinant viruses. The reverse-genetics systems for the generation of recombinant influenza WSN (8) and PR/8 (23) viruses have been described previously. For the generation of recombinant WSN viruses, 293T and MDCK cells were cotransfected with eight pPOLI vectors expressing negative-sense vRNA and four pCAGGS protein expression vectors encoding the viral polymerase subunits and NP. After 48 h of incubation, recombinant viruses were passaged to fresh MDCK cells for amplification. For the generation of recombinant PR/8 viruses, 293T cells were transfected with eight pDZ plasmids, which utilize bidirectional transcription to both encode the viral genomic RNA and express viral protein. Cells were collected 24 h posttransfection and inoculated into the allantoic cavities of 10-day-old embryonated chicken eggs. Rescue of recombinant viruses was assessed by hemagglutination activity. Each of the newly generated viruses was further plaque purified, and mutations were confirmed by sequencing of mutated genes.
Viral growth kinetics. The growth of mutant recombinant WSN viruses was evaluated by inoculation of MDCK cells at a multiplicity of infection (MOI) of 0.001. Cells were incubated at 37°C for 48 h, at which time the virus titer in the supernatant was determined by plaque assay of MDCK cells. Growth kinetics of PR/8 viruses was determined by inoculation of 10-day-old eggs with 100 PFU of virus. At 48 h postinoculation, the virus titer in the allantoic fluid was determined by titration of plaques on MDCK cells.
Isolation of packaged vRNAs. Packaged vRNA for WSN recombinants was analyzed by purification of virus from the supernatants of three 150-mm dishes of MDCK cells. Cells were inoculated at an MOI of 0.001. Cells were incubated for 36 to 48 h until maximal cytopathic effect was visible. Supernatant was clarified by low-speed centrifugation and then further clarified by centrifugation at 10,000 rpm using a Beckman SW20 rotor (Beckman Coulter, Fullerton, CA). Clarified supernatant was then layered on a 30% sucrose cushion and further centrifuged at 25,000 rpm for 2.5 h. Pelleted virus was resuspended in TMK (10 mM Tris-HCl, pH 7.5, 1.5 mM MgCl2, 10 mM KCl) and vRNA extracted using TRIzol LS reagent (Invitrogen). Precipitated vRNA was resuspended in a final volume of 15 µl of 10 mM Tris-HCl (pH 8.0) and stored at –80°C. To analyze packaged vRNA for PR/8 mutant viruses, 10-day-old eggs were inoculated with approximately 1,000 PFU and incubated for 2 days. Allantoic fluid was harvested and processed as described above.
qPCR analysis of packaged vRNAs. Extracted vRNAs (approximately 200 ng) were reverse transcribed using a universal 3' primer (5'-AGGGCTCTTCGGCCAGCRAAAGCAGG) and Superscript II reverse transcriptase (RT) (Invitrogen). The RT product was then diluted 10,000-fold and used as a template for quantitative PCR (qPCR). Separate PCRs were then carried out as previously reported (14) with segment-specific primers utilizing a LightCyler 480 apparatus (Roche, Nutley, NJ). The 10-µl reaction mixture contained 1 µl of diluted RT product, SYBR green I (Molecular Probes), a 0.5 µM concentration of each primer, a 200 µM concentration of each of the deoxynucleoside triphosphates, 3 mM MgCl2, 2 µl of 10x PCR buffer II, and 1 U of AmpliTaq Gold enzyme (Applied Biosystems, Foster City, CA). Relative concentrations of vRNA were determined on the basis of an analysis of cycle threshold values (14), normalization for the total vRNA amount by equalization of the level of PB2 vRNA, and then calculation of the percentages of incorporation relative to the levels of wt vRNA packaging. Results are presented as the average levels of incorporation of vRNA ± standard deviations, with results derived from two independent virus purifications, with vRNA levels quantified in triplicate.
|
|
|---|
To determine whether sequence conservation was at the RNA level, we have to factor out amino acid conservation and the population structure. Viral sequences are highly related to each other due to the evolution of the virus and the sampling bias. This correlation in the sequences could produce false positives if a few sequences are overweighed due to the population of the virus or a sampling bias.
How can we factor out amino acid conservation and population structure? The first step we followed was to consider the conservation in third-codon positions, eliminating most of the conservation due to the protein structure. If the conservation is due to population structure (null hypothesis), we would expect the conserved positions to be independent from each other, i.e., for them to be randomly distributed along the genome of the virus and not clustered together. Fixing a conservation cutoff, c (in our analysis, we used a conservation cutoff of 98%; i.e., in at least 98% of the particular nucleotides, the position is conserved), we computed the probability that we would get by chance a position that is more conserved than the cutoff value. This probability is just the number of positions that are more conserved than the cutoff value, n(c), over the total number of positions, n(0). The quotient of these two quantities [Pc = n(c)/n(0)] defines the probability of finding a particular well-conserved position, i.e., of finding a nucleotide that is more conserved than the particular cutoff value. Now, the probability that there are N conserved third-codon positions together, assuming the accuracy of the null hypothesis (these are independent events), is PcN. After correcting for multiple hypotheses (Bonferroni correction, with a P of <0.01), we found the lengths, N, of nearby third-codon positions. In particular, for 98% conservation, we hold that a string of highly conserved third-codon positions with a length of >6 cannot be explained by the stratification of the data. Table 1 summarizes the conserved sequences found within the polymerase vRNAs. The results are very robust with any change in the conservation cutoff from 95 to 99%.
|
View this table: [in a new window] |
TABLE 1. Conserved sequences identified within the three polymerase subunit vRNAs for avian and human influenza A virusesa
|
Mutational analysis of the conserved residues of the PB2 gene and their role in genome packaging. The highly conserved region of the PB2 gene, nucleotides 2209 to 2304 (numbering corresponds to that of the positive-sense strand), is contained within the 5' packaging region of the PB2 gene identified by Muramoto et al. (15), which includes 120 nucleotides at the 5' end of the PB2 gene and is slightly longer than the 80 nucleotides that Liang et al. (12) reported to be essential for packaging. In order to determine the role of this highly conserved region in packaging, sequential regions of 15 to 20 nucleotides were mutated (Fig. 1).
![]() View larger version (14K): [in a new window] |
FIG. 1. Mutational analysis of the highly conserved region of the WSN PB2 vRNA. (Left) Schematic representations of the regions of synonymous mutations (white boxes) introduced into the WSN PB2 vRNA. nt, nucleotide; ORF, open reading frame. (Right) Synonymous nucleotide changes introduced for each construct. The upper line is the parental WSN virus PB2 sequence, and nucleotide changes are presented in bold on the lower line. The numbering of nucleotides is based on that of the positive-sense RNA sequence.
|
Recombinant viruses were then rescued using each of the mutated PB2 constructs. The replication of each of these recombinant viruses was assessed by infecting MDCK cells at an MOI of 0.001 and incubating them at 37°C for 48 h. Virus titers from the supernatants at 48 h are shown in Table 2, with PB2-51 growing to a level similar to that of wt WSN. Mutant viruses PB2-52, PB2-53, PB2-54, and PB2-55 all showed an approximately 2-log reduction in titer compared with wt WSN, suggesting a role for this highly conserved region in the replication cycle.
|
View this table: [in a new window] |
TABLE 2. Effect on replication in MDCK cells of recombinant WSN viruses with the introduction of synonymous changes within the highly conserved region of the PB2 vRNA
|
|
View this table: [in a new window] |
TABLE 3. Effect of synonymous changes within the highly conserved region of PB2 vRNA on the packaging of individual vRNAs into progeny WSN virions
|
In an attempt to select for mutations restoring the packaging of the PB2 vRNA, PB2-52, PB2-53, PB2-54, and PB2-55 viruses were sequentially passaged at a low MOI 8 to 10 times in MDCK cells. Even after the final passage in MDCK cells, no restoration of virus replication to wt levels was achieved. Additionally, no changes in the reduced packaging of the PB2 vRNA were observed. This suggests that more than one nucleotide is important for packaging and that multiple changes are necessary in order to restore packaging to the correct balance.
Mutational analysis of the single conserved residues of the PB2 gene and their role in vRNA packaging. While this work was being completed, Gog et al. (11) reported a similar study identifying and examining the same highly conserved sequence of the PB2 vRNA. In that study, Gog et al. created a packaging construct with a GFP gene flanked by 159 coding nucleotides of the 3' end and 166 coding nucleotides of the 5' end of the PB2 coding nucleotides, along with the untranslated regions. When this construct was transfected into cells, followed by infection at a high MOI with influenza A virus WSN, they were able to passage the GFP construct onto fresh MDCK cells. By making changes to single conserved amino acid codons, they were able to reduce the transduction of this GFP construct from 20- to 50-fold. To examine the effects of these nucleotide changes in the context of a complete PB2 vRNA, we engineered viruses with the same changes (Fig. 2).
![]() View larger version (22K): [in a new window] |
FIG. 2. Mutation of single highly conserved codons of the WSN PB2 vRNA. Synonymous nucleotide changes were introduced for each mutated codon. The upper line is the parental WSN virus PB2 sequence, and nucleotide changes are presented in bold on the lower line. Numbering of residues and changes introduced are based on the report by Gog et al. (11). mut, mutation.
|
|
View this table: [in a new window] |
TABLE 4. Effect of single residue synonymous changes in WSN PB2 mutants on the packaging of individual vRNAs into progeny virions
|
![]() View larger version (13K): [in a new window] |
FIG. 3. Mutational analysis of the short highly conserved region of the WSN PB1 vRNA. (Top) Schematic representations of the region of synonymous mutations (white boxes) introduced into the WSN PB1 vRNA. nt, nucleotide; ORF, open reading frame. (Bottom) Synonymous nucleotide changes introduced. The upper line is the parental WSN virus PB1 sequence, and nucleotide changes are presented in bold on the lower line. The numbering of nucleotides is based on that of the positive-sense RNA sequence.
|
|
View this table: [in a new window] |
TABLE 5. Effect of synonymous changes within the highly conserved region of WSN PB1 on the packaging of individual vRNAs into progeny virions
|
![]() View larger version (13K): [in a new window] |
FIG. 4. Mutational analysis of the highly conserved region of the WSN PA vRNA. (Left) Schematic representation of the regions of synonymous mutations (white boxes) introduced into the WSN PA vRNA. nt, nucleotide; ORF, open reading frame. (Right) Synonymous nucleotide changes introduced for each construct. The upper line is the parental WSN virus PA sequence, and nucleotide changes are presented in bold on the lower line. The numbering of nucleotides is based on that of the positive-sense RNA sequence.
|
|
View this table: [in a new window] |
TABLE 6. Effect on replication in MDCK cells of recombinant WSN viruses with the introduction of synonymous changes within the highly conserved region of the PA vRNA
|
|
View this table: [in a new window] |
TABLE 7. Effect of synonymous changes within the highly conserved region of PA on the packaging of individual vRNAs into progeny WSN virions
|
![]() View larger version (16K): [in a new window] |
FIG. 5. Mutational analysis of the conserved region of the PR/8 PA vRNA. (Left) Schematic representations of the regions of synonymous mutations (white boxes) introduced into the PR/8 PA vRNA. (Right) Synonymous nucleotide changes introduced at the 3' and 5' ends of the vRNA for each construct. The upper line is the parental PR/8 virus PA sequence, and nucleotide changes are presented in bold on the lower line. The numbering of nucleotides is based on that of the positive-sense RNA sequence.
|
|
View this table: [in a new window] |
TABLE 8. Effect on replication in embryonated eggs of recombinant PR/8 viruses with the introduction of synonymous changes within the highly conserved region of the PA vRNA
|
|
View this table: [in a new window] |
TABLE 9. Effect of synonymous changes within the highly conserved region of PA on the packaging of individual vRNAs into progeny PR/8 virions
|
|
|
|---|
To date there have been two different approaches to identify and map the cis-acting RNA packaging sequences in the eight segments of the influenza virus: (i) deletion of nonessential RNA sequences from a gene segment and replacement with a gene (GFP) that permits the detection of this segment when the gene is packaged into a virus particle and delivered to a new host cell (9, 10, 12, 15, 17, 25) and (ii) mutation of the essential RNA sequences for packaging the RNA followed by a quantitative measurement of the packaging of each of the eight RNA species in a virus preparation (11, 15). In the present study, we first introduced mutations (that did not change an amino acid) into the third position of a codon located in the conserved avian and human influenza virus RNA sequences of the PB1, PB2, and PA genes of the WSN and PR/8 viruses (Tables 2 to 9). These mutant viruses were replicated in cells and the virus progeny purified. Employing one RNA segment to normalize the RNA levels in the virions (the NA segment), we quantitated the levels of the other seven RNA segments in this virus preparation. Several interesting conclusions and patterns emerged from these studies. (i) Selected mutations in PB2 reduced the packaging of PB2 by up to 93% of the normalized level of NA. (ii) Mutations in the PB1 segment reduced the level of packaging of PB1 by up to 90% compared with that of a normalized NA control. (iii) Mutations in the PA segment reduced the packaging of the PA segment by up to 81% compared with that of the NA segment. Clearly, these RNA sequences play a role in segment packaging into influenza virus. However, it appears that not all of the conserved bases are involved in packaging, and the reason for this conservation is still unknown.
Perhaps more interesting were the observations that demonstrate that mutations in one segment of RNA also can have a dramatic impact upon the packaging of another segment into a virion, and this pattern is clearly nonrandom; i.e., changes in one segment affect the other segments in a nonuniform fashion. For example, selected mutations in the conserved sequence of PB2 (PB2-55) (Table 3) reduced the packaging of PB1 (86%), PA (89%), NP (88%), M (91%), and NS (92%) vRNAs. These same mutations failed to have any impact upon the levels of either HA or NA RNA segments in virions. In contrast, mutations in the conserved region of PB1 (PB2-51) that drastically reduced PB1 vRNA segment packaging (90%) had only a modest or no effect upon the reduction of packaging of PB2 (41%), PA (2%), NP (33%), or M (49%) vRNA and no impact upon HA or NA vRNA. Mutations in PA (PB2-53) reduced PA vRNA packaging by 81%, PB2 vRNA packaging by 64%, and NP vRNA packaging by 85% but had no impact upon the packaging of PB1, HA, NA, or M vRNA. These results both confirm and extend the work of Gog et al. (11), who focused upon PB2 and employed an approach different from that used in the studies reported here. The fact that different segments of the influenza virus RNA impact the packaging of other RNA segments suggests an ordered, nonrandom process of packaging where RNA-RNA and possibly RNA-protein interactions guide into the virion the eight different segments from a pool of mixed (able to undergo reassortment) RNA molecules. It is of some interest that none of the mutations in PA, PB1, or PB2 had any impact upon the packaging of the HA and NA segments and that no conserved RNA sequences that span both avian and human NA and HA sequences have been identified. This suggests that the packaging of these two segments (and therefore the reassortment of these two segments) may be accomplished by a mechanism different than the ones described here. It should be noted, however, that a virus lacking the HA RNA segment and grown in a complementing cell line expressing an HA protein also showed differential packaging of the other vRNAs (11). Thus, the HA RNA is also packaged by a mechanism which is dependent on or affects the incorporation of the other segments.
Just as the packaging of the eight RNA segments of the influenza virus appears to be nonrandom and has certain asymmetries, we have previously shown that the reassortment of the eight influenza virus gene segments does not fit a binomial distribution, which would be expected if all of the eight segments randomly reassorted into progeny viruses (13, 22). Lubeck et al. (13) first showed that a mixed infection of cells in culture by two different viruses gave rise to progeny viruses that did not fit a random segregation of all eight segments (a binomial distribution) and that would give rise to 256 different viruses. This experiment was carried out in cell culture so that no immunoselective forces acted upon the progeny obtained. The nonrandom nature of the progeny gave rise to patterns. (i) The most common reassorted virus had one new segment and seven segments from the other virus. Most often the new segment was the HA or NA segment. (ii) When the PB2 segment reassorted from one virus, it commonly brought along several other segments (PB1, PA, NP, M, and NS). These patterns are similar to those identified as important in the packaging process, and these data suggest that the packaging process guides the reassortment of segments into progeny virus in a nonrandom fashion (not a binomial distribution). A recent analysis (22) of two different human epidemics in New York and New Zealand has demonstrated a nonrandom pattern of reassortment of different vRNA segments in wild populations of viruses. In spite of the fact that both reassortment and immunoselection act upon this process, the results of that analysis in vivo (22) were similar to the results of the in vitro study of Lubeck et al. (13). Most commonly, when a single segment reassorts, it is usually the HA or NA segment. When several segments reassort, there is a complex pattern.
In conclusion, we have rescued and passaged recombinant influenza viruses with mutations within highly conserved regions located within previously identified packaging signals of the polymerase genes. These mutations resulted in a
10-fold reduction in the packaging level of that particular vRNA and also reductions in the levels of packaging of other vRNAs. This work complements other studies using GFP reporter packaging constructs to identify these regions as important for packaging, demonstrating in the context of a fully infectious eight-segmented virus the importance of these sequences. Further studies of the packaging regions in other RNA segments of a variety of subtype viruses will be necessary in order to better understand the precise molecular mechanism by which vRNAs are incorporated into influenza virions.
Partial support of this work was provided by NIH grants PO1AI058113, Northeast Biodefense Center grant U54 AI05718, and UO1AI070469 and Center for Research on Influenza Pathogenesis grant HHSN2662000700010C. The work of A.J.L. and R.R. was supported by the Simons Foundation, the Ambrose Monell Foundation, and the Leon Levy Foundation.
Published ahead of print on 19 December 2007. ![]()
These two authors contributed equally to this work. ![]()
# Present address: CSIRO Livestock Industries, CSIRO Australian Animal Health Laboratory, P.O. Bag 24, Geelong, Victoria 3220, Australia. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»