Previous Article | Next Article ![]()
Journal of Virology, January 2002, p. 259-268, Vol. 76, No. 1
0022-538X/01/$04.00+0 DOI: 10.1128/JVI.76.1.259-268.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Unité dEpidémiologie et Physiopathologie des Virus Oncogènes, Département du SIDA et des Rétrovirus, Institut Pasteur, 75724 Paris Cedex 15, France,1 Centre Pasteur du Cameroun, Yaoundé, Cameroon,2 International Zoo Veterinary Group, Keigthley, West Yorkshire, United Kingdom3
Received 11 July 2001/ Accepted 25 September 2001
|
|
|---|
|
|
|---|
Most of the viruses belonging to the PTLV-1 lineage, which comprises HTLV-1 and STLV-1, cannot be separated into distinct phylogenetic lineages according to their species of origin. Their phylogenetic intermixing has been interpreted as evidence for past and recent interspecies transmission episodes (5, 6, 24, 29, 34). However, with regard to the viral transmission from monkeys to humans, this hypothesis is supported by a still-limited but increasing number of observations (13, 14, 17, 35, 3739). The situation for PTLV-2, which comprises HTLV-2 and STLV-2, is different, since HTLV-2 and STLV-2, while clustering in the same large phylogenetic clade, are distantly related, with no evidence for recent interspecies transmissions (32). Regarding PTLV-L, the only known strain (STLV-L/PH-969) was isolated in 1994 from an Eritrean baboon (Papio hamadryas) kept in a captive colony in Leuven, Belgium. STLV-3/PH-969 (formerly STLV-L/PH-969), which remains the unique prototype of its type, exhibits 40 and 38% divergence at the nucleotide level from HTLV-1 and HTLV-2, respectively (30). The goals of this study were to search for highly divergent PTLV strains among African primates in order to gain new insights into the origin, evolution, and modes of dissemination of such viruses and their hosts.
We report here the isolation, molecular characterization (complete nucleotide sequence), and phylogenetic analysis of a novel STLV subtype infecting two wild-caught red-capped mangabeys (Cercocebus torquatus torquatus) (CTO-604 and CTO-602) originating from southern Cameroon. These viral strains, named STLV-3/CTO-604 and STLV-3/CTO-602, which elicit in their host an HTLV-2-like serology as determined by Western blotting, are genetically highly divergent from the PTLV-1 and PTLV-2 strains and distantly related to STLV-3/PH-969. The finding of this viral subtype in Cameroon greatly enlarges the geographical distribution of this PTLV type in the African continent. In addition, the presence of two highly divergent (compared to PTLV-1 and PTLV-2) but related viruses in two evolutionarily distantly related African monkeys species (i.e., mangabeys and Eritrean baboons) living in two opposite ecosystems (rain forest versus desert) reinforces the possible African origin of PTLV. It also opens new avenues regarding the search for a possible human counterpart of these viruses in individuals exhibiting such HTLV-2-like seroreactivities, especially in the African continent.
|
|
|---|
Serological tests. Serum or plasma samples were screened for the presence of HTLV and STLV antibodies by enzyme-linked immunosorbent assay (Sanofi Diagnostics Pasteur, France) and by particle agglutination (PA) tests (Fujirebio Japan), as well as with an in-house immunofluorescence assay (IFA) using HTLV-1-producing (MT2) or HTLV-2-producing (C19) cell lines (16). IFA and PA tests were also used to determine the titers of HTLV and STLV antibodies. All positive or borderline samples were then tested with a Western blot assay (HTLV 2.4; Diagnostic Biotechnology, Singapore), which contains HTLV-1-purified virions enriched with a gp21 recombinant protein (RGD21) that reacts with sera containing HTLV-1 or HTLV-2 antibodies and two gp46 Env synthetic peptides specific either for HTLV-1 (MTA1) or HTLV-2 (K55) (36).
Cell culture and virus isolation. Heparinized blood specimens were drawn from the two animals (CTO-602 and CTO-604) identified as having HTLV-2-like seropositivity in the pilot survey and then rushed to our unit, where the peripheral blood mononuclear cells (PBMCs) were separated with Ficoll-Hypaque (Eurobio, Les Ulis, France). Ten million cells were placed in culture in RPMI 1640 medium with 20% heat-inactivated fetal calf serum, 1% L-glutamine, and 1% penicillin-streptomycin. During the subsequent 3 days, the cells were stimulated with phytohemagglutinin (Difco, Detroit, Mich.) at 2 µg/106 cells. The cells were then cultivated in a humidified 5% CO2 atmosphere in the same medium as described above in the presence of 10% interleukin-2 (Boehringer, Mannheim, Germany); the medium was changed twice per week.
Indirect immunofluorescence and antigen detection. Indirect immunofluorescence was performed on cultured cells in order to detect the expression of viral antigens, using either mouse monoclonal antibodies directed against HTLV-1 p19 or p24 (Cambridge Biotech, Cambridge, Mass.) or human HTLV-1 and HTLV-2 polyclonal sera. Sera from the two animals studied (CTO/604 and CTO/602) were also used.
PCR. PCR was implemented using previously described conditions (7) on high-molecular-weight DNA extracted from PBMCs before and after 2 months of culture. A first round of PCR was performed with the primer sets GAG1-GAG2 (gag region of HTLV-1), SK110-SK111 (pol region of HTLV-1 and HTLV-2), and KKPX1-KKPX2, SK43-SK44, and TR101-TR102 (tax region of HTLV-1 and HTLV-2) as described previously (17, 21). In order to obtain the complete sequence of this novel isolate, successive PCRs were performed using nine different primers sets designated A to I (Fig. 1). For each set, the first primer was specific for the previously obtained sequence and the second primer was a consensus sequence of all previously known HTLV and STLV prototype strains (Fig. 1 and Table 1). For all of the PCRs, the amplification conditions were as follows: denaturation at 94°C for 9 min and then 40 cycles of 94°C for 30 s, annealing at a primer-specific temperature (Table 1) for 30 s, and extension at 72°C for 30 s per 500 bp. An extension of 10 min was performed after the last cycle. Reaction tubes, prepared in a room physically separate from the laboratory, contained 1 µg of DNA, 0.2 mM deoxynucleoside triphosphate mix (Boehringer), 5 µl of 10x reaction buffer, 1.5 to 2.5 mM MgCl2, and 2.5 U of Taq Gold DNA polymerase (Perkin-Elmer) in a total volume of 50 µl. All PCR products were purified on a 1% agarose gel by gel extraction using the QIAquick gel extraction kit (Qiagen GmbH, Hilden, Germany). Purified DNA was then cloned in the pCR2.1 vector of the TA cloning kit (Invitrogen, Carlsbad, Calif.), sequenced, and verified on both DNA strains.
![]() View larger version (17K): [in a new window] |
FIG. 1. (A) PCR strategy for amplifying the entire STLV-3/CTO-604 proviral genome. The nine proviral fragments which were amplified by PCR, cloned, and sequenced are shown (black bars). (B) Schematic representation of the STLV-3/CTO-604 proviral genome (top) and of the resulting viral messengers (bottom). The start codon used for the translation of the precursor protein (asterisks), the primers used for detection of singly spliced or doubly spliced messengers (horizontal arrows), and the positions and designations of the spliced sites (vertical arrows) identified in the STLV-3/CTO-604 genome are indicated. Nucleotide numbering is according to the STLV-3/CTO604 proviral genome. sa, splice acceptor; sd, splice donor.
|
|
View this table: [in a new window] |
TABLE 1. Sequences of primer sets used for amplifying the complete proviral genome of STLV-3/CTO-604
|
Phylogenetic analyses. Multiple nucleotide and amino acid sequence alignments were performed with the ClustalW algorithm implemented in MacVector 6.5 (Oxford Molecular). The 16 PTLV Tax amino acid sequences (331 amino acids [aa]; MO Tax aa 1 to 331) were aligned using a PAM 250 matrix, and the 16 Env polyprotein PTLV nucleotide sequences were aligned to the ATK Env polyprotein sequence (nucleotides [nt] 1 to 1467). Phylogenetic analyses were performed using the PHYLIP package with two different methods: neighbor joining (NJ) and maximum parsimony (MP). The SEQBOOT program generated 1,000 data sets that are randomly resampled versions of the aligned sequences. A distance matrix was calculated for each data set using the PROTDIST and DNADIST programs with the Kimura two-parameter model, and an empirical transition/transversion ratio for the Env polyprotein (2.23) was used. This ratio was estimated from the data set with the Treepuzzle 5.0 program. The NEIGHBOR program generated a tree for each data set, and a consensus tree was constructed by using the CONSENSE program with the majority-rule criteria. The same data sets were examined with the PROTPARS and DNAPARS programs, based on the MP method, with the same parameters used in NJ to test the robustness of the phylogeny.
Nucleotide sequence accession numbers. The GenBank accession numbers for the sequences determined in this study are AF391797 for the complete STLV-3/CTO-604 sequence and AF391796 for the complete tax sequence of STLV-3 CTO/602.
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 2. Western blot serological patterns of the two C. torquatus animals from Cameroon, using the Western blot from Diagnostic Biotechnology (HTLV blot version 2.4). Lane 1, HTLV-1 positive control; lane 2, HTLV-2 positive control; lane 3, HTLV-1 and -2 negative control; lane 4, plasma from CTO-604; lane 5, plasma from CTO-602.
|
In order to confirm the presence of HTLV-related viruses in these two animals, high-molecular-weight DNA, extracted from uncultured and cultured PBMCs, was subjected to PCR screening with gag, pol, and tax primer pairs. The amplified products were analyzed by Southern blotting with HTLV-1-specific (gag, pol, and tax), HTLV-2-specific (pol), or both HTLV-1- and -2-specific (tax)
-32P-radiolabeled probes. The CTO-602 and CTO-604 DNA samples scored positive only with a highly conserved tax probe (SK45) that can detect both HTLV-1 and HTLV-2. Such data strongly suggested that the two red-capped mangabeys were infected with a highly divergent STLV strain. Using highly conserved primers (KKPX1-MacArc4S) (18), a 433-bp sequence (tax region) was obtained and found to be nearly identical (only one nucleotide difference) for the four uncultured and cultured CTO-602 and CTO-604 samples. Comparison with all of the available PTLV prototype sequences indicated that this 433-bp fragment exhibited a high similarity (95%) to the STLV-3/PH-969 sequence, while it was much more divergent but roughly equidistant from all of the other PTLV prototypes, including HTLV-1 ATK (78%), HTLV-2 MO (77%), and STLV-2 PP1664 (78%).
The complete sequence of one of these two strains (STLV-3/CTO-604) was then obtained by cloning and sequencing nine proviral fragments amplified by successive PCR on cellular DNA from cultured PBMCs of the CTO-604 animal (Fig. 1). The comparison of this complete sequence (8,919 bp) with the other PTLV prototypes indicated that although it was unique, this new strain was more related to STLV-3/PH-969 (87% similarity for the complete sequence) than to HTLV-1 (60%), HTLV-2 (62%), or STLV-2 (62%) prototypes (Table 2). Furthermore, the complete STLV-3/CTO-602 tax cDNA sequence obtained from short-term-cultured PBMCs demonstrated that these two novel strains were almost identical (99.9% for a 1,053-bp tax sequence).
|
View this table: [in a new window] |
TABLE 2. Nucleotide sequence comparisons between the PTLV prototype subtype strains and STLV-3/CTO-604
|
Interestingly, the long terminal repeat (LTR) of STLV-3/CTO-604 has high similarity (86.2%) with the STLV-3/PH-969 LTR and a similar length (694 and 695 bp, respectively). Its overall organization is also identical to the one present in the HTLV-bovine leukemia virus genus, with highly conserved regions (positions of U3-R and R-U5 boundaries, TATA box, polyadenylation site and signal, and potential splice donor site). Interestingly, the STLV-3/CTO-604 LTR was shorter than the HTLV-1 (756 bp) and HTLV-2 (764 bp) LTRs. This was mainly due to the presence of only two 21-bp repeats (the middle and the proximal ones) in the STLV-3/CTO-604 U3 region. This was also previously reported for STLV-3/PH-969 (30).
Although the STLV-3/CTO-604 sequence was related to STLV-3/PH-969, the overall divergence between the two strains in the major ORFs was equivalent to or higher than (8 to 13%) those between the different subtypes of HTLV-1 (A, B, C, and D) or HTLV-2 (A, B, and D) (7 to 9%) (Table 2).
A comparison of the protein sequences of the different PTLV prototypes is given in Table 3. p19, p24, and p15, in analogy to the HTLV-1 proteins, were defined in the 422-aa Gag precursor. The p24 protein is the most conserved among the different types (around 85% similarity) and subtypes (range, 96 to 98%). This is in agreement with the strong reactivities observed against the conformational p24 epitope by Western blotting. The highest divergence among the Gag proteins is observed in p19, with an important deletion in the carboxy-terminal part of the protein in the new STLV-3/CTO-604. Thus, the p19 immunodominant epitope in HTLV-1 ATK (residues 102 to 117, PPSSPTHDPPDSDPQI) shows only 18% amino acid similarity with STLV-3/CTO-604.
|
View this table: [in a new window] |
TABLE 3. Protein sequence comparisons between the PTLV prototype subtype strains and STLV-3/CTO-604
|
Although previous studies (30) have shown the existence of four putative ORFs in the pX region of STLV3-PH969, only three alternatively spliced messengers could be detected by RT-PCR. These doubly spliced messengers encoded the Tax and Rex proteins and a putative protein called RORFI. This protein is related at the amino acid level to p12I and p10I of HTLV-1 and HTLV-2, respectively. The sequence analyses of the new STLV-3/CTO-604 and CTO-602 Px regions show only two ORFs that correspond to Tax and Rex ORFs (Table 3). Indeed, regarding the possible equivalent of the STLV-3/PH-969 RORFI sequence in the novel strain, there are three mutations that replace two glutamines (CAG) with stop codons (TAG and TAA), and another mutation in the splice acceptor site replacing AG with AA, that probably eliminate the splice junction.
To gain new insights into the genomic organization of the Px region and to search for the presence of singly or doubly spliced viral messengers previously described for the other HTLV and STLV retroviruses, we performed a series of RT-PCR experiments on the total RNA extracted from a 2-month culture of the CTO-602 PBMCs. RT-PCR using primer pair 602 RT LTR and 602 RT env (Fig. 1) amplified a cDNA (950 bp) that encodes the major Env polyprotein, consisting of exon 1 (splice donor [sd]-LTR at nt 409) spliced directly to exon 2 (splice acceptor [sa]-env at nt 4903). A doubly spliced mRNA (424 bp long) was also detected using primer pair 602 RT LTR and MVB Rex ex2 with a second splicing event between the exon 2 env splice donor (sd-env nt 5070) and exon 3 Tax-Rex splice acceptor (sa-T/R at nt 7249). The splice junction between exon 2 and exon 3 is located at a nucleotide position very similar to that found in STLV-3/PH-969, while the one located between exon 1 and exon 2 occurs 16 nt downstream from the equivalent splice of the STLV-3/PH-969 strain. As previously described for HTLV-1 (3), we found at least one alternative splice acceptor for exon 2 (sa-env2 nt 4959) that could be used to generate the env and all doubly spliced mRNAs.
To detect other potential doubly spliced viral messengers, we hybridized the RT-PCR products with two specific
-32P-labeled probes, one in the LTR (604 RT 328S) and the other in the Px proximal region (604 RT ORF1). None of the several sequenced clones that doubly hybridized correspond to a new cDNA.
Although the MVB Rex ex2 primer could amplify all of the potential singly and doubly spliced mRNAs, our method probably lacks sensitivity. This results in the preferential amplification of the Tax-Rex doubly spliced mRNA. To improve our RT-PCR sensitivity, we used two primers specific for the proximal Px region: 602 RT-ORFI (antisense, 6900-GCTAAGCTATTGGCGAGAGCG-6880) and 602 RT Px2 (antisense, 7025-GATCAGGTGGACATGCTTCAG-7005). We did not obtain any single or doubly spliced mRNAs.
Phylogenetic analyses were performed with the STLV-3/CTO-604 strain compared to representatives of the PTLV-1 and PTLV-2 subtypes and to STLV-3/PH-969. Due to the high sequence divergence observed between the PTLV sequences in noncoding regions, we aligned both the Tax amino acid sequence and the Env polyprotein nucleotide sequence rather than the complete genomes. Unrooted trees generated three distinct PTLV groups, supported by high bootstrap values (100% by NJ and MP), in both env nucleotide and Tax amino acid analyses (Fig. 3 and 4). Two of them, i.e., PTLV-1 and PTLV-2, comprised human and simian viral strains, while the third one included only the new simian strains (STLV-3/CTO-604 and STLV-3/CTO-602) described in this study and the previously described STLV-L/PH-969. Phylogenetic clustering of all of these strains in three clearly distinct groups strongly suggests a long independent evolution of these viruses.
![]() View larger version (17K): [in a new window] |
FIG. 3. Unrooted phylogenetic tree generated by the NJ method on the Tax amino acid sequences (aa 1 to 331 of the HTLV-2A MO prototype sequence). Bootstrap support (1,000 replicates) for the NJ tree is noted on the branches of the tree. The STLV-3/CTO-604 and STLV-3/CTO-602 sequences were analyzed with HTLV and STLV prototype sequences available from the GenBank database (26). Branch lengths are proportional to the evolutionary distance (scale bar) between the taxa.
|
![]() View larger version (18K): [in a new window] |
FIG. 4. Unrooted phylogenetic tree generated by the NJ method on the Env polyprotein nucleotide sequence (nt 1 to 1467 of the HTLV-1A ATK prototype sequence). Bootstrap support (1,000 replicates) for the NJ tree is noted on the branches of the tree. The STLV-3/CTO-604 strain was analyzed with HTLV and STLV prototype sequences available from the GenBank database (26). Branch lengths are proportional to the evolutionary distance (scale bar) between the taxa.
|
|
|
|---|
Originally, this third type of PTLV was provisionally designated L, which stands for Leuven (9). From now on, we propose to designate this third type PTLV-3. Thus, PTLV-3 would provisionally comprise two subtypes designated STLV-3 subtype A, with STLV-3/PH-969 as the prototype, and STLV-3 subtype B, with STLV-3/CTO-604 as the prototype and STLV-3/CTO-602 as a second strain of subtype B.
Recent data (26) based on molecular clock analysis of PTLV evolution, using the third codon position, estimated the evolutionary rate of PTLV to be not higher than 1.67 x 10-6 ± 0.17 x 10-6 nucleotide substitution per site per year. The clock was calibrated by employing the earliest human migration from Asia to Melanesia 60,000 years ago as the lower limit for the node separating HTLV-1 Melanesian subtype C from the other HTLV-1 subtypes. The separation between the ancient African HTLV-2 subtype D and HTLV-2 subtypes A and B was also estimated to occur around 58,000 years ago, and that between the STLV-3/PH969 and PTLV-2 strains was estimated to occur around 1,026,000 ± 110,000 years ago. Assuming that PTLV-3 evolves at the same evolutionary rate as other PTLVs, we tested the molecular clock hypothesis for PTLV (without the IIA and IIB drug user strains) on the third codon position of the Env polyprotein and Env gp21. In both cases the clock is indeed valid, and we calibrated the evolutionary rate at around 1.9 x 10-6 ± 0.1 x 10-6 nucleotide substitution per site per year. Using this value, we estimated that the separation between STLV-3/PH-969 and STLV-3/CTO-604 occurred 200,000 ± 30,000 years ago. These provisional values would probably be modified with the characterization of other PTLV-3 strains.
At the molecular level, HTLV-1 and HTLV-2 promoters consist of three 21-bp repeated elements which contain a core element essential for the LTR transcriptional activation (28). The analysis of the STLV-3/CTO-604 LTR showed a U3 region smaller than that of HTLV-1 or HTLV-2, due to the deletion of the TATA-distal 21-bp repeat. These HTLV-1 and -2 promoter elements play an important role in basal transcription in absence of Tax, but mutating one of them unequally reduces the basal HTLV-1 transcription level (1). Although only two 21-bp repeats seem to be sufficient for a high level of Tax activation in HTLV-1-infected cells (2), the absence of one 21-bp repeat may suggest that the STLV-3/CTO-604 LTR transactivation mediated by Tax could be different than that in HTLV-1 or -2. Recent studies demonstrated that HTLV-1 and HTLV-2 subtype A Tax proteins repress the p53 transcriptional activity differently (19). These data suggest that the transactivations mediated by different PTLV Tax proteins are nonequivalent. Therefore, the potential transformant capacity of STLV-3/CTO-604 Tax is now under investigation.
These two new strains (STLV-3/CTO-604 and CTO-602) were present in two red-capped mangabeys (C. torquatus torquatus), small primates whose habitat is mostly restricted to valley forests and swamps of the western part of Central Africa (12). When we initiated the study, we did not detect any other monkeys (all of wild-caught origin) infected by such a variant strain in the wildlife center studied. Indeed, among the 9 HTLV- or STLV-seropositive monkeys out the 65 animals tested, 7 were infected by a classical African STLV-1 strain (15, 21), while the two remaining HTLV-seropositive animals were CTO-602 and the CTO-604. These data suggest that these novel viruses naturally infect C. torquatus torquatus in the wild. Furthermore, to our knowledge there are no other data reporting HTLV or STLV infection in C. torquatus torquatus.
HTLV-2-like serology (strong p24 associated with RGD21 and low or no p19 seroreactivities) seems to be very rare in African monkeys. It has been reported for only a few Pan paniscus animals infected by STLV-2 strains (8, 33), in rare Papio hamadryas animals infected by STLV-L strains (9), and in one Papio anubis animal (Bab 503) infected with an African STLV-1 strain (15). It is worthwhile to note, however, that most of the seroepidemiological surveys conducted with monkeys cannot be considered representative of the situation in the wild, since most of the samples were obtained from captive animals. The biodiversity of such viruses in the wild, especially in central Africa, is thus far from being known. The recent findings of several new and highly divergent simian immunodeficiency viruses illustrate this fact well (4, 20, 40).
Multiple episodes of interspecies transmission of PTLV-1 (STLV-1 or HTLV-1) have occurred between different primates, including humans, in Central Africa (13, 14, 17, 34). It is thus tempting to speculate that some other STLV-3 strains or related viruses may exist in other monkeys species but also that HTLV strains related to STLV-3 may exist in human populations living in such areas. Searches for a possible human counterpart (possible HTLV-3) using specific consensus or degenerated primers are ongoing, especially in individuals living in the same rain forest areas as monkeys infected by STLV-3 and whose sera exhibit these HTLV-2-like seroreactivities.
Regarding the name of the new PTLV subtype described in this paper, we have tentatively and provisionally named it STLV-3 subtype B. However, among the specialists in the field, a new proposal for PTLV nomenclature is being discussed and debated in order to clarify the situation regarding these new strains and other recent findings (W. M. Switzer, V. Shanmugam, S. Van Dooren, A.-M. Vandamme, V. Bhullar, B. Parekh, and W. Heneine, Abstr. 10th Int. Conf. Hum. Retroviruses [AIDS Res. Hum. Retroviruses, Suppl. 17], abstr. O-26, 2001; S. Van Dooren, X. Pourrut, M. Peeters, E. Delaporte, and A.-.M Vandamme, Abstr. 10th Int. Conf. Hum. Retroviruses [AIDS Res. Hum. Retroviruses, Suppl. 17], abstr. P-045, 2001; T. Takemura, M. Yamashita, M. K. Shimada, T. Shotake, T. Miura, and M. Hayami, Abstr. 10th Int. Conf. Hum. Retroviruses [AIDS Res. Hum. Retroviruses, Suppl. 17], P-Add-10, 2001). When new names are approved by the consensus of such specialist groups, we will, of course, modify the names of these new PTLVs in our papers.
We acknowledge the Pandrillus (Cameroon/Nigeria) directors, Elizabeth Gadsby and Peter Jenkins, for their great help in obtaining the blood samples studied. We thank Marco Salemi for helpful discussion.
|
|
|---|
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»