- Open Access
Screening and analysis of RNAs associated with activated memory CD4 and CD8 T cells in liver cancer
World Journal of Surgical Oncology volume 20, Article number: 2 (2022)
Liver cancer is one of the most common malignant tumors in the world. T cell-mediated antitumor immune response is the basis of liver cancer immunotherapy.
To screen and analyze the RNAs associated with activated memory CD4 T cells and CD8 T cells in liver cancer.
ESTIMATE was used to calculate the stromal and immune scores of tumor samples, which were downloaded from The Cancer Genome Atlas (TCGA). The differentially expressed genes (DEGs) in high and low stromal and immune scores were screened, followed by functional enrichment of overlapped DEGs. We then conducted a survival analysis to identify immune-related prognostic indicators and constructed protein-protein interaction (PPI) networks and ceRNA networks. Finally, chemical small-molecule–target interaction pairs associated with liver cancer were screened.
A total of 55,955 stromal-related DEGs and 1811 immune-related DEGs were obtained. The 1238 overlapped DEGs were enriched in 1457 biological process terms and 74 KEGG pathways. In addition, a total of 120 activated memory CD4 T cell-related genes and 309 CD8 T cell-related genes were identified. The survival analysis revealed that upregulated expression of T cell-related genes including EOMES, CST7, and CD5L indicated the favorable prognosis of liver cancer. EOMES was regulated by has-miR-23b-3p and has-miR-23b-3p was regulated by lncRNA AC104820.2 in the ceRNA network of activated memory CD4 T cell-related genes. In addition, EOMES was regulated by has-miR-23a-3p and has-miR-23a-3p was regulated by lncRNA AC000476.1 in the ceRNA network of CD8 T cells.
T cell-related RNAs EOMES, CST7, CD5L, has-miR-23b-3p, and has-miR-23a-3p may be associated with the prognosis of liver cancer. And the molecular characteristics of these T cell-related genes were plotted.
A total of 309 CD8 T cell-related genes and 120 activated memory CD4 T cell-related genes were screened in liver cancer tumor samples.
Forty-four chemical small-molecule–target interaction pairs associated with activated memory CD4 T cells and 276 pairs associated with CD8 T cells were screened.
Upregulated expression of T cell-related genes including EOMES, CST7, and CD5L indicated the favorable prognosis of liver cancer.
Liver cancer is one of the ordinary malignant tumors, with high morbidity and mortality [1, 2]. Although the incidence of liver cancer has decreased, there are still approximately 840,000 new cases worldwide every year . Surgical resection and liver transplantation are currently recognized as effective methods for the treatment of liver cancer; however, the 5-year survival rate is only approximately 50%, and the 5-year recurrence rate can reach 60~70% . Epigenetic-related molecules (such as SNRPC) , non-coding RNA (such as hsa-mir-221) , and immune-related molecules (such as DCK)  could be used as potential biomarkers for diagnosis, treatment, and prognosis monitoring in liver cancer. Therefore, it is critically necessary to further explore the molecular pathogenesis of this disease by screening for new molecular biomarkers that could enhance early diagnosis and improve current therapies for liver cancer.
According to the competing endogenous RNA (ceRNA) hypothesis, different RNA transcripts participate in pathological processes, primarily by competing for the binding sites of shared miRNAs. Li et al. have suggested that long noncoding RNA LINC01139 promotes the progression of liver cancer by upregulating MYBL2 via competitively binding to members of the miR-30 family . In addition, tumor infiltration of immune cells is closely related to the prognosis of tumors and the determination of immunotherapeutic targets . Tumor immune escape is an important characteristic of tumor formation and is related to the decline of the T cell response-ability . T cells are key mediators of tumor destruction and are critical to the specificity of tumor antigen expression . Moreover, numerous studies have shown that the T cell-mediated antitumor immune response is the basis of tumor immunotherapy, which is associated with a favorable prognosis . Parriott et al. found that T cells expressing the chimeric-PD1-Dap10-CD3zeta receptor can reduce the tumor burden in a variety of mouse solid cancer syngeneic models . CD8+ T lymphocytes are essential factors affecting the efficacy of immunotherapy. Donghua et al. identified CD8+ T cell co-expression genes (C1QC, CD3D, GZMA, and PSMB9 ) that promoted infiltration of CD8+ T cells in liver cancer . Immunotherapy is widely used to treat advanced cancer. However, immunotherapy has many limitations, such as low success rate, many complications, and rapid progression. Studies have shown that the low success rate of immunotherapy is related to less T lymphocyte infiltration. Analysis of infiltrating T cell related molecules may be helpful for immunotherapy.
However, there are few reports on the immune-related prognostic indicators of liver cancer. In this study, we aimed to analyze T cells related to tumor immunity and investigate RNAs associated with T cells in liver cancer to identify useful molecular markers related to the prognosis of liver cancer. In the tumor microenvironment, immune cells and stromal cells are two main types of non-tumor components. They are of great significance for the diagnosis and prognosis evaluation. The immune and stroma scores which was based on the ESTIMATE algorithm are the quantification of the immune and stroma components in the tumor microenvironment. The differential genes were divided into upregulated group and downregulated group according to the intermediate value of quantitative scores. In order to find more representative genes, common stroma and immune genes were screened in the upregulated group and downregulated group, respectively. Then, we analyzed the molecular mechanisms associated with T cells in liver cancer based on the representative genes.
Materials and methods
Data sources and data preprocessing
The RNA expression RNAseq sequencing data (Counts) and clinical data (phenotype) of LIHC (Version: 07-19-2019) were downloaded from the TCGA database  (https://xenabrowser.net/). Fragments Per Kilobase of exon model per Million mapped fragments (FPKM) was used to standardize the original data. The FPKM of a gene in sample was equal to the ratio of the total exon fragments to the reading falling in the genome (map reading (Millions)) and the length of the gene (exon length (KB)). FPKM = Total exon fragments/(Mapped reads×Exon length). The RNA expression RNAseq sequencing data (FPKM) was also downloaded from the TCGA database. A total of 58,387*369 Counts and FPKM representation matrices were obtained. The detection platform of sequencing data was the Illumina HiSeq 2 000 RNA Sequencing platform.
The RNA-Seq was annotated based on the annotation file of the Gencode database  (https://www.gencodegenes.org/). Then, genes with “protein_coding” annotations were extracted as mRNA. Meanwhile, RNAs with “TEC,” “known_ncrna,” “macro_lncRNA,” “bidirectional_promoter_lncrna,” and “lncRNA” annotations were extracted as lncRNA, and other genes were defined as “undefined.” Then, the sequencing data log2 (count+1) values of Counts were restored to the count value. Reads and the sequencing data log2 (count+1) values of FPKM were restored to the fpkm value. In the expression profile analysis of the Ensembl_ID, the mapping probe was used to calculate the gene expression value (obtained from the annotation files of the chip platform and microarray dataset) to Symbol_ID. The average value was taken as the level of Ensembl_ID expression when multiple probes matched one Symbol_ID. Lastly, 369 tumor samples, 60,483 ENSEMBLE IDs, and 58,387 RNA symbols were obtained in Table 1. The clinical information related to prognosis in TCGA was also collected in Table 1, containing overall survival (OS) and OS status, and the lifetime was converted from days to months (days/30).
Immune score calculation
The stromal and immune scores of all tumor samples were calculated with the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE algorithm in R .
Differential expression analysis
The tumor samples were split into high- and low-score groups based on the median value of stromal and immune scores of all tumor samples. The typical Bayesian modified t test in limma package  (Version 3.40.6) was applied to analyze differentially expressed RNA (DERNA) between high and low immune-score groups with a cutoff of |logFC| > 1 and P < 0.05. Lastly, the ggscatter function of the ggpubr package  in the R language (Version: 0.2.2) was used to draw a volcano plot of the DERNAs. The DEGs of stromal and immune-score groups were selected as the DEGs related to the immune microenvironment.
Based on up- and downregulated DEGs, the clusterProfiler package  in R was used to perform Gene Ontology (GO)  and Kyoto Encyclopedia of Genes and Genomes (KEGG)  pathway enrichment analysis with a cutoff of P < 0.05 and count ≥ 2.
Screening of T cell-related genes
RNA-Seq expression profile data were used to target DEGs, and the abundance matrix of the immune cells of the samples was evaluated using the CIBERSORT deconvolution algorithm  to analyze the abundance of infiltration of 22 immune cells in tumor samples. The LM22 dataset provided by the Cibersort website was used as the RNA expression signature template following the analysis of the abundance of infiltration of immune cells, and parameters were set as perm = 50 and QN = N. The landscape of the immune cells was then plotted using the ggplot2 package in R. The Pearson correlation coefficient between the expression value of the DEGs and the abundance of immune infiltration in activated memory CD4 T cells and CD8 T cells was calculated, and the activated memory CD4 T and CD8 T cell-related genes were corrected for BH with a threshold value of |r| > 0.3 and P < 0.05. Moreover, KEGG pathway analysis on T cell-related genes was carried out based on the carcinoma- and hepatocellular-related pathways in the Comparative Toxicogenomics Database (CTD)  (http://ctdbase.org/).
Based on the median expression in T cell-related genes, samples were divided into high- and low-expression groups. The survival analysis was then carried out using a Kaplan-Meier plot. Finally, a log-rank statistical test was conducted with a significance threshold of P < 0.05.
PPI network analysis
The STRING database  was used to analyze the protein-protein interactions encoded by T cells. The PPI score was set as 0.7 (high-confidence value). Afterward, the PPI network of T cell-related genes was constructed using Cytoscape software .
Prediction of miRNA-target interaction pairs
The miRNAs of the T cell-related genes were predicted using miRWalk 3.0 . Then, the miRNA-target interaction pairs in the TargetScan , MiRDB , and MirTarBase  databases were also obtained using a threshold of score > 0.95. In addition, the HMDD V3.2 database  was used to further validate and screen the predicted miRNA with the keywords “Carcinoma, Hepatocellular”.
Construction of the ceRNA network
The lncRNAs-miRNAs relationship between T cell-related genes was predicted using the DIANA-LncBase database . Finally, the T cell-related gene lncRNA-miRNA-mRNA network was constructed using Cytoscape software.
Chemical small-molecule–target network analysis
To search for liver cancer-related genes and chemicals, the Comparative Toxicogenomics Database (http://ctd.mdibl.org/) was searched using “Carcinoma, Hepatocellular” as keywords. Genes that were both associated with liver cancer and belonged to the T cell-related genes ceRNA network were used to screen chemical-target pairs. The T cell-related genes chemical small-molecule–target network was obtained utilizing the Cytoscape software.
The DEGs between the high and low stromal and immune score groups
The aim was to screen immune-related differential genes in liver cancer. In total, 10,221 genes from the TCGA-LIHC dataset were matched by the ESTIMATE algorithm, and the number of unmatched genes was 191. There were 1811 DEGs (of which 1744 and 67 were up- and downregulated, respectively) were selected according to the high and low immune score groups (Fig. 1A and Supplementary file 1). Meanwhile, a total of 55,955 DEGs (of which 2279 and 153 were up- and downregulated, respectively) were selected between high and low stromal score groups (Fig. 1B and Supplementary file 2). In addition, 1211 overlapped upregulated DEGs and 27 overlapped downregulated DEGs were selected in the stromal score and immune score groups (Fig. 1C, D). The names of DEGs are shown in Supplementary file 3. In the enrichment analysis of 1238 overlapped DEGs, the GO analysis showed that the 1238 DEGs were mainly enriched in 1457 GO-biological processes (BPs; e.g., regulation of lymphocyte activation, leukocyte migration, and lymphocyte differentiation) and 74 KEGG pathways [e.g., cell-adhesion molecules (CAMs)] (Fig. 1E, F).
Identification of T cell-related genes
T cells played an important role in liver cancer occurrence and development. The aim was to further analyze T cell-related genes in liver cancer. The abundance of immune-cell infiltration in tumor samples was estimated and the immune-cell landscape was shown in Fig. 2. Meanwhile, according to the Pearson correlation coefficient, a total of 120 activated memory CD4 T cell-related genes (Supplementary file 4) and 309 CD8 T cell-related genes were identified (Supplementary file 5). In the enrichment analysis of activated memory CD4 T cells and CD8 T cell-related genes, activated memory CD4 T cell-related genes were enriched in 406 GO-BPs and 35 KEGG pathways [e.g., CAMs, Th1 and Th2 cell differentiation, and Hematopoietic cell lineage], the Top 8 GO terms and Top 10 KEGG pathways were shown in Fig. 3A and B, respectively. In addition, CD8 T cell-related genes were involved in 596 GO-BPs and 42 KEGG pathways [e.g., Th17 cell differentiation, Th1 and Th2 cell differentiation, and CAMs], the Top 8 GO terms and Top 10 KEGG pathways were shown in Fig. 3C and D, respectively. In addition, based on the CTD database, activated memory CD4 T cell-related genes were involved in 31 liver cancer-related pathways [e.g., CAMs, Th1 and Th2 cell differentiation, hematopoietic cell lineage] (Table 2), and CD8 T cell-related genes were involved in 36 liver cancer-related pathways [e.g., Th17 cell differentiation, Th1 and Th2 cell differentiation, and CAMs] (Table 3).
T cells associated with prognosis of liver cancer
Liver cancer was an often fatal malignant with poor prognosis. The purpose was to further analyze the role of T cell-related genes in the prognosis of hepatocellular carcinoma. A total of 363 sets of survival information from patients with liver cancer was obtained. Moreover, nine activated memory CD4 T cell-related genes were associated with liver cancer prognosis [e.g., eomesodermin (EOMES), glutathione S-transferase (CST7), and adhesion G protein-coupled receptor E2 (EMR2)] (Table 4). In addition, there were 30 CD8 T cell-related genes associated with liver cancer prognosis [e.g., CD5 molecules like CD5L, EOMES, and CST7] (Table 5). Upregulated expression of T cell-related genes including EOMES, CST7, and CD5L indicated the favorable prognosis of liver cancer. Downregulated expression of T cell-related genes including NCF2 and HTAR3 indicated the poor prognosis of liver cancer.
PPI network of the T cell-related genes
The PPI network was constructed to plot the characteristics of these molecules. PPI network analysis of activated memory CD4 T cell-related genes revealed 53 nodes, including one survival-related gene (EOMES), and 162 interaction pairs (Fig. 4A). The CD8 T cell-related gene PPI network contained 127 nodes, including 11 survival-related genes [e.g., EOMES, CD69 molecule (CD69), and zeta chain of T cell receptor-associated protein kinase 70 (ZAP70)], and 613 interaction pairs (Fig. 4B).
CeRNA network of the T cell-related genes
The ceRNA network was constructed to describe the regulatory effect of non-coding RNA on T cell-related differential molecules. Based on the HMDD database, ten miRNA-mRNA relationships of activated memory CD4 T cell-related genes were obtained (10 miRNAs and three target genes). Also, 26 miRNA-mRNA relationships of CD8 T cell-related genes were obtained (22 miRNAs and 10 target genes). Four miRNA-lncRNA relationships of activated memory CD4 T cell-related genes were obtained (four miRNAs and one lncRNA). Furthermore, 21 miRNA-lncRNA relationships of CD8 T cell-related genes were obtained (21 miRNAs and 13 lncRNAs). Based on the ten miRNA-mRNA relationships and four miRNA-lncRNA relationships of activated memory CD4 T cell-related genes, a ceRNA network of activated memory CD4 T cell-related genes was constructed (Fig. 5A). Here, EOMES was regulated by has-miR-23b-3p and has-miR-23b-3p was regulated by lncRNA AC104820.2. In addition, the ceRNA network of CD8 T cell-related genes was constructed among 26 miRNA-mRNA relationships and 21 miRNA-lncRNA relationships (Fig. 5B). Here, EOMES was regulated by has-miR-23a-3p and has-miR-23a-3p was regulated by lncRNA AC000476.1.
Chemical small-molecule–target network analysis of T cell-related genes
The chemical small-molecule–target network was constructed to find agents that regulate these differentially expressed genes. There were 44 chemical small-molecule–target interaction pairs associated with activated memory CD4 T cells (Fig. 6A), including five mRNAs and 26 chemical small molecules. In addition, there were 276 CD8 T cell-associated chemical small-molecule–target interaction pairs, containing 19 mRNAs and 110 chemical small molecules (Fig. 6B).
In this study, a total of 55,955 stromal-related DEGs and 1811 immune-related DEGs were obtained. Then, the 1238 overlapped DEGs were enriched in 1457 BPs and 74 KEGG pathways. In addition, a total of 120 activated memory CD4 T cell-related genes and 309 CD8 T cell-related genes were identified. The survival analysis showed that the T cell-related genes EOMES, CST7, CD5L, and EMR2 were associated with the prognosis of liver cancer. Our analysis recommended the AC104820.2-has-miR-23b-3p-EOMES axis and the AC000476.1-has-miR-23a-3p-EOMES axis from activated memory CD4 T cell- and CD8 T cell-related ceRNAs, respectively.
A total of 120 activated memory CD4 T cell-related genes and 309 CD8 T cell-related genes were identified in this study. Moreover, T cell-related genes were involved in Th1 and Th2 cell differentiation and CAM signaling pathways. Wang et al. found that abnormal expression of CD4 T cell subsets in peripheral blood of patients with ovarian cancer and an imbalance of Th1/Th2 and Treg/Th17 provide a basis for clinical immunotherapy of ovarian cancer . In addition, Su et al. identified lncRNAs and genes related to the prognosis of gastric cancer and showed that prognostic genes positively related to mortality were enriched in the CAM signaling pathway by performing a functional enrichment analysis . Furthermore, Talia et al. showed that GPER levels are related to metastasis and premigration genes belonging to the CAM signaling pathway in breast cancer , which further suggested that T cell-related genes may be related to the progression of liver cancer via Th1- and Th2-cell differentiation and CAM signaling pathways.
EOMES is an important transcription factor that regulates the differentiation and function of effector T cells, and previous studies have shown that the survival of tumor patients is closely related to EOMES expression. For instance, Gao et al. suggested that the potential anticancer functions of EOMES, ATF5, and ECM1 were confirmed by siRNA experiments . The CD8 T cells were divided into three distinct subpopulations: PD1Hi, PD1Int, and PD1-. Ma et al. showed that compared with adjacent non-tumor liver tissues, the PD1Hi CD8 T cells were significantly enriched in tumors, and the PD1Hi CD8 T cells in liver cancer were highly expressed with depletion-related inhibitory receptors (TIM3, ctla-4, and so on) and transcription factors (EOMES, BATF, and so on) . In the present study, EOMES was a T cell-related gene associated with liver cancer prognosis. Moreover, EOMES was regulated by has-miR-23b-3p and has-miR-23b-3p was regulated by lncRNA AC104820.2 in activated memory CD4 T cell-related genes ceRNA network, and EOMES was regulated by has-miR-23a-3p and has-miR-23a-3p was regulated by lncRNA AC000476.1 in CD8 T cell-related genes ceRNA network. Zaman et al. found an inhibitory effect of miR-23b-3p on the expression of the PTEN gene in renal carcinoma . Meanwhile, Wen et al. showed that osthole inhibited EMT-mediated metastasis of prostate cancer by inhibiting snail signaling and miR-23a-3p . Considering these results together, we speculate that EOMES may be the potential target of has-miR-23b-3p and has-miR-23a-3p in liver cancer and that they play important roles in the progression of the disease.
Moreover, in the present study, CST7 was a T cell-related gene associated with liver cancer prognosis. CST7, also known as CMAP, has been reported to have close connections with liver cancer. For instance, Zhou et al. suggested that CST7 and CSTB genes may serve as potential prognostic and diagnostic biomarkers for liver cancer . Interestingly, CMAP is a novel cystatin-like gene involved in liver metastasis . CD5L is a soluble scavenger cysteine-rich protein that regulates the inflammatory response. Aran et al. suggest that CD5L is upregulated in hepatocellular carcinoma and promotes the proliferation and anti-apoptotic response of liver cancer cells by binding to HSPA5 (GRP78) . Thus, we speculate that CST7 and CD5L contribute to liver cancer progression.
In the present study, immune-related prognostic indicators and chemical small-molecule–target interaction pairs associated with liver cancer were screened, which could guide liver cancer clinical decision-making. Moreover, we constructed a ceRNA network of immune-related genes associated with liver cancer prognosis and predicted miRNA targets, which provided a basis for the study of tumor-induced T cell-mediated immune escape. However, we did not use molecular biological experiments to verify our results. The molecular mechanisms that T cell-related genes influence the prognosis of liver cancer have not been investigated, and further molecular biological experiments are needed to explore the specific roles of these genes.
T cell-related RNAs EOMES, CST7, CD5L, has-miR-23b-3p, and has-miR-23a-3p may be associated with the prognosis of liver cancer. And the molecular characteristics of these T cell-related genes were plotted by PPI network, ceRNA network, chemical small-molecule–target network.
Availability of data and materials
The Cancer Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data
- CD5L :
CD5 molecule like
- EOMES :
- CST7 :
- CD69 :
- ZAP70 :
Zeta chain of T cell receptor-associated protein kinase 70
Lee S, et al. Diagnostic performance of CT/MRI liver imaging reporting and data system v2017 for hepatocellular carcinoma: a systematic review and meta-analysis. Liver Int. 2020;40(6):1488–97.
Zhiyong D. Cinobufacini injection for moderate and advanced primary liver cancer: a systematic review and meta-analysis. J Chin Pharm Sci. 2019;28(4):264–75.
Valery PC, et al. Projections of primary liver cancer to 2030 in 30 countries worldwide. Hepatology. 2018;67(2):600–11.
Margonis GA, et al. Association of BRAF mutations with survival and recurrence in surgically treated patients with metastatic colorectal liver cancer. JAMA Surg. 2018;153(7):e180996.
Jianxin X, et al. N6-methyladenosine (m6A) RNA methylation regulator SNRPC is a prognostic biomarker and is correlated with immunotherapy in hepatocellular carcinoma. World J Surg Oncol. 2021;19(1):241.
Zhongjun W, et al. Identification and interaction analysis of key genes and microRNAs in hepatocellular carcinoma by bioinformatics analysis. World J Surg Oncol. 2017;15(1):63.
Wang X, et al. DCK is a promising prognostic biomarker and correlated with immune infiltrates in hepatocellular carcinoma. World J Surg Oncol. 2020;18(1):176.
Li ZB, et al. Long noncoding RNA LINC01139 promotes the progression of hepatocellular carcinoma by upregulating MYBL2 via competitively binding to miR-30 family. Biochem Biophys Res Commun. 2020;525(3):581–8.
Rohr-Udilova N, et al. Deviations of the immune cell landscape between healthy liver and hepatocellular carcinoma. Sci Rep. 2018;8(1):6220.
Dmitry G. Mechanisms and functional significance of tumour-induced dendritic-cell defects. Nat Rev Immunol. 2004;4(12):941–52.
Kishton RJ, Sukumar M, Restifo NP. Metabolic regulation of T cell longevity and function in tumor immunotherapy. Cell Metab. 2017;26(1):94–109.
Ribas A, Wolchok JD. Cancer immunotherapy using checkpoint blockade. Science. 2018;359(6382):1350–5.
Parriott G, et al. T cells expressing a chimeric-PD1-Dap10-CD3zeta receptor reduce tumor burden in multiple murine syngeneic models of solid cancer. Immunology. 2020;160(3):280–94.
Qi Pan Y, et al. Identification of CD8+ T cell-related genes: correlations with immune phenotypes and outcomes of liver cancer. J Immunol Res. 2021;2021:9960905.
Tomczak K, Czerwińska P, Wiznerowicz M. The cancer genome atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol. 2015;19(1A):A68.
Harrow J, et al. GENCODE: the reference human genome annotation for the ENCODE project. Genome Res. 2012;22(9):1760–74.
Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4(1):1–11.
Smyth GK, et al. LIMMA: linear models for microarray data. In: Bioinformatics and computational biology solutions using R and bioconductor. Statistics for biology and health; 2005.
Kassambara A. Ggpubr:“ggplot2” based publication ready plots. R package version 0.1; 2017. p. 6.
Yu G, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.
Ashburner M, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Newman AM, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Davis AP, et al. The comparative toxicogenomics database: update 2019. Nucleic Acids Res. 2019;47(D1):D948–54.
Szklarczyk D, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2016;45(D1):D362–8.
Shannon P, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Dweep H, et al. miRWalk–database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011;44(5):839–47.
Agarwal V, et al. Predicting effective microRNA target sites in mammalian mRNAs. elife. 2015;4:e05005.
Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43(D1):D146–52.
Chou C-H, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302.
Huang Z, et al. HMDD v3. 0: a database for experimentally supported human microRNA–disease associations. Nucleic Acids Res. 2019;47(D1):D1013–7.
Paraskevopoulou MD, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016;44(D1):D231–8.
Wang LH, et al. Th1/Th2 and Treg/Th17 cell balance in peripheral blood of patients with ovarian cancer. Nan Fang Yi Ke Da Xue Xue Bao. 2017;37(8):1066–70.
Su X, et al. Identification of the prognosis-related lncRNAs and genes in gastric cancer. Front Genet. 2020;11:27.
Talia M, et al. The G protein-coupled estrogen receptor (GPER) expression correlates with pro-metastatic pathways in ER-negative breast cancer: a bioinformatics analysis. Cells. 2020;9(3):622.
Gao F, et al. Integrated analyses of DNA methylation and hydroxymethylation reveal tumor suppressive roles of ECM1, ATF5, and EOMES in human hepatocellular carcinoma. Genome Biol. 2014;15(12):533.
Ma J, et al. PD1(hi) CD8(+) T cells correlate with exhausted signature and poor clinical outcome in hepatocellular carcinoma. J Immunother Cancer. 2019;7(1):331.
Zaman MS, et al. Inhibition of PTEN gene expression by oncogenic miR-23b-3p in renal cancer. PLoS One. 2012;7(11):e50203.
Wen YC, et al. By inhibiting snail signaling and miR-23a-3p, osthole suppresses the EMT-mediated metastatic ability in prostate cancer. Oncotarget. 2015;6(25):21120–36.
Zhou X, et al. Investigation of the clinical significance and prospective molecular mechanisms of cystatin genes in patients with hepatitis B virusrelated hepatocellular carcinoma. Oncol Rep. 2019;42(1):189–201.
Morita M, et al. CMAP: a novel cystatin-like gene involved in liver metastasis. Cancer Res. 1999;59(1):151–8.
Aran G, et al. CD5L is upregulated in hepatocellular carcinoma and promotes liver cancer cell proliferation and antiapoptotic responses by binding to HSPA5 (GRP78). FASEB J. 2018;32(7):3878–91.
The authors gratefully acknowledge the multiple databases, which made the data available.
This work was supported by the Welfare Technology Application Research Program of Huzhou (No. 2019GY01).
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.
About this article
Cite this article
Yan, Z., Lijuan, Y., Yinhang, W. et al. Screening and analysis of RNAs associated with activated memory CD4 and CD8 T cells in liver cancer. World J Surg Onc 20, 2 (2022). https://doi.org/10.1186/s12957-021-02461-6
- Liver cancer
- Differentially expressed analysis
- Immune-related genes
- T cells
- Survival analysis
- CeRNA network