Previous Article | Next Article ![]()
Journal of Virology, November 2003, p. 12310-12318, Vol. 77, No. 22
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.22.12310-12318.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Satish Pillai,7 Robert W. Doms,2 Preston Marx,8 and Steven M. Wolinsky1*
Departments of Medicine,1 Cell and Molecular Biology, The Feinberg School of Medicine, Northwestern University, Chicago, Illinois 60611,6 Department of Pathology, University of Pennsylvania, Philadelphia, Pennsylvania 19104,2 Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545,3 National Primate Research Center, University of Washington, Seattle, Washington 98195,4 The Laboratory of Molecular Microbiology, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, Maryland 20892,5 Departments of Zoology,7 Tropical Medicine, Tulane Regional Primate Center, Tulane University Medical Center, New Orleans, Louisiana 704338
Received 19 December 2002/ Accepted 12 August 2003
|
|
|---|
|
|
|---|
Whereas the envelope glycoproteins of the primate immunodeficiency viruses share many common features, including similar locations of cysteine residues, conserved and variable regions, and CD4 binding domains, significant differences exist, especially between SIV and HIV-1 envelope glycoproteins (29). The external envelope glycoproteins of SIV and HIV-2 have approximately 70% amino acid sequence similarity to each other, and their nucleotide sequences have less than 40% similarity to those of selected HIV-1 envelope genes (29). Nonetheless, the divergent strains of the distantly related primate lentiviruses use a chemokine receptor from the seven-transmembrane-domain G-protein-coupled receptor superfamily for virus entry, and CCR5 use is indeed almost universal (4, 28, 49). Knowledge of CCR5 use by SIV per se is restricted to receptors that were derived from humans and selected nonhuman primates, however (4, 15, 49). Therefore, we cloned and analyzed CCR5 sequences derived from 24 representative primate species and subspecies of the suborders Prosimii and Anthropoidea. Phylogenetic reconstruction of the nonhuman CCR5 homologue sequences by maximum-likelihood and neighbor-joining methods recapitulated the evolutionary relationships for these representative primate species and subspecies. Despite a number of amino acid substitutions in multiple extracellular domains of the primate CCR5 homologues, many of which were highly conserved in primates within a given species, the primate CCR5 homologues supported SIV and HIV-1 entry.
|
|
|---|
PCR amplification of nonhuman primate CCR5 homologues. We isolated total cellular DNA by using a Puregene DNA isolation kit (Gentra Systems, Minneapolis, Minn.). DNA pellets were resuspended in water, the nucleic acid content was quantified, and the DNA was used directly for in vitro amplification. To amplify a functional CCR5 gene from each species, we designed primers based on the conserved 5' (CCR5A; 5'-GGAGGGCAACTAAATACATTCTAGG-3') and 3' (CCR5B; 5'-GACTGGTCACCAGCCCACTTGAGTCC-3') untranslated regions that flank the entire human CCR5 gene (GenBank accession number U57840). We were not able to amplify the entire CCR5 gene from cells obtained from primates of the families Lorisidae, Tarsiidae, and Lemuridae, a finding that suggests a nucleotide sequence mismatch within a region that flanks the gene. We were able to amplify the CCR5 DNA of the family Lorisidae by using PCR with oligonucleotide primer pairs located within the CCR5 5' (ATGGATTATCAAGTGTCAAG) and 3' (TCACAAGCCCACAGATATTT) sequences. Five hundred nanograms of input cell DNA was used in each reaction. PCR was performed as described previously (27) with a Perkin-Elmer Cetus 9600 automated thermal cycler programmed for 5 cycles at 98°C for 10 s, 55°C for 30 s, and 72°C for 45 s and then for 35 cycles at 98°C for 10 s, 60°C for 30 s, and 72°C for 45 s. Reagent controls were run in parallel. PCR product DNAs were resolved on a 1.0% agarose gel by electrophoresis. The appropriately sized band was excised, and the DNA was extracted by using the QIAquick gel extraction method (Qiagen, Valencia, Calif.). Purified PCR product DNA was inserted into vector pCR2.1 by the TA cloning method. All polynucleotide extractions and amplification reactions were done with appropriate negative controls to detect contamination at each step of the procedure.
DNA sequencing. For each species, we obtained CCR5 sequences from at least one representative animal. Six clones containing the appropriately sized CCR5 homologue insert were randomly selected from each animal to ensure a 97% probability of analyzing both alleles [1 - 2(1 -n), where n is the number of independent clones]. The PCR product DNA inserts were sequenced in both directions by using ABI dye-terminator chemistry. Extension reactions were resolved on an ABI 377 automated sequencer, and the resulting chromatograms were assembled by using AutoAssembler (Applied BioSystems, Foster City, Calif.). Stop codons were found in two of the clones sequenced: 1 of 12 clones from two C. jacchus animals and 1 of 18 clones from six M. mulatta animals. Because these stop codons were not found in other clones, they likely represent PCR-associated errors. There were no frameshift mutations, in-frame insertions, or deletions in any of the clones sequenced. A consensus sequence was computed for each animal and was used in all DNA sequence analyses. We numbered the CCR5 homologue sequences according to the human CCR5 sequence under GenBank accession number U57840.
Phylogenetic sequence analysis of primate CCR5 homologue sequences. Comparisons were based on various breakdowns of the available PCR product DNA sequence data. Both DNA and protein alignments were hand edited by using a multiple-alignment sequence editor (18). Analysis of the number of synonymous mutations per potential synonymous site and of the number of nonsynonymous mutations per potential nonsynonymous site was performed by using the formula of Nei and Gojobori (34). Because the CCR5 nucleotide sequences of primates in the same species proved to be very similar or identical to each other, we selected one representative CCR5 sequence for analysis. We were not able to amplify the full-length CCR5 homologue sequence from primates of the taxonomic suborder Prosimii; we therefore used the 1,019 nucleotides that we could evaluate for all nonhuman primates for all phylogenetic analyses.
The CCR5 sequence alignments were subjected to three different tree-building methods: maximum parsimony, neighbor joining (a distance method), and maximum likelihood. These methods were implemented with the computer programs PAUP, version 3.1.1 (48); neighbor in PHYLIP, version 3.5c (48); and RevDNArates, a modification of DNArates code that incorporates a general reversible base substitution model, and RevDNAml, a modification of fastDNAml code that also incorporates a reversible model for rapid implementation of PHYLIP maximum-likelihood code (26, 35). We used the exhaustive branch-and-bound search method to find all possible shortest trees under parsimony. We compared 57 models of nucleotide substitution for neighbor-joining distance by using hierarchical likelihood ratio tests (22, 23). Of these 57 models, 56 are implemented in Modeltest, version 3.06 (37). The general time-reversible model with site-specific rates is implemented in PAUP, RevDNAml, and RevDNArates; in Modeltest, it is computed by hand. We tested all 57 models by using a single neighbor-joining tree based on the assumptions and methods implemented in Modeltest. This strategy includes the Jukes-Cantor model (equal base frequencies, equal rates of substitution), minimum evolution, no invariant sites, negative branch lengths set to zero, and random breaking of ties. We used RevDNArates to optimize parameters for the model with the general time-reversible model with site-specific rates that include base frequencies, substitution rates, and four site-specific categories of rates. After determining the best model, we used RevDNAml to search for the best tree. For the final tree, we used PAUP with 500 neighbor-joining search replicates to compute bootstrap percentages for each branch.
Expression vector constructs. We selected a representative CCR5 homologue insert from each species to transfer into mammalian expression vector pCDNA3.1. Briefly, clone pCR2.1 was digested with the appropriate flanking restriction enzymes to remove the CCR5 insert. The insert was inserted into pCDNA3.1 and then sequenced in its entirety to ensure that no mutations had been introduced in the cloning process. We performed in vitro transcription-translation reactions (Promega, Madison, Wis.) by using the subcloned genes to ensure that a protein product of the approximately correct molecular weight could be made from each clone. Translatable clones were used to assess the role of primate CCR5 genes as cofactors for virus entry. Luciferase reporter viruses were prepared by cotransfecting 293T cells with pNL-Luc-Env or pSIV-Luc-R-E and vectors expressing different SIV or HIV-1 envelopes as described previously (11, 14).
Primary SIV strains. Primary SIV strains (SIVsmSL92b and SIVsmFNJ) were isolated from cells obtained from two sooty mangabeys by cocultivation with uninfected cells from the native host species as described previously (4).
Virus entry assay. HOS.CD4 cells were transiently transfected with pcDNA3.1-CCR5 constructs containing the CCR5 homologues of the various primate species by using a CellPhect kit (Pharmacia Biotech). To test for virus entry, cells (2 x 104 per well) in culture medium (Dulbecco's modified Eagle's medium containing 10% fetal calf serum and 1.0 µg of puromycin/ml) were seeded in 24-well plates 24 h posttransfection. After an additional 24 h of culturing, cells were infected with luciferase reporter viruses (40 ng of p24 or p27) in a total of 300 µl of culture medium containing 1 ng of Polybrene/ml. The medium from the reporter virus-infected cells was replaced with fresh medium (1 ml per well) after overnight incubation. After an additional 3 days of culturing, we prepared 100 µl of cell lysate and measured the luciferase activity in 20 µl by using commercial reagents (Promega).
Cell-cell fusion assay. The cell-cell fusion assay has been described in detail elsewhere (41). Briefly, effector QT6 cells were infected with vaccinia virus expressing T7 polymerase (vTF1.1) and then transfected with plasmid pcDNA3.1 containing the designated envelope gene via the calcium phosphate transfection method. Target QT6 quail cells were transfected with plasmids expressing CD4, coreceptor, and T7 luciferase. At 24 h after infection or transfection, effector cells were added to target cells, and fusion was quantified 8 h after mixing by lysing cells with 0.5% Triton X-100. An aliquot of the lysate was mixed with an equal volume of luciferase assay reagent (Promega), and luminescence was read in a Luminometer (42).
Nucleotide sequence accession numbers. After rigorous quality control testing (30), the 169 human and nonhuman CCR5 homologue sequences were submitted to GenBank under accession numbers AF161887 to AF162055.
|
|
|---|
![]() View larger version (23K): [in a new window] |
FIG. 1. Phylogenetic reconstruction of human and nonhuman primate CCR5 homologue sequences by using the neighbor-joining algorithm. We used one representative CCR5 sequence per species (Table 1). The numbers are bootstrap values (percentages); each value refers to the branch to its right or left. Branches with no bootstrap values had either a bootstrap value not greater than 50% (rare) or a branch length near or equal to zero (common). Apes are highlighted in blue, Old World monkeys are highlighted in magenta, New World monkeys are highlighted in green, and prosimians are highlighted in red.
|
|
View this table: [in a new window] |
TABLE 1. List of species (or subspecies which taxonomists may treat as species) studied
|
![]() View larger version (60K): [in a new window] |
FIG. 2. Alignment of deduced human and nonhuman primate CCR5 homologue amino acid sequences. Directly below the human CCR5 consensus sequence are the nonhuman primate CCR5 homologue sequences for apes (P. troglodytes, G. gorilla, P. pygmaeus, and H. concolor), Old World monkeys (P. cynocephalus papio, M. sphinx, M. fascicularis, M. fuscata, M. mulatta, M. nemestrina, C. guereza, C. diana, C. aethiops sabaeus, C. pygerythrus, C. mona, C. nictitans, and E. patas), New World monkeys (A. caraya, A. trivirgatus, Saguinus sp., and C. jacchus), and prosimians (L. catta and V. variegata rubra). Each nonhuman primate CCR5 homologue sequence is the consensus sequence of six clones obtained from each animal to ensure a 97% probability that both alleles were represented. Amino acids in the alignment that share similarity with the human CCR5 consensus sequence shown at the top are indicated by a dash. The truncated 5' and 3' ends of the prosimian CCR5 homologues are indicated by asterisks. The amino acid alignment is annotated for the amino terminus, extracellular loops (ECL), intracellular loops (ICL), transmembrane regions (TM), and carboxy terminus. Apes are highlighted in blue, Old World monkeys are highlighted in magenta, New World monkeys are highlighted in green, and prosimians are highlighted in red.
|
Among the species of the superfamily Hominoidea, we found an additional amino acid substitution at position 9 (I to T) for H. leucogenys, G. gorilla, and P. pygmaeus but not for P. trogodytes. Among the species of the superfamily Cercopithecoidea, we commonly found amino acid substitutions in the amino terminus (position 9 [I to T]), the first intracellular loop (positions 58 [M to I] and 61 [I to V]), the second transmembrane domain (position 78 [F to L]), and the second extracellular loop (position 171 [K to R]). We found an additional conserved substitution in the second extracellular loop (position 198 [I to M]) for M. facicularis, M. fuscata, M. mulatta, M. tonkeana, and M. nemestrina. Among the species of the superfamily Ceboidea, we found amino acid substitutions in the amino terminus (positions 15 [Y to G], 16 [T to A and T to P], 24 [N to D], 27 [I to M], 28 [A to G], and 30 [R to H]), the first transmembrane domain (positions 38 [L to M and L to I] and 40 [I to L]), the first intracellular loop (position 61 [I to V]), and the second transmembrane domain (position 81 [L to F]). Additional amino acid substitutions were found in the first extracellular loop (position 92 [A to G]), the third transmembrane domain (position103 [L to F]), the second extracellular loop (positions 174 [L to Y], 180 [S to P], 184 [Y to F], 194 [Q to E], and 198 [I to M]), the sixth transmembrane domain (position 260 [F to Y]), the seventh transmembrane domain (position 292 [I to V]), and the carboxy terminus (positions 329 [E to K] and 335 [S to N]). Last, among the species of the suborder Prosimii, from which we were able to obtain representative clones of the CCR5 homologue with species-specific amino and carboxy termini, we found additional amino acid substitutions in positions throughout the molecule. The limited number of species of the suborder Prosimii precludes a formal analysis of deduced amino acid sequence conservation. Nevertheless, analysis of the number of synonymous mutations per site and the number of nonsynonymous mutations per site (24) suggests that the nonhuman CCR5 homologues underwent purifying selection, as witnessed by the high ratio of nonsynonymous to synonymous substitutions per site (average over all species, 6.62), especially in the carboxy-terminal third of the protein, the region beyond codon 200.
Diverse primate CCR5 homologues support membrane fusion by virus strains. To determine whether the diverse primate CCR5 homologues function as coreceptors for SIV, we tested the recombinant proteins in a cell-cell fusion assay by using SIVmac Env 239.2 and SIVmac Env 316. No fusion was observed when Env-negative cells were used. Figure 3 shows that expression of the diverse nonhuman primate CCR5 homologues in QT6 cells along with human CD4 allowed fusion with QT6 effector cells expressing either of the two different envelope proteins.
![]() View larger version (37K): [in a new window] |
FIG. 3. Human and nonhuman primate CCR5 homologue proteins support membrane fusion. QT6 cells expressing the SIVmac Env 239.2 and SIVmac Env 316 proteins and T7 polymerase were mixed with QT6 cells expressing CD4, luciferase under the control of the T7 promoter, and the different primate CCR5 homologues. Cell-cell fusion was measured 8 to 10 h later by quantification of luciferase activity in relative light units. All results are the averages of at least three experiments with the different CCR5 homologues. Values are normalized to the luciferase activity for the M. fascicularis CCR5 homologue. Apes are highlighted in blue, Old World monkeys are highlighted in magenta, and New World monkeys are highlighted in green.
|
|
|
|---|
Alignment with and comparison to the human CCR5 amino acid sequence of the deduced nonhuman primate CCR5 homologue amino acid sequences showed species-specific changes among the diverse species, with greater amino acid differences for prosimians than for Old World monkeys, New World monkeys, and apes. Specifically, we found conserved amino acid substitutions at positions 33 (N to D) and 129 (V to I) among all nonhuman primate CCR5 homologue sequences relative to the human CCR5 consensus sequence. Additionally, we found highly conserved amino acid substitutions at position 9 (I to T) for apes and at positions 9 (I to T), 58 (M to I), 61 (I to V), 78 (F to L), and 171 (K to R) for Old World monkeys and an additional conserved substitution at position 198 (I to M) for macaques. Similar amino acid changes were observed in limited data sets of nonhuman primate CCR5 homologue sequences derived from representative species of the superfamilies Ceboidea and Cercopithecoidea (33). We observed larger numbers of relatively conserved amino acid substitutions among New World monkeys and prosimians.
CCR5 polymorphisms were also found within a given species, suggesting that, like humans (3, 38), nonhuman primates have polymorphisms (5, 28). Most of these chemokine receptors functioned as coreceptors for SIV and HIV-1 entry. Other than the stop codon found in two alleles, we identified only one mutation, an N13D change for the vervet CCR5 homologue, that adversely affected virus entry (17). Although we did not find any deletions within the nonhuman primate CCR5 homologue coding regions sequenced here, a 32-bp deletion allele exists in the human population (31, 38) and a 24-bp deletion allele exists in the C. torquatus torquatus population (36). Such polymorphisms may play a role in the limited pathogenicity of SIV in naturally and experimentally infected C. atys hosts.
Each of the nonhuman primate CCR5 homologues tested here permitted the entry of genetically diverse SIV and HIV-1 Env constructs, including SIVmac1A11 and HIV-1 laboratory strains JR-FL, ADA, and HXB2. Thus, whereas simian-human immunodeficiency virus constructs allow HIV-1 envelope fusion and entry in primate cells expressing CD4, the SIV accessory genes overcome the restrictions to replication at the postentry level. These findings indicate that the nonhuman primate CCR5 homologue is not a species-specific restriction for primate immunodeficiency virus entry and support the hypothesis that HIV-1 originated as a cross-species transmission event (4, 20).
The topology of the phylogenetic tree based on the nonhuman primate CCR5 homologues sequenced here agrees with that constructed by others using limited numbers of nonhuman primate CCR5 sequences from phylogenetically distinct species of Anthropoidea and Hominoidea (33, 43, 46). We generally found excellent agreement between the phylogenetic reconstructions of nonhuman primate CCR5 homologue sequences and available COii and CD4 gene sequences (12, 21, 52) created by using the maximum-likelihood method (data not shown). The branching patterns of the nonhuman primate CCR5 homologue sequence sets estimated by the maximum-likelihood and neighbor-joining methods revealed distinct genetic lineages, with quantitatively greater branch lengths for prosimians than for Old World monkeys, New World monkeys, and apes. The branching order was also consistent with the presumed historical diversification of the nonhuman primates (21, 46, 52). Specifically, closer linkages were observed among species representing Ceboidea, Cercopithecoidea, and Hominoidea. The only difference between the phylogenetic reconstructions created by using the maximum-likelihood and neighbor-joining methods was the location of Colobus, which clustered in an ancestral location relative to the C. aethiops cluster in the maximum-likelihood tree but in an ancestral location relative to all cercopithecines in the neighbor-joining tree (data not shown). This discrepancy, noted by others (52), suggests that the nonhuman primate CCR5 homologue may not be a suitable marker for low-level phylogenetic analysis. Analysis of the number of synonymous mutations per site and the number of nonsynonymous mutations per site (24, 34) suggests that the nonhuman primate CCR5 homologues underwent purifying selection. Furthermore, a difference in the number of synonymous mutations per site and the number of nonsynonymous mutations per site for cercopithecines and Colobus implies a difference in selective pressures for these two lineages (52).
In summary, nonhuman primate CCR5 homologue sequences are highly homologous, with amino acid sequence variations in the amino-terminal and second extracellular domains showing evolutionary conservation among the different species. Except for an N13D change for the C. aethiops CCR5 homologue (17), all other nonhuman primate CCR5 homologues supported SIV and HIV-1 entry. The topologies of the nonhuman primate CCR5 homologue phylogenetic trees recapitulated the presumed evolutionary relationship among the nonhuman primates, suggesting that this phylogenetic information is useful for resolving the relationships among the taxonomic suborders.
Support for this work was provided by grants from the Public Health Service (HD37356, AI41420, and AI40880), the Elizabeth Glaser Pediatric AIDS Foundation, and the Burroughs Wellcome Fund and by a gift from an anonymous foundation.
Present address: Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 06511. ![]()
|
|
|---|
ccr5 mutation. J. Virol. 71:3219-3227.[Abstract]
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»