Tax Gene Characterization of Human T-Lymphotropic Virus Type 1 Strains from Brazilian HIV-Coinfected Patients

Abstract

The tax gene of human T-lymphotropic virus type 1 (HTLV-1) diverges among isolates according to geographic regions and has been classified into two genotypes: taxA and taxB. In Brazil, taxA is the most prevalent genotype in symptomatic and asymptomatic carriers. Few studies have been conducted in HIV-infected patients. The present study characterized the tax gene (1059 bp) in 13 Brazilian HIV-1/HTLV-1-coinfected patients from the south and southeast regions. The results confirmed the transcontinental HTLV-1 subgroup A of the Cosmopolitan subtype and showed high nucleotide similarity both among Brazilian sequences and in relation to the ATK prototype (99.5% and 99.2%, respectively). Six nucleotide substitutions were highly conserved among isolates, ranging from 76.9% to 100%: C7401T, T7914C, C7920T, C7982T, G8231A, and A8367C. The presence of the Brazilian molecular signature of genotype taxA was confirmed in all of the isolates, and they clustered into two Latin American clusters, which confirms the double introduction of HTLV-1 in Brazil.

H uman T-lymphotropic virus type 1 (HTLV-1) is the etiological agent of the hematological disorder adult T cell leukemia/lymphoma (ATLL) and the neurologic syndrome HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP).^1,2 It is estimated that 20 million individuals are infected by HTLV-1 worldwide, mainly in southwestern Japan, the Caribbean area, some regions of Africa, and Central and South America.^1,2 Based on the phylogenetic analysis of the long terminal repeat (LTR) and env regions of the genome of HTLV-1, different subtypes and subgroups of HTLV-1 have been described around the world, mostly related to human migratory routes and geographic areas.^1,2

HTLV-1 is widespread in the northeast region of Brazil,^3,4 and studies on its origin and dissemination showed that HTLV-1 was recently introduced in this region, most likely through the post-Columbian migrations of the African slave trade between the sixteenth and nineteenth centuries.^5,6 HTLV-1 has also been detected in other regions of the country, but with a smaller percentage of cases than in the northeast.⁴ The molecular characterization of HTLV-1 in Brazil showed that the Cosmopolitan subtype (a) of transcontinental subgroup A was prevalent, although historical data indicated that the slave trade movement originated from west Africa, where only the western African subgroup (C) circulates.^7,8 One study on the genetic background of HTLV-1-infected individuals from Brazil confirmed the introduction of HTLV-1 by infected individuals from the southern region of Africa.⁷

Brazil also has the largest Japanese population outside of Japan, almost 1.5 million, and 1 million live in São Paulo (southeastern Brazil).⁹ Therefore, the introduction of the Japanese HTLV-1 subgroup B in this geographic area could be hypothesized. In fact, this subgroup was detected in some individuals from southeastern Brazil.¹⁰

On the basis of the molecular analysis of the tax gene, HTLV-1 strains were classified into two genotypes, taxA and taxB, whose association with clinical disease was suggested.¹¹

Although there are few studies on tax molecular characterization in Brazil, one study with a large number of individuals (symptomatic and asymptomatic carriers) from São Paulo and Rio de Janeiro (southeast) characterized the entire tax gene and demonstrated a Brazilian molecular signature (five nucleotide substitutions in relation to the ATK prototype: C7401T, T7914C, C7920T, C7982T, and G8231A) characteristic of the taxA genotype.¹⁰ Corroborating these data, another study from São Paulo showed the same pattern of mutations and a lack of association between the tax genotype and disease progression.¹²

In addition to some regions of Brazil being considered endemic for HTLV-1 infection,⁴ Brazil is considered epidemic for HIV/AIDS, with prevalence rates of 0.3% to 0.6%¹³; as a result, HIV/HTLV-1 coinfection has been described in this country, with prevalence rates varying according to the geographic area. For example, in the cities of Salvador (northeastern Brazil), Santos (southeastern Brazil), and Londrina (southern Brazil), HTLV-1 prevalence rates of 17.1%, 6.0%, and 0.8%, respectively, were detected in HIV-coinfected individuals.^14
–16

To gain further insight into the history of HTLV-1 in Brazil, especially in the at-risk population infected by HIV-1 in the south and southeast regions, we conducted the present and preliminary studies. A few years ago, we determined the prevalence of HTLV-1 and HTLV-2 in patients with HIV/AIDS from Londrina and its vicinities in the south of Brazil and detected HIV/HTLV-1/2 coinfection in 6.4% of cases¹⁶; subsequently, we characterized HTLV-1 isolates from southern and southeastern Brazil using the LTR and env regions¹⁷ and confirmed the HTLV-1a Cosmopolitan subtype, subgroup A, which clustered into two Latin American clusters.

In the present study, another genomic region (tax) was analyzed for the molecular signature described in Brazilian isolates. The study population consisted of 13 patients with HIV/AIDS who were coinfected with HTLV-1 from Londrina and its vicinities (BRLO) and São Paulo (BRSP).¹⁷ Five patients are from Londrina (five men) and eight are from São Paulo (five women and three men). DNA samples were extracted from peripheral blood leukocytes, and the tax region was amplified using a previously described modified protocol,¹¹ in which PCR was optimized using the GoTaq Colorless Master Mix (Promega Corporation, Madison, WI) to improve stability. Fragments of 1119 bp were obtained, and 1059 bp were sequenced with an ABI 3130 Genetic Analyzer (Applied Biosystems, Foster, CA). All of the sequencing chromatograms were assembled and 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 (http://www.ncbi.nlm.nih.gov/genbank), in which the nucleotide and amino acid substitutions were searched.

HTLV-1 subtyping was screened with the NCBI Genotyping (http://www.ncbi.nlm.nih.gov/projects/genotyping/formpage.cgi) tool website. Neighbor-joining (NJ) and maximum likelihood (ML) phylogenetic trees were constructed based on an appropriate nucleotide substitution models determined by Modeltest v3.7 (TrN+G model) using PAUP v4b10 software. Bootstrapping was performed with the stepwise addition algorithm for 1000 replicates. The EL (HTLV-1b) sequence was used as the outgroup. Nucleotide distances were estimated using MEGA4 software. Sequences of 1059 bp (nt 7324–8382 in relation to the ATK prototype) were initially analyzed using the NCBI Genotyping tool, and the presence of the HTLV-1a subtype was found in all of the study samples. The phylogenetic analysis of a 977-bp segment of 49 HTLV-1 tax sequences available in the GenBank database, including the 13 new sequences, was used to construct the tax phylogenetic tree (Fig. 1). The phylogenetic tree confirmed that all of the study sequences belong to the HTLV-1a subtype, genotype taxA. Moreover, the 13 new tax sequences clustered into two different clusters inside the genotype taxA.

FIG. 1.

Phylogenetic tree constructed by the neighbor-joining method using the PAUP v4b10 software for the partial tax region of 977 bp (position 7378–8354 relative to the ATK prototype) of human T-lymphotropic virus type 1 (HTLV-1) of 49 isolates, including sequences from the south and southeast regions of Brazil (GenBank AN JN887698–JN887710) in bold, generated with the TrN+G model. Bootstrap values above 65% and zero length using the likelihood ratio test with p<0.001 (**) and p≤0.05 (*) in key branches are depicted. The HTLV-1 EL isolate was used as the outgroup.

Of note, the previous molecular analysis of the LTR and env regions of such isolates characterized these strains as belonging to the HTLV-1a Cosmopolitan subtype, transcontinental subgroup A, clustering into two branches of Latin America,¹⁷ the same branches detected in the present study (Fig. 1). These findings strengthen the suggested concept of the different introduction of HTLV-1 in Brazil^{3,5

–8,10,17} and among HIV-coinfected patients.

Nucleotide similarity estimations showed a similarity of 99.5% among all of the study sequences and when they were stratified into two groups: one from São Paulo and the other from Londrina and its vicinities. The analysis of sequences in relation to the ATK prototype showed a nucleotide similarity of 99.2%. Interestingly, the comparative analysis of HTLV-1 isolates from Chile with the ATK prototype showed 99.6% nucleotide similarity in the tax region.¹⁸ Because Brazil and Chile are located in Latin America, these results could suggest an earlier introduction of HTLV-1 in Brazil.

Six nucleotide substitutions were highly conserved among the 13 new Brazilian sequences in relation to ATK: C7401T and T7914C in 76.9% of the sequences (except for Latin American cluster B) and C7920T, C7982T, G8231A, and A8367C in 100% of the sequences. These data are presented in Table 1 and confirm the molecular signature previously described in Brazilian sequences of HTLV-1.^10,12 Two of the nucleotide substitutions at positions C7982T and G8231A resulted in amino acid changes of A2661V and S2744N.

Table 1.

Comparison of tax Region Sequences of Human T-Lymphotropic Virus Type 1, Showing the Nucleotide Substitutions, Amino Acid Changes, and Their Positions in the 13 Samples from the South and Southeast Regions of Brazil, in Relation to the International ATK Prototype

	Nucleotide change at position
Sequence	7401	7410	7431	7476	7527	7593	7722	7908	7914	7920	7953	7974	7980	7982	8001	8002	8066	8109	8128	8184	8191	8231	8232	8288	8295	8313	8366	8367
ATK	C	C	G	T	G	C	A	T	T	C	C	C	T	C	A	A	C	A	G	C	C	G	T	A	G	T	G	A
BRLO14-02	T								C	T		T		T								A				C		C
BRLO15-02					A					T				T	G							A			A			C
BRLO34-02										T				T	G							A			A			C
BRLO37-02	T						G		C	T				T				G				A						C
BRLO48-02	T						G		C	T				T				G	A			A						C
BRSP134-08	T		A	C					C	T		T		T								A				C		C
BRSP145-08	T								C	T	T			T								A		G				C
BRSP206-08	T								C	T				T							T	A	C					C
BRSP232-08	T	A							C	T			C	T								A					A	C
BRSP42-09	T								C	T		T		T								A				C		C
BRSP205-09	T								C	T				T		G						A						C
BRSP320-09						T		C		T				T	G		G					A			A			C
BRSP414-09	T								C	T				T						T		A						C
aa change														A→V		T→A	T→R		V→I			S→N		E→G			R→H

aa, amino acid.

It is important to note that the protocols selected and optimized for the PCR, nested PCR, and sequencing in the present study amplified a large segment of the tax gene, which only lacks the three initial tax nucleotides. These three nucleotides characterize the start codon of Tax (amino acid methionine) and are located at the beginning of the viral envelope gene.

In conclusion, despite the continental size of Brazil, the different social habits and risk factors of its population, and the regional differences in the genetic background of individuals, the same molecular signature of genotype taxA was detected in HTLV-1 strains isolated from HIV-1-coinfected patients from the south and southeast regions. The hypothesis of multiple introductions of HTLV-1 in Brazil was also confirmed in the present study.

Sequence Data

The GenBank accession numbers for the HTLV-1 tax sequences included in the phylogenetic study are as follows: ATK (J02029), TSP-1 (M86840), BRRJ017_95 (DQ323833), BRRJl69_00 (DQ323845), BRRJc127_97 (DQ323840), BRRJ3476 (DQ323834), BRRJa107_96 (DQ323835), BRRP388 (DQ323855), BRRP165 (DQ323875), BRRJa156_971 (DQ323836), BRRP004 (DQ323870), BRRP438 (DQ323862), BRRP050 (DQ323850), BRRP464 (DQ323883), ArJ13-1(DQ227134), ArJ2-1 (DQ227130), ArJ16-1 (DQ227135), ArJ11-5 (DQ227158), ArJ55-2 (DQ227189), BOI (L36905), WHP (AF259264), htlvtax5 (AF226593), ATL150 (AB036376), RKI3-Ger (AF042071), RKI2-Rum (AF003012), HC34-13 (AB045549), MT-2 (AF292000), HC139 (AB045520), ATL-YS (U19949), HAM50 (AB045481), (AF485380), (AF485381), 31/98 (AF133525), 122/97 (AF133524), EL (M67514).

The GenBank accession numbers of the HTLV-1 tax sequences obtained in our laboratory and included in the phylogenetic analysis are as follows: BRLO14-02 (JN887698), BRLO15-02 (JN887699), BRLO34-02 (JN887700), BRLO37-02 (JN887701), BRLO48-02 (JN887702), BRSP134-08 (JN887703), BRSP145-08 (JN887704), BRSP206-08 (JN887705), BRSP232-08 (JN887706), BRSP42-09 (JN887707), BRSP205-09 (JN887708), BRSP320-09 (JN887709), BRSP414-09 (JN887710).

Footnotes

Acknowledgments

This study was supported by the 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). The authors thank Karoline Rodrigues Campos for her technical support.

Author Disclosure Statement

No competing financial interests exist.

References

Gonçalves

, Proietti

, Ribas

et al. Epidemiology, treatment, and prevention of human T-cell leukemia virus type 1-associated diseases. Clin Microbiol Rev, 2010; 23:577–589.

Proietti

, Carneiro-Proietti

ABF

, Catalan-Soares

, Murphy

. Global epidemiology of HTLV-1 infection, associated diseases. Oncogene, 2005; 24:6058–6068.

Dourado

, Alcantara

, Barreto

, Teixeira

, Galvão-Castro

. HTLV-I in the general population of Salvador, Brazil: A city with African ethnic and sociodemographic characteristics. J Acquir Immume Defic Syndr, 2003; 34:527–531.

Catalan-Soares

, Carneiro-Proietti

ABF

, Proietti

. the Interdisciplinary HTLV Research Group: Heterogeneous geographic distribution of human T-cell lymphotropic viruses I and II (HTLV-I/II): Serological screening prevalence rates in blood donors from large urban areas in Brazil. Cad Saúde Púb (RJ), 2005; 21:926–931.

Rego

FFA

, Alcantara

LCJ

, Neto

JPM

, Mota-Miranda

ACA

, Pereira

, Gonçalves

, Galvão-Castro

. HTLV type 1 molecular study in Brazilian villages with African characteristics giving support to the post-Columbian introduction hypothesis. AIDS Res Hum Retroviruses, 2008; 24,5:673–677.

Magalhães

, Mota-Miranda

, Alcantara

LCJ

, Olavarria

, Galvão-Castro

, Rios-Grassi

. Phylogenetic and molecular analysis of HTLV-1 isolates from a medium sized town in northern of Brazil: Tracing a common origin of the virus from the most endemic city in the country. J Med Virol, 2008; 80:2040–2045.

Alcantara

LCJ

, Oliveira

, Gordon

et al. Tracing the origin of Brazilian HTLV-1 as determined by analysis of host, viral genes. AIDS, 2006; 20,5:780–782.

Mota

ACA

, Van Dooren

, Fernandes

FMC

et al. The close relationship between South African and Latin American HTLV type 1 strains corroborated in a molecular epidemiological study of the HTLV type 1 isolates from a blood donor cohort. AIDS Res Hum Retroviruses, 2007; 23,4:503–507.

Governo do Estado de São Paulo. 1908–2008 100 Anos. Imigração Japão Brasil. http://www.saopaulo.sp.gov.br/imigracaojaponesa/historia.php. 2012 March 16.

10.

Kashima

, Alcantara

, Takayanagui

et al. Distribution of human T cell lymphotropic virus type 1 (HTLV-1) subtypes in Brazil: Genetic characterization of LTR and tax region. AIDS Res Hum Retroviruses, 2006; 22,10:953–959.

11.

Furukawa

, Yamashita

, Usuku

, Izumo

, Nakagawa

, Osame

. Phylogenetic subgroups of human T cell lymphotropic virus (HTLV) Type I in the tax gene and their association with different risks for HTLV-I-associated myelopathy/tropical spastic paraparesis. J Infect Dis, 2000; 182:1343–1349.

12.

Casseb

, Fukumori

LMI

, Vergara

MPP

et al. Lack of tax diversity for tropical spastic paraparesis/human T-cell lymphotropic virus type 1 (HTLV-I) associated myelopathy development in HTLV-I-infected subjects in São Paulo, Brazil. Mem Inst Oswaldo Cruz, 2006; 101,3:273–276.

13.

UNAIDS. Global Report–UNAIDS report for the global AIDS epidemic–2010. http://www.unaids.org/globalreport/HIV_prevalence_map.htm. 2012 January 19.

14.

Brites

, Goyanna

, França

et al. Coinfection by HTLV-I/II is associated with an increased risk of strongyloidiasis and delay in starting antiretroviral therapy for AIDS patients. Braz J Infect Dis, 2011; 15,1:6–11.

15.

Etzel

, Shibata

, Rozman

, Jorge

MLSG

, Damas

, Segurado

AAC

. HTLV-1 and HTLV-2 infections in HIV-infected individuals from Santos, Brazil: Seroprevalence and risk factors. J Acquir Immune Defic Syndr, 2001; 26:185–190.

16.

Morimoto

, Caterino-de-Araujo

, Morimoto

et al. Seroprevalence and risk factors for human T-cell lymphotropic virus type 1 and 2 infection in human immunodeficiency virus (HIV)-infected patients attending AIDS Referral Center Health Units in Londrina and other communities in Paraná, Brazil. AIDS Res Hum Retroviruses, 2005; 21,4:256–262.

17.

Magri

, Brigido

LFM

, Rodrigues

, Morimoto

, Ferreira

JLP

, Caterino-de-Araujo

. Phylogenetic and similarity analysis of HTLV-1 isolates from HIV-coinfected patients from the south and southeast regions of Brazil. AIDS Res Hum Retrovir, 2012; 28,1:110–114.

18.

Ramirez

, Cartier

, Villota

, Fernandez

. Genetic characterization and phylogeny of human T-cell lymphotropic virus type I from Chile. Virus Res, 2002; 84,1–2:135–149.