Skip to main content

HTLV-1-infected CD4+ T-cells display alternative exon usages that culminate in adult T-cell leukemia



Reprogramming cellular gene transcription sustains HTLV-1 viral persistence that ultimately leads to the development of adult T-cell leukemia/lymphoma (ATLL). We hypothesized that besides these quantitative transcriptional effects, HTLV-1 qualitatively modifies the pattern of cellular gene expression.


Exon expression analysis shows that patients’ untransformed and malignant HTLV-1+ CD4+ T-cells exhibit multiple alternate exon usage (AEU) events. These affect either transcriptionally modified or unmodified genes, culminate in ATLL, and unveil new functional pathways involved in cancer and cell cycle. Unsupervised hierarchical clustering of array data permitted to isolate exon expression patterns of 3977 exons that discriminate uninfected, infected, and transformed CD4+ T-cells. Furthermore, untransformed infected CD4+ clones and ATLL samples shared 486 exon modifications distributed in 320 genes, thereby indicating a role of AEUs in HTLV-1 leukemogenesis. Exposing cells to splicing modulators revealed that Sudemycin E reduces cell viability of HTLV-1 transformed cells without affecting primary control CD4+ cells and HTLV-1 negative cell lines, suggesting that the huge excess of AEU might provide news targets for treating ATLL.


Taken together, these data reveal that HTLV-1 significantly modifies the structure of cellular transcripts and unmask new putative leukemogenic pathways and possible therapeutic targets.


HTLV-1 causes chronic infection that relies on the persistent clonal expansion of infected CD4+ and CD8+ T-cells. Reprogramming of gene transcription sustains HTLV-1 viral persistence and ultimately leads to the development of adult T-cell leukemia/lymphoma (ATLL) from a CD4+-infected clone in a minority of carriers after a prolonged latency [1]. Transcriptome analyses using microarray technology has provided key insights into the biological processes involved in the pre-malignant and malignant expansion of infected cells [2]-[14]. HTLV-1 has been found to deregulate the expression of numerous genes involved in key cellular pathways such as the cell cycle, apoptosis, telomeres/telomerase, DNA repair, and immune and inflammatory responses through expression of the Tax oncoprotein [2]-[6].

Recent works have highlighted that in addition to their quantitative effects on gene expression, numerous pathogenic processes, such as persistent viral infections [15] or tumor development [16],[17], rely on acquired alternate exon usage (AEU) events. As for other retroviruses, alternative splicing plays a pivotal role in HTLV-1 expression. From its 5’LTR, HTLV-1 transcribes a single polycistronic pre-mRNA that codes for structural and enzymatic proteins required for viral particle production. This pre-mRNA also undergoes multiple alternative splicing events that generate mono-spliced transcripts coding for the regulatory proteins p21Rex, p12 and p13, and double-spliced transcripts coding for Tax, p27Rex and p30 [18],[19]. Similarly, minus-strand transcription initiated from the 3’LTR generates spliced and unspliced RNA isoforms of HBZ that synthesize HBZ proteins with distinct properties on cell proliferation [20]-[22]. How HTLV-1 intervenes in alternative splicing processes is still incompletely understood. It has been reported that the RNA-binding protein Rex regulates viral splicing through interacting with host splicing machinery to inhibit viral RNA splicing and export unspliced and single-spliced transcripts [23]-[25]. Gene-by-gene analyses have further shown that HTLV-1 may affect alternative splicing of cellular genes including CD44 [26] and IL-6- and IL-2-receptors [27],[28], proposing that HTLV-1-induced AEUs might contribute to molecular mechanisms that underlie the clonal expansion and the malignant transformation of infected cells. However, no systematic study has been hitherto conducted to ascertain the extent of alternative splicing modifications upon HTLV-1 infection. Here, by using integrative analysis of exon expression profiles and gene ontology of CD4+ T-cells derived from infected individuals with and without malignancy, we show that HTLV-1 induces multiple AEU alterations that unmask new putative leukemogenic pathways and possible therapeutic targets.

Results and discussion

Comparative microarray analysis of exon expression profiles was performed with three ATLL samples and 12 untransformed CD4+ T-cell clones (six infected) derived from HTLV-1-infected individuals with no clinical signs of malignancy. ATLL samples were obtained from patients with an acute form of ATLL (>95% circulating malignant cells). CD4+ clones were obtained through cloning by limiting dilution of peripheral blood mononuclear cells (PBMCs) derived from three HTLV-1-infected individuals with tropical spastic paraparesis/HTLV-1-associated myelopathy with a disease duration of 6, 11, and >26 years. This study was conducted according to the principles outlined in the Declaration of Helsinki, and approved by the Institutional Review Board of the Hospices Civils de Lyon (France). As previously described [29],[30], CD4+ clones were submitted to ex vivo culture for one month prior to HTLV-1 screening and RNA extraction. At this time, infected and uninfected cells remained non-immortalized and required IL-2 for continued growth. Given that in vitro T-cell activation is known to modify AEU [31],[32], exon array analysis of untransformed CD4+ clones was carried out with RNAs extracted from unstimulated and phytohemagglutinin (PHA)-stimulated CD4+ T-cells (Additional file 1: Table S1). By this approach that take into account the in vitro cell culture, we assumed that significant changes in exon expression between infected and uninfected clones mainly resulted from the global impact of HTLV-1 infection, irrespective of activation status. The microarray data have been deposited in toto into the Gene Expression Omnibus database and are available under record number GSE52244.

HTLV-1 expression was assessed by quantitative RT-PCR (qRT-PCR) of Tax transcripts. In accordance with previous reports, Tax transcripts were weakly detected in 2 out of 3 ATLL samples (Additional file 2: Figure S1) [33]-[35]. In comparison, there was a wide range of Tax expression between clones (0.16-2.20; median, 1.2; mean ± standard deviation [SD], 1.172 ± 0.799). From these data, we assumed that HTLV-1 samples exhibited typical features of untransformed- and transformed-infected CD4+ T-cells in vivo. Computational analyses of 15 exon arrays and annotation of exon events were achieved as detailed in Figure 1. Figure 1A represents the distribution of quantitative and qualitative gene modifications in untransformed HTLV-1-infected CD4+ clones and ATLL cells. A total of 18/20 (90%) array-predicted exon usages were validated by exon-specific RT-PCR (Figures 1C and D). Among those, we found alternative splicing changes in IL2-R and CD44 transcripts (Figure 1C). These data confirmed and extended previous observations indicating that IL2-R and CD44 variant isoforms are detected in the PBMC of infected individuals [26],[28]. Array-predicted changes in transcript abundance were validated by qRT-PCR for 8/9 (88%) genes (Additional file 3: Figure S2).

Figure 1
figure 1

Distribution of AEU in ATLL cells and cloned CD4 +T-cells derived from HTLV-1-infected individuals. (A) Distribution of quantitative and qualitative HTLV-1 modifications in untransformed HTLV-1-positive CD4 +clones and ATLL samples. The Venn-diagrams show the distribution of genes modified at the whole gene expression level (white), AEU (dark grey), or both (light grey) in ATLL samples (bottom) and untransformed infected clones (top; six independent data sets including three PHA-activated and three unactivated clones) versus uninfected clones (six independent data sets including three PHA-activated and three unactivated clones). Graphs present the distribution of each class of AEU annotated according to the FASTERDB database [41]. Acceptor: alternative acceptor site; deletion: exon deletion; donor: alternative donor site; exon: exon skipping; intron-ret: intron retention; polya: alternative polyadenylation site; prom: alternative promoter. The total number of each class of exon event is indicated on the histogram bar. (B) GO analysis of genes (n = 329) presenting common AEUs in untransformed-infected clones and in ATLL samples (compared with uninfected clones). The complete set of genes featured in microarrays (42304 Ensembl geneIDs) was used as a reference background. Gene annotations are presented in Additional file 5: Table S3. (C) RT-PCR validation of microarray-predicted exon events in CD4 +T-cells derived from HTLV-1 carriers compared with uninfected T-cell clones. Exon-specific RT-PCR was performed with pooled RNA samples derived from six infected (I) and six uninfected (UI) clones, thereby reflecting the distribution of splice variants in all clones irrespective of activation status. Numbers indicate the expected band (bp) of PCR products, and the SI and p-value are indicated for each exon event; a SI ≥ 1.2 was considered to be a significant change in exon expression and used for comparisons. (D) RT-PCR validation of microarray-predicted exon events in ATLL cells (array) and exon-specific RT-PCR of six additional ATLL samples compared to uninfected T-cell clones.

To examine the global impact of HTLV-1 on AEU of CD4+ T-cells, we investigated AEU in untransformed CD4+ T-cells by comparing the pattern of exon expression in infected versus uninfected clones (Figure 1A; Additional file 1: Table S1). Overall, 2201 AEU events were identified in 1271 genes, of which only 656 (51%) were found altered at the level of whole gene expression. This demonstrated that in untransformed cells, HTLV-1 significantly modifies the exon content of multiple transcripts, considering half of these genes were transcriptionally unmodified. Secondly, when ATLL samples were compared with uninfected clones, the overall number of AEU events was found to be 11-fold (23504 versus 2201) higher than that distinguishing HTLV-1-positive from -negative untransformed CD4+ T-cells (Figure 1A). In addition to this quantitative difference, the proportion of AEU events found in differentially expressed genes was significantly higher in ATLL samples (72% versus 51% in uninfected clones; p < 0.0001, Pearson’s Chi-squared test; Figure 1A). Together these results indicate that the distribution of AEU alterations was significantly different between ATLL and HTLV-1-positive untransformed cells.

The comparative analysis of AEU between the two categories of infected cells revealed that infected CD4+ clones and ATLL samples shared 486 AEUs in 329 genes (p = 1.83 × 10−21, Hypergeometric test), suggesting that 22% of AEUs arising at the chronic stage of infection might pertain to ATLL development. The number of AEUs per gene ranged from one to 20 (median, 1; mean ± SD, 1.47 ± 1.42; Additional file 4: Table S2). GO enrichment analysis (DAVID resources, performed with this set of genes identified several enriched terms that fit well with HTLV-1 infection and related diseases, such as pathways in cancer (fold enrichment [FE]: 2.71), endocytosis (FE: 4.45), antigen processing and presentation (FE: 7.23), cell adhesion molecules (CAMs; FE: 4.56), and natural killer cell-mediated cytotoxicity (FE: 4.46; Figure 1B; Additional file 5: Table S3).

Hierarchical clustering of exon array data showed that alternative expression of 3977 exons (2218 genes) permitted sample gathering in three distinct clusters, including untransformed-uninfected and infected CD4+ clones, as well as ATLL samples (p < 0.05, Kruskal-Wallis ANOVA; Figure 2A). For untransformed cells, these exon expression profiles were irrespective of PHA stimulation, implying a direct role of HTLV-1 in these recurrent alternative splicing modifications. ATLL cells exhibited marked changes in AEU compared to untransformed HTLV-1-positive or -negative CD4+ clones. Such distribution suggests that some AEU events could represent an ATLL molecular signature. Therefore, we used exon-specific RT-PCR to analyze six additional samples harvested from patients with acute forms of ATLL. Over the 11 events identified in Figure 1D, the expression of 4 AEUs was found to be recurrent in ATLL (Figure 1D). Whether the AEUs shown in Figures 1B, D, and 2A are involved in ATLL development and possess diagnostic and/or prognostic implications deserves further investigation.

Figure 2
figure 2

Exon-based hierarchical clustering. Hierarchical clustering analysis was performed with Mev4.0 software ( using the gene-normalized exon intensities. (A) The selection of exons useful for differentiating between the 15 cell samples was statistically analyzed using a Kruskal-Wallis ANOVA with p < 0.05 representing statistical significance. (B) GO analysis. Top pathways are presented. The gene set corresponded to 2000 genes that displayed the highest splicing index (SI) values. The complete set of genes featured in microarrays (42304 Ensembl geneIDs) was used as a reference background. Official gene symbol annotations are presented in Additional file 5: Table S3.

A GO analysis was performed to gain insight into the functional significance of the exon expression profiles (3977 exons) that discriminate between the three groups of samples presented in Figure 2A. Significant enrichment was detected for genes involved in important biological processes, such as those related to cancer (50 genes, FE = 1.9; p = 6.6 × 10−6), MAPK signaling (46 genes, FE: 2.2; p = 6.1 × 10−7), focal adhesion (34 genes, FE: 2.2; p = 5.5 × 10−5), cytokine-cytokine receptor interactions (34 genes, FE: 1.7; p = 0.003), and endocytosis (32 genes, FE: 2.2; p = 3.2 × 10−5; Figure 2B). Official symbol annotation of genes with AEUs that are enriched in these pathways are presented in (Additional file 6: Table S4). Secondly, GO analysis of AEU events detected in ATLL samples (7616 genes) and untransformed HTLV-1-positive CD4+ T-cell clones (1271 genes) showed distinct biological changes (Additional file 7: Figure S3). For example, genes encoding spliceosome components appeared exclusively enriched in ATLL cells, while cell cycle-related genes were found enriched in untransformed and transformed HTLV-1 positive samples. Moreover, AEU ontological enrichments were different from those of differentially expressed genes (Additional file 7: Figure S3). Together, these results suggest that HTLV-1-related transcriptional and post-transcriptional changes might pose distinct functional impacts on untransformed and malignant cells.

The plethora of splicing events found in transformed and untransformed-infected cells suggests that spliceosome components might constitute new therapeutic targets in the management of HTLV-1 infection. Thus, we examined whether the spliceosome modulators Sudemycin E and D1 [36] exerted anti-proliferative activity against HTLV-1 T-cell lines. To this end, cell viability upon drug exposure was compared between HTLV-1-transformed cell lines (MT2 and HUT-102) and uninfected cells, including primary CD4+ T-cells and the HTLV-1-negative acute lymphoblastic leukemia cell line, MOLT4. Figure 3A shows that increasing the concentration of Sudemycin D1 led to a rapid decrease in the viability of human cell lines whether they carried HTLV-1 proviruses or not, while growth of primary CD4+ T-cells was inhibited at higher drug concentrations. In contrast, a similar Sudemycin E concentration range was found to significantly decrease the viability of HTLV-1 cell lines, while displaying weak effects on primary CD4+ and MOLT4 cells. High concentration of Sudemycin E (>5 μM) affected cell viability of all tested cells, carrying or not HTLV-1 (not shown). These data suggest that, compared to both transformed and untransformed-uninfected cells, HTLV-1 tumoral cells are highly sensitive to drugs that target splicing machinery. Time-course analysis confirmed these results and further revealed that the anti-proliferative effect of Sudemycin E was persistent in MT2 and HUT-102 HTLV-1-positive cells relative to HTLV-1-negative CD4+ and MOLT4 cells (Figure 3B).

Figure 3
figure 3

Anti-proliferative effect of Sudemycin E on HTLV-1-transformed cell lines. (A) Cell viability was determined by MTT assay (Cell Titer 96® Non-Radioactive Cell Proliferation Assay, Promega) 72 h after Sudemycin (D1 or E) exposure. Data represent normalized mean ± SD; *p < 0.05 was considered statistically significant using a Mann–Whitney U-test. (B) Time-course analysis over 96 h using 1 μM of Sudemycin. Experiments were performed twice, in triplicate for each cell type, with DMSO used as control. Data represent the normalized mean ± SD. (C) Sudemycin E-induced expression of short spliceoforms encoding caspases 2 and 9 and MDM2. Exon-specific RT-PCR was carried out 72 h after Sudemycin E exposure. The oligonucleotides used have been previously described [36]. (D) Sudemycin E reduces Tax expression. Western blot analysis of Tax expression was carried out with a mouse monoclonal anti-Tax antibody (clone 474) 72 h after Sudemycin E exposure. Histograms represent the quantification of Tax signals normalized against beta-actin. (E) Sudemycin E decreases gene expression level of HBZ and the relative ratio sHBZ/unHBZ. Isoform-specific qRT-PCR analysis was carried out after 72 h of Sudemycin exposure. HBZ isoform expressions were normalized to hprt-1 transcripts used as internal control.

Further analyses examined the impact of Sudemycin E on the expression of viral and cellular genes. Figure 3C shows that Sudemycin E dose-dependently induced the expression of short spliceoforms of MDM2 and caspases 2 and 9 that have been previously correlated with apoptosis induction [37],[38]. In contrast, the expression of ubiquitin, resulting from constitutive splicing of intron 3, remained unaffected. This indicates that Sudemycin E only affected alternative splicing, rather than whole gene transcription, of apoptotic regulatory genes. In parallel, as Tax is synthetized by double-spliced viral transcripts, we examined whether Sudemycin E could affect Tax protein levels. As shown in Figure 3D, western blot analysis of MT2 and HuT102 cell lines revealed that increasing the concentration of Sudemycin E led to reduced Tax expression. Because ATLL cells barely express Tax but regularly show high levels of HBZ mRNA, we next aimed at testing the effects of Sudemycin E on HBZ expression. In contrast to Tax, HBZ expression is not strictly dependent on splicing events. In fact, HTLV-1 minus strand encodes spliced (sHBZ) and unspliced isoforms (unHBZ) of HBZ, yet sHBZ is the most abundant in HTLV-1 transformed cell lines and in fresh ATLL samples [20],[39]. In contrast to its unspliced counterpart, sHBZ promotes cell proliferation in vitro thereby likely explaining the positive correlation between sHBZ expression and HTLV-1 proviral loads in vivo [21],[39]. These data prompted us to assess the expression levels of both sHBZ and usHBZ RNA isoforms in Sudemycin E exposed cells using exon-specific quantitative RT-PCR analysis. The results illustrated in the Figure 3E show that both MT2 and HUT102 cells exhibited reduced ratios of sHBZ/unHBZ upon Sudemycin E together with a dramatic decrease in the overall level of unHBZ expression, which is consistent with the transactivating functions of sHBZ on HBZ gene expression [40]. Given the crucial role of both Tax and sHBZ in conferring proliferative and anti-apoptotic properties, these data help explain the heightened sensitivity of HTLV-1-infected cells to the anti-proliferative properties of Sudemycin E compared to their uninfected counterparts. Taken together, our results indicate that by impairing splicing and, in turn, gene expression of HTLV-1 oncogenes, splicing modulators such as Sudemycin E might constitute promising tools for ATLL treatment.


Our investigation revealed that HTLV-1 significantly qualitatively modifies gene transcripts, as well as unmasks new putative leukemogenic pathways and possible therapeutic targets. Because alternative splicing can result in mRNA isoforms coding for proteins with different biological activity, our results suggest that excessive splicing alterations in HTLV-1-positive cells might underlie their phenotypic plasticity. Future studies are needed to address the molecular mechanisms underlying HTLV-1-induced AEU and their involvement in clonal persistence, immune escape, and/or cellular transformation of infected CD4+ cells, as well as their putative role in treatment resistance and relapse of ATLL.


T cell limiting cloning

Peripheral blood mononuclear cells (PBMCs) were obtained by Ficoll separation of whole blood of HTLV-1 infected individuals. For T-cell limiting dilution cloning, PBMCs were seeded at 0.1 cell/well in Terasaki plates after removal of adherent cells. T-lymphocytes were cultured in RPMI 1640 containing penicillin and streptomycin, sodium pyruvate, non-essential amino acids, 2-mercaptoethanol, 10% filtered human AB serum, 100 U/mL recombinant IL-2 (Chiron Corporation), and 75 μM HTLV-1 integrase inhibitor L-731,988 [30]. T-lymphocytes were re-stimulated every 14 days with PHA (1 μg/mL) and fresh feeder cells (5.105 cells/mL) that were composed of lethally irradiated PBMCs from three distinct allogenic HTLV-1-negative donors. Clones were phenotyped by flow cytometry using antibodies against CD4 and CD8 (DakoCytomation) and isotype-matched controls on a FACScan system using CellQuest software (Becton Dickinson). HTLV-1-positive clones were assessed by PCR using long-terminal repeat region–specific primers as previously described [29]. For each clone, RNA extraction was performed before and 24 h after PHA stimulation.

Microarray analysis

Affymetrix exon array hybridization was completed by labeling 1 μg of TRIzol-purified total RNA with Affymetrix reagents according to the manufacturer’s instructions. Hybridization cocktails containing 5–5.5 μg cDNA were prepared and hybridized to Affymetrix-GeneChip Human Exon 1.0 ST arrays (Affymetrix). Affymetrix Expression Console Software was used for quality assessment, while exon array data were normalized using quintile normalization and analyzed using FasterDB annotation ( [41]. Background correction and probe selection were performed as described previously [42]. Statistical analyses were performed using a Student’s t-test on the splicing index (SI) that corresponds to comparison of gene-normalized exon intensity values between two experimental conditions [30],[42]. According to exon-specific RT-PCR (Figures 1C and D), microarray data were considered statistically significant for p < 0.05 and SI ≥ 1.2.

Polymerase chain reaction

PCR was performed using the Herculase II Fusion DNA Polymerase (Agilent Technologies) in 30–35 cycles. PCR products were analyzed by agarose gel electrophoresis with ethidium bromide under UV light. Primers were designed to encompass alternative splicing events. Primer sequences are available upon request.

Quantitative RT-PCR

Expression of Tax mRNA was quantified with Rotor Gene (Qiagen) as previously described [29],[43]. Expression of sHBZ and unHBZ was quantified using primers previously described by Murata et al. [44]. qRT-PCR was performed in triplicate and relative quantitation (RQ) was calculated by the 2ddCT method to normalize gene expression to the endogenous control U6 (ENSG00000206625) or hprt1 (ENSG00000165704).

Drug exposure and cell viability

Cell viability was determined by MTT assay (Cell Titer 96® Non-Radioactive Cell Proliferation Assay, Promega) 72 h after Sudemycin (D1 or E) exposure. The MTT dye was added during the last hour of incubation. After, supernatants were removed and cells were incubated with 100 μL/well Solubilization Solution/Stop Mix for 1 h at room temperature. The OD was measured using an ELISA reader at 570 nm, with 650 nm as a reference. Experiments were completed twice, in triplicate for each cell type, with DMSO used as control.

Additional files


  1. Mortreux F, Gabet AS, Wattel E: Molecular and cellular aspects of HTLV-1 associated leukemogenesis in vivo. Leukemia. 2003, 17: 26-38. 10.1038/sj.leu.2402777.

    Article  CAS  PubMed  Google Scholar 

  2. Ng PW, Iha H, Iwanaga Y, Bittner M, Chen Y, Jiang Y, Gooden G, Trent JM, Meltzer P, Jeang KT, Zeichner SL: Genome-wide expression changes induced by HTLV-1 Tax: evidence for MLK-3 mixed lineage kinase involvement in Tax-mediated NF-kappaB activation. Oncogene. 2001, 20: 4484-4496. 10.1038/sj.onc.1204513.

    Article  CAS  PubMed  Google Scholar 

  3. Pise-Masison CA, Radonovich M, Mahieux R, Chatterjee P, Whiteford C, Duvall J, Guillerm C, Gessain A, Brady JN: Transcription profile of cells infected with human T-cell leukemia virus type I compared with activated lymphocytes. Cancer Res. 2002, 62: 3562-3571.

    CAS  PubMed  Google Scholar 

  4. Taylor JM, Ghorbel S, Nicot C: Genome wide analysis of human genes transcriptionally and post-transcriptionally regulated by the HTLV-I protein p30. BMC Genomics. 2009, 10: 311-10.1186/1471-2164-10-311.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Tsukasaki K, Tanosaki S, DeVos S, Hofmann WK, Wachsman W, Gombart AF, Krebs J, Jauch A, Bartram CR, Nagai K, Tomonaga M, Said JW, Koeffler HP: Identifying progression-associated genes in adult T-cell leukemia/lymphoma by using oligonucleotide microarrays. Int J Cancer. 2004, 109: 875-881. 10.1002/ijc.20028.

    Article  CAS  PubMed  Google Scholar 

  6. Zane L, Sibon D, Legras C, Lachuer J, Wierinckx A, Mehlen P, Delfau-Larue MH, Gessain A, Gout O, Pinatel C, Lancon A, Mortreux F, Wattel E: Clonal expansion of HTLV-1 positive CD8+ cells relies on cIAP-2 but not on c-FLIP expression. Virology. 2010, 407: 341-351. 10.1016/j.virol.2010.07.023.

    Article  CAS  PubMed  Google Scholar 

  7. Chevalier SA, Durand S, Dasgupta A, Radonovich M, Cimarelli A, Brady JN, Mahieux R, Pise-Masison CA: The transcription profile of Tax-3 is more similar to Tax-1 than Tax-2: insights into HTLV-3 potential leukemogenic properties. PLoS One. 2012, 7: e41003-10.1371/journal.pone.0041003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. de La Fuente C, Deng L, Santiago F, Arce L, Wang L, Kashanchi F: Gene expression array of HTLV type 1-infected T cells: Up-regulation of transcription factors and cell cycle genes. AIDS Res Hum Retroviruses. 2000, 16: 1695-1700. 10.1089/08892220050193164.

    Article  CAS  PubMed  Google Scholar 

  9. Harhaj EW, Good L, Xiao G, Sun SC: Gene expression profiles in HTLV-I-immortalized T cells: deregulated expression of genes involved in apoptosis regulation. Oncogene. 1999, 18: 1341-1349. 10.1038/sj.onc.1202405.

    Article  CAS  PubMed  Google Scholar 

  10. Mok SC, Bonome T, Vathipadiekal V, Bell A, Johnson ME, Wong KK, Park DC, Hao K, Yip DK, Donninger H, Ozbun L, Samimi G, Brady J, Randonovich M, Pise-Masison CA, Barrett JC, Wong WH, Welch WR, Berkowitz RS, Birrer MJ: A gene signature predictive for outcome in advanced ovarian cancer identifies a survival factor: microfibril-associated glycoprotein 2. Cancer Cell. 2009, 16: 521-532. 10.1016/j.ccr.2009.10.018.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Pise-Masison CA, Radonovich M, Dohoney K, Morris JC, O’Mahony D, Lee MJ, Trepel J, Waldmann TA, Janik JE, Brady JN: Gene expression profiling of ATL patients: compilation of disease-related genes and evidence for TCF4 involvement in BIRC5 gene expression and cell viability. Blood. 2009, 113: 4016-4026. 10.1182/blood-2008-08-175901.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Ruckes T, Saul D, Van Snick J, Hermine O, Grassmann R: Autocrine antiapoptotic stimulation of cultured adult T-cell leukemia cells by overexpression of the chemokine I-309. Blood. 2001, 98: 1150-1159. 10.1182/blood.V98.4.1150.

    Article  CAS  PubMed  Google Scholar 

  13. Sasaki H, Nishikata I, Shiraga T, Akamatsu E, Fukami T, Hidaka T, Kubuki Y, Okayama A, Hamada K, Okabe H, Murakami Y, Tsubouchi H, Morishita K: Overexpression of a cell adhesion molecule, TSLC1, as a possible molecular marker for acute-type adult T-cell leukemia. Blood. 2005, 105: 1204-1213. 10.1182/blood-2004-03-1222.

    Article  CAS  PubMed  Google Scholar 

  14. Tattermusch S, Skinner JA, Chaussabel D, Banchereau J, Berry MP, McNab FW, O’Garra A, Taylor GP, Bangham CR: Systems biology approaches reveal a specific interferon-inducible signature in HTLV-1 associated myelopathy. PLoS Pathog. 2012, 8: e1002480-10.1371/journal.ppat.1002480.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Dowling D, Nasr-Esfahani S, Tan CH, O’Brien K, Howard JL, Jans DA, Purcell DF, Stoltzfus CM, Sonza S: HIV-1 infection induces changes in expression of cellular splicing factors that regulate alternative viral splicing and virus production in macrophages. Retrovirology. 2008, 5: 18-10.1186/1742-4690-5-18.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Przychodzen B, Jerez A, Guinta K, Sekeres MA, Padgett R, Maciejewski JP, Makishima H: Patterns of missplicing due to somatic U2AF1 mutations in myeloid neoplasms. Blood. 2013, 122: 999-1006. 10.1182/blood-2013-01-480970.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Germann S, Gratadou L, Dutertre M, Auboeuf D: Splicing programs and cancer. J Nucleic Acids. 2012, 2012: 269570-10.1155/2012/269570.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Cavallari I, Rende F, D’Agostino DM, Ciminale V: Converging strategies in expression of human complex retroviruses. Viruses. 2011, 3: 1395-1414. 10.3390/v3081395.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Ciminale V, Pavlakis GN, Derse D, Cunningham CP, Felber BK: Complex splicing in the human T-cell leukemia virus (HTLV) family of retroviruses: novel mRNAs and proteins produced by HTLV type I. J Virol. 1992, 66: 1737-1745.

    PubMed Central  CAS  PubMed  Google Scholar 

  20. Usui T, Yanagihara K, Tsukasaki K, Murata K, Hasegawa H, Yamada Y, Kamihira S: Characteristic expression of HTLV-1 basic zipper factor (HBZ) transcripts in HTLV-1 provirus-positive cells. Retrovirology. 2008, 5: 34-10.1186/1742-4690-5-34.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Yoshida M, Satou Y, Yasunaga J, Fujisawa J, Matsuoka M: Transcriptional control of spliced and unspliced human T-cell leukemia virus type 1 bZIP factor (HBZ) gene. J Virol. 2008, 82: 9359-9368. 10.1128/JVI.00242-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Bakker A, Li X, Ruland CT, Stephens DW, Black AC, Rosenblatt JD: Human T-cell leukemia virus type 2 Rex inhibits pre-mRNA splicing in vitro at an early stage of spliceosome formation. J Virol. 1996, 70: 5511-5518.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Li M, Kannian P, Yin H, Kesic M, Green PL: Human T lymphotropic virus type 1 regulatory and accessory gene transcript expression and export are not rex dependent. AIDS Res Hum Retroviruses. 2012, 28: 405-410. 10.1089/aid.2011.0130.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Nakano K, Watanabe T: HTLV-1 Rex: the courier of viral messages making use of the host vehicle. Front Microbiol. 2012, 3: 330-10.3389/fmicb.2012.00330.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Matsuoka E, Usuku K, Jonosono M, Takenouchi N, Izumo S, Osame M: CD44 splice variant involvement in the chronic inflammatory disease of the spinal cord: HAM/TSP. J Neuroimmunol. 2000, 102: 1-7. 10.1016/S0165-5728(99)00155-1.

    Article  CAS  PubMed  Google Scholar 

  27. Horiuchi S, Yamamoto N, Dewan MZ, Takahashi Y, Yamashita A, Yoshida T, Nowell MA, Richards PJ, Jones SA, Yamamoto N: Human T-cell leukemia virus type-I Tax induces expression of interleukin-6 receptor (IL-6R): Shedding of soluble IL-6R and activation of STAT3 signaling. Int J Cancer. 2006, 119: 823-830. 10.1002/ijc.21918.

    Article  CAS  PubMed  Google Scholar 

  28. Horiuchi S, Koyanagi Y, Tanaka Y, Waki M, Matsumoto A, Zhou YW, Yamamoto M, Yamamoto N: Altered interleukin-2 receptor alpha-chain is expressed in human T-cell leukaemia virus type-I-infected T-cell lines and human peripheral blood mononuclear cells of adult T-cell leukaemia patients through an alternative splicing mechanism. Immunology. 1997, 91: 28-34. 10.1046/j.1365-2567.1997.00236.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Sibon D, Gabet AS, Zandecki M, Pinatel C, Thete J, Delfau-Larue MH, Rabaaoui S, Gessain A, Gout O, Jacobson S, Mortreux F, Wattel E: HTLV-1 propels untransformed CD4 lymphocytes into the cell cycle while protecting CD8 cells from death. J Clin Invest. 2006, 116: 974-983. 10.1172/JCI27198.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Rabaaoui S, Zouhiri F, Lancon A, Leh H, d’Angelo J, Wattel E: Inhibitors of strand transfer that prevent integration and inhibit human T-cell leukemia virus type 1 early replication. Antimicrob Agents Chemother. 2008, 52: 3532-3541. 10.1128/AAC.01361-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Ip JY, Tong A, Pan Q, Topp JD, Blencowe BJ, Lynch KW: Global analysis of alternative splicing during T-cell activation. RNA. 2007, 13: 563-572. 10.1261/rna.457207.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Whistler T, Chiang CF, Lonergan W, Hollier M, Unger ER: Implementation of exon arrays: alternative splicing during T-cell proliferation as determined by whole genome analysis. BMC Genomics. 2010, 11: 496-

    PubMed Central  PubMed  Google Scholar 

  33. Kannagi M, Sugamura K, Sato H, Okochi K, Uchino H, Hinuma Y: Establishment of human cytotoxic T cell lines specific for human adult T cell leukemia virus-bearing cells. J Immunol. 1983, 130: 2942-2946.

    CAS  PubMed  Google Scholar 

  34. Kinoshita T, Shimoyama M, Tobinai K, Ito M, Ito S, Ikeda S, Tajima K, Shimotohno K, Sugimura T: Detection of mRNA for the tax1/rex1 gene of human T-cell leukemia virus type I in fresh peripheral blood mononuclear cells of adult T-cell leukemia patients and viral carriers by using the polymerase chain reaction. Proc Natl Acad Sci U S A. 1989, 86: 5620-5624. 10.1073/pnas.86.14.5620.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Yoshida M, Seiki M, Yamaguchi K, Takatsuki K: Monoclonal integration of human T-cell leukemia provirus in all primary tumors of adult T-cell leukemia suggests causative role of human T-cell leukemia virus in the disease. Proc Natl Acad Sci U S A. 1984, 81: 2534-2537. 10.1073/pnas.81.8.2534.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Fan L, Lagisetti C, Edwards CC, Webb TR, Potter PM: Sudemycins, novel small molecule analogues of FR901464, induce alternative gene splicing. ACS Chem Biol. 2011, 6: 582-589. 10.1021/cb100356k.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Fushimi K, Ray P, Kar A, Wang L, Sutherland LC, Wu JY: Up-regulation of the proapoptotic caspase 2 splicing isoform by a candidate tumor suppressor, RBM5. Proc Natl Acad Sci U S A. 2008, 105: 15708-15713. 10.1073/pnas.0805569105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Logette E, Wotawa A, Solier S, Desoche L, Solary E, Corcos L: The human caspase-2 gene: alternative promoters, pre-mRNA splicing and AUG usage direct isoform-specific expression. Oncogene. 2003, 22: 935-946. 10.1038/sj.onc.1206172.

    Article  CAS  PubMed  Google Scholar 

  39. Satou Y, Yasunaga J, Yoshida M, Matsuoka M: HTLV-I basic leucine zipper factor gene mRNA supports proliferation of adult T cell leukemia cells. Proc Natl Acad Sci U S A. 2006, 103: 720-725. 10.1073/pnas.0507631103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Gazon H, Lemasson I, Polakowski N, Cesaire R, Matsuoka M, Barbeau B, Mesnard JM, Peloponese JM: Human T-cell leukemia virus type 1 (HTLV-1) bZIP factor requires cellular transcription factor JunD to upregulate HTLV-1 antisense transcription from the 3’ long terminal repeat. J Virol. 2012, 86: 9070-9078. 10.1128/JVI.00661-12.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Mallinjoud P, Villemin JP, Mortada H, Polay Espinoza M, Desmet FO, Samaan S, Chautard E, Tranchevent LC, Auboeuf D: Endothelial, epithelial, and fibroblast cells exhibit specific splicing programs independently of their tissue of origin. Genome Res. 2014, 24: 511-521. 10.1101/gr.162933.113.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Dutertre M, Sanchez G, De Cian MC, Barbier J, Dardenne E, Gratadou L, Dujardin G, Le Jossic-Corcos C, Corcos L, Auboeuf D: Cotranscriptional exon skipping in the genotoxic stress response. Nat Struct Mol Biol. 2010, 17: 1358-1366. 10.1038/nsmb.1912.

    Article  CAS  PubMed  Google Scholar 

  43. Satou Y, Yasunaga J, Zhao T, Yoshida M, Miyazato P, Takai K, Shimizu K, Ohshima K, Green PL, Ohkura N, Yamaguchi T, Ono M, Sakaguchi S, Matsuoka M: HTLV-1 bZIP factor induces T-cell lymphoma and systemic inflammation in vivo. PLoS Pathog. 2011, 7: e1001274-10.1371/journal.ppat.1001274.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Murata K, Hayashibara T, Sugahara K, Uemura A, Yamaguchi T, Harasawa H, Hasegawa H, Tsuruda K, Okazaki T, Koji T, Miyanishi T, Yamada Y, Kamihira S: A novel alternative splicing isoform of human T-cell leukemia virus type 1 bZIP factor (HBZ-SI) targets distinct subnuclear localization. J Virol. 2006, 80: 2495-2505. 10.1128/JVI.80.5.2495-2505.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


The authors would like to thank Joel Lachuer and Nicolas Nazaret (ProfilExpert, Lyon, France) for providing transcriptomic facilities. This work was supported by the Ligue Nationale Contre le Cancer (Comités de l’Ain, de la Drome, de la Saone et Loire, and du Rhône), Fondation de France, Association Laurette Fugain, Association pour la Recherche sur le Cancer (ARC), Association Guillaume Espoir, and the Agence Nationale pour la Recherche (EPIVIR). M.T. and C.V. were supported by bursaries from the French Ministry of Higher Education and Science. F.M. was supported by INSERM and by the Hospices Civils de Lyon (AVIESAN CHRT program 2010). E.W. was supported by Hospices Civils de Lyon and Lyon I University (France). T.W. was supported in part by a US NIH Grant (CA140474), the American Lebanese-Syrian Associated Charities (USA), and St. Jude Children’s Research Hospital (USA).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Eric Wattel or Franck Mortreux.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MT, CV, HM, MK, and CP performed experiments and analyzed data. AG provided patient samples and analyzed data. TW provided Sudemycin E and analyzed data. EW and DA analyzed data and participated in writing the manuscript. FM designed the study, analyzed data, and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1.: List of samples. (XLSX 32 KB)


Additional file 2: Figure S1.: Tax expression in infected clones and ATLL samples. Tax mRNA was quantified by real-time qRT-PCR as described previously [29],[43]. *Statistically significant differences by Mann–Whitney U-test. (PPTX 68 KB)


Additional file 3: Figure S2.: qRT-PCR confirmation of overall gene expression. cDNA was amplified by qRT-PCR in 20 μl reactions using a QuantiFast SYBR Green PCR Kit (Qiagen) with 10 μM of each primer. Reactions were run on a Rotor-Gene 3000 (Corbett Life Science, Australia). HRPT1 gene expression was used as an internal control. A melting curve (57–95°C) was generated at the end of each run to verify primer specificity. The Pfaffl method was used for relative quantification. Primer sequences are available upon request. FC, indicates microarray-predicted positive or negative fold-changes in gene expression between infected and uninfected CD4+ cells; CD4UI, uninfected CD4+ clones; CD4I, infected CD4+ clones. (PPTX 140 KB)


Additional file 4: Table S2.: Common AEUs between untransformed-infected CD4+ clones and ATLL samples (compared with uninfected clones). (XLSX 79 KB)


Additional file 5: Table S3.: GO terms enriched by genes with common AEUs in untransformed-infected clones and ATLL samples (compared with uninfected clones). (XLSX 11 KB)


Additional file 6: Table S4.: GO terms enriched by genes displaying AEUs that discriminate between uninfected and infected clones and ATLL samples. (XLSX 34 KB)


Additional file 7: Figure S3.: GO analysis in ATLL samples and untransformed clones according to whole gene expression and AEU. The Venn-diagrams show the distribution of genes modified at the whole gene expression level (white; FC ≥ 1.2; p < 0.05), AEU level (dark grey; SI ≥ 1.2; p < 0.05), or both (light grey). Each gene set was analyzed using the DAVID browser (; KEGG pathways). Each gene set was limited to 2000 genes that displayed the highest fold-change in expression and splicing index (SI) values. The complete set of genes featured in microarrays was used as reference background. Top pathways are presented. (PPTX 62 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Thénoz, M., Vernin, C., Mortada, H. et al. HTLV-1-infected CD4+ T-cells display alternative exon usages that culminate in adult T-cell leukemia. Retrovirology 11, 119 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: