Abstract
HTLV-1 is endemic in Brazil and HIV/HTLV-1 coinfection has been detected, mostly in the northeast region. Cosmopolitan HTLV-1a is the main subtype that circulates in Brazil. This study characterized 17 HTLV-1 isolates from HIV coinfected patients of southern (n=7) and southeastern (n=10) Brazil. HTLV-1 provirus DNA was amplified by nested PCR (env and LTR) and sequenced. Env sequences (705 bp) from 15 isolates and LTR sequences (731 bp) from 17 isolates showed 99.5% and 98.8% similarity among sequences, respectively. Comparing these sequences with ATK (HTLV-1a) and Mel5 (HTLV-1c) prototypes, similarities of 99% and 97.4%, respectively, for env and LTR with ATK, and 91.6% and 90.3% with Mel5, were detected. Phylogenetic analysis showed that all sequences belonged to the transcontinental subgroup A of the Cosmopolitan subtype, clustering in two Latin American clusters.
T
The HTLV-1 retrovirus has a conserved genome, with only 0.1% to 6.9% genetic diversity among strains. The env and long terminal repeat (LTR) regions have the highest variability, 4 and are useful for viral subtype characterization. Mostly individuals remain asymptomatic, but this infection may cause a wide spectrum of diseases, including leukemia (adult T cell leukemia/lymphoma, ATL) and neurodegenerative disorders (HTLV-1-associated myelopathy/tropical spastic paraparesis, HAM/TSP). 1 In Brazil, the majority of HTLV-1 molecular characterization studies were conducted in patients presenting ATL and HAM/TSP and in asymptomatic blood donors. Studies of HTLV-1 in HIV-coinfected patients are limited and are attributed to HTLV-1a, the main strain also circulating in HIV/HTLV-1-coinfected individuals. 3,5 –8
The present study was conducted to add information concerning HTLV-1 strains that circulate in HIV-coinfected patients from the south and southeast regions of Brazil, thought to have more cases of HIV/HTLV-2 coinfection. 9,10 Seventeen HIV/HTLV-1-coinfected patients attending Aids Reference Centers, seven (one woman and six men) from Londrina-LO (south), and 10 (seven women and three men) from São Paulo-SP (southeast), were enrolled in this analysis. HTLV genomic DNA was extracted from peripheral blood leukocytes (PBLs). Nested polymerase chain reaction (n-PCR) was performed to amplify a 763-bp env and 766-bp LTR region of HTLV-1, using protocols previously described (Tables 1 and 2) 6,11 in which some primers have undergone minor nucleotide modifications and thermal cycling conditions were adjusted. Env and LTR amplified products were sequenced in an ABI 3100 Genetic Analyser (Applied Biosystems, USA), with the same primers used in the n-PCR.
Primer nucleotide position is provided as aligned with ATK (HTLV-1-infected cell line; GenBank Accession number J02029). bp, base pairs; n-PCR, nested polymerase chain reaction.
Adapted from Caterino-de-Araujo et al. 11
Adapted from Laurentino et al. 6
Primers were employed to detect HTLV-1 provirus DNA in peripheral blood leucocytes of HIV coinfected patients from the south (Londrina-LO) and the southeast (São Paulo-SP) regions of Brazil.
Adapted from Caterino-de-Araujo et al. 11
Adapted from Laurentino et al. 6
Protocol for detecting HTLV-1 provirus DNA in peripheral blood leukocytes of HIV co-infected patients from the south (Londrina-LO) and the southeast (São Paulo-SP) regions of Brazil.
All sequencing chromatograms obtained were assembled and manually edited with Sequencher 4.7 software. Multiple alignments were performed using the Clustal W multiple-sequence alignment tool from BioEdit Sequence Alignment Editor version 7.0.5.3 software, with a reference set available in the GenBank database (
Env sequences (705 bp) from 15 isolates and LTR sequences (731 bp) from 17 isolates were obtained. Molecular analysis of env and LTR sequences disclosed, respectively, nucleotide similarities of 99.5% and 98.8% among sequences, 99% and 97.4% with ATK, and 91.6% and 90.3% with Mel5.
Phylogenetic analysis of a segment of 705 bp of 44 HTLV-1 env sequences available in the GenBank database including the 15 new sequences from the south and southeast regions of Brazil were used to construct the env phylogenetic tree (Fig. 1). A segment of 692 bp of 59 HTLV-1 LTR sequences available in GenBank including the 17 new sequences of this study was used to construct the LTR phylogenetic tree (Fig. 2).
Phylogenetic tree constructed by the maximum-likelihood method using the PAUP v4b10 software for the partial env region of 705 bp (position 5614–6318 relative to the ATK prototype) of HTLV-1 of 44 isolates, including sequences from the south and southeast regions of Brazil (GenBank AN HM770426-HM770440) in bold, generated with the TrN + G model. Bootstrap values above 70% and zero length using the likelihood ratio test with p<0.001 (**) and p≤0.05 (*) in key branches are depicted. The HTLV-1 Mel5 isolate was used as the outgroup.
Phylogenetic tree constructed by the maximum-likelihood method using the PAUP v4b10 software for the partial LTR region of 692 bp (position 75–766 relative to the ATK prototype) of HTLV-1 of 59 isolates, including sequences from the south and southeast regions of Brazil (GenBank AN JF271836-JF271852) in bold, generated with the GTR + G model. Bootstrap values above 70% and zero length using the likelihood ratio test with p<0.001 (**) and p≤0.05 (*) in key branches are depicted. The HTLV-1 Mel5 isolate was used as the outgroup.
Four sequences that showed the highest similarities with ATK clustered separately in the trees (p≤0.001 for env and LTR, bootstrap=96% for LTR); they had C386T, T594C and C595T nucleotide substitutions in the LTR region (BRLO14-02, BRLO30-02, BRLO34-02, and SP320-09). All except these four sequences had mutations T5637G, T5730C, and A6120G in env and T500A, A549G, G676A, T766C, and T767G in the LTR regions. Six sequences presented the amino acid change V1981I and clustered together in the env tree (BRLO48-02, BRSP145-08, BRSP205-09, BRSP232-08, BRSP414-09, and BRSP01-10).
According to the NCBI-Genotyping and the REGA-Subtyping websites all isolates from the present study belonged to the Cosmopolitan HTLV-1a subtype, the same subtype detected by phylogenetic analysis (bootstrap value of 80% for env and 87% for LTR, and a p<0.001 for both). In addition, all isolates belonged to the Transcontinental subgroup A (supported by a p<0.001) and clustered in two Latin American clusters (A and B) in the LTR tree (supported by a p<0.001 in ML analyses).
These phylogenetic results are in accordance with the results obtained in HTLV-1 strains characterized from the north and northeast regions of Brazil. In fact, nucleotide similarities in LTR sequences of 99.2% and 98.9%, respectively, were detected in strains isolated from individuals from Cruz das Almas and Salvador cities (southeastern Brazil), and of 99.3% and 97.3%, respectively, in strains isolated from individuals from Fortaleza and Rio Branco cities (northern Brazil). 12,13
Molecular analysis of HTLV-1 isolates from HIV-1/HTLV-1-coinfected individuals from northern Brazil disclosed the Cosmopolitan HTLV-1a subtype circulating in four women from Feira de Santana, state of Bahia, 6 and in two men from Belém, state of Para. 7 Corroborating these and the present study, phylogenetic analysis of HTLV-1 isolates from several populations in Brazil confirms HTLV-1a, Transcontinental subgroup A as the most prevalent in the country, 8,12,13 and suggests a post-Columbian introduction of HTLV-1 into the Brazilian population in the state of Bahia and maybe in other places in Latin America. 3,7,12,14
Indeed, the majority of the sequences from this study clustered into Latin American Cluster A, although some sequences clustered in Latin American Cluster B. The well-divided clusters of Latin American could suggest different introductions of HTLV-1 in Brazil, and also in HIV/HTLV-1-coinfected individuals (both from the south and southeast regions). Still, the high level of similarities among sequences leads us to speculate on a recent introduction of HTLV-1 in Brazil.
In conclusion, these data show that the Cosmopolitan HTLV-1a subtype is also frequent in HIV coinfected patients from southern and southeastern Brazil, suggesting spreading of HTLV-1 in the country. The mutations observed should be monitored in the context of molecular epidemiology. Moreover, this is the first study of the molecular characterization of HTLV-1 sequences from Londrina, southern Brazil, providing insight into the HTLV-1 dynamic epidemic in this country.
Sequence Data
The GenBank accession numbers for the HTLV-1 sequences included in the phylogenetic study are as follows:
The GenBank accession numbers of the HTLV-1 fragments sequenced in our laboratory and included in the phylogenetic analysis are as follows:
Footnotes
Acknowledgments
This study was supported by Ministério da Ciência e Tecnologia/Conselho Nacional de Desenvolvimento Científico e Tecnológico (MCT/CNPq), Brazil (Universal Grant 481040/2007-2), a fellowship to A.C.A. (Grant 303328/2009-6), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Brazil (Ph.D. fellowship to M.C.M.), and Instituto Adolfo Lutz (Grant 39/07).
Author Disclosure Statement
No competing financial interests exist.