Genetic characterization of the complete genome of a highly divergent simian T-lymphotropic virus (STLV) type 3 from a wild Cercopithecus mona monkey
© Sintasath et al; licensee BioMed Central Ltd. 2009
Received: 17 August 2009
Accepted: 27 October 2009
Published: 27 October 2009
The recent discoveries of novel human T-lymphotropic virus type 3 (HTLV-3) and highly divergent simian T-lymphotropic virus type 3 (STLV-3) subtype D viruses from two different monkey species in southern Cameroon suggest that the diversity and cross-species transmission of these retroviruses are much greater than currently appreciated.
We describe here the first full-length sequence of a highly divergent STLV-3d(Cmo8699AB) virus obtained by PCR-based genome walking using DNA from two dried blood spots (DBS) collected from a wild-caught Cercopithecus mona monkey. The genome of STLV-3d(Cmo8699AB) is 8913-bp long and shares only 77% identity to other PTLV-3s. Phylogenetic analyses using Bayesian and maximum likelihood inference clearly show that this highly divergent virus forms an independent lineage with high posterior probability and bootstrap support within the diversity of PTLV-3. Molecular dating of concatenated gag-pol-env-tax sequences inferred a divergence date of about 115,117 years ago for STLV-3d(Cmo8699AB) indicating an ancient origin for this newly identified lineage. Major structural, enzymatic, and regulatory gene regions of STLV-3d(Cmo8699AB) are intact and suggest viral replication and a predicted pathogenic potential comparable to other PTLV-3s.
When taken together, the inferred ancient origin of STLV-3d(Cmo8699AB), the presence of this highly divergent virus in two primate species from the same geographical region, and the ease with which STLVs can be transmitted across species boundaries all suggest that STLV-3d may be more prevalent and widespread. Given the high human exposure to nonhuman primates in this region and the unknown pathogenicity of this divergent PTLV-3, increased surveillance and expanded prevention activities are necessary. Our ability to obtain the complete viral genome from DBS also highlights further the utility of this method for molecular-based epidemiologic studies.
Simian and human T-lymphotropic viruses (STLV and HTLV, respectively) are diverse deltaretroviruses now consisting of four broad primate T-lymphotropic virus (PTLV) groups. PTLV-1, PTLV-2 and PTLV-3 include human (HTLV-1, HTLV-2, and HTLV-3) and simian (STLV-1, STLV-2, and STLV-3) viruses, respectively [1–8]. To date, a total of three individuals from southern Cameroon with reported nonhuman primate (NHP) exposures were found to be infected with the recently identified HTLV-3 [1, 7, 8]. PTLV-4 consists of only HTLV-4 which was reported from one individual in Cameroon with known exposure to NHPs . A simian counterpart of this virus has yet to be identified. Moreover, recent phylogenetic analyses of a highly divergent STLV-1-like virus from a captive Macaca arctoides suggest the possibility of a fifth group, tentatively referred to as PTLV-5 . There is currently no evidence that STLV-5 has crossed into humans. These recent discoveries of novel HTLVs and STLVs suggest a greater diversity of PTLVs than is currently appreciated.
Both HTLV-1 and -2 have spread globally and are pathogenic human viruses [10–13]. HTLV-1 causes adult T-cell leukemia/lymphoma (ATL), HTLV-1 associated myelopathy/tropical spastic paraparesis (HAM/TSP), and other inflammatory diseases in less than 5% of those infected [2, 11, 13]. HTLV-2 is less pathogenic than HTLV-1, but has been associated with a neurologic disease similar to HAM/TSP [10, 12]. The recently discovered HTLV-3 and HTLV-4 viruses have not yet been associated with any diseases, but molecular analyses of the full-length genomes have identified functional motifs important for viral expression and possibly oncogenesis [14, 15].
STLVs have been identified in diverse Old World monkeys and apes. STLV-1 has been found in at least 20 different Old World primate species in Africa and Asia, and phylogenetic analysis shows that STLV-1s cluster by geography rather than by host species suggesting they are easily transmitted among NHPs [2, 3, 5, 16, 17]. There are currently seven recognized PTLV-1 subtypes (A to G) that are comprised of genetically related HTLV-1 and STLV-1 strains from different primate species. The close relatedness and clustering of the various HTLV-1s and STLV-1s into distinct subtypes suggests that at least seven independent cross-species transmission events formed the genetic diversity of HTLV-1. Currently STLV-2 is comprised of only two strains, STLV-2(PP1664) and STLV-2(PanP), both of which were identified in two different troops of captive bonobos (Pan paniscus) .
Like STLV-1, STLV-3 has a wide geographic distribution amongst NHPs in Africa [18–27]. Because of the phylogeographical clustering of STLV-3 into distinct clades, four separate molecular subtypes have been proposed: East African (subtype A), West and Central African (subtype B), and West African (subtype C and D) clades . STLV-3 infection has been identified in captive Ethiopian gelada baboons (Theropithecus gelada) , wild sacred baboons (Papio hamadryas) , wild hybrid baboons (P. hamadryas X P. anubis hybrid) [25, 27], and captive Eritrean hamadryas baboons (P. hamadryas) , which together comprise the STLV-3 East African (subtype A) clade. The STLV-3 West and Central African (subtype B) clade is made up strains found among Senegalese olive baboons (P. papio) , Cameroonian and Nigerian red-capped mangabeys (Cercocebus torquatus torquatus), and Cameroonian agile mangabeys (Cercocebus agilis) [18, 22, 23]. Somewhat divergent subtype B STLV-3s have also been recently identified in grey- cheeked mangabeys (Lophocebus albigena) and moustached monkeys (Cercopithecus cephus) in Cameroon although the phylogeny of these viruses was inferred using relatively short tax and LTR sequences [20, 24]. That all three HTLV-3 strains which have been recently discovered in Cameroon [1, 7, 8] cluster within the STLV-3 subtype B clade is of phylogenetic significance. STLV-3 subtype C consists of divergent viruses found in Cameroonian spot-nosed guenons (Cercopithecus nictitans) though phylogenetic inference of this particular clade is limited by analysis of only very short tax-rex sequences [20, 26]. Full-length genomes of STLV-3 subtype C are currently not available. More recently, we identified a highly divergent STLV-3 strain in Cameroon from two different primate species, C. mona (Cmo8699AB) and C. nictitans (Cni78676AB) . Based on preliminary analysis of partial gene regions, these new STLVs formed a possible fourth STLV-3 lineage outside all PTLV-3 subtypes but within the diversity of the PTLV-3 group that we tentatively called STLV-3 subtype D . Both STLV-3(Cmo8699AB) and STLV-3(Cni7867AB) share 99% sequence homology in the pol, tax, and LTR regions and cluster together with high bootstrap support within the STLV-3 subtype D clade . Together, these findings demonstrate the broad range of NHP host species susceptible to STLV infection and that STLV diversity is driven more by phylogeography than by co-divergence with host species, illustrating the ease with which STLV is transmitted across species barriers [28, 29].
Here, we report the first full-length genome sequence of STLV-3(Cmo8699AB) from a wild C. mona monkey. We confirm that this virus is a highly divergent and novel STLV-3. Across the genome, we found evidence that STLV-3d(Cmo8699AB) is unique from other PTLVs. Robust phylogenetic analysis of major gene regions of STLV-3d(Cmo8699AB) as well as new tax sequences from the divergent STLV-3d(Cni3034) and STLV-3d(Cni3038) viruses demonstrate that STLV-3d(Cmo8699AB) is a novel and ancient lineage outside the diversity of all known PTLV-3, thus strongly supporting its subtype D designation. Detailed examination of the complete genome predicted that all enzymatic, structural, and regulatory genes were intact. Viral replication and pathogenic potential shown or hypothesized for other PTLV-3s have yet to be determined [14, 15, 30]. Given the inferred ancient origin of STLV-3d(Cmo8699AB), its prevalence in two primate species from the same geographical region, and the documented propensity for STLVs to cross species boundaries, STLV-3d may be more widespread than currently realized. These results underscore an unknown public health concern for STLV-3d, particularly in a region with frequent exposure to NHPs through hunting and butchering.
DNA preparation and PCR-based genome walking
Using the NucliSens nucleic acid isolation kits (Biomérieux, Durham, NC) as previously described , nucleic acids were extracted from two dried blood spots (DBS) each collected by two different hunters from a wild-caught C. mona monkey (Cmo8699AB) and a C. nictitans monkey (Cni7867AB). Due to the limited DBS material available, we successfully maximized DNA yield through additional elution of nucleic acids from the silica beads with water. DNA from Cni3034 and Cni3038 were prepared from whole blood using the Qiagen DNA extraction protocol (Valencia, CA). DNA quality and yield were evaluated in a semi-quantitative PCR amplification of the β-actin gene as previously described [31, 32] and confirmed with the QuantiT dsDNA HS Assay kit (Invitrogen, Carlsbad, CA). A minimum total input of 10 ng of DNA was used in each reaction mixture with standard PCR conditions. DNA preparation and PCR assays were performed in different laboratories specifically equipped for the processing and testing of only NHP samples according to established precautions to prevent contamination.
PCR primer pairs1,2 used to amplify overlapping regions of the STLV-3d(Cmo8699AB) genome
TCA ACC TTT TCT CCC CAA CGC CCT
AYT GGR GGC TRC CWG GGG CGG AAG
GCA CCT TCG CTT CTC CTG TCC TGG
GRT AGG GYG GAG GCT TTT GRG GGT
GTG CCG CCA ACC CCA TCC CCA AGG
GGY RTG IAR CCA RRC IAG KGG CCA
AAA GGG CTA GCA ATT CAC CAC TGG
GAT AGG GTT ATT GCC TGG TCC TTG ATA
ACC CCC CCA GTA AGC ATC CAG GCG
GGY RTG IAR CCA RRC IAG KGG CCA
AGA TGT CCT CCA GCA ATG CCA AAG
GRY RGG IGT ICC TTT IGA GAC CCA
TCC ACA GAA AAA ACC CAA TCC ACT
GGG CAG TAG CAA TGG GAC CAA GGA
CAC TCC TGG TCC CAT ACA CTT TCT CGG
GGT GGG GCC TGT GTA GTT TGG GAG
AAA GTC TAA ACC CTC CAT GCC CAG
TTT GGT AGG GAT TTT TGT TAG GAA GG
TCC TTG TAT CTT TTT CCC CAT TGG
AAG GTA TTG TAG AGG CGA GCT GAC
PCR amplicons were purified with Qiaquick PCR or gel purification kits (QIAGEN, Valencia, CA) and sequenced directly using ABI PRISM Big Dye terminator kits (Foster City, CA) on an ABI 3130xl sequencer or after cloning into a TOPO vector (Invitrogen, Carlsbad, CA).
Sequence and phylogenetic analysis and dating the origin of STLV-3d(Cmo8699AB)
Comparison of the full-length, gap-stripped PTLV-3 genomes was performed with the SimPlot program (Version 3.5.1) where STLV-3d(Cmo8699AB) was the query sequence using the F84 (ML) model and a transition/transversion ratio of 2.0 . RNA secondary structure of the LTR region was predicted using the mfold web server program  found at http://mfold.bioinfo.rpi.edu/. Prediction of splice acceptor (sa) and splice donor (sd) sites was performed using the NetGene2 program available at the web server http://www.cbs.dtu.dk/services/NetGene2/. Identification and analysis of ORFs were performed using the ORF Finder program available at http://www.ncbi.nlm.nih.gov/projects/gorf/.
Percent nucleotide divergence was calculated using the DNASTAR MegAlign 7.2 software (http://www.DNASTAR.com). For phylogenetic analysis two datasets were used. To investigate the phylogenetic relationship between PTLV, the first dataset included tax sequences from complete PTLV genomes available at GenBank and the new STLV-3 tax sequences from Cmo8699AB, Cni7867AB, Cni3034, Cni3038, and Lal9859 obtained in the current study, respectively. For further phylogenetic resolution of STLV-3d among PTLV, a larger dataset was used and included concatenated gag, pol, env, and tax sequences from complete PTLV genomes available at GenBank and the complete genome of STLV-3d(Cmo8699AB) determined here. Sequences were aligned using the Clustal W program, followed by manual editing and removal of indels. Nucleotide substitution saturation was assessed using pair-wise transition and transversion versus divergence plots using the DAMBE program . Unequal nucleotide composition was measured by using the TREE-PUZZLE program . Nucleotide substitution models and parameters were estimated from the edited Clustal W sequence alignments by using Modeltest v3.7 . A variant of the general time reversible (GTR) model, which allows six different substitution rate categories (rA ↔ C = 2.62, rA ↔ G = 13.07, rA ↔ T = 2.79, rC ↔ G = 2.26, rC ↔ T = 4.54, rG ↔ T = 1) with gamma-distributed rate heterogeneity (α = 0.7071) and an estimated proportion of invariable sites (0.3436) was determined to best fit the data for the tax only alignments. The best model for the concatenated gag-pol-env-tax alignment was GTR+G, with six different rate substitutions (rA ↔ C = 2.53, rA ↔ G = 11.47, rA ↔ T = 2.58, rC ↔ G = 2.15, rC ↔ T = 4.3, rG ↔ T = 1) and gamma-distributed rate heterogeneity (α = 0.366). Phylogenetic trees were inferred using Bayesian analysis implemented in the BEAST software package  and with maximum likelihood (ML) using the PhyML program available online at the webserver http://atgc.lirmm.fr/phyml/. Support for branching order of the ML-inferred trees was evaluated using 500 bootstraps. Two independent BEAST runs consisting of 10 - 100 million Markov Chain Monte Carlo (MCMC) generations for the tax only and PTLV concatamer alignments, respectively, with a sampling every 1,000 generations, an uncorrelated log-normal relaxed molecular clock, and a burn-in of 100,000 to 1 million generations. Both the constant coalescent and the Yule process of speciation were used as tree priors to infer the viral tree topologies. Convergence of the MCMC was assessed by calculating the effective sampling size (ESS) of the runs using the program Tracer (v1.4; http://beast.bio.ed.ac.uk/Tracer). All parameter estimates showed significant ESSs (> 300). The tree with the maximum product of the posterior clade probabilities (maximum clade credibility tree) was chosen from the posterior distribution of 9,001 sampled trees (after burning in the first 1,000 sampled trees) with the program TreeAnnotator version 1.4.6 included in the BEAST software package . Trees were viewed and edited using FigTree v1.1.2 http://tree.bio.ed.ac.uk/software/figtree.
Divergence dates for the most recent common ancestor (MRCA) of STLV-3d(Cmo8699AB) were obtained by using both the tax only and the concatenated gag-pol-env-tax alignments, using Bayesian inference and using a relaxed molecular clock in the BEAST program. The PTLV evolutionary rate assumed a global molecular clock model and was estimated according to the formula: evolutionary rate (r) = branch length (bl)/divergence time (t) . Divergence dates were obtained from well-established genetic and archaeological evidence for the timing of migration of the ancestors of indigenous Melanesians and Australians from Southeast Asia [14, 16, 29, 41]. The PTLV evolutionary rate was estimated by using the divergence time of 40,000 - 60,000 years ago (ya) for the Melanesian HTLV-1 lineage (HTLV-1mel) and 15,000-30,000 ya for the most recent common ancestor of HTLV-2a/HTLV-2b native American strains as strong priors in a Bayesian MCMC relaxed molecular clock method implemented in the BEAST software package . The use of two calibration points has previously been shown to provide more reliable estimates of PTLV substitution rates than a single calibration date [41, 42]. The upper and lower divergence times estimated from anthropological data were used to define the interval of a strong uniform prior distribution from which the MCMC sampler would sample possible divergence times for the corresponding node in the tree.
Nucleotide accession numbers
The STLV-3d(Cmo8699AB) complete proviral genome has the GenBank accession number EU231644. Partial STLV-3d genomic sequences obtained from monkey Cni7867AB were assigned the GenBank accession numbers FJ957879 (LTR-partial gag) and FJ957880 (pol-partial env). Longer tax sequences obtained from STLV-3d(Cni7867AB), STLV-3c(Cni3034), STLV-3c(Cni3038), and STLV-3b(Lal9589NL) have the GenBank accession numbers EU152281, FJ957877, FJ957878, and GQ241937, respectively.
Comparison of the STLV-3d(Cmo8699AB) proviral genome with prototypical PTLVs
Percent nucleotide and amino acid identity of STLV-3d(Cmo8699AB) with other prototypical PTLVs1
PTLV-3 (subtype A)
PTLV-3 (subtype B)
Evolutionary relationship of STLV-3d to other PTLVs
Divergence dates for the most recent common ancestor of STLV-3d(Cmo8699AB)
PTLV evolutionary rate and time-scale calculated with a Bayesian relaxed molecular clock using 1st + 2nd codon positions of concatenated gag-pol-env-tax genes and tax only1.
6.29 × 10-7
(3.29 × 10-7 - 9.53 × 10-7)
5.36 × 10-7
(3.21 × 10-7 - 8.1 × 10-7)
(147,042 - 529,980)
(88,914 - 299,436)
(58,833 - 109,552)
(45,899 - 118,645)
(38,355 - 76,651)
(39,783 - 59,155)
(77,653 - 305,591)
(46,324 - 180,712)
(41,349 - 182,273)
(29,660 - 111,773)
(11,650 - 87,100)
(8,744 - 56,742)
(14,419 - 40,104)
(13,591 - 27,792)
HTLV-2a, b 4
(14,426 - 28,212)
(13,703 - 27,783)
(52,894 - 201,260)
(34,342 - 127,209)
(26,648 - 102,445)
(17,149 - 68,097)
(4,386 - 36,666)
(52,822 - 200,926)
Genomic organization of STLV-3d(Cmo8699AB) and identification of conserved functional motifs
STLV-3d proteome analysis
The predicted protein translation of the STLV-3d(Cmo8699AB) genome revealed all major structural and enzymatic (Gag, Pro, Pol, and Env) and regulatory proteins (Tax and Rex) (Fig. 1 and Table 2). Analysis of the overlapping open reading frames (ORFs) of gag and pro and pro and pol predicts that translation occurs by one or more successive -1 ribosomal frameshifts that align different ORFs. The conserved-slippage sequence (6(A)-8 nt-6(G)-11 nt-6(C) can be found in the gag-pro overlap of STLV-3d(Cmo8699AB). The pro-pol overlap slippage sequence has the same point mutation found among the other prototypical PTLV-3s (GTTAAAC versus TTTAAAC in HTLV-1, HTLV-2, and HTLV-4). Comparable to other PTLV-3s, the Gag protein of STLV-3d(Cmo8699AB) is composed of 420 amino acids (aa) and is predicted to cleave into three core protein products: p19 (matrix), p24 (capsid), and p15 (nucleocapsid). One of the most highly conserved PTLV domains, the Gag protein of STLV-3d(Cmo8699AB) has > 88% similarity to that of prototypical PTLV-3 subtypes (Table 2). The highest amino acid similarity to other PTLV-3 subtypes is found in the p24 capsid protein (94 - 96%), while the p15 nucleocapsid protein was the most divergent (80 - 83%).
The predicted length of the STLV-3d(Cmo8699AB) Env glycoprotein is 493 aa, similar to the Env protein of STLV-3b(CtoNG604) and HTLV-3b(Pyl43) (10, 35). The surface (SU) and transmembrane (TM) proteins are comparable to all other PTLV-3 subtypes at 315 aa and 178 aa, respectively. The TM protein is highly conserved across PTLV-3 subtypes (90 - 92% similarity) including STLV-3d(Cmo8699AB). The high aa identity of the Gag p24 and Env proteins suggests that this divergent virus would be cross-reactive on standard HTLV-1/2 Western blot (WB) assays. Unfortunately, serum or plasma was not available from animals Cmo8699AB and Cni7867AB to confirm this hypothesis. The STLV-3d(Cmo8699AB) SU also contains highly conserved residues believed important for viral entry (data not shown) similar to those described recently for HTLV-3b(Pyl43) .
PTLV Tax proteins are important for the trans-activation of viral gene expression, viral replication and viral pathogenesis. Comparison of the Tax proteins of prototypical PTLVs and STLV-3d revealed the conservation of critical functional motifs including the nuclear localization signal (NLS), cAMP response element (CREB) binding protein (CBP)/P300 binding motifs, and nuclear export signal (NES) motifs (data not shown). Amino acid sequences (M1, M22, and M47) that are important for Tax1 transactivation and activation of κβ (NF-κβ) pathway  are also preserved in the STLV-3d(Cmo8699AB) and STLV-3d(Cni7867AB) Tax proteins (data not shown). The C-terminal transcriptional activating domain (CR2) at positions 313 - 318 of the protein is important for CBP/P300 binding and up-regulation of transcription and is also present. The CR2 motif [(S/T)T(V/I)PFS] is conserved among all PTLV-3 subtypes and is identical to those found in STLV-3a subtypes. In the Tax C-terminus, STLV-3d also possesses a conserved PDZ-binding motif present in PTLV-1 and PTLV-3 Tax but not in PTLV-2 or HTLV-4 [14, 30, 42, 49, 50]. The PDZ domain has been shown to be an important binding site for Tax in mediating signal transduction and interleukin-2-independent growth induction for T-cell transformation [50, 51]. Taken together, preservation of the predicted STLV-3d(Cmo8699AB) Tax protein sequence motifs suggests Tax interactions with cellular regulatory pathways similar to those of both PTLV-1 and PTLV-3. All functional motifs, including a potential PDZ domain, are present in the STLV-3d Tax, although 12 aa residues are missing from the N-terminus of the Tax proteins of STLV-3c(Cni3034 and Cni3038) obtained in the current study. This suggests that these motifs are highly conserved among the very divergent PTLV-3 group (data not shown).
While comparable in length to other PTLV-3s, the 182 aa Rex protein of STLV-3d(Cmo8699AB) is the most divergent viral protein sharing only about 70% similarity. The activation domain/NES (DALSARLYNTLSLDSPP) (aa 81 - 97) is an important motif conserved in all PTLVs for shuttling unspliced viral RNA transcripts from the nucleus to the cytoplasm [52, 53]. However, like HTLV-2 and HTLV-4, the Rex protein of STLV-3d(Cmo8699AB) has an alanine (A) at position 94 instead of the glycine (G) found in HTLV-3b(Pyl43) and STLV-3b(CTONG409) or aspartic acid (D) present in all other PTLV-3s.
The pX encoding region between env and the 3' LTR contains multiple coding regions shown to be important for HTLV-1 viral replication T-cell activation, and cellular gene expression with two of the open reading frames (ORFs) encoding for the ubiquitous Tax and Rex proteins . Two putative splice donor sites with high confidence were predicted at positions 414 and 5058 in the LTR (sd-LTR) and Env (sd-Env), respectively, that code for the Env protein (Fig. 1). A conserved splice acceptor site is located at position 7552 that with the sd-Env site code for the singly spliced Tax and Rex proteins (Fig. 1). The positions of these putative splice junction sites are similar to those of other PTLV-3s [14, 15, 21–23, 27]. Analysis of the pX region of STLV-3d(Cmo8699AB) revealed only a single additional ORF (ORFI that begins with a methionine and is predicted to code for a proteins of 131 aa in length (Fig. 1)), in contrast to other PTLV-3s which have been predicted to have at least two additional ORFs in the pX region. BLAST analysis of the ORFI protein resulted in matches to miscellaneous fungal and mammalian proteins with very low identity (< 30%). Further studies are required to evaluate the function of the ORFI viral protein.
Screening of human populations with high exposure to NHPs has resulted in the successful discovery of novel retroviruses, including HTLV-3, HTLV-4, and simian foamy virus (SFV) [1, 7, 8, 32]. We have previously demonstrated that hunter-collected DBS specimens from wild-caught NHPs are not only an effective collection strategy to demonstrate STLV diversity but also allow for monitoring of retroviral cross-species transmission events at the primate-hunter interface . Using these primate DBS specimens, we recently identified novel STLV-3s in wild-caught C. mona and C. nictitans monkeys by analysis of partial gene sequences . To characterize this new virus, we obtained its complete proviral genome using nucleic acids extracted entirely from two DBS, collected by two hunters in the field. To our knowledge, this is the first full-length genome of a simian retrovirus obtained entirely from DBS. The ability to generate a complete viral genome from the equivalent of about 0.25 ml whole blood demonstrates further the utility of this collection strategy for monitoring and characterizing viral diversity.
Robust phylogenetic analysis of both the conserved tax region and gag-pol-env-tax concatenated sequences inferred a novel lineage with high statistical support within the PTLV-3 clade that is highly divergent. The formation of a fourth lineage within the diversity of PTLV-3, containing STLV-3 sequences from two distinct primate species (C. mona and C. nictitans), strongly supports the proposed nomenclature and classification of this new virus as STLV-3 subtype D. The discovery of nearly identical STLV-3d(Cmo8699AB and Cni7867AB) viruses in two different primate species within the same region of Cameroon and the inferred ancient divergence of STLV-3d about 115,000 ya also suggests a higher prevalence and a more widespread distribution for this virus.
Our results also suggest that taxonomic subdivision within the current subtype B strains is warranted. PTLV-3 from Cameroon (STLV-3(Cto604), HTLV-3(Pyl43), and HTLV-3(Lobak18)) are distantly related to other PTLV-3 subtype B strains from West-Central Africa ((HTLV-3(Cam2026ND), STLV-3(NG409), STLV-3(PPAF3), and STLV-3(Lal9589NL) sharing < 90% nucleotide identity. Since both PTLV-3 lineages are presently known as B subtypes, we propose re-naming them as subtypes B1 (from Cameroon only) and B2 (from West-Central Africa). A similar nomenclature strategy has been previously adopted for subtyping HTLV-2 . This re-classification and the categorization of the new STLV-3 subtype D are based upon highly supported phylogenetic division of these subtypes and genetic distances of at least 5% across the genome, as recently proposed for PTLV classification .
PTLVs have an ancient evolutionary history with the ancestral HTLVs being inferred to have first occurred many thousands of years ago following zoonotic transmission from STLV-infected NHPs [14, 24, 29, 42, 61]. This finding contrasts with the relatively recent emergence of the human immunodeficiency virus (HIV) from simian immunodeficiency virus-infected NHPs in the last century [62, 63]. The recent discovery of HTLV-3 and HTLV-4 and novel STLV-1-like viruses among people who hunt and butcher NHPs suggests that these interspecies transmission events are not rare and are most likely contemporaneous [1, 7]. From phylogenetic analysis it has been inferred that STLV-1 may have crossed species boundaries to humans on at least seven separate occasions resulting in the multiple HTLV-1 subtypes . Given the inferred ancient origin of STLV-3d(Cmo8699AB) and PTLV-3, the wide geographic distribution of STLV-3 across Africa, the long history of human exposure to simians in Africa and the lack of screening for HTLV in blood banks in Africa, human infections with STLV-3d-like viruses might be expected to occur there. Thus, although HTLV-3 so far has only been identified in three persons from Cameroon and all three are subtype B viruses, it is tempting to speculate that like HTLV-1 diversity, HTLV-3 diversity will be driven by transmission of each of the four STLV-3 subtypes to humans. More surveillance studies at the NHP-human interface are needed to determine the prevalence, diversity, and epidemiology of STLV-3d(Cmo8699AB) and HTLV-3.
Molecular differences between HTLV-1 and HTLV-2 Tax proteins have been proposed to modulate function, transmissibility, and pathogenesis . We therefore examined the predicted protein sequences of STLV-3d(Cmo8699AB) to determine whether important functional and regulatory motifs were present to infer the replication-competency and pathogenic potential for this divergent viral subtype. All enzymatic, structural, and regulatory proteins were preserved in STLV-3d(Cmo8699AB), including the ubiquitous Tax binding domains CBP/P300, NES, and CR2, which are all important for viral transcription and transformation [64–66]. In addition, the presence of a PDZ-binding motif in the STLV-3d Tax, which has been shown to be critical in signal transduction and T-cell transformation of HTLV-1-infected cells [50, 51], suggests that the STLV-3d(Cmo8699AB) Tax is more similar to the Tax of PTLV-1 and other PTLV-3s, than it is to the Tax of PTLV-2 which lacks a PDZ motif [14, 42]. Furthermore, as has been demonstrated with all PTLVs, STLV-3d also possesses a conserved ASP basic leucine zipper motif in the antisense strand between the env and tax/rex gene regions. ASP has been shown to participate in regulation of viral replication and possibly oncogenesis [50, 51]. Combined, these findings show that the STLV-3d(Cmo8699AB) genome is intact, is likely to be replication competent, and may have a pathogenic potential similar to HTLV-1 which is also predicted for HTLV-3 subtype B; however, further studies are required to validate this hypothesis.
The LTR region of STLV-3d(Cmo8699AB) has two of the three 21-bp repeat Tax-responsive elements (TRE) typically found in the HTLV-1 and HTLV-2 LTRs. The three TREs (distal, central, and proximal) are involved in basal transcription and have been shown to confer Tax1, Tax2, and Tax3 responsiveness . Studies have also shown that mutations in the central TRE compared to the distal or proximal TRE-1 result in the greatest loss of basal transcription levels . As with all PTLV-3s, the STLV-3d(Cmo8699AB) LTR lacks only the distal TRE, which does not appear to have deleterious effects on gene expression and viral replication [30, 69]. Nonetheless, more studies are necessary to determine if these differences will affect the transcriptional activity of STLV-3d(Cmo8699AB).
Another notable difference of STLV-3d from other PTLV-3s was observed in the leucine-rich activation region of the putative NES domain of the Rex protein involved in regulation of viral expression. STLV-3d(Cmo8699AB) has a single aa mutation from aspartic acid or glycine to alanine at position 94 similar to that seen in the HTLV-2 Rex protein (Rex2). Mutagenesis studies substituting alanine for serine residues in this region have demonstrated a significant reduction in the phosphorylation activation required for efficient RNA binding of Rex-2 . These results suggest that the alanine mutation at aa position 94 of the STLV-3d(Cmo8699AB) Rex may also have a similar loss of biologic activity. The effects of these changes on the processing of viral transcripts and regulation of viral replication by the STLV-3d Rex will require further investigation.
In summary, complete genome analysis of STLV-3d(Cmo8699) reveals this novel virus is a highly divergent member of the PTLV-3 group that we name subtype D. We show by robust genetic analysis that STLV-3d(Cmo8699AB) has an ancient origin and an intact genome. Furthermore, we demonstrate that complete viral genomes can be obtained using limited amounts of genomic material extracted from DBS collected in the field. This collection strategy will facilitate the monitoring of viral diversity and cross-species transmission at the human-primate interface. Expanded surveillance will help us to better understand the epidemiology and public health importance of STLV zoonoses.
DMS was funded through a National Science Foundation Graduate Research Fellowship and the Edward and Kathy Ludwig Scholarship. NDW was supported by awards from the National Institutes of Health Director's Pioneer Award (Grant DP1-OD000370), the WW Smith Charitable Trust, the US Military HIV Research Program, and grants from the NIH Fogarty International Center (International Research Scientist Development Award Grant 5 K01 TW000003-05), AIDS International Training and Research Program (Grant 2 D 43 TW000010-17-AITRP), and the National Geographic Society Committee for Research and Exploration (Grant #7762-04). This research was supported in part by the Global Viral Forecasting Initiative. The entire staff of GVFI-Cameroon is thanked for their support and assistance. The collaboration of numerous hunters participating voluntarily in the GVFI surveillance program is also appreciated. The Cameroon Ministry of Defense, Ministry of Scientific Research and Innovation and Ministry of Forestry and Fauna provided authorizations and support for this work. Use of trade names is for identification only and does not imply endorsement by the U.S. Department of Health and Human Services, the Public Health Service, or the Centers for Disease Control and Prevention. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention.
- Calattini S, Chevalier SA, Duprez R, Bassot S, Froment A, Mahieux R, Gessain A: Discovery of a new human T-cell lymphotropic virus (HTLV-3) in Central Africa. Retrovirology. 2005, 2: 30-10.1186/1742-4690-2-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Gessain A, Mahieux R: [Epidemiology, origin and genetic diversity of HTLV-1 retrovirus and STLV-1 simian affiliated retrovirus]. Bull Soc Pathol Exot. 2000, 93: 163-171.PubMedGoogle Scholar
- Mahieux R, Chappey C, Georges-Courbot MC, Dubreuil G, Mauclere P, Georges A, Gessain A: Simian T-cell lymphotropic virus type 1 from Mandrillus sphinx as a simian counterpart of human T-cell lymphotropic virus type 1 subtype D. J Virol. 1998, 72: 10316-10322.PubMed CentralPubMedGoogle Scholar
- Meertens L, Rigoulet J, Mauclere P, Van Beveren M, Chen GM, Diop O, Dubreuil G, Georges-Goubot MC, Berthier JL, Lewis J, Gessain A: Molecular and phylogenetic analyses of 16 novel simian T cell leukemia virus type 1 from Africa: close relationship of STLV-1 from Allenopithecus nigroviridis to HTLV-1 subtype B strains. Virology. 2001, 287: 275-285. 10.1006/viro.2001.1018.View ArticlePubMedGoogle Scholar
- Slattery JP, Franchini G, Gessain A: Genomic evolution, patterns of global dissemination, and interspecies transmission of human and simian T-cell leukemia/lymphotropic viruses. Genome Res. 1999, 9: 525-540.PubMedGoogle Scholar
- Van Brussel M, Salemi M, Liu HF, Goubau P, Desmyter J, Vandamme AM: The discovery of two new divergent STLVs has implications for the evolution and epidemiology of HTLVs. Rev Med Virol. 1999, 9: 155-170. 10.1002/(SICI)1099-1654(199907/09)9:3<155::AID-RMV242>3.0.CO;2-3.View ArticlePubMedGoogle Scholar
- Wolfe ND, Heneine W, Carr JK, Garcia AD, Shanmugam V, Tamoufe U, Torimiro JN, Prosser AT, Lebreton M, Mpoudi-Ngole E, et al: Emergence of unique primate T-lymphotropic viruses among central African bushmeat hunters. Proc Natl Acad Sci USA. 2005, 102: 7994-7999. 10.1073/pnas.0501734102.PubMed CentralView ArticlePubMedGoogle Scholar
- Calattini S, Betsem E, Froment A, Bassot S, Chevalier SA, Mahieux R, Gessain A: New strain of human T lymphotropic virus (HTLV) type 3 in a Pygmy from Cameroon with peculiar HTLV serologic results. J Infect Dis. 2009, 199: 561-564. 10.1086/596206.View ArticlePubMedGoogle Scholar
- Van Dooren S, Meertens L, Lemey P, Gessain A, Vandamme AM: Full-genome analysis of a highly divergent simian T-cell lymphotropic virus type 1 strain in Macaca arctoides. J Gen Virol. 2005, 86: 1953-1959. 10.1099/vir.0.80520-0.View ArticlePubMedGoogle Scholar
- Araujo A, Hall WW: Human T-lymphotropic virus type II and neurological disease. Ann Neurol. 2004, 56: 10-19. 10.1002/ana.20126.View ArticlePubMedGoogle Scholar
- Proietti FA, Carneiro-Proietti AB, Catalan-Soares BC, Murphy EL: Global epidemiology of HTLV-I infection and associated diseases. Oncogene. 2005, 24: 6058-6068. 10.1038/sj.onc.1208968.View ArticlePubMedGoogle Scholar
- Roucoux DF, Murphy EL: The epidemiology and disease outcomes of human T-lymphotropic virus type II. AIDS Rev. 2004, 6: 144-154.PubMedGoogle Scholar
- Yamashita M, Ido E, Miura T, Hayami M: Molecular epidemiology of HTLV-I in the world. J Acquir Immune Defic Syndr Hum Retrovirol. 1996, 13 (Suppl 1): S124-131. 10.1097/00042560-199600001-00021.View ArticlePubMedGoogle Scholar
- Switzer WM, Qari SH, Wolfe ND, Burke DS, Folks TM, Heneine W: Ancient origin and molecular features of the novel human T-lymphotropic virus type 3 revealed by complete genome analysis. J Virol. 2006, 80: 7427-7438. 10.1128/JVI.00690-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Calattini S, Chevalier SA, Duprez R, Afonso P, Froment A, Gessain A, Mahieux R: Human T-cell lymphotropic virus type 3: complete nucleotide sequence and characterization of the human tax3 protein. J Virol. 2006, 80: 9876-9888. 10.1128/JVI.00799-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Song KJ, Nerurkar VR, Saitou N, Lazo A, Blakeslee JR, Miyoshi I, Yanagihara R: Genetic analysis and molecular phylogeny of simian T-cell lymphotropic virus type I: evidence for independent virus evolution in Asia and Africa. Virology. 1994, 199: 56-66. 10.1006/viro.1994.1097.View ArticlePubMedGoogle Scholar
- Van Dooren S, Verschoor EJ, Fagrouch Z, Vandamme AM: Phylogeny of primate T lymphotropic virus type 1 (PTLV-1) including various new Asian and African non-human primate strains. Infect Genet Evol. 2007, 7: 374-381. 10.1016/j.meegid.2006.06.003.View ArticlePubMedGoogle Scholar
- Courgnaud V, Van Dooren S, Liegeois F, Pourrut X, Abela B, Loul S, Mpoudi-Ngole E, Vandamme A, Delaporte E, Peeters M: Simian T-cell leukemia virus (STLV) infection in wild primate populations in Cameroon: evidence for dual STLV type 1 and type 3 infection in agile mangabeys (Cercocebus agilis). J Virol. 2004, 78: 4700-4709. 10.1128/JVI.78.9.4700-4709.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Goubau P, Van Brussel M, Vandamme AM, Liu HF, Desmyter J: A primate T-lymphotropic virus, PTLV-L, different from human T-lymphotropic viruses types I and II, in a wild-caught baboon (Papio hamadryas). Proc Natl Acad Sci USA. 1994, 91: 2848-2852. 10.1073/pnas.91.7.2848.PubMed CentralView ArticlePubMedGoogle Scholar
- Liegeois F, Lafay B, Switzer WM, Locatelli S, Mpoudi-Ngole E, Loul S, Heneine W, Delaporte E, Peeters M: Identification and molecular characterization of new STLV-1 and STLV-3 strains in wild-caught nonhuman primates in Cameroon. Virology. 2008, 371: 405-17. 10.1016/j.virol.2007.09.037.View ArticlePubMedGoogle Scholar
- Meertens L, Gessain A: Divergent simian T-cell lymphotropic virus type 3 (STLV-3) in wild-caught Papio hamadryas papio from Senegal: widespread distribution of STLV-3 in Africa. J Virol. 2003, 77: 782-789. 10.1128/JVI.77.1.782-789.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Meertens L, Mahieux R, Mauclere P, Lewis J, Gessain A: Complete sequence of a novel highly divergent simian T-cell lymphotropic virus from wild-caught red-capped mangabeys (Cercocebus torquatus) from Cameroon: a new primate T-lymphotropic virus type 3 subtype. J Virol. 2002, 76: 259-268. 10.1128/JVI.76.1.259-268.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Meertens L, Shanmugam V, Gessain A, Beer BE, Tooze Z, Heneine W, Switzer WM: A novel, divergent simian T-cell lymphotropic virus type 3 in a wild-caught red-capped mangabey (Cercocebus torquatus torquatus) from Nigeria. J Gen Virol. 2003, 84: 2723-2727. 10.1099/vir.0.19253-0.View ArticlePubMedGoogle Scholar
- Sintasath DM, Wolfe ND, Lebreton M, Jia H, Garcia AD, Le Doux-Diffo J, Tamoufe U, Carr JK, Folks TM, Mpoudi-Ngole E, et al: Simian T-lymphotropic virus diversity among nonhuman primates, Cameroon. Emerg Infect Dis. 2009, 15: 175-184. 10.3201/eid1502.080584.PubMed CentralView ArticlePubMedGoogle Scholar
- Takemura T, Yamashita M, Shimada MK, Ohkura S, Shotake T, Ikeda M, Miura T, Hayami M: High prevalence of simian T-lymphotropic virus type L in wild Ethiopian baboons. J Virol. 2002, 76: 1642-1648. 10.1128/JVI.76.4.1642-1648.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Dooren S, Salemi M, Pourrut X, Peeters M, Delaporte E, Van Ranst M, Vandamme AM: Evidence for a second simian T-cell lymphotropic virus type 3 in Cercopithecus nictitans from Cameroon. J Virol. 2001, 75: 11939-11941. 10.1128/JVI.75.23.11939-11941.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Dooren S, Shanmugam V, Bhullar V, Parekh B, Vandamme AM, Heneine W, Switzer WM: Identification in gelada baboons (Theropithecus gelada) of a distinct simian T-cell lymphotropic virus type 3 with a broad range of Western blot reactivity. J Gen Virol. 2004, 85: 507-519. 10.1099/vir.0.19630-0.View ArticlePubMedGoogle Scholar
- Koralnik IJ, Boeri E, Saxinger WC, Monico AL, Fullen J, Gessain A, Guo HG, Gallo RC, Markham P, Kalyanaraman V, et al: Phylogenetic associations of human and simian T-cell leukemia/lymphotropic virus type I strains: evidence for interspecies transmission. J Virol. 1994, 68: 2693-2707.PubMed CentralPubMedGoogle Scholar
- Salemi M, Desmyter J, Vandamme AM: Tempo and mode of human and simian T-lymphotropic virus (HTLV/STLV) evolution revealed by analyses of full-genome sequences. Mol Biol Evol. 2000, 17: 374-386.View ArticlePubMedGoogle Scholar
- Chevalier SA, Meertens L, Pise-Masison C, Calattini S, Park H, Alhaj AA, Zhou M, Gessain A, Kashanchi F, Brady JN, Mahieux R: The tax protein from the primate T-cell lymphotropic virus type 3 is expressed in vivo and is functionally related to HTLV-1 Tax rather than HTLV-2 Tax. Oncogene. 2006, 25: 4470-4482. 10.1038/sj.onc.1209472.View ArticlePubMedGoogle Scholar
- Busch MP, Switzer WM, Murphy EL, Thomson R, Heneine W: Absence of evidence of infection with divergent primate T-lymphotropic viruses in United States blood donors who have seroindeterminate HTLV test results. Transfusion. 2000, 40: 443-449. 10.1046/j.1537-2995.2000.40040443.x.View ArticlePubMedGoogle Scholar
- Wolfe ND, Switzer WM, Carr JK, Bhullar VB, Shanmugam V, Tamoufe U, Prosser AT, Torimiro JN, Wright A, Mpoudi-Ngole E, et al: Naturally acquired simian retrovirus infections in central African hunters. Lancet. 2004, 363: 932-937. 10.1016/S0140-6736(04)15787-5.View ArticlePubMedGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999, 73: 152-160.PubMed CentralPubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMed CentralView ArticlePubMedGoogle Scholar
- Hebsgaard SM, Korning PG, Tolstrup N, Engelbrecht J, Rouze P, Brunak S: Splice site prediction in Arabidopsis thaliana pre-mRNA by combining local and global sequence information. Nucleic Acids Res. 1996, 24: 3439-3452. 10.1093/nar/24.17.3439.PubMed CentralView ArticlePubMedGoogle Scholar
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online--a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33: W557-559. 10.1093/nar/gki352.PubMed CentralView ArticlePubMedGoogle Scholar
- Lemey P, Pybus OG, Van Dooren S, Vandamme AM: A Bayesian statistical analysis of human T-cell lymphotropic virus evolutionary rates. Infect Genet Evol. 2005, 5: 291-298. 10.1016/j.meegid.2004.04.005.View ArticlePubMedGoogle Scholar
- Switzer WM, Salemi M, Qari SH, Jia H, Gray RR, Katzourakis A, Marriott SJ, Pryor KN, Wolfe ND, Burke DS, et al: Ancient, independent evolution and distinct molecular features of the novel human T-lymphotropic virus type 4. Retrovirology. 2009, 6: 9-10.1186/1742-4690-6-9.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki T, Fujisawa JI, Toita M, Yoshida M: The trans-activator tax of human T-cell leukemia virus type 1 (HTLV-1) interacts with cAMP-responsive element (CRE) binding and CRE modulator proteins that bind to the 21-base-pair enhancer of HTLV-1. Proc Natl Acad Sci USA. 1993, 90: 610-614. 10.1073/pnas.90.2.610.PubMed CentralView ArticlePubMedGoogle Scholar
- Yin MJ, Paulssen EJ, Seeler JS, Gaynor RB: Protein domains involved in both in vivo and in vitro interactions between human T-cell leukemia virus type I tax and CREB. J Virol. 1995, 69: 3420-3432.PubMed CentralPubMedGoogle Scholar
- Kress E, Baydoun HH, Bex F, Gazzolo L, Duc Dodon M: Critical role of hnRNP A1 in HTLV-1 replication in human transformed T lymphocytes. Retrovirology. 2005, 2: 8-10.1186/1742-4690-2-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Lieu YK, Kumar A, Pajerowski AG, Rogers TJ, Reddy EP: Requirement of c-myb in T cell development and in mature T cell function. Proceedings of the National Academy of Sciences. 2004, 101: 14853-14858. 10.1073/pnas.0405338101.View ArticleGoogle Scholar
- Jones KS, Huang YK, Chevalier SA, Afonso PV, Petrow-Sadowski C, Bertolette DC, Gessain A, Ruscetti FW, Mahieux R: The receptor complex associated with human T-cell lymphotropic virus type 3 (HTLV-3) Env-mediated binding and entry is distinct from, but overlaps with, the receptor complexes of HTLV-1 and HTLV-2. J Virol. 2009, 83: 5244-5255. 10.1128/JVI.02285-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith MR, Greene WC: Identification of HTLV-I tax trans-activator mutants exhibiting novel transcriptional phenotypes. Genes Dev. 1990, 4: 1875-1885. 10.1101/gad.4.11.1875.View ArticlePubMedGoogle Scholar
- Hirata A, Higuchi M, Niinuma A, Ohashi M, Fukushi M, Oie M, Akiyama T, Tanaka Y, Gejyo F, Fujii M: PDZ domain-binding motif of human T-cell leukemia virus type 1 Tax oncoprotein augments the transforming activity in a rat fibroblast cell line. Virology. 2004, 318: 327-336. 10.1016/j.virol.2003.10.006.View ArticlePubMedGoogle Scholar
- Rousset R, Fabre S, Desbois C, Bantignies F, Jalinot P: The C-terminus of the HTLV-1 Tax oncoprotein mediates interaction with the PDZ domain of cellular proteins. Oncogene. 1998, 16: 643-654. 10.1038/sj.onc.1201567.View ArticlePubMedGoogle Scholar
- Tsubata C, Higuchi M, Takahashi M, Oie M, Tanaka Y, Gejyo F, Fujii M: PDZ domain-binding motif of human T-cell leukemia virus type 1 Tax oncoprotein is essential for the interleukin 2 independent growth induction of a T-cell line. Retrovirology. 2005, 2: 46-10.1186/1742-4690-2-46.PubMed CentralView ArticlePubMedGoogle Scholar
- Ballaun C, Farrington GK, Dobrovnik M, Rusche J, Hauber J, Bohnlein E: Functional analysis of human T-cell leukemia virus type I rex-response element: direct RNA binding of Rex protein correlates with in vivo activity. J Virol. 1991, 65: 4408-4413.PubMed CentralPubMedGoogle Scholar
- Kusuhara K, Anderson M, Pettiford SM, Green PL: Human T-cell leukemia virus type 2 Rex protein increases stability and promotes nuclear to cytoplasmic transport of gag/pol and env RNAs. J Virol. 1999, 73: 8112-8119.PubMed CentralPubMedGoogle Scholar
- Bindhu M, Nair A, Lairmore MD: Role of accessory proteins of HTLV-1 in viral replication, T cell activation, and cellular gene expression. Front Biosci. 2004, 9: 2556-2576. 10.2741/1417.View ArticlePubMedGoogle Scholar
- Gaudray G, Gachon F, Basbous J, Biard-Piechaczyk M, Devaux C, Mesnard JM: The complementary strand of the human T-cell leukemia virus type 1 RNA genome encodes a bZIP transcription factor that down-regulates viral transcription. J Virol. 2002, 76: 12813-12822. 10.1128/JVI.76.24.12813-12822.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Arnold J, Yamamoto B, Li M, Phipps AJ, Younis I, Lairmore MD, Green PL: Enhancement of infectivity and persistence in vivo by HBZ, a natural antisense coded protein of HTLV-1. Blood. 2006, 107: 3976-3982. 10.1182/blood-2005-11-4551.PubMed CentralView ArticlePubMedGoogle Scholar
- Landry S, Halin M, Lefort S, Audet B, Vaquero C, Mesnard JM, Barbeau B: Detection, characterization and regulation of antisense transcripts in HIV-1. Retrovirology. 2007, 4: 71-10.1186/1742-4690-4-71.PubMed CentralView ArticlePubMedGoogle Scholar
- Basbous J, Arpin C, Gaudray G, Piechaczyk M, Devaux C, Mesnard JM: The HBZ factor of human T-cell leukemia virus type I dimerizes with transcription factors JunB and c-Jun and modulates their transcriptional activity. J Biol Chem. 2003, 278: 43620-43627. 10.1074/jbc.M307275200.View ArticlePubMedGoogle Scholar
- Shaulian E, Karin M: AP-1 as a regulator of cell life and death. Nat Cell Biol. 2002, 4: E131-136. 10.1038/ncb0502-e131.View ArticlePubMedGoogle Scholar
- Switzer WM, Pieniazek D, Swanson P, Samdal HH, Soriano V, Khabbaz RF, Kaplan JE, Lal RB, Heneine W: Phylogenetic relationship and geographic distribution of multiple human T-cell lymphotropic virus type II subtypes. J Virol. 1995, 69: 621-632.PubMed CentralPubMedGoogle Scholar
- Feuer G, Green PL: Comparative biology of human T-cell lymphotropic virus type 1 (HTLV-1) and HTLV-2. Oncogene. 2005, 24: 5996-6004. 10.1038/sj.onc.1208971.PubMed CentralView ArticlePubMedGoogle Scholar
- Hahn BH, Shaw GM, De Cock KM, Sharp PM: AIDS as a zoonosis: scientific and public health implications. Science. 2000, 287: 607-614. 10.1126/science.287.5453.607.View ArticlePubMedGoogle Scholar
- Sharp PM, Bailes E, Gao F, Beer BE, Hirsch VM, Hahn BH: Origins and evolution of AIDS viruses: estimating the time-scale. Biochem Soc Trans. 2000, 28: 275-282.View ArticlePubMedGoogle Scholar
- Alefantis T, Barmak K, Harhaj EW, Grant C, Wigdahl B: Characterization of a nuclear export signal within the human T cell leukemia virus type I transactivator protein Tax. J Biol Chem. 2003, 278: 21814-21822. 10.1074/jbc.M211576200.View ArticlePubMedGoogle Scholar
- Kashanchi F, Brady JN: Transcriptional and post-transcriptional gene regulation of HTLV-1. Oncogene. 2005, 24: 5938-5951. 10.1038/sj.onc.1208973.View ArticlePubMedGoogle Scholar
- Scoggin KE, Ulloa A, Nyborg JK: The oncoprotein Tax binds the SRC-1-interacting domain of CBP/p300 to mediate transcriptional activation. Mol Cell Biol. 2001, 21: 5520-5530. 10.1128/MCB.21.16.5520-5530.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Shimotohno K, Takano M, Teruuchi T, Miwa M: Requirement of multiple copies of a 21-nucleotide sequence in the U3 regions of human T-cell leukemia virus type I and type II long terminal repeats for trans-acting activation of transcription. Proc Natl Acad Sci USA. 1986, 83: 8112-8116. 10.1073/pnas.83.21.8112.PubMed CentralView ArticlePubMedGoogle Scholar
- Barnhart MK, Connor LM, Marriott SJ: Function of the human T-cell leukemia virus type 1 21-base-pair repeats in basal transcription. J Virol. 1997, 71: 337-344.PubMed CentralPubMedGoogle Scholar
- Van Brussel M, Goubau P, Rousseau R, Desmyter J, Vandamme AM: Complete nucleotide sequence of the new simian T-lymphotropic virus, STLV-PH969 from a Hamadryas baboon, and unusual features of its long terminal repeat. J Virol. 1997, 71: 5464-5472.PubMed CentralPubMedGoogle Scholar
- Narayan M, Younis I, D'Agostino DM, Green PL: Functional domain structure of human T-cell leukemia virus type 2 rex. J Virol. 2003, 77: 12829-12840. 10.1128/JVI.77.23.12829-12840.2003.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.