HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) versus adult T-cell leukemia/lymphoma (ATLL)
BMC Research Notes volume 14, Article number: 109 (2021)
Human T cell leukemia virus-1 (HTLV-1) infection may lead to one or both diseases including HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) or adult T cell leukemia lymphoma (ATLL). The complete interactions of the virus with host cells in both diseases is yet to be determined. This study aims to construct an interaction network for distinct signaling pathways in these diseases based on finding differentially expressed genes (DEGs) between HAM/TSP and ATLL.
We identified 57 hub genes with higher criteria scores in the primary protein–protein interaction network (PPIN). The ontology-based enrichment analysis revealed following important terms: positive regulation of transcription from RNA polymerase II promoter, positive regulation of transcription from RNA polymerase II promoter involved in meiotic cell cycle and positive regulation of transcription from RNA polymerase II promoter by histone modification. The upregulated genes TNF, PIK3R1, HGF, NFKBIA, CTNNB1, ESR1, SMAD2, PPARG and downregulated genes VEGFA, TLR2, STAT3, TLR4, TP53, CHUK, SERPINE1, CREB1 and BRCA1 were commonly observed in all the three enriched terms in HAM/TSP vs. ATLL. The constructed interaction network was then visualized inside a mirrored map of signaling pathways for ATLL and HAM/TSP, so that the functions of hub genes were specified in both diseases.
Infection with Human T cell leukemia virus-1 (HTLV-1) which normally takes an innocuous and insidious course in humans can rarely manifest as adult T cell leukemia/lymphoma (ATLL) and/or HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). HAM/TSP results in lower back pain, limb spasticity, progressive neurological decline, and urinary disturbances . Contrastingly, ATLL is characterized by unrestrained growth of T-cell precursors in blood, bone marrow, thymus, or lymph nodes. Despite thorough investigations, it is still unclear what triggers HTLV-1 infection to remain innocuous or to progress to either of these complications [2, 3]. A thorough examination of the pathogenic pathways and the disruptions in host proteins and their interactions might therefore reveal novel aspects of the disease through identification of unbalanced pathways in HTLV-1 related diseases. Identifying the key dysregulated pathways could be utilized to screen individuals with higher risk and reveal novel therapeutic approaches to prevent the progression to either of these diseases. Analyzing the cellular gene expression signature in both of these complications highlights common and specific deregulated pathways in HTLV-1 related diseases which are likely to provide novel insight into the pathogenesis and the management of HTLV-1 related diseases.
Several analytical approaches and softwares have been developed (https://pubs.acs.org/doi/abs/10.1021/acs.jproteome.5b01080, https://academic.oup.com/bioinformatics/article-abstract/35/8/1436/5102873, https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0189922)  as means to further analyze the high-throughput data derived from previous genomic studies to extract and conclude additional results which would only be available when the data from different studies were compared with each other. A well-known method to analyze previous datasets is to search for differentially expressed genes (DEGs) between normal individuals and those harboring the disease or between several complications associated with a single pathogen [5, 6].
In this study we performed the protein–protein interaction network-based analysis to determine DEGs between HAM/TSP and ATLL samples. The resulting networks and related hub genes were enriched in gene ontology for biological processes (BP) and the results and their implication is discussed.
The gene expression profile GSE19080 was acquired from the public repository database gene expression omnibus (GEO) (www.ncbi.nlm.nih.gov/geo) which includes an individual platform, GPL9686. The dataset contains the results of microarray experiments using the human ImmuneArray cDNA array. The data related to ATLL (7 specimens) and HAM/TSP (12 specimens) patients were extracted and then analyzed.
Exploration of differentially expressed genes
The Data was first normalized and preprocessed with log2 transformation. Then, the differentially expressed genes (DEGs) and their value of fold changes (FC) were acquired by GEO2R, which is an interactive web tool based on GEOquery and limma packages in R computing language. The adjusted P-value < 0.05 (calculated by FDR) was selected as the criterion for selection of DEGs. The direction of dysregulation of each for DEGs was reported as upregulation (positive logFC) and downregulation (negative logFC) compared to baseline.
Construction of protein–protein interaction (PPI) network
In order to construct the PPI network, the online STRING (Search Tool for the Retrieval of Interacting Genes) database version 10.5 was employed . The information and interactions from various biological sources, including physical interactions, functional association, high-throughput experiments, genomic context, co-expression, databases, and text-mining were considered. The cut-off criterion was set at combined score > 0.4 to analyze the PPINs.
Identification of hub genes
The network was analyzed by Network Analyzer app in Cytoscape (3.5.1) to calculate “degree” and “betweenness” centrality measures. The number of edges of a node is assigned as the degree  and the number of node visiting during moving all shortest paths is defined as betweenness centrality . The genes with higher degree and betweenness scores were selected as hub genes. These genes and their associated PPIN were visualized using Gephi version 0.9.1 .
Gene ontology analysis
Gene ontology BP was assessed by Enrichr website . Top ten major functional terms were selected based on z-scores for further analysis.
Signaling network analysis
The HTLV-1-implicated signaling network was depicted based on the KEGG and WikiPathway databases. The upregulated and downregulated genes were presented with red and blue coloring in visualized pathways, respectively.
Identification of DEGs and hub genes
The number of 1116 DEGs was recognized according to FDR < 0.05 by Benjamini–Hochberg procedure. Following the construction of the primary network by STRING, the network was analyzed by means of degree and betweenness centrality measures. 57 genes (Table 1) were selected as hub genes based on the aforementioned criteria. The logFC for each hub gene is reported in Additional file 1: Figure S1.
The relationship between hub genes was specified using STRING database. The network was constructed, as shown in Additional file 2: Figure S2. The network consists of 57 nodes and 716 edges. The size of each node was determined based on its degree value and the color of each node was specified according to their direction of deregulation with higher values of logFC represented by colors closer to red and lower values by blue.
Gene ontology biological processes
The upregulated and downregulated hub genes were individually enriched in gene ontology BP and the top functional terms were selected. Additional file 3: Figure S3a, b reveals the BP of hub genes associated with upregulated and downregulated genes, respectively. The comparison between the Gene Ontology BP of the enriched upregulated and downregulated genes revealed that positive regulation of transcription from RNA polymerase II promoter (GO: 0045944), positive regulation of transcription from RNA polymerase II promoter involved in meiotic cell cycle (GO: 0010673) and positive regulation of transcription from RNA polymerase II promoter by histone modification (GO: 1903757) are common in the samples of this analysis. Figure 1 is a circos plot demonstrating the contribution amount of proteins in each of the GO BPs (GO: 0045944, GO: 0010673 and GO: 1903757). The upregulated genes including TNF, PIK3R1, HGF, NFKBIA, CTNNB1, ESR1, SMAD2, PPARG and downregulated genes comprised of VEGFA, TLR2, STAT3, TLR4, TP53, CHUK, SERPINE1, CREB1 and BRCA1 were commonly observed in all of the three enriched terms in TSP vs. ATLL.
Signaling network displaying differentiation between HAM/TSP and ATLL
The related pathways and the connections between them are depicted in Fig. 2. The signaling network mapping for ATLL and HAM/TSP is illustrated mirror wise to provide a visual comparison of hub genes were specified in both diseases. The expression level of the following genes is increased in ATLL; TNF, AP1, NFAT, CLAM1, PI3K-AKT, PRKCB, NGF which is accompanied by spontaneous promotion of NF-κB pathway and persistent lymphocyte activation. In addition, the upregulation of PRKCB, ITGA2, IL8, and NOS2 genes is associated with inflammation, angiogenesis, cell survival, and migration pathways. Contrarily, promotion of pathways related to apoptosis and immune dysregulation are noticeable in HAM/TSP. The upregulation of genes including TP53, EGR1, Serpine1, IGFR1 can induce apoptosis, while the increase in STAT3, TLR2/4, MAP3K1, CREB1, and APP lead to disruption of the immune response.
The results of this analysis indicate high dissimalirities in regulation of cell proliferation and inflammatory pathways in ATLL and HAM/TSP. While the enriched pathways regarding RNA polymerase II  and DNA replication in this study just allude to a surge in clonal expansion or increased viral transcription of TATA-box containing 5′LTR of HTLV-1  in the provirus-harboring cells in ATLL, the identified hub genes provide further insight into the pathogenesis of both diseases.
The SRC gene is highly overexpressed in a variety of human cancers . In this study, however, ATLL was associated with lower levels of SRC expression. Counterintuitively, the loss of src homology 2 containing tyrosine phosphatase (SHP-1) is associated with spontaneous activation in HTLV-1 infected T-cells . Furthermore, another src-like tyrosine kinase LYN, which is normally upregulated via Tax in HTLV-1 cell lines , is downregulated in this study compared to HAM/TSP. This suggests that the increased LYN expression is not maintained with the loss of Tax expression in later stages of ATLL development . Further studies are needed to determine the role of src-like tyrosine kinases in ATLL and HAM/TSP.
One of the major pathways in ATLL is lymphocyte activation, which is of particular importance in ATLL compared to ACs and HAM/TSP . Transduction of lymphocyte activation signals occurs through the following ways: (i) activation of NF-kappa B, (ii) increase the AP1 gene expression via up-regulation of TNF , and (iii) enhancement the expression of NFAT gene due to up-regulation of CALM1 . The diminished expression of the aforementioned genes in HAM/TSP indicates their lesser significance compared to ATLL. Therefore, it could be speculated that NF-κB pathway is likely the major signaling pathways in ATLL compared to other forms of HTLV-1 infection, as inhibition of NF-κB by a super-repressor form of IκBα (SR-IκBα) in infected T-cells in ATLL results in cell death regardless of Tax expression . Furthermore, there was a significant disruption of genes related to phosphatidylinositol 3-kinase-Protein kinase B (PI3K-AKT) pathway in ATLL (PIK3CD, PIK3R1, and IRS1) which has been demonstrated to have a role in activation of NF-κB pathway .
Several genes like FN1, ITGA, PTK2, HGF, and CBL can affect PI3K-AKT pathway in the focal adhesion and ECM Receptor interaction pathways [22,23,24]. Moreover, PI3K-AKT itself is indirectly activated through the ROS pathway. This pathway can be stimulated through upregulation of PRKCB due to functions of NGF and HGF which are in turn mediated by PLC-γ and PLC-β, respectively [25,26,27]. The increase in NF-κB action is also important in upregulation of PTGS2, TNF, NGF, IL8, and NFKBIA in ATLL [25, 28,29,30,31], as well as NOS2 [32, 33]. The NFKBIA increase in ATLL cells however, supposedly inhibit NF-κB pathway . However, pathways such as PARP-1 that lead to NFKBIA degradation, are activated in ATLL to overcome the inhibition. The NF-κB pathway is less prominent in HAM/TSP, as the effective genes mentioned above are inversely expressed in it.
Inhibition of apoptosis is another pivotal distinction in signaling pathways between ATLL and HAM/TSP, marked by significant dowregulation of TP53 in ATLL . TP53 loss or mutation in ATLL cells reflects their aggressive proliferation and poor prognosis as described in many other cancers [34, 35]. The upregulation of SMAD2 in ATLL facilitates metastasis and is associated with TP53 mutation .
Furthermore, MDM2 can increase the expression of CDKN2A, which in turn increases TP53 . Although, upregulation of CDKN2A secondary to the induction of oncogenes is unavoidable, EZH2 functionally inhibits the role of CDKN2A in cell cycle . The decrease in expression level of TP53 is accompanied by downregulation of its target genes including Serpin1 in ATLL. Serpin1 has been linked to reduction of tumor invasion and growth . In addition, the downregulation of EGR1 tumor repressor is also observed in ATLL in this study. The increase in the expression level of EGR1 in HAM/TSP along with the higher levels of TP53 and Serpin1 compared to ATLL restricts the uncontrolled proliferation of infected cells as seen in ATLL . Furthermore, the analysis revealed PRKCB, ITGA2, IL8, and NOS2 genes to be significantly upregulated in ATLL. The afforementioned genes are associated with cell motility and angiogenesis and could be surmised to contribute to the metastasic properties of leukemic cells in ATLL.
Viral persistence cannot continue without significant dysregulation of the host immune response. The balance of immune response can tip to either immunosuppression or hypersensitivity in HTLV-1 infections, resulting in ATLL and HAM/TSP, respectively. The mild immunodeficiency almost exclusively associated with ATLL is partly the consequence of disruption of regulatory genes implicated in the immune response. In this study, the expression levels of STAT3, TLR2/4, MAP3K1, CREB1, and APP genes in ATLL have been reportedly decreased, leading to reduced antigen presentation, T-cell costimulation, and immune response against the virus [41,42,43,44,45]. The different patterns of gene expression described throughout the article further elaborates the unique immunophenotype observed in each of these two diseases.
This study provides a novel approach to gene expression analyses regarding HTLV-1 related diseases by comparing the gene expression signature between samples from HAM/TSP and ATLL patients. The results revealed distinct patterns of gene expression, especially in cell cycle regulation and immune response between the two diseases.
Further detailed studies help us understand other functions of the involved genes in the pathogenesis of HTLV-1.
Availability of data and materials
The dataset analyzed in this study [GSE19080] was accessed through gene expression omnibus public repository (www.ncbi.nlm.nih.gov/geo).
Human T-lymphotropic virus type 1
HTLV-1-associated myelopathy/tropical spastic paraparesis
Adult T-cell leukemia/lymphoma
Gene Expression Omnibus
Differentially expressed genes
Search Tool for the Retrieval of Interacting Genes
Protein–protein interaction network
Kyoto Encyclopedia of Genes and Genomes
Phosphatidylinositol 3-kinase-Protein kinase B
Bangham CR, Araujo A, Yamano Y, Taylor GP. HTLV-1-associated myelopathy/tropical spastic paraparesis. Nat Rev Dis Prim. 2015;1:15012.
Bangham CR, Ratner L. How does HTLV-1 cause adult T-cell leukaemia/lymphoma (ATL)? Curr Opin Virol. 2015;14:93–100.
Bangham CRM, Human T. Cell leukemia virus type 1: persistence and pathogenesis. Annu Rev Immunol. 2018;36:43–71.
Leung MK, Delong A, Alipanahi B, Frey BJ. Machine learning in genomic medicine: a review of computational problems and data sets. Proc IEEE. 2016;104(1):176–97.
Azimzadeh Jamalkandi S, Mozhgani SH, Gholami Pourbadie H, Mirzaie M, Noorbakhsh F, Vaziri B, et al. Systems biomedicine of rabies delineates the affected signaling pathways. Front Microbiol. 2016;7:1688.
Mozhgani SH, Zarei-Ghobadi M, Teymoori-Rad M, Mokhtari-Azad T, Mirzaie M, Sheikhi M, et al. Human T-lymphotropic virus 1 (HTLV-1) pathogenesis: a systems virology study. J Cell Biochem. 2018;119(5):3968–79.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2014;43(D1):D447–52.
Diestel R. Graph theory. Grad Texts in Math. 2005;101.
Yoon J, Blumer A, Lee K. An algorithm for modularity analysis of directed and weighted biological networks based on edge-betweenness centrality. Bioinformatics. 2006;22(24):3106–8.
Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. Icwsm. 2009;8:361–2.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.
Butler JE, Kadonaga JT. The RNA polymerase II core promoter: a key component in the regulation of gene expression. Genes Dev. 2002;16(20):2583–92.
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(19):9359–68.
Irby RB, Yeatman TJ. Role of Src expression and activation in human cancer. Oncogene. 2000;19(49):5636–42.
Migone TS, Cacalano NA, Taylor N, Yi T, Waldmann TA, Johnston JA. Recruitment of SH2-containing protein tyrosine phosphatase SHP-1 to the interleukin 2 receptor; loss of SHP-1 expression in human T-lymphotropic virus type I-transformed T cells. Proc Natl Acad Sci USA. 1998;95(7):3845–50.
Uchiumi F, Semba K, Yamanashi Y, Fujisawa J, Yoshida M, Inoue K, et al. Characterization of the promoter region of the src family gene lyn and its trans activation by human T-cell leukemia virus type I-encoded p40tax. Mol Cell Biol. 1992;12(9):3784.
Grassmann R, Aboud M, Jeang K-T. Molecular mechanisms of cellular transformation by HTLV-1 Tax. Oncogene. 2005;24(39):5976–85.
Qiao Y, He H, Jonsson P, Sinha I, Zhao C, Dahlman-Wright K. AP-1 is a key regulator of proinflammatory cytokine TNFalpha-mediated triple-negative breast cancer progression. J Biol Chem. 2016;291(10):5068–79.
Boczek T, Lisek M, Ferenc B, Zylinska L. Cross talk among PMCA, calcineurin and NFAT transcription factors in control of calmodulin gene expression in differentiating PC12 cells. Biochem Biophys Acta. 2017;1860(4):502–15.
Hironaka N, Mochida K, Mori N, Maeda M, Yamamoto N, Yamaoka S. Tax-independent constitutive IkappaB kinase activation in adult T-cell leukemia cells. Neoplasia. 2004;6(3):266–78.
Bai D, Ueno L, Vogt PK. Akt-mediated regulation of NFkappaB and the essentialness of NFkappaB for the oncogenicity of PI3K and Akt. Int J Cancer. 2009;125(12):2863–70.
Xia H, Nho RS, Kahm J, Kleidon J, Henke CA. Focal adhesion kinase is upstream of phosphatidylinositol 3-kinase/Akt in regulating fibroblast survival in response to contraction of type I collagen matrices via a beta 1 integrin viability signaling pathway. J Biol Chem. 2004;279(31):33024–34.
Liyasova MS, Ma K, Lipkowitz S. Molecular pathways: cbl proteins in tumorigenesis and antitumor immunity-opportunities for cancer treatment. Clin Cancer Res. 2015;21(8):1789–94.
Kim MS, Shutov LP, Gnanasekaran A, Lin Z, Rysted JE, Ulrich JD, et al. Nerve growth factor (NGF) regulates activity of nuclear factor of activated T-cells (NFAT) in neurons via the phosphatidylinositol 3-kinase (PI3K)-Akt-glycogen synthase kinase 3beta (GSK3beta) pathway. J Biol Chem. 2014;289(45):31349–60.
Chou WY, Chuang KH, Sun D, Lee YH, Kao PH, Lin YY, et al. Inhibition of PKC-induced COX-2 and IL-8 expression in human breast cancer cells by glucosamine. J Cell Physiol. 2015;230(9):2240–51.
Kondratyev A, Popova O, Severin S, Choladze M, Shmyrev I, Tubasheva I, et al. Nerve growth factor stimulates protein kinase C translocation in PC12 cells. FEBS Lett. 1990;264(1):75–7.
Rusanescu G, Gotoh T, Tian X, Feig LA. Regulation of Ras signaling specificity by protein kinase C. Mol Cell Biol. 2001;21(8):2650–8.
Hoesel B, Schmid JA. The complexity of NF-kappaB signaling in inflammation and cancer. Mol Cancer. 2013;12:86.
Manna SK, Ramesh GT. Interleukin-8 induces nuclear transcription factor-kappaB through a TRAF6-dependent pathway. J Biol Chem. 2005;280(8):7010–21.
Ahn SH, Park H, Ahn YH, Kim S, Cho MS, Kang JL, et al. Necrotic cells influence migration and invasion of glioblastoma via NF-kappaB/AP-1-mediated IL-8 regulation. Sci Rep. 2016;6:24552.
Lim JW, Kim H, Kim KH. Nuclear factor-kappaB regulates cyclooxygenase-2 expression and cell proliferation in human gastric cancer cells. Lab Investig J Tech Methods Pathol. 2001;81(3):349–60.
Arias-Salvatierra D, Silbergeld EK, Acosta-Saavedra LC, Calderon-Aranda ES. Role of nitric oxide produced by iNOS through NF-kappaB pathway in migration of cerebellar granule neurons induced by Lipopolysaccharide. Cell Signal. 2011;23(2):425–35.
Jafarian M, Mozhgani SH, Patrad E, Vaziri H, Rezaee SA, Akbarin MM, et al. Evaluation of INOS, ICAM-1, and VCAM-1 gene expression: a study of adult T cell leukemia malignancy associated with HTLV-1. Adv Virol. 2017;162(4):1009–15.
Whibley C, Pharoah PDP, Hollstein M. p53 polymorphisms: cancer implications. Nat Rev Cancer. 2009;9(2):95–107.
Tawara M, Hogerzeil SJ, Yamada Y, Takasaki Y, Soda H, Hasegawa H, et al. Impact of p53 aberration on the progression of adult T-cell leukemia/lymphoma. Cancer Lett. 2006;234(2):249–55.
Adorno M, Cordenonsi M, Montagner M, Dupont S, Wong C, Hann B, et al. A mutant-p53/Smad complex opposes p63 to empower TGFbeta-induced metastasis. Cell. 2009;137(1):87–98.
Ichimura K, Bolin MB, Goike HM, Schmidt EE, Moshref A, Collins VP. Deregulation of the p14ARF/MDM2/p53 pathway is a prerequisite for human astrocytic gliomas with G1-S transition control gene abnormalities. Can Res. 2000;60(2):417–24.
Jacobsen JA, Woodard J, Kee BL. EZH2 regulates cell cycle and survival in B and T lymphocyte progenitors. Am Assoc Immnol. 2016;196:122.4.
Shetty S, Shetty P, Idell S, Velusamy T, Bhandary YP, Shetty RS. Regulation of plasminogen activator inhibitor-1 expression by tumor suppressor protein p53. J Biol Chem. 2008;283(28):19570–80.
Baron V, Adamson ED, Calogero A, Ragona G, Mercola D. The transcription factor Egr1 is a direct regulator of multiple tumor suppressors including TGFbeta1, PTEN, p53, and fibronectin. Cancer Gene Ther. 2006;13(2):115–24.
Liu B, Bao P, Zhou B, Chen K, Wang R, Hu T, et al. New role of Tlr4 as a mediator for inflammation-apoptosis signaling pathway in necrotizing enterocolitis. Pediatr Res. 2011;70(S5):137.
Avalle L, Camporeale A, Camperi A, Poli V. STAT3 in cancer: a double edged sword. Cytokine. 2017;98:42–50.
van den Elsen PJ, Holling TM, Kuipers HF, van der Stoep N. Transcriptional regulation of antigen presentation. Curr Opin Immunol. 2004;16(1):67–75.
Aliprantis AO, Yang RB, Weiss DS, Godowski P, Zychlinsky A. The apoptotic signaling pathway activated by Toll-like receptor-2. EMBO J. 2000;19(13):3325–36.
Li C, Che LH, Ji TF, Shi L, Yu JL. Effects of the TLR4 signaling pathway on apoptosis of neuronal cells in diabetes mellitus complicated with cerebral infarction in a rat model. Sci Rep. 2017;7:43834.
Many thanks to the Vice Chancellor for Research, Alborz University of Medical for supporting the study.
This study was funded by Alborz University of Medical Sciences (ABZUMS).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of the up-regulated (positive logFC) and down-regulated (negative logFC) hub genes in TSP vs. ATLL and vice versa.
The PPINs formed between the identified hub genes of ATLL vs. TSP group. The node size is representative of the degree of nodes and the node color is indicative of the expression level of nodes ranging from red (upregulated genes) to blue (downregulated genes).
The most meaningful GO biological process terms in top ranks of z-score were specified for two series of upregulated and downregulated genes.
About this article
Cite this article
Zarei-Ghobadi, M., Sheikhi, M., Teymoori-Rad, M. et al. HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) versus adult T-cell leukemia/lymphoma (ATLL). BMC Res Notes 14, 109 (2021). https://doi.org/10.1186/s13104-021-05521-y