Previous Article | Next Article ![]()
Journal of Virology, February 2004, p. 2158-2163, Vol. 78, No. 4
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.4.2158-2163.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Ellen W. Collisson,1* Margaret Slater,2 and Guan Zhu1
Department of Veterinary Pathobiology,1 Department of Public Health and Anatomy, College of Veterinary Medicine, Texas A&M University, College Station, Texas 778432
Received 1 August 2003/ Accepted 3 November 2003
|
|
|---|
|
|
|---|
A number of FIV strains have been identified throughout the world, and their corresponding genomes have been partially sequenced (1, 5, 6, 10, 13, 15, 18, 19, 22, 27, 32). Studies developed on the basis of the nucleotide sequences from the envelope gene separated FIV strains into five distinct phylogenetic subtypes designated A to E (1, 2, 10, 13, 22, 27). Although geographic isolation could be a major factor in the evolution of FIV, individual subtypes are found in more than one part of the world (1, 2, 5, 13, 19, 27). For example, subtypes A and B have been isolated from cats in the United States, Europe, Japan, and Australia, and subtype C-infected cats have been identified in North America, Europe, and Taiwan. Consequently, recombinant strains have also been described (notably between subtype A and B subtypes and between B and D subtypes) (1, 2). The isolation and characterization of new FIV isolates with a broad feline representation, including high-risk feral cats, are critical for better understanding of ongoing genetic diversity.
In the present study, FIV was isolated from 11 feral cats that were trapped by the Aggie Feral Cat Alliance in Brazos County, Tex., and surrounding areas. Rigorous phylogenetic analyses were used to determine nucleotide sequences from the V3-V4 region of the envelope gene and from the whole capsid gene of the new Texas FIV strains and to compare those sequences with existing sequences of other FIV strains.
Feline peripheral blood mononuclear cells (PBMCs) were isolated from EDTA (K3)-treated whole-blood samples by Histopaque-1077 (Sigma Chemicals, St. Louis, Mo.) density gradient centrifugation. PBMCs were grown ex vivo for 7 to 10 days as previously described (27). Cells were grown at 37°C in a humidified atmosphere containing 5% CO2, with a change of medium every 3 to 4 days. FIV replication in the culture was determined by detection of capsid antigen in the supernatants with an antigen capture enzyme-linked immunosorbent assay (IDEXX, Portland, Maine) as described elsewhere (25).
Total RNA was extracted (using TRIzol reagent [Invitrogen, Carlsbad, Calif.]) from virus-infected PBMCs and stored at -70°C. Viral cDNA was synthesized using SuperScript II RNase H reverse transcriptase (Invitrogen) primed with oligo(dT). The V3-V4 sequence and capsid gene were separately amplified by nested PCR. For the V3-V4 sequence, the primary PCR used primers 6785F (5'-GCG CAA GTA GTG TGG AG-3') and 8842R (5'-GCT TCA TCA TTC CTC CTC TT-3'); the nested, secondary amplification employed primers 7316F (5'-ATA CCA AAA TGT GGA TGG TG-3') and 7866R (5'-CAA GAC CAA TTT CCA GCA AT-3'). The following two sets of primers were used for amplifying the capsid genes: 1053F (5'-GAA GAT CTC ATG GGG AAC GGA CAG GGG CGA-3') and 2403R (5'-GTC GAC TTA TAA ATC CAA TAG TCT CTC CTC-3') for the primary PCR and 1458F (5'-GAA GAT CTC CCT ATT CAA ACA GCA AAT GGA GCA-3') and 2127R (5'-CCA TCG ATG GCA AGA GTT GCA TTT TAT ATC CTG G-3') for the nested, secondary PCR. Primary reactions (50 µl) contained 100 ng of cDNA, a 100 nM concentration of each primer, 200 µM concentrations of each of the four deoxynucleoside triphosphates, and 1.5 U of TaqDNA polymerase. After incubating at 94°C for 5 min, reactions were subjected to five cycles of amplification at 94°C for 60 s, 53°C for 60 s, and 72°C for 2 min, an additional 30 cycles at 94°C for 15 s and 53°C for 15 s (with a 0.1°C increase after each cycle), and 72°C for 2 min, which was followed by 72°C for 15 min of incubation to ensure the completion of elongation. The same conditions and parameters were applied to the second set of amplifications, except that 10 µl of each of the primary PCR products was used as a template. The PCR products were purified using a Genelute PCR clean-up kit (Sigma Chemical Co.) and directly used as templates for DNA sequencing using a dye terminator method (Gene Technologies Core Facility, Texas A&M University).
Phylogenetic relationships of these new isolates were inferred from analyses of the V3-V4 and capsid sequences (7). To determine the relationships of the V3-V4 region, nucleotide sequences from 11 Texas isolates were aligned with 43 homologues retrieved (using a MacVector v7.0 program [Genetics Computer Group, Madison, Wis]) from the GenBank database. All aligned sequences were also inspected manually. Positions containing gaps or ambiguously aligned positions were removed from the data set. Finally, a total of 501 nucleotide (nt) positions from 54 taxa that can be unambiguously aligned by codons were included in the final phylogenetic analyses. Similarly, DNA sequences of capsid genes from 27 isolates (including 7 Texas isolates) were aligned; a total of 606 nt positions were used for phylogenetic reconstructions.
Phylogenetic analyses were performed using Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) approaches for both V3-V4 and capsid datasets. For BI phylogeny, nucleotide sequences were partitioned by codons. The MrBayes v3.0 program (9) was used to perform a total of 1,000,000 generations of searches with a Markov chain Monte Carlo sampling method with four chains simultaneously running. Start trees for all chains were randomly generated by the MrBayes program. Posterior possibility (PP) values were obtained using an extended majority ruling method to calculate the consensus tree from the best BI trees. ML analysis employed a quartet-puzzling method using a HKY85 nucleotide substitution model with the consideration of two rates of heterogeneity (one invariable plus one variable). Best quartet-likelihood trees were constructed using the Treepuzzle v5.0 program with 1,000 puzzling steps (30). All parameters, including transition-transversion ratio, nucleotide frequency, and fraction of invariable sites, were estimated from the datasets on the basis of neighbor-joining topologies. For MP bootstrapping analysis, a full heuristic algorithm and the PAUP* v4.0b10 program (31) were used to search for the best trees from 1,000 replicates. Sequences were added randomly, and tree-bisection reconnection was applied for branch swapping during the MP tree searches.
During BI analysis, the best-likelihood values were reached before 40,000 generations of Markov chain Monte Carlo searches for both ENV-V3-V4 and the capsid datasets, indicating that 1,000,000 generations were sufficient for the chains. The BI consensus trees inferred from all alignments clearly separated FIV viral strains into six major clusters, five of which were the previously established subtypes (i.e., A, B, C, D, and E); the sixth cluster consisted of the 11 newly cloned Texas (TX) isolates plus the previously reported Txmtex03 strain (Fig. 1 and 2). The separation of these major subtypes and the cluster was robustly supported by the PP values (i.e., PP in all major nodes = 97 to 100%). The branch lengths between the B subtype and the Texas cluster were better separated in V3-V4 trees than in the capsid gene trees. The B subtype and Texas cluster are clearly related to each other, indicating that these two clusters evolved from a common FIV viral ancestor. When the two alignments were subjected to the quartet-puzzling ML and MP analyses, both methods again divided all viral isolates into six major clusters. However, relationships of isolates among these major FIV clusters were not well resolved by quartet likelihood and MP analyses (Fig. 1 and 2; collapsed branches in ML/MP analyses are marked by an x). In all trees, the ML puzzling indices and MP bootstrapping supporting values were lower than the BI PP values at corresponding nodes with both puzzling and bootstrapping methods (indicative of rapidly evolving genomes), suggesting an overconservative nature.
![]() View larger version (29K): [in a new window] |
FIG.1. (A) Phylogenetic relationships among 54 FIV isolates determined (using BI, MP, and quartet-puzzling methods) on the basis of envelope V3-V4 region nucleotide sequences. A total of 501 nt positions were included for the tree reconstructions. The tree shown here is a consensus tree calculated from the best 9,000 BI trees. Statistical supporting values shown at major branches are PP (BI [from 9,000 best trees]), bootstrapping proportions (MP [from 1,000 replicates]), and puzzling indices (ML [with 10,000 puzzling steps]). Collapsed (unsupported) branches (as determined by MP and/or ML analysis) are marked by an x at corresponding positions. The strain names (GenBank accession numbers) for the V3-V4 sequences are as follows: EngUK8 (X69496); ScotUK2 (X69494); Wales UK14 (X69497); Wo (L06312); Dutch19K32 (M73965); USCAzepy01A (U02417); USCAhnky12A (U02403); USCAsam01A (U02410); USCAtt09A (U02413); SAP03 (AB010404); Sendai2 (D37814); Aomori 2 (D37817); Aomori 1 (D37816); Yokohama (D37812); AICO2 (AB010397); TY1 (D67064); USMAboy03B (U02419); USILbrny03B (U02418); USMOglwd03B (U02420); USTXmtex03B (U02422); M88 (Y13868); M2 (X69501); M3 (X69502); Lp9 (D84497); CABCpady02C (U02392); CABCpbar02C (U02394); CABCpbar03C (U02395); CABCpbar01C (U02393); CABCpbar07C (U02397); DEBAfredC (U57020); TI1 (AB016025); TI2 (AB016026); TI3 (AB016027); TI4 (AB016028); Fukuoka (D37815); Shizuoka (D37811); OKA01D (AB010400); Lp3 (D84496); Lp20 (D84498); Lp24 (D84500); TM2 (M59418); Usil2489_7B (U11820); Petaluma (M25381); PPR (M36968): TX125 (AY139094); TX120 (AY139095); TX200 (AY139096); TXMK (AY139097); TX133 (AY139098); TX132(AY139099); TX84 (AY139100); TXTG (AY139101); TX78 (AY139102); TX109 (AY139102); and TX77 (AY139103). (B) Consensus BI tree inferred from 27 capsid gene sequences with a total of 606 nt positions. PP values (BI [from 9,000 best trees]), bootstrapping proportions (MP [from 1,000 replicates]), and puzzling indices (ML [with 10,000 puzzling steps]) are shown at major branches. Collapsed branches (as determined by MP and/or ML analysis) are marked by an x at corresponding positions. The strain names (GenBank accession numbers) for the capsid sequences are as follows: Wo (L06136); Sendai1 (D37820); Z1 (X57001); Aomori2 (D37824); Sendai2 (D37821); Yokohama (D37819); Aomori1 (D37820); ItM3 (Y13866); ItM2 (Y13867); TI3 (AB027300); TI4 (AB027301); TI2 (AB027299); TI1 (AB027298); Shizuoka (D37822); Fukuoka (D37818); Lp24 (AB027304); Lp20 (AB027303); Lp3 (AB027302); TM2 (M59418); Usil2489_7B (U11820); Petaluma (M25381); PPR (M36968); TX120 (AY139105); TX77 (AY139106); TX78 (AY139107); TX84 (AY139108); TX53 (AY139109); TX200 (AY139110); TXTG (AY139111); and TX132 (AY139112).
|
![]() View larger version (26K): [in a new window] |
FIG. 2. Phylogenetic tree determined on the basis of the concatenated V3-V4 and capsid sequences (total positions = 1,107 nt). Statistically supporting values, which included PP values (BI [from 9,000 best trees]), bootstrapping proportions (MP [from 1,000 replicates]), and puzzling indices (ML [with 10,000 puzzling steps]), are again shown at major branches. Collapsed branches determined by MP and/or ML analysis are marked by an x at corresponding positions. The GenBank accession numbers are listed in the legend for Fig. 1.
|
The divergence of Texas strains (as a unique cluster) from subtype B viruses was also confirmed by the genetic distances of the V3-V4 regions. The means of the distances of the FIV isolates within and among various subtypes were calculated using a Jukes-Cantor method (Table 1). The results revealed that the Texas cluster was most divergent from subtype A and least divergent from subtype B. The genetic distances between sequences of isolates within the Texas group (mean of distances = 6.4) were similar to those of subtype B (i.e., 6.2), which were less than those within subtypes A and C (i.e., 10.1) but more than those within subtypes D and E (i.e., 4.4 to 4.9). Although double-crossover events could have occurred in regions not sequenced to maintain the observed coevolution of V3-V4 and the capsid sequences, genetic shifts were not identified and the similarities were maintained. Therefore, observed divergence would seem to have occurred through point mutations rather than recombination (1, 2).
|
View this table: [in a new window] |
TABLE 1. Means of Jukes-Cantor distances between isolates within and among various subtypesa
|
The divergence of FIV in the United States raises broader questions relative to vaccine design (1, 4). Although the only commercial vaccine for FIV is comprised of subtypes A and D viruses, this vaccine may not be effective against natural challenge strains (12, 29). While there are no reports of significant differences in the biological impact of FIV isolated from Texas, the phylogenetic analyses presented in this study suggested that a possible new subtype (or subtype cluster) should be taken into consideration, especially with regard to vaccine design and testing. The amount of variation within the United States warrants additional ongoing comprehensive analyses to truly understand the breadth of FIV diversity throughout the world.
We thank the Department of Veterinary Pathobiology (department head, Ann Kier) for financial support. This work was further funded by the National Institute of Allergy and Infectious Diseases (grants AI32360 [E.W.C.] and AI44594 [G.Z.]), by the Center for AIDS Research at Baylor (grant 502342), and by the Morris Animal Health Foundation (grant 96FE-09).
Present address: Human Vaccine Institute, Duke University Medical Center, Durham, NC 27710. ![]()
|
|
|---|
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»