Previous Article | Next Article ![]()
Journal of Virology, May 2002, p. 4750-4763, Vol. 76, No. 10
0022-538X/02/$04.00+0 DOI: 10.1128/JVI.76.10.4750-4763.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Department of Veterinary Diagnostic and Production Animal Medicine,1 Department of Statistics, Iowa State University, Ames, Iowa 50010,2 Department of Veterinary PathoBiology, University of Minnesota, St. Paul, Minnesota 551083
Received 6 November 2001/ Accepted 12 February 2002
|
|
|---|
|
|
|---|
PRRS virus (PRRSV) is a member of the family Arteriviridae, along with equine arteritis virus, lactate dehydrogenase-elevating virus of mice, and simian hemorrhagic fever virus (9, 48). PRRSV is enveloped and has a polyadenylated, single-stranded, nonsegmented, positive-sense RNA genome 15 kb in size (7, 14, 34, 35). The genome consists of eight open reading frames (ORFs) that are expressed through the production of a nested set of subgenomic 3'-coterminal mRNAs (14, 34, 35). ORF 1, which comprises two-thirds of the genome, encodes the viral RNA-dependent RNA polymerase (14, 37). ORFs 2 to 7 are postulated to encode structural proteins, but only three proteins have been consistently identified in virions and/or lysates of virus-infected cells. These three are the 15-kDa nucleocapsid (N), 19-kDa matrix (M), and 25-kDa envelope (E) glycoproteins that are encoded by ORFs 7, 6, and 5, respectively (41, 42). Proteins encoded by ORFs 2 to 4 are designated GP2, GP3, and GP4, where GP indicates glycoprotein and the number designates the ORF from which it is derived (36, 56).
PRRSV infection is prevalent in swine-producing regions throughout the world. Infection with PRRSV produces viremia with subsequent dissemination to and viral replication in multiple organs (50, 52). The virus replicates predominantly in cells of the monocyte/macrophage lineage (47, 50, 52, 57). Clinically, infection may range from totally inapparent to causing severe disease. Generally, clinical signs involve reproductive disorders in pregnant animals and/or respiratory disease in pigs of all ages (11, 20, 21, 23, 47, 51, 63). Reproductive disease in pregnant animals is manifested as late-term abortions or premature farrowings. Affected litters have higher proportions of stillbirths and piglets born weak and increased preweaning mortality (11, 23, 63).
PRRSV produces a persistent infection despite an active immune response (1, 6, 20, 27, 60). The PRRSV carrier state was first recognized in an experiment that documented the transmission of virus from animals infected 99 days earlier to commingled sentinel pigs (66). Subsequently, Wills et al. (60) reported the isolation of virus at up to 157 days postinoculation (p.i.). Persistence is an important epidemiological feature because it provides a ready means for PRRSV to perpetuate itself through a cycle of transmission from carrier to susceptible animals. As a consequence, the elimination of PRRSV from herds is difficult, and cyclic bouts of PRRSV-associated health problems are commonplace.
The mechanism(s) by which the virus persists in the host is not known, but PRRSV isolates are characterized by a high degree of genetic and antigenic variability. Persistence in the host and viral diversity could be two manifestations of the same function. Genetic analyses have shown the existence of at least two major virus genotypes, the European and the North American, with extensive genetic variability both within and between these genotypes. Nelsen et al. (40) found important differences between the prototypic North American (VR-2332) and European (Lelystad) viruses in the 5' leader sequence and parts of ORF 1a. Marked differences were also found between European and North American isolates in some structural genes (24, 39). ORF 7 is highly conserved among North American isolates, with 95 to 100% amino acid homology, but a comparison of North American viruses and the Lelystad virus revealed only 57 to 59% amino acid homology (32, 39). ORF 6 is the most conserved gene among North American isolates, with up to 100% amino acid identity, and is the most conserved gene between North American and European isolates, with 70 to 81% identity (24, 33, 39). The amino acid sequence homology for ORF 5 varies from 88 to 97% among North American isolates and from 51 to 59% when North American viruses are compared to the Lelystad virus (4, 24, 31, 39). A comparison of the Lelystad virus with isolate ATCC VR-2332, the North American prototype virus, revealed amino acid identities of 63, 58, and 68% for ORFs 2, 3, and 4, respectively (39). Similar results were reported when the Lelystad virus was compared with U.S. isolate VR-2385 (38).
Genetic diversity is mirrored in antigenic diversity among PRRSV isolates. Antigenic variation was initially demonstrated in a comparison of European and North American isolates. Using an immunoperoxidase monolayer assay, Wensvoort et al. (58) evaluated the reactivity of porcine polyclonal antibodies raised against either the Lelystad virus or North American isolate ATCC VR-2332 with PRRSV isolates from around the world. The investigators were able to differentiate European from North American isolates on the basis of differences in immunoperoxidase monolayer assay antibody titers. That is, significantly higher antibody titers were obtained in the homologous assay system. Later studies found even greater antigenic diversity among PRRSV isolates than initially suspected (5, 15, 17, 28, 64). Yoon et al. (64) examined 22 PRRSV isolates recovered from samples collected between 1989 and 1993 in eight different U.S. states. With a panel of five monoclonal antibodies (MAbs) specific for the N protein, the 22 virus isolates fell into one of three groups based on their reactivity patterns. Yang et al. (61) expanded this study to 70 North American isolates from samples collected between 1989 and 1995 and a panel of 23 MAbs against the N protein. Their results showed five antigenic groups, with the European Lelystad virus representing an antigenic group distinct from any of the North American groups identified. Furthermore, with antibodies against discontinuous epitopes of the N and M proteins and continuous epitopes of the E protein and GP3, the 65 North American isolates in the first and second antigenic groups (I15 and II15) were further subdivided into nine and four antigenic subgroups, respectively (62).
That field isolates of PRRSV show a remarkable degree of genetic and antigenic variability has become abundantly evident. Understanding the source and degree of virus diversity is important because of its potential impact on viral virulence, persistence in the host, protective immunity, and even on the ability to diagnose the infection. The assumption is that these changes occur during the course of in vivo replication in swine and arise, in large part, because of errors that occur during RNA replication. However, the degree and rate of mutation of PRRSV in infected pigs over time are not known. To address this question, we examined the genetic and antigenic changes that occurred in a well-characterized PRRSV isolate during the course of multiple pig passages over a period of 367 days.
|
|
|---|
![]() View larger version (22K): [in a new window] |
FIG. 1. Experimental design. The diagram illustrates three independent lines (A, B, and C) and seven pig passages (P1 to P7) of PRRSV. Pigs 1A, 1B, and 1C were inoculated with plaque-cloned PRRSV. Pigs in P2 to P7 were inoculated with tissue filtrates from the corresponding pig in the previous passage.
|
Viruses, cells, and media. The North American prototype PRRSV ATCC VR-2332 (American Type Culture Collection, Manassas, Va.) was used in the study. VR-2332 was selected because the complete genomic sequence has been published (GenBank accession no. PRU87392) (40) and the epitopic profile has been described (61, 62). The virus was propagated on MARC-145 cells, a clone of the African monkey kidney cell line MA-104 that is considered highly permissive to PRRSV (25). The cells were cultured and maintained in culture medium composed of Dulbecco's modified Eagle's medium (DMEM; Sigma Chemical Co., St Louis, Mo.) supplemented with 10% fetal bovine serum (HyClone Laboratories, Inc., Logan, Utah), 50 µg of gentamicin (Sigma)/ml, and 0.25 µg of amphotericin B (Fungizone; Sigma)/ml.
To produce a highly homologous challenge virus, the stock virus was subjected to three rounds of plaquing on MARC-145 cells. Briefly, the stock virus was 10-fold serially diluted (i.e., 10-1 to 10-6) with culture medium. Two milliliters of each dilution was inoculated into each well of a six-well plate (Corning Inc., Corning, N.Y.) containing 24-h-old confluent MARC-145 cells. Inoculated cells were incubated for 2 h at 37°C in a humidified 5% CO2 incubator. At the end of the incubation period, the inocula were discarded and the cells were rinsed twice with culture medium and then covered with 4 ml of overlay medium composed of 0.2% agarose (FMC Bioproducts, Rockland, Maine), 10% fetal bovine serum, 50 µg of gentamicin/ml, and 0.25 µg of amphotericin B/ml in DMEM. The cells were incubated for 3 to 7 days at 37°C in a humidified 5% CO2 atmosphere. During the incubation period, well-demarcated plaques were selected, suspended in 2 ml of culture medium, and then propagated once on MARC-145 cells. This process was repeated three times. Following the third round of plaquing, one virus clone, designated VR-2332-CC-01 (hereafter referred to as CC-01), was selected and propagated once on MARC-145 cells to obtain an amount sufficient for inoculating the pigs in P1. Surplus CC-01 was divided into aliquots and stored at -80°C for future use.
The same procedure was used to recover plaque-cloned virus isolates from pigs. That is, 15 plaque-cloned virus isolates were recovered from each of the three pigs in each of the seven passages by directly diluting the serum samples collected on day 7 p.i. with the procedure described above. The use of viremic serum samples made it possible to recover viruses directly from pigs and eliminated the requirement for virus isolation in cell cultures prior to plaque cloning. This process also served to reduce the selection effect of cell culturing on the recovery of plaque-cloned viruses. It was not expected that sequencing 15 clones from each pig and passage would identify all the unique nucleotide and amino acid sequences present in each pig. A rare sequence was unlikely to be included in the sample. However, if the true prevalence of a sequence exceeded 18%, then it was presumed with 95% or higher probability that we would be able to observe at least one clone of that sequence in our sample. Selected plaque-cloned viruses were propagated once on MARC-145 cells, divided into aliquots, and stored at -80°C until used. Fifteen plaque-cloned viruses were also obtained from the original virus inoculum to document the genetic and antigenic characteristics of CC-01, assess the variability of the virus population within the inoculum, and provide a baseline for comparisons.
Animals and animal care. Throughout the study, animals were housed and cared for in compliance with the requirements given in Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching (17a).
For the study, crossbred pigs were obtained from a closed specific-pathogen-free herd known to be free of PRRSV. Pigs were weaned at 10 to 14 days of age, ear tagged, and randomly assigned to a treatment. Weaned pigs were group housed in HEPA-filtered isolation units for several days prior to exposure to virus and then individually housed in the isolation units from the time of exposure until the end of the observation period. The isolation units were equipped with internal flush mechanisms that allowed for disposal of waste products while maintaining a biosecure environment. Likewise, a sealed system permitted feeding without jeopardizing biosecurity. Contact with the external environment and other animals was eliminated to the greatest degree possible. Inoculations and sample collections were scheduled so that only one pig was removed from its isolation unit in any 24-h period. Immediately following inoculation or sample collection, the environment was cleaned and then disinfected with chlorhexidine diacetate disinfectant (Nolvasan Solution; Fort Dodge Laboratories, Fort Dodge, Iowa).
Biological samples. Blood samples were collected from all pigs on days 0, 7, 14, 21, 28, 35, and 60 p.i. with a single-use blood collection system (Vacutainer; Becton Dickinson, Franklin Lakes, N.J.). After 30 min at room temperature, the serum was harvested by centrifugation of the sample at 1,000 x g for 10 min. At day 60 p.i., pigs were euthanatized and necropsied. Whole blood, tonsils, tracheobronchial and medial iliac lymph nodes, spleen, lungs, and bronchoalveolar lavage fluid were collected from each animal. Whole blood samples were collected with a single-use blood collection system containing EDTA (Vacutainer) and processed immediately after sampling to obtain peripheral blood leukocytes (PBLs) as previously described (43). Briefly, whole blood in EDTA was centrifuged at 800 x g and 4°C for 30 min. The PBLs at the interface of the plasma and red blood cells were collected and resuspended in 3 ml of Hanks' balanced salt solution (HBSS; Sigma). The cell suspension was then layered over an equal volume of Histopaque-1077 (Sigma) and centrifuged at 1,200 x g and 4°C for 20 min. The PBLs at the interface of the upper medium layer and the lower erythrocyte layer were collected and then washed twice by resuspension in 10 ml of HBSS and centrifugation at 200 x g for 10 min. All samples were stored at -80°C.
To prepare inocula for P2 through P7, biological samples (tonsils, lymph nodes, spleen, lungs, PBLs, and lung lavage fluid) were suspended in cold HBSS at approximately 20% (wt/vol). The pool was homogenized in a Stomacher 80 (Seward, London, United Kingdom) for 2 min and then centrifuged at 4,000 x g for 30 min at 4°C. The supernatant was passed through a 0.2-µm-pore-size nonpyrogenic nitrocellulose membrane filter (Corning) and used to inoculate the next pig.
ELISA.
A commercial enzyme-linked immunosorbent assay (ELISA) kit (IDEXX Laboratories, Inc., Westbrook, Mass.) for the detection of antibody specific for PRRSV was used by following the directions supplied by the manufacturer. According to the manufacturer, samples for which the ratio of net optical density of test sample to net optical density of positive control (S/P) is
0.4 are positive for antibody against PRRSV.
Virus titration. A microtitration infectivity assay was performed to assess the virus titer of the original inoculum, as well as the levels of PRRSV in serum samples collected over time from pigs in P1 through P7. Briefly, samples were 10-fold serially diluted (i.e., 100 to 10-6) in culture medium. One hundred microliters of each dilution was added to each of 3 wells of a 96-well microtitration plate (Corning) containing 24-h-old confluent MARC-145 cell monolayers. Inoculated cells were incubated at 37°C in a humidified 5% CO2 incubator. Each sample was run in duplicate. The cells were monitored daily for cytopathic effects (CPE) for up to 7 days. If CPE were not evident, the cells were fixed with aqueous 80% acetone solution, dried, stained with fluorescein isothiocyanate-conjugated MAb SDOW17 (Rural Technologies, Brookings, S.Dak.), which is specific for the N protein of PRRSV, and visualized with fluorescence microscopy. The presence of PRRSV was determined based on the observation of virus-specific CPE and/or fluorescence reaction. Virus titers were determined by using the method described by Reed and Muench (49) and were expressed as TCID50 per milliliter.
VN assay. A one-way virus neutralization (VN) assay with constant antibody and various virus concentrations was performed as described previously (62) to compare the relative susceptibilities of pig-derived plaque-cloned virus isolates to the neutralizing activity of serum and to screen for potential escape mutant viruses. A total of 105 plaque-cloned viruses were examined in the one-way VN assay, i.e., 15 clones from the CC-01 inoculum, 45 clones from P1 (15 clones from each of three pigs), and 45 clones from P7 (15 clones from each of three pigs). The P1 clones represented 7 days and the P7 clones represented 367 days of in vivo replication.
To perform the assay, virus clones were 10-fold serially diluted in DMEM supplemented with 50 µg of gentamicin/ml and 0.25 µg of amphotericin B/ml. Serum collected from each of three pigs (1A, 1B, and 1C) at the termination of P1 (day 60 p.i.) served as the source of antibody. The assay was performed in duplicate for each combination of antibody source and virus clone. Serum samples were heat inactivated at 56°C for 45 min and then diluted 1:2. Fifty microliters of each diluted virus clone was mixed with an equal volume of diluted serum. An additional 50 µl of each diluted virus clone mixed with an equal volume of cell culture medium instead of antiserum served as an untreated control. All mixtures were incubated for 1 h at 37°C, and then 100 µl of each mixture was added to 96-well plates containing 24-h-old confluent MARC-145 cell monolayers in the same manner as that described above for virus titration. Inoculated cells were incubated for up to 7 days. At the end of the 7-day incubation period, the virus titer of each clone was calculated for each plaque clone in the presence and absence of antiserum and was expressed as TCID50 in log10 units. The relative susceptibility of plaque-cloned viruses was calculated with the following equation and was expressed as the VN index: [(TCID50 without antiserum - TCID50 with antiserum)/TCID50 without antiserum] x 100.
MAb analysis. The epitopic profile of each of 150 plaque-cloned viruses, i.e., 15 clones from the CC-01 inoculum, 45 clones from P1 (day 7 of in vivo replication), 45 clones from P2 (day 67 of in vivo replication), and 45 clones from P7 (day 367 of in vivo replication) was determined with a panel of MAbs in an indirect immunofluorescence assay (61). The panel consisted of five MAbs against the N protein and two against the E protein. The production and characterization of the MAbs have been described elsewhere (61, 62).
To prepare PRRSV-infected cell monolayers, 24-h-old confluent MARC-145 cells in 96-well plates were inoculated with optimally diluted virus. Procedures for control wells were identical except that the cells were not inoculated with virus. At 24, 48, or 72 h, plates were fixed by immersion in aqueous 80% acetone solution for 10 min, air dried, and stored at -20°C until used. MAbs used in the assay were optimally diluted with 0.01 M phosphate-buffered saline (PBS, pH 7.2). Individual wells received 40 µl of each MAb dilution, after which plates were incubated at 37°C for 45 min in a humid environment. After the MAb solution was discarded, plates were washed three times with PBS. To visualize the presence of antigen-antibody complexes, each well received 40 µl of 1:100-diluted goat anti-mouse immunoglobulin G conjugated to fluorescein isothiocyanate (Kirkegaard & Perry Laboratories, Gaithersburg, Md.), after which plates were incubated for 45 min in a humid environment at 37°C. Plates were then washed three times with PBS and air dried. The reactivity of each MAb with each virus clone was evaluated by using fluorescence microscopy.
RT-PCR and sequencing. A total of 150 plaque-cloned viruses were sequenced for ORFs 1b and 7, i.e., 15 clones from the CC-01 inoculum, 45 clones from P1 (day 7 of in vivo replication), 45 clones from P2 (day 67 of in vivo replication), and 45 clones from P7 (day 367 of in vivo replication). A total of 330 virus clones were sequenced for ORF 5, including 15 clones from the CC-01 inoculum and 45 plaque-cloned viruses collected from each of the seven passages (n = 315). Fifty plaque-cloned viruses from the CC-01 inoculum were included in each sequencing reaction to assess the homogeneity of the inoculum and any random mutation(s) which could be introduced during reverse transcription (RT)-PCRs or sequencing reactions.
Viral RNA for RT-PCR amplification and sequencing was extracted from each plaque-cloned virus with a QIAamp viral RNA minikit (Qiagen Inc., Valencia, Calif.) by following the protocols recommended by the manufacturer. Viral RNA was collected and stored at -80°C until used.
The RT-PCR of ORF 1b was performed with 5 µl of extracted RNA by using a Qiagen OneStep RT-PCR kit according to the manufacturer's protocol and primers 8100 and VRBP2 (Table 1). RT was performed for 30 min at 50°C. Reverse transcriptase was inactivated and Taq polymerase was activated by raising the temperature to 95°C for 15 min. PCR amplification was achieved with 35 cycles of annealing at 57°C for 30 s, extension at 72°C for 45 s, and denaturation at 94°C for 45 s. Products were polished by incubation at 72°C for 7 min.
|
View this table: [in a new window] |
TABLE 1. Primers used for PCR and sequencing
|
The ORF 7 RT-PCR procedure was identical to that described above for ORF 5 except for the use of a pair of sense (P7F) and antisense (P7R) primers (Table 1), designed as previously reported (12), and 58°C instead of 50°C for the annealing process.
In every round of RT-PCR, RNA extracted from reference PRRSV was used as a positive control and mock-infected cell culture fluid was used as a negative control. The final RT-PCR products were purified with either QIAquick or 96-well PCR purification kits (Qiagen) by following the procedures recommended by the manufacturer. The final purified product was collected and stored at 4°C until sequenced.
Sequencing reactions for ORF 1b were carried out at the University of Minnesota Advanced Genetic Analysis Center (St. Paul) by using 4 µl of purified PCR product and 3.2 pmol each of primers RESP1BP.1 and VRBP2 (Table 1). Sequencing of ORF 5 and ORF 7 was performed at the Iowa State University Nucleic Acid Facility (Ames) by using the appropriate amount of purified PCR product (30 µg/µl) and 5 pmol of each primer. Primers P5F, P5R, PS1, and PS1R were used for ORF 5 and primers P7F and P7R were used for ORF 7 (Table 1). For all three ORFs, pairs of primers were primed from the positions before and after the target sequences. However, for ORF 5, PS1 and PS1R were also used to increase sequence accuracy and were primed from the middle of ORF 5 in both directions.
Analysis of data.
ELISA S/P values were analyzed by fitting a regression model to longitudinal data (16). An appropriate model for the time course in each pig was chosen by using plots of the results and observations of antibody responses. A simple model that fit the data from day 7 through day 60 was ELISA S/P =
+ ß log (day) +
. Differences among lines and generations were tested separately for each parameter of the model.
Viremia data from days 7 to 35 were analyzed by analysis of variance and by longitudinal regression methods. Data for days 0 and 60 were excluded because no virus was detected at those times. All analyses were done with the average of two replicate determinations of log10 TCID50 for each pig and day. The analysis of variance was done with a split-plot analysis to account for the repeated measures across days. However, because there was a significant three-way interaction among passage, lines (A, B, and C), and day p.i., longitudinal data analysis (16) was used to find a simpler model for the time course of viremia in each pig. A quadratic regression equation, log10 TCID50 =
+ ß day +
day2 +
, was fit to data from each pig. In the model,
was an intercept that was the same for all groups, ß was the slope of the line,
was the error associated with each observation, and
was the quadratic coefficient. The day of peak viremia (Dpeak) was estimated for each pig from the quadratic regression coefficients as follows: Dpeak = -ß/2
. Differences among lines and generations were tested separately for each parameter.
Data from the VN assay were summarized by using the mean VN index for each combination of isolate and P1 antiserum. An escape mutant was considered to be a virus clone with a significantly smaller mean VN index.
Sequence alignment, comparisons, and phylogenetic analyses were carried out by using computer software from DNASTAR Inc., Madison, Wis. Nucleotide sequences of VR-2332 (40) and the RespPRRS vaccine strain (GenBank accession no. AF159149) (2) were obtained and included in the sequence analyses for comparative purposes. Like the RespPRRS vaccine strain, CC-01 is derived from VR-2332. Nucleotide mutation rates and amino acid sequence changes over the course of the experiment were compared to those of CC-01. The mutation rates were calculated as the proportion of substitutions and were expressed as mean percent change separately for ORFs 1b, 5, and 7 per plaque-cloned virus. Phylogenetic trees representing amino acid variations in the E protein (ORF 5) and the polymerase (ORF 1b), GP5, and the N protein combined were computed by using the Hein method with the PAM250 weight table (22).
|
|
|---|
ELISA-detectable antibody responses and PRRS viremia in successive passages (P2 to P7) showed that inoculation of pigs with tissue homogenate filtrates resulted in transmission of the infection. All pigs in lines A, B, and C showed similar humoral immune responses, as measured by the ELISA (Fig. 2). No statistically significant differences were found in the ELISA responses among lines A, B, and C for any sampling time.
![]() View larger version (15K): [in a new window] |
FIG. 2. Antibody response to PRRSV infection. Values represent means and standard deviations of ELISA S/P ratios in line A (), line B ( ), and line C ( ) over all passages (P1 to P7).
|
![]() View larger version (18K): [in a new window] |
FIG. 3. Level and duration of viremia. Values represent means and standard deviations in pigs (P1) inoculated with cell culture-derived PRRSV ( ) and pigs (P2 to P7) inoculated with pig-passaged virus in line A (), line B ( ), and line C ( ).
|
|
View this table: [in a new window] |
TABLE 2. Estimated day of peak viremia by line and passage
|
Genetic characterization of the CC-01 inoculum. Partial ORF 1b sequences (435 bases), complete ORF 5 nucleotide sequences (603 bases), and complete ORF 7 nucleotide sequences (372 bases) were obtained for the 15 plaque-cloned viruses derived from the CC-01 inoculum. The ORF 1b, 5, and 7 sequences of all 15 clones were identical during repeated sequencing attempts. This result demonstrated the homogeneity of the original virus population, proved that no random mutation was introduced during PCR and sequencing, and provided a baseline for comparing and contrasting virus clones recovered in later passages.
Relative to VR-2332, the ORF 5 nucleotide sequence of CC-01 differed at three positions: 38 (G to A), 252 (C to T), and 451 (A to G). These changes resulted in amino acid substitutions at two residues: 13 (arginine to glutamine) and 151 (arginine to glycine). Compared to the RespPRRS vaccine virus, the ORF 5 sequence of CC-01 had one synonymous nucleotide change at position 252 (C to T). In contrast, the ORF 1b nucleotide sequence of CC-01 was identical to those of VR-2332 and the RespPRRS vaccine virus. The ORF 7 nucleotide sequences of VR-2332 and the RespPRRS vaccine virus were identical, but CC-01 differed in two synonymous mutations at positions 30 (G to T) and 345 (C to T).
Assessment of genetic changes. (i) Sequence analysis of PRRSV ORF 5. Complete ORF 5 nucleotide sequences were obtained for 315 swine-derived virus clones (15 clones per pig, three pigs per passage, and seven passages). These clones represented 7, 67, 127, 187, 247, 307, and 367 days of in vivo replication in P1 through P7, respectively. As shown in Table 3, three or more ORF 5 nucleotide variants (5NVs) were detected in every pig beginning as early as day 7 p.i. Six 5NVs were recovered from one pig. Among the 315 virus clones, 291 differed from CC-01 in the ORF 5 nucleotide sequence. In fact, no clone with 100% homology to the CC-01 ORF 5 sequence was detected after P1. In total, 48 5NVs, each containing from one to eight nucleotide substitutions, were observed. No nucleotide insertions or deletions were detected.
|
View this table: [in a new window] |
TABLE 3. Chronological appearance of 48 5NVs within three independent lines and over seven pig passages
|
|
View this table: [in a new window] |
TABLE 4. Chronological appearance of 22 5AVs within three independent lines and over seven pig passages
|
The ORF 5 sequences of CC-01, VR-2332, the RespPRRS vaccine virus, and the 22 amino acid variants were aligned and compared (Fig. 4). Among the 22 variants, no specific "hot spots" were identified, except at position 151 (glycine to arginine). Most amino acid substitutions due to nucleotide changes were observed in amino acid residues located between residues 10 and 34. Eleven of the 22 5AVs had an amino acid substitution at either residue 33 (asparagine to serine) or residue 34 (asparatic acid to asparagine or serine), which is part of the ectodomain of GP5. An amino acid substitution at residue 33 abolished a postulated variable N-linked glycosylation site. Phylogenetic analysis (Fig. 5) revealed two distinct branches of evolution among the 22 variants, with 8 of 9 variants from line B on one branch and 13 of 14 variants from lines A and C on the other branch.
![]() View larger version (27K): [in a new window] |
FIG. 4. GP5 amino acid sequence comparisons of CC-01, 22 amino acid variants, VR-2332, and the RespPRRS vaccine virus. Dots indicate bases identical to those in CC-01, and letters indicate amino acid substitutions. Positions 41 to 90 and 101 to 120 have been truncated.
|
![]() View larger version (17K): [in a new window] |
FIG. 5. Phylogenetic relationship of 22 PRRSV 5AVs detected during 367 days of in vivo replication. Letters in parentheses represent the lines (A, B, and C) of pig passages in which the variant was recovered. Lines represent the distances of sequence divergence from the ancestor (CC-01).
|
Relative to CC-01, a total of 16 ORF 1b nucleotide variants (1NVs) (Table 5) were observed, each containing from one to four nucleotide changes. In P1, only one variant was detected among 45 virus clones, but all 45 virus clones from P7 were variants. All nucleotide changes were randomly distributed, and only one ORF 1b amino acid variant (1AV) was found (Table 6). Specifically, the mutation at position 207 (T to C) in one virus clone recovered from line A in P7 resulted in an amino acid substitution (valine to alanine). The remaining nucleotide mutations were synonymous.
|
View this table: [in a new window] |
TABLE 5. Appearance of 16 1NVs within three independent lines and in three pig passages
|
|
View this table: [in a new window] |
TABLE 6. Appearance of one 1AV within three independent lines and in three pig passages
|
|
View this table: [in a new window] |
TABLE 7. Appearance of nine 7NVs within three independent lines and in three pig passages
|
|
View this table: [in a new window] |
TABLE 8. Appearance of three 7AVs within three independent lines and in three pig passages
|
|
View this table: [in a new window] |
TABLE 9. Appearance of 40 CNVs within three independent lines and in three pig passages
|
|
View this table: [in a new window] |
TABLE 10. Appearance of 19 CAVs within three independent lines and in three pig passages
|
![]() View larger version (16K): [in a new window] |
FIG. 6. Phylogenetic relationship of 19 PRRSV CAVs. Letters (A, B, and C) represent the lines of pig passages in which the variant was detected during 367 days of in vivo replication. Lines represent the distances of sequence divergence from the ancestor (CC-01).
|
![]() View larger version (18K): [in a new window] |
FIG. 7. Mutation rates for PRRSV ORFs 1b, 5, and 7 over seven animal passages (367 days of in vivo replication). The mutation rates were calculated as the proportions of substitutions in each plaque-cloned virus and are represented as the mean percent change on the nucleotide () and amino acid ( ) levels. Each symbol represents the mean of mutation rates for 45 cloned isolates, except for 44 cloned viruses in P2 for ORF 1b. Error bars are the standard errors of the mean.
|
The types of mutations occurring in ORFs 1b, 5, and 7 during animal passages are summarized in Tables 11 and 12. Overall, transitions (G 171 A or T 171 C) were more common than transversions (GA 171 TC). The numbers of transitions per plaque-cloned virus during passage were consistent with the mutation rates of the nucleotide sequences in the corresponding genes, since transversions occurred only in a few isolates. While ORFs 1b and 7 were relatively conserved during serial animal passages, ORF 5 showed stable increases in substitutions per clone from 0.58 to 3.87 for P1 through P7.
|
View this table: [in a new window] |
TABLE 11. Mean number of nucleotide substitutions per plaque-cloned virus in ORF 5 (603 bases) of PRRSV over seven animal passages
|
|
View this table: [in a new window] |
TABLE 12. Mean number of nucleotide substitution per plaque-cloned virus in ORF 1b (435 bases) and ORF 7 (372 bases) of PRRSV in three animal passages
|
|
|
|---|
Several strengths of the experimental design are worth noting. First, a plaque-cloned virus inoculum, i.e., a uniform virus population, rather than a whole virus isolate composed of multiple variants (quasispecies) was used to initiate the infection in P1. This strategy provided a definitive ancestor against which to compare and contrast descendant isolates. Second, in P2 to P7, the use of tissue homogenates to pass virus directly from pig to pig circumvented the selection pressures arising from the in vitro isolation and propagation of viruses. Third, the 60-day passage allowed the immune system to respond and exert selective pressure on the virus population. Fourth, the use of HEPA-filtered housing systems and adherence to strict biosecurity procedures made it possible to maintain virus evolution in three independent lines (replicates) of pig passages. Fifth, the recovery of 15 plaque-cloned isolates directly from serum made it possible to randomly sample the population of PRRSV variants circulating in each pig and characterize their distribution.
For the purposes of data interpretation, it should be recognized that the sample size of 15 clones did not identify all the unique nucleotide and protein sequences present. Specifically, if the prevalence of a variant in the population of viruses circulating in a pig exceeded 18%, the probability of recovering at least one clone in the sample was greater than 0.95. As a corollary, the proportion of any sequence among the 15 clones was only an estimate of the prevalence in the pig. The exact prevalence values should not be overinterpreted because there is a moderate amount of uncertainty associated with these estimates. For example, if a sequence was observed in 12 of the 15 clones (prevalence of 80%), the 95% confidence interval for that prevalence estimate would be 52 to 95%.
Successful transmission of PRRSV via tissue samples collected at 60 days p.i. provided additional evidence that PRRSV infection persists in swine. However, it was observed that cell-adapted virus behaved much differently in pigs than pig-adapted virus. That is, inoculation of pigs in P1 with a cell-adapted clone derived from ATCC VR-2332 produced a pattern of viremia strikingly different, both in virus titer and in duration of viremia, from that seen with such inoculation of pigs in P2 to P7 (Fig. 3 and Table 2). These data are important because they suggest that the results of experiments may differ significantly depending on whether cell-adapted or pig-adapted viruses are used. Although no genetic basis for this biological difference was discovered, it is recommended that future studies with direct application to the field use PRRSV with a minimum amount of cell culture adaptation.
To date, no specific genetic markers have been linked to either the attenuation or the virulence of PRRSV. However, it is interesting to note the consistent substitution of a glycine residue with arginine at position 151 of ORF 5 in all virus clones after just one passage in pigs (Fig. 4). This change represents a reversion to the amino acid residue of the parental strain, VR-2332, a field isolate. Both PRRSV CC-01 and the RespPRRS vaccine virus were derived from VR-2332 through cell culture passage and, for the ORFs sequenced, shared 100% amino acid homology. Reversion of the amino acid residue at position 151 of ORF 5 has been reported for the RespPRRS vaccine virus and VR-2332 (2, 3, 40, 54, 59). Cumulatively, these results lead us to speculate that the glycine residue at position 151 of ORF 5 may be a marker for cell culture adaptation. Whether this specific change affects the virulence or in vivo replication of PRRSV remains to be determined.
Virus evolution in the three independent lines of pig passages was independent, i.e., highly random. Although numerous genotypes were identified, with few exceptions, a single variant was dominant in each pig at each passage. Unexpectedly, the dominant variant differed from one pig to the next, even within the same line. The minor variants within a population were often different from the dominant form by more than one mutation. These observations suggested that a "purifying selection" may have occurred. This was especially true for line C, where after seven passages and a total of 367 days of in vivo replication, there was essentially only one kind of virus, based on ORF 5. No evidence was obtained for radiation, for an increase in diversity, or for the establishment of a new, dominant mutation in the population. Rather, it appeared that a single genotype gained a clear competitive advantage and became the dominant form within a specific pig or passage.
The results of this study demonstrated that PRRSV evolves continuously in infected pigs, with different genes of the viral genome undergoing various degrees of change. As anticipated for a housekeeping gene, PRRSV ORF 1b, which encodes the RNA polymerase domain, was highly conserved. ORF 7 (N protein) was also conserved during animal passages relative to ORF 5 (E protein), suggesting that the internal structural protein is genetically more stable than the external structural proteins. This scenario may be due to the immune pressures to which the external structural proteins are exposed. The relatively low rates of mutation observed in ORFs 1b and 7 are consistent with field observations (4, 18, 55).
Relative to ORFs 1b and 7, ORF 5 mutated at a much higher rate; numerous genotypes (i.e., quasispecies) were identified, with mutiple genotypes being present concurrently in all pigs. Based on an analysis of algorithms for hydrophobicity and hydrophilicity and surface probability, six potential antigenic sites were identified on the surface of GP5 of both North American and European isolates (4, 45). Sequence analyses of field isolates of PRRSV have revealed that variable regions are located among amino acid residues near either the N terminus (3 to 39) or the C terminus (164 to 200), whereas the most conserved regions (II to IV) are located between amino acid residues near 40 to 163 (4, 45). In agreement with previous reports, our study found a high degree of amino acid substitutions between residues 10 and 34, a region which contains the signal peptide (1 to 31) and one variable N-linked glycosylation site, whereas no or a few inconsistent changes were observed in amino acid residues 40 through 160. Many ORF 5 variants had amino acid substitutions at residue 33 (asparagine to arginine) or 34 (aspartic acid to asparagine or arginine). Since amino acids at positions 33 and 34 are known to be part of the ectodomain of GP5, any change in this region may affect the N-linked glycosylation or the conformation of GP5, which is important for interactions with the M protein (19, 24, 29, 30, 32, 45, 44, 46, 53). Theoretically, such a change could result in the loss of an epitope or an alteration in the immune response, but no such changes were found in this study.
The appearance of ORF 5 quasispecies during in vivo replication suggests a potential role for mutation in viral persistence. However, the overall mutation rate for ORF 5 was lower than expected on the basis of the degree of diversity seen in field isolates (24). The low rate of mutation in ORF 5 suggests that ongoing mutation is not likely to be the mechanism by which PRRSV eludes the immune system and persists in pigs.
It is believed that North American and European PRRSV isolates were derived from a common lactate dehydrogenase-elevating virus-like ancestor (40). Based on the relatively low rate of mutation observed in this study and the 25 and 36% amino acid divergences for ORFs 1b and 7 between strain VR-2332 and the Lelystad virus (39, 40), we conclude that these viruses diverged in the very distant past and experienced independent evolutionary events. Given that the earliest evidence places the introduction of PRRSV into domestic swine in 1979 (8), the slow rate of mutation documented in this study suggests that the extensive antigenic and genetic diversity observed in PRRSV field isolates cannot be fully explained simply as a function of mutations. Alternate hypotheses to explain the variability observed in field isolates would include recombination, as shown among field isolates (24) or in vitro (65), and possible ongoing evolution in alternate host species (67). At present, the factors and conditions contributing to the diversity of PRRSV remain to be established.
This study was supported in part by funding from Check-Off dollars on behalf of the National Pork Board, an Iowa Healthy Livestock Initiative research grant, and the Iowa Livestock Health Advisory Council.
|
|
|---|
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»