- Research Note
- Open access
- Published:
The mirror like expression of genes involved in the FOXO signaling pathway could be effective in the pathogenesis of human lymphotropic virus type 1 (HTLV-1) through disruption of the downstream pathways
BMC Research Notes volume 16, Article number: 147 (2023)
Abstract
Objectives
Human lymphotropic virus type 1 (HTLV-1) is the cause of two major diseases, ATLL and HAM/TSP in a percentage of carriers. Despite progress in understanding the pathogenesis of these two diseases, the exact pathogenesis mechanism is still not well understood. High-throughput technologies have revolutionized medical research. This study aims to investigate the mechanism of pathogenesis of these two diseases using the results of high-throughput analysis of microarray datasets.
Results
A total of 100 differentially expressed genes were found between ATLL and HAM/TSP. After constructing protein-protein network and further analyzing, proteins including ATM, CD8, CXCR4, PIK3R1 and CD2 were found as the hub ones between ATLL and HAM/TSP. Finding the modules of the subnetwork revealed the enrichment of two common pathways including FOXO signaling pathway and Cell cycle with two common genes including ATM and CDKN2D. Unlike ATLL, ATM gene had higher expressions in HAM/TSP patients. The expression of CDKN2D was increased in ATLL patients. The results of this study could be helpful for understanding the pathogenic mechanism of these two diseases in the same signaling pathways.
Introduction
Human lymphotropic virus type I (HTLV-1) is a member of the Retroviridae family, which is in the genus Deltavirus. This virus does not cause clinical symptoms in 95% of infected people, but the remaining 5% progress to adult T-cell leukemia/lymphoma (ATLL) or/and (HTLV-1)-associated myelopathy/tropical spastic paraparesis (HAM/TSP). It is estimated that between 15 and 20 million people in the world are infected with the virus. This virus is more prevalent in Japan, Iran, South Africa, and South America [1,2,3].
High-throughput studies provide possibility to find the simultaneous expression of thousands of gene. The further analyses help understand the affected pathways resulting from the changes in the expression value of the involved genes. In these types of studies, the whole transcriptome is evaluated and the differentially expressed genes (DEGs) or co-expressed genes could be identified. Therefore, the molecular disorders leading to development of a disease could be more clearly understood [4,5,6,7]. In the HTLV-1 associated diseases (ATLL and HAM/TSP), in addition to the viral mechanisms, diverse cellular signaling pathways are also involved [8, 9]. Common genes with different and even reverse expressions are involved in the development of these diseases with two different clinical signs. It is not yet clear why some people develop these diseases and others do not. A series of previous high-throughput studies compared ATLL and HAM/TSP with healthy groups as well as asymptomatic carriers (ACs), but few studies have compared these two diseases with each other [6, 10,11,12,13,14]. Since the molecular mechanism leading to the progression of the ATLL and HAM/TSP has not yet been well clarified, this study aims to identify signaling pathways and genes with significant expression differences between ATLL and HAM/TSP patients using the analysis of the microarray data. The identified genes are introduced as the possible molecular players implicated in the fate of HTLV-1 infection toward two possible diseases.
Materials and methods
Gene expression microarray dataset
In this study, gene expression profile with accession number GSE19080 and platform number 9686 through gene expression omnibus public repository (www.ncbi.nlm.nih.gov/ geo) was used. The authors performed the microarray experiments using the human ImmuneArray cDNA array. This dataset includes a total of 38 samples, of which 7 samples belong to ATLL patients, 12 samples belong to HAM/TSP patients, 11 samples included healthy carriers and 8 samples belong to normal people. In accordance with the objectives of this survey, only ATLL and HAM/TSP samples were used.
Gene expression dataset and differential expression analysis
The GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/) was employed to perform log2 transformation, recognition of DEGs, and calculation of fold change (FC). In addition, differential expression analysis between two groups of ATLL and HAM/TSP patients was visualized using package pheatmap in R 3.2.5. Based on variance of expression, first 100 genes of each group with highest values were selected to generate heatmap plot.
Protein-protein interaction network (PPIN)
The online STRING database was performed to construct the PPIN. The interactions presented in this database are based on genomic context, high-throughput experiments, co-expression, and previous knowledge (databases and text-mining) [15]. The combined score higher than 0.4 was considered as cut- off to analyze the PPIN.
Reconstruction of PPIN and centrality analysis
The PPIN obtained from the STRING database was analyzed using the Network Analyzer app in the Cytoscape (3.5.1) software, and the degree centrality criteria was calculated. These criteria reflect the amount of extension from each node to other accessible nodes. Using the information obtained from these calculations, the nodes with the highest score in the degree value were used to construct the subnetwork. Gephi version 0.9.1, an open-source network visualization and manipulation software was used for further analysis and visualization of the obtained subnetwork [16].
The identification of functional modules and gene enrichment analysis
To divide the obtained subnetwork into the corresponding modules, the fast-unfolding clustering algorithm was executed in Gephi version 0.9.1. Then, each of the obtained modules was enriched using the Enricher web tool. Significantly Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways terms were taken based on the top ten combined scores.
Study participants and quantitative real-time PCR
For data validation quantitative real-time polymerase chain reaction (RT-qPCR) was performed on the cDNA samples of 10 HAM/TSP and 10 ATLL patients to measure the expression of ATM, CDKN2D and RPLP0 using the SYBR Green qPCR Master Mix (TaKaRa, Otsu, Japan). The nucleotide sequences of designed primers were as follow: forward primer of ATM 5- CTGCTGCCGTCAACTAGAAC-3, reverse primer of ATM 5-AGGCTTGTGTTGAGGCTGAT-3, forward primer of CDKN2D 5-TGATGTCAACGTGCCTGATG-3’, reverse primer of CDKN2D 5- AGCTCCAAGGGTGTGAGAC-3. RPLP0 was utilized as a housekeeping gene in order to normalize the mRNA expression levels, as well as control error between samples [17]. The relative expressions were calculated based on the expression ratio of ATM/RPLP0 and CDKN2D /RPLP0. Data were analyzed by Prism GraphPad Software Version 8.0.2 (GraphPad software, Inc., San Diego, CA, USA). This study was approved by the medical ethics committee of Alborz University of Medical Sciences (IR.ABZUMS.REC.1398.108).
Results
Differential expression analysis (DEGs)
The heatmap plot of the first 100 genes belonging to ATLL and HAM/TSP samples were constructed. To this purpose, the mean expression level of genes was calculated for each group. Then, the top 100 genes among two sample groups were selected based on variances difference. In the constructed heatmap, the color gradient from red to green indicates the expression level of genes from the highest to the lowest expression level. The down-regulated genes are identified as green and up-regulated genes are specified as red (Fig. 1a). Table 1 shows DEGs and their expression values by LogFC. A positive LogFC indicates an increase in gene expression in HAM/TSP patients compared to ATLL patients. The statistically significant value between the two mentioned groups is showed considering adjusted p < 0.05. value. The five genes with the highest expression value were PRMT1, S100B, RAD1, APP, and SPN, and the five genes with the highest decrease in the expression value were CD48, SELPLG, KLRB1, GTF2F2, and CDC7.
Construction of PPI network and extraction of the PPI subnetwork
The PPIN was used in order to inspect the stock relationship between DEGs. In order to construct the PPIN with the highest degree of connectivity, the primary network was constructed using STRING database. Afterward, the subnetwork including the proteins with the highest degree of connectivity in the primary network was build and visualized (Fig. 1b). The subnetwork consisted of 54 nodes and 212 edges for ATLL vs. HAM/TSP patients. In the obtained subnetwork, the size of the nodes indicates the degree of each node and the color of the nodes shows the expression level of each protein. The red color is indicative of the proteins with higher expression and the blue color is representative of the proteins with higher negative expression. The five proteins with the highest degree value are ATM, CD8, CXCR4, PIK3R1 and CD2.
The identification of functional modules and pathway enrichment analysis
To find the modules of the subnetwork, the fast-unfolding clustering algorithm was applied. As a result, four modules (M1, M2, M3, M4) were obtained. The up-regulated and down-regulated genes are tagged by color. To undertake enrichment analysis among the detected modules, EnrichR web-based tool was used. The most significant GO signaling pathways in top ranks of combined score were FOXO signaling pathway and Cell cycle (Fig. 2a,b).
The identification of common genes
Figure 2c shows the involve genes in the abovementioned pathways as well as common genes between them. ATM and CDKN2D genes are common in both pathways. Unlike ATLL, ATM gene had higher expressions in HAM/TSP patients but the expression of CDKN2D was increased in ATLL patients.
ATM and CDKN2D gene expression
The mean ATM gene expression in HAM/TSP and ATLL reported as 1.21 ± 0.24 and 0.35 ± 0.17, respectively. A significant increase observed in HAM/TSP compared to ATLL (P = 0.004). The mean CDKN2D gene expression in HAM/TSP participants and ATLL illustrated to be 0.26 ± 0.25 and 1.40 ± 0.16, respectively. In ATLL, CDKN2D expressed at a significantly higher level than in HAM/TSP individuals (P = 0.003). (Supplementary 1).
Discussion
The obtained results in this study, which is based on high-throughput analysis of substantial data of the whole transcriptome, revealed that PRMT1, S100B, RAD1, APP, and SPN genes of patients with HAM/TSP have the highest expression value. Also, CD48, SELPLG, KLRB1, GTF2F2, and CDC7 genes have the lowest expression in patients with HAM/TSP. One of the advantages of this study, which has rarely been mentioned in previous studies [8, 9], is the simultaneous comparison of genes with differential expression in patients with HAM/TSP compared to patients with ATLL. This study specifically identifies that if a gene has an increased expression in a HAM/TSP disease, it has a decreased expression in ATLL disease and vice versa. This inverse or mirror like expression can be of particular interest to researchers who work on the molecular pathogenesis of the mentioned associated diseases.
In this study based on variance of expression, the first 100 genes with highest values were selected and visualized by heatmap. These 100 genes were further analyzed and their interaction whit each other were shown in a PPIN. The five genes with the highest degree value in the drawn PPIN were ATM, CD8, CXCR4, PIK3R1 and CD2 that introduced as hub genes in the PPIN. Therefore, regardless of the increase or decrease in the expression, are of potential importance for further investigation.
In this study, four modules were obtained from the bioinformatic analysis of PPIN, and each module contains its own genes depending on the physical interaction of them together. Enrichment of genes in the obtained modules determined that the cell cycle signaling pathway and the FOXO signaling pathway are common among them. It is interesting to note that ATM and CDKN2D genes are common in both pathways.
FOXO1 activation plays a role in cell cycle progression regulation. The transcription and half-life of cyclin-dependent kinase inhibitor p27KIP1 rises when FOXO1 is active and could impact cell cycle progression [18]. ATM serine/threonine kinase, symbol ATM, is a serine/threonine protein kinase that is recruited and activated by DNA double-strand breaks [19]. The activation of this protein plays an effective role in activation of the DNA damage checkpoint, leading to DNA repair, apoptosis or cell cycle arrest [20]. CDKN2D is a member of the INK4 family CKIs. This protein regulates the G1-to-S phase transition by specifically inhibiting the activity of CDK4 and CDK6 [21]. In previous studies, the role of this protein in the induction of malignancy has not been well clarified [22]. The results of this study can potentially confirm the role of this protein in the development of ATLL.
Conclusion
In this study, it was found that dysregulation of FOXO signaling pathways and Cell cycle can be effective in the development of the ATLL or HAM/TSP diseases. In addition, it was found that ATM and CDKN2D genes are common in the mentioned pathways and their expression disorder is effective in the development of these diseases. The obtained results could be helpful for understanding the pathogenic mechanism of these two diseases in the same signaling pathways.
Limitations
Further detailed studies help us understand other functions of the involved genes in the pathogenesis of HTLV-1.
Data availability
In this study, gene expression profile with accession number GSE19080 and platform number 9686 through gene expression omnibus public repository (www.ncbi.nlm.nih.gov/ geo) was used.
Abbreviations
- HTLV‑1:
-
Human T‑lymphotropic virus type 1
- HAM/TSP:
-
HTLV‑1‑associated myelopathy/tropical spastic paraparesis
- ATLL:
-
Adult T‑cell leukemia/lymphoma
- GEO:
-
Gene Expression Omnibus
- DEGs:
-
Differentially expressed genes
- STRING:
-
Search Tool for the Retrieval of Interacting Genes
- PPIN:
-
Protein–protein interaction network
References
Karimi G, Zadsar M, Pourfathollah AA. Seroprevalence and geographical distribution of human T-lymphotropic virus type 1 among volunteer blood donors in endemic areas of Iran. Virol J. 2017;14(1):14.
Cook LB, Elemans M, Rowan AG, Asquith B. HTLV-1: persistence and pathogenesis. Virology. 2013;435(1):131–40.
Verdonck K, Gonzalez E, Van Dooren S, Vandamme AM, Vanham G, Gotuzzo E. Human T-lymphotropic virus 1: recent knowledge about an ancient infection. Lancet Infect Dis. 2007;7(4):266–81.
Bebek G, Koyuturk M, Price ND, Chance MR. Network biology methods integrating biological data for translational science. Brief Bioinform. 2012;13(4):446–59.
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.
Zarei Ghobadi M, Emamzadeh R, Teymoori-Rad M, Mozhgani SH. Decoding pathogenesis factors involved in the progression of ATLL or HAM/TSP after infection by HTLV-1 through a systems virology study. Virol J. 2021;18(1):175.
Ghobadi MZ, Emamzadeh R, Afsaneh E. Exploration of mRNAs and miRNA classifiers for various ATLL cancer subtypes using machine learning. BMC Cancer. 2022;22(1):433.
Zarei-Ghobadi M, Sheikhi M, Teymoori-Rad M, Yaslianifard S, Norouzi M, Yaslianifard S, et al. HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) versus adult T-cell leukemia/lymphoma (ATLL). BMC Res Notes. 2021;14(1):109.
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.
Mozhgani SH, Zarei Ghobadi M, Norouzi M, Rahimi H, Valizadeh N, Teymoori-Rad M, et al. Signaling factors potentially associated to the pathogenesis of adult T-cell leukemia /lymphoma: a network-analysis and novel findings assessment. Virus Res. 2022;319:198875.
Zarei Ghobadi M, Mozhgani SH, Erfani Y. Identification of dysregulated pathways underlying HTLV-1-associated myelopathy/tropical spastic paraparesis through co-expression network analysis. J Neurovirol. 2021;27(6):820–30.
Shadabi S, Delrish N, Norouzi M, Ehteshami M, Habibian-Sezavar F, Pourrezaei S, et al. Comprehensive high-throughput meta-analysis of differentially expressed microRNAs in transcriptomic datasets reveals significant disruption of MAPK/JNK signal transduction pathway in adult T-cell leukemia/lymphoma. Infect Agent Cancer. 2021;16(1):49.
Ghobadi MZ, Emamzadeh R, Mozhgani SH. Deciphering microRNA-mRNA regulatory network in adult T-cell leukemia/lymphoma; the battle between oncogenes and anti-oncogenes. PLoS ONE. 2021;16(2):e0247713.
Mozhgani SH, Piran M, Zarei-Ghobadi M, Jafari M, Jazayeri SM, Mokhtari-Azad T, et al. An insight to HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) pathogenesis; evidence from high-throughput data integration and meta-analysis. Retrovirology. 2019;16(1):46.
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. 2015;43(Database issue):D447–52.
Bastian M. In: Heymann S, Jacomy M, editors. Gephi: an open source software for exploring and manipulating networks. Proceedings of the international AAAI conference on web and social media; 2009.
Aghamahdi F, Shafiee A, Rostami S, Mokhames Z, Safavi M, Yaslianifard S, et al. Comparative study of CNR1 and CNR2 cannabinoid receptors expression levels in COVID-19 patients with and without diabetes mellitus: recommendations for future research targets. Diabetes Metab Syndr. 2022;16(5):102499.
Lu H, Huang H. FOXO1: a potential target for human diseases. Curr Drug Targets. 2011;12(9):1235–44.
Canman CE, Lim DS. The role of ATM in DNA damage responses and cancer. Oncogene. 1998;17(25):3301–8.
Kurz EU, Lees-Miller SP. DNA damage-induced activation of ATM and ATM-dependent signaling pathways. DNA Repair (Amst). 2004;3(8–9):889–900.
Kim WY, Sharpless NE. The regulation of INK4/ARF in cancer and aging. Cell. 2006;127(2):265–75.
Wang Y, Jin W, Jia X, Luo R, Tan Y, Zhu X, et al. Transcriptional repression of CDKN2D by PML/RARalpha contributes to the altered proliferation and differentiation block of acute promyelocytic leukemia cells. Cell Death Dis. 2014;5:e1431.
Acknowledgements
Not applicable.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
“SY1, SY2, MM and S‑ HM designed the project. S‑ HM interpreted and wrote the manuscript. SY1, SY2, did systematic search and literature review. All authors read and approved the final manuscript.“
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was performed in accordance with the Declaration of Helsinki for all human or animal experimental investigations. This study was approved by the medical ethics committee of Alborz University of Medical Sciences (IR.ABZUMS.REC.1398.108) and informed consent obtained from participants.
Consent for publication
Not applicable.
Competing interest
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary 1
. The expression levels of ATM (a) and CDKN2D (B) in HAM/TSP and ATLL individuals.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yaslianifard, S., Movahedi, M., Yaslianifard, S. et al. The mirror like expression of genes involved in the FOXO signaling pathway could be effective in the pathogenesis of human lymphotropic virus type 1 (HTLV-1) through disruption of the downstream pathways. BMC Res Notes 16, 147 (2023). https://doi.org/10.1186/s13104-023-06423-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104-023-06423-x