- Research
- Open access
- Published:
Comprehensive bioinformatic analysis of the expression and prognostic significance of TSC22D domain family genes in adult acute myeloid leukemia
BMC Medical Genomics volume 16, Article number: 117 (2023)
Abstract
Background
TSC22D domain family genes, including TSC22D1-4, play a principal role in cancer progression. However, their expression profiles and prognostic significance in adult acute myeloid leukemia (AML) remain unknown.
Methods
The online databases, including HPA, CCLE, EMBL-EBI, GEPIA2, BloodSpot, GENT2, UCSCXenaShiny, GSCALite, cBioportal, and GenomicScape, utilized the data of TCGA and GEO to investigate gene expression, mutation, copy number variation (CNV), and prognostic significance of the TSC22D domain family in adult AML. Computational analysis of resistance (CARE) was used to explore the effect of TSC22D3 expression on drug response. Functional enrichment analysis of TSC22D3 was performed in the TRRUST Version 2 database. The STRING, Pathway Commons, and AnimalTFDB3.0 databases were used to investigate the protein–protein interaction (PPI) network of TSC22D3. Harmonizome was used to predict target genes and kinases regulated by TSC22D3. The StarBase v2.0 and CancermiRNome databases were used to predict miRNAs regulated by TSC22D3. UCSCXenaShiny was used to investigate the correlation between TSC22D3 expression and immune infiltration.
Results
Compared with normal adult hematopoietic stem cells (HSCs), the expression of TSC22D3 and TSC22D4 in adult AML tissues was markedly up-regulated, whereas TSC22D1 expression was markedly down-regulated. The expression of TSC22D1 and TSC22D3 was significantly increased in adult AML tissues compared to normal adult tissues. High TSC22D3 expression was significantly associated with poor overall survival (OS) and event-free survival (EFS) in adult AML patients. Univariate and multivariate Cox analysis showed that overexpression of TSC22D3 was independently associated with adverse OS of adult AML patients. High TSC22D3 expression had a adverse impact on OS and EFS of adult AML patients in the chemotherapy group. TSC22D3 expression correlated with drug resistance to BCL2 inhibitors. Functional enrichment analysis indicated that TSC22D3 might promote AML progression. MIR143-3p sponging TSC22D3 might have anti-leukemia effect in adult AML.
Conclusions
A significant increase in TSC22D3 expression was observed in adult AML tissues compared to normal adult HSCs and tissues. The prognosis of adult AML patients with high TSC22D3 expression was unfavorable, which could severe as a new prognostic biomarker and potential target for adult AML.
Introduction
Acute myeloid leukemia (AML) is an aggressive hematopoietic malignancy with high biological and clinical heterogeneity [1]. Despite advances made in the diagnosis and treatment of AML, the increased risk of relapse and low 5-year survival rate after diagnosis remain significant challenges [2]. Authentication of new AML biomarkers can help to clarify the pathogenesis of the disease and guide the diagnosis, treatment, and prognosis evaluation of AML [3]. TSC22D domain family genes have been extensively reported to play an essential role in tumors [4,5,6,7]. Nonetheless, their expression profiles and prognosis in adult AML remain unclear. Herein, we conducted an integrated analysis of the expression and prognostic value of TSC22D domain family genes in adult AML by using data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The flow chart of our study was shown in Fig. 1.
Materials and methods
Data retrieval and processing
Gene expression analysis
Gene expression of the TSC22D domain family in AML cell lines
Human Protein Atlas (HPA, https://www.proteinatlas.org) [8] is a comprehensive database of proteomics, transcriptomics, and systems biology data. The expression of the TSC22D domain family genes in 88 leukemia cell lines (including 38 AML cell lines) was determined and visualized using “HPA”.
Cancer Cell Line Encyclopedia (CCLE, https://www.broadinstitute.org/ccle) [9] is a multiomics online database that provides a large-scale transcriptome sequencing data for the study of human cancer cell lines. The expression data of TSC22D domain family genes in 43 AML cell lines was downloaded from the “Expression 22Q4 Public” dataset of the CCLE database and visualized by the cluster heatmap tool.
EMBL's European Bioinformatics Institute (EMBL-EBI, https://www.ebi.ac.uk) [10] is an integrated bioinformatics research database. The expression of the TSC22D domain family genes in 16 AML cell lines was determined and visualized using “EMBL-EBI”.
Gene expression of the TSC22D domain family in adult AML tissues and CD34 positive hematopoietic stem cells (HSCs) from normal adult bone marrow tissues
BloodSpot (http://servers.binf.ku.dk/bloodspot/) [11] is an online open data analysis platform that provides gene expression and survival prognosis data from TCGA and GEO databases. Gene expression data of the TSC22D domain family in adult AML tissues and CD34 positive HSCs from normal adult bone marrow tissues was downloaded from the “Normal hematopoiesis with AMLs” dataset and the “Bloodpool: AML samples with normal cells” dataset of the BloodSpot database.
Gene Expression database of Normal and Tumor tissues 2 (GENT2, http://gent2.appex.kr) [12] integrates publicly available expression profile microarray data from the GEO database to compare and analyze gene expression in normal and cancer patient tissues. Gene expression data of the TSC22D domain family in 2802 adult AML tissues and 17 CD34 positive HSCs from normal adult bone marrow tissues was downloaded from the “ GPL570 platform (HG-U133_Plus_2)” of the GENT2 database (See excel sheet 1 in the Additional file 1).
Gene expression of the TSC22D domain family in adult AML tissues and normal adult tissues
Gene Expression Profling Interactive Analysis 2 (GEPIA2, http://gepia2.cancer-pku.cn/) [13] is an updated and enhanced online publicly accessible database based on TCGA and Genotype-Tissue Expression (GTEx) databases for tumor and normal samples for gene expression analysis. The expression of the TSC22D domain family genes in 173 TCGA-LAML tissues and 70 GTEx-Normal tissues was compared and visualized using “GEPIA2”.
Gene expression data of the TSC22D domain family in 542 adult AML tissues and 73 normal adult.bone marrow tissues was downloaded from the “ Leukemia MILE Study” dataset (GSE13159) of the BloodSpot database.
Gene expression data of the TSC22D domain family in 2802 adult AML tissues and 134 normal adult bone marrow tissues was downloaded from the “ GPL570 platform (HG-U133_Plus_2)” of the GENT2 database (See excel sheet 1 in the Additional file 1).
Survival analysis
We retrieved and analyzed the RNAseq gene expression data of the TSC22D domain family and the corresponding clinical prognostic data in the GEPIA2, Bloodspot, GSCALite ( http://bioinfo.life.hust.edu.cn/web/GSCALite/) [14], UCSCXenaShiny (https://shiny.hiplot-academic.com/ucsc-xena-shiny) [15], cBioportal (https://www.cbioportal.org) [16], and GenomicScape (http://genomicscape.com/) [17] databases. Furthermore, RNASeq (RNA-seq V2 RSEM) gene expression data of TSC22D3 and the corresponding clinical prognostic data was downloaded from the “ TCGA-LAML, NEJM 2013” [18] dataset of the cBioPortal database (See excel sheet 2 in the Additional file 2). Then adult patients with AML were stratified into a low expression group and a high expression group based on TSC22D3 mRNA median expression. We explored the relationship between TSC22D3 expression and clinical parameters and performed the analyses of OS, EFS, and univariate and multivariate Cox OS.
Effect of TSC22D3 expression on drug response
Computational analysis of resistance [19] (CARE, http://care.dfci.harvard.edu/) is used to identify genomes and biomarkers of response to targeted therapies. A positive CARE score indicated that gene expression was associated with drug sensitivity, whereas a negative CARE score indicated drug resistance.
Data of the correlation between TSC22D3 expression and drug response was downloaded from the “Cancer Genome Project (CGP)” dataset and the “Cancer Therapeutics Response Portal (CTRP)” dataset of the CARE database and visualized by the arc link tool.
Gene mutation and copy number variation (CNV) analysis of TSC22D3
Gene mutation data and the corresponding survival data of of TSC22D3 in adult AML from the “ TCGA-LAML, PanCancer Atlas” dataset was analyzed and visualize using “cBioPortal”. And CNV data and the corresponding survival data of TSC22D3 in adult AML was analyzed and visualized using “GSCALite”.
Functional enrichment analysis of TSC22D3
Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining Version 2 (TRRUST Version 2, http://www.grnpedia.org/trrust/) [20] is an online, open database of human and mouse transcriptional regulatory networks. Gene ontology (GO) biological process, disease ontology (DO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway data associated with human TSC22D3 transcription factor (TF) was downloaded from “TRRUST Version 2” and then visualized by the bar with a color gredient tool.
Protein–protein interaction (PPI) analysis of TSC22D3
STRING (https://string-db.org) [21] is an online open database aimed at providing customized protein–protein networks.
Pathway Commons (http://www.pathwaycommons.org) [22] is an integrated platform of multiple database for predicting protein–protein interactions.
AnimalTFDB3.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB/) [23] is an online database aimed at providing the most comprehensive and accurate information for animal (including human) TFs and cofactors.
The relationship between TSC22D3 and other proteins was predicted and visualized using “String,” “Pathway Commons,” and “AnimalTFDB3.0.” The expression of potential protein in 173 TCGA-LAML tissues and 70 GTEx-Normal tissues was analyzed using “UCSCXenaShiny”. The correlation between TSC22D3 protein and potential protein was analyzed and visualized using “UCSCXenaShiny”. The effect of potential protein on OS of adult AML patients was analyzed and visualized using “GenomicScape”.
Analysis of TSC22D3 regulated target genes and kinases
Harmonizome (http://amp.pharm.mssm.edu/Harmonizome) [24] integrates many publicly available online databases to predict the functions of genes or proteins.
Data of TSC22D3 regulated target genes was downloaded from the “CHEA Transcription Factor Targets” dataset, the “ENCODE Transcription Factor Targets” dataset, and the “JASPAR Predicted Transcription Factor Targets” dataset of the Harmonizome database and visualized using the jvenn tool.
The expression of potential target gene in 173 TCGA-LAML tissues and 70 GTEx-Normal tissues was analyzed and visualized using “UCSCXenaShiny”. The correlation between TSC22D3 and potential target gene was analyzed and visualized using “UCSCXenaShiny”. The effect of potential target gene on OS of adult AML patients was analyzed and visualized using “GenomicScape”.
Data of the top 20 predicted kinases with a high Z score regulated by TSC22D3 was downloaded from the Harmonizome database and visualized using the circular heatmap tool. The expression of predicted kinases in 173 TCGA-LAML tissues and 70 GTEx-Normal tissues was analyzed and visualized using “UCSCXenaShiny”. The correlation between TSC22D3 and predicted kinases was analyzed and visualized using “UCSCXenaShiny”. The effect of predicted kinases on OS of adult AML patients was analyzed and visualized using “UCSCXenaShiny”.
Analysis of TSC22D3 regulated miRNAs
StarBase v2.0 ( https://starbase.sysu.edu.cn/) [25] integrates multiple online microRNA (miRNA) databases to explore miRNA interactions.
CancerMIRNome ( http://bioinfo.jialab-ucr.org/CancerMIRNome) [26] is an online database for interactive analysis and visualization of the miRNome spectrum in human cancer.
Data of the predicted miRNAs regulated by TSC22D3 was downloaded from the “PITA” dataset, the “microT” dataset, the “miRmap”dataset, the “miRanda”dataset, the “PicTar”dataset, and the “TargetScan” dataset of the StarBase v2.0 database and visualized using the jvenn tool. The correlation between the TSC22D3 expression and potential miRNA was analyzed and visualized using “StarBase v2.0” and “UCSCXenaShiny”. The effect of potential miRNA on OS of adult AML patients was analyzed and visualized using “UCSCXenaShiny” and “CancerMIRNome”. DO and KEGG pathway analysis of potential miRNA was performed and visualized using “CancerMIRNome”.
Immune infiltration analysis of TSC22D3
Data of the correlation between the TSC22D3 expression and immune cell infiltration in adult AML by using the “CIBERSORT” algorithm, the “QUANTISEQ” algorithm, the “MCPCOUNTER” algorithm, the “EPIC” algorithm, and the “XCELL” algorithm was downloaded from the UCSCXenaShiny database and visualized using the correlation analysis tool.
Data analysis and visualization
Wilcoxon rank-sum test was used for comparative analysis of gene expression. Kaplan Meier survival analysis (including OS and EFS) was performed using the log-rank test. Univariate and multivariate Cox's survival analysis was performed using SPSS software version 21.0. Graph Pad Prism version 8 was used for chi-square test analysis. P value < 0.05 indicated a significance level.
Data visualization was performed using the cluster heatmap tool, the circular heatmap tool, box tool, KM survival curve tool, jvenn tool, bar with color gredient tool, arc link tool, and correlation analysis tool from the website (http://www.bioinformatics.com.cn).
Results
Analysis of the expression of the TSC22D domain family genes in AML cell lines, normal adult HSCs, adult AML tissues and normal adult tissues
Three different databases, including “HPA”, “CCLE”, and EMBL-EBI”, demonstrated that TSC22D domain family genes were abnormally expressed in AML cell lines at different levels (Fig. 2A–C).
Then we explored the expression of the TSC22D domain family genes in adult AML tissues and normal adult HSCs using TCGA and GEO data from the BloodSpot and GENT2 databases. These results showed that the expression of TSC22D3 and TSC22D4 was significantly upregulated in adult AML tissues relative to normal adult HSCs, whereas the expression trend of TSC22D1 was the opposite (P < 0.05) (Fig. 3A–C).
We subsequently investigated gene expression of the TSC22D domain family in adult AML tissues and normal adult tissues utilizing TCGA data from the GEPIA2 database, and GEO data from the BloodSpot and GENT2 databases. The results revealed that the expression of TSC22D1 and TSC22D3 was markedly increased in adult AML tissues compared with normal adult tissues (P < 0.05) (Fig. 4A–C).
Survival analysis according to the expression of the TSC22D domain family genes in adult AML
Survival analysis was performed using TCGA data from the GEPIA2, Bloodspot, GSCALite, UCSCXenaShiny, and cBioportal databases and GEO data from the GenomicScape database. Amusingly, only TSC22D3 expression was of survival prognostic significance in adult AML. However, other members of the TSC22D family genes had little effect on OS of adult AML patients (See Table 1).
Morever, we found that high TSC22D3 expression was significantly correlated with white blood cell (WBC) counts (> 20 × 10^9/L), bone marrow (BM) blasts (> 70%), FAB M1 subtype, FAB M5 subtype, and positive NPM1 mutation (P < 0.05). Low TSC22D3 expression was significantly correlated with the FAB M2 subtype and the FAB M3 subtype (P < 0.05) (See Table 2).
We found that high TSC22D3 expression significantly affected OS and EFS of adult AML patients (P < 0.05) (Fig. 5A, D). Both univariate and multivariate COX regression analysis showed that the increased mortality in adult AML patients was significantly associated with over 60 years old, cytogenetics, DNMT3A positive mutation, TP53 positive mutation, treatment type (chemotherapy), and high TSC22D3 expression (P < 0.05) (See Table 3). High TSC22D3 expression had a detrimental effect on OS and EFS of adult AML patients in the chemotherapy group (P < 0.05) (Fig. 5B, E). However, high TSC22D3 expression had no effect on OS and EFS of adult AML patients in the transplantation group (P > 0.05) (Fig. 5C, F).
Analysis of the effect of TSC22D3 expression on drug response
Analysis of the effect of TSC22D3 expression on drug response using the “CGP” dataset and the “CTRP” dataset demonstrated that TSC22D3 expression was significantly associated with drug resistance to BCL2 inhibitors (Fig. 5G, H).
Gene mutation and CNV analysis of TSC22D3
Gene mutation rate of TSC22D3 was 8%, and TSC22D3 gene mutation did not affect the OS of adult AML patients (Fig. 6A, C). An the incidence of CNV of TSC22D3 was low in adult AML and did not affect the OS of adult AML patients (Fig. 6B, D).
Functional enrichment analysis of TSC22D3
The results showed that TSC22D3 has many biological functions, including response to DNA damage stimulus, G1 phase of mitotic cell cycle, regulation of cell proliferation, cell cycle arrest, and response to drug, etc. (Fig. 7A). DO analysis revealed that TSC22D3 was involved in tumors, including myeloid leukemia (Fig. 7A). Furthermore, KEGG pathway analysis indicated that TSC22D3 was involved in the regulation of multiple signaling pathways (Fig. 7A).
PPI analysis of TSC22D3
The consistent analysis of the STRING, Pathway Commons, and AnimalTFDB3.0 databases indicated that TSC22D3 interacted with FOS and SCNN1B (Fig. 7B–D). And FOS has been extensively reported in tumor progression. Then the analysis of gene expression, correlation, and survival prognosis of FOS in adult AML was performed in TCGA and GEO datasets. The results revealed that FOS was significantly increased in adult AML (P < 0.05) (Fig. 7E). TSC22D3 was positively correlated with FOS (P < 0.05) (Fig. 7F). High FOS expression was associated with unfavorable OS of 78 adult AML patients (P < 0.05) (Fig. 7G).
Analysis of TSC22D3 regulated target genes and kinases
The results showed that TSC22D3 might regulate CREB1 (Fig. 8A). Gene expression analysis indicated that CREB1 was significantly elevated in adult AML (P < 0.05) (Fig. 8B). Unexpectedly, TSC22D3 had statistical no correlation with CREB1 (P > 0.05) (Fig. 8C). However, high CREB1 expression had an adverse impact on OS of 162 adult AML patients (P < 0.05) (Fig. 8D).
We investigated the top 20 kinases with high Z score regulated by TSC22D3 (Fig. 8E) and analyzed their expression, as well as the correlation of TSC22D3, and survival prognosis in adult TCGA-LAML (Additional file 3: See Table S1). The results revealed that the expression of MAP4K1, MAP2K3, TYK2, and STK10 was markedly up-regulated in adult AML and significantly positively correlated with TSC22D3 (P < 0.05) (Fig. 8F, G). And These kinases had an unfavorable effect on OS of 161 adult AML patients (P < 0.05) (Fig. 8H).
Analysis of miRNAs regulated by TSC22D3
Analysis of six different miRNA datasets revealed that TSC22D3 might be a possible target gene of MIR143-3p (Fig. 9A). TSC22D3 was negatively correlated with MIR143-3p (P < 0.05) (Fig. 9B, C). High expression of MIR143-3p was a favorable prognostic factor for OS of adult AML patients (P < 0.05) (Fig. 9D, E). DO and KEGG analysis indicated that MIR143-3p was involved in bone marrow cancer, including myeloid leukemia (Fig. 9F, G).
Immune infiltration analysis of TSC22D3
Analysis of the correlation of TSC22D3 expression and immune cell infiltration in adult AML by using five different algorithms showed that TSC22D3 expression was significantly associated with monocyte/macrophage (P < 0.05) (Fig. 10A–E).
Discussion
TSC22D domain family genes, including TSC22D1-4, belong to the leucine zipper TF family and have been reported to be involved in regulating cell proliferation and differentiation [27]. TSC22D1, also called transforming growth factor-β-stimulated clone-22, was reported to play a tumor suppressor role in tumors [28]. TSC22D2 depends on the TSC22D2-PKM2-CyclinD1 regulatory axis to inhibit tumor cell growth in colorectal cancer [29]. TSC22D3, also known as glucocorticoid-induced leucine zipper (GILZ), can promote or suppress tumor growth, depending on the type of tumor and its microenvironment. TSC22D3 plays a dual role in tumors: it not only exerts a tumor-promoting effect by influencing the immune system and tumor microenvironment but also inhibits tumor growth by inducing apoptosis or suppressing the proliferation of cancer cells [30]. TSC22D4, also known as THG-1, was reported to promote esophageal squamous cell carcinoma cell tumorsphere growth [31].
In our study, TCGA and GEO data was used to investigate the expression of TSC22D domain family genes and their prognostic significance in adult AML. These results showed that the expression of TSC22D1 and TSC22D3 was markedly increased in adult AML tissues. Stunningly, it was TSC22D3, not other TSC22D family genes, that had prognostic significance for OS of adult AML patients. Therefore, we focused on the possible role of TSC22D3 in adult AML. Our study revealed that adult AML patients with high expression of TSC22D3 had adverse OS and EFS. And overexpression of TSC22D3 was an independently survival prognostic factor in adult AML patients. Subgroup survival analysis according to treatment type also showed that high TSC22D3 expression was significantly associated with unfavorable OS and EFS in the chemotherapy group. However, we found no effect of TSC22D3 expression on OS and EFS of adult AML patients in the transplantation group. This suggested that transplantation might overcome the disadvantages of TSC22D3 expression. Furthermore, we found that high TSC22D3 expression was significantly associated with high WBC counts, high BM blasts, FAB M1 subtype, FAB M5 subtype, and positive mutation of NPM1. This partly explained why high expression of TSC22D3 was associated with a poor survival prognosis in adult AML.
Hyperactivation of BCL2 is associated with the development, progression, prognosis, and resistance to chemotherapy in AML. BCL2 inhibitors including Venetoclax have been applied in the clinical treatment of AML [32]. However, with the widespread use of Venetoclax, drug resistance has gradually emerged in AML patients, especially in the relapsed/refractory AML patients. Preclinical and clinical studies have partially unraveled the mechanism of drug resistance to Venetoclax [33]. The clinical application of BCL2 inhibitors still faces many challenges, which may be relevant to the fact that the complex mechanism of drug resistance has not been fully unraveled. Fascinatingly, our study showed that TSC22D3 expression was significantly correlated with resistance to BCL2 inhibitors. This might be one of the reasons for drug resistance of BCL2 inhibitors, and it was worth further exploring the underlying mechanism of drug resistance mediated by TSC22D3.c-Fos has been reported to play crucial parts in the maintenance and proliferation of AML [34]. Interestingly, our study indicated that TSC22D3 might transcriptionally up-regulate the expression of FOS, which might play a certain role in AML progression. TSC22D3 promoted tumor cell proliferation by regulating AKT kinase [35]. And hyperactivity of the kinases was involved in cancer progression. Therefore, we analyzed the kinases regulated by TSC22D3. Our study indicated that TSC22D3 might transcriptionally activate the kinases of MAP4K1, MAP2K3, TYK2, and STK10. MAP4K1, as an oncogene, promoted AML progression by regulating the cell cycle through the MAPK pathway [36]. MAP2K3 promoted tumor progression by regulating tumor cell migration and invasion through the JNK signaling pathway [37]. Dysregulated activation of TYK2 in cancers may lead to hyperactive JAK/STATs signal, which may play an important role in the occurrence and development of cancers [38]. The prognosis of AML patients with expressing high levels of STK10 was poor, which could severe as a new prognostic biomarker for AML [39].
MIR143-3p has been reported to function as a tumor suppressor [40]. Our study indicated that MIR143-3p might exhibit anti-leukemic effect by downregulating the expression of TSC22D3. TSC22D3 has been reported to be involved in the supervision of the cell cycle, differentiation, and apoptosis of immune cells [41]. TSC22D3 may play an anti-inflammatory and immunosuppressive role in tumor development. Activation of the immunosuppressive TSC22D3 TF in dendritic cells can result in treatment failure [42]. Overexpression of TSC22D3 subverted therapy-induced anticancer immuno-surveillance [43]. As a TF, TSC22D3 may mediate the immunosuppressive and anti-inflammatory effects of T cells and macrophages by inhibiting nuclear factor-κB (NF-κB)-dependent transcription [44, 45]. Furthermore, TSC22D3 played a significant role in tumor progression by mediating the increase in cell quantity and activity of Treg cells through the TGF-β signaling pathway [46, 47]. TSC22D3 could play an indispensable role in the tumor microenvironment by influencing all immune system cells that infiltrated the tumor microenvironment [30]. In addition, TSC22D3 may serve as a pivotal regulator of T cell predysfunction [48]. Recent research shows that the proliferation, survival, and drug resistance of AML cells may be sustained and modulated by the bone marrow immunosuppressive microenvironment [49]. Our study showed a significantly positive correlation between monocyte/macrophage and TSC22D3 expression. How did TSC22D3 regulate monocyte/macrophage needed further study in AML immune microenvironment.
To sum up, TSC22D3 might be involved in AML progression through multiple mechanisms, including the regulation of target genes, kinases, signaling pathways, drug resistance, and immune cell infiltration. MIR143-3p sponging TSC22D3 might exhibit anti-leukemic effect in adult AML. Our study extended our understanding of TSC22D3 as a novel prognostic factor in adult AML and its potential role in AML.
Availability of data and materials
The datasets provided for this study can be found and accessed in online databases. These online databases were accessible from the following addresses. HPA, https://www.proteinatlas.org; CCLE, https://www.broadinstitute.org/ccle; EMBL-EBI, https://www.ebi.ac.uk; BloodSpot, http://servers.binf.ku.dk/bloodspot/; GENT2, http://gent2.appex.kr; GEPIA2, http://gepia2.cancer-pku.cn/; GSCALite, http://bioinfo.life.hust.edu.cn/web/GSCALite/; UCSCXenaShiny, https://shiny.hiplot-academic.com/ucsc-xena-shiny; cBioportal, https://www.cbioportal.org; GenomicScape, http://genomicscape.com/; CARE, http://care.dfci.harvard.edu/; TRRUST Version 2, http://www.grnpedia.org/trrust/; STRING, https://string-db.org; Pathway Commons, http://www.pathwaycommons.org; AnimalTFDB3.0, http://bioinfo.life.hust.edu.cn/AnimalTFDB/; Harmonizome, http://amp.pharm.mssm.edu/Harmonizome; StarBase v2.0, https://starbase.sysu.edu.cn/; CancerMIRNome, http://bioinfo.jialab-ucr.org/CancerMIRNome.
Abbreviations
- AML or LAML:
-
Acute myeloid leukemia
- HPA:
-
Human Protein Atlas
- CCLE:
-
Cancer Cell Line Encyclopedia
- EMBL-EBI:
-
EMBL’s European Bioinformatics Institute
- GEPIA2:
-
Gene Expression Profling Interactive Analysis 2
- GENT2:
-
Gene Expression database of Normal and Tumor tissues 2
- CARE:
-
Computational analysis of resistance
- CGP:
-
Cancer Genome Project
- CTRP:
-
Cancer Therapeutics Response Portal
- TRRUST:
-
Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining
- CNV:
-
Copy number variation
- TCGA:
-
The Cancer Genome Atlas
- GEO:
-
Gene Expression Omnibus
- PPI:
-
Protein–protein interaction
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- GTEx:
-
Genotype-Tissue Expression
- TF:
-
Transcription factor
- OS:
-
Overall survival
- EFS:
-
Event-free survival
- GILZ:
-
Glucocorticoid-induced leucine zipper
- HSC:
-
Hematopoietic stem cell
- GO:
-
Gene ontology
- DO:
-
Disease ontology
- miRNA:
-
MicroRNA
- WBC:
-
White blood cell
- BM:
-
Bone marrow
- FAB:
-
French–American–British
References
Döhner H, Weisdorf DJ, Bloomfield CD. Acute myeloid leukemia. N Engl J Med. 2015;373(12):1136–52. https://doi.org/10.1056/NEJMra1406184.
Doucette K, Karp J, Lai C. Advances in therapeutic options for newly diagnosed, high-risk AML patients. Ther Adv Hematol. 2021;12:20406207211001136. https://doi.org/10.1177/20406207211001138.
Prada-Arismendy J, Arroyave JC, Röthlisberger S. Molecular biomarkers in acute myeloid leukemia. Blood Rev. 2017;31(1):63–76. https://doi.org/10.1016/j.blre.2016.08.005.
Meijer D, Jansen MP, Look MP, et al. TSC22D1 and PSAP predict clinical outcome of tamoxifen treatment in patients with recurrent breast cancer. Breast Cancer Res Treat. 2009;113(2):253–60. https://doi.org/10.1007/s10549-008-9934-3.
Xiao L, Wei F, Liang F, et al. TSC22D2 identified as a candidate susceptibility gene of multi-cancer pedigree using genome-wide linkage analysis and whole-exome sequencing. Carcinogenesis. 2019;40(7):819–27. https://doi.org/10.1093/carcin/bgz095.
Qadir F, Aziz MA, Sari CP, et al. Transcriptome reprogramming by cancer exosomes: identification of novel molecular targets in matrix and immune modulation. Mol Cancer. 2018;17(1):97. https://doi.org/10.1186/s12943-018-0846-5.
Zhang X, Koga N, Suzuki H, et al. Promotion of cellular senescence by THG-1/TSC22D4 knockout through activation of JUNB. Biochem Biophys Res Commun. 2020;522(4):897–902. https://doi.org/10.1016/j.bbrc.2019.11.145.
Thul PJ, Lindskog C. The human protein atlas: a spatial map of the human proteome. Protein Sci. 2018;27(1):233–44. https://doi.org/10.1002/pro.3307.
Ghandi M, Huang FW, Jané-Valbuena J, et al. Next-generation characterization of the Cancer Cell Line Encyclopedia. Nature. 2019;569(7757):503–8. https://doi.org/10.1038/s41586-019-1186-3.
Madeira F, Park YM, Lee J, et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019;47(W1):W636–41. https://doi.org/10.1093/nar/gkz268.
Bagger FO, Kinalis S, Rapin N. BloodSpot: a database of healthy and malignant haematopoiesis updated with purified and single cell mRNA sequencing profiles. Nucleic Acids Res. 2019;47(D1):D881–5. https://doi.org/10.1093/nar/gky1076.
Park SJ, Yoon BH, Kim SK, et al. GENT2: an updated gene expression database for normal and tumor tissues. BMC Med Genom. 2019;12(Suppl 5):101. https://doi.org/10.1186/s12920-019-0514-7.
Tang Z, Kang B, Li C, et al. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–60. https://doi.org/10.1093/nar/gkz430.
Liu CJ, Hu FF, Xia MX, et al. GSCALite: a web server for gene set cancer analysis. Bioinformatics. 2018;34(21):3771–2. https://doi.org/10.1093/bioinformatics/bty411.
Wang S, Xiong Y, Zhao L, et al. UCSCXenaShiny: an R/CRAN package for interactive analysis of UCSC Xena data. Bioinformatics. 2022;38(2):527–9. https://doi.org/10.1093/bioinformatics/btab561.
Gao J, Aksoy BA, Dogrusoz U, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1. https://doi.org/10.1126/scisignal.2004088.
Kassambara A, Rème T, Jourdan M, et al. GenomicScape: an easy-to-use web tool for gene expression data analysis. Application to investigate the molecular events in the differentiation of B cells into plasma cells. PLoS Comput Biol. 2015;11(1):e1004077. https://doi.org/10.1371/journal.pcbi.1004077.
Cancer Genome Atlas Research Network; Ley TJ, Miller C, Ding L, et al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013;368(22):2059–2074. https://doi.org/10.1056/NEJMoa1301689.
Jiang P, Lee W, Li X, et al. Genome-scale signatures of gene interaction from compound screens predict clinical efficacy of targeted cancer therapies. Cell Syst. 2018;6(3):343-354.e5. https://doi.org/10.1016/j.cels.2018.01.009.
Han H, Cho JW, Lee S, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380-D386. https://doi.org/10.1093/nar/gkx1013.
Szklarczyk D, Kirsch R, Koutrouli M, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023;51(D1):D638–46. https://doi.org/10.1093/nar/gkac1000.
Cerami EG, Gross BE, Demir E, et al. Pathway Commons, a web resource for biological pathway data. Nucleic Acids Res. 2011;39(Database issue):D685–690. https://doi.org/10.1093/nar/gkq1039.
Hu H, Miao YR, Jia LH, et al. AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors. Nucleic Acids Res. 2019;47(D1):D33–D38. https://doi.org/10.1093/nar/gky822.
Rouillard AD, Gundersen GW, Fernandez NF, et al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database (Oxford).2016;2016:baw100. https://doi.org/10.1093/database/baw100.
Li JH, Liu S, Zhou H, et al. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7. https://doi.org/10.1093/nar/gkt1248.
Li R, Qu H, Wang S, et al. CancerMIRNome: an interactive analysis and visualization database for miRNome profiles of human cancer. Nucleic Acids Res. 2022;50(D1):D1139–46. https://doi.org/10.1093/nar/gkab784.
Dragotto J, Canterini S, Del Porto P, et al. The interplay between TGF-β-stimulated TSC22 domain family proteins regulates cell-cycle dynamics in medulloblastoma cells. J Cell Physiol. 2019;234(10):18349–60. https://doi.org/10.1002/jcp.28468.
Nakamura M, Kitaura J, Enomoto Y, et al. Transforming growth factor-β-stimulated clone-22 is a negative-feedback regulator of Ras/Raf signaling: implications for tumorigenesis. Cancer Sci. 2012;103(1):26–33. https://doi.org/10.1111/j.1349-7006.2011.02108.x.
Liang F, Li Q, Li X, et al. TSC22D2 interacts with PKM2 and inhibits cell growth in colorectal cancer. Int J Oncol. 2016;49(3):1046–56. https://doi.org/10.3892/ijo.2016.3599.
Ayroldi E, Cannarile L, Delfino DV, et al. A dual role for glucocorticoid-induced leucine zipper in glucocorticoid function: Tumor growth promotion or suppression? Cell Death Dis. 2018;9(5):463. https://doi.org/10.1038/s41419-018-0558-1.
Hwang J, Haque MA, Suzuki H, et al. THG-1 suppresses SALL4 degradation to induce stemness genes and tumorsphere formation through antagonizing NRBP1 in squamous cell carcinoma cells. Biochem Biophys Res Commun. 2020;523(2):307–14. https://doi.org/10.1016/j.bbrc.2019.11.149.
Roberts AW, Wei AH, Huang DCS. BCL2 and MCL1 inhibitors for hematologic malignancies. Blood. 2021;138(13):1120–36. https://doi.org/10.1182/blood.2020006785.
Ong F, Kim K, Konopleva MY. Venetoclax resistance: mechanistic insights and future strategies. Cancer Drug Resist. 2022;5(2):380–400. https://doi.org/10.20517/cdr.2021.125.
Xie X, Yang W, Zhang W, et al. Tegaserod maleate exhibits antileukemic activity by targeting TRPM8. Biomed Pharmacother. 2022;154:113566. https://doi.org/10.1016/j.biopha.2022.113566.
Redjimi N, Gaudin F, Touboul C, et al. Identification of glucocorticoid-induced leucine zipper as a key regulator of tumor cell proliferation in epithelial ovarian cancer. Mol Cancer. 2009;8:83. https://doi.org/10.1186/1476-4598-8-83.
Ling Q, Li F, Zhang X, et al. MAP4K1 functions as a tumor promotor and drug mediator for AML via modulation of DNA damage/repair system and MAPK pathway. EBioMedicine. 2021;69:103441. https://doi.org/10.1016/j.ebiom.2021.103441.
Sun Y, Zhang D, Guo X, et al. MKK3 modulates JNK-dependent cell migration and invasion. Cell Death Dis. 2019;10(3):149. https://doi.org/10.1038/s41419-019-1350-6.
Wöss K, Simonović N, Strobl B, et al. TYK2: an upstream kinase of STATs in cancer. Cancers (Basel). 2019;11(11):1728. https://doi.org/10.3390/cancers11111728.
Bi L, Jia S, Hu W, et al. Systematic analysis of prognostic significance, functional enrichment and immune implication of STK10 in acute myeloid leukemia. BMC Med Genom. 2022;15(1):101. https://doi.org/10.1186/s12920-022-01251-7.
Sugito N, Heishima K, Akao Y. Chemically modified MIR143-3p exhibited anti-cancer effects by impairing the KRAS network in colorectal cancer cells. Mol Ther Nucleic Acids. 2022;30:49–61. https://doi.org/10.1016/j.omtn.2022.09.001.
Ronchetti S, Migliorati G, Riccardi C. GILZ as a mediator of the anti-inflammatory effects of glucocorticoids. Front Endocrinol (Lausanne). 2015;6:170. https://doi.org/10.3389/fendo.2015.00170.
Ma Y, Yang H, Kroemer G. Endogenous and exogenous glucocorticoids abolish the efficacy of immune-dependent cancer therapies. Oncoimmunology. 2019;9(1):1673635. https://doi.org/10.1080/2162402X.2019.1673635.
Yang H, Xia L, Chen J, et al. Stress-glucocorticoid-TSC22D3 axis compromises therapy-induced antitumor immunity. Nat Med. 2019;25(9):1428–41. https://doi.org/10.1038/s41591-019-0566-4.
Ayroldi E, Migliorati G, Bruscoli S, et al. Modulation of T-cell activation by the glucocorticoid-induced leucine zipper factor via inhibition of nuclear factor kappaB. Blood. 2001;98(3):743–53. https://doi.org/10.1182/blood.v98.3.743.
Berrebi D, Bruscoli S, Cohen N, et al. Synthesis of glucocorticoid-induced leucine zipper (GILZ) by macrophages: an anti-inflammatory and immunosuppressive mechanism shared by glucocorticoids and IL-10. Blood. 2003;101(2):729–38. https://doi.org/10.1182/blood-2002-02-0538.
Lebson L, Wang T, Jiang Q, et al. Induction of the glucocorticoid-induced leucine zipper gene limits the efficacy of dendritic cell vaccines. Cancer Gene Ther. 2011;18:563–70. https://doi.org/10.1038/cgt.2011.23.
Bereshchenko O, Coppo M, Bruscoli S, et al. GILZ promotes production of peripherally induced Treg cells and mediates the crosstalk between glucocorticoids and TGF-beta signaling. Cell Rep. 2014;7(2):464–75. https://doi.org/10.1016/j.celrep.2014.03.004.
Yan M, Hu J, Yuan H, et al. Dynamic regulatory networks of T cell trajectory dissect transcriptional control of T cell state transition. Mol Ther Nucleic Acids. 2021;26:1115–29. https://doi.org/10.1016/j.omtn.2021.10.011.
Tabe Y, Konopleva M. Role of microenvironment in resistance to therapy in AML. Curr Hematol Malig Rep. 2015;10(2):96–103. https://doi.org/10.1007/s11899-015-0253-6.
Acknowledgements
We acknowledge TCGA and GEO database for providing their platforms and contributors for uploading their meaningful datasets. We thank this website (http://www.bioinformatics.com.cn) for the visualization of the pictures. And we thank all the authors for participating and editing this article.
Authors’ Information
XQX is currently a doctoral student in the First Central Clinical School of Tianjin Medical University and works at Shanxi Fenyang Hospital.
Funding
This work was funded by the General Project of National Natural Science Foundation of China (81970180), the Key Science and Technology Support Project of Tianjin Science and Technology Bureau (20YFZCSY00800), Tianjin Key Medical Discipline(Specialty) Construction Project (TJYXZDXK-056B), Key projects of Tianjin Applied Basic Research and Multi-Investment Fund (21JCZDJC01240), Science and Technology Project of Tianjin Municipal Health Committee (TJWJ2022XK018), Science and Technology Project of Tianjin Municipal Health Committee (TJWJ2022QN030), and Tianjin Municipal Natural Science Foundation (22JCQNJC00820).
Author information
Authors and Affiliations
Contributions
MFZ and XQX contributed to study design; XQX, YZL, JXW, and XJ retrieved and analyzed the data; XQX and RS wrote the main manuscript; RS, YZL, MZ, XX, and DNX prepared the figures; MFZ contributed to supervision, review, and revision of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
TCGA and GEO belong to public databases. The patients involved in the database have obtained ethical approval. Users can download relevant data for free for research and publish relevant articles. Our study is based on open source data, so there are no ethical issues and other conflicts of interest.
Consent for publication
Not applicable.
Competing interests
All authors declare that there are no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
. The expression of the TSC22D domain family genes in normal adult HSCs, normal adult tissues, and adult AML tissues using the Gent2 database.
Additional file 2
. The RNASeqexpression data of TSC22D3 and the corresponding clinical prognostic data using the “ TCGA-LAML, NEJM 2013” dataset of the cBioPortal database.
Additional file 3: Table S1
. Analysis of the expression, correlation, and survival prognosis of predicted kinases regulated by TSC22D3 using the UCSCXenaShiny and Harmonizome database.
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
Xu, X., Sun, R., Li, Y. et al. Comprehensive bioinformatic analysis of the expression and prognostic significance of TSC22D domain family genes in adult acute myeloid leukemia. BMC Med Genomics 16, 117 (2023). https://doi.org/10.1186/s12920-023-01550-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-023-01550-7